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

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

基于Boltzmann模型方程的气体运动论HPF并行算法

从修正的BGK-Blotzmann模型方程出发，引入离散速度坐标技术对气体分子速度分量进行离散降维，基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散差分格式，发展直接求解气体分子速度分布函数的气体运动论有限差分数值格式，提出了一套能有效模拟各流域三维绕流问题的气体运动论统一算法，在分析研究气体运动论数值算法内在并行度的基础上，开展各流域三维绕流问题统一算法的HPF（高性能FORTRAN）并行化程度设计，建立一套能有效模拟各流域复杂外形体绕流的HPF并行算法软件，并进行了不同Knudsen(克努森）数下三维球体绕流及类“神舟号”返回舱外形体绕流的初步数值实验，将计算结果与过渡区有关实验数据及各流域气体绕流现象进行比较分析，证实了发展的气体运动论HPF并行算法在求解稀薄流到连续流不同流域复杂绕流问题方面的可行性。     

激波对超音速平板发汗冷却的影响

本文分析了激波影响下,不同激波强度,冷却气体和注入率对发汗冷却效果的影响。马赫数3的超音速主流遇到流道内的楔形激波发生器产生斜激波入射到多孔平板表面。楔块楔角φ分别为0°、4°、8°、12°模拟不同的入射激波强度,冷却介质分别为空气、甲烷与氢气。计算结果表明,激波使多孔区域出口表面静压上升,冷却流体流出受阻而破坏冷却效果,且冷却效率随激波强度增强而下降,随冷却气体分子量增大,冷却效率下降幅度减小;但在激波强度较强时,激波在多孔表面形成逆压梯度,迫使冷却流体流向入射点上游区域而使该区域冷却效果得到恢复。提高冷却剂注入率可以减弱激波对冷却效率的破坏作用,但使壁面温度不均匀性增加。     

高精度非定常激波装配法

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

微管道气体流动的蒙特卡洛直接模拟

采用蒙特卡洛直接模拟(DSMC)方法,数值模拟了微管道中压力驱动的气体流动,结果表明固壁边界存在速度滑移,稀薄气体效应明显;整个流场温度变化很小,流动马赫数很小,密度、压力流向变化非常大而横向几乎不变;可压缩性导致压力随流向的非线性分布,但这种效应随Knudsen数增大而减弱.     

基于第二粘性的Navier-Stokes方程组求解正激波结构

为考察气体第二粘性（体积粘性）对正激波内部流动的影响机制，数值求解含第二粘性的一维Navier-Stokes方程组.结果表明：第二粘性对激波内部的密度、热流和能量分布等物理量具有抹平效应，导致热流和熵流的峰值减小、激波厚度增加，体积粘性耗散的增加使得一部分机械能转化为内能；考虑第二粘性所计算的密度分布和激波厚度大为改善，与实验数据吻合较好；当马赫数为1.2≤Ma≤10，激波内部的Knudsen数满足0.12≤Kn≤0.4，对于马赫数Ma≤4.0的激波内部流动，考虑第二粘性的连续流Navier-Stokes方程组能够准确地模拟正激波结构.     

考虑转动能的一维/二维Boltzmann-Rykov模型方程数值算法

研究考虑转动能的Boltzmann-Rykov模型方程,基于转动自由度对气体分子速度分布函数矩积分,引入约化速度分布函数,应用离散速度坐标法与数值积分技术,将气体运动论模型方程化为在离散速度坐标点处关于三个约化速度分布函数的联立方程组.应用拓展计算流体力学有限差分方法,数值计算考虑转动自由度的双原子气体一维、二维Boltzmann模型方程,得到高、低Knudsen数一维激波管内流动和二维竖直平板绕流问题的流场,分析验证考虑转动能的Boltzmann-Rykov模型方程全流域统一算法求解一维/二维气体流动问题的可靠性.结果表明,气体稀薄程度与分子内自由度对流场具有较大影响,且Knudsen数较高的稀薄气体流动呈现严重的非平衡流动特点.     

高重复频率激光能量沉积减小超声速波阻的数值研究

 采用数值模拟的方法研究了激光重复频率、点火位置及来流马赫数等参数对激光能量沉积减小超声速钝头体波阻的影响。数值模拟结果表明,由于激光能量的沉积产生的低密度区与弓形激波相互作用,在钝头体前形成了类似虚拟尖锥的回流区,使原弓形激波逐渐向阻力较小的斜激波转变。阻力随着频率的增加而减小,当频率增加到200 kHz时,阻力减小到约为原来的17%,能量效率的最大值出现在频率为50 kHz处。说明控制参数的选择对减阻性能起着关键的作用。     

Laval喷管内激波位置的计算及制冷性能分析

明确激波在Laval喷管内的发生位置,能够为设计及优化Laval喷管线型提供理论依据,从而提高Laval喷管的制冷性能和整个天然气超音速分离器的分离效率。文中对Laval喷管进行了结构设计,对不同背压条件下Laval喷管内激波位置进行了理论分析与计算,并利用FLUENT软件进行了数值模拟,通过研究不同背压条件下Laval喷管内气体马赫数、压力和温度分布,对喷管制冷性能进行了对比分析。结果表明:保持Laval喷管入口压力不变,随出口背压增大,激波位置逐渐从喷管出口向喷管入口方向移动,气体受到激波的影响,在喷管所能达到的最大马赫数不断降低,所能产生的最低温度不断上升,喷管的制冷性能越差;理论计算与数值模拟结果基本一致,数值模拟验证了理论计算的正确性。     

不同来流马赫数下激光能量形成高超声速进气道虚拟唇口的特性

用基于三维非定常可压缩雷诺平均Navier-Stokes方程的有限体积法计算了马赫数低于设计值6时一种高超声速进气道的性能参数,发现其性能存在明显下降。为提升进气道性能,将功率为15kW的激光能量注入进气道固体唇口前的流场中,形成虚拟唇口,马赫数为4.5,5.0和5.5时,计算得到来流捕获率分别提高了34%,20.6%和15.6%。绘制了不同马赫数下来流捕获率达到峰值时的流场压强云图,说明了虚拟唇口的特性及形成机制。结果表明:来流马赫数越低,来流捕获率越小,但相对于无能量注入时的来流捕获率的提升程度越明显;在不同来流马赫数条件下,通过改变激光能量引致的激波结构和位置,可达到最优状态,即激波与进气道前缘斜激波相交后的透射波打在进气道肩部位置的状态。     

Gasdynamic design of a two-dimensional supersonic inlet with the increased flow rate factor

A two-dimensional inlet of external compression with the increased flow rate factor at high supersonic velocities is constructed by the method of gasdynamic design. Its feature is that a flow with the initial oblique shock wave and the subsequent centered isentropic compression wave is formed over the external compression ramp of the inlet. These waves interact with one another so that a resulting stronger oblique shock wave and a velocity discontinuity arise in front of the entrance to the inlet internal duct. An example of an inlet configuration with the design flow regime corresponding to the Mach number M_d = 7 is considered. The characteristics of this inlet were obtained in the range of the free-stream Mach numbers M = 4–7 with the use of a Navier—Stokes code for turbulent flow. They are compared with characteristics of an equivalent conventional shocked inlet. As computations have shown, the inlet with the isentropic compression wave has much higher values of flow rate factor φ at Mach numbers M < M_d. So, for example, at M = 4 the value φ ≈ 0.72 for it is by 33 % higher in comparison with φ ≈ 0.54 for the equivalent shocked inlet.     

Acoustic disturbance from gas non-uniformities convected through a nozzle

The non-steady flow generated by convection of gas containing non-uniform temperature regions or “entropy spots” through a nozzle is examined analytically as a source of acoustic disturbance. The first portion of the investigation treats the “compact nozzle”, the case where all wave lengths are much longer than the nozzle. Strengths of transmitted and reflected one-dimensional waves are given for supersonic and subsonic nozzles and for one configuration of supersonic nozzle with normal shock at the outlet. In addition to a wave reflected from the nozzle inlet, the supersonic nozzle discharges two waves, one facing upstream and the other facing downstream. For reasonable values of the nozzle inlet Mach number, the pressure amplitude of each wave increases directly as the discharge Mach number.The acoustic perturbations from a supercritical nozzle of finite length, in which the undisturbed gas velocity increases linearly through the nozzle, are analyzed for several inlet and discharge Mach number values and over a wide frequency range. The results which agree with the compact analysis for low frequency, deviate considerably as the frequency rises, achieving pressure fluctuation levels of several times the compact values. It is shown that this result originates in a phase shift between the two waves emitted downstream and that the pressure fluctuations for moderate frequencies may be approximated from the compact analysis with an appropriate phase shift.In all cases, the pressure fluctuations caused by a 2% fluctuation in absolute inlet temperature are large enough to require consideration in acoustic analysis of nozzles or turbine blade channels.     

弱激波冲击无膜重气柱和气帘界面的实验研究

研制了一个改进的激波管设备,对马赫数为1.2的弱激波冲击作用下空气中SF6气柱和气帘界面的演变过程进行了初步的实验研究。通过设计激波管实验段、烟雾发生器、气体箱、进气吸气系统和激波管尾段,控制混合气体中SF6的峰值浓度和初始气流速度,建立了稳定、可重复的无膜气柱和气帘初始界面形成技术。利用高速摄影技术,在水平面内观测了气柱和气帘的初始界面图像,沿垂直方向观测了界面RM(Richtmyer-Meshkov)不稳定性的演变过程。气柱演变图像显示了典型的对涡结构,气帘演变图像显示了早期的多蘑菇形结构和后期的相邻波长干扰效应。图像后处理表明,气柱的高度和宽度、气帘的宽度均随时间单调增加,且宽度比高度增加快得多。从二维涡量动力学方程出发,对图像中涡的演变过程进行了初步解释。     

Structure and parameters of the shock layer formed when a supersonic underexpanded jet interacts with a counterpropagating hypersonic flow in the transition regime

The method of molecular dynamics (the Bird system) has been used to mathematically model a planar, strongly underexpanded supersonic jet that encounters a hypersonic flow of rarefied gas. Particular attention is paid to the structure and parameters of the shock layer close to the plane of symmetry. The results of calculations are presented for currents of a monatomic gas, simulating argon, with a Mach number of the external flow of M_∞=5.48, a Mach number at the nozzle edge of M_a=1, a ratio of the density at the nozzle edge to the density of the unperturbed flow equal to 130, and various stagnation temperatures of the external flow and of the jet. The evolution of the structure and the parameters of the shock layer as the Knudsen number Kn_∞ varies from 0.02 to 0.35 is considered. The results are compared with the data calculated for the shock layer when argon flows around thermally insulated cylinders. The main features and regularities of the relaxation of the translational degrees of freedom of the gas for external and jet flows are considered. Data are presented on the form of the distribution function over velocities and its evolution as gas moves through the shock layers. Zh. Tekh. Fiz. 68, 13–18 (July 1998)     

Local effect of electric and magnetic fields on the position of an attached shock in a supersonic diffuser

Effective ways for controlling shock wave configurations by means of external actions are sought. One such way is a local effect of electric and magnetic fields. In this paper, the local effect of external fields is implemented by current localization in a limited region of a diffuser. The experiment is carried out in a diffuser providing the complete internal compression of the gas with a Mach number at the inlet M=4.3. As a working medium, a xenon plasma is used. The plasma flow is formed in a shock tube equipped with an accelerating nozzle. Two ways of current localization are tested. In the first one, the diffuser inlet is a short channel of Faraday generator type. In this case, the ponderomotive force basically decelerates or accelerates the flow depending on the direction of the electric field. In the second way, the current flows through a narrow near-wall region between adjacent electrodes. In this case, the ponderomotive force compresses or expands the gas. In both cases, it is shown that the angle of an attached shock due to MHD interaction can be both decreased and increased. The central problem with the MHD control of shock waves is near-electrode and near-wall phenomena.     

球头高超音速热绕流的研究

本文研究了加热（冷却）对流动的影响。通过计算分析了在不同来流马赫数和加热量下，加热对球头高超音速绕流的影响规律。研究表明，加热对流动有着显著的影响，将使流动出现“热斥冷吸”、表面压力和激波距离增大、加热与冷却影响的非线性、加热使超音速区内马赫数下降趋向于１等现象。     

A unified gas kinetic scheme with moving mesh and velocity space adaptation

There is great difficulty for direct Boltzmann solvers to simulate high Knudsen number flow due to the severe steep slope and high concentration of the gas distribution function in a local particle velocity space. Local mesh adaptation becomes necessary in order to make the Boltzmann solver to be a practical tool in aerospace applications. The present research improves the unified gas-kinetic scheme (UGKS) in the following two aspects. First, the UGKS is extended in a physical space with moving mesh. This technique is important to study a freely flying object in a rarefied environment. Second, the adaptive quadtree method in the particle velocity space is implemented in the UGKS. Due to the new improvements in the discretization of a gas distribution function in the six dimensional phase space, the adaptive unified gas kinetic scheme (AUGKS) is able to deal with a wide range of flow problems under extreme flying conditions, such as the whole unsteady flying process of an object from a highly rarefied to a continuum flow regime. After validating the scheme, the capability of AUGKS is demonstrated in the following two challenge test cases. The first case is about the free movement of an ellipse flying at initial Mach number 5 in a rarefied flow at different Knudsen numbers. The force on the ellipse and the unsteady trajectory of the ellipse movement are fully captured. The gas distribution function around the ellipse is analyzed. The second case is about the study of unsteady flight of a nozzle under a bursting process of the compressed gas expanding into a rarefied environment. Due to the strong expansion wave and the huge density difference between interior and exterior regions around the nozzle, the particle distribution function changes dramatically in the particle velocity space. The use of an adaptive velocity space in the AUGKS becomes necessary to simulate such a flow and to control the computational cost to a tolerable level. The second test is a challenge problem for any existing rarefied flow solver.     

测量喷气Z箍缩负载的气流马赫数

 利用双压力传感器的方法测量了喷气Z箍缩负载的气流马赫数。测量得到的气流马赫数最高可达到4.25，它和定常计算的结果4.8很接近，这表明可以用定常计算来估算喷气Z箍缩中超声速喷嘴产生的气流马赫数。     

Numerical study of shock wave interaction in steady flows of a viscous heat-conducting gas with a low ratio of specific heats

Specific features of shock wave interaction in a viscous heat-conducting gas with a low ratio of specific heats are numerically studied. The case of the Mach reflection of shock waves with a negative angle of the reflected wave with respect to the free-stream velocity vector is considered, and the influence of viscosity on the flow structure is analyzed. Various issues of nonuniqueness of the shock wave configuration for different Reynolds numbers are discussed. Depending on the initial conditions and Reynolds numbers, two different shock wave configurations may exist: regular configuration interacting with an expansion fan and Mach configuration. In the dual solution domain, a possibility of the transition from regular to the Mach reflection of shock waves is considered.     

喷嘴直径对喷射器性能的影响分析

利用二维轴对称,真实气体模型对喷射式制冷空调系统的喷射器进行CFD计算。真实气体模型占用较多的计算机资源,但它的计算结果更加真实。计算的目的在于得到在不同喷嘴直径条件下喷射器性能的变化和制冷系统性能的变化情况。在计算工况条件下,存在一最优的喷嘴直径Dn*(5mm),使得此时的喷射系数和系统COP最大。当Dn>Dn*时,喷射系数和COP随着喷嘴直径的增加而减少;当Dn相似文献