求解一维理想磁流体方程的5阶紧凑CWENO中心迎风格式

通过结合求解双曲型守恒律的5阶紧凑CWENO格式和半离散中心迎风格式,推广应用于求解一维理想磁流体力学方程,得到计算一维理想磁流体力学方程的5阶紧凑CWENO中心迎风格式.     

高阶半离散中心迎风方法的研究与应用

结合四阶Central Weighted Essentially Non-Oscillatory格式、三阶Central-Upwind格式构造了一种新的四阶半离散中心迎风差分方法求解双曲守恒律、浅水波方程及有关问题.而且由于数值粘性与时间步长无关,从而在涉及到对流扩散方程的求解时时间步长可根据稳定性需要尽可能的小.     

五阶高分辨率熵稳定算法

针对双曲守恒律方程的数值求解问题,构造一种新型的熵稳定算法.新算法空间方向采用五阶中心加权基本无振荡(CWENO)重构格式,时间方向采用四阶强稳定龙格-库塔(Runge-Kutta)方法.将新算法应用于若干一维Burgers方程和Euler方程组问题数值算例的求解.结果表明:新算法精度高,有效抑制了伪振荡的产生,与理论分析的结果一致.     

一维WENO格式及其在标量守恒方程数值计算中的应用

研究高阶精度加权本质无振荡(WENO)格式及其在标量守恒律方程中的应用.应用五阶WENO空间离散格式和三阶TVD Runger-kutta时间离散格式对一维标量守恒方程以及二维标量守恒方程进行了数值模拟.数值结果表明该格式具有高精度性和本质无振荡性.     

双曲方程的一种二阶TVD差分格式构造方法

本文采用一阶迎风差分格式作Ｔａｙｌｏｒ展开,消去低阶项,给出求解守恒型双曲方程初值问题的一种二阶ＴＶＤ 差分格式的构造方法,并导出可能形式上用于非守恒型双曲方程的差分公式     

W-B-K方程的多辛Preissmann格式

引入正则动量,验证了W-B-K方程具有Hamilton系统多辛格式,并证实此格式具有多辛守恒律、局部能量守恒律和动量守恒律.基于Hamilton空间体系的多辛理论研究了W-B-K方程的数值解法,利用中心Preissmann方法构造离散多辛格式的途径,并构造了一种典型的半隐式的多辛格式,该格式满足多辛守恒律.数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性.     

Level Set方法的研究与应用

本文把与四阶CWENO格式、四阶NCE（Natural Continuous Extensions）Runge-Kutta方法结合之后的Level Set方法应用到一维、二维双曲守恒律标量方程的求解。将所得的数值解与高阶激波捕捉方法所得的数值解进行比较,说明了Level Set方法能很好的处理标量双曲守恒律标量方程的激波追踪问题。     

利用Hamilton-Jacobi方程求解双曲守恒律组的有限元法

将Galerkin高次有限元应用于双曲守恒律组的Hamilton Jacobi方程形式 ,得到了求解一维双曲守恒律组的数值格式。对于标量守恒律方程以及线性双曲方程组 ,这类计算格式具有TVD性质。非线性方程组的计算结果表明该方法具有较好的收敛性     

利用Hamilton-Jacobi方程求解双曲守恒律组的有限元法

将Galerkin高次有限元应用于双曲守恒律组的Hamilton-Jacobi方程形式，得到了求解一维双曲守恒律组的数值格式。对于标量守恒律方程以及线性双曲方程组，这类计算格式具有TVD性质。非线性方程组的计算结果表明该方法具有较好的收敛性。     

具波动算子非线性Schr(o)dinger方程的一种守恒差分格式

研究了一类具波动算子的非线性Schr?dinger方程的数值计算问题.给出了该方程的两个守恒律,构造了求解该方程近似解的一种守恒差分格式,使该差分格式的精度在时间和空间上均达到二阶精度,并对该格式的收敛性及稳定性进行了证明.数值实验与理论结果相一致,很好地验证了本文提出的离散格式.     

求解双曲守恒律方程的高次有限元方法

将Galerkin高次有限元应用于双曲守恒律的Hamilton Jacobi方程 ,得到了求解双曲守恒律的一类数值格式。这类格式在CFL条件下具有TVD性质 ,在更强的条件下 ,其半离散格式的数值解收敛于守恒律方程的熵解。数值结果表明 ,这类格式具有较高的分辨激波的能力。     

求解双曲守恒律方程的高次有限元方法

将Galerkin高次有限元应用于双曲守恒律的Hamilton-Jacobi方程,得到了求解双曲守恒律的一类数值格式,这类格式在CFL条件下具有TVD性质,在更强的条件下,其半离散格式的数值解收敛于守恒律方程的熵解,数值结构表明,这类格式具有较高的分辨激波的能力。     

带浸入边界法的新型五阶WENO格式求解双曲守恒律方程

采用一种带浸入边界法的新型五阶有限差分WENO(weighted essentially non-oscillatory)格式在笛卡尔网格上求解含有复杂物面的双曲型守恒律方程。这种结构网格上的新型WENO格式因对计算网格质量依赖性较高,故一般不能直接应用于上述问题的数值模拟。而浸入边界法是一种能较好处理复杂物面边界的方法。将两种方法结合起来,可在笛卡尔网格上数值解决跨音速复杂流动问题,并用四个经典算例验证新型五阶WENO方法的有效性。     

二维定常对流占优扩散方程的指数变换及相应的有限体积法

利用指数变换消除了二维对流占优扩散方程中源于不对称算子的迎风效应,将对流扩散方程等价转化为守恒型扩散方程,然后通过有限体积法对守恒型扩散方程离散,构造了一种新型的差分格式,而此格式也可理解为利用广义差分法得到,故收敛阶也可以得到证明.数值实验表明,此格式优于以往的几种差分格式,数值解的收敛效果令人满意.     

三维飞船绕流问题的高精度迎风紧致差分格式

从迎风紧致逼近出发,给出了一个求解三维可压Navier-Stokes方程的一种高精度的数值方法,利用Steger-Warming的通量分裂技术将守恒型方程中的流通向量分裂成2部分,在此基础上据风向构造逼近于无粘项的三阶迎风紧致有限差分格式,对方程中的粘性部分采用通常的二阶差分逼近,所建立的差分格式被用来数值求解了三维飞船绕流问题。     

满足3个守恒律的Godunov型格式

为将已有的一维守恒律方程满足多个守恒律的Godunov型格式推广到高维守恒律方程中,对二维的线性传输方程设计了一个满足3个守恒律的Godunov型格式．数值试验表明,该格式具有长时间的保结构性．     

满足两个守恒律的Godunov型格式

茅德康等对一维的守恒律方程设计了满足多个守恒律的Godunov型格式。此格式具有超收敛性和长时间保结构性。为了把这种数值模拟方法推广到高维的守恒律方程中,先考虑二维的线性传输方程,对其设计了一个满足两个守恒律的Godunov型格式。从数值试验可看出,该格式也具有一定的保结构性。     

基函数法在黏性可压缩流动中的研究与应用

研究基函数法在二维和三维轴对称情形下可压缩黏性流动中的应用, 并利用此法数值地计算了二维圆柱的超音速黏性绕流和三维球头的轴对称超音速黏性绕流两个算例。取一阶三角函数为基函数, 构造出导数的中心格式和迎风格式。对于可压缩N-S 方程中的对流项, 采用通量分裂法及中心格式与迎风格式相结合的技术;对于可压缩N-S方程中的黏性项, 则采用中心格式进行处理。由此, 构造出了数值求解黏性可压缩流动的一阶三角函数类型的基函数格式。两个算例的数值计算结果表明, 基函数法不仅在处理无黏可压缩流动时, 是一种高精度、高分辨率的新型计算方法, 而且在处理黏性可压缩流动问题时也一样行之有效。     

双曲型守恒律的一个4阶精度差分格式

基于TVD限制器函数方法选取数值导数,在空间方向用分段3次多项式进行重构,对时间积分用Simpson求积公式,并用四阶Runge-Kutta NCE方法求中间时间点的值,得到求解一维非线性双曲型守恒律方程的4阶精度差分格式;之后给出2个经典数值算例,以验证格式的高精度高分辨率优点。     

时间分数阶扩散方程的一种数值解法

将block-by-block法扩展到分数阶偏微分方程的求解中,即采用block-by-block法离散时间分数阶扩散方程的时间方向,同时采用经典的二阶中心差分格式处理空间方向,得到了新的求解时间分数阶扩散方程的数值格式。数值实验表明,该格式能有效地数值求解一类时间分数阶扩散方程的初边值问题。