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

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

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

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

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

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

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

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

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

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

对流占优扩散问题的一种特征差分方法

用基于一般的 L agrange插值的特征差分方法求解对流占优扩散问题 ,会出现较大的数值扩散或者数值振荡等困难 ,高阶单调插值又计算复杂。该文采用 A.A .Sam arskii构造差分格式的方法 ,建立了一种新的特征差分方法。先对对流扩散方程的扩散项进行修改 ,然后再进行特征差分。此方法具有较高精度 ,并消除了非物理振荡。证明了方法的无条件稳定性。数值结果表明 ,该方法可成功求解对流占优扩散问题。     

对流占优扩散方程一种新的特征差分算法

将特征线法和有限差分法相结合，借助于双线性插值，给出了求解对流占优扩散方程数值解的一种新的特征差分格式，并研究了算法的收敛性。该算法的优点是特别适用于求解变系数的对流占优扩散方程，能更有效地消除数值震荡现象。     

对流占优扩散方程的一种基于斜线性插值的特征差分算法

将特征线法和有限差分法相结合，借助于斜线性插值。给出了求解对流占优扩散方程数值解的一种新的特征差分格式。并研究了算法的收敛性。该算法的优点是特别适用于求解变系数的对流占优扩散方程。能更有效地消除数值振荡现象。     

QUICK格式在对流占优扩散问题中的应用分析

用基于一般有限差分方法迎风格式求解土壤水盐对流占优扩散问题,会出现较大的数值扩散或者数值振荡现象。该文采用有限体积数值计算方法,将QUICK格式引入上述问题的求解过程中。数值结果表明,该方法可以有效的求解对流占优扩散问题,为土壤水盐运移问题求解提供理论基础。     

求解二维非定常对流扩散方程的高精度指数型差分方法

针对二维非定常对流扩散方程,提出了一种高精度指数型差分方法,证明了所构造差分格式的无条件稳定性.通过数值算例验证了差分格式的有效性和合理性,并且对于对流占优问题的求解该方法更优越.     

三维区域上的多重网格算法

对三维椭圆型偏微分方程边值问题,设定其求解区域为曲边六面体,在非均匀剖分条件下,使用多重网格算法,在规则区域和均匀网格下,多重网格方法的实施有其标准的算法流程,而支工程计算中常见的任意几何区域和非均匀网格剖分、多重网格方法的应用相对困难。此时可施行一坐标变换（等参变换）,把物理空间中曲边六面体上的非均匀网格,映射到计算空间中长方体区域上的均匀网格,然后在计算空间中求解相应的偏微分方程边值问题。这种     

三维抛物型流动问题的非错列网格计算方法

用控制容积有限差分法计算三维抛物型流动，通常要把三个速度分量分别定义在三套相互错列的网格上，以消除不皮形场对解的真空性所产生的影响。网格的错列使得计算过程和程度变得相当复杂，造成了应用上的困难。本文利用虚拟网格和插值技术，在网格不错列的的前提下，建立了三维抛物型流动的新的计算方法。对典型流动问题的计算结果与实验数据、近似分析解以及用错列网格方法计算的结果一致。这种方法对简化述化床悬浮空间内的气－…     

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

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

非等距网格上奇异扰动问题的有限差分格式

考虑一维定常对流扩散方程的Dirichlet边值问题,利用Taylor级数构造一个基于非等距网格的有限差分格式,给出了格式的截断误差估计,并分析了其稳定性.采用网格生成函数构造非等距网格,并与一些已有的差分格式对比,数值实验表明该格式可以得到更为精确的数值结果,能很好地模拟边界层效应.     

新型高分辨率格式及其在CFD的应用

 利用小波奇异分析构造了一种新型的高分辨率高效计算格式,借助小波奇异分析所找到的流场中奇异点区域和光滑区域,合理选择高分辨率、高精度WENO格式和高精度中心差分格式进行流场的离散求解;利用该格式对超声速二维前台阶流动问题、跨声速RAE 2822典型翼型的二维外流场、跨声速VKI-LS59涡轮叶栅的二维流场以及NASA Rotor 37跨声速压气机转子和NASA Rotor 67跨声速风扇转子的三维流场进行数值模拟,计算结果表明,采用这种新型格式可以用于内外流复杂流场计算。数值计算还发现,在采用相同网格的情况下,本文格式不仅具有通常WENO格式高分辨率的特点,而且计算效率较WENO格式高出3～5倍。     

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

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