辐射平衡壁温下有限催化模型数值模拟

文章发展了高超声速飞行器辐射平衡壁温下有限催化的数值方法,将数值模拟结果与基于返回舱外形的风洞实验数据进行了对比,并进一步针对典型高超声速飞行器钝双锥研究了辐射平衡壁温下有限催化对气动热环境的影响规律.针对返回舱外形的数值实验表明,完全催化与完全非催化边界条件下壁面热流密度均与风洞实验结果偏差较大,而采用合适的有限催化模型获得的壁面热流密度与风洞实验结果符合良好.针对典型高超声速飞行器钝双锥的研究表明,在辐射平衡温度边界条件下驻点附近氧原子的催化复合系数约为0.17,氮原子的催化复合系数约为0.026,大面积区则分别降为0.005 3和0.01.在驻点热流密度方面,完全催化的壁面热流峰值比有限催化高约21%,而完全非催化的壁面热流峰值比有限催化低约29%.      

临近空间高超声速飞行器气动力及气动热研究现状

临近空间高超声速飞行器具有速度快、突防能力强、杀伤力大等特点,是当今世界各军事强国新型武器的重点发展方向.其中,气动力和气动热是高超声速飞行器的两项重要指标,也是高超声速技术研究的重点内容.文章综述了国内外临近空间高超声速飞行器气动力及气动热研究现状,分析了研究的发展趋势,并分别从工程计算、数值仿真以及实验研究3个方面介绍了高超声速飞行器气动力及气动热的研究技术和方法.      

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

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

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

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

高超声速可压缩流中粗糙壁热流研究

用计算流体力学（CFD）数值模拟和理论方法对高超声速可压缩湍流中粗糙壁面热增量进行研究. 着重考虑粗糙单元密度和粗糙单元形状对粗糙表面热流的影响. 结果表明： ①粗糙单元密度变化时, CFD数值方法计算所得的粗糙单元等效热流随着粗糙单元密度的降低而增加, 理论方法预测结果的规律随方法不同而不同;②在相同粗糙单元密度和高度时,如果粗糙单元形状改变,CFD计算结果也随之发生改变. 理论预测方法所得结果不发生变化.     

基于遗传算法的升力体外形优化设计

升力体由于低热流率再入物理特性和高效的内部容积利用率，是高超声速飞行器气动外形的一种典型布局.文章对升力体飞行器进行参数化数值建模，并提取其表征外形的参数作为设计变量，综合考虑飞行器再入过程中的气动力、气动热、容积利用率及稳定性等性能指标.运用多目标混合遗传算法对升力体进行了多变量、多约束下的气动外形优化设计，获得了再入飞行器外形的最优Pareto解.数值模拟结果表明，典型状态下最优Pareto解与CFD结果相差12%，验证了优化结果的准确性.      

CFD在高速飞行器热气弹问题中的应用

采用气动力/热/结构耦合的方法对高速细长体飞行器结构热静气动弹性问题进行了研究.为保证耦合计算精度，达到准确预测热气动弹性特性的能力，气动力和气动热计算采用CFD数值模拟方法，热应力和热变形计算采用有限元方法并通过热考核试验验证.以该简单细长体飞行器模型为研究对象，对其热静气动弹性特性进行了计算与分析，计算结果表明:CFD/CSD耦合可准确模拟热气弹问题，且气动加热造成结构温升不均衡是结构变形的主导因素，力热耦合静气弹变形与单纯受力分析变形形式不同，对飞行器气动特性影响规律不同.准确预测飞行器热气动弹性特性对飞行器结构设计十分必要.      

基于Ludwieg管的高超声速边界层转捩实验

高超声速边界层层/湍流转捩是高超声速飞行器气动力和气动热设计中的难点和热点问题.为了降低开展高超声速边界层不稳定性与转捩实验研究的门槛,研究基于Ludwieg管原理设计并建造了一座Mach 6高超声速管风洞,重点对Ludwieg管风洞的启动和运行过程开展了数值模拟,分析了储气段弯管布局对试验段流场的影响;之后,对该高超声速风洞的自由来流品质进行了静态和动态的标定,验证了风洞的设计Mach数,并给出了流场的动态扰动特征;最后,基于7°半张角尖锥标模开展了高超声速边界层转捩实验,通过表面齐平式安装的高频PCB传感器获得边界层不稳定波,分析了高超声速边界层不稳定波的演化特征.以上工作表明,Ludwieg管相对常规高超声速风洞具有建设和运行成本低、运行效率高、流场品质好等优点,适合开展高超声速边界层转捩等基础实验研究.      

面向高超声速飞行器的激波智能预测方法

高超声速飞行器激波位置的准确预测能够有效提升数值模拟的精度和效率。一方面, 对高超声速飞行器激波附近网格进行正交和加密处理, 可有效提升数值计算精度; 另一方面, 使用高超声速飞行器激波位置对计算网格进行修正, 能够加速CFD计算收敛过程。提出了一种基于机器学习的高超声速飞行器激波智能预测方法, 对典型高超声速飞行器外形进行激波位置的高效准确预测。首先, 针对典型高超声速飞行器外形和典型飞行状态, 使用数值模拟方法获得收敛的流场, 并采用基于Mach数等值线的激波提取方法, 从流场中判别激波面并提取构成激波面的关键点位置, 形成训练数据; 然后采用有监督学习算法, 学习关键点位置, 并利用二次曲线沿流向拟合关键点形成初步的激波线族; 最后, 基于剖面压力云图, 构造基于投影压力图像的智能预测神经网络, 对初步形成的激波线族进行修正, 并获得三维激波面。大量的实验结果表明, 激波预测模型能够对高超声速飞行器激波位置做出准确预测, 预测的激波面与CFD数值计算结果中提取的激波面误差在10-4量级。      

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

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

求解飞机蒙皮耦合热效应的壁面热流函数法

分析了飞机内部换热与蒙皮外辐射换热和气动换热的耦合作用,提出了一种求解飞机内外耦合热效应的工程计算方法。通过构造反映蒙皮外辐射换热和气动换热的壁面热流函数,将蒙皮外热环境转化为浮动的热边界条件,实现了飞机蒙皮外部热环境求解与蒙皮内换热求解的解耦处理,避免了难以实现的直接耦合求解。在验证该方法可靠性的基础上, 采用该方法,并结合热网络法与FLUENT软件对某型飞机外部热环境及结构热响应问题进行了计算。     

DNS of turbulent flow with hemispherical wall roughness

The present study of the effect of roughness density on the mean flow turbulence parameters is motivated by the need for new generation of boundary conditions for multiphase computational multiphase fluid dynamics (CMFD) models applied to boiling flows. Effect of roughness element density on the turbulent flow in a channel is quantified through direct numerical simulations (DNSs). The Navier--Stokes equations are solved using finite element method and bubbles are approximated as rigid near-hemispherical obstacles at the wall. Six different cases were analysed including channel flow with smooth wall and channel flow with rough wall for five different bubble nucleation site densities. Friction factor and the law of the wall was calculated and compared with the previously published results. Existing correlations for nucleating bubble site density dependency on a wall heat flux were used to obtain a relation between the heat flux and the friction factor, leading to the law of the wall dependency on the heat flux. This separate effect study provides new guidelines on how the heat flux in subcooled boiling regime affects the turbulence behaviour near the wall and guides the computational fluid dynamics model development for boiling two-phase flows.     

交叉三角形波纹板流道在过渡流状态下的传热与阻力特性

交叉三角形波纹板流道中流动常处于过渡态。本文研究了交叉波纹板中的周期性完全发展流动及热传递。利用周期性降低几何流道的复杂性以及简化模拟对象。为了模拟这个拓扑结构中的过渡流,利用了已经验证的低雷诺数k-ε,湍流模型来说明流动中的湍流流动。得到了三维复杂计算区的温度、速度以及湍流场。计算了在恒壁温和恒热流密度两种边界条件下的摩擦系数和平均Nusselt数及其与雷诺数的关系。湍流中心从上层壁面的波纹处移向下层壁面的波纹处并逐步增强。     

Numerical investigation on properties of attack angle for an opposing jet thermal protection system

The three-dimensional Navier-Stokes equation and the k-ε viscous model are used to simulate the attack angle characteristics of a hemisphere nose-tip with an opposing jet thermal protection system in supersonic flow conditions.The numerical method is validated by the relevant experiment.The flow field parameters,aerodynamic forces,and surface heat flux distributions for attack angles of 0°,2°,5°,7°,and 10° are obtained.The detailed numerical results show that the cruise attack angle has a great influence on the flow field parameters,aerodynamic force,and surface heat flux distribution of the supersonic vehicle nose-tip with an opposing jet thermal protection system.When the attack angle reaches 10°,the heat flux on the windward generatrix is close to the maximal heat flux on the wall surface of the nose-tip without thermal protection system,thus the thermal protection has failed.     

Numerical investigation on properties of attack angle for opposing jet thermal protection system

The three-dimensional Navier-Stokes equation and the k-ε viscous model are used to simulate the attack angle characteristics of a hemisphere nose-tip with an opposing jet thermal protection system in supersonic flow condition. The numerical method is validated by the relevant experiment. The flow field parameters, aerodynamic forces, and surface heat flux distributions for attack angles of 0°, 2°, 5°, 7°, and 10° are obtained. The detailed numerical results show that the cruise attack angle has a great influence on the flow field parameters, aerodynamic force, and surface heat flux distribution of the supersonic vehicle nose-tip with opposing jet thermal protection system. When the attack angle reaches 10°, the heat flux on the windward generatrix is close to the maximal heat flux on the wall surface of the nose-tip without thermal protection system, thus the thermal protection is failure.     

A numerical analysis on flow field and heat transfer by interaction between a pair of vortices in rectangular channel flow

This paper is a numerical study of the effect of flow field and heat transfer created by interactions between a pair of vortices generated by a vortex generator in a rectangular channel flow. In order to analyze the vortices produced by the vortex generator, the pseudo-compressibility method is introduced into the Navier–Strokes (NS) equation of a three-dimensional unsteady, incompressible viscous flow. A two-layer k–ε turbulence model is used on the flat plate three-dimensional turbulence boundary to predict the turbulence characteristics of the vortices. The computational results accurately predict the vortex characteristics, which are related to Reynolds stress, turbulent kinetic energy, and flow field. Also, in the prediction of thermal boundary layers, skin friction characteristics, and heat transfers, the present results are reasonably close to the experimental results obtained by other researchers.     

Double-diffusive radiative magnetic mixed convective slip flow with Biot and Richardson number effects

This paper investigates combined heat and mass transfer by mixed magneto-convective flow of an electrically conducting flow along a moving radiating vertical flat plate with hydrodynamic slip and thermal convective boundary conditions. The governing transport equations are converted into a system of coupled nonlinear ordinary differential equations with prescribed boundary conditions using similarity variables developed by Lie group theory. The transformed nondimensional boundary value problem is then solved numerically with MAPLE13 quadrature. Excellent correlation with previous nonmagnetic, no-slip studies is achieved. Surface shear stress function and local Nusselt number (heat transfer gradient at the wall) are increased with Richardson number, whereas local Sherwood number is found to initially decrease then subsequently increase. The “thermally thick” scenario (Biot number > 0.1) is investigated and increasing Biot number is observed to enhance shear stress function (skin friction), local Nusselt number, and local Sherwood number. Increasing thermal radiation flux increases thermal boundary layer thickness as does increasing the magnetic field effect. Increasing hydrodynamic slip parameter reduces skin friction but enhances local Nusselt and Sherwood numbers. The study has applications in high-temperature polymeric synthesis and magnetic field flow control.     

高超声速飞行器前缘流固耦合计算方法研究

对高超声速流场和结构温度场进行了耦合计算分析, 同时基于准静态假设对结构应力进行了分析. 流场部分采用基于非定常Navier-Stokes (N-S)方程的有限体积法, 湍流模型采用SST k-ω 模型, 固体部分采用基于非稳态热传导方程的有限元法, 同时基于准静态假设对固体结构的应力应变进行了分析. 在流固交界面处, 高速流体从固体结构得到温度边界条件, 固体结构从高速流体得到热流边界条件, 从而实现了流场和固体温度场的紧耦合计算.通过与超声速无限长圆管绕流试验结果进行对比, 验证了该方法的可靠性. 同时对二维圆管结构在气动加热过程中的温度、应力等的变化进行了比较详细的分析. 研究结果表明: 随着气动加热时间的推进, 由于圆管结构的高温区在不断扩大, 导致了结构的热变形在不断地增大; 圆管最小变形区出现在θ为60°处; 同时研究发现在计算时间内圆管热变形对外部流场的影响可以忽略不计.     

Experimental characterization of the convective heat transfer in a vortex-wall interaction

The development of turbulence models and wall laws for the numerical simulation of flows in complex geometries requires a detailed experimental analysis of turbulence and of the phenomena that appear in turbulent boundary layers. There is a strong need to develop new measurement systems allowing the determination of unsteady wall heat transfer coefficients. In order to improve the knowledge of the unsteady phenomena occurring in perturbed boundary layers, a fundamental study is conducted on the interaction of a single vortex with a flat plate. An experimental methodology using a specific thermal sensor whose surface temperature is measured by an infrared thermography system is presented. It allows the characterization of the unsteady convective heat transfer coefficient whose evolution is compared with the fluctuations of the wall friction coefficient, calculated from velocity profiles measured by laser Doppler velocimetry.