粘性流场三维NS方程有限体积法数值解

提出了有限体积法中的隐式ＮＮＤ格式数值求解三维全ＮＳ方程的薄层三维ＮＳ方程的方法，对绕带座舱的再入飞行器高超声速粘性流动进行数值模拟。     

超声速二维湍流分离流动的数值解法

本文采用层流中引导简化N-S方程的方法,建立了简化的湍流运动方程组,分别用零方程、一方程模型计算了绕二维压缩拐角的湍流分离流动。数值求解用显隐式的反扩散格式。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

绕任意刚性翼型,机翼非定常流动的跨声速欧拉方程解

本文在有限体积方法的基础上引入变系数隐式残值光顺及固接网格法求解了绕二维翼型及三维机翼非定常流动的跨声速欧拉方程,计算结果与同类结果及实验结果吻合。     

用三维N—S方程计算钝体绕流

本文采用有限体积方法求解三维低速Ｎ－Ｓ方程，研究了三维钝体和类汽车物体的绕流。用中心差分和迎风差分的混合差分格式构造代数方程组，采用了隐式和显式相结合的方法使主对角元素占优以提高迭代求解过程的稳定性。利用压力校正算法，采用面松驰迭代来实现方程的求解，湍流的计算采用了标准的双方程模型，本文给出了钝体和类小汽车外形的计算结果，钝体的计算结果和实验结果进行了对比，与其它方法的计算结果也进行了对比，它们之     

用三维N－S方程计算钝体绕流

本文采用有限体积方法求解三维低速Ｎ－Ｓ方程，研究了三维钝体和类汽车物体的绕流。用中心差分和迎风差分的混合差分格式构造代数方程组，采用了隐式和显式相结合的方法使主对角元素占优以提高迭代求解过程的稳定性。利用压力校正算法，采用面松弛迭代来实现方程的求解。湍流的计算采用了标准的双方程模型。本文给出了钝体和类小汽车外形的计算结果。钝体的计算结果和实验结果进行了对比，与其它方法的计算结果也进行了对比，它们之间具有合理的一致性。     

N-S方程数值模拟弹体不可压脱体涡绕流

应用在连续方程中增加压力对时间导数项的拟压缩性方法 ,数值求解三维定常不可压N- S方程。采用隐式 Beam- Warming近似因式分解格式求解。以修正的 Baldwin- Lomax代数模型模拟湍流流动。对绕正切拱形旋成柱体的大迎角脱体涡流动进行了数值模拟 ,并分析了层流粘性和湍流粘性对脱体涡及二次分离涡的影响。在头部施加微小的非对称射流扰动可得到与实验相符的非对称脱体涡     

有限体积格子Boltzmann方法用于近空间连续流区绕流模拟

针对航天器再入解体形成残骸碎片的近空间绕流计算问题,拓展了格子Boltzmann方法在可压缩流动模拟的能力。引入有限体积隐式格式求解耦合双分布函数格子Boltzmann模型方程及采用圆函数为基础构造的D2Q13离散速度模型;引入IMEX-RK格式进行了时间项离散解决模型方程的源项刚性问题;对Riemann问题、平板双马赫反射问题、RAE2822翼型跨声速绕流等近空间连续流区、可压缩典型案例进行了数值模拟。通过比较分析,初步验证了耦合双分布函数有限体积格子Boltzmann方法对连续流区、可压缩流动的模拟能力;进一步开展了方柱形解体残骸的超声速绕流模拟,得到与N-S方程计算结果吻合一致的绕流结果,证实经改进的耦合双分布函数有限体积格子Boltzmann方法对解体残骸碎片绕流问题具有较好的模拟能力。     

拟压缩性方法数值模拟绕弹体的不可压缩层流脱体涡流动

本文应用在连续生方程中增加压力对时间导数项的拟压缩性方法，数值求解三维定常不可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程。数值求解过程采用隐式的Ｂｅａｍ－Ｗａｒｍｉｎｇ近似因式分解格式。计算绕尖头旋成柱体的大迎角脱体涡流动，得到与实验相符的结果。     

不同马赫数的无粘和粘性流动高阶精度隐式计算方法

 应用隐式时间推进法对不同马赫数的无粘和粘性流动进行数值分析,给出了基于预处理方法的高阶精度隐式求解方法。方程离散采用改进的高阶精度对流迎风分裂格式。该方法通过求解曲线坐标系可压缩Euler和 Navier-Stokes方程,对低速到超音速范围内的无粘和粘性流动的典型问题进行了数值分析。计算结果与文献的数值结果或实验数据对比分析表明,该方法对低速到超音速范围内的无粘和粘性流动问题进行数值分析是可行而有效的。     

大展弦比飞机翼身Euler绕流计算

给出了用Ｅｕｌｅｒ方程对大展弦比翼身组合体绕流的数值模拟，将Ｋｒｏｌｌ的求解欧拉方程的有限体积法，推广到求解机翼机身流动问题，基本思想是使用保角变换、代数变换和椭圆方程变换混合方法生成机翼机身贴体网格。     

基于粘性非结构网格的三维增升装置N-S方程数值模拟

应用基于粘性非结构网格求解N—S方程方法数值模拟了某型飞机三维增升装置的绕流。在物面附近粘性主导作用区域内，采用三棱柱型非结构网格，在其他流场区域，采用传统四面体非结构网格的混合非结构网格生成技术。湍流模型采用Spalart—Allmaras一方程模型，运用格心格式有限体积法，对控制方程和湍流模型方程进行数值离散。用四步龙格-库塔方法作显式时间推进，采用当地时间步长、隐式残值光顺等加速收敛措施．求解某型飞机全机带增升装置构型的粘性绕流流场，其数值计算与风洞试验结果显示了良好的吻合度。     

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

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

湍流附面层方程的有限差分解

本文用三点隐式有限差分法数值求解了轴对称流动的层流、转捩和湍流可压缩附面层方程,湍流附面层用两层概念处理,而每一层用适当的涡粘性模型代替Reynolds应力项,用常值湍流Prandtl数将湍流热流项与Reynolds应力关联起来。支配方程用迭代法求解。计算结果表明,采用本文改进的涡粘性模型和涡粘性分布的平滑办法,对文中所计算的几种情况,都给出了满意的结果。     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

二维粘性不可压缩流场的分块耦合求解方法

本文用分块耦合的方法，并采用S.Abdallah^[1,2]和M.L.Mansour^[3]的在非交错网络上隐式求解原始变量的Navier-Stokes方程的提法和离散格式，成功地解决了二维复杂流场的数值求解问题，获得了空腔流、T型叉管流、NACA0020翼型绕流和方块钝体绕流的稳态解，并进行了讨论。此方法还可推广到三维复杂流场的计算问题。     

三维层流分离的N—S方程隐式上风格式数值模拟

数值模拟了绕半圆球柱M∞＝1．2，a＝19°的粘性分离流和涡流．用隐式LU－SGS及Roe的通量差分分裂格式数值求解了薄层N－S方程，分析计算了所得的流场结构，包括物体背风面的涡流及三维分离方式，计算结果同实验结果定性及定量吻合很好，另外，还用流动拓扑理论分析了计算所得的流场形态。     

二维带动力吸气式高超声速飞行器绕流的PNS-NS混合求解

为了提高吸气式高超声速飞行器绕流的求解效率,采用空间推进求解抛物化PNS(parabolizedNavier-Stokes)方程和时间迭代求解Navier-Stokes(N-S)方程的混合计算流体动力学(CFD)方法来求解高超声速飞行器整机绕流.在超声速占主导的流动区域采用空间推进求解抛物化N-S方程的方法,在亚声速和分离区采用时间迭代求解N-S方程的方法.对于求解二维带化学反应的吸气式高超声速飞行器绕流,混合CFD方法和完全时间迭代方法相比,可得到同等准确的数值模拟结果,并且求解效率提高了数倍.      

高速列车绕流的数值模拟

给出了用低速Ｎ－Ｓ方程求解高速列车绕流的数值计算方法，采用显隐式结合的方法来加快收敛速度，基于该方法编制了计算软件系统Ｔｒａｉｎ ３Ｄ，用该软件系统对高速列车外形的外流场进行了数值模拟，获得了该外形的空气动力特性和流场结构。     

跨音速喷管化学反应流数值模拟

用时间相关的半隐格式有限差分数值方法求解了化学非平衡反应跨音速喷管流场，在喷管收敛段，流动接近化学平衡状态，控制方程的刚性问题严重，数值积分困难。通过对时间差分项隐式离散、对空间差分项显式离散，流场边界采用参考平面上的特征线法计算等，成功地解决了由于化学反应有限速率带来的数值解不稳定问题。该格式简单、需要计算机存贮空间少。本文完成了一维和轴对称非平衡化学反应喷管流动计算，并与化学平衡流和冻结流的计算结果做了比较。