Numerical investigation of hypersonic step-flows

Hypersonic aerospace vehicles are exposed to extreme flight conditions with heavy contour loads during their mission. Especially at ridges and sharp corners, the wall heat flux and pressure may cause serious damage to the body. Sometimes, the surface material cannot resist the high loading and fails completely. In this work the laminar hypersonic flow over forward and backward facing steps is investigated by CFD techniques and the results are compared with experimental data. The selected flow conditions correspond to cold hypersonic flow according to the availability of experimental data. The Navier-Stokes equations in the high temperature gas approximation of a thermally perfect gas in local equilibrium serve as the model for the physical problem. A multiblock finite-volume method is used to discretize consistently all spatial derivatives appearing in the balance equations. A second order in space Godunov-type method is utilized for the non-diffusive part of the governing equations whereas centered differences are used for the diffusive part. Time integration is performed by a second order implicit scheme. In each time step, the resulting nonlinear system of equations is solved by Newton's method employing a relaxation scheme based on conjugate gradients for the linear equation system. The results obtained permit a close insight into the physics of the flow problems under consideration and by this provide valuable information for construction concepts of hypersonic vehicles. Besides a careful comparison of the numerical results with experimental data, numerical aspects like the grid influence are addressed. Received 9 November 1998 / Accepted 2 December 1999     

A hybrid DSMC/Navier–Stokes frame to solve mixed rarefied/nonrarefied hypersonic flows over nano‐plate and micro‐cylinder

We extend a hybrid DSMC/Navier–Stokes (NS) approach to unify the DSMC and the NS simulators in one framework capable of solving the mixed non‐equilibrium and near‐equilibrium flow regions efficiently. Furthermore, we use a one‐way state‐based coupling (Dirichlet–Dirichlet boundary‐condition coupling) to transfer the required information from the continuum region to the rarefied one. The current hybrid DSMC–NS frame is applied to the hypersonic flows over nanoflat plate and microcylinder cases. The achieved solutions are compared with the pure DSMC and NS solutions. The results show that the current hybrid approach predicts the surface heat transfer rate and shear stress magnitudes very accurately. Some important conclusions can be drawn from this study. For example, although the shock wave region would be a non‐equilibrium region, it is not necessary to use a pure DSMC simulator to solve it entirely. This is important when the researchers wish to predict the surface properties such as velocity slip, temperature jump, wall heat flux rate, and friction drag magnitudes accurately. Our investigation showed that our hybrid solution time would be at least 40% (for the flat plate) and 35% (for the cylinder) of the time that must be spent by a pure DSMC solver to attain the same accuracy.Copyright © 2013 John Wiley & Sons, Ltd.     

高超音速湍流分离表面热流率的脉动特性

本文给出高超音速湍流分离不稳定特性的实验研究结果。试验条件是:自由流马赫数为 7.8,单位长度雷诺数为 3.5×10~7/米。分离流场由有限展长前向台阶产生,并用有高空间分辨率和快速响应的一列平齐安装的铂膜电阻温度计和多通道系统测量其表面热流率脉动。信号的条件采样分析结果表明:分离激波的根部由一束压缩波构成,流向展长约二分之一来流边界层厚度,在边界层外汇聚成单一主激波。这种激波结构极其不稳定,出现大尺度运动,流向运动的尺度约为分离激波上游影响区域长度的22％。激波振荡频率为一宽频带,主要集中在 1～3 千赫。在分离激波运动区域,热流脉动呈间歇性,在无扰动和激波扰动间跳跃。可以认为这种间歇性是分离激波系统大尺度振荡的结果。在激波运动区域的下游为分离区,流体继续压缩,热流脉动无间歇。     

Vortex structure generation on the frontal surface of a cylinder set transversely in a hypersonic flow

The problem of formation of spatially periodic structures on the frontal surface of a cylindrically blunted body set transversely in a hypersonic flow is studied. Within the framework of the model adopted, a possible mechanism of vortex structure generation on the frontal surface of the blunt body is proposed and confirmed by calculations; in this mechanism, the curved bow shock produces a vortex flow, while in its turn the vortex, which persists under weak dissipation, acts on the shock thus maintaining its curved shape. It is shown that the spatially periodic mode of hypersonic flow past a cylinder can exist in the case of a uniform incident flow and under homogeneous boundary conditions on the body surface.     

论粘性激波层理论的理论基础及其进一步发展

本文将Davis的量级分析方法改用匹配渐近展开方法,作为一级近似推导出了高超音速化学反应粘性激波层方程。证明了粘性激波层方程是NS方程在匹配渐近意义上的一级近似方程。进一步讨论了这一方程的基本假设条件。本文首次推导出的二级近似方程,是对Davis粘性激波层方程的修正。这种修正可以提高数值解的精度,有助于对问题获得更全面的了解,对进一步发展与完善高超音速钝头体绕流问题的数值求解方法起一定的作用。     

Magnetohydrodynamic interaction in hypersonic air flow past a blunt body

Hypersonic MHD air flow past a blunt body in the presence of an external magnetic field is considered. The MHD effect on the flow consists in a significant increase in the shock wave stand-off from the body surface and a significant reduction in the heat flux to the wall (up to 50%). It is shown that even in the presence of a strong Hall effect the intensity of the magnetohydrodynamic interaction in the plasma behind the shock wave remains at a high level commensurable with the ideal case of the absence of a Hall effect.     

Hypersonic Flow Past the Spherical Nose of a Body in the Presence of a Magnetic Field

The hypersonic flow past the nose of a spherical body containing current sources generating a magnetic field is investigated theoretically and numerically. The magnetohydrodynamic (MHD) flow is analyzed on the basis of the complete system of Navier-Stokes equations containing the force and thermal MHD terms and the electrodynamic equations. Local and integral thermal and aerodynamic characteristics of the body are found. It is shown that the presence of a magnetic field makes it possible to reduce the heat flow to the body in the neighborhood of the stagnation point by several times. However, in this case the total body drag increases.     

返回舱三维高超声速绕流及近尾迹流场数值模拟

为使返回舱安全、稳定、可靠地飞行，准确地计算其周围的复杂绕流流场，对飞船的初步设计是十分必要的。用Ｈａｒｔｅｎ－Ｙｅｅ的二阶迎风ＴＶＤ有限差分格式求解薄层Ｎ－Ｓ方程，模拟了返回舱三维高超声速流场，Ｍ＿∞＝７．３５，Ｒｅ＿∞＝７．５×１０￣５，α＝１０°、２０°。给出了详细的绕流结构，不同攻角、不同子午面上的物面压力分布与Ｍｏｓｅｌｅｙ和ｗｅｌｌｓ的实验数据进行了比较，符合较好。通过分析表明，在一定的攻角下，倒锥体上低压区压力的计算精度，对力矩系数及压心位置仍有明显的影响。     

Heating characteristics of blunt swept fin-induced shock wave turbulent boundary layer interaction

An experimental study was conducted on shock wave turbulent boundary layer interactions caused by a blunt swept fin-plate configuration at Mach numbers of 5.0, 7.8, 9.9 for a Reynolds number range of (1.0∼4.7)×10⁷/m. Detailed heat transfer and pressure distributions were measured at fin deflection angles of up to 30° for a sweepback angle of 67.6°. Surface oil flow patterns and liquid crystal thermograms as well as schlieren pictures of fin shock shape were taken. The study shows that the flow was separated at deflection of 10° and secondary separation were detected at deflection of ϑ≥20°. The heat transfer and pressure distributions on flat plate showed an extensive plateau region followed by a distinct dip and local peak close to the fin foot. Measurements of the plateau pressure and heat transfer were in good agreement with existing prediction methods, but pressure and heating peak measurements atM≥6 were significantly lower than predicted by the simple prediction techniques at lower Mach numbers. The project supported by China Academy of Launch Vehicle Technology     

三维N-S方程计算隐式有限体积法

利用NND有限差分格式，发展了一种新的完全隐式的有限体积数值方法，以求解与时间相关的N－S方程．对通过单元体界面的无粘流和粘性流通量均作隐式处理．对绕流钝锥体和不同攻角的气动辅助实验飞行器的高超声速粘性流和化学反应流获得了定常数值解．对流加热率和流场电子密度的计算值与实验数据进行了比较，符合较好，证实了本方法的精确性．