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

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

时间发展可压化学反应混合层直接数值模拟

基于一个简化的单步不可逆燃烧动力学模型,对时间发展二维可压化学反应平面混合流场,实现了高精度迎风紧致有限差分直接数值模拟。研究了化学反应放热与流体混合运动间的耦合作用,分析了前者对可压混合层厚度、拟序结构及产物生成率等方面的影响。     

超声速粘性剪切流空间稳定性对称紧致差分数值分析

用一类高精度对称紧致差分格式数值离散一阶改型三维可压粘性扰动方程,对导出的非线性离散特征值问题采用二阶修正Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ边值迭代局部解法,实现了超声速剪切流的线性空间稳定性分析。基于对流Ｍ数、Ｒｅ数、速度比、密度比等参数研究,讨论了压缩性效应、粘性效应、超声速扰动快／慢模态等,结果显示超声速粘性剪切流的弱不稳定性同多种影响因素密切相关。     

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

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

三元Navier-Stokes方程的数值计算——激波与附面层的相互干扰

前言 对于空气动力学所提出的问题,通常所需要的只是其定常解。稳定法被广泛地用来解决这类问题。显式格式更为简单,但稳定性要求很小的时间步长,因而计算机时间题是一个突出问题。为克服这一问题,近年来人们作了很多工作。本文采用了一个单步差分格式来计算三元激波与附面层干扰问题。此方法简单,在解的光滑区,格式具有形式上的二阶精度,并且格式是高度稳定的。计算中另一个问题是,在激波附近的     

一个三元分离流动的数值分析

1.前言 弹体飞行过程中会遇到各种分离流动,数值求解Navier-Stokes方程的办法是研究分离流动的重要手段,二元问题的数值模拟已有很大进展,但解决三元流动仍有很大困难,除流动本身的复杂性外,计算机的容量和速度也给以很大限制,计算方法有待于不断改进。 在[1]中利用一单步差分格式数值计算了一简单的三元分离流动,这里将对结果进行分析。本文第一部分对物理模型做一简述,第二部分分析计算结果并讨论三元分离流动特性。     

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

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

空间发展混合流的大涡模拟

采用大涡模拟方法数值模拟了空间发展混合流，通过在平均流中引入小扰动，研究了混合流大涡结构的产生和演化过程，捕捉了展向涡的卷起、配对、合并，以及二次流向涡的出现等大尺度的三维拟序结构。与实验及直接数值模拟的结果相比，吻合程度较好，表明本文所采用的LES方法对湍流的数值模拟是切实可行的。     

后向台阶横向喷流干扰流场的数值研究

 用有限差分法求解层流全N-S方程,对超音速后向台阶绕流与台阶后横向喷流干扰流场进行了数值模拟。不带喷流的超音速后向台阶绕流数值模拟结果与实验及其他计算结果吻合较好。最后给出横向喷流压力比p_j/p_∞=3.033～7.033和喷口宽度w=0.1h～0.2h(h为台阶高度)时的数值模拟结果。得到了合理的流场结构;并可看到喷流压力比和喷口宽度对流场特性有显著影响,在计算参数范围内,底部压力与系数Q(=wp_j/hp_∞)几乎成线性关系。     

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

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