二维高超声速后台阶表面传热特性实验研究

高超声速后台阶流动是大气层内高速飞行器发动机设计、表面热防护以及高超声速拦截器红外成像窗口气动光学效应校正等诸多先进高超声速技术研发过程中所涉及的一类基础流动问题. 研究高超声速后台阶流动特性对有效提升飞行器综合性能, 进一步掌握高超声速流动机理具有重大基础 意义. 本文以二维高超声速后台阶流动为研究对象, 在KD-01高超声速激波风洞中测量了二维后台阶模型表面传热系数和表面静压, 并将实测台阶下游表面传热系数分布同采用高超声速边界层理论所得估计值进行了比较. 为进一步验证实验结果, 使用NPLS技术测量了其中一种实验状态下台阶周围流动结构. 研究发现, 对于二维高超声速后台阶流动, 台阶下游表面传热分布受台阶处边界层外缘流动特性的直接影响; 在台阶下游分离区和再附区内, 气体黏性占主导作用; 在台阶下游远场区域, 边界层流动特性趋同于平板边界层; 下游边界层基本结构取决于台阶处边界层相对厚度. 对高超声速后台阶流动, 若使用数值模拟方法研究气动热问题, 应当使用湍流模型.     

跨流域高超声速绕流环境Boltzmann模型方程统一算法研究

如何准确可靠地模拟从外层空间高稀薄流到近地面连续流的航天器高超声速绕流环境与复杂流动变化机理是流体物理的前沿基础科学问题. 基于对Boltzmann方程碰撞积分的物理分析与可计算建模, 确立了可描述自由分子流到连续流区各流域不同马赫数复杂流动输运现象统一的Boltzmann模型速度分布函数方程, 发展了适于高、低不同马赫数绕流问题的离散速度坐标法和直接求解分子速度分布函数演化更新的气体动理论数值格式, 建立了模拟复杂飞行器跨流域高超声速飞行热环境绕流问题的气体动理论统一算法. 对稀薄流到连续流不同Knudsen数0.002 ≤Kn_∞ ≤1.618、不同马赫数下可重复使用卫星体再入过程(110–70 km)中高超声速绕流问题进行算法验证分析, 计算结果与典型文献的Monte Carlo直接模拟值及相关理论分析符合得较好. 研究揭示了飞行器跨流域不同高度高超声速复杂流动机理、绕流现象与气动力/热变化规律, 提出了一个通过数值求解介观Boltzmann模型方程, 可靠模拟高稀薄自由分子流到连续流跨流域高超声速气动力/热绕流特性统一算法.     

高超声速三维热化学非平衡电离流动计算研究

研究发展了一种计算高超声速三维热化学非平衡电离流动的数值方法。真实气体热力学状态是通过平动-转动温度和电子-振动温度来模拟。对零攻角轴对称钝锥体和有攻角钝锥体的高超声速电离流动,采用隐式有限体积法NND格式,数值求解NS方程。     

高速动能弹温度场数值模拟

随着兵器发射技术和空气动力学技术的发展,动能弹的发射初速和飞行状态正从超声速向高超声速发展,由此产生了气动热问题.准确预测动能弹温度场是其气动力和热防护设计的关键技术.采用CFD预测温度场的方法,包括平衡流流动控制方程及差分格式,构造平衡流通量Jacob矩阵,在差分格式矢通量分裂过程中嵌入平衡流真实气体模型模拟温度场,获得平衡流气体状态方程.对典型高速动能弹热环境进行验证,考察方法的合理性.对设计的一种新型高超声速动能弹温度场进行数值模拟,为其气动设计及热防护提供了较可靠的数据.     

具有局部稀薄气体效应的高超声速尖锥气动加热特征研究

以微钝尖锥为飞行器前缘模型,采用基于分子运动论的DSMC方法模拟不同前缘曲率半径的尖锥在高超声速来流下的气动热环境,计算驻点热流率,并与Fay-Riddell公式和其他修正理论作对比,研究具有局部稀薄气体效应的高超声速尖锥气动加热特征及其变化规律.发现修正的Cheng参数适合作为工程上判断驻点区域稀薄气体效应影响大小的判据.     

壁面催化对再入飞行器等离子体鞘套及电磁参数的影响

针对超高声速流动中的高温真实气体效应,采用数值模拟求解考虑化学非平衡的三维Navier-Stokes(N-S)方程,研究了壁面催化对典型再入飞行器等离子体鞘套及电磁参数的影响规律.研究发现:(1)完全催化壁条件下等离子体密度计算结果与飞行试验符合较好;(2)完全催化壁条件下,离解原子在壁面的复合导致波后气体可压缩性增强...     

热化学模型对高超声速磁流体控制数值模拟影响分析

针对等离子体流场的模拟准确性问题及其对高超声速磁流体控制的影响,通过数值求解三维非平衡Navier-Stokes流场控制方程和Maxwell电磁场控制方程,建立了三维低磁雷诺数磁流体数值模拟方法及程序,分析了不同空气组分化学反应模型和壁面有限催化效率等因素对高超声速磁流体控制的影响.研究表明:不同空气组分化学反应模型对高超声速磁流体流场结构、气动力/热特性控制的影响不容忽视;对于本文计算条件,Park化学反应模型在组分模型一致性、等离子体模拟准确性等方面具有一定优势;磁控热防护效果,受壁面有限催化复合系数影响较大,两者呈非线性关系,不同表面区域差异较大;磁场对磁阻力伞及其磁阻力特性影响,受壁面催化效应的影响相对较小.     

吸气式高超声速飞行器气动力气动热的数值模拟方法及应用

对吸气式高超声速飞行器而言,物面热流和摩阻的准确预测对飞行器设计及安全十分关键.介绍采用CFD准确预测气动力和气动热的方法,包括流动的控制方程、湍流模型及湍流的先进壁面函数边界条件,介绍流动的数值求解方法.对典型超声速层流和湍流流动的摩擦阻力和热流进行详细的验证与确认,考察CFD工具在使用先进壁面函数边界条件后,湍流计算的法向网格无关性能力.对设计的一种吸气式高超声速飞行器的气动力和气动热进行数值模拟,为飞行器的气动设计及热防护提供了可靠的数据.     

颗粒增强湍流的新模型和气粒流动的湍流变动

用雷诺应力方程模型和极细的网格系对单个颗粒受湍流气体绕流进行了数值模拟,研究了改变颗粒直径和气体相对速度时颗粒增强气体湍流的规律.据此构造了颗粒尾涡增强气体湍流的新模型.将此子模型加入到两相流动模型中,对竖直和水平通道内气粒两相流动进行了数值模拟,和实验结果的对照表明,考虑颗粒尾涡增强气体湍流效应得到的气体湍流脉动速度的模拟结果比不考虑此效应的模拟结果好得多.     

高精度非定常激波装配法

为正确模拟高超声速绕流中,来流小扰动与弓形激波之间的干扰对流动特征的影响,将弓形激波作为动边界,利用非定常特征关系处理激波处的边界条件.应用五阶精度迎风紧致格式和六阶精度的对称格式与三阶精度的R-K方法相结合,建立高精度非定常激波装配方法.采用该方法数值模拟钝锥高超声速定常流场和二维抛物外形高超声速边界层流动的感受性问题,数值模拟来流小扰动与弓形激波干扰激波后非定常扰动流场,研究扰动波进入边界层产生边界层不稳定波的特征.     

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.     

Gas-kinetic numerical studies of three-dimensional complex flows on spacecraft re-entry

The gas-kinetic numerical algorithm solving the Boltzmann model equation is extended and developed to study the three-dimensional hypersonic flows of spacecraft re-entry into the atmosphere in perfect gas. In this study, the simplified velocity distribution function equation for various flow regimes is presented on the basis of the kinetic Boltzmann–Shakhov model. The discrete velocity ordinate technique and numerical quadrature methods, such as the Gauss quadrature formulas with the weight function 2/π^1/2exp(?V²) and the Gauss–Legendre numerical quadrature rule, are studied to resolve the barrier in simulating complex flows from low Mach numbers to hypersonic problems. Specially, the gas-kinetic finite-difference scheme is constructed for the computation of three-dimensional flow problems, which directly captures the time evolution of the molecular velocity distribution function. The gas-kinetic boundary conditions and numerical procedures are studied and implemented by directly acting on the velocity distribution function. The HPF (high performance fortran) parallel implementation technique for the gas-kinetic numerical method is developed and applied to study the hypersonic flows around three-dimensional complex bodies. The main purpose of the current research is to provide a way to extend the gas-kinetic numerical algorithm to the flow computation of three-dimensional complex hypersonic problems with high Mach numbers. To verify the current method and simulate gas transport phenomena covering various flow regimes, the three-dimensional hypersonic flows around sphere and spacecraft shape with different Knudsen numbers and Mach numbers are studied by HPF parallel computing. Excellent results have been obtained for all examples computed.     

高超声速湍流流场高折射率梯度区域气动光学畸变仿真研究

基于折射率界面厚度的描述建立了一种高折射率梯度门限的数学模型,在此梯度门限下,研究了高超声速流场中高折射率梯度区域的气动光学传输效应.提出了一种用折射率梯度的调和平均值描述高折射率梯度门限的方法.采用高超声速流场的计算流体力学结果作为分析折射率梯度和进行气动光学传输仿真的源数据,忽略绝对值低于该门限的梯度值重构折射率场,并采用变折射率介质中光线追迹算法仿真其气动光学传输畸变.不同流场状况、不同位置截面的仿真结果表明,采用本门限,重构折射率场和原折射率场的相关性达0.9以上,仿真光程差均方根的相对误差不超过±5%,验证了该高折射率梯度门限模型的有效性和适用性,同时从数值角度证实了高超声速湍流流场中高折射率梯度区域是气动光学传输畸变的主要成因.     

Study on the unified algorithm for three-dimensional complex problems covering various flow regimes using Boltzmann model equation

The Boltzmann simplified velocity distribution function equation describing the gas transfer phenomena from various flow regimes will be explored and solved numerically in this study. The discrete velocity ordinate method of the gas kinetic theory is studied and applied to simulate the complex multi-scale flows. Based on the uncoupling technique on molecular movement and colliding in the DSMC method, the gas-kinetic finite difference scheme is constructed to directly solve the discrete velocity distribution functions by extending and applying the unsteady time-splitting method from computational fluid dynamics. The Gauss-type discrete velocity numerical quadrature technique for different Mach number flows is developed to evaluate the macroscopic flow parameters in the physical space. As a result, the gas-kinetic numerical algorithm is established to study the three-dimensional complex flows from rarefied transition to continuum regimes. The parallel strategy adapted to the gas-kinetic numerical algorithm is investigated by analyzing the inner parallel degree of the algorithm, and then the HPF parallel processing program is developed. To test the reliability of the present gas-kinetic numerical method, the three-dimensional complex flows around sphere and spacecraft shape with various Knudsen numbers are simulated by HPF parallel computing. The computational results are found in high resolution of the flow fields and good agreement with the theoretical and experimental data. The computing practice has confirmed that the present gas-kinetic algorithm probably provides a promising approach to resolve the hypersonic aerothermodynamic problems with the complete spectrum of flow regimes from the gas-kinetic point of view of solving the Boltzmann model equation. Supported by the National Natural Science Foundation of China (Grant Nos. 90205009 and 10321002) and the National Parallel Computing Center     

Numerical Simulations of Unsteady Flows from Rarefied Transition to Continuum Using Gas-Kinetic Unified Algorithm

Numerical simulations of unsteady gas flows are studied on the basis of Gas-Kinetic Unified Algorithm (GKUA) from rarefied transition to continuum flow regimes. Several typical examples are adopted. An unsteady flow solver is developed by solving the Boltzmann model equations, including the Shakhov model and the Rykov model etc. The Rykov kinetic equation involving the effect of rotational energy can be transformed into two kinetic governing equations with inelastic and elastic collisions by integrating the molecular velocity distribution function with the weight factor on the energy of rotational motion. Then, the reduced velocity distribution functions are devised to further simplify the governing equation for one- and two-dimensional flows. The simultaneous equations are numerically solved by the discrete velocity ordinate (DVO) method in velocity space and the finite-difference schemes in physical space. The time-explicit operator-splitting scheme is constructed, and numerical stability conditions to ascertain the time step are discussed. As the application of the newly developed GKUA, several unsteady varying processes of one- and two-dimensional flows with different Knudsen number are simulated, and the unsteady transport phenomena and rarefied effects are revealed and analyzed. It is validated that the GKUA solver is competent for simulations of unsteady gas dynamics covering various flow regimes.     

Lattice Boltzmann computational fluid dynamics in three dimensions

The recent development of the lattice gas method and its extension to the lattice Boltzmann method have provided new computational schemes for fluid dynamics. Both methods are fully paralleled and can easily model many different physical problems, including flows with complicated boundary conditions. In this paper, basic principles of a lattice Boltzmann computational method are described and applied to several three-dimensional benchmark problems. In most previous lattice gas and lattice Boltzmann methods, a face-centered-hyper-cubic lattice in four-dimensional space was used to obtain an isotropic stress tensor. To conserve computer memory, we develop a model which requires 14 moving directions instead of the usual 24 directions. Lattice Boltzmann models, describing two-phase fluid flows and magnetohydrodynamics, can be developed based on this simpler 14-directional lattice. Comparisons between three-dimensional spectral code results and results using our method are given for simple periodic geometries. An important property of the lattice Boltzmann method is that simulations for flow in simple and complex geometries have the same speed and efficiency, while all other methods, including the spectral method, are unable to model complicated geometries efficiently.     

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

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

Variation character of stagnation point heat flux for hypersonic pointed bodies from continuum to rarefied flow states and its bridge function study

This paper is a research on the variation character of stagnation point heat flux for hypersonic pointed bodies from continuum to rarefied flow states by using theoretical analysis and numerical simulation methods. The newly developed near space hypersonic cruise vehicles have sharp noses and wingtips, which desires exact and relatively simple methods to estimate the stagnation point heat flux. With the decrease of the curvature radius of the leading edge, the flow becomes rarefied gradually, and viscous interaction effects and rarefied gas effects come forth successively, which results in that the classical Fay-Riddell equation under continuum hypothesis will become invalid and the variation of stagnation point heat flux is characterized by a new trend. The heat flux approaches the free molecular flow limit instead of an infinite value when the curvature radius of the leading edge tends to 0. The physical mechanism behind this phenomenon remains in need of theoretical study. Firstly, due to the fact that the whole flow regime can be described by Boltzmann equation, the continuum and rarefied flow are analyzed under a uniform framework. A relationship is established between the molecular collision insufficiency in rarefied flow and the failure of Fourier’s heat conduction law along with the increasing significance of the nonlinear heat flux. Then based on an inspiration drew from Burnett approximation, control factors are grasped and a specific heat flux expression containing the nonlinear term is designed in the stagnation region of hypersonic leading edge. Together with flow pattern analysis, the ratio of nonlinear to linear heat flux W _r is theoretically obtained as a parameter which reflects the influence of nonlinear factors, i.e. a criterion to classify the hypersonic rarefied flows. Ultimately, based on the characteristic parameter W _r , a bridge function with physical background is constructed, which predicts comparative reasonable results in coincidence well with DSMC and experimental data in the whole flow regime.     

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

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