探路者号火星探测器气动热和传热耦合分析

本文建立高超声速气动热和结构传热的松耦合计算方法,以探路者号火星探测器为研究对象,开展了探测器高超声速气动热和结构传热的耦合计算研究,分析了探测器进入条件下某轨道点上气动加热和结构传热机理。耦合计算表明,随着时间的推进,表面结构温度逐渐升高,壁面热流相对降低,表面趋向辐射平衡温度。因热防护系统结构传热时间尺度长于飞行器进入过程总时间,探测器进入的真实过程滞后于辐射平衡过程。耦合研究表明,使用耦合计算方法能较好地再现真实的气动加热和结构传热过程。     

吸气式高超声速飞行器机体/推进一体化内外流非定常耦合方法

基于吸气式高超声速飞行器机体/推进一体化的气动布局设计方式，文章提出了一种内外流一体化流场的耦合求解方法，其中燃烧室内流场采用考虑有限速率化学反应动力学模型的一维非稳态方法求解，进气道和尾喷管外流场采用二维CFD软件计算，进气道与燃烧室在耦合界面处通过一维平均方法实现静温、静压和Mach数等参数传递.并分别以日本国家航空与航天实验室（NAL）的氢燃料燃烧室模型作为内流场验证算例，以某典型高超声速飞行器一体化模型作为内外耦合流场验证算例.研究结果表明:有限速率化学反应准一维方法能较为准确地模拟燃烧室内燃烧流场，提出的内外流场耦合方法能够有效地计算出内外流耦合效应，计算后体压力分布与理论值较接近.该方法可为超燃冲压发动机的性能快速分析和吸气式高超声速飞行器机体/推进一体化的初步分析设计提供重要参考.      

激光辐照高超声速舵面热力响应研究

建立了基于流固耦合的高超声速飞行器舵面结构在激光辐照下静气动弹性模型,流体控制方程为三维雷诺平均N-S方程,分别采用了中心格式和AUSM+up格式对粘性项和对流项进行空间离散,时间推进采用了高斯-塞德尔隐式推进方法,湍流粘性系数求解使用Menter SST模型。利用冯·卡门研究所高超风洞实验结果对模型进行了校核,预估了激光辐照对高超舵面的热力影响。结果表明,气动力/热计算模型与实验数据符合较好,能够准确模拟高超飞行器的热力参数,根据模型外推结果,激光在较低功率下加热高超飞行器舵面可能导致舵面材料弹性模量大大降低,继而发生弯曲发散而折断,高超飞行器可能因此发生气动失稳而坠毁。     

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

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

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

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

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

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

气动加热对高超声速飞行器激光毁伤效应影响

通过采用工程计算方法求解高超声速飞行器碳-碳复合材料分别在气动热、激光单独作用以及气动热/激光耦合作用下的热化学烧蚀。计算分析表明:激光单独作用下,碳-碳复合材料的烧蚀速率较小;随激光能量的增加,碳-碳复合材料的烧蚀速率增加;气动加热条件下激光对高超声速飞行器碳-碳复合材料的烧蚀毁伤效应会明显增强;沿弹道的气动加热累积效应对碳-碳复合材料气动热/激光耦合烧蚀作用不明显。     

低超声速飞行器瞬态热状况解耦算法

针对低超声速飞行器非稳态飞行条件下内外流固耦合一体化计算的复杂性,将飞行器外部流场的实时气动热转化为浮动的第三类边界条件进行解耦.以加速俯冲的超声速三维头锥体为例,分别采用浮动温差法和辐射平衡法提取表面对流换热系数进行解耦计算,并与直接耦合计算结果进行比较,验证两种解耦算法的可靠性.结果表明,将非稳态飞行过程离散为不同飞行状态点,通过提取对流换热系数解耦计算得到的不同状态点的锥体表面温度分布与直接耦合计算得到的结果吻合较好.两种解耦算法在计算效率方面均要优于耦合计算方法;在外界气动环境发生剧烈变化的过程中,最大相对误差均不超过2%.     

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

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

高超声速飞行器热气动弹性仿真计算方法综述

高超声速飞行器遭遇的热气动弹性问题具有危险性和复杂性，而型号设计又对计算环境下的热气动弹性仿真提出了效率和精度的双重要求.针对此问题，首先回顾了热气动弹性分析中气动力、气动热、结构热传导、热结构动力学等子学科的仿真计算方法，包括工程方法、数值模拟及降阶建模方法.其中，降阶建模方法有望解决目前计算效率和精度的矛盾，是调研的重点.之后，从时空耦合的角度总结了热气动弹性分析的多场耦合计算架构，包括单向、双向空间耦合策略，以及全耦合、松耦合、时域-频域耦合等时间耦合策略.在结论中，提出了基于子学科降阶模型进行多场耦合仿真的思路，并对热气动弹性分析的发展趋势进行了展望.      

高超声速飞行器热管冷却前缘结构一体化建模分析

针对高超声速飞行器工作时前缘恶劣的气动加热环境, 为了保证飞行器前缘的尖锐外形, 提出内嵌高温热管前缘结构. 针对热管内部复杂流动与换热情况, 对内嵌高温热管前缘结构进行一体化建模, 将模型的核心部件液态金属热管工作状况的计算与实验进行对比以验证模型的可靠性. 本文还分析了给定工况下内嵌高温热管前缘结构的热防护效果, 其中壁面最高温度下降了11.6%, 最低温度上升了8%, 高温区和低温区均被封闭在前缘外层区域, 内层温度更加均匀, 实现了热流由高温区向低温区的转移, 削弱了高温区的热负荷. 本文还分析了接触热阻对热管冷却前缘结构效果的影响. 关键词： 热管 前缘 疏导式热防护 气动热     

A theoretical study on the unsteady aerothermodynamics for attached flow models

The principle of the unsteady aerothermodynamics was theoretically investigated for the attached flow. Firstly, two simplified models with analytic solutions to the N-S equations were selected for the research, namely the compressible unsteady flows on the infinite flat plate with both time-varying wall velocity and time-varying wall temperature boundary conditions. The unsteady temperature field and the unsteady wall heat flux (heat flow) were analytically solved for the second model. Then, the interaction characteristic of the unsteady temperature field and the unsteady velocity field in the simplified models and the effects of the interaction on the transient wall heat transfer were studied by these two analytic solutions. The unsteady heat flux, which is governed by the energy equation, is directly related to the unsteady compression work and viscous dissipation which originates from the velocity field governed by the momentum equation. The main parameters and their roles in how the unsteady velocity field affects the unsteady heat flux were discussed for the simplified models. Lastly, the similarity criteria of the unsteady aerothermodynamics were derived based on the compressible boundary layer equations. Along with the Strouhal number Stu, the unsteadiness criterion of the velocity field, StT number, the unsteadiness criterion of the temperature field was proposed for the first time. Different from the traditional method used in unsteady aerodynamics which measures the flow unsteadiness only by the Stu number, present results show that the flow unsteadiness in unsteady aerothermodynamics should be comprehensively estimated by comparing the relative magnitudes of the temperature field unsteadiness criterion StT number with the coefficients of other terms in the dimensionless energy equation.     

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

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

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

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

Investigation on Rijke pipe＇s acoustic characteristics by numerical simulation： modeling the pulsing flow field coupled the inner of pipe with its outer space

Based on the results of fluid dynamics, heat transfer and acoustics, a Computational Fluid Dynamics （CFD） method was utilized to study the acoustic characteristics and self-excited pulsation mechanism inside a Rijke pipe. To avoid settling the irrational boundary conditions of the finite-amplitude standing wave in the Rijke thermo-acoustic system, the simulation modeling in the flow field, which coupled the inner of pipe with its outer space, was carried out to replace the traditional way in form of internal flow field numerical investigations. A hypothesis for heat source in energy equation including the relationship on unsteady heat of air around heat source, oscillation pressure and oscillation velocity was presented. To reflect the essence of Rijke pipe, simulation on self-excited oscillation was conducted by means of its own pulsation of pressure, velocity and temperature. This method can make the convergence process steady and effectively avoid divergence. The physical phenomenon of the self-excited Rijke pipe was analyzed. Moreover, the mechanisms on the Rijke pipe＇s self-excited oscillation were explained. Based on this method, comparative researches on the acoustic characteristic of the Rijke pipe with different size and different shape of nozzle were performed. The simulation results agreed with the experimental data satisfactorily. The results show that this numerical simulation can be used to study the sound pressure of nozzle for the engineering application of Rijke pipes.     

类前缘防热层流场与热响应耦合计算研究

本文在以往对类前缘防热层热响应计算分析基础上,进一步研究实现了外流场高超音速NS方程数值计算表面气动加热与防热层结构热响应的耦合计算,这对于常用的非耦合计算方法来说是一进步,也为进一步开展外流场/结构热响应/热应力全耦合一体化计算研究和防热层表面吹气强化传热问题的场协同研究打下了基础。     

碳氢燃料流动换热与裂解反应的建模及仿真

碳氢燃料的化学热沉的释放特性不仅与燃料自身的化学性质有关,还受燃料温度和驻留时间等流动换热因素的影响,是高超声速热防护技术中的重要研究课题。考虑流动换热对裂解反应的影响及裂解反应所引起的组分变化对流动换热过程的影响,本文提出了一个一维非稳态的模型来研究碳氢燃料的裂解反应与流动换热过程的耦合特性.并对不同的温度、工作压力...     

Numerical investigation of thermal protection of hypersonic flying vehicles

The problem of spatial flow around a hypersonic flying vehicle is considered for trajectories with different attack angles for flight through air with chemical equilibrium. The conjugate problem statement gives solutions for gas state in the boundary layer, thermal regime of streamlined body (made of different composites), and the rate of mass loss for heat protecting material. Physical processes in the condensed phase of carbon-containing coatings have complex nature: processes of heating, pyrolysis, heterogeneous oxidation, and sublimation. These processes result in protective material destruction. It was shown that using different materials for passive protection can be beneficial in reduction of the surface temperature, characteristics of thermochemical degradation, and this allows a control over heat and mass transfer for a flying body.     

HEAT TRANSFER IN HIGH-PRESSURE COMPRESSOR GAS TURBINE INTERNAL AIR SYSTEMS: A ROTATING DISC-CONE CAVITY WITH AXIAL THROUGHFLOW

This article reports on heat transfer measurements made on a rotating test rig representing the internal disc-cone cavity of a gas turbine high-pressure (H.P.) compressor stack. Tests were carried out for a range of flow rates and rotational speeds at engine representative nondimensional conditions. The rig also had a central drive shaft, which could rotate in the same direction as the discs, contrarotate relative to the discs, or remain static. Measurements of heat transfer were obtained from a conduction solution method using measured surface temperatures as boundary conditions. Results from the outer surface of the cone are in reasonable agreement with theoretical predictions for the heat transfer from a free cone in turbulent flow. The heat transfer measurements from the inner surface of the cone reveal two regimes of heat transfer: one governed by rotation, the other by action of the throughflow. In the rotationally dominated regime, the heat transfer from the inner surface of the cone is higher for a co-rotating shaft than for either a static or contra-rotating shaft. In the throughflow-dominated regime the heat transfer shows little consistent dependence on the direction of shaft rotation. Tests carried out at different values of surface-to-fluid temperature difference add support to the hypothesis that in the rotationally dominated regime the heat transfer occurs through a process of free convection, where the buoyancy force is induced by rotation. The heat transfer from the disc is significantly lower than that from the inner surface of the cone and more or less insensitive to the sense of shaft rotation. The disc average Nusselt numbers show similar behavior to those from the inner surface of the cone and suggest that the disc heat transfer too is governed either by rotationally induced buoyancy or by the axial throughflow.