扫描有限元法求解二元喷管的三维可压粘性流场

发展了一种二元喷管流场计算的三维Ｎ－Ｓ方程扫描有限元方法,给出了Ｎ－Ｓ方程中二阶导数项新的离散方法,以及更完善、简洁的扫描有限元方程的表达式。开发了相应的计算程序,用所发展的方法对二维平面收扩喷管、两种大宽高比的二元喷管及其圆变方转接段的三维粘性流场进行了计算,计算结果表明本方法具有较宽广的适用范围,计算结果良好     

计算矢量喷管流场的有限元法

本文通过采用扫描有限元法及多步龙格－库塔时间推进格式求解非定常平均Ｎ－Ｓ方程来模拟二元矢量喷管的内部跨速流场,成功地计算了二维有激波的跨声速喷管内部流场及仰 仰角３０°的三维矢量喷管内部流场。     

二元喷管跨音速流场的数值计算

本采用有限体积法进行空间离散，Runge-Kutta法进行时间推进的办法，通过求解Euler方程来模拟二元收敛－扩散喷管跨音速流场。运用隐式残差平均技术和当地时间步长法加速计算收敛。为了消除气流参数的波动和激波前后的振荡，在方程中添加了自适应耗散项。在100个时间步以内便可获得稳态解，计算结果与试验数据和其他数值方法吻合良好。     

多重网格加速的LUSGS算法用于喷管流场数值模拟

采用有限体积法,二阶精度插值的Roe格式对三维流场N-S方程进行空间离散;采用LUSGS隐式时间推进法和多重网格技术加快计算收敛速度;采用独特的显式算法克服了两方程湍流模型刚性对流场计算收敛速度的影响,并给出了该方法的理论依据。通过轴对称喷管内流场、二元喷管内外流场的数值模拟验证了采用以上方法在计算效率上的提高。最后将RNG模型计算结果与SA模型计算结果及实验数据做对比分析。     

姿控发动机喷管内流场及推力计算研究

运用时间相关法数值求解定常Ｎ－Ｓ方程和采用有限体积方法离散控制方程,对姿控发动机喷管内流场及其产生的推力进行了计算,得到了较为精细的流场图象。计算结果表明,不同壁面温度条件可以显著改变喷管壁面附近的流场结构及其产生的推力。     

二元矢量喷管流场及性能数值模拟

用有限差分法对二元矢量喷管的内流场及推力性能进行数值模拟。流场计算采用ＭａｃＣｏｒ－ｍａｃｋ两步格式，运用两层代数紊流模型，对初场、网格进行处理，采用加快收敛的当地时间步长，使计算精度和时间都能较好地满足工程应用的要求，得到了与实验数据较一致的喷管内特性。     

矢量二元收扩喷管热射流复杂流场特征研究

提出了一种二元收－扩矢量喷管的概念,用热射流实验台对二元矢量喷管进行了过膨胀状态下的喷管模型热态实验。测量了喷管壁面静压和喷管出口总压分布,用红外热成像技术给出了二元喷管热射流场的清晰画面,并对流场特征进行了描述和分析。研究表明,在非设计状态,二元喷管管内有明显的流动分离现象,喷口射流总压严重畸变,在热射流流场中,射流流场呈现双势核结构。     

二元喷管无粘流的有限体积法分析

本介绍了一种计算收敛－扩散喷管内流特性的有限体积法，它将求解域分成许多有限体积单元，将欧拉方程在小体积单元上积分，从而得到半离散方程，然后采用Mac Cormack格式进行时间推进，求得定常解。计算结果与试验值吻合良好，是一种有效的计算具有喉道小曲率半径不连续的收敛－扩散喷管特性的方法。     

次流喷射控制推力矢量喷管的流场数值模拟

为了探讨次流喷射控制二维推力矢量喷管的性能,用时间推进求解Ｎ－Ｓ方程的方法数值模拟了喷管的流场,数值格式用中心有限体积格式,借助于当地时间步长和隐式残差光顺技术加速收敛,得到了在不同二次喷流压力下推力矢量喷管内流场的变化特征。计算结果与实验结果比较有满意的一致性。研究表明,应用次流喷射控制主流流动可以实现较大的推力矢量转折,但是,二次喷流必须具有足够的压力值。     

轴对称矢量喷管内流特性的数值模拟研究

利用矢通量分裂有限体积格式和B-L湍流模型,采用TTM网格生成和局部网格加密技术,求解三维N-S方程,对轴对称矢量喷管的内流流场进行数值模拟,并将计算结果与模型试验结果进行了分析比较。研究结果表明,在矢量状态下,计算结果和试验结果相比,壁面静压分布的相对误差不大于13%,气动矢量角的相对误差不大于10%,推力系数的相对误差不大于2%,所开发的计算方法和计算程序可用于轴对称矢量喷管的工程设计。      

强制偏流喷管的数值仿真

采用薄层Ｎ－Ｓ方程成功地求解了强制流喷管内的三维复杂流场，对通量项和压力项用Ｖａｎ Ｌｅｅｒ发展的矢通量分裂法进行了离散，对离散方程采用线超松弛方法迭代求解，在亚音速与超音速气流共同存在的出口边界分情况进行了处理。计算结果验证了强制偏喷管具有自动补偿外界压力变化的能力，与实验结果基本吻合。     

超音进气道流场的数值计算

用Jameson的有限体积法,对某型号飞机在不同飞行状态下的进气道流场进行了数值模拟。求解的控制方程为三维非定常N-S方程。在对流项的计算上用二阶中心差分有限体积法,粘性项的计算则采用作者在［1］中发表的一种有限体积框架下的新的离散化技术,时间方向用多步Runge-Kutta方法进行推进,以增大CFL条件数。湍流模型则采用Baldwin-Lomax代数模型。求解过程中还采用了当地时间步长,隐式残量平均、人工粘性等加速收敛技巧。最后给出了各种飞行状态下进气逼出口截面的平均总压恢复与总压畸变,得出了比较满意的结果。     

用N－S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型￣［１］。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

固体火箭发动机内流场分区耦合数值计算

采用分区耦合方法计算固体火箭轴对称燃烧室与喷管流场。对于低速的燃烧室流场，选用不可压流的Ｎ－Ｓ方程描述并用ＳＩＭＰＬＥＣ方法数值求解；对于高Ｒｅ数的喷管流场，则采用Ｅｕ－ｌｅｒ方程描述并用ＳＣＭ方法求解。计算时用燃烧室出流为喷管流场提供入口参数，同时用喷管流场压强分布反馈影响燃烧室流动状况。对耦合边界条件处理方法进行了探讨。对典型的侧壁加质燃烧室与喷管流场进行了计算，计算结果揭示了单独喷管流场计算难以反映的喷管收敛段近壁区的低速区域，与已有的燃烧室流场实验结果一致并反映了燃烧室与喷管流场之间的联系，较好地模拟了流动中的物理现象。     

旋转固体火箭发动机喷管热结构分析

以某固体火箭发动机喷管为例,研究其旋转情况下的流动、传热及热结构响应性能。利用雷诺应力湍流模型,结合增强型壁面函数,求解N-S方程。以喷管内壁面温度分布为边界条件,求解二维轴对称瞬态热传导方程。将瞬态温度场按时间步长加载,进行瞬态静力学分析,得到不同转速下的流场、温度场及应力场。结果表明,喷管在旋转情况下,喷管内壁面尤其是喉部及扩张段温度分布在转速为100 r/min及600 r/min左右时变化剧烈;旋转内流场与喷管结构的耦合作用加剧了喷管的传热,尤其是喉衬的烧蚀,特征点温度值随转速增大而升高;最大热应力位于喷管最外层尾端,整体热应力在转速低于100 r/min时得到释放,随着转速的不断增大,喷管整体最大热应力及扩张段特征点应力随之增大,而喉衬特征点由于旋转导致了温度梯度降低,其应力值随之减小。     

高速粘性内流的高分辨率高精度迎风型杂交格式

本文在有限体积离散和ＬＵ分解的基础上，构造出一个新的隐式迎风型杂交格式，并用于求解定常流动的稳态解。     

螺旋式发射管内火箭燃气冲击流场数值模拟

用改进的MUSCL格式解三维、可压缩平均雷诺纳维尔-斯托克斯方程组,湍流模型为Spalart-Allmaras代数模型,计算和分析了火箭发射管内燃气冲击流场.在建立差分算法时,将有限体积离散和黎曼解算器相结合,简化了计算工作.首先,以小喉截面超声速喷管流动问题开展数值实验,数值结果与实验数据吻合很好;其次,计算了火箭发射管内燃气流场,获得发射管内气流压强和温度等物理参数分布.     

Simulation of low-pressure water flow through the Laval nozzle taking into account evaporation on the basis of a two-temperature model of two-phase flow

A water flow through the Laval nozzle is studied using the model of the one-velocity two-temperature two-phase medium motion written on the basis of the Navier-Stokes equations. As an equation of state for the carrying phase, Tait??s equation is used. The vapor phase dynamics is described by the conservation equation for the concentration of bubbles having the critical radius, the Rayleigh-Lamb equation, and the equation of vapor pressure inside a bubble. The equations for the one-velocity two-temperature medium dynamics are solved by the Beam-Warming method. The explicit contraflow scheme is used to solve the Rayleigh-Lamb equation and the conservation equation for the bubble concentration. The solutions of the quasistationary type obtained for the plane Laval nozzle make it possible to estimate the local flow parameters and the rate of evaporation at with different values of counterpressure. We compare the solutions for the models of independent phase motion and those taking into account the interphase exchange of mass, impulse, and energy.     

新型斜波瓣超声速混合器设计研究及其数值计算

根据超声速两股流掺混的特点, 在波瓣喷管的基础上提出了一种新的斜波瓣型超声速混合器设计方法, 设计的原则是在较小流动损失的前提下尽可能增加混合程度.利用有限体积法对三维可压粘流N-S方程进行离散, 紊流模型为k-ε模型, 分别对欠膨胀、接近完全膨胀和过膨胀3种工况进行了计算研究.计算结果表明, 当内通道瓣腔顶端与外通道相交处成一定夹角, 射流产生鞍形激波和旋涡可以增加掺混效率, 混合器内两股超声速流的流态基本符合设计要求.在欠膨胀、接近完全膨胀的工况, 效果是令人满意的, 但对严重过膨胀的工况, 到出口截面的混合效果不理想.