高超声速非平衡流动-辐射特性数值模拟研究

飞行器高超声速飞行过程中所承受对流加热和辐射加热可能具有相当的量级,因此合理准确预测气动加热需要将二者进行综合考虑.文章发展了具有非玻尔兹曼电子能级分布和振动能级分布的高温空气碰撞辐射模型,并耦合一维激波后流动方程计算不同飞行条件下激波后的非平衡流动特性,采用逐线辐射输运模型计算获得激波后非平衡辐射特性、辐射强度和辐射输运通量,深入比较分析了不同飞行高度和马赫数对非平衡流动和辐射输运过程的影响.计算结果表明对于高空高马赫飞行条件,其波后流动存在显著的热力学非平衡、化学非平衡和能级非平衡特征,在近激波区域高振动能级和原子高束缚电子激发态明显低于玻尔兹曼分布.在高空高马赫条件下真空紫外辐射占据主导地位,主要是由高能原子束缚-束缚跃迁造成的.随着高度和马赫数的下降,激波层内气体解离和电离程度降低,原子辐射贡献下降,分子辐射贡献增加,导致红外、可见光和紫外波段的辐射输运增强,真空紫外辐射输运过程减弱.     

基于格子Boltzmann方法的非平衡气动加热算法

提出了一种新的Euler与Boltzmann方程结合的气动加热数值计算方法。该方法首先采用CFD软件中的FLUENT来获取三维复杂外形边界层外缘气动参数,通过可压缩格子Boltzmann方法(LBM)的微观非平衡粒子的流动和碰撞来模拟边界层中气体热流动过程,解决了工程方法和N-S方法的适用范围不足和计算效率低下的问题。在每个LBM的流动与碰撞的过程中通过更新密度分布函数来使得宏观的速度和温度计算得以耦合。同时,对热LBM模型还做了湍流比拟和改进来扩展其应用范围。最后,通过一组高超声速气动加热试验数据证明了该结合方案的正确性和有效性,并讨论了该结合方案较其它常用方法的优缺点。     

k-ω SST湍流模式三维激波分离流动修正

“激波?边界层分离”是航空气动领域的典型湍流非平衡流动问题, 准确模拟激波分离对于跨声速飞行器气动性能评估和优化设计具有重要意义. 然而传统涡黏性湍流模式中涡黏性系数的定义方式并不适用于非平衡流动, k-ω SST湍流模式为此引入的Bradshaw假设在应用于三维强逆压梯度和较大分离流动时反而限制了雷诺应力的生成, 导致包括k-ω SST在内的常用涡黏性湍流模式均无法对此类流动进行准确模拟. 同时, 现有的非线性雷诺应力本构关系也并不能有效提高模拟精度. 为此, 针对k-ω SST模式分别提出了基于Bradshaw假设和基于长度尺度的两种激波分离流动修正方法. 前者通过提高Bradshaw常数的方式放宽了对雷诺应力生成的限制, 后者则从湍流长度尺度概念出发, 利用混合长度理论、湍动能生成/耗散之比和一种新定义的长度尺度之比构造了ω方程耗散项修正函数, 提高了模式在三维激波分离流动中的建模长度尺度. 两种方法对ONERA M6机翼跨声速大攻角流动均能得到较雷诺应力模式更好的模拟结果. 进一步的雷诺应力分析表明, 三维激波分离流动中“主雷诺应力分量”的概念不再成立, 各雷诺应力分量大小接近. 网格收敛性分析、对其他攻角状态的验证以及湍流平板边界层壁面律验证进一步确认了所提出的两种修正方法的合理性、有效性和通用性.      

高温激波管化学非平衡流动数值模拟研究

为了预测氢氧定容燃烧驱动的高温激波管性能,需要准确分析激波管非定常化学非平衡流动过程.本文在破膜前的驱动段定容燃烧以及破膜后的化学非平衡流动数值模拟中,引入双时间步长方法,发展高温激波管化学非平衡流动数值模拟方法,该方法在时间上具有二阶精度.计算结果与目前存在的激波管流动解析解以及零维化学反应系统的数值解进行了比较,吻合较好.对于典型高温激波管状态,采用有限体积方法离散准一维流动Euler控制方程,并通过将流动过程和化学反应动力学过程耦合求解,获得了激波管内部的化学非平衡流动特征.     

风力机翼型动态失速等离子体流动控制数值研究

针对动态失速引起的风力机翼型气动性能恶化的问题,本文基于动网格和滑移网格技术, 开展了大涡模拟数值计算研究,探索了非定常脉冲等离子体的动态流动控制机理. 结果表明,等离子体气动激励能够有效控制翼型动态失速, 改善平均和瞬态气动力,减小力矩负峰值和迟滞环面积. 压力分布在等离子体施加范围内出现了负压"凸起",上翼面吸力峰值明显增大.脉冲频率和占空比这两个非定常控制参数对流动控制影响显著,无因次脉冲频率为1.5时等离子体控制效果较好,占空比为0.8时即可接近连续工作模式下的气动收益. 翼型深失速状态,等离子体促使流动分离位置明显向后缘移动, 抵抗了大尺度动态失速涡的发生,分离涡结构破碎耗散、重新附着, 涡流影响范围减小; 浅失速状态,等离子体激励具有较强的剪切层操纵能力, 诱导了翼型边界层提前转捩,促进了与主流的动量掺混. 等离子体气动激励诱导出前缘附近贴体翼面"涡簇",起到了虚拟气动外形的作用.不同尺度、频域的动态涡结构与等离子体气动激励的非线性、强耦合作用导致了气动力/力矩的谐波振荡.      

基于滑移速度壁模型的复杂边界湍流大涡模拟

采用滑移速度壁模型实现了浸入边界方法与壁模型相结合的大涡模拟.本文首先分别采用平衡层模型和非平衡壁模型对周期山状流进行数值模拟,以考查在壁模型中考虑切向压力梯度的作用.数值结果表明,流场的压力对本文所采用的壁模型形式并不敏感,但是考虑切向压力梯度可以显著改进壁面摩擦力的计算结果,并且能够准确的预测强压力梯度区以及分离区内的流动平均统计特性.不考虑压力梯度效应的平衡层模型显著低估了壁面摩擦力的分布,同时无法准确预测分离区内的平均速度剖面.非平衡模型的修正项正比于切向压力梯度和壁面法向距离,因此在强压力梯度区或者网格较粗时,计算得到的平均压力和摩擦力分布以及流动的低阶统计量均与参考的实验和计算结果吻合.在此基础上,通过回转体绕流的大涡模拟考查了该方法用于模拟高雷诺数壁湍流的适用性,非平衡壁模型可以准确地捕捉流动的物理结构并较准确地预测其水动力学特性.结果表明,将浸入边界方法与非平衡滑移速度壁模型相结合的大涡模拟,有望成为数值模拟复杂边界高雷诺数壁湍流的工具.      

NUMERICAL STUDY ON DYNAMIC STALL FLOW CONTROL FOR WIND TURBINE AIRFOIL USING PLASMA ACTUATOR 1)

针对动态失速引起的风力机翼型气动性能恶化的问题,本文基于动网格和滑移网格技术, 开展了大涡模拟数值计算研究,探索了非定常脉冲等离子体的动态流动控制机理. 结果表明,等离子体气动激励能够有效控制翼型动态失速, 改善平均和瞬态气动力,减小力矩负峰值和迟滞环面积. 压力分布在等离子体施加范围内出现了负压"凸起",上翼面吸力峰值明显增大.脉冲频率和占空比这两个非定常控制参数对流动控制影响显著,无因次脉冲频率为1.5时等离子体控制效果较好,占空比为0.8时即可接近连续工作模式下的气动收益. 翼型深失速状态,等离子体促使流动分离位置明显向后缘移动, 抵抗了大尺度动态失速涡的发生,分离涡结构破碎耗散、重新附着, 涡流影响范围减小; 浅失速状态,等离子体激励具有较强的剪切层操纵能力, 诱导了翼型边界层提前转捩,促进了与主流的动量掺混. 等离子体气动激励诱导出前缘附近贴体翼面"涡簇",起到了虚拟气动外形的作用.不同尺度、频域的动态涡结构与等离子体气动激励的非线性、强耦合作用导致了气动力/力矩的谐波振荡.     

HLLE++格式在高马赫数空腔流动模拟中的应用

空腔流动存在剪切层运动、涡脱落与破裂,以及激波与激波、激波与剪切层、激波与膨胀波和激波/涡/剪切层相互干扰等现象,流动非常复杂,特别是高马赫数(M>2)时,剪切层和激波更强,激波与激波干扰更严重,对数值格式的要求更高,既需要格式耗散小,对分离涡等有很高的模拟精度,又需要格式在激波附近具有较大的耗散,可以很好地捕捉激波,防止非物理解的出现。Roe和HLLC等近似Riemann解格式在高马赫数强激波处可能会出现红玉现象,而HLLE++格式大大改善了这种缺陷,在捕捉高超声速激波时避免了红玉现象的发生,同时还保持在光滑区域的低数值耗散特性。本文在结构网格下HLLE++格式的基础上,通过改进激波探测的求解,建立了基于非结构混合网格的HLLE++计算方法,通过无粘斜坡算例,验证了HLLE++格式模拟高马赫数流动的能力,并应用于高马赫数空腔流动的数值模拟,开展了网格和湍流模型影响研究,验证了方法模拟高马赫数空腔流动的可靠性和有效性。     

基于滑移速度壁模型的复杂边界湍流大涡模拟

采用滑移速度壁模型实现了浸入边界方法与壁模型相结合的大涡模拟.本文首先分别采用平衡层模型和非平衡壁模型对周期山状流进行数值模拟,以考查在壁模型中考虑切向压力梯度的作用.数值结果表明,流场的压力对本文所采用的壁模型形式并不敏感,但是考虑切向压力梯度可以显著改进壁面摩擦力的计算结果,并且能够准确的预测强压力梯度区以及分离区内的流动平均统计特性.不考虑压力梯度效应的平衡层模型显著低估了壁面摩擦力的分布,同时无法准确预测分离区内的平均速度剖面.非平衡模型的修正项正比于切向压力梯度和壁面法向距离,因此在强压力梯度区或者网格较粗时,计算得到的平均压力和摩擦力分布以及流动的低阶统计量均与参考的实验和计算结果吻合.在此基础上,通过回转体绕流的大涡模拟考查了该方法用于模拟高雷诺数壁湍流的适用性,非平衡壁模型可以准确地捕捉流动的物理结构并较准确地预测其水动力学特性.结果表明,将浸入边界方法与非平衡滑移速度壁模型相结合的大涡模拟,有望成为数值模拟复杂边界高雷诺数壁湍流的工具.     

高温热化学非平衡气动热试验与仿真技术研究进展

临近空间新型飞行器向全空域、更高马赫数发展,面临的气动热环境会越发恶劣,高温流场气动热预测技术是该类飞行器发展的关键技术之一.高超声速气流通过激波压缩或黏性阻滞减速,分子动能转化为内能,产生了高温.高温引起体分子振动、电子激发,伴随离解、电离反应等一系列复杂气动物理现象,其流场气动热预测面临诸多挑战.文章对高温热化学非平衡气动热预测技术的发展情况进行了分析探讨.首先,阐述了国内外高温气动热地面试验技术的发展历程,重点介绍分析了气动热风洞试验设备的模拟能力及目前试验测试技术的研究水平;然后,调研和讨论了高温气动热数值模拟研究现状,分别从热化学模型、辐射输运和壁面催化/烧蚀等多个角度探讨了热化学非平衡流场气动热数值模拟规律;最后,对气动热预测技术的发展趋势进行了讨论,提出了高温气动热试验与仿真技术后续应重点解决的问题.     

DSMC (direct simulation Monte Carlo) analysis of microchannel flow with gas–liquid boundary

High heat capacity and constant operation temperature make a 2-phase heat remover tool promising for solving high heat dissipation problems in MEMS devices. However, microscale analysis of the flow with the conventional Navier–Stokes equation is inadequate, because the non-continuum effect is important when the characteristic dimension is comparable to the local mean free path. DSMC is a direct, particle-based numerical simulation method that uses no continuum assumption. In this paper, the gas–liquid boundary effects in microchannel flow are studied using this method. Modified DSMC code is used to simulate low-speed flow—under which viscous heating produces no significant temperature change—and MD results are incorporated into the DSMC boundary condition. Steady Couette flow simulation results show that the gas–liquid boundary affects the density distribution and the temperature dependence of the slip velocity. Unsteady simulation results show that mass transfer by diffusion is faster than momentum transfer by collision.     

Generalised Couette flow of two immiscible viscous fluids with heat transfer using Brinkman model

The present study deals with generalised Couette flow of two viscous, incompressible, immiscible fluids with heat transfer in presence of heat source through two straight parallel horizontal walls. The lower wall is bounded below, by a naturally permeable material of high porosity and the flow inside the porous medium is assumed to be of moderate permeability, modelled by Brinkman equation. The flow domain is divided into three zones to obtain analytical solutions of the momentum and energy equations. To link various flow regions, appropriate matching conditions have been used. The effects of permeability parameter, Reynolds number and viscous parameter on velocity field and the effects of Reynolds number, viscous parameter, permeability parameter, constant heat source and Brinkman number on temperature distribution in different zones are discussed graphically. The mass flow rate, skin-friction factor and rates of heat transfer at the upper boundary and porous interface are discussed with the help of tables.     

Heat transfer to flowing polymer melts

A theoretical and experimental study of heat transfer to polymer melts flowing through ducts is presented. The mathematical model provides for shear heating and expansion cooling effects for various wall boundary conditions.Experimental results, obtained in tube flow, show reproducible temperature and velocity profiles and the data confirm the predictions of the magnitude of the shear heating and expansion cooling effects.It is concluded that the method can be justified for predicting heat transfer in more complex geometries and some preliminary results are presented.     

Cross-flow heat transfer in fixed bed

Radial flow reactor operated at cross-flow heat transfer is focused for large scale methanol synthesis. The effects of operating conditions including the reactor inlet air temperature, the heating pipe temperature and the air flow rate on the cross-flow heat transfer were investigated and results show that the temperature profile of the area in front of the heating pipe is slightly affected by all the operating conditions. The main area whose temperature profile is influenced is located behind the heating pipe. The heat transfer direction is related to the direction of the flow. In order to obtain the basic parameters for radial flow reactor designing calculation, the dimensionless number group method was used for data fitting of the bed effective thermal conductivity and the wall heat transfer coefficient which were calculated by the mathematical model with the product of Reynolds number and Prandtl number. The comparison of experimental data and calculated values shows that the calculated values fit the experimental data satisfactorily and the formulas can be used for reactor designing calculation.     

Comparison of velocity and temperature measurements with simulations in a hypersonic wake flow

A hypersonic shock-tunnel flow around an axisymmetric model of a planetary entry probe is analyzed. Planar laser-induced fluorescence is applied to measure both the velocity and the rotational temperature everywhere in the central plane of the flow field. The experimental test case is compared to simulations using the direct simulation Monte Carlo (DSMC) method. While the Mach 9.7 flow at a freestream Reynolds number based on the model diameter of 35,000 is chemically frozen, effects of thermal non-equilibrium and localized rarefaction cannot be neglected. DSMC and measurements agree well within the outer wake, but disagree close to the centerline, where in particular the measured velocity is higher than values predicted by the simulations. The experimental results indicated a shorter recirculation region and increased local fluctuations in the free shear layer upstream of the wake recompression shock when compared to the simulation. These effects are attributed to incipient transition, which is not observed in the simulations, as the simulations did not model the effects of freestream fluctuations. Furthermore, measured and simulated vorticities are compared with theoretical predictions.     

Heat and mass transfer to a reacting particle in a gas flow for arbitrary temperature dependence of the transport coefficients

A study is made of simultaneous heat and mass transfer to a reacting particle of any shape in a translational (and shear) flow of a viscous heat conducting compressible gas for which the thermal conductivity and diffusion coefficient, and also the specific heat depend on the temperature. The first two terms of the asymptotic expansion with respect to the small Reynolds number are obtained for the mean Sherwood and Nusselt numbers. The case of a power-law temperature dependence of the gas viscosity is considered in detail.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 111–119, January–February, 1984.     

Direct numerical simulation of turbulent heat transfer on the Reynolds analogy over irregular rough surfaces

The effect of rough surface topography on heat and momentum transfer is studied by direct numerical simulations of turbulent heat transfer over uniformly heated three-dimensional irregular rough surfaces, where the effective slope and skewness values are systematically varied while maintaining a fixed root-mean-square roughness. The friction Reynolds number is fixed at 450, and the temperature is treated as a passive scalar with a Prandtl number of unity. Both the skin friction coefficient and Stanton number are enhanced by the wall roughness. However, the Reynolds analogy factor for the rough surface is lower than that for the smooth surface. The semi-analytical expression for the Reynolds analogy factor suggests that the Reynolds analogy factor is related to the skin friction coefficient and the difference between the temperature and velocity roughness functions, and the Reynolds analogy factor for the present rough surfaces is found to be predicted solely based on the equivalent sand-grain roughness. This suggests that the relationship between the Reynolds analogy factor and the equivalent sand-grain roughness is not affected by the effective slope and skewness values. Analysis of the heat and momentum transfer mechanisms based on the spatial- and time-averaged equations suggests that two factors decrease the Reynolds analogy factor. One is the increased effective Prandtl number within the rough surface in which the momentum diffusivity due to the combined effects of turbulence and dispersion is larger than the corresponding thermal diffusivity. The other is the significant increase in the pressure drag force term above the mean roughness height.     

Energy Estimate of the Critical Reynolds Numbers in a Compressible Couette Flow. Effect of Bulk Viscosity

A variational problem of determining the critical Reynolds number of the laminar-turbulent transition is numerically solved within the framework of the nonlinear energy theory of stability of compressible flows. Stability of various modes in the Couette flow of a compressible gas is estimated by the method of collocations. It is demonstrated that the minimum critical Reynolds numbers in the range of the ratio of the bulk viscosity ηb to the shear viscosity η, which is realistic for diatomic gases, are reached for modes of streamwise disturbances. The critical Reynolds numbers increase as the bulk viscosity is increased in the interval η_b = 0-2η, with the maximum increase in the limit being approximately 30%.     

Oscillating plate temperature effects on a flow past an infinite vertical porous plate with constant suction and embedded in a porous medium

 An approximate solution to the problem of flow of a viscous incompressible dissipative fluid past an infinite vertical porous plate embedded in a porous medium is presented here. The plate temperature is assumed to be oscillating about a constant mean temperature. Mean velocity and mean temperature, the transient velocity and temperature profiles are shown graphically. The mean skin-friction and the mean rate of heat transfer are also shown graphically. The expressions for the amplitude and the phase of the skin-friction and the rate of heat transfer are derived and their numerical values are listed in Tables. The effects of different parameters governing the unsteady flow are discussed. Received on 23 November 1998     

Linear stability of supersonic Couette flow of a molecular gas under the conditions of viscous stratification and excitation of the vibrational mode

The linear stability of viscous two-dimensional perturbations in the supersonic plane Couette flow of perfect and vibrationally excited gases is investigated. In both cases an alternative is considered so that the transport coefficients were taken either constant or dependent on the static flow temperature. The Sutherland model is used to take the temperature dependence of the shear viscosity into account. It is shown that “viscous” stratification increases considerably the flow stability as compared with the case of constant viscosity. At the same time, the simple constant viscosity model conserves all characteristic features of the development of viscous perturbations in the Sutherland model. The dissipation effect of excitation of the vibrational mode is conserved in taking the temperature dependence of the transport coefficients into account. For both models the corresponding increase in the critical Reynolds number is of approximately 12%.