论文部分内容阅读
本文主要研究各向异性扩散问题的多边形/多面体上的节点型保正有限体积格式.首先,我们构造了扩散问题的一种任意多边形上的非线性节点型保正格式.该格式以节点未知量为主未知量,单元中心未知量和边中点未知量为辅助未知量.它是国际上第一个纯节点型的非线性保正有限体积格式.与绝大多数保正格式采用动态模板不同,该格式的构造基于一个固定模板,不需要联合法向量的凸分解,避免了繁琐的搜索过程,而且辅助未知量的插值非常简单.此外,该格式不需特殊处理就能天然地克服所谓的数值热障问题.我们严格分析了格式的保正性和截断误差.为了提高计算效率,我们还讨论了 Picard迭代方法及其Anderson加速算法.数值实验表明,该格式在应用于一些大变形网格和间断/各向异性张量问题时具有二阶精度和良好的保正性.此外,我们也验证了Anderson加速的有效性.其次,我们将上述节点型格式推广到三维问题,给出了三维扩散问题的一般多面体网格上的非线性节点型保正格式.该格式继承了二维格式优点,且适用于单元面为非平面的网格.由于对偶面上的单侧流的离散是在主单元的一个固定的四面体子单元上完成的,使得算法在多面体网格上实现非常简单.随后,我们将节点型保正格式应用到二维二/三温辐射扩散系统的数值模拟中.数值结果表明该格式能够很好地模拟具有强间断性,强非线性,强耦合性的辐射扩散系统.最后,我们构造了扩散问题的多边形网格上的一类解耦的保正的离散对偶有限体积(DDFV)格式.与已有的离散对偶有限体积格式相同,我们分别在主网格和对偶网格上建立两套有限体积方程组.节点未知量方程组的系数矩阵是对称正定的,而中心未知量方程组的系数矩阵的转置是一个M-矩阵.对求得的节点未知量采取某种截断技巧可使得两组未知量都满足保正性质.该格式在对偶网格上的局部守恒性是有条件成立的,而在主单元上是无条件成立的.与已有的非线性保正格式不同,该格式对于线性问题不需要非线性迭代.对于非线性问题,由于保正机制和非线性迭代方法解耦,可以自由选取非线性求解器.此外,我们严格证明了格式的保正性以及线性问题的适定性.为了数值对比,我们还构造了一种基于插值截断的非线性中心型保正格式.数值实验表明,该解耦格式具有很高的计算效率;它在保正的同时还能够有效地避免数值热障问题.