Investigation of hypersonic flows through a cavity with sweepback angle in near space using the DSMC method

Near space has been paid more and more attentionin recent years due to its military application value.However,flow characteristics of some fundamental configurations(e.g.,the cavity)in near space have rarely been investigated due to rarefied gas effects,which make the numerical simulation methods based on continuous flow hypothesis lose validity.In this work,the direct simulation Monte Carlo(DSMC),one of the most successful particle simulation methods in treating rarefied gas dynamics,is employed to explore flow characteristics of a hypersonic cavity with sweepback angle in near space by considering a variety of cases,such as the cavity at a wide range of altitudes 20-60 km,the cavity at freestream Mach numbers of 6-20,and the cavity with a sweepback angle of 30°-90°.By analyzing the simulation results,flow characteristics are obtained and meanwhile some interesting phenomena are also found.The primary recirculation region,which occupies the most area of the cavity,causes pressure and temperature stratification due to rotational motion of fluid inside it,whereas the pressure and temperature in the secondary recirculation region,which is a small vortex and locates at the lower left corner of the cavity,change slightly due to low-speed movement of fluid inside it.With the increase of altitude,both the primary and secondary recirculation regions contract greatly and it causes them to separate.A notable finding is that rotation direction of the secondary recirculation region would be reversed at a higher altitude.The overall effect of increasing the Mach number is that the velocity,pressure,and temperature within the cavity increase uniformly.The maximum pressure nearby the trailing edge of the cavity decreases rapidly as the sweepback angle increases,whereas the influence of sweepback angle on velocity distribution and maximum temperature within the cavity is slight.

NON-STATIONARY EFFECTS IN HYPERSONIC NONUNIFORM DUSTY-GAS FLOW PAST A BLUNT BODY

In the framework of the two-fluid model, a hypersonic flow of a nonuniform dusty gas with low inertial (non-depositing) particles around a blunt body is considered. The particle mass concentration is assumed to be small, so that the effect of particles on the carrier phase is significant only inside the boundary layer where the particles accumulate. Stepshaped and harmonic nonuniformities of the particle concentration ahead of the bow shock wave are considered and the corresponding nonstationary distributions of the particle concentration in the shock layer are studied. On the basis of numerical study of nonstationary two-phase boundary layer equations derived by the matched asymptotic expansion method, the effects of free-stream particle concentration nonuniformities on the thermal flux, and the friction coefficient in the neighborhood of stagnation point are investigated, in particular, the most “dangerous” nonuniformity periods are found.

Stable Runge-Kutta discontinuous Galerkin solver for hypersonic rarefied gaseous flow based on 2D Boltzmann kinetic model equations

A stable high-order Runge-Kutta discontinuous Galerkin(RKDG) scheme that strictly preserves positivity of the solution is designed to solve the Boltzmann kinetic equation with model collision integrals. Stability is kept by accuracy of velocity discretization, conservative calculation of the discrete collision relaxation term, and a limiter. By keeping the time step smaller than the local mean collision time and forcing positivity values of velocity distribution functions on certain points, the limiter can preserve positivity of solutions to the cell average velocity distribution functions. Verification is performed with a normal shock wave at a Mach number 2.05, a hypersonic flow about a two-dimensional(2D) cylinder at Mach numbers 6.0 and 12.0, and an unsteady shock tube flow. The results show that, the scheme is stable and accurate to capture shock structures in steady and unsteady hypersonic rarefied gaseous flows. Compared with two widely used limiters, the current limiter has the advantage of easy implementation and ability of minimizing the influence of accuracy of the original RKDG method.

Prediction of hypersonic boundary layer transition on sharp wedge flow considering variable specific heat

When the air temperature reaches 600 K or higher, vibration is excited. The specific heat is not a constant but a function of temperature. Under this condition, the transition position of hypersonic sharp wedge boundary layer is predicted by using the improved eN method considering variable specific heat. The transition positions with different Mach numbers of oncoming flow, half wedge angles, and wall conditions are computed condition, the nearer to the Mach number The results show that for the same oncoming flow condition and wall transition positions of hypersonic sharp wedge boundary layer move much leading edge than those of the flat plate. The greater the oncoming flow the closer the transition position to the leading edge.

SEPARATED FLOW AROUND BENT-NOSE BICONIC IN HYPERSONIC FLOW

This article investigates the separated flow around the bent-nose biconic, in which the second ordered TVD(Total Variation Diminishing) method is applied. The governing cqua-tion is a thin layer Navier-Stokes equation. The surface pressure distribution at the leeward side is a little higher in the vicinity of the symmetrical plane. This seems to be caused by the development of the vortex at the leeward side due to the separated flow. The location of a separation line and the size are in very good agreement with the experiment. Also the secondary separation is captured very clearly.

AN ANALYTIC SOLUTION ON HYPERSONIC FLOW OVER AN ARBTTRARY SLENDER BODY WITH NEAR POWER-LAW PROFILE

On the basis of a self-similar solution as well as of the assumption of the 'Transserse Motion',a general linear theory on hypersonic flow over a general slender body is set up in this paper By means of this theory, the problem concerned can he put into a universal system of O.D Eqs .which can be integrated manerically in adyance

An Empirical Method for Prediction of Hypersonic Rarefied Flow-Field Structure

Numerical simulations are presented about the effects of gas rarefaction on hypersonic flow field.Due to the extremely difficult experiment,limited wind-tunnel conditions and high cost,most problems in rarefied flow regime are investigated through numerical methods,in which the direct simulation Monte-Carlo(DSMC)method is widely adopted.And the unstructured DSMC method is employed here.Flows around a vertical plate at a given velocity 7 500 m/s are simulated.For gas rarefaction is judged by the free-stream Knudsen number(Kn),two vital factors are considered:molecular number density and the plate′s length.Cases in which Kn varies from 0.035 to13.36 are simulated.Flow characters in the whole rarefied regime are described,and flow-field structure affected by Knis analyzed.Then,the dimensionless position D*of a certain velocity in the stagnation line is chosen as the marker of flow field to measure its variation.Through flow-field tracing and least-square numerical method analyzing,it is proved that hypersonic rarefied flow field expands outward linearly with the increase of 1/2Kn.An empirical method is proposed,which can be used for the prediction of the hypersonic flow-field structure at a given inflow velocity,especially the shock wave position.

NATURE OF THE SURFACE HEAT TRANSFER FLUCTUATION IN A HYPERSONIC SEPARATED TURBULENT FLOW

This paper presents the results of an experimental study of the unsteady nature of a hypersonic sepa- rated turbulent flow.The nominal test conditions were a freestream Mach number of 7.8 and a unit Reynolds number of 3.5x10^7/m.The separated flow was generated using finite span forward facing steps.An array of flush mounted high spatial resolution and fast response platinum film resistance thermometers was used to make mul- ti-channel measurements of the fluctuating surface heat trtansfer within the separated flow.Conditional sampling ana- lysis of the signals shows that the root of separation shock wave consists of a series of compression wave extending over a streamwise length about one half of the incoming boundary layer thickness.The compression waves con- verge into a single leading shock beyond the boundary layer.The shock structure is unsteady and undergoes large-scale motion in the streamwise direction.The length scale of the motion is about 22 percent of the upstream influence length of the separation shock wave.There exists a wide band of frequency of oscillations of the shock system.Most of the frequencies are in the range of 1-3 kHz.The heat transfer fluctuates intermittently between the undisturbed level and the disturbed level within the range of motion of the separation shock wave.This inter mittent phenomenon is considered as the consequence of the large-scale shock system oscillations.Downstream of the range of shock wave motion there is a separated region where the flow experiences continuous compression and no intermittency phenomenon is observed.

A reduced-order model for fast predicting ionized flows of hypersonic vehicles along flight trajectory

An improved Reduced-Order Model(ROM)is proposed based on a flow-solution preprocessing operation and a fast sampling strategy to efficiently and accurately predict ionized hypersonic flows.This ROM is generated in low-dimensional space by performing the Proper Orthogonal Decomposition(POD)on snapshots and is coupled with the Radial Basis Function(RBF)to achieve fast prediction speed.However,due to the disparate scales in the ionized flow field,the conventional ROM usually generates spurious negative errors.Here,this issue is addressed by performing flow-solution preprocessing in logarithmic space to improve the conventional ROM.Then,extra orthogonal polynomials are introduced in the RBF interpolation to achieve additional improvement of the prediction accuracy.In addition,to construct high-efficiency snapshots,a trajectory-constrained adaptive sampling strategy based on convex hull optimization is developed.To evaluate the performance of the proposed fast prediction method,two hypersonic vehicles with classic configurations,i.e.a wave-rider and a reentry capsule,are used to validate the proposed method.Both two cases show that the proposed fast prediction method has high accuracy near the vehicle surface and the free-stream region where the flow field is smooth.Compared with the conventional ROM prediction,the prediction results are significantly improved by the proposed method around the discontinuities,e.g.the shock wave and the ionized layer.As a result,the proposed fast prediction method reduces the error of the conventional ROM by at least 45%,with a speedup of approximately 2.0×105compared to the Computational Fluid Dynamic(CFD)simulations.These test cases demonstrate that the method developed here is efficient and accurate for predicting ionized hypersonic flows.

Experimental study on shock interaction control of double wedge in high-enthalpy hypersonic flow subject to plasma synthetic jet

The hypersonic shock-shock interaction flow field at double-wedge geometries controlled by plasma synthetic jet actuator is experimentally studied in a Ma = 8 high-enthalpy shock tunnel with the purpose of exploring a novel technique for reducing surface heat flux in a real flight environment. The results demonstrate that increasing the discharge energy is advantageous in eliminating the shock wave, shifting the shock wave interaction point, and shortening the control response time. The oblique shock wave can be completely removed when the actuator's discharge energy grows from 0.4 J to 11.5 J, and the displacement of the shock wave interaction point increases by 124.56%, while the controlled response time is shortened by 30 μs. Besides, the reduction in diameter of the jet exit is firstly proved to have a negative impact on energy deposition in a working environment with incoming flow, which reduces the discharge energy and hence decreases the control effect. The shock wave control response time lengthens when the jet exits away from the second wedge. Along with comparing the change in wall heat flux at the second wedge over time, the control effect of plasma synthetic jet actuator with and without inflation is also analyzed. When plasma synthetic jet works in inflatable mode, both the ability to eliminate shock waves and the shifting effect of the shock wave interaction point are increased significantly, and the wall heat flux is also reduced.

Study on the characteristics of interaction flowfields induced by supersonic jet on a revolution body

The paper focuses on the triple jets interaction with a hypersonic external flow on a revolution body. The experimental model is a ogive-cylinder body with three supersonic nozzles, which are aligned along the flow direction. The freestream Mach numbers are 5 and 6. The spatial and surface flow characteristics are illustrated by the schlieren photographs and the typical pressure distribution. The results show that there are multi-wave system, separation, reattachment, multi-peak pressure, high-pressure and low-pressure zone boundaries obvious distinction in tri-jets interference flowfield. The present paper also analyzes how do the pressure ratio, the angle of attack, and Mach number effect on tri-jets interaction characteristics.

Hypersonic Shock Wave/Boundary Layer Interactions by a Third-Order Optimized Symmetric WENO Scheme

A novel third-order optimized symmetric weighted essentially non-oscillatory(WENO-OS3)scheme is used to simulate the hypersonic shock wave/boundary layer interactions.Firstly,the scheme is presented with the achievement of low dissipation in smooth region and robust shock-capturing capabilities in discontinuities.The Maxwell slip boundary conditions are employed to consider the rarefied effect near the surface.Secondly,several validating tests are given to show the good resolution of the WENO-OS3 scheme and the feasibility of the Maxwell slip boundary conditions.Finally,hypersonic flows around the hollow cylinder truncated flare(HCTF)and the25°/55°sharp double cone are studied.Discussions are made on the characteristics of the hypersonic shock wave/boundary layer interactions with and without the consideration of the slip effect.The results indicate that the scheme has a good capability in predicting heat transfer with a high resolution for describing fluid structures.With the slip boundary conditions,the separation region at the corner is smaller and the prediction is more accurate than that with no-slip boundary conditions.

Gao's interacting shear flows( ISF) theory and its inferences and their applications

Gao＇s viscous/in-viscid interacting shear flows （ISF） theory, proposed by professor Gao Zhi in Institute of Mechanics, China Academy of Science, and its inferences and their applications in computational fluid dynamics （CFD） are reviewed and some subjects worthy to be studied are pro- posed in this paper. The flow-field and motion law of ISF, mathematics definition of strong viscous shear layer flow in ISF, ISF equations, wall-surface compatibility criteria （Gao＇s criteria ）, space scale variety law of strong viscous shear layer reveals flow mechanism and local space small scale triggered by strong interaction that cause some abnormal severe local pneumatic heating phenomenon in hypersonic flow. Gao＇s ISF theory was used in near wall flow, free ISF flow simulation and design of computing grids, Gao＇s wall-surface criteria were used to verify calculation reliability and accuracy of near wall flows, ISF theory approximate analytical result of shock waves-boundary layer interac- tion and ISF equations were used to obtain the numerical exact solution of local area flow （ such as stationary point flow）. Some new subjects, such as, improving near-wall turbulent models according to the turbulent flow simulation satisfying the wall-criteria and illustrating relation between grid-con- vergence based on the wall criteria and other convergence tactics, are suggested. The necessity of applying Gao＇s ISF theory and wall criteria is revealed. Difficulties and importance of hypersonic vis- cous/in-viscid interaction phenomenon were also emphasized.

Fluctuating Pressure Inside/Outside the Flow Separation Region in High Speed Flowfield

An experimental study was conducted on the interactions of shock wave/turbulence or laminar boundary layer caused by fin-type protuberance, as the lack of detailed understanding of fluctuating pressure loads inside and outside the laminar or turbulence boundary layer separation region in hypersonic flow. The changes of fluctuating pressure in separation region were focused on in this paper. The study shows that the existence of fin changes flowfiled on the plate significantly. The laminar boundary layer separation occurs earlier and the separation region is more extensive. Similar flow is observed between a couple of measurement points outside the laminar separation region. However, there are significant differences between the flow inside and outside the separation region. The level of fluctuating pressure of laminar boundary layer is smaller than that in turbulent case. Even so, in laminar case, the peak fluctuating pressure still reaches a high level. Therefore, the structural influence （damage and/or early fatigue） of fluctuating pressure loads caused by the laminar boundary layer separation should not be ignored.

Physical criterion study on forward stagnation point heat flux CFD computations at hypersonic speeds

In order to evaluate uncertainties in computational fluid dynamics （CFD） computations of the stagnation point heat flux, a physical criterion is developed. Based on a quasi-one-dimensional hypothesis along the stagnation line, a new stagnation flow model is applied to obtain the governing equations of the flow near the stagnation point at hypersonic speeds. From the above equations, the compatibility relations are given at the stagnation point and along the stagnation line, which consist of the physical criterion for checking the accuracy in the stagnation point heat flux computations. The verification of the criterion is made with various numerical results.

On transition of type V interaction in double-wedge flow with non-equilibrium effects

The transition between regular reflection （RR） and Mach reflection （MR） of type V shock-shock interaction on a double-wedge geometry with high temperature non-equilibrium effects is investigated by extended shock-polar method and numerical simulation. First, the critical angles of transition from detachment criterion and yon Neumann criterion are determined by the extended shock-polar method considering the non-equilibrium effects. Then wave patterns and the transition process are numerically obtained. Results of the critical transition angles from shock-polar calculation and numerical simulation show evident disagreement, indicating transition mechanism between RR and MR of type V interaction is changed. By comparing with the frozen counterpart, it is also found that non-equilibrium effects lead to a larger critical wedge angle and a larger hysteresis interval.

Effects of mesh resolution on hypersonic heating prediction

Aeroheating prediction is a challenging and critical problem for the design and optimization of hypersonic vehicles.One challenge is that the solution of the Navier-Stokes equations strongly depends on the computational mesh.In this letter,the effect of mesh resolution on heat flux prediction is studied.It is found that mesh-independent solutions can be obtained using fine mesh,whose accuracy is confirmed by results from kinetic particle simulation.It is analyzed that mesh-induced numerical error comes mainly from the flux calculation in the boundary layer whereas the temperature gradient on the surface can be evaluated using a wall function.Numerical schemes having strong capability of boundary layer capture are therefore recommended for hypersonic heating prediction.

Numerical simulation of hypersonic thermochemical nonequilibrium flows using nonlinear coupled constitutive relations

To predict aeroheating performance of hypersonic vehicles accurately in thermochemical nonequilibrium flows accompanied by rarefaction effect,a Nonlinear Coupled Constitutive Relations(NCCR)model coupled with Gupta’s chemical models and Park’s two-temperature model is firstly proposed in this paper.Three typical cases are intensively investigated for further validation,including hypersonic flows over a two-dimensional cylinder,a RAM-C II flight vehicle and a type HTV-2 flight vehicle.The results predicted by NCCR solution,such as heat flux coefficient and electron number densities,are in better agreement with those of direct simulation Monte Carlo or flight data than Navier-Stokes equations,especially in the extremely nonequilibrium regions,which indicates the potential of the newly-developed solution to capture both thermochemical and rarefied nonequilibrium effects.The comparisons between the present solver and NCCR model without a two-temperature model are also conducted to demonstrate the significance of vibrational energy source term in the accurate simulation of high-Mach flows.

Prediction of hypersonic boundary layer transition with variable specific heat on plane flow

When an aircraft flies at a hypersonic speed,the temperature of gas inner boundary layer near the wall is so high that the specific heat is no longer a constant but dependent upon the temperature.It is necessary to consider its effect on transition location.In this paper,the transition locations of hypersonic plane boundary layer are predicted with the improved e N method,and the results of the specific heat dependent upon temperature are compared with those of constant specific heat.The flow parameters are taken as those corresponding to the condition at a height of 40 km and the Mach numbers of oncoming flow are 6,7,and 8,respectively.It is found that the transition locations calculated by the variable specific heat are closer to the leading edge than those by the constant specific heat.The deviations in most cases are around 30 percent.All the results prove that the real gas effect should be taken into consideration when one predicts transition location for hypersonic flow.Whether the first or second mode wave determines the transition location relies on the oncoming flow Mach number and the wall condition.

A hybrid CFD/characteristics method for fast characterization of hypersonic blunt forebody/inlet flow

A hybrid CFD/characteristic method(CCM) was proposed for fast design and evaluation of hypersonic inlet flow with nose bluntness, which targets the combined advantages of CFD and method of characteristics. Both the accuracy and efficiency of the developed CCM were verified reliably, and it was well demonstrated for the external surfaces design of a hypersonic forebody/inlet with nose bluntness. With the help of CCM method, effects of nose bluntness on forebody shock shapes and the flowfield qualities which dominate inlet performance were examined and analyzed on the two-dimensional and axisymmetric configurations. The results showed that blunt effects of a wedge forebody are more substantial than that of related cone cases. For a conical forebody with a properly blunted nose, a recovery of the shock front back to that of corresponding sharp nose is exhibited, accompanied with a gradually fading out of entropy layer effects. Consequently a simplification is thought to be reasonable for an axisymmetric inlet with a proper compression angle, and a blunt nose of limited radius can be idealized as a sharp nose, as the spillage and flow variations at the entrance are negligible, even though the nose scale increases to 10% cowl lip radius. Whereas for two-dimensional inlets, the blunt effects are substantial since not only the inlet capturing/starting capabilities, but also the flow uniformities are obviously degraded.