二维Euler方程的非结构网格再生成自适应方法

给出了提高二维Ｅｕｌｅｒ方程定常解质量的非结构网格自适应方法和模拟结果。流场算法为结点型有限体积法和Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式，网格用阵面推进法生成，生成初始网格时使用均匀的直角坐标背景网格。     

高精度差分格式及多尺度流场特性的数值模拟

简要介绍了作者近年发展的高精度差分格式,给出了两类格式精度的分析方法,分析了粘性面差分逼近所引入的耗散量低亏效应。文中按群速度将格式分为快型,慢型和混合型三类。对每类格式分析了高频波的传播特性。     

不可压扰动方程高精度对称紧致差分数值解法及应用

从粘性不可压扰动方程一阶改型形式出发,对其实现了高精度对称紧致差分离散,就导出的扰动线性特征值问题给出了一个高铲双重迭代局部解法,以相同精度将特征值和特征函数同时得到。通过不可平面Ｐｏｉｓｅｕｉｌｌｅ流时间稳定性算例详细对比显示了算法良好的谱分辨能力和较弱的网格依赖性,并结合复矩阵广义特征隐式移ＱＺ算法获取了一个 扰动特征值谱计算结果。     

捕捉激波的群速度控制方法

回顾及综述了群速度控制方法基本原理、发展过程及最新进展。群速度控制方法是通过控制色散误差以捕捉间断的一种数值方法。该方法通过在间断两侧采用不同色散类型的格式,将数值波的传播压制在间断附近,避免数值振荡向周围扩散,从而达到抑制数值振荡的效果。该方法主体上采用数值色散抑制振荡,实际使用中仅需补充较低水平的数值耗散就可将振荡压制在可接受的范围内。因而群速度控制格式的总体耗散水平较低,易于实现小尺度波的高分辨率捕捉,是复杂流动高分辨率数值模拟的有效方法。     

紧致格式数值模拟超音速粘性绕流问题

通过数值求解可压Navier-Stokes方程的方法计算了球锥粘性绕流问题。方程中粘性项按通常办法逼近,为改善收敛速度利用了文[1,2]中之算子附加修正方法。方程中之流向导数利用了四阶精度的紧致差分。这是由于该方向上流动参量光滑和边界条件易于处理的缘故。在法向方向上利用了中心差分。对差分方程之隐式部分利用了近似因式分解法。由于选用了文[3]中之特殊Jacobiall系数矩阵分裂法而使计算工作量大为减少。用这种方法计算了超音速粘性绕流问题。计算结果与实验结果相符很好。     

高超音速平板前缘干扰数值解

本文采用了文献〔1〕中提出的一个数值求解Navier-Stokes方程的空间二阶精度单步差分格式,数值计算了超音速和高超音速平板前缘干扰问题。通过简单的模型方程对格式进行了描述和分析。用这一方法计算了M_∞=20,Tw/T_0=0.06,Re_(∞L)=4×10~3和M_∞=3,Tw=T_0,Re_(∞L)=10~3的平板前缘干扰问题。计算结果与实验值比较符合的很好,通过这些计算结果可使我们更清楚地了解平板前缘干扰的流动特性。     

钝体绕流流场分离涡的数值模拟

 采用近似因式分解和特殊Jacobin系数矩阵分裂法,求解可压缩Navier-Stokes方程,数值模拟矩形柱体绕流流场中分离涡的形成和发展。给出了Ma_∞=0.3,0.8;Re=10~4,10~6的计算结果。从结果中可清楚地看到柱体非对称涡形成的非定常过程。