稳定性可保证二阶格式在多重网格中的有效性和经济性

基于一种稳定性可保证的二阶差分格式（SGSD），对SIMPLE算法实施了完全多重网格循环以加速外迭代的收敛．采用规正变量的方法实施了SGSD．通过对二维顶盖驱动流动的计算，分析了多重网格在SIMPLE算法中的收敛特性．计算结果表明：SGSD格式具有与其他高阶格式及高阶组合格式相同的计算精度，且收敛速度优于其他高阶格式，在雷诺数较高时（Re=3000），其收敛速度是二阶迎风格式的1．77倍，是QUICK格式的1．37陪，同时在疏密网格层次上均可以保证计算的稳定性；采用多重网格加速SIMPLE算法的迭代时，不仅要考虑多重网格的循环方式，还要考虑对流项的离散格式，在计算中SGSD格式具有明显的优势。     

QUICK格式在非结构化网格中的应用

在有限容积法的基础上发展了非结构化网格的对流项二次迎风插值（QUICK）格式．详细推导了扩散项采用格林函数法，对流项采用改进的QUICK格式离散方程，对项盖驱动流和圆柱绕流问题进行了计算，讨论了不同Re下计算的准确性和格式的收敛性，并与高精度结构化网格计算结果进行对比分析．结果表明，该格式的临界网格Peclet数为8／3左右，与中心差分相比较，该格式计算精度相当，对流稳定性好，收敛速度高．同等条件下较结构化网格对复杂区域的模拟更接近实际测量结果，是一种对复杂区域计算有应用前景的对流格式．     

对流项二次迎风插值格式在非结构化网格中的应用

在有限容积法的基础上发展了非结构化网格的对流项二次迎风插值（QUICK）格式．详细推导了扩散项采用格林函数法、对流项采用改进的QUICK格式的离散方程，对顶盖驱动流和圆柱绕流问题进行了计算，讨论了不同Re下计算的准确性和格式的收敛性，并与高精度结构化网格计算结果进行对比分析．结果表明，该格式的临界网格Peclet数为8／3左右，与中心差分相比较，该格式的计算精度与其相当，对流稳定性好，收敛速度高．同等条件下较结构化网格对复杂区域的模拟更接近实际测量结果，是一种对复杂区域计算有应用前景的对流格式．     

对流扩散方程不同格式的数值稳定性分析

利用数值计算方法对一维、线性、无源、两点边值问题的对流扩散方程在均分网格下不同格式的稳定性进行了分析，通过理论推导和实际计算，得出了上述条件下中心差分格式(CD)、QUICK格式和稳定性可控(SCSD)格式的P△cr数，并对稳定性可以保证的(SGSD)差分格式进行分析探讨，验证了数值稳定性是格式固有的属性。在此基础上，二维问题进行数值计算，以资为复杂的多维问题对流离散格式的稳定性分析提供依据。     

网格Peclet数和网格尺寸对对流扩散方程差分格式精度的影响

用泰勒级数展开法对对流扩散方程中的扩散项进行了理论分析,从而证明对于给定的差分格式,不仅网格Ｐｅｃｌｅｔ数不能任意选取,而且网格尺寸也不能随意设定。用四种格式对一维有源对流扩散方程进行了计算,结果表明网格尺寸对差分格式精度的影响比网格Ｐｅｃｌｅｔ数更明显。为了得到真实可靠的结果,所用差分格式的阶数愈高,相应的网格尺寸就必须愈小。二阶以上的高精度迎风差分格式与二阶迎风格式相比,并无明显的优势,建议一般不必要采用二阶以上的高精度差分格式。     

网格Peclet数和网格尺寸对流—扩散方程差分格式精度的影响

用泰勒级数展开法对流－扩散方程中的扩散项进行了理论分析,从而证明对于给定的差分格式,不仅网格Ｐｅｃｌｅｔ数不能任意选取,而且网格尺寸也不能随着设定,用四种格式对一维有源对流－扩散方程进行了计算,结果表明网格尺寸对差分格式精度的影响比网格Ｐｅｃｌｅｔ数更明显,为了得到真实可靠的结果,所用差分格式的阶数愈高,相应的网格尺寸就必须愈小,二阶以上的高精度迎风差分格式与二阶迎风格式相比,并无明显的优势,建议     

QUICK和乘方格式在顶盖驱动方腔流动数值计算中的比较

推导了非均分网络系统中QUICK格式附加源项的表达式，并用QUICK和乘方两种差分格式，分别对不同网格数下的二维、三维顶盖驱动方腔流动进行了数值计算。通过与文献中基准解的比较，考察了两种对流差分格式的数值预测性能。     

二阶迎风有限体积法方腔流数值模拟

介绍了一种基于非结构网格的二阶迎风有限体积离散格式，在对流项的离散过程中，为了达到二阶精度，在界面上对物理量咖作Taylor展开，在处理展开项中的梯度时经入NND格式的优点，克服了中心差分格式不稳定的缺点，并对方腔流动进行了系统的数值模拟．计算网格采用三角形网格，节点数为12960，单元数为25600．压强修正基于SIMPLEC方法．给出了雷诺数Re达到5000时的定常流动结果，与以往方腔流计算的标准解非常吻合，但所用的网格数要少．计算结果表明该二阶迎风有限体积算法有很好的收敛性和稳定性，且松弛因子对计算结果影响很小．     

扭曲网格上扩散方程的一个线性精确差分格式

本文讨论了构造扩散方程差分格式的线性精确方法,它要求当控制方程的解析解是差于自变量的线性函数且扩散系数是常数时所构造的差分格式是精确的.在此基础上推导了在结构四边形网格上求解扩散方程的线性精确差分格式.数值算例表明我们构造的格式在许多扭曲严重的网格上获得了二阶或接近二阶的计算精度.     

不同差分格式在同位网格系统中的计算效果比较

针对流体流动数值计算的有限差分法,系统地研究了离散对流项的6种差分格式:CDS、FUDS、HDS、PLDS、SUDS和QUICK·比较计算采用同位网格系统·采用两个有分析解或基准解的算例,就不同格式对数值求解NS方程的精度、稳定性和收敛特性的影响进行了分析比较·计算结果表明,当扩散项占主导地位时,所有格式在同位网格中几乎具有相同的计算精度·随着对流项的增加直到占主导地位,FUDS、HDS和PLDS的在同位网格中具有相同的精度,而SUDS和QUICK的精度比前三种高,CDS次之·对于相同的速度、压力松弛因子和收敛准则,各种格式在同位网格中的收敛速度相差甚微·     

流动模拟二阶精度的简易保持研究

首先分析了各类低阶和高阶差分格式,指出二阶精度格式具有综合优势．接着以二维对流扩散方程为例,推导出延迟修正的CDS格式,并定义为DCDS格式．通过标量绕静点输运、自然对流及混合对流等实例,并与UDS,CDS及QUICK格式比较,论证了DCDS格式不仅能维持计算结果的二阶精度,而且极大地约束了CDS格式所固有的振荡和越界行为,边界处理也很简单．总之,在程序的效率、调试、边界处理、计算精度及稳定性等因素权衡下．DCDS格式是保持二阶数值精度的简易途径。     

二维非定常对流扩散方程的对角占优隐格式及多重网格算法

提出了一种数值求解二维非定常变系数对流扩散方程的对角占优、空间为二阶精度的隐格式。利用Fourier分析方法证明了该格式是无条件稳定的．由于格式具有对角占优性,因此适用于大梯度（高雷诺数）问题的数值求解．另一方面,为了克服传统迭代法在求解隐格式时收敛速度慢的缺陷,采用了多重网格加速技术,大大加快了迭代收敛速度,提高了求解效率．数值实验结果证明了该格式的精确性、稳定性和对高网格雷诺数问题的强适应性．     

求解三维非定常对流扩散方程的隐式差分方法

提出了一种数值求解三维非定常变系数对流扩散方程,对角占优、空间为二阶精度的隐格式,利用Fourier分析方法证明了该格式是无条件稳定的,并且由于格式具有对角占优性,因此适合于大梯度(高雷诺数)问题的数值求解.另外,为了克服传统迭代法在求解隐格式时收敛速度慢的缺陷,采用了多重网格加速技术,大大加快了迭代收敛速度,提高了求解效率.数值实验结果证明了该方法的精确性、稳定性和对高网格雷诺数问题的强适应性.     

非均匀网格上求解对流扩散问题的高阶紧致差分方法

基于非均网格上函数的泰勒级数展开,推导出求解一维对流扩散问题的高阶紧致差分格式.对于离散化得到的代数方程组,采用BiCGStab（2）迭代法求解.数值实验表明,该格式对于扩散占优、对流占优及边界层问题都有很好的适应性,对于数值模拟待求物理量的大梯度变化具有很高的分辨率,计算结果明显优于传统的均匀网格上的差分格式.在具体的数值模拟中,可根据实际物理量的变化规律,选取适当的网格生成变换函数,合理地调整非均匀网格的疏密分布,从而获得比在含相同结点数的均匀网络系统中更为精确的数值结果.