论文部分内容阅读
数值模拟流体流动时,常常需要同时求解两个未知函数,并且往往是一个标量函数、一个向量函数,这类问题通常采用混合有限元方法处理。混合有限元方法的理论与应用已经有比较系统的研究,并获得了丰硕的成果。混合有限元方法的突出优点是满足局部守恒性质,其缺点是必须处理源于鞍点方程的、系数矩阵非对称正定的线性方程组(通常的数值代数(例如共轭梯度法)方法不易处理),此外,混合有限元方法计算量较大。本文以一维发展型反应扩散方程为模型,系统研究发展型问题积分形式的流量重构算法。我们将分别研究半离散和全离散的计算格式。在两种计算格式中,空间变量均采用有限元方法逼近,在半离散格式中时间变量不离散,在全离散计算格式中,时间变量采用有限差分方法逼近。本文算法分别计算逼近压力Ph和逼近流量uh。首先,采用标准的Galerkin有限元格式求得逼近压力Ph;然后,在每一个空间剖分区域上,直接构造一个极其简单的逼近流量的局部计算公式,直接得到精确流量u的有效逼近uh。此流量重构算法不仅有效地克服了通常混合有限元方法的上述缺点,大大节省了计算工作量,而且保持了混合有限元方法满足局部守恒性质的优越性。本文严格分析了计算格式的收敛性,数值计算结果显示了该算法的优越性。