进入轨道偏差对火星科学实验室气动力特性的影响

针对火星科学实验室(MSL)高超声速进入过程,利用三维并行程序求解流体力学Navier-Stokes方程,耦合真实气体模型,分析火星大气中真实气体效应对进入器气动力特性的影响量在进入轨道发生偏差时的变化规律.结果表明:对海盗号的计算结果与飞行数据符合很好,验证了火星大气真实气体模型和计算方法;真实气体效应影响下,激波层厚度大为减小,温度下降明显,进入器阻力系数明显增加,升力系数变化不大,俯仰力矩系数增加,基准状态下配平攻角较完全气体减小约2.2°;高度不变,Ma数增加导致阻力系数和俯仰力矩系数增大,配平攻角和完全气体的差值由1.6°增加到2.6°,表明Ma数变大时真实气体效应引起的气动力变化增强;Ma数不变,高度增加略微减弱波后化学反应,对进入器气动力特性基本没有影响.     

双锥/双楔流动中的高温气体效应仿真模拟

在高超声速飞行和再入地球大气过程中, 气体分子的振动、电子态激发, 伴随离解、电离反应, 从而产生高温真实气体效应。不同数值方法对高温真实气体效应的模化会造成气体热物性参数的差异, 从而对流场模拟引入不确定度。以高超声速的双锥/双楔流动为例, 通过计算流体力学方法和直接模拟Monte Carlo (DSMC)方法, 研究高温真实气体模型对复杂干扰流动的预测能力。结果表明, 有别于量热完全气体, 若考虑真实物理过程的热化学非平衡过程带来气体热力学性质、输运特性的变化, 会导致激波角、边界层厚度、分离区尺寸等流动结构的改变。因此, 在研究高超声速模拟中应注意数值模型的正确应用。      

标模外形化学非平衡流场数值模拟

利用三维并行计算代码求解Navier-Stokes方程,数值模拟标模(ELECTRE)化学非平衡绕流,研究真实气体效应对标模气动热特性的影响,反应模型为Dunn和Kang的7组元7反应化学动力学模型.利用典型弹道点的飞行试验数据验证化学非平衡流计算程序的可靠性.在此基础上,研究不同壁面催化条件下攻角和高度变化对热流的影响.计算表明:真实气体效应主要发生在物面附近很薄的激波层内,并使激波脱体距离减小;完全催化壁驻点热流值高于非催化壁热流值;随着攻角增大,热流分布差异明显,而且攻角越大时,物面电子数密度越小;飞行高度越高,O₂和N₂离解程度越低,驻点热流越低.     

超临界水在倾斜上升管内的非均匀传热特性数值模拟

对超临界压力水在管径为Φ32mm×3mm、长度为8000mm、倾角为20°的倾斜光管内的流动与传热特性进行数值模拟研究,重点考察超临界水在大比热区内的奇异物性变化对倾斜管内的传热特性的影响.通过与实验数据的对比,验证计算模型的正确性;计算压力P=26,34 MPa时,不同质量流速和热负荷下倾斜光管内壁温随超临界水焓值增大而变化的规律,以及倾斜光管内壁周向温度及热负荷分布的不均匀性;分析大比热区管内上下母线处内壁温度差随工质焓值变化的特性及机理,讨论大比热区水的物性变化对倾斜光管内传热不均匀性的影响;引入截面横向动能与截面相对横向动能两个变量对二次流的强度进行描述,重点考察超临界水在大比热区内的二次流的流动特点及其对倾斜管内传热特性影响的机理.并利用截面中垂线上的密度梯度定量分析二次流变化的规律,讨论超临界压力下大比热区内倾斜光管内壁温分布异常的机理.     

飞行器大攻角非定常气动特性神经网络建模

大攻角气动特性预测与气动建模是新型飞行器提升飞行性能的重要内容.以轴对称导弹简化模型为研究对象，首先采用计算流体力学方法，对70°大攻角状态的非定常气动特性进行数值模拟，计算方法基于RANS的N-S方程，湍流模型采用SA模型，对流场采用有限体积法离散，无黏项采用Roe通量差分分裂格式，黏性项采用中心差分，时间推进采用LU-SGS格式的双时间步法.飞行器运动模式采用强迫振荡的方式，对5种不同振荡频率进行了非定常数值计算，并记录每一内迭代周期最终的气动力和力矩数值.其次，以CFD预测结果作为气动建模的样本，采用动导数模型、多项式模型等传统方法，进行气动建模，并分析其有效性和精度.最后采用神经网络方法对大攻角非定常气动力进行建模，并和动导数模型、多项式模型进行精度对比.结果表明，基于神经网络的人工智能气动建模方法具有较高的精度和适应性.该方法为飞行器大攻角非定常非线性气动建模，大攻角飞行稳定性分析与控制提供理论参考.      

高超声速层流绕三维压缩拐角的数值模拟

通过数值模拟研究了高超声速来流绕过压缩拐角的层流分离三维流动特性.数值方法采用三维N-S方程,结合2阶精度Roe格式以及分区结构网格有限体积法进行离散.数值模拟的空间激波结构与实验纹影结果符合较好;激波/边界层干扰区内3条纵向线上的计算压力分布与实验结果进行了对比分析,计算获得在三维楔侧面存在低压力区,与实验结果反映的规律一致,计算结果表明低压力区是由楔体侧缘尖端发起的二次涡的抽吸作用造成的.此外,在楔体后端尾流区的低压沿边界层内的亚声速区往上游传递了一定距离.      

基于汽相/液相界面跟踪方法的弯管内空化流动的数值预测

采用发展的基于Reynolds-Averaged Navier-Stokes方程和液相/气相界面跟踪方法的单相空化模型和数值迭代算法,数值预测了二维弯管内空化卒泡的形状和相应的压力系数分布.将发展的二维空化模型和算法推广到i维空化流动的数值模拟,预测了三维非对称弯管内的三维空泡形状和流动特性.数值预测三维非对称弯管内空化流动时的卒泡位置和相应压力系数分布符合理论分析结论.研究结果证明了所发展的空化模型和算法能够应用于三维空化流动时的数值预测.     

波状前缘对风力机翼型失速性能影响数值分析

为研究仿生波状前缘对翼型失速性能的影响,本文采用S-A湍流模型,对风力机翼型NACA634-021(光滑前缘)以及对应的正弦波状前缘仿生翼型的绕流流场进行了数值模拟。结果表明,光滑翼型在20°攻角附近发生深度失速,升力系数骤然下降;而波状前缘仿生翼型有效改善了失速特性,升力系数变化较平稳,在大攻角下高于光滑翼型。通过流场分析发现光滑翼型失速前后升力系数骤然下降的主要原因在于前缘压力面和吸力面的压差大幅度下降,而仿生翼型改变了前缘的压力分布特性,进而改变了大攻角下的分离特性,促进流向涡对的产生和发展,使得凸峰附近保持附着流动,进而提高升力。     

燃料空气供给条件对热声不稳定燃烧影响的研究

对Solar低排放预混燃烧系统的燃烧稳定性进行了数值研究.应用非定常N-S方程、雷诺应力紊流模型及涡团耗散燃烧模型,数值模拟了该类型燃烧器在不同的燃料空气供给条件下的气流流动特性和压力振荡特性,并给出了不稳定发生时压力和速度振荡的幅值和频率.根据供给条件的不同,燃烧可以是稳定的或是不稳定的,取决于燃料到火焰前沿的迟滞时间.采用CFD方法,可精确地获得燃料到火焰前沿的迟滞时间,证实了所采用的模型能够精确预测不稳定燃烧的出现及振荡特性.通过调整燃料与空气的供给条件,可使振荡激励或阻尼.     

重费米子材料CeCuGa的低温物性研究

利用溶融和热扩散方法,我们合成了金属化合物CeCuGa.X射线衍射表明样品为单相的正交结构,分别利用稳定绝热脉冲法和四引线法测量了样品的低温比热及电阻率,电阻率在低温区出现Kondo极小,10K以上的比热数据与电子项γT和声子项βT~3符合的很好,在2K左右出现Kondo宽峰,2K以下的数据外推到绝对零度得到γ=1.46J/K~2mole,表明CeCuGa是一种典型的重费米子材料。     

Influence of gas law on ultrasonic behaviour of porous media under pressure

This paper deals with the influence of gas law on ultrasonic behaviour of porous media when the saturating fluid is high pressured. Previous works have demonstrated that ultrasonic transmission through a porous sample with variations of the static pressure (up to 18 bars) of the saturating fluid allows the characterization of high damping materials. In these studies, the perfect gas law was used to link static pressure and density, which is disputable for high pressures. This paper compares the effects of real and perfect gas laws on modeled transmission coefficient for porous foams at these pressures. Direct simulations and a mechanical parameters estimation from minimization show that results are very similar in both cases. The real gas law is thus not necessary to describe the acoustic behaviour of porous media at low ultrasonic frequencies (100 kHz) up to 20 bars.     

On the effect of viscosity and thermal conductivity on sound propagation in ducts: A re-visit to the classical theory with extensions for higher order modes and presence of mean flow

The dispersion equation for the axisymmetric modes of viscothermal acoustic wave propagation in uniform hard-walled circular ducts containing a quiescent perfect gas is classical. This has been extended to cover the non-axisymmetric modes and real fluids in contemporary studies. The fundamental axisymmetric mode has been the subject of a large number of studies proposing approximate solutions and the characteristics of the propagation constants for narrow and wide ducts with or without mean flow is well understood. In contrast, there are only few publications on the higher order modes and the current knowledge about their propagation characteristics is rather poor. On the other hand, there is a void of papers in the literature on the effect of the mean flow on the quiescent modes of propagation. The present paper aims to contribute to the filling of these gaps to some extent. The classical theory is re-considered with a view to cover all modes of acoustic propagation in circular ducts carrying a real fluid moving axially with a uniform subsonic velocity. The analysis reveals a new branch of propagation constants for the axisymmetric modes, which appears to have escaped attention hitherto. The solution of the governing wave equation is expressed in a modal transfer matrix form in frequency domain and numerical results are presented to show the effects over wide ranges of frequency, viscosity and mean flow parameters on the propagation constants. The theoretical formulation allows for the duct walls to have finite impedance, but no numerical results are presented for lined ducts or ducts carrying a sheared mean flow.     

Numerical analysis of self-excited tangential combustion instability for an MMH/NTO rocket combustor

This study presented the numerical simulation of the tangential combustion instability in a hypergolic liquid bipropellant rocket thrust chamber, which applied fuel liquid film cooling method and unlike impinging injectors. The liquid spray was modeled using Lagrangian approach, while the gas was regarded as Euler phase. Stress-blended eddy simulation and finite rate/eddy–dissipation model were adopted to simulate the turbulent combustion process. Consistent with the experiment results, this work successfully simulated the transformation of tangential combustion instability from standing mode to spinning mode. The mean pressure, amplitude and frequency of limit cycle oscillation were in good agreement with the experiment. There was a detailed analysis about the flow field, Rayleigh index, and driving mechanism of the combustion instability. It was found that the oscillation began with hot spots of heat release rate due to the interaction between the spray of impinging injectors and cooling fuel jet. More than that, cooling fuel jet also contributed to drive the oscillation. In the standing mode, injectors in the inner and outer rings drive the oscillation together, while the spinning mode is mainly driven by injectors in the outer ring. The pressure wave is subsonic and its Mach number is close to 1. It was shown that the pressure wave contained a complex structure divided into three parts. This led to the in-phase of the pressure along the axial direction and the double-peak characteristic of the downstream pressure signal. Besides, a positive feedback closed-loop system associated with periodic oxidizer/fuel ratio was believed to sustain the combustion instability. The oscillation can be maintained when pressure, heat release and oxidizer/fuel ratio are coupled together. The analysis results indicate that rotating detonation is an implication to tangential combustion instability.     

Pressure based finite volume method for calculation of compressible viscous gas flows

A pressure based, iterative finite volume method is developed for calculation of compressible, viscous, heat conductive gas flows at all speeds. The method does not need the use of under-relaxation coefficient in order to ensure a convergence of the iterative process. The method is derived from a general form of system of equations describing the motion of compressible, viscous gas. An emphasis is done on the calculation of gaseous microfluidic problems. A fast transient process of gas wave propagation in a two-dimensional microchannel is used as a benchmark problem. The results obtained by using the new method are compared with the numerical solution obtained by using SIMPLE (iterative) and PISO (non-iterative) methods. It is shown that the new iterative method is faster than SIMPLE. For the considered problem the new method is slightly faster than PISO as well. Calculated are also some typical microfluidic subsonic and supersonic flows, and the Rayleigh–Bénard convection of a rarefied gas in continuum limit. The numerical results are compared with other analytical and numerical solutions.     

基于真实气体效应的湍流标度律特性

在量热完全气体、热完全气体和化学反应完全气体等3种气体模型假设下,利用Mach数为4.05、壁温为1 300 K的超声速槽道湍流的直接数值模拟（direct numerical simulation,DNS）结果,对标度律和自相似性做了详细分析.结果表明,不仅在量热完全气体模型下存在标度律和扩展自相似性,而且在热完全气体和化学反应完全气体模型下标度律和扩展自相似性仍然成立.压缩性的影响使得速度结构函数通过Favre平均获得更为合适.与热完全气体模型的结果相比,化学反应完全气体和量热完全气体模型的结果吻合更好.      

High-order shock-fitting methods for direct numerical simulation of hypersonic flow with chemical and thermal nonequilibrium

In recent years, much progress has been made in the direct numerical simulation of laminar-turbulent transition of hypersonic boundary layer flow. However, most of the efforts at the direct numerical simulation of transition previously have been focused on the idealized perfect gas flow or “cold” hypersonic flows. For practical problems in hypersonic flows, high-temperature effects of thermal and chemical nonequilibrium are important and cannot be modeled by a perfect gas model. Therefore, it is necessary to include the real gas models in the numerical simulation of hypersonic boundary layer transition in order to accurately predict flow field parameters. Currently most numerical methods for hypersonic flow with thermo-chemical nonequilibrium are based on shock-capturing approach at relatively low order of accuracy. Shock capturing schemes reduce to first-order accuracy near the shock and have been shown to produce spurious oscillations behind curved strong shocks. There is a need to develop new methods capable of simulating nonequilibrium hypersonic flow fields with uniformly high-order accuracy and avoid spurious oscillations near the shock. This paper presents a fifth-order shock-fitting method for numerical simulation of thermal and chemical nonequilibrium in hypersonic flows. The method is developed based on the state-of-the-art real gas models for thermo-chemical nonequilibrium and transport phenomena. Shock-fitting approach is used because it has the advantage of capturing the entire flow field with high-order accuracy and without any oscillations near the shock. The new method has been tested and validated for a number of test cases over a wide span of free stream conditions. The developed method is applied for the study of receptivity of free stream acoustic waves over a blunt cone for hypervelocity flow. Some preliminary results of the computations of the high order shock fitting method for the above mentioned study have also been presented.     

Numerical investigation of hydrodynamics behaviour of melt layer during laser cutting of steel

Understanding of the melt layer hydrodynamic behaviour during laser-cutting process under gas jet assistance is of high importance for cut quality control. In the present work, a numerical model is developed to calculate the three-dimensional behaviour of the melt flow on the kerf front, while an inert gas jet interacts with the melt film. Fluent CFD code is used to solve the governing hydrodynamic equations by finite volume method. The results show that the melt flow on the kerf front reveals a strong instability, which depends on the cutting speed and on the gas jet velocity. Global flow behaviour (gas and molten metal flows) computed using a laminar model, reveals oscillations of the gas–metal liquid interface, which is assimilated to Kelvin–Helmholtz instability. The origin of this instability is discussed in terms of instabilities in thermal dynamics and hydrodynamics. Instability in thermal dynamics is related to the localized melting, while the instability in hydrodynamics is governed by forces balance between gas and resistant surface tension.     

PIV技术在气体界面不稳定性实验研究中的应用

采用粒子图像测速(Particle Image Velocimetry,PIV)技术,以乙二醇烟雾作为示踪粒子,实验测量了平面激波作用下SF6气柱-空气界面Richtmyer-Meshkov不稳定性演化图像和二维速度场.测量结果揭示了不稳定性流场的典型特征和细微结构,与高速摄影和数值模拟结果符合,说明建立的PIV技术适...