Numerical Simulation of Two-Dimensional Shock/Boundary-Layer Interaction between a Rocket and Booster

A two-dimensional Reynolds-averaged Navier-Stokes solver is applied to analyze the aerodynamic behavior of the Shock/Boundary-Layer interaction of rocket with a boosted The K-ε turbulence model and a finite volume method in a unstructured body-fitted curvilinear coordinates have been used. The results indicate that the separation and the reattachment occur in the Boundary-Layer of the main rocket because of the shock interaction. The shape of the booster nose effects the flow field obviously. In the case of the hemisphere booster nose the pressure has complicate distributions and the separation is very clear. The distance between the booster and main rocket has the evident effect on the flow field. If the distance is smaller the pressure coefficient is bigger the separation zone even the separation bubble occurs.

Interaction between internal solitary waves and the seafloor in the deep sea

Internal solitary wave(ISW),as a typical marine dynamic process in the deep sea,widely exists in oceans and marginal seas worldwide.The interaction between ISW and the seafloor mainly occurs in the bottom boundary layer.For the seabed boundary layer of the deep sea,ISW is the most important dynamic process.This study analyzed the current status,hotspots,and frontiers of research on the interaction between ISW and the seafloor by CiteSpace.Focusing on the action of ISW on the seabed,such as transformation and reaction,a large amount of research work and results were systematically analyzed and summarized.On this basis,this study analyzed the wave–wave interaction and interaction between ISW and the bedform or slope of the seabed,which provided a new perspective for an in‐depth understanding of the interaction between ISW and the seafloor.Finally,the latest research results of the bottom boundary layer and marine engineering stability by ISW were introduced,and the unresolved problems in the current research work were summarized.This study provides a valuable reference for further research on the hazards of ISW to marine engineering geology.

Model for Asymmetry of Shock/Boundary Layer Interactions in Nozzle Flows

The reason for the asymmetry phenomenon of shock/boundary layer interactions(SBLI)in a completely symmetric nozzle with symmetric flow conditions is still an open question.A model for the asymmetry of nozzle flows was proposed based on the properties of fluid entrainment in the mixing layer and momentum conservation.The asymmetry model is deduced based on the nozzle flow with restricted shock separation,and is still applicable for free shock separation.Flow deflection angle at nozzle exit is deduced from this model.Steady numerical simulations are conducted to model the asymmetry of the SBLIs in a planar convergent-divergent nozzle tested by previous researchers.The obtained values of deflection angle based on the numerical results of forced symmetric nozzle flows can judge the asymmetry of flows in a nozzle at some operations.It shows that the entrainment of shear layer on the separation induced by SBLTs is one of the reasons for the asymmetry in the confined SBLIs.

Numerical evaluation of passive control of shock wave/boundary layer interaction on NACA0012 airfoil using jagged wall

Shock formation due to flow compressibility and its interaction with boundary layers has adverse effects on aerodynamic characteristics, such as drag increase and flow separation. The objective of this paper is to appraise the practicability of weakening shock waves and, hence, reducing the wave drag in transonic flight regime using a two-dimensional jagged wall and thereby to gain an appropriate jagged wall shape for future empirical study. Different shapes of the jagged wall, including rectangular, circular, and triangular shapes, were employed. The numerical method was validated by experimental and numerical studies involving transonic flow over the NACA0012 airfoil, and the results presented here closely match previous experimental and numerical results. The impact of parameters, including shape and the length-to-spacing ratio of a jagged wall, was studied on aerodynamic forces and flow field. The results revealed that applying a jagged wall method on the upper surface of an airfoil changes the shock structure significantly and disintegrates it, which in turn leads to a decrease in wave drag. It was also found that the maximum drag coefficient decrease of around 17 % occurs with a triangular shape, while the maximum increase in aerodynamic efficiency（lift-to-drag ratio）of around 10 % happens with a rectangular shape at an angle of attack of 2.26？.

Receptivity of hypersonic boundary layer due to fast-slow acoustics interaction

The objective of receptivity is to investigate the mechanisms by which external disturbances generate unsta- ble waves. In hypersonic boundary layers, a new receptivity process is revealed, which is that fast and slow acoustics through nonlinear interaction can excite the second mode near the lower-branch of the second mode. They can generate a sum-frequency disturbance though nonlinear interaction, which can excite the second mode. This receptivity process is generated by the nonlinear interaction and the nonparal- lel nature of the boundary layer. The receptivity coefficient is sensitive to the wavenumber difference between the sumfrequency disturbance and the lower-branch second mode. When the wavenumber difference is zero, the receptivity coefficient is maximum. The receptivity coefficient decreases with the increase of the wavenumber difference. It is also found that the evolution of the sum-frequency disturbance grows linearly in the beginning. It indicates that the forced term generated by the sum-frequency disturbance resonates with the second mode.

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.

INTERACTION BETWEEN COHERENT STRUCTURES IN WALL REGION OF A TURBULENT BOUNDARY LAYER

Using the idea of general resonant triad of the hydrodynamic stability, the theoretical models for coherent structures in the wall region of a turbulent boundary layer is proposed. The interaction between coherent structures in the wall region of a turbulent boundary layer is studied by combining the compact finite differences of high numerical accuracy and the Fourier spectral hybrid method for solving the three dimensional Navier Stokes equations. In this method, the third order mixed explicit implicit scheme is employed for the time integration. The fifth order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth order center compact schemes for the derivatives in spectral space are descried, respectively. The fourth order compact schemes satisfied by the velocities and pressure in spectral space is derived. As an application, the method is implemented to the wall region of a turbulent boundary to study the interaction between coherent structures. It is found that the numerical results are satisfactory.

STUDY OF SWEPT SHOCK WAVE AND BOUNDARY LAYER INTERACTIONS

This paper presents briefly the recent progress on study of swept shock wave/boundary layer interactions with emphasis on application of zonal analysis and correlation analysis to them. Based on the zonal analysis an overall framework of complicated interaction flow structure including both surface flowfield and space flowfield is discussed. Based on correlation analysis the conical interactions induced by four families of shock wave generators have been discussed in detail. Some control parameter and physical mechanism of conical interaction have been revealed. Finally some aspects of the problem and the prospects for future work are suggested.

Study of interaction between shock wave and unsteady boundary layer

This paper reports theoretical and experimental study of a new type of interaction of a moving shock wave with an unsteady boundary layer. This type of shock wave-boundary layer interaction describes a moving shock wave interaction with an unsteady boundary layer induced by another shock wave and a rarefaction wave. So it is different from the interaction of a stationary shock wave with steady boundary layer, also different from the interaction of a reflected moving shock wave at the end of a shock tube with unsteady boundary layer induced by an incident shock. Geometrical shock dynamics is used for the theoretical analysis of the shock wave-unsteady boundary layer interaction, and a double-driver shock tube with a rarefaction wave bursting diaphragm is used for the experimental investigation in this work.

Entropy Generation through the Interaction of Laminar Boundary-Layer Flows: Sensitivity to Initial Conditions

A modified form of the Townsend equations for the fluctuating velocity wave vectors is applied to the interaction of a longitudinal vortex with a laminar boundary-layer flow. These three-dimensional equations are cast into a Lorenz-format system of equations for the spectral velocity component solutions. Tsallis-form empirical entropic indices are obtained from the solutions of the modified Lorenz equations. These solutions are sensitive to the initial conditions applied to the time-dependent coupled, non-linear differential equations for the spectral velocity components. Eighteen sets of initial conditions for these solutions are examined. The empirical entropic indices yield corresponding intermittency exponents which then yield the entropy generation rates for each set of initial conditions. The flow environment consists of the flow of hydrogen gas with impurities at a given temperature and pressure in the interaction of a longitudinal vortex with a laminar boundary layer flow. Results are presented that indicate a strong correlation of predicted entropy generation rates and the corresponding applied initial conditions. These initial conditions may be ascribed to the turbulence levels in the boundary layer, thus indicating a source for the subsequent entropy generation rates by the interactive instabilities.

THE INTERACTION OF A SHOCK WAVE WITH THE BOUNDARY LAYER IN A REFLECTED SHOCK TUNNEL

The influence of a nontotal reflection on the interaction of a reflected shock wave with the boundary layer in a reflected shock tunnel has been investigated. The calculating method of the velocity, the temperature and the Mach number profiles in the boundary layer in reflected shock fixed coordinates has been obtained. To account for equilibrium real gas effects of nitrogen, the numerical results show that the minimum Mach number in the boundary layer has been moved from the wall into the boundary layer with the increasing of the incident shock Mach number. The minimum Mach number, the shock angle in the bifurcated foot and the jet velocity along the wall to the end plate are reduced owing to the Increasing of the area of nozzle throat. The numerical results are in good agreement with measurements.

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.similar to 4.7) x 10(7)/m. Detailed heat transfer and pressure distributions were measured at fin deflection angles of up to 30 degrees for a sweepback angle of 67.6 degrees. 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 degrees and secondary separation were detected at deflection of theta greater than or equal to 20 degrees. 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 at M greater than or equal to 6 were significantly lower than predicted by the simple prediction techniques at lower Mach numbers.

Structural uncertainty quantification of Reynolds-Averaged Navier–Stokes closures for various shock-wave/boundary layer interaction flows

Accurate prediction of Shock-Wave/Boundary Layer Interaction(SWBLI)flows has been a persistent challenge for linear eddy viscosity models.A major limitation lies in the isotropic representation of the Reynolds stress,as assumed under the Boussinesq approximation.Recent studies have shown promise in improving the prediction capability for incompressible separation flows by perturbing the Reynolds-stress anisotropy tensor.However,it remains uncertain whether this approach is effective for SWBLI flows,which involve compressibility and discontinuity.To address this issue,this study systematically quantifies the structural uncertainty of the anisotropy for oblique SWBLI flows.The eigenspace perturbation method is applied to perturb the anisotropy tensor predicted by the Menter Shear–Stress Transport(SST)model and reveal the impacts of anisotropy on the prediction of quantities of interest,such as separation and reattachment positions,wall static pressure,skin friction,and heat flux.The results demonstrate the potential and reveal the challenges of eigenspace perturbation in improving the SST model.Furthermore,a detailed analysis of turbulent characteristics is performed to identify the source of uncertainty.The findings indicate that eigenspace perturbation primarily affects turbulent shear stress,while the prediction error of the SST model is more related to turbulent kinetic energy.

Effects of Number of Bleed Holes on Shock-Wave/Boundary-Layer Interactions in a Transonic Compressor Stator

An extensive numerical investigation is conducted to characterize the flow separation control in a transonic compressor cascade with a porous bleed.The bleed holes are arranged on the suction surface in a single row,two staggered rows and three staggered rows.For each bleed scheme,five bleed pressure ratios are examined at an inlet Mach number of 1.0.The results indicate that the aerodynamic performance of the cascade is significantly improved by the porous bleed.For the single-row scheme,the maximum reduction in total pressure losses is 57%.For the two-staggered-row and three-staggered-row schemes,there is an optimal bleed pressure ratio of 1.0,and the maximum reductions in total pressure loss are 68% and 75%,respectively.The low loss in the cascade is due to the well-controlled boundary layer.The new local supersonic region created by the bleed hole is the key reason for the improved boundary layer.The vortex induced by side bleeding provides another mechanism for delaying flow separation.Increasing the bleed holes could create multiple local supersonic regions,which reduce the range of the adverse pressure gradient that the boundary layer needs to withstand.This is the reason why cascades with more bleed holes perform better.

Dynamic soil-tunnel interaction in layered half-space for incident plane SH waves

The dynamic soil-tunnel interaction is studied by indirect boundary element method （IBEM）, using the model of a rigid tunnel in layered half-space, which is simplified to a single soil layer on elastic bedrock, subjected to incident plane SH waves. The accuracy of the results is verified through comparison with the analytical solution. It is shown that soil-tunnel interaction in layered half-space is larger than that in homogeneous half-space and this interaction mechanism is essentially different from that of soil-foundation-superstructure interaction.

Dynamic soil-tunnel interaction in layered half-space for incident P-and SV-waves

The dynamic soil-tunnel interaction is studied by the model of a rigid tunnel embedded in layered half-space, which is simplified as a single soil layer on elastic bedrock to the excitation of P- and SV-waves. The indirect boundary element method is used, combined with the Green＇ s function of distributed loads acting on inclined lines. It is shown that the dynamic characteristics of soil-tunnel interaction in layered half-space are different much from that in homoge- neous half-space, and that the mechanism of soil-tunnel interaction is also different much from that of soil-founda- tion-superstructure interaction. For oblique incidence, the tunnel response for in-plane incident SV-waves is com- pletely different from that for incident SH-waves, while the tunnel response for vertically incident SV-wave is very similar to that of vertically incident SH-wave.

Efficient seismic analysis for nonlinear soil-structure interaction with a thick soil layer

The direct finite element method is a type commonly used for nonlinear seismic soil-structure interaction(SSI)analysis.This method introduces a truncated boundary referred to as an artificial boundary meant to divide the soil-structure system into finite and infinite domains.An artificial boundary condition is used on a truncated boundary to achieve seismic input and simulate the wave radiation effect of infinite domain.When the soil layer is particularly thick,especially for a three-dimensional problem,the computational efficiency of seismic SSI analysis is very low due to the large size of the finite element model,which contains an whole thick soil layer.In this paper,an accurate and efficient scheme is developed to solve the nonlinear seismic SSI problem regarding thick soil layers.The process consists of nonlinear site response and SSI analysis.The nonlinear site response analysis is still performed for the whole thick soil layer.The artificial boundary at the bottom of the SSI analysis model is subsequently relocated upward from the bottom of the soil layer(bedrock surface)to the location nearest to the structure as possible.Finally,three types of typical sites and underground structures are adopted with seismic SSI analysis to evaluate the accuracy and efficiency of the proposed efficient analysis scheme.

Recent progress in conical shock wave/boundary layer interaction with spanwise pressure gradient

A common denominator between conical-symmetry and conical shock interaction is the spanwise pressure gradient,which perform more non-uniformity and its interactionflow is more complicated than the spanwise-homogeneous planar shock wave.Recent advances in conical-symmetry and conical shock interactions with turbulent boundary layer are reviewed in specific areas:(i)quasi-conical swept interactions due to compression ramps and sharpfins,(ii)impinging conical shock wave with interactions of plate wall,(iii)laminar double cone in-teractions with consideration of real-gas effects.Substantial success has been achieved in describing the phenomena of the time averaged and instantaneousflow features and the low-frequency unsteadiness,including correlations and coherent structures in the separation bubble,through complementary experimental and numerical studies of swept shock interactions.All available observations are here scrutinized to infer underlying mechanisms of interactions in conicalflow,and provide theoretical foundation and hints forfluidic control techniques.Com-parison with high-fidelity direct numerical simulations is used to quantified the uncertainty of RANS turbulence models in complex interactions.Regarding heat transfer,extensive studies of hypersonicflow over double cone geometries have shown that those can be predicted with reasonable accuracy,even in the presence of high-temperature effects.

Direct numerical simulation of shock wave/boundary layer interaction controlled by steady microjet in a compression ramp

Shock wave/boundary layer interaction in a 24°turning angle of the compression ramp at Mach number 2.9 controlled by steady microjet is investigated using direct numerical simulation.Three different jet spacings which are termed as sparse,moderate and dense are considered,and the induced vortex system and shock structures are compared.A moderate jet spacing configuration is found to generate counter-rotating vortex pairs that transport high-momentum fluid towards the vicinity of wall and strengthen the boundary layer to resist separation,reducing the separation region.The dense jet spacing configuration creates a larger momentum deficit region,reducing the friction downstream of the corner.Analysis of pressure and pressure gradient reveals that dense jet spacing configuration reduces the intensity of separation shock.The impact of varying jet spacings on the turbulent kinetic energy transport mechanism is also investigated by decomposing the budget terms in the transport equation.Furthermore,the spectral characteristics of the separation region are studied using power spectral density and dynamic mode decomposition methods,revealing that moderate jet spacing configuration suppresses low-frequency fluctuations in the separation region.

Selective enhancement of oblique waves caused by finite amplitude second mode in supersonic boundary layer

Nonlinear interactions of the two-dimensional （2D） second mode with oblique modes are studied numerically in a Mach 6.0 fiat-plate boundary layer, focusing on its selective enhancement effect on amplification of different oblique waves. Evolution of oblique modes with various frequencies and spanwise wavenumbers in the presence of 2D second mode is simulated successively, using a modified parabolized stability equation （PSE） method, which is able to simulate interaction of two modes with different frequen- cies efficiently. Numerical results show that oblique modes in a broad band of frequencies and spanwise wavenumbers can be enhanced by the finite amplitude 2D second mode instability wave. The enhancement effect is accomplished by interaction of the 2D second mode, the oblique mode, and a forced mode with difference frequency. Two types of oblique modes are found to be more amplified, i.e., oblique modes with frequency close to that of the 2D second mode and low-frequency first mode oblique waves. Each of them may correspond to one type of transition routes found in transition experiments. The spanwise wavenumber of the oblique wave preferred by the nonlinear interaction is also determined by numerical simulations.