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

通过数值求解可压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的计算结果。从结果中可清楚地看到柱体非对称涡形成的非定常过程。     

气动计算中一个简单有效的差分格式

本文采用算子附加修正方法,提出了一个气动计算中简单而有效的差分格式。此差分格式的特点是保持原显式格式的精度,计算过程有显式格式的简单,又有隐式格式的特点,即对于任意步长此格式都是稳定的。本文利用此格式求解了N-S方程,给出了二维激波与边界层干扰的分离流计算结果,并与MacCormack的隐式格式进行了比较。计算结果表明,本文提出的格式与MacCormack隐式格式相比,大大简化了计算过程,节省了计算机时。     

用空间推进法来解钝锥三维粘性绕流

本文将文献［１］中提出的方法推广应用于三维可压缩Ｎ－Ｓ方程，文中给出了钝锥超声速三维粘性绕流的结果。得到了背风面分离的流动图画，计算结果与实验相吻合。     

用CSCM格式求解激波反射问题

本文对Euler方程采用CSCM格式,数值求解了激波反射问题,与目前流行的交替方向隐式(ADI)格式以及总变差递减(TVD)格式的计算结果进行了比较,并作了简要分析讨论。     

可压混合层涡并放热效应的涡动力分析

 基于双元单步不可逆燃烧动力学简化模型 ,采用紧致差分算法直接求解二维可压 Navier-Stokes方程 ,研究了时间发展非预混扩散火焰平面反应混合层流动。拟序结构由数值线性稳定性分析方法设定的粘性可压最不稳定小扰动初场促生。随放热强度增大 ,得到了基频涡卷撕裂现象。涡动力分析表明 ,放热主要通过膨胀和斜压这两种强非线性过程抑制大尺度剪切涡的演化     

超音速钝体绕流数值计算

本文提出了关于超音速钝体绕流数值求解的一种简便方法——非定常直线法。该方法在一般直线法的基础上利用具有时间项的双曲型方程组,将原求解一阶拟线性混合型方程组的边值问题,化为求解拟线性双曲型方程组的初边值问题。本方法的特点是对初始值要求不高,可克服一般直线法所遇到的奇异性问题,同非定常差分法比较,本方法的算法简单,花费机时少,计算结果也较为满意。用本方法给出了物形有曲率间断的超音速钝体绕流数值计算结果。     

六阶精度的群速度直接控制紧致格式及其应用

差分数值解产生非物理振荡的直接原因是在于非均一的波群的群速度,因此本文利用直接群速度控制的方法重新构造具有六阶精度的紧致型差分格式,以达到改善激波数值解的目的。对所构造的格式的精度及数值解随控制参数变化的行为进行了分析。最后文中给出了一些典型算例,证明了文中所构造的格式具有方法简便,物理含义清楚,精度高,网格基架点小,和捕捉激波能力较强的优点。      

分离流动的特性分析

本文分析了三维分离流的特性,且从基本物理现象出发,以数值求解NavierStokes方程所给出的三维粘性分离的结果和三维分离流的实验结果为依据,给出了三维粘性绕流产生分离的必要条件和物面上分离线的定义。