带配平翼钝体高超声速粘性绕流的数值模拟

采用ＮＮＤ差分格式，通过求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，数值模拟了带配平翼钝体的高超声速粘性绕流。文中首先以球头驻点压强和热流为参数，讨论了在差分方程左端采用ＬＵ－ＳＧＳ、ＬＵ－ＡＤＩ和对角化ＡＤＩ三种不同隐式算法收敛效率的异同。然后以带配平翼钝体的高超声速粘性绕流为模型，对所研制的程序进行了计算验证。在计算中采用了代数方法和求解椭圆型方程方法相结合的网格生成技术，针对配平翼外形给出了贴体性、     

化学非平衡高超声速电离流动数值研究

本文对非平衡化学区应高超声速绕流零攻角轴对称钝头体电离空气粘性流动进行了研究，通过带有化学源项的Ｎａｖｉｃｒ－Ｓｔｏｋｅｓ方程对流动作数值模拟，总的连续方程由单个组元守恒方程代替。高温电离空气中存在７个主要的组元，它们依次为Ｎ＿２，Ｏ＿２，ＮＯ，ＮＯ￣＋，Ｎ，Ｏ，和ｅ￣－。化学动力学模型由７个组元和这些组元之间６个可能的反应来描述。发展了一种数值方法以求解这组完全耦合的与时间有关的偏微分方程组。利用一类新的有限差分方法来求解这组以守恒形式表示的控制方程。这是一种以ＴＶＤ格式来离散粘性流通量的方法。在数值实验中采用了ＴＶＤＹｅｅ－Ｈａｒｔｅｎ格式，对全部算例都是用激波捕捉法。关于化学组元，我们试图处理完全催化壁、非催化壁和有限催化壁等三种不同的表面边界层条件。计算是对ＲＡＭＣ－Ⅱ试验条件进行的，得出了７１ｋｍ高空情形的典型流场数值解。对表面压力，Ｓｔａｎｔｏｎ数和热流率给出了化学反应气体结果与完全气体结果的比较。计算的流场电子密度与ＲＡＭＣ－Ⅱ实验数据进行了比较，计算结果与实验数据符合是相当好的。     

Navier—Stokes方程粘性项的有限体积离散法研究

Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程粘性项的有限体积离散法研究封建湖西北工业大学应用数学系［摘要］：三维Ｎ－Ｓ方程是描述真实气体流动的非线性偏微分方程组，如何快速准确地计算其粘性项至关重要。这不仅因为粘性项的计算量非常巨大，而且粘性项对急剧变化的流动梯度的...     

激波运动与反射流场的ENO算法

在ＴＶＤ格式的基础上，以修正限制函数的方法，将通常的ＴＶＤ格式改造为性能更为优越的ＥＮＯ格式，并利用这种有限差分格式求解欧拉方程和薄层Ｎ－Ｓ方程。分别数值模拟了激波管中运动激波与三种楔形物体相互作用流场；给出了激波在不同物体上的反射及波涡干扰过程中非定常流场的变化情况；算出的波系结构在与ＣＡＲＤＣ小型激波管中做的相关实验结果一致。     

求解Euler／Navier—Stokes方程的有限体积龙格库塔方法

提出了一个求解定常Ｅｕｌｅｒ／Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的计算方法,求解的主管方程是积分形式的Ｅｕｌｅｒ方程或雷诺平均的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,计算方法是有限体积多步龙格库塔方法,通过在格式中加入自适应的二阶加四阶人工粘性项以保证数值计算的稳定性,用空间变时间步长、焓阻尼来加速定常问题的收敛速度。     

旋转矩形通道中粘性流动的数值模拟

通过有限差分方法（ＭＡＣ－ＡＦ１）求解时间相关的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，对旋转矩形通道中粘性流动进行了数值研究，旋转轴方向与主流方向垂直，给出了充分发展时矩形通道中的二次流结构及主流速度分布，得到的数值解与实验吻合良好。     

静止叶栅三维粘性数值模拟

应用Ｎ－ Ｓ方程压力修正算法和ｋ－ ε两方程湍流模型,对端弯叶片和基准（无端弯）叶片的平面叶栅进行了三维粘性数值模拟。在采用简单交错网格的情况下,通过对计算控制体通量的速度进行修正的方法解决了Ｎ－ Ｓ方程压力修法算中遇到的压力场和速度场呈锯齿状分布的问题。用实验测量值对本数值模拟方法的精度和可靠性进行了检验     

粘性不可压流动压力修正的有限元计算

将ＳＩＭＰＬＥ方法压力修正思想应用于有限元方法，对粘性不可压缩流动的ＳＵＰＧ有限元方程，建立将压力与速度分离，进行压力修正迭代的计算格式；对典型算例进行多个Ｒｅ数的层流，湍流计算，显示此方法具有良好的收敛性和精度。     

飞船高超声速粘性绕流的数值模拟

本文采用稳式ＮＮＤ格式，通过求解三维薄层近似Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，数值模拟了飞船高超声速粘性流场。在计算中充分考虑了网格雷诺数和网格分布对计算结果的影响，不但得到了满意的表面压力、热流分布和流场的物理量分布，而且得到与实验观察一致的表面油流图像和分离结构。     

跨声速机翼流场的N—S方程计算

本文给出了三维机翼粘性流场的数值模拟过程，控制方程为时均化的三维可压缩薄层Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，湍流模型采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数湍流模型，空间离散采用中心有限体积格式，时间采用Ｒｕｎｇｅ－Ｋｕｔｔａ多齿格式进行式积分。     

微波电热推力器钝头体方案的性能计算

进行了微波电热推力器（ＭＥＴ）钝头体方案的性能计算,对微波作了谐振模式ＴＭ０１１和波导模式ＴＭ０１两种选择。将计入电磁场作用的轴对称钝头体粘性流动方程,结合ＴＭ０１１或ＴＭ０１模式,与等离体区耦合求解,获得了ＴＥＴ比冲和推力等性能,结果与文献的实验及计算结果趋向一致。     

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。     

径向开槽药柱发动机内流场数值模拟

采用时间相关的显示ＭａｃＣｏｒｍａｃｋ预估－校正两步格式和可压缩全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，计算了轴对称径向开槽药柱发动机内流场。计算中采用了贴体坐标系统和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型。边界参数计算采用参考平面上的特征线法和物理边界条件。通过大量算例的分析表明，方法有效，计算结果合理。当有集中质量加入时，燃烧室头部压强升高，喷射强度和开槽口的倾斜方向对头部压强也有较大的影响。     

跨声速和超声速流中激波/边界层干扰数值模拟

对SST湍流模型中的Bradshaw常数a1进行了修正,并对跨声速和超声速流中激波/边界层干扰进行了数值模拟研究,空间离散采用二阶精度差值的低耗散通量分裂格式(LDFSS),时间离散采用对称高斯-赛德尔(SGS)算法。结果表明:在跨声速流动中,计算得到的壁面压力分布、分离区长度和速度剖面都与实验值吻合较好,而且很好地模拟了典型的λ激波结构;在超声速流动中,修正后模型的计算精度较原始模型有了较大改善,计算得到壁面压力分布和分离点的位置都和实验值吻合较好。     

高马赫数下激波湍流边界层干扰数值研究

应用GAO—YONG可压缩湍流方程组数值模拟了入射斜激波／平板湍流边界层相互干扰现象,计算了来流马赫数为5．0,激波入射角度分别为15．876°、23．287°两种不同激波干扰强度下的流场。计算程序中的对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge—Kutta显式时间推进法求解空间离散后的控制方程。计算较好地模拟了高马赫数下的激波／湍流边界层干扰的流场结构,位移边界层厚度,动量损失厚度等,也比较准确地预测了平板壁面压力、摩阻系数等气动力参数的分布。     

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

微可压缩模型(SCM)与可压缩N-S方程数值计算对比研究

本文开展了微可压缩模型与可压缩N-S方程之间的对比研究.通过对层流平板边界层和椭球高雷诺数湍流绕流的数值模拟,比较了微可压缩模型和可压缩N-S方程之间的计算精度和效率.计算结果表明,微可压缩模型和可压缩N-S方程具有相同的计算精度,但是微可压缩模型具有更高的计算效率,内存和CPU时间都节省20%左右.     

Gao-Yong理性湍流方程

对近十年来Gao-Yong理性湍流方程的研究进行综述,给出矢量形式的可压湍流封闭方程组的推导过程。侧偏统计平均保留了湍流脉动量的一阶统计信息,在引入加权漂移速度对称性及正交各向异性后,导出了漂移流的连续方程、动量方程及机械能方程,最后依据湍流物理的唯象论,使用动量传输链概念模化封闭了整个方程组。方程组不含任何经验系数,不使用壁面函数,保留了NS方程的均化的非线性特性。其级数形式的能量方程与非线性现象多尺度层次现象相对应,具备了描述湍流统计平均流动及拟序结构流动的双重功能。已进行了大量算例验证,验证结果证明了Gao-Yong湍流方程对广泛范围的复杂湍流问题的适定性。     

不同轴长比椭球体的绕流计算与分离形态分析

用拟压缩性方法求解不可压Ｎ－Ｓ方程,对４∶１旋成体椭球和其他轴长比椭球的层流绕流进行数值模拟。分别用二阶中心差分格式和三阶迎风紧致差分格式来逼近无粘项。从计算所得的表面摩擦应力线、空间流线和横截面流场特性,分析了１０°≤α≤７０°迎角范围内的分离和脱体涡系形态,算出了相应的气动特性。计算结果表明,三阶迎风紧致差分比二阶中心差分可更好地模拟分离涡结构。迎角直到４０°时,算出的流态仍然对称,而迎角为７０°时,算出的流态为非对称非定常,与实验流态基本相符。椭球体横截面形状越扁平,开式分离就越提前产生。椭球体横截面形状横向越窄,分离区就越小。由于改变椭球体横截面形状具有改变分离和脱体涡系的效果,本文的研究结果有助于改善三维物体的气动力特性,为改进机体设计提供理论依据     

MHD控制激波诱导湍流边界层分离的机理分析

为了研究磁流体动力学(Magnetohydrodynamics:MHD)加速边界层对激波-湍流边界层相互作用的影响,用高阶有限差分法求解了小磁雷诺数近似的MHD湍流方程。其中,无粘通量采用WENN格式离散、粘性通量采用Roe平均中心差分离散,时间采用半隐式推进,并采取追赶法求解。计算给出了湍流、电场、磁场和电导率等参数对边界层分离的影响,数值结果显示:在同样的逆压梯度下,湍流边界层分离能更快地趋于稳态流场,且分离区比层流小;通过施加洛仑兹力加速,边界层速度型面变得更加饱满、位移厚度减小、分离点和再附点向激波与固壁的交点靠近,分离区尺寸减小甚至最终被消除。