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

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

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

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

水气两相均质瞬变流二阶有限体积法数学模型

采用有限体积二阶Godunov格式对水气两相均质流瞬变现象进行建模和模拟研究.在均质流控制方程中引入动态摩阻,采用二阶Godunov格式求解,运用MUSCL(守恒定律的单调上游中心格式)方法进行二阶线性重构,同时为避免虚拟振荡引入MINMOD(斜率限制器函数);提出双虚拟边界法进行边界计算,实现了计算区域的所有控制体包括边界同时达到二阶精度且计算更简便.将所建模型计算结果与已有模型计算结果、试验数据进行对比,并对模型各参数进行敏感性分析,结果表明:本文模型能更好地模拟试验压力,并有效避免虚假数值振荡,提高了计算结果稳定性和精确度;随着初始气体含气率的增加,波速降低,水气两相均质流的瞬变压力峰值逐渐减小、波动周期延长.     

2类高阶格式数值测试和比较

通过求解二维Burgers方程初边值问题、二维周期漩涡问题、Richtmyer Meshkov(简称RM)不稳定性问题和Rayleigh-Taylor(简称RT)不稳定性问题,对五阶FD-WENO格式（简称WENO5）及二阶Godunov格式MUSCL进行了数值测试和比较.所得结果对于求解气体动力学问题及辐射流体力学问题时,怎样恰当地选择数值方法具有一定参考作用.     

求解一阶线性双曲型偏微分方程组的一个差分格式

利用有限差分法,构造一个求解一阶线性双曲型偏微分方程组的高精度隐式差分格式.该格式的相容性和收敛性被证明.进一步,Fourier级数法分析表明该格式无条件稳定.另外,该隐式差分格式可以转化为显格式求解,较同类格式的计算量小.通过实例验证,该格式数值结果绘制的图形稳定、光滑、效果良好.     

对流扩散方程差分格式稳定性分析

用Fourier方法分析了离散线性对流扩散方程一些差分格式的稳定性和其截断误差.在这些格式的基础上,给出一个新的跳点格式,该格式具有更优的计算效率,数值实验结果与理论分析结果一致.     

满足最大值原理的熵格式计算线性传输方程

茅德康等发展了熵格式计算一维双曲守恒型方程,熵格式具有超收敛性并且适合长时间计算.但是熵格式不满足最大值原理,在最值点处会出现过高或者过低现象.发展了满足最大值原理的熵格式并且对一维和二维线性传输方程进行了数值模拟,数值结果表明该格式在最值点不会出现过高过低现象而且不会发生非物理振荡.     

线性与非线性发展方程差分格式计算稳定性的比较分析

针对线性与非线性发展方程的几种差分格式,以一维线性和非线性平流方程为例,对线性与非线性发展方程差分格式的计算稳定性进行了比较分析,揭示了差分格式结构和初值形式与计算稳定性的关系.理论分析和数值试验证明,线性与非线性发展方程差分格式计算稳定性在本质上是完全不同的.     

双曲型线性方程三阶和四阶TVD格式的新构造

利用Taylor级数理论和总变差减小(TVD)格式的充分条件构造了时间二阶、空间五点三阶和四阶新TVD格式．给出了新TVD格式与传统TVD格式及近期建立的二阶新TVD格式用于线性双曲型方程的计算结果，表明本文新格式特别是四阶TVD格式具有比二阶新TVD格式和传统TVD格式峰值衰减更慢、间断更陡，而计算工作量具有与传统二阶TVD格式相当的良好数值性能．     

线性传输方程满足3个守恒律的差分格式

考虑单个线性传输方程,对其设计了一种满足3个守恒律的差分格式.此格式为3阶Godunov型的,用的是分片2次重构,重构函数的系数由3个守恒量来确定.虽然微分方程是线性的,但所设计的格式是非线性的.数值实验结果表明,此格式是非线性稳定的,并且对长时间的数值模拟有很好的保结构性质.     

一维双曲系统的网格自适应算法

文章基于一类Godunov型的中心-迎风算法构造了根据解的局部性质调整网格疏密程度的自适应算法;提出了一个新的光滑探测器用以检测光滑区域和间断区域;该探测器简单实用,数值实验证明光滑探测器能很好地捕获间断区域;同时,分析了自适应算法的误差以及提供了该算法相关步骤的理论证明,理论分析表明,对于阶数r≥2的算法,网格自适应算法能减小误差。     

一类求解浅水波方程的基本无振荡熵稳定格式

针对浅水波方程,提出了一类低耗散基本无振荡熵稳定格式.在Roe型熵稳定通量中添加熵守恒格式的熵数值黏性绝对值的量来抵消解在跨过激波时所产生的熵增,从而抑制伪振荡;并且,利用通量限制器函数构造出相应的高分辨率熵稳定格式.利用新格式模拟一维和二维经典问题,数值结果表明,该格式具有低耗散、高分辨率、基本无振荡性等特点,是求解浅水波方程较为理想的方法.     

一种计算激波的守恒差分格式

介绍了一种处理激波间断解的差分格式,此格式是基于在间断解的两侧使用不同计算公式的守恒差分格式;证明了该格式满足离散熵条件。     

Riemann问题基于双激波近似的激波捕捉格式

Riem ann问题迭代求解的计算量很大 ,为了提高效率 ,减小计算量 ,提出一种 Riemann问题的近似解法。把Riem ann问题中的膨胀波看成是“膨胀激波”,认为左右状态之间存在由两个激波围成的相似解区域。该方法的优点是保留了 Riem ann问题的非线性特征 ,且保证左右波之间的自相似解区域的熵不减少。该方法在一维及多维无粘可压缩流动求解中取得了较好的结果。     

低耗散TVD格式的构造

目的　为进一步提高双曲型方程数值解的精度,减少数值耗散,同时保证稳定性,处理有强间断的流动问题; 方法　从总变差减小（ＴＶＤ）定义出发,分析保证数值解ＴＶＤ的条件,找出满足研究目的的可能性和途径; 结果　首次提出和证明双曲型方程数值解ＴＶＤ的充分必要条件,并据此简化导出适用于时间显示格式的ＴＶＤ充分条件,并据此构造了一种低耗散ＴＶＤ 格式; 结论　进行的数值实验证明,高分辨格式可以进一步减小数值耗散,提高数值解精度;