Xinliang Li

Direct numerical simulation of hypersonic boundary layer transition over a blunt cone with a small angle of attack

The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [Li et al., AIAA J. 46, 2899 (2008)], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 20°≤θ≤30° (θ=0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20°...

Direct Numerical Simulation Of Hypersonic Boundary-Layer Transition Over A Blunt Cone

Direct numerical simulation of transition How over a blunt cone with a freestream Mach number of 6, Reynolds number of 10,000 based on the nose radius, and a 1-deg angle of attack is performed by using a seventh-order weighted essentially nonoscillatory scheme for the convection terms of the Navier-Stokes equations, together with an eighth-order central finite difference scheme for the viscous terms. The wall blow-and-suction perturbations, including random perturbation and multifrequency perturbation, are used to trigger the transition. The maximum amplitude of the wall-normal velocity disturbance is set to 1% of the freestream velocity. The obtained transition locations on the cone surface agree well with each other far both cases. Transition onset is located at about 500 times the nose radius in the leeward section and 750 times the nose radius in the windward section. The frequency spectrum of velocity and pressure fluctuations at different streamwise locations are analyzed and compared with the linear stability theory. The second-mode disturbance wave is deemed to be the dominating disturbance because the growth rate of the second mode is much higher than the first mode. The reason why transition in the leeward section occurs earlier than that in the windward section is analyzed. It is not because of higher local growth rate of disturbance waves in the leeward section, but because the growth start location of the dominating second-mode wave in the leeward section is much earlier than that in the windward section.

Evolution of mixing width induced by general Rayleigh-Taylor instability

A theory determining the evolution of general Rayleigh-Taylor mixing fronts is established to reproduce firstly all of the documented experiments conducted for diverse acceleration histories and all density ratios. The theory is established in terms of the fundamental conservation and symmetry principles, with special consideration given to the symmetry breaking of the density fields occurring in actual flows. The results reveal the sensitivity/insensitivity of the evolution of a mixing front neighbouring light/heavy fluid to the degree of symmetry breaking, and also explain the distinct evolutions in two experiments with the same configurations.Turbulent mixing induced by Rayleigh-Taylor (RT) instability occurs ubiquitously in many natural phenomena and engineering applications. As the simplest and primary descriptor of the mixing process, the evolution of mixing width of the mixing zone plays a notable role in the flows. The flows generally involve complex varying acceleration histories and widely varying density ratios, two dominant factors affecting the evolution of mixing width. However, no satisfactory theory for predicting the evolution has yet been established. Here a theory determining the evolution of mixing width in general RT flows is established to reproduce, first, all of the documented experiments conducted for diverse (i.e., constant, impulsive, oscillating, decreasing, increasing, and complex) acceleration histories and all density ratios. The theory is established in terms of the conservation principle, with special consideration given to the asymmetry of the volume-averaged density fields occurring in actual flows. The results reveal the sensitivity or insensitivity of the evolution of a mixing front of a neighboring light or heavy fluid to the degree of asymmetry and thus explain the distinct evolutions in two experiments with the same configurations.

DNS of shock/boundary layer interaction flow in a supersonic compression ramp

A direct numerical simulation of the shock/turbulent boundary layer interaction flow in a supersonic 24-degree compression ramp is conducted with the free stream Mach number 2.9. The blow-and-suction disturbance in the upstream wall boundary is used to trigger the transition. Both the mean wall pressure and the velocity profiles agree with those of the experimental data, which validates the simulation. The turbulent kinetic energy budget in the separation region is analyzed. Results show that the turbulent production term increases fast in the separation region, while the turbulent dissipation term reaches its peak in the near-wall region. The turbulent transport term contributes to the balance of the turbulent conduction and turbulent dissipation. Based on the analysis of instantaneous pressure in the downstream region of the mean shock and that in the separation bubble, the authors suggest that the low frequency oscillation of the shock is not caused by the upstream turbulent disturbance, but rather the instability of separation bubble.

An investigation of electromagnetic wave propagation in plasma by shock tube

This paper presents the electromagnetic wave propagation characteristics in plasma and the change of the attenuation coefficients of the microwave with some parameters by using #800mm high temperature shock tube. In order to get the attenuation of the electromagnetic wave through the plasma behind shock wave, the microwave transmission is used to measure the relative change of the wave power. The experimental results reveal that the dependence of the attenuation coefficients on ne, v, w, and L are in good agreement with the theoretical predictions of electromagnetic wave propagations in the uniform infinite plasma when the plasma layer is thick (such as L/λ ≥ 10). We also measure the electromagnetic wave power attenuate value with and without the magnetic field.

Difference schemes on non-uniform mesh and their application

High order accurate schemes are needed to simulate the multi-scale complex flow fields to get fine structures in simulation of the complex flows with large gradient of fluid parameters near the wall, and schemes on non-uniform mesh are desirable for many CFD (computational fluid dynamics) workers. The construction methods of difference approximations and several difference approximations on non-uniform mesh are presented. The accuracy of the methods and the influence of stretch ratio of the neighbor mesh increment on accuracy are discussed. Some comments on these methods are given, and comparison of the accuracy of the results obtained by schemes based on both non-uniform mesh and coordinate transformation is made, and some numerical examples with non-uniform mesh are presented.

Scale-adaptive subgrid-scale modelling for large-eddy simulation of turbulent flows

The proportionality between the subgrid-scale (SGS) drain rate of kinetic energy and the viscous dissipation rate of the resolved motions is studied a priori by filtering a given fully resolved field and evaluating a generic form of the hypothesized energy spectrum. The ratio of the SGS drain to the resolved dissipation, on which a balance condition for the SGS dissipation across an arbitrary grid scale is founded, is shown to be independent of the turbulence Reynolds number, and can be described by a function in terms of the averaged mesh Reynolds number. Such a balance condition can serve as a physical constraint in the SGS modeling to account for the scale effects of the model coefficient(s). Scale-adaptive dynamic Smagorinsky-Lilly model and mixed nonlinear model are formulated for large-eddy simulation of transitional and/or turbulent flows in such a way that the constraint is satisfied. The newly proposed scale-adaptive dynamic SGS models are validated in simulations of homogeneous isotropic turbule...

Resolution-optimised nonlinear scheme for secondary derivatives

A 5-point-stencil optimised nonlinear scheme with spectral-like resolution within the whole wave number range for secondary derivatives is devised. The proposed scheme can compensate for the dissipation deficiency of traditional linear schemes and suppress the spurious energy accumulation that occurs at high wave numbers, both of which are frequently encountered in large eddy simulation. The new scheme is composed of a linear fourth-order central scheme term and an artificial viscosity term. These two terms are connected by a nonlinear weight. The proposed nonlinear weight is designed based on Fourier analysis, rather than Taylor analysis, to guarantee a spectral-like resolution. Moreover, the accuracy is not affected by the optimisation, and the new scheme reaches fourth-order accuracy. The new scheme is tested numerically using the one-dimensional diffusion problem, one-dimensional steady viscous Burger’s shock, two-dimensional vortex decaying, three-dimensional isotropic decaying turbulence and fully developed turbulent channel flow. All the tests confirm that the new scheme has spectral-like resolution and can improve the accuracy of the energy spectrum, dissipation rate and high-order statistics of turbulent flows.

Direct numerical simulation of supersonic turbulent boundary layer subjected to a curved compression ramp

Numerical investigations on a supersonic turbulent boundary layer over a longitudinal curved compression ramp are conducted using direct numerical simulation for a free stream Mach number M∞ = 2.9 and Reynolds number Reθ = 2300. The total turning angle is 24°, and the concave curvature radius is 15 times the thickness of the incoming turbulent boundary layer. Under the selected conditions, the shock foot is transferred to a fan of the compression wave because of the weaker adverse pressure gradient. The time-averaged flow-field in the curved ramp is statistically attached where the instantaneous flow-field is close to the intermittent transitory detachment state. Studies on coherent vortex structures have shown that large-scale vortex packets are enhanced significantly when the concave curvature is aligned in the spanwise direction. Consistent with findings of previous experiments, the effect of the concave curvature on the logarithmic region of the mean velocity profiles is found to be small. The intensi...

A New High Order Accurate Shock Capture Method with Wave Booster

A new kind of shock capturing method is developed. Before applying the high order accurate traditional scheme which is called as base scheme in this paper the fluid parameters are preconditioned in order to control the group velocity. The newly constructed scheme is high order accurate, simple, has high resolution of the shock, and less computer time consumed.