解析方法与数值模拟相结合求解炸药在容器内部爆炸产生的壁面载荷

本文采用总变差递减（TVD）数值方法与强爆炸解析解相结合的手段,求解了炸药在爆炸容器内部产生的作用载荷。计算结果与实验值在几个典型位置进行了比较,二者有良好的符合。与单纯的数值方法相比,本文所探索的方法能有效地克服对称轴附近的数值振荡,不失为确定爆炸载荷的一种新途径。     

用高阶高精度 WENO 格式求解二维激波附面层干扰流场

采用时间相关法求解二维 Navier-Stokes方程 ,数值模拟二维平板层流附面层与激波干扰流场 ,给出了物面压力分布和应力分布。计算中 ,对流项空间导数的差分离散采用高阶高精度 WENO格式 ,时间方向采用具有 TVD性质的 Runge-Kutta方法 ,粘性项采用二阶中心差分。所得压力分布和应力分布与国外实验结果吻合较好 ,计算实践表明高阶 WENO格式具有优异的性能 ,应用前景广阔     

用多重网格TVD方法求解三维高超音速粘性流场

本文采用时间相关法求解了三维Ｎ－Ｓ方程。在计算中，对于压力和对流项采用Ｒｏｅ的通量差分分裂技术和Ｏｓｈｅｒ－Ｃｈａｋｒａｖａｒｔｈｙ的ＴＶＤ格式，对于粘性项，采用了中心差分格式。通过数值模拟得到了钝头双锥体高超音速层流流场在不同攻角下的数值解。计算的重点放在壁面热流率的预测。在攻角为４°和１２°的状态下，Ｓｔａｎｔｏｎ数的计算值与实验值进行了比较，二者吻合较好。为了加快收敛速率，计算时采用了多重网格方法。     

高超音速流场的隐式TVD格式并行计算

本文采用区域分割技术和拼接网格的并行策略，发展了一个适合于分布式存贮多机系统的ＴＶＤ隐式有限体积并行算法；并在ＰＶＭ并行环境下，对三维高超音速绕流流场实现了多机并行计算，通过负载平衡等方法得到了较高的加速比（在二处理机系统上加速比为１∶８４，在四处理机系统上为３∶４４）。     

基于TVD思想的高阶迎风紧致格式

紧致格式具有模板紧凑、精度高的特点,然而直接应用该格式捕捉激波会产生虚假振荡。解决该问题有多种思路,本文着重研究基于TVD思想的紧致方法。首先推导了两种5阶迎风紧致格式,并采用Fourier分析比较其耗散与色散特性。然后,采用保单调保精度方法和限制器两种不同的TVD方法,进行算例验证。计算结果表明,应用TVD方法后,格式对激波、接触间断分辨率好且鲁棒性增强,但不同的TVD方法仍存在精度降低、耗散不足等问题。     

高超音速粘性流场数值计算中差分格式的对比分析

本文通过三维高超音速粘性流场的求解对三种常用的差分格式： Ｈａｒｔｅｎ 的 Ｔ Ｖ Ｄ 格式、修正后的 Ｏｓｈｅｒ Ｃｈａｋｒａｖａｒｔｈｙ 的 Ｔ Ｖ Ｄ 格式以及 Ｊａｎｅｓｏｎ 的当地极值递减 （ Ｌ Ｅ Ｄ） 格式对激波的分辨率、边界层计算的准确性和格式的计算效率等三个方面进行了较细致的对比分析， 为三维高超音速复杂流场的数值模拟打下基础。