非结构四面体网格上三维非线性扩散方程的有限体积法

为了有效地数值模拟科学和工程中有广泛应用的非线性扩散方程,在三维线性扩散方程非结构四面体网格的有限体积法的基础上,提出了一个计算非结构四面体网格上非线性扩散方程的有限体积法。方法采用网格单元中心作为计算节点,相对于网格点的方法,计算量减少了一半。用L agrange因子法得到网格点上的值,考虑了网格中心点和网格点的相对位置,更适应大变形的网格。利用算子分裂,使计算更加简单。用N ew ton-B iCG STAB法来求解得到非线性方程组。数值结果表明:该方法具有二阶精度、保持通量守恒、对大变形的网格适应性强。