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

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

超声速层流/湍流压缩拐角流动结构的实验研究

在Ma=3.0的超声速风洞中, 分别对上游边界层为超声速层流和湍流, 压缩角度为25°和28°的压缩拐角流动进行了实验研究. 采用纳米粒子示踪平面激光散射(NPLS)技术获得了流场整体和局部区域的精细结构, 边界层、剪切层、分离激波、回流区和再附激波等典型结构清晰可见, 测量了超声速层流压缩拐角壁面的压力系数. 从时间平均的流场结构中测量出分离激波、再附激波的角度和再附后重新发展的边界层的增长情况, 通过分析时间相关的流场NPLS图像, 可以发现流场结构随时间的演化特性. 实验结果表明: 在25°的压缩角度下, 超声速层流压缩拐角流动发生了典型的分离, 边界层迅速增长失稳转捩, 并引起一道诱导激波, 流场中出现了K-H涡、剪切层和微弱压缩波结构, 而超声速湍流压缩拐角流动没有出现分离, 湍流边界层始终表现为附着状态; 在28° 的压缩角度下, 超声速层流压缩拐角流动进一步分离, 回流区范围明显扩大, 诱导激波、分离激波向上游移动, 再附激波向下游移动, 分离区流动结构复杂, 相比之下, 超声速湍流压缩拐角流动的回流区范围明显较小, 边界层增长缓慢, 流场中没有出现诱导激波、K-H涡和压缩波, 流动分离区域的结构也相对简单, 但分离激波的强度则明显更强. 关键词： 压缩拐角 层流 湍流 流动结构     

Nonlinear Modeling Study of Aerodynamic Characteristics of an X38-like Vehicle at Strong Viscous Interaction Regions

Strong viscous interaction and multiple flow regimes exist when vehicles fly at high altitude and high Mach number conditions. The Navier–Stokes(NS) solver is no longer applicable in the above situation. Instead, the direct simulation Monte Carlo (DSMC) method or Boltzmann model equation solvers are usually needed. However, they are computationally more expensive than the NS solver. Therefore, it is of great engineering value to establish the aerodynamic prediction model of vehicles at high altitude and high Mach number conditions. In this paper, the hypersonic aerodynamic characteristics of an X38-like vehicle in typical conditions from 70 km to 110 km are simulated using the unified gas kinetic scheme (UGKS), which is applicable for all flow regimes. The contributions of pressure and viscous stress on the force coefficients are analyzed. The viscous interaction parameters, Mach number, and angle of attack are used as independent variables, and the difference between the force coefficients calculated by UGKS and the Euler solver is used as a dependent variable to establish a nonlinear viscous interaction model between them in the range of 70–110 km. The evaluation of the model is completed using the correlation coefficient and the relative orthogonal distance. The conventional viscous interaction effect and rarefied effect are both taken into account in the model. The model can be used to quickly obtain the hypersonic aerodynamic characteristics of X38-like vehicle in a wide range, which is meaningful for engineering design.     

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

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

Numerical simulation of the thermal non-equilibrium flow-field characteristics of a hypersonic Apollo-like vehicle

In order to investigate the relationship between the flow-field parameters outside the vehicle and the altitude, this paper takes the Atmospheric Reentry Demonstrator (ARD) with an angle of attack of -20° as the research object and adopts a two-temperature model coupled with the shear-stress transport k-ω turbulence model to focus on the variation of flow-field parameters including flow-field pressure, Mach number and temperature with the reentry altitude. It is found that the flow-field high-pressure region and low-Mach region both appear in the shock layer near the head of the ARD, while the maximum pressure of the surface appears on the windward side of the ARD's head with a toroidal distribution, and the numerical magnitude is inversely proportional to the radius of the torus. With fluid through the shoulder of the ARD flow expansion plays a dominant role, the airflow velocity increases, the Mach number of the windward side of the rear cone increases and the flow-field pressure and surface pressure rapidly decrease. When the fluid passes through the shock layer, the translational-rotation temperature will increase before the vibration-electron temperature, there is a thermal non-equilibrium effect and the two temperatures will rapidly decrease again when approaching the surface of the ARD due to the existence of temperature gradient. At the same time, both the windward side of the shoulder and the back cover of the ARD suffer from a large thermal load and require thermal protection.     

无缓冲气COIL扩压器流场数值模拟

 以无缓冲气化学氧碘激光器(COIL)实验器件的数据计算得到的混合喷管出口参数平均值作为光腔入口条件,对6种不同构型的扩压器从光腔入口至扩压器出口的流场进行了数值模拟,得出了各流场参数分布;对不同构型扩压器的流场特点、总压恢复性能进行了分析;研究了扩压器出口背压对流场参数的影响。结果表明：对于主流无缓冲气的COIL,等截面扩压器具有较好的压力恢复性能;增大扩压器出口背压可以使扩压器的压力恢复性能提高,然而,较高的背压使激波串向光腔方向移动,从而使光腔流场受到干扰,影响光腔的光束质量。     

超音速等离子体炬的数值模拟

介绍了描述大气压下超音速等离子体炬内等离子体特性的磁流体力学模型，在二维近似下，对会聚-扩展型喷口等离子体炬进行了数值模拟，获得了等离子体炬内等离子体速度、温度、压力以及马赫数的分布.结果表明超音速等离子体炬内的流场特性可以分为亚音速、跨音速和超音速三个明显的区域. 关键词： 等离子体炬 磁流体力学 数值模拟     

Spatio-temporal evolution of cavity ignition in supersonic flow

Successful ignition in the recirculating flow of a scramjet flame holder can be highly dependent upon the location of energy deposition because of the spatial variation of fuel concentration and flow properties. The current work experimentally investigated ignition processes when energy was deposited (～100 mJ) via a spark discharge at four locations in the base of a cavity or by laser-induced breakdown in a Mach 2 flow with a stagnation temperature and pressure of 590 K and 483 kPa, respectively. The cavity was directly fueled with ethylene injection. The time dependent heat release was imaged at 40,000 frames per second and fuel concentration and distribution measurements were taken in the cavity prior to ignition. The average fuel concentration at the lean and rich ignition limits near the energy deposition locations measured 4.4–9.3% (Φ= 0.75 to 1.47). Energy deposition near the cavity step resulted in near immediate ignition kernel development and rapid achievement of self-sustained flame propagation in the front of the cavity, often faster than the bulk recirculation time of the cavity, leading to a spike in heat release. Energy deposition away from the cavity step region led to competition between local flow velocity, fuel concentration, and flame propagation rates. Ignition kernels formed along the floor of the cavity towards the closeout ramp and were rapidly advected towards the cavity step region before flame propagation could ensue. The fastest and most robust ignition events for all fueling cases showed rapid spanwise flame propagation near the cavity step.     

压缩湍流剪切层-空腔流体气动光学效应的计算机仿真

对亚音速机载光电系统由于剪切层-腔流体而导致的气动光学效应进行了数值模拟.剪切层内的压缩效应用对流马赫数0.8表示,腔体的几何形状为L/H=6.7,平均腔内温度为230 K.应用商用CFD软件对这个非稳定二维流体进行了仿真.结果表明Strehl比随着波长的增加而增大,但是模糊角却变得相对复杂.证明了波长效应并给出了优化的波长模型λ*=2πσt.     

Flow characteristics of supersonic gas passing through a circular micro-channel under different inflow conditions

Gas flow in a micro-channel usually has a high Knudsen number. The predominant predictive tool for such a microflow is the direct simulation Monte Carlo(DSMC) method, which is used in this paper to investigate primary flow properties of supersonic gas in a circular micro-channel for different inflow conditions, such as free stream at different altitudes, with different incoming Mach numbers, and with different angles of attack. Simulation results indicate that the altitude and free stream incoming Mach number have a significant effect on the whole micro-channel flow field, whereas the angle of attack mainly affects the entrance part of micro-channel flow field. The fundamental mechanism behind the simulation results is also presented. With the increase of altitude, thr free stream would be partly prevented from entering into micro-channel.Meanwhile, the gas flow in micro-channel is decelerated, and the increase in the angle of attack also decelerates the gas flow. In contrast, gas flow in micro-channel is accelerated as free stream incoming Mach number increases. A noteworthy finding is that the rarefaction effects can become very dominant when the free stream incoming Mach number is low. In other words, a free stream with a larger incoming velocity is able to reduce the influence of the rarefaction effects on gas flow in the micro-channel.     

Aerothermodynamic Analysis of a Reentry Brazilian Satellite

This work deals with a computational investigation on the small ballistic reentry Brazilian vehicle SAtélite de Reentrada Atmosférica (SARA). Hypersonic flows over the vehicle SARA at zero-degree angle of attack in chemical equilibrium and thermal nonequilibrium are modeled by the direct simulation Monte Carlo method, which has become the main technique for studying complex multidimensional rarefied flows, and which properly accounts for the nonequilibrium aspects of the flows. The emphasis of this paper is to examine the behavior of the primary properties during the high-altitude portion of SARA reentry. In this way, velocity, density, pressure, and temperature field are investigated for altitudes of 100, 95, 90, 85, and 80?km. In addition, comparisons based on geometry are made between axisymmetric and planar two-dimensional configurations. Some significant differences between these configurations were noted on the flowfield structure in the reentry trajectory. The analysis showed that the flow disturbances have different influence on velocity, density, pressure, and temperature along the stagnation streamline ahead of the capsule nose. It was found that the stagnation region is a thermally stressed zone. It was also found that the stagnation region is a zone of strong compression, high wall pressure. Wall pressure distributions are compared with those of available experimental data, and good agreement is found along the spherical nose for the altitude range investigated.     

光腔与扩压器的一体化优化数值模拟

采用数值模拟方法对化学氧碘化学激光器光腔通道、超声速扩压器一体化方案的优化展开研究,对扩压器的角度、构型、背压等参数对扩压性能的影响以及对光腔内流场的影响进行计算和分析。研究结果表明：传统的直接扩散型以及平直段+扩散段型的超声速扩压器,抵抗背压影响的能力较弱,且光腔出口处静压急剧升高,影响了光腔内的流场;通过在平直段+扩散段型的超声速扩压器的平直段部分,插入数片楔形体,可以将扩压器的工作背压提升33%以上,且可以有效地隔绝扩压器对光腔内流场的不利影响,从而使光腔下游的逆压梯度大大降低;同时,由于缩短了扩压器的长度,扩压器的总压损失明显降低,冷流状态下的总压恢复系数达到0.484。     

稀薄气体二维外部柱体绕流特性研究

本文采用直接模拟蒙特卡罗方法对稀薄气体二维外部柱体绕流问题进行了数值模拟。结果表明:外部绕流问题,在特定情况下会产生激波,激波的产生,不仅与气体的稀薄程度有关,还与来流马赫数有关。而气流与壁面之间的换热,随来流马赫数增加而增加,随气体稀薄程度增加而减小。     

Numerical investigation of flow over two sweepback wedges at M = 4 and 6

Results of numerical simulation are discussed: simulation was carried out for a configuration of two wedges with sweepback leading edges placed on a pre-compression ramp in a way that skewed surfaces of the wedges deflect the compressed flows in the opposite directions. It was demonstrated that this configuration produces a flow with irregular interaction in the plane of symmetry for shock waves produced by sweepback wedges. The shock waves formed by the skew wedges induce 3D boundary layer separations along sweepback leading edges of the wedges. Flows in the separation zones are directed toward the plane of symmetry of this configuration; they interact and produce in the central part a “swollen” zone of separation flow with a typical S-shaped profile of velocity. Simulation data was obtained for the free stream flow with Mach number M = 4 and 6 and based upon Navier—Stokes equations and k-ωSST turbulence model using FLUENT computation code. Inviscid flow described by Euler equations was considered as well.     

预混段长度对旋流燃烧器气固两相流动特性的影响

本文利用三维PDA测量了预混段减小的斗山巴布科克旋流燃烧器气、固两相流动特性,并与设计的斗山巴布科克旋流燃烧器PDA测量结果对比,得出:预混段减小的燃烧器在x/d=0.1,0.3截面存在回流区,并将设计燃烧器回流区的位置由一二次风之间推迟至外二次风区域;预混段减后的燃烧器与设计燃烧器相比,轴向,径向,切向平均速度在二次风区域产生的峰值远离燃烧器轴线;在x/d=0.1～0.5截面,预混段减小的燃烧器与设计燃烧器相比一次风区域颗粒体积流量的峰值变大,二次风区域颗粒体积流量的峰值变小,在燃烧器出口附近回流区内的颗粒体积流量明显降低。     

偏滤器运行模式对托卡马克边缘区等离子体平行流的影响

我们使用一维流体模型,根据在不同偏滤器运行模式下静压强沿着磁力线方向的分布变化,讨论了偏滤器运行模式对托卡马克边缘区等离子体平行流的影响.低再循环模式下,静压强从X点 (X-point)附近的刮削层区域开始明显下降,变化趋势与密度变化趋势一致;等离子体平行流的马赫数在偏滤器区域逐步变大,变化从平缓到迅速.高再循环模式下,静压强在靶板附近的区域迅速下降,在其他区域变化非常小;等离子体平行流的马赫数仅在靠近靶板附近的区域迅速变大,在其他区域变化平缓.在弱脱靶模式下的静压强变化与高再循环模式下类似,不过静压强在X-point附近的刮削层区域开始出现下降的趋势,导致等离子体平行流的马赫数在X-point处的值比在高再循环模式下大.强脱靶模式下,静压强在刮削层区域开始明显下降,在远离靶板的偏滤器区域,静压强迅速下降的地方,观察到高马赫数等离子体平行流.静压强迅速下降引起动压强迅速上升来维持总的压强守恒是在强脱靶状态 下产生高马赫数平行流的一种可能驱动机理.     

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.     

基于数值模拟的飞行员高速气流吹袭损伤分析

研究飞行员高速气流吹袭的损伤机理,应用尺度适应模型对人椅系统外流场进行数值模拟.在马赫数为0.6,迎角为-90°~90°条件下的数值计算结果与风洞试验结果吻合较好.对不同马赫数、不同迎角的无防护状态下的人椅系统进行数值模拟,研究发现高速气流吹袭下飞行员面颈部、胸腹部、上臂内侧、小腿迎风侧压力较大;增大人椅系统的迎角可以改善飞行员身体各部分的受力情况.     

一种宽速域乘波体的设计及气动特性研究

为了解决乘波体偏离设计条件下气动特性会恶化,特别在低速时,升力严重不足这个问题,提出了通过增大后掠角生成前缘涡,增加背风面的升力,以改善乘波体低速气动性能.首先使用Visual Basic编程语言,并通过CATIA软件二次开发技术,实现了锥导乘波体的参数化设计和自动生成.再通过控制圆锥角和流场长度这两个设计参数,获得了大后掠乘波体构型.最后,运用剪切应力输运（shear-stress-transport,SST）模型,计算了所得乘波体的气动特性,并分析了流场变化,发现乘波体在设计状态下激波能很好附着在前缘上,在小的正攻角下,乘波体可获得比设计状态更高的升阻比,满足巡航要求.运用k-ω模型计算了乘波体的低速气动特性,得到了不同攻角下升力、阻力和升阻比的变化规律.研究结果发现,乘波体在低速下产生了明显的涡结构,在合适攻角下,能产生数量可观的附加升力,提高了乘波体的水平起降性能.      

An experimental study of pressure fluctuations generated by an open shallow cavity performed using joint time-frequency techniques of data analysis

An experimental study of the flow over a single open shallow cavity is reported. Pressure fluctuations were measured using two cavity models at free-stream Mach numbers 0.3, 0.5, and 0.6. At certain conditions, resonance of acoustic waves inside cavity with disturbances in the shear layer over cavity was observed. Pressure fluctuation characteristics were analysed using Fourier, wavelet, and Hilbert — Huang transforms. Those methods have allowed us to accurately examine the transition of the flow over cavity into resonant state.