基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

铣槽结构液体火箭发动机推力室壳体热应力分析

介绍了液体火箭发动机推力室铣槽结构热应力的数值分析方法,通过建立液体火箭发动机推力室的流场燃烧和导热理论模型,运用有限体积法考虑液膜冷却计算出发动机工作时的燃气、燃烧室壳体和冷却工质的温度场,将得出的结果作为壳体热应力计算模型的边界条件进行热应力场有限元分析。内、外壁温度的计算数据与实验结果基本相符。     

Investigation on the Thermal Conductivity of 3-Dimensional and 4-Directional Braided Composites

It is vital to choose a factual and reasonable micro-structural model of braided composites for improving the calculating precision of thermal property of 3-D braided composites by finite element method (FEM). On the basis of new microstructure model of braided composites proposed recently, the model of FEM calculation for thermal conductivity of 3-dimennsional and 4-directional braided composites is set up in this paper. The curves of coefficient of effective thermal conductivity versus fiber volume ratio and interior braiding angle are obtained. Furthermore, comparing the results of FEM with the available experimental data,the reasonability and veracity of calculation are confirmed at the same time.     

用非结构直角网格和欧拉方程计算飞机绕流

 采用八叉树结构 ,生成复杂外形绕流计算的非结构直角网格。物面附近用投影方法 ,使网格贴体。并将Jameson的有限体积法推广用于这种网格的欧拉方程计算。对歼击机模型的绕流计算表明 ,网格生成的机时花费很少 ,总体质量好 ,因而欧拉方程解算的收敛质量也好。     

细观力学有限元法预测复合材料宏观有效弹性模量

基于能量等效原理提出了复合材料有效弹性模量的定义，并指出了该定义的基础及前提条件。为从理论上计算复合材料宏观有效弹性模量，建立了通过细观力学有限元法计算复合材料有效弹性模量的方法。复合材料宏观弹性模量，是通过对复合材料细观结构代表性体积元的力学响应的计算来得到，在该计算方法中，给出了施加简便的边界载荷以及恰当的边界变形约束条件的方法。数值计算结果与部分试验结果具有较好的一致性，表明所提出的方法能够较好地计算复合材料的宏观有效弹性模量。     

分区和有限容积法计算长尾喷管中的湍流流动

为了给长尾喷管热防护结构设计提供热边界条件，从轴对称二维雷诺平均N－S方程出发，采用分区算法和有限容积法计算了长尾喷管中的湍流流动，分析了流场的结构特性，提出了长尾喷管中燃气温度，压强和马赫数分布的某些规律以及燃气温度沿壁面的分布，计算和分析结果提出了长尾喷管热防护结构设计应注意的事项。     

基于Euler／N—S方程的跨音速非线性静气动弹性问题研究

在C-H网格的基础上,采用Jameson的中心差分有限体积法求解Euler/N-S方程,采用结构影响系数法计算结构的弹性变形,用三角元面积加权法和常体积转换法(CVT)实现流固耦合,计算了弹性后掠机翼在跨音速状态的静气动弹性问题     

求解Euler方程的无网格与有限体积混合算法研究

研究了无网格方法与有限体积法相结合的、求解Euler方程的混合算法.与单一的无网格方法相比较,由于在大部分计算区域采用了有限体积法,使得发展的算法在运算效率上能与有限体积法相当;同时在物体附近嵌入了无网格区,使得在处理几何外形上更具灵活性.有关算法涉及的无网格与网格区域搭接处理,通过在交界面处引入辅助卫星点和网格单元构成边界信息,实现了区域间的流动信息传递.Euler方程的空间导数分别在两区域用有限体积法和无网格法离散近似,时间方向都采用四步显式Runge-Kutta格式推进求解,数值模拟了喷管内流和绕翼型外流,并分别与整体有限体积法和整体无网格方法进行了比较.算例展示出,提出的混合算法能有效捕捉激波间断,且两区域等值线过渡光滑,算法效率如预期与有限体积法相当,表明该方法是可行的.     

替代燃料贫油熄火边界影响因素

为了研究燃油以及入口空气压力对于贫油熄火(LBO)边界的影响大小及规律,采用航空煤油(RP-3)、高沸点费托油(FT)和柴油进行了三种不同燃烧室入口压力工况下的贫油熄火实验并进行规律分析。分析结果表明:入口压力对贫油熄火边界的影响(19.17%)要大于燃油性质造成的影响(6.26%)。导致熄火油气比变化的主要因素包括入口空气压力,火焰体积,燃烧室温度,燃油雾化直径以及燃油的热值和密度,其中火焰体积和燃油雾化直径主要受燃油性质影响,而燃烧室温度则与入口压力有很大关系。入口压力影响的贫油熄火油气比变化会受其影响的火焰体积和燃烧室温度变化而削弱。碳数高支链烷烃含量少的燃油可能会提高火焰体积对熄火油气比的影响,使其在低入口压力下有更好的贫油熄火边界。      

化学非平衡流动解耦算法的回顾与新进展

超声速化学非平衡流动的数值模拟一直是计算流体力学领域的难点,主要体现在如下几个方面：物理过程非常复杂,存在着激波、燃烧波等各种复杂波系的相互作用;超声速化学非平衡流动属于典型的时空多尺度物理问题,其控制方程存在严重的刚性,给数值求解带来了很大困难。对国内外的解耦算法研究现状简单综述后,主要介绍1993年刘君提出的解耦算法的理论基础,流动方程采用冻结流模型,源项方程模拟流体微团在当地绝热、定容的热力学系统内发生的化学反应过程。通过引入两个中间变量,即等效内能和等效比热比,将与温度无关的生成焓从流动方程组能量项中分离出去,源项方程组中包含等效内能,使用不同算子对流动方程和源项方程解耦求解。与传统解耦算法相比,源项方程的求解过程中包含状态参数和组元同时变化。结合刘君解耦算法机理和有限体积法空间平均特性,介绍近期在提高算法计算效率方面的研究进展,包括流动方程优化算法和耦合过程优化。采用优化算法对经典的激波诱导燃烧算例进行数值模拟,与不同文献结果进行对比,验证了优化算法的时空精度。通过总结经验发现化学非平衡反应仅发生在流场局部区域,提出质量生成率判据,结合相应模拟结果验证了方法的可行性,可以进一步提高化学反应算子的计算效率。