边界层对流扩散方程在自适应网格上的高精度紧致格式

对流扩散方程在非均匀网格上的高精度紧致格式具有精度高、模板小等优势,然而现有方法往往需要事先指定边界层或大梯度的位置,利用网格分布函数生成非均匀网格并调整网格分布参数,这严重影响方法的适用性.本文提出了一种在二维计算区域上生成正交网格的自适应方法,根据物理解的特征对网格的分布进行自适应调整,并结合非均匀网格上的高精度紧致格式对一维及二维边界层对流扩散问题进行求解.数值结果表明,自适应网格方法结合高阶紧致格式可以有效求解边界层问题,提高数值解的精度,减少计算所需的网格和计算量.     

一种非均匀网格上的高精度紧致差分格式

提出了一种形式简单、网格剖分灵活、具有一定通用性的非均匀网格上的三点四阶紧致差分格式,对格式的截断误差进行了分析.采用文中提出的格式对Burgers方程和对流方程进行数值求解,并与均匀网格上的三点四阶紧致差分格式所得数值解对比,结果证明本文提出的格式对于大梯度问题的数值模拟有更高的精度.     

对流扩散反应方程基于坐标变换的高阶紧致差分格式

基于非均匀网格,提出了一种求解一维定常对流扩散反应方程的高精度紧致差分格式。首先采用坐标变换方法将原方程由物理空间的非均匀网格转换为计算空间的均匀网格,然后给出一阶导数和二阶导数在均匀网格上的中心差分逼近式,并结合变换后的方程,得到了定常对流扩散反应方程具有四阶精度的紧致差分格式。最后,通过数值算例验证了该方法的精确性和高分辨率的特点。数值实验结果表明,对于所研究问题,该方法较不进行坐标变换而直接在物理域上建立的非均匀网格上的高阶紧致格式具有更高精度。     

非均匀网格上三维对流扩散方程高精度紧致差分方法

利用降维法推导出非均匀网格上三维对流扩散方程的高精度紧致差分格式，对于离散得到的代数方程组采用BiCGStab（2）迭代法求解．数值算例表明，在网格节点数相同的情况下，基于非均匀网格的计算格式较均匀网格格式具有高精度、高分辨率的优点，对于合边界层的对流扩散问题有很好的适应性．     

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

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

三维双调和方程的高阶紧致差分格式及其多重网格方法

提出一类求解三维双调和方程的高精度紧致差分格式.该类格式是以泊松方程的高精度格式为基础的四阶精度19点紧致差分格式和六阶精度27点紧致差分格式.采用多重网格方法求解由高精度紧致差分格式所形成的代数方程组,并与低精度方法进行比较.讨论多重网格方法中不同松驰算子的迭代收敛效果.数值实验结果验证四阶紧致差分格式和六阶紧致差分格式的精度以及多重网格方法的可靠性和高效性.     

对流扩散方程非均匀网格上四阶紧致格式的预条件迭代法

利用带填补数的不完全LU分解(ILUT(τ,s))作预处理器以及FGMRES(20)作迭代加速器,对非均匀网格上二维对流扩散方程的高精度紧致差分格式进行数值实验,并与均匀网格上的计算结果进行对比,数值结果显示出非均匀网格上本文方法的优越性,在合适的网格伸缩系数下,本文方法不仅能够保证格式的四阶精度,而且降低了误差的数量级.同时,比较了预条件迭代法与传统迭代法的求解效率,结果表明预条件方法的单位对数残差几乎成直线下降,相比传统迭代法有明显的计算优势.     

紧致差分格式求解二维非稳态不可压N-S方程

采用泰勒展式系数匹配的方法构造基于非等距网格的紧致差分格式并得出了它的截断误差.紧致差分格式能够很好的模拟不同时刻流场的变化情况,网格系统的选择对精度的影响很大,基于非等距网格的紧致差分方法是一种比经典差分方法精度更高的求解非稳态纳维斯托克斯方程的有效算法.     

求解热传导方程的四阶有限容积紧致格式

对非齐次热传导方程提出了一种数值求解的有限容积紧致格式,该格式具有空间上的四阶精度,且与有限差分紧致格式相比,更好地保持了问题的物理守恒性.数值算例表明,在相同的结点下,有限容积紧致格式比有限差分非紧致格式的精度更高.     

构造对流扩散方程高精度非均匀网格格式的指数变换方法

提出了一种基于非均匀差分网格,构造求解对流扩散方程的高精度格式的指数变换方法.引入指数函数,将对流扩散方程变换为扩散反应方程,消除了数值求解中较难处理的对流项.采用优化差分方法推导出扩散反应方程基于非均匀网格的高精度差分格式,进而通过逆变换得到对流扩散方程的高精度格式.理论分析表明,该方法具有3至4阶精度,当计算区域为均匀网格时取得4阶精度.数值实例表明,在相同的非均匀网格系统中,此方法的计算精度明显优于传统的隐式差分方法.在水环境的实际模拟计算中,根据物理量的变化规律灵活地调整非均匀网格的间距,不仅能增强高精度差分方法的实用性,而且可以取得比均匀网格方法更为精确的计算结果.     

奇性方腔黏性流动的完全高精度紧致差分方法

以涡量流函数形式的Navier-Stokes(N-S)方程为例,详细介绍了构造完全高精度紧致差分格式的一般方法.所建立的高精度差分格式,无论是在计算区域的内点还是在边界点上均可以达到4阶精度,且具有紧致性,与已有数值实验结果相比只需要用很少的网格(61×61)就可以求得较高计算精度的数值解,从而大大节省了计算时间,提高了计算效率.     

紧致差分格式的分辨率与精度的实例比较与讨论

通过比较先前建立的4阶最优紧致差分格式,以及传统的6阶和8阶紧致差分格式,来研究精度和分辨率之间的关系,主要比较了空间离散格式的有效波数范围、实际数值计算精度、以及对小尺度波动的模拟能力.数值试验结果表明:(1)这3种格式的计算精度都可以达到理论精度,并且此时精度越高,误差越小;(2)对于小尺度波动,最优4阶紧致格式比6阶和8阶紧致格式具有更高的分辨率;(3)对于行波问题,最优4阶紧致格式能够更加准确地模拟波动的传播行为.理论分析和数值算例的比较结果均表明,数值格式的精度和分辨率并不能互相替代,而是要根据计算问题的需要选择具有合适的精度和分辨率的数值格式.     

GAMMA格式及DATE方法在SIMPLE算法中的应用

为了提高N-S方程数值计算的精度并解决压力-速度的耦合问题,将GAMMA格式、DATE方法与SIMPLE算法相结合,得到了一种高精度非结构同位网格下的SIMPLE算法.该算法用GAMMA格式离散对流项,用DATE方法处理方程求解过程中压力-速度的耦合问题,可以弥补传统动量插值法不能处理强非线性压力场的不足,保证对流项离散的稳定性及高精度特性.     

三对角四阶跳点紧致格式优化和初步应用

 根据修正波数应在充分大的波数范围内接近准确波数的思想,构造了优化的3对角4阶跳点紧致差分格式及插值格式.优化跳点紧致格式仍然具有4阶精度,但提高了分辨率,能够在更大的波数范围内保持群速度特性.通过实验数据表明,优化跳点紧致差分格式的分辨率可达到0.86π,优化紧致插值格式可达0.63π,可较好保持群速度的最大波数为0.75π,均大于标准4阶和6阶跳点紧致格式.分别用优化格式,标准4阶和6阶跳点紧致格式计算小尺度波动的性能,结果表明优化格式在模拟小尺度波动,在减小计算误差并保持群速度方面具有明显优势.     

三对角四阶紧致差分格式的优化和初步应用

 提高数值解的精度和分辨率,有助于更精确地求解日趋复杂的工程问题。本文依据差分格式的伪波数应该在尽可能大的波数范围内接近物理波数的思想,构造了满足四阶精度的具有高分辨率的三对角紧致差分格式。一方面,它可以与近些年发展的求解(循环)三对角方程组的高效算法相结合,以更高的分辨率、更小的计算量来计算一阶导数;另一方面,与传统格式相比,该格式的最大精确求解波数可以达到2.5761,大于传统格式的1.13097。因此,优化格式更适合模拟小尺度波动。数值计算结果表明：(1) 虽然优化格式仍然是四阶精度,但要比传统四阶紧致差分格式的计算误差小,尤其对于小尺度波动,优化格式的计算误差会更小;(2) 对于行波问题,优化格式能够更加准确地模拟波动的传播行为,其优势也更加明显。理论分析和数值算例的比较结果均表明,优化的紧致差分格式更适合求解小尺度波动问题。     

基于高阶紧致格式的二维不可压NS方程求解

本文介绍了一种基于原始变量的用于求解二维非定常不可压Navier-Stokes方程的高阶紧致格式。这种紧致格式最初是用于计算声学（CAA）的高精度格式,相对于传统的紧致格式,使用该格式的优点在于减少计算量的同时降低了边界模板的处理难度。这种方法建立在非交错网格上,空间离散具有六阶精度。压力Poisson方程基于九基点模板的四阶紧致格式进行离散,超松弛迭代进行求解。时间推进上采用四阶Runge-Kutta方法。为验证该方法的精度和有效性,利用该格式计算了一个具有解析解的问题,以及二维非定常情况下的方腔驱动流动问题,并且和传统的紧致格式进行了计算时间的对比。     

通风分裂算法在有限元非结构网格中的推广

将在结构网格中发展起来的通风分裂算法推广到有限体积非结构网格中,并成功地进行了定常跨声速二维欧拉方程的数值计算．此格式具有流通量矢量分裂格式的简单性及通量差分分裂格式的高阶精度．最后给出的ＮＡＣＡ００１２翼型绕流的结果与实验结果符合的较好．     

满足熵增原则的高精度高分辨率格式

提出了构造高精度、高分辨率格式的一条新途径。用此途径构造的格式能满足保证格式具有某些所希望性能的一些原则和要求（例如熵增原则,稳定性原则,精度阶数,格式点数等）。作为例子,文中提出了一个三点三阶精度的紧致格式,它满足熵增原则和稳定性原则。３个例子的数值结果表明该格式能无振荡地捕捉激波并在光滑区具有高的精度。