基于物面弱边界条件的并行全隐式流场求解方法

提出了一种基于物面弱边界条件的Navier-Stokes方程并行全隐式求解方法。首先,实现了非结构网格格点格式有限体积法的粘性物面弱边界条件。以此为基础,推导了相应的粘性Jacobian矩阵。相比传统的物面强边界条件,弱边界条件的应用简化了该推导过程。接着,为了解决Spalart-Allmaras(S-A)模型的收敛和保正性问题,提出了基于物面弱边界条件的S-A模型无条件保正性收敛方法,通过构造特别的隐式矩阵同时保证S-A模型的收敛和湍流工作变量的正值。流场方程全隐式时间推进和无条件保正性收敛方法所得到的线性方程组采用多色高斯-塞德尔迭代法进行求解,并由MPI实现并行化。通过若干二维和三维算例对本文所提出的方法进行了测试,其结果表明:1)随着罚系数的增大,基于弱边界条件的结果趋向于基于强边界条件的结果;2)对于力的计算,弱边界条件自动满足通量一致关系而无需特别处理。

二维非定常对流扩散方程的隐式多重网格方法

提出了数值求解二维非定常对流扩散方程的一种新的加权平均隐格式,利用Fourier分析方法证明了该格式是无条件稳定的.利用多重网格加速技术,提出了基于时间修正的多重网格的全近似格式(FAS),从而克服了传统迭代法在求解隐格式时收敛速度慢的缺陷,大大加快了迭代收敛速度,提高了问题的求解效率.数值计算结果表明,修正的多重网格FAS格式较传统的多重网格粗网格校正格式(CS)具有更好的收敛效率.并且随着σ和s的增大,它可以将传统迭代法的收敛速度提高几十倍,甚至几百倍.

一维非定常对流扩散方程非均匀网格上的高精度紧致差分格式

本文在非均匀网格上给出了求解非定常对流扩散方程的一种高精度紧致差分格式,特别适合边界层和大梯度等问题的求解.从稳态对流扩散方程入手,首先,基于非均匀网格上的泰勒级数展开对空间导数项进行离散,然后对时间项采用二阶向后欧拉差分公式,从而得到一维非定常对流扩散方程在非均匀网格上的三层全隐式紧致差分格式.新格式在时间具有二阶精度,空间具有三到四阶精度,并且是无条件稳定的.最后,通过数值实验验证了本文格式的精确性,以及在处理诸如边界层和大梯度问题上的优势.

带有一个噪声的随机微分方程的全隐式格式（英文）

通过将It型随机微分方程转换为等价的Stratonovich型和向后It型随机微分方程,构造出分别与Stratonovich型和向后It型随机微分方程相容的两种全隐式随机数值格式.这两种数值格式应用于带一个噪声的随机微分方程,且均为均方1阶收敛的.

二维非稳态对流扩散方程的高阶紧致差分格式

针对二维非稳态变系数对流扩散方程,对时间的离散分别采用二阶和三阶向后差分公式,对空间的离散分别采用四阶紧致差分和六阶紧致差分方法,提出了两种高精度紧致差分格式,两种格式的截断误差分别为O(τ^(2)+h^(4)_(x)+h^(4)_(y))和O(τ^(3)+h^(6)_(x)+h^(6)_(y)),并且均是无条件稳定的,最后给出了数值算例验证了理论结果.

Cahn-Hilliard方程的隐显BDF2方法

Cahn-Hilliard方程作为一类重要的四阶扩散方程已成为偏微分方程研究领域一个倍受关注的问题.本文考虑带有Neumann边界的Cahn-Hilliard方程的隐显BDF2半离散格式和全离散格式,并证明了该格式是质量守恒的.