迎风TVD格式在粘性流计算中的应用研究与改进

首先以高超声速钝楔绕汉热流计算和Ｍ∞＝２．０、激波角θ＝３２．５８°的平板激波附面层干扰分离流动计算,研究了迎风ＴＶＤ格式的粘性特性、发现即使在很小的物面网格下,原有格式计算驻点热流仍较实验和文献结果低,激波附面层干扰的分离区小于实验值,原因主要在于熵修正公式。针对这些不足提出了新的熵修正公式,用改进后的格式重新计算,得到了与实验符合较好的结果。最后用改进的格式对分离流动的壁面温度控制效应进行了研     

非线性Schrodinger方程的守恒数值格式

对非线性Ｓｃｈｒｏｄｉｎｇｅｒ方程提出了一种新的守恒差分格式,并证明了该格式的收性与稳定性。通过数值计算,对非线性Ｓｃｈｒｏｄｉｎｇｅｒ方程中非线性项的离散进行了讨论,获得如下结论在取适当的参数后,所提出孤差分格式工上好于作为该格式特例的文（７）中的格式。     

激光烧蚀瑞利泰勒不稳定性数值研究

 介绍了激光烧蚀流体不稳定性计算程序EUL3D,其计算结果与Takabe 公式、FAST2D程序和LASNEX程序的结果以及日本大阪大学激光烧蚀瑞利—泰勒(RT)不稳定性实验,都较好符合,发现了横向电子热传导烧蚀在长波长扰动的非线性瑞利—泰勒不稳定性演变中起重要作用。在合理近似下,得到了烧蚀RT不稳定性线性增长率的预热致稳公式,此公式除包含了烧蚀对流致稳和密度梯度致稳因素外,还包含了Atwood数变小致稳因素,因此与各种情况的二维计算值都很好符合。     

色散方程的高稳定性两层四点显格式的单点精细积分法

基于单点精细积分的思想,对色散方程Ｕｔ＝ａＵｘｘｘ构造了一类高稳定性的两层四点显式差分格式,其局部截断误差为Ｏ（τ＋ｈ）稳定性条件为│Ｒ│＝│ａτ／ｈ＾３│≤ｆ（β）,对任意正实数β为单调递增函数,它们不仅显著地改善了同类格式的稳定性条件│Ｒ│≤０．２５而且也优于众多三层多点（５点或５点以上）显格式的稳定性条件。     

新型上风格式(AUSM+)的性能分析及其应用

研究了新型上风格式AUSM+的分辨率、效率等性能,并用它与Roe、vanLeer上风格式数值模拟了前向台阶激波反射流动,通过对激波、膨胀波、接触间断及其间相互干扰的复杂波系的模拟对比,分析探讨了AUSM+格式的低数值耗散、间断高分辨率等特性。     

高超声速三维热化学非平衡电离流动计算研究

研究发展了一种计算高超声速三维热化学非平衡电离流动的数值方法。真实气体热力学状态是通过平动-转动温度和电子-振动温度来模拟。对零攻角轴对称钝锥体和有攻角钝锥体的高超声速电离流动,采用隐式有限体积法NND格式,数值求解NS方程。     

二维三温能量方程的九点差分格式及其迭代解法

给出了激光驱动内爆数值模拟中二维三温能量方程的九点差分格式及其迭代解法,并给出了九点差分格式与五点差分格式对比计算结果,对一维和二维物理模型进行了数值模拟,得到了令人满意的数值结果。     

7阶WENO-S格式的计算效率研究

WENO-S格式是一类适合于含间断问题数值模拟的加权本质无振荡格式.这类格式的光滑因子满足对单频波为常数,这使得其近似色散关系与线性基底格式一致,并且具有良好的小尺度波动模拟能力．计算效率是数值方法性能指标的一个重要方面.由于WENO-S格式的光滑因子在各子模板上的计算公式除下标不同外形式一致,在计算线性对流方程相邻数值通量时,部分光滑因子完全相同．为此提出一种消除WENO-S格式冗余光滑因子计算的方法.该方法要求一条网格线上用于重构或插值的量可以表示为一个序列.基于此要求分析其对于几种不同物理问题的可行性和使用方法.以7阶WENO-S格式为例介绍了格式性质和去除冗余光滑因子计算的方法.该方法中预先计算和存储一条网格线上的所有光滑因子,在网格点较多的情况下,光滑因子计算次数约为原7阶WENO-S格式的1/4.对一维对流问题、球面波传播问题、二维旋转问题、二维小扰动传播问题及一维和二维无黏流动问题进行了数值模拟.结果表明该格式对多种流动结构具有良好的捕捉能力,并且同时具有良好的计算效率,去除冗余计算后又降低了约20%的计算时间.     

Vlasov-Poisson方程半拉格朗日守恒算法

利用三阶迎风插值多项式结合限制子方法,构造Vlasov-Poisson方程的半拉格朗日守恒型格式,可保持Vlasov-Poisson方程解的正性.计算朗道阻尼,双束不稳定性等典型问题,并与样条插值方法、UGKS方法进行比较,模拟结果表明半拉格朗日守恒性格式在Vlasov-Poisson方程求解中具有较高分辨率.     

基于二阶格式的有限体积保正格式

从二阶线性格式出发,通过对法向通量进行重构,得到非线性两点通量,获得四面体网格上的单元中心型有限体积保正格式.该格式适用于求解间断和各向异性扩散系数问题.无需假设辅助未知量非负,避免了辅助未知量计算出负时"遇负置零"的人为处理方式;并且证明该格式在每个非线性Picard迭代步具有强保正性,即当源项和边界条件非负时,线性化格式的非平凡解是严格大于零的.数值算例验证该格式具有二阶收敛性且是保正的.