二维泊松方程非均匀网格上的高精度紧致差分格式

提出了数值求解二维泊松方程基于非均匀网格的高阶紧致差分格式,通过选取合适的网格分布参数求解具有边界层的数值算例,空间可以达到四阶精度.并与均匀网格上的计算结果进行比较,充分验证了本文非均匀网格高精度紧致格式的精确性和优越性.     

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

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

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

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

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

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

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

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

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

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

对流扩散问题非均匀网格上的部分半粗化多重网格方法

结合非均匀网格上的 HOC 格式与部分半粗化的多重网格方法对具有边界层的2维对流扩散问题进行了求解,并基于面积率构造了部分半粗化多重网格方法的插值算子和限制算子。数值实验表明:对于只需要在1个方向进行网格加密的边界层问题,基于部分半粗化的网格分布策略及多重网格算法可以大大减少无边界层方向的网格数,从而较完全粗化的网格分布策略及多重网格算法具有更高的计算精度和求解效率。     

一种求解反应对流扩散方程的稳定化谱元方法

针对谱元方法求解二维非稳态反应对流扩散方程中出现的稳定性问题,提出了一种稳定的高精度数值方法。该方法在空间上将Chebyshev谱元方法和一致逼近迎风方法相结合,时间上采用分步θ-格式。通过解析解算例验证了该方法的精度及数值稳定性,并对含有不同类型边界层的反应对流扩散问题进行了求解。研究表明:一致逼近迎风项的增加扩大了谱元方法求解反应对流扩散方程的稳定域,在对流项及反应项占优时保持了数值解的高精度;对于含有边界层的复杂反应对流扩散问题,数值解在整个计算区域内获得了一致收敛的结果。研究工作对谱元方法在反应对流扩散问题高精度数值求解中的应用提供参考。     

基于非均匀网格求解非线性对流扩散问题的一种高精度差分格式

将作者基于均匀网格提出的优化差分法和反演差分法推广到非均匀网格中，提出了一种有效求解定常非线性对流扩散问题的高精度差分格式，在此基础上进一步发展了相应的非定常非线性对流扩散问题的高精度格式。数值实例表明，该格式对对流占优和扩散占优问题均具有较好的适应性，对待求量的大梯度变化有极高的分辨能力，计算结果明显优于传统的差分格式。此格式亦可方例地应用于非均匀网络在计算区域内取所有空间步长相等时的特例－－均匀网络中。在水环境模拟的实际计算中，根据待求量的变化规律合理地调整非均匀网络的疏密分布，不仅增强了高精度差分格式的实用效果，而且可使该格式获得比在含相同结点数的均匀网络系统中更为精确的数值结果。     

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

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

自适应的移动单位分解法求解边界层问题

提出一种新的无网格方法,即自适应移动单位分解法来处理边界层问题,该方法的主要思想是将移动网格技巧与单位分解法结合起来.利用等弧长分配自适应地移动节点分布,在此基础上使用单位分解格式求解一维对流扩散模型,并给出方法的理论基础与算法流程.文中数值实验结果验证了该自适应移动方法的有效性和优越性.     

非线性椭圆问题的紧致差分格式及瀑布两网格法

为了讨论来源于科学工程问题的二维非线性椭圆问题的离散格式及其数值解法。首先,将泊松方程的四阶紧致差分格式推广至二维非线性椭圆问题,提出了紧致差分(CFD)格式,基于CFD格式,选取合适的步长,形成粗网格层和细网格层。在粗网格层上,使用牛顿法求得对应的非线性方程的高精度数值解;在细网格层上,运用插值算子将粗网格上的数值解进行插值,得到细层上较好的初始值,并再次使用牛顿法进行求解,提出了CFD格式下的瀑布两网格(CTG)法。数值实验表明:提出的CFD格式具有四阶计算精度,CTG法迭代步数少、计算时间短。     

随机波动下障碍期权定价的有限差分方法

为有效求解随机波动影响下,障碍期权定价的二维对流扩散方程的初值边值问题,采用非均匀有限差分近似方法,构造了非均匀空间网格,利用泰勒级数展开式导出了非均匀网格上的一阶偏导、二阶偏导以及混合偏导项的差分格式,对离散得到的常微分方程组采用Craig-Sneyd格法迭代求解,通过数值实验将所得结果同蒙特卡洛方法进行了比较.研究结果表明,非均匀有限差分方法是求解障碍期权定价问题的一种稳健、有效的数值方法.     

1维非定常对流扩散方程的有理型高阶紧致差分格式

针对1维非定常对流扩散方程,首先建立了1种2层有理型高阶紧致差分格式,其局部截断误差为O（h4+τ2）。然后采用 von Neumann 分析方法证明了该格式是无条件稳定的。由于在每个时间层上只涉及到3个网格点,因此可直接采用追赶法求解此差分方程。最后通过3个数值算例验证了方法的精确性和可靠性。数值结果表明:所述格式不仅能够适用于非定常对流扩散问题,而且能够较好地求解非定常纯对流问题或纯扩散问题,并且其计算效果均优于 Crank-Nicolson（C-N）格式和指数型高阶紧致（EHOC）差分格式。     

封闭方腔内自然对流问题的高精度紧致差分格式

提出数值求解二维非定常不可压涡量-流函数Navier-Stokes/Boussinesq方程组的高精度紧致差分格式,格式空间为四阶精度,时间为二阶精度,并且是无条件稳定的.为了验证高精度紧致差分格式的精确性和可靠性,对有解析解的二维非定常不可压Navier-Stokes/Boussinesq方程组的Dirichlet问题和典型的封闭方腔自然对流问题进行数值模拟.     

求解对流扩散反应方程的一种高精度紧致差分方法

首先,针对一维对流扩散反应方程,借助截断误差余项修正的方法,将中心差分格式余项中未知函数的三阶和四阶导数项利用一阶导数的表达式来代替,从而提出一种新的紧致差分格式,具有四阶精度.然后,为了简化计算,对格式常系数形式的耗散误差和色散误差进行分析,证实该格式的低耗散性.接着,将该方法推广到二维,运用降维的思想转化成2个一维形式的定常对流扩散反应方程,并用求解一维方程的方法,离散后相加即得二维对流扩散反应方程的紧致差分格式.最后,通过数值实验验证本文格式的精确性和可靠性.     

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

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

基于弹簧近似的非结构化网格自适应处理方法

首先对弹簧近似法进行了说明,并对节点弹簧法和边弹簧法进行了对比.在边弹簧法的基础上,从调节弹簧平衡长度的角度出发,提出了一种非结构化网格的自适应处理方法,并结合有限单元法,给出了网格自适应处理的流程.最后,对一维对流扩散方程、二维线性抛物型方程以及欧拉方程进行了求解.计算结果表明:文中提出的网格自适应处理方法可以使几个问题的网格节点布置更为合理,有效地降低数值误差,提高计算精度,特别对于算例中所涉及的二维非结构化网格,自适应化过程中并没有出现网格重叠现象.     

求解非定常不可压Navier-Stokes方程的一种高精度并行算法

采用一阶投影法,建立了一种基于MPI求解非定常不可压N-S方程的高精度并行算法.该算法在空间上可达到4阶精度,其中,对流项中的1阶导数和粘性项中的2阶导数分别采用WENO格式和4阶对称型宽格式进行离散,而Poisson方程则采用4阶精度的紧致格式进行迭代求解.通过对2维Taylor涡列和双周期双剪切边界层流动问题及3维回转体绕流问题的数值计算,验证了算法的可靠性及其并行效率.     

一维和二维溃坝波的混合有限分析解

针对水深平均二维浅水方程高分辨率模型，采用非交错网格下的二维对流扩散方程的混合有限分析法，运用SIMPLE算法结合Rhie and Chow的动量插值技术计算二维溃坝水流．一维理论解和其他格式的二维数值结果符合良好，表明非交错网格下的混合有限分析法能正确有效地求解溃坝问题．进而对部分溃决情况下的溃坝波传播和绕流问题进行了数值模拟，给出了精细的数值结果，揭示了复杂的运动特性，进一步表明混合有限分析法具有较强的溃坝波的捕捉能力，是求解溃坝流动的有效方法之一．