As a multidisciplinary phenomenon,panel aeroelasticity in shock-dominated flow is featured by two primary interactions:Fluid-Structure Interactions(FSIs)and Shock-Boundary Layer Interactions(SBLIs).The former raises s...As a multidisciplinary phenomenon,panel aeroelasticity in shock-dominated flow is featured by two primary interactions:Fluid-Structure Interactions(FSIs)and Shock-Boundary Layer Interactions(SBLIs).The former raises structural concerns,and the latter is of aerodynamic interest.Thus,panel aeroelasticity in shock-dominated flow represents a vital topic for the development and optimization of supersonic vehicles and propulsion systems.This review systematically summarizes recent advances in the methodologies applied to capture structural and fluid dynamics,including theoretical models,numerical simulations,and wind tunnel experiments.The application of data-driven modal decomposition,an advanced technique to extract physically crucial features,on the topic is introduced.From the perspective of FSIs,the distinctive aeroelastic behaviors in shock-dominated flow,including hysteresis phenomena and nonlinear responses,are highlighted.From the perspective of SBLIs,the modifications in their spatial and temporal characteristics imposed by the aeroelastic responses are emphasized.Motivated by the interaction between the shock waves and structural response,different strategies have been proposed to implement aeroelastic suppression and shock control,which have the potential to enhance structural safety and aerodynamic performance in the next generation of high-speed flight vehicles.展开更多
The stability of supersonic inlets faces challenges due to various changes in flight conditions,and flow control methods that address shock wave/boundary layer interactions under only one set of conditions cannot meet...The stability of supersonic inlets faces challenges due to various changes in flight conditions,and flow control methods that address shock wave/boundary layer interactions under only one set of conditions cannot meet developmental requirements.This paper proposes an adaptive bump control scheme and employs dynamic mesh technology for numerical simulation to investigate the unsteady control effects of adaptive bumps.The obtained results indicate that the use of moving bumps to control shock wave/boundary layer interactions is feasible.The adaptive control effects of five different bump speeds are evaluated.Within the range of bump speeds studied,the analysis of the flow field structure reveals the patterns of change in the separation zone area during the control process,as well as the relationship between the bump motion speed and the control effect on the separation zone.It is concluded that the moving bump endows the boundary layer with additional energy.展开更多
Three-dimensional curved shock wave/boundary layer interaction with streamwise and spanwise curvatures widely exists in practical aerodynamic design.To explore the effects of composite shock curvatures on boundary lay...Three-dimensional curved shock wave/boundary layer interaction with streamwise and spanwise curvatures widely exists in practical aerodynamic design.To explore the effects of composite shock curvatures on boundary layer separation,a canonical model with a cone placed above plate was utilized as a reference.Configurations of straight,convex,and concave conical shock waves inducing the curved conical shock wave/boundary layer interactions were studied,using CFD based on Reynolds-averaged numerical simulation method.The flow structure and separation region of each case were discussed quantitively on the symmetry plane,flat plate,and plane perpendicular to flow direction,respectively.The focus of the analysis was on the characteristic patterns of separation scale variation in the streamwise and spanwise directions,which were observed to consistently change with respect to both directions with alterations in the incident shock wave shape.A simplified control volume model was established to qualitatively discuss the influence source of curved shock waves on separation scales,based on mass conservation equations.The results suggest that the curved shock wave has a holistic effect on separation,which is not solely dependent on the shock foot strength.展开更多
A Discrete Boltzmann Method(DBM)with a Maxwell-type boundary condition is constructed to investigate the influence of rarefaction on laminar Shock Wave/Boundary Layer Interaction(SWBLI).Due to the complexity of compre...A Discrete Boltzmann Method(DBM)with a Maxwell-type boundary condition is constructed to investigate the influence of rarefaction on laminar Shock Wave/Boundary Layer Interaction(SWBLI).Due to the complexity of compressible flow,a Knudsen number vector Kn,whose components include the local Knudsen numbers such as Kn_(ρ)and Kn_(U),is introduced to characterize the local structures,where Kn_(ρ)and Kn_(U)are Knudsen numbers defined in terms of the density and velocity interfaces,respectively.Since first focusing on the steady state of SWBLI,the DBM considers up to the second-order Kn_(ρ)(rarefaction/non-equilibrium)effects.The model is validated using Mach number 2 SWBLI and the necessity of using DBM with sufficient physical accuracy is confirmed by the shock collision problem.Key findings include the following:the leading-edge shock wave increases the local density Knudsen number Kn_(ρ)and eventually leads to the failure of linear constitutive relations in the Navier-Stokes(N-S)model and surely also in the lower-order DBM;the non-equilibrium effect differences in regions behind the leading-edge shock wave are primarily correlated with Kn_(ρ),while in the separation region are primarily correlated with Kn_(U);the non-equilibrium quantities D_(2)and D_(4,2),as well as the viscous entropy production rate S_(NOMF)can be used to identify the separation zone.The findings clarify various effects and main mechanisms in different regions associated with SWBLI,which are concealed in N-S model.展开更多
Cowl-induced incident Shock Wave/Boundary Layer Interactions (SWBLI) under the influence of gradual expansion waves are frequently observed in supersonic inlets. However, the analysis and prediction of interaction len...Cowl-induced incident Shock Wave/Boundary Layer Interactions (SWBLI) under the influence of gradual expansion waves are frequently observed in supersonic inlets. However, the analysis and prediction of interaction lengths have not been sufficiently investigated. First, this study presents a theoretical scaling analysis and validates it through wind tunnel experiments. It conducts detailed control volume analysis of mass conservation, considering the differences between inviscid and viscous cases. Then, three models for analysing interaction length under gradual expansion waves are derived. Related experiments using schlieren photography are conducted to validate the models in a Mach 2.73 flow. The interaction scales are captured at various relative distances between the shock impingement location and the expansion regions with wedge angles ranging from 12° to 15° and expansion angles of 9°, 12°, and 15°. Three trend lines are plotted based on different expansion angles to depict the relationship between normalised interaction length and normalised interaction strength metric. In addition, the relationship between the coefficients of the trend line and the expansion angles is introduced to predict the interaction length influenced by gradual expansion waves. Finally, the estimation of normalised interaction length is derived for various coefficients within a unified form.展开更多
It is of great significance to improve the accuracy of turbulence models in shock-wave/ boundary layer interaction flow. The relationship between the pressure gradient, as well as the shear layer, and the development ...It is of great significance to improve the accuracy of turbulence models in shock-wave/ boundary layer interaction flow. The relationship between the pressure gradient, as well as the shear layer, and the development of turbulent kinetic energy in impinging shock-wave/turbulent bound- ary layer interaction flow at Mach 2.25 is analyzed based on the data of direct numerical simulation （DNS）. It is found that the turbulent kinetic energy is amplified by strong shear in the separation zone and the adverse pressure gradient near the separation point. The pressure gradient was non-dimensionalised with local density, velocity, and viscosity. Spalart Allmaras （S A） model is modified by introducing the non-dimensional pressure gradient into the production term of the eddy viscosity transportation equation. Simulation results show that the production and dissipation of eddy viscosity are strongly enhanced by the modification of S-A model. Compared with DNS and experimental data, the wall pressure and the wall skin friction coefficient as well as the velocity profile of the modified S-A model are obviously improved. Thus it can be concluded that the mod- ification of S-A model with the pressure gradient can improve the predictive accuracy for simulat- ing the shock-wave/turbulent boundary laver interaction.展开更多
The effect of magnetohydrodynamic(MHD)plasma actuators on the control of hypersonic shock wave/turbulent boundary layer interactions is investigated here using Reynolds-averaged Navier-Stokes calculations with low mag...The effect of magnetohydrodynamic(MHD)plasma actuators on the control of hypersonic shock wave/turbulent boundary layer interactions is investigated here using Reynolds-averaged Navier-Stokes calculations with low magnetic Reynolds number approximation.A Mach 5 oblique shock/turbulent boundary layer interaction was adopted as the basic configuration in this numerical study in order to assess the effects of flow control using different combinations of magnetic field and plasma.Results show that just the thermal effect of plasma under experimental actuator parameters has no significant impact on the flow field and can therefore be neglected.On the basis of the relative position of control area and separation point,MHD control can be divided into four types and so effects and mechanisms might be different.Amongst these,D-type control leads to the largest reduction in separation length using magnetically-accelerated plasma inside an isobaric dead-air region.A novel parameter for predicting the shock wave/turbulent boundary layer interaction control based on Lorentz force acceleration is then proposed and the controllability of MHD plasma actuators under different MHD interaction parameters is studied.The results of this study will be insightful for the further design of MHD control in hypersonic vehicle inlets.展开更多
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...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.展开更多
The interaction length induced by Shock Wave/Turbulent Boundary-Layer Interactions(SWTBLIs)in the hypersonic flow was investigated using a scaling analysis,in which the interaction length normalized by the displacemen...The interaction length induced by Shock Wave/Turbulent Boundary-Layer Interactions(SWTBLIs)in the hypersonic flow was investigated using a scaling analysis,in which the interaction length normalized by the displacement thickness of boundary layer was correlated with a corrected non-dimensional separation criterion across the interaction after accounting for the wall temperature effects.A large number of hypersonic SWTBLIs were compiled to examine the scaling analysis over a wide range of Mach numbers,Reynolds numbers,and wall temperatures.The results indicate that the hypersonic SWTBLIs with low Reynolds numbers collapse on the supersonic SWTBLIs,while the hypersonic cases with high Reynolds numbers show a more rapid growth of the interaction length than that with low Reynolds numbers.Thus,two scaling relationships are identified according to different Reynolds numbers for the hypersonic SWTBLIs.The scaling analysis provides valuable guidelines for engineering prediction of the interaction length,and thus,enriches the knowledge of hypersonic SWTBLIs.展开更多
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 app...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？.展开更多
Flow separation due to shock wave/boundary layer interaction is dominated in blade passage with supersonic relative incoming flow,which always accompanies aerodynamic performance penalties.A loss reduction method for ...Flow separation due to shock wave/boundary layer interaction is dominated in blade passage with supersonic relative incoming flow,which always accompanies aerodynamic performance penalties.A loss reduction method for smearing the passage shock foot via Shock Control Bump(SCB)located on transonic compressor rotor blade suction side is implemented to shrink the region of boundary layer separation.The curved windward section of SCB with constant adverse pressure gradient is constructed ahead of passage shock-impingement point at design rotor speed of Rotor 37 to get the improved model.Numerical investigations on both two models have been conducted employing Reynolds-Averaged Navier-Stokes(RANS)method to reveal flow physics of SCB.Comparisons and analyses on simulation results have also been carried out,showing that passage shock foot of baseline is replaced with a family of compression waves and a weaker shock foot for moderate adverse pressure gradient as well as suppression of boundary layer separations and secondary flow of low-momentum fluid within boundary layer.It is found that adiabatic efficiency and total pressure ratio of improved blade exceeds those of baseline at 95%-100%design rotor speed,and then slightly worsens with decrease of rotatory speed till both equal below 60%rated speed.The investigated conclusion implies a potential promise for future practical applications of SCB in both transonic and supersonic compressors.展开更多
In shock tube experiments,the interaction between the reflected shock and boundary layer can induce shock bifurcation and weak ignition.The weak ignition can greatly affect the ignition delay time measurement in a sho...In shock tube experiments,the interaction between the reflected shock and boundary layer can induce shock bifurcation and weak ignition.The weak ignition can greatly affect the ignition delay time measurement in a shock tube experiment.In this work,two-dimensional simulations considering detailed chemistry and transport are conducted to investigate the shock bifurcation and non-uniform ignition behind a retlected shock.The objectives are to interpret the formation of shock bifurcation induced by the reflected shock and boundary layer interaction and to investigate the weak ignition and its transition to strong ignition for both hydrogen and dimethyl ether.It is found that the non-uniform reflection of the incident shock at the end wall produces a wedge-shaped oblique shock foot at the wall.The wedge-shaped structure results in strong interactions between reflected shock and boundary layer,which induces the shock bifurcation.It is demonstrated that the local high-temperature spots at the foot of the bifurcated shock is caused by viscous dissipation and pressure work.As the post-reflected shock temperature increases,the transition from weak ignition to strong ignition in a stoichiometric hydrogen/oxygen mixture is observed.The relative sensitivity of ignition delay time to the post-rellected shock temperature is introduced to characterize the appearance of weak ignition behind the reflected shock.Unlike in the hydrogen/oxygen mixture,weak ignition is not observed in the stoichiometric dimethyl-ether/oxygen mixture since it has a relatively longer ignition delay time and smaller relative sensitivity.展开更多
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 achieveme...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.展开更多
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...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.展开更多
This study investigates the three-dimensional(3D)effects introduced by the end walls for an aspect ratio of1 in ramp-induced shock wave boundary layer interactions.The simulations are performed using a symmetry bounda...This study investigates the three-dimensional(3D)effects introduced by the end walls for an aspect ratio of1 in ramp-induced shock wave boundary layer interactions.The simulations are performed using a symmetry boundary condition in the spanwise direction at free-stream Mach numbers in 3D.The simulations are performed using an in-house compressible supersonic solver“Open SBLIFVM”.Two free stream Mach numbers 2.5,and3 are used in the current work,and the simulated results are compared with the aspect ratio 1 simulations by Mangalagiri and Jammy.The inflow is initialized with a similarity solution;its Reynolds number based on the boundary layer thickness is adjusted such that the Reynolds number at the start of the ramp is kept at 3×10^(5)for all simulations.From the results,it is evident that the introduction of sidewalls resulted in a shorter centerline separation length when compared with the two-dimensional(2D)simulations.This contradicts the results at Mach 2 by Mangalgiri and Jammy where the vortex observed at Mach 2 in the central separation region disappeared with increasing free-stream Mach number.Additionally,the topology of interaction shifted from owl-like separation of the second kind to the first kind when the freestream Mach number increased from2 to 2.5.It can be concluded that the interaction topology is crucial to the increase or decrease of the central separation length when compared to 2D simulations.展开更多
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,...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.展开更多
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 ...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.展开更多
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 ...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.展开更多
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 an...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.展开更多
Numerical simulations and theoretical models are developed in this paper for the Detonation-Wave/Boundary-Layer Interactions(DWBLIs)under reflections.Transient flow fields demonstrate the highly non-stationarity of th...Numerical simulations and theoretical models are developed in this paper for the Detonation-Wave/Boundary-Layer Interactions(DWBLIs)under reflections.Transient flow fields demonstrate the highly non-stationarity of the DWBLIs when Mach Reflection(MR)occur,and subsequent analyses show that the subsonic region introduced by the boundary layer exacerbates the instability.Further quantitative analyses show that viscosity has little effect on propulsive performance and the separation wave can be considered as an oblique detonation wave.Influence parameters to DWBLIs such as combustion chamber height,incoming Mach number,equivalence ratio,and inlet channel length are categorized and studied.Besides simulations,theoretical analytical modeling is established for Regular Reflection(RR)and MR of DWBLIs.Multiple formulas for the separation zone length are obtained according to the mass conservation under different transformation type between inviscid and viscid reflections.Comparison with the numerical simulations verifies the validity of the model and it can be further generalized to the curved DWBLIs.The developed model makes the theoretical solution process of DWBLIs possible and provides the key foundation for further analysis and solution.展开更多
基金supported by the National Natural Science Foundation of China(No.12372233)the Innovation Foundation for Doctor Dissertation of Northwestern Polytechnical University,China(No.25GH01020005)the“111 Project”of China(No.B17037)。
文摘As a multidisciplinary phenomenon,panel aeroelasticity in shock-dominated flow is featured by two primary interactions:Fluid-Structure Interactions(FSIs)and Shock-Boundary Layer Interactions(SBLIs).The former raises structural concerns,and the latter is of aerodynamic interest.Thus,panel aeroelasticity in shock-dominated flow represents a vital topic for the development and optimization of supersonic vehicles and propulsion systems.This review systematically summarizes recent advances in the methodologies applied to capture structural and fluid dynamics,including theoretical models,numerical simulations,and wind tunnel experiments.The application of data-driven modal decomposition,an advanced technique to extract physically crucial features,on the topic is introduced.From the perspective of FSIs,the distinctive aeroelastic behaviors in shock-dominated flow,including hysteresis phenomena and nonlinear responses,are highlighted.From the perspective of SBLIs,the modifications in their spatial and temporal characteristics imposed by the aeroelastic responses are emphasized.Motivated by the interaction between the shock waves and structural response,different strategies have been proposed to implement aeroelastic suppression and shock control,which have the potential to enhance structural safety and aerodynamic performance in the next generation of high-speed flight vehicles.
基金supported by the National Key R&D Program of China(Grant No.2019YFA0405300)the National Natural Science Foundation of China(Grant No.11972368)the Natural Science Foundation of Hunan Province(Grant No.2021JJ10045).
文摘The stability of supersonic inlets faces challenges due to various changes in flight conditions,and flow control methods that address shock wave/boundary layer interactions under only one set of conditions cannot meet developmental requirements.This paper proposes an adaptive bump control scheme and employs dynamic mesh technology for numerical simulation to investigate the unsteady control effects of adaptive bumps.The obtained results indicate that the use of moving bumps to control shock wave/boundary layer interactions is feasible.The adaptive control effects of five different bump speeds are evaluated.Within the range of bump speeds studied,the analysis of the flow field structure reveals the patterns of change in the separation zone area during the control process,as well as the relationship between the bump motion speed and the control effect on the separation zone.It is concluded that the moving bump endows the boundary layer with additional energy.
基金the support of the National Natural Science Foundation of China(Nos.12372295,U21B6003,U20A2069,12302389 and 123B2037)。
文摘Three-dimensional curved shock wave/boundary layer interaction with streamwise and spanwise curvatures widely exists in practical aerodynamic design.To explore the effects of composite shock curvatures on boundary layer separation,a canonical model with a cone placed above plate was utilized as a reference.Configurations of straight,convex,and concave conical shock waves inducing the curved conical shock wave/boundary layer interactions were studied,using CFD based on Reynolds-averaged numerical simulation method.The flow structure and separation region of each case were discussed quantitively on the symmetry plane,flat plate,and plane perpendicular to flow direction,respectively.The focus of the analysis was on the characteristic patterns of separation scale variation in the streamwise and spanwise directions,which were observed to consistently change with respect to both directions with alterations in the incident shock wave shape.A simplified control volume model was established to qualitatively discuss the influence source of curved shock waves on separation scales,based on mass conservation equations.The results suggest that the curved shock wave has a holistic effect on separation,which is not solely dependent on the shock foot strength.
基金support from the National Key R&D Program of China(No.2020YFC2201100)the Foundation of National Key Laboratory of Shock Wave and Detonation Physics,China(No.JCKYS2023212003)+1 种基金the National Natural Science Foundation of China(No.12172061)the Opening Project of State Key Laboratory of Explosion Science and Safety Protection(Beijing Institute of Technology)(No.KFJJ25-02M).
文摘A Discrete Boltzmann Method(DBM)with a Maxwell-type boundary condition is constructed to investigate the influence of rarefaction on laminar Shock Wave/Boundary Layer Interaction(SWBLI).Due to the complexity of compressible flow,a Knudsen number vector Kn,whose components include the local Knudsen numbers such as Kn_(ρ)and Kn_(U),is introduced to characterize the local structures,where Kn_(ρ)and Kn_(U)are Knudsen numbers defined in terms of the density and velocity interfaces,respectively.Since first focusing on the steady state of SWBLI,the DBM considers up to the second-order Kn_(ρ)(rarefaction/non-equilibrium)effects.The model is validated using Mach number 2 SWBLI and the necessity of using DBM with sufficient physical accuracy is confirmed by the shock collision problem.Key findings include the following:the leading-edge shock wave increases the local density Knudsen number Kn_(ρ)and eventually leads to the failure of linear constitutive relations in the Navier-Stokes(N-S)model and surely also in the lower-order DBM;the non-equilibrium effect differences in regions behind the leading-edge shock wave are primarily correlated with Kn_(ρ),while in the separation region are primarily correlated with Kn_(U);the non-equilibrium quantities D_(2)and D_(4,2),as well as the viscous entropy production rate S_(NOMF)can be used to identify the separation zone.The findings clarify various effects and main mechanisms in different regions associated with SWBLI,which are concealed in N-S model.
基金co-supported by the National Natural Science Foundation of China (No. 12172175)the National Science and Technology Major Project, China (No. J2019-II0014-0035)the Science Center for Gas Turbine Project, China (Nos. P2022-C-II-002-001, P2022-A-II-002-001)
文摘Cowl-induced incident Shock Wave/Boundary Layer Interactions (SWBLI) under the influence of gradual expansion waves are frequently observed in supersonic inlets. However, the analysis and prediction of interaction lengths have not been sufficiently investigated. First, this study presents a theoretical scaling analysis and validates it through wind tunnel experiments. It conducts detailed control volume analysis of mass conservation, considering the differences between inviscid and viscous cases. Then, three models for analysing interaction length under gradual expansion waves are derived. Related experiments using schlieren photography are conducted to validate the models in a Mach 2.73 flow. The interaction scales are captured at various relative distances between the shock impingement location and the expansion regions with wedge angles ranging from 12° to 15° and expansion angles of 9°, 12°, and 15°. Three trend lines are plotted based on different expansion angles to depict the relationship between normalised interaction length and normalised interaction strength metric. In addition, the relationship between the coefficients of the trend line and the expansion angles is introduced to predict the interaction length influenced by gradual expansion waves. Finally, the estimation of normalised interaction length is derived for various coefficients within a unified form.
基金supported by the National Natural Science Foundation of China (No.11302012,51376001,51136003)the National Basic Research Program of China (No.2012CB720205)+3 种基金the National Magnetic Confinement Fusion Research Program of China (No.2012GB102006)the Aeronautical Science Foundation of China (No.2012ZB51014)the ‘‘111’’ Project(No.B08009)the Astronautical Technology Innovation Foundation of China
文摘It is of great significance to improve the accuracy of turbulence models in shock-wave/ boundary layer interaction flow. The relationship between the pressure gradient, as well as the shear layer, and the development of turbulent kinetic energy in impinging shock-wave/turbulent bound- ary layer interaction flow at Mach 2.25 is analyzed based on the data of direct numerical simulation （DNS）. It is found that the turbulent kinetic energy is amplified by strong shear in the separation zone and the adverse pressure gradient near the separation point. The pressure gradient was non-dimensionalised with local density, velocity, and viscosity. Spalart Allmaras （S A） model is modified by introducing the non-dimensional pressure gradient into the production term of the eddy viscosity transportation equation. Simulation results show that the production and dissipation of eddy viscosity are strongly enhanced by the modification of S-A model. Compared with DNS and experimental data, the wall pressure and the wall skin friction coefficient as well as the velocity profile of the modified S-A model are obviously improved. Thus it can be concluded that the mod- ification of S-A model with the pressure gradient can improve the predictive accuracy for simulat- ing the shock-wave/turbulent boundary laver interaction.
基金Project supported by the National Key R&D Program of China(Nos.2019YFA0405300 and 2019YFA0405203)the Chinese Scholarship Council(CSC)(No.201903170195)。
文摘The effect of magnetohydrodynamic(MHD)plasma actuators on the control of hypersonic shock wave/turbulent boundary layer interactions is investigated here using Reynolds-averaged Navier-Stokes calculations with low magnetic Reynolds number approximation.A Mach 5 oblique shock/turbulent boundary layer interaction was adopted as the basic configuration in this numerical study in order to assess the effects of flow control using different combinations of magnetic field and plasma.Results show that just the thermal effect of plasma under experimental actuator parameters has no significant impact on the flow field and can therefore be neglected.On the basis of the relative position of control area and separation point,MHD control can be divided into four types and so effects and mechanisms might be different.Amongst these,D-type control leads to the largest reduction in separation length using magnetically-accelerated plasma inside an isobaric dead-air region.A novel parameter for predicting the shock wave/turbulent boundary layer interaction control based on Lorentz force acceleration is then proposed and the controllability of MHD plasma actuators under different MHD interaction parameters is studied.The results of this study will be insightful for the further design of MHD control in hypersonic vehicle inlets.
基金The project supported by China Academy of Launch Vehicle Technology
文摘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.
基金supported by the National Natural Science Foundation of China(Nos.11772325 and 11621202)。
文摘The interaction length induced by Shock Wave/Turbulent Boundary-Layer Interactions(SWTBLIs)in the hypersonic flow was investigated using a scaling analysis,in which the interaction length normalized by the displacement thickness of boundary layer was correlated with a corrected non-dimensional separation criterion across the interaction after accounting for the wall temperature effects.A large number of hypersonic SWTBLIs were compiled to examine the scaling analysis over a wide range of Mach numbers,Reynolds numbers,and wall temperatures.The results indicate that the hypersonic SWTBLIs with low Reynolds numbers collapse on the supersonic SWTBLIs,while the hypersonic cases with high Reynolds numbers show a more rapid growth of the interaction length than that with low Reynolds numbers.Thus,two scaling relationships are identified according to different Reynolds numbers for the hypersonic SWTBLIs.The scaling analysis provides valuable guidelines for engineering prediction of the interaction length,and thus,enriches the knowledge of hypersonic SWTBLIs.
文摘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？.
基金the funding from the National Key Research and Development Program of China(No.2016YFB0901402)the Key Project of National Natural Science Foundation of China(No.51790513)。
文摘Flow separation due to shock wave/boundary layer interaction is dominated in blade passage with supersonic relative incoming flow,which always accompanies aerodynamic performance penalties.A loss reduction method for smearing the passage shock foot via Shock Control Bump(SCB)located on transonic compressor rotor blade suction side is implemented to shrink the region of boundary layer separation.The curved windward section of SCB with constant adverse pressure gradient is constructed ahead of passage shock-impingement point at design rotor speed of Rotor 37 to get the improved model.Numerical investigations on both two models have been conducted employing Reynolds-Averaged Navier-Stokes(RANS)method to reveal flow physics of SCB.Comparisons and analyses on simulation results have also been carried out,showing that passage shock foot of baseline is replaced with a family of compression waves and a weaker shock foot for moderate adverse pressure gradient as well as suppression of boundary layer separations and secondary flow of low-momentum fluid within boundary layer.It is found that adiabatic efficiency and total pressure ratio of improved blade exceeds those of baseline at 95%-100%design rotor speed,and then slightly worsens with decrease of rotatory speed till both equal below 60%rated speed.The investigated conclusion implies a potential promise for future practical applications of SCB in both transonic and supersonic compressors.
基金supported by the National Natural Science Foundation of China(Grant Nos.52006001,and 52176096).
文摘In shock tube experiments,the interaction between the reflected shock and boundary layer can induce shock bifurcation and weak ignition.The weak ignition can greatly affect the ignition delay time measurement in a shock tube experiment.In this work,two-dimensional simulations considering detailed chemistry and transport are conducted to investigate the shock bifurcation and non-uniform ignition behind a retlected shock.The objectives are to interpret the formation of shock bifurcation induced by the reflected shock and boundary layer interaction and to investigate the weak ignition and its transition to strong ignition for both hydrogen and dimethyl ether.It is found that the non-uniform reflection of the incident shock at the end wall produces a wedge-shaped oblique shock foot at the wall.The wedge-shaped structure results in strong interactions between reflected shock and boundary layer,which induces the shock bifurcation.It is demonstrated that the local high-temperature spots at the foot of the bifurcated shock is caused by viscous dissipation and pressure work.As the post-reflected shock temperature increases,the transition from weak ignition to strong ignition in a stoichiometric hydrogen/oxygen mixture is observed.The relative sensitivity of ignition delay time to the post-rellected shock temperature is introduced to characterize the appearance of weak ignition behind the reflected shock.Unlike in the hydrogen/oxygen mixture,weak ignition is not observed in the stoichiometric dimethyl-ether/oxygen mixture since it has a relatively longer ignition delay time and smaller relative sensitivity.
基金supported by the National Key Basic Research and Development Program (No.2014CB744100)
文摘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.
基金supported by the National Natural Science Foundation of China (Grants 11332007 and 11202147)the Specialized Research Fund for the Doctoral Program of Higher Education (Grants 20120032120007)
文摘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.
基金sponsored by the Department of Science and Technology,Science and Engineering Research Board(SERB),Core Research(Grant No.CRG/2020/03859)。
文摘This study investigates the three-dimensional(3D)effects introduced by the end walls for an aspect ratio of1 in ramp-induced shock wave boundary layer interactions.The simulations are performed using a symmetry boundary condition in the spanwise direction at free-stream Mach numbers in 3D.The simulations are performed using an in-house compressible supersonic solver“Open SBLIFVM”.Two free stream Mach numbers 2.5,and3 are used in the current work,and the simulated results are compared with the aspect ratio 1 simulations by Mangalagiri and Jammy.The inflow is initialized with a similarity solution;its Reynolds number based on the boundary layer thickness is adjusted such that the Reynolds number at the start of the ramp is kept at 3×10^(5)for all simulations.From the results,it is evident that the introduction of sidewalls resulted in a shorter centerline separation length when compared with the two-dimensional(2D)simulations.This contradicts the results at Mach 2 by Mangalgiri and Jammy where the vortex observed at Mach 2 in the central separation region disappeared with increasing free-stream Mach number.Additionally,the topology of interaction shifted from owl-like separation of the second kind to the first kind when the freestream Mach number increased from2 to 2.5.It can be concluded that the interaction topology is crucial to the increase or decrease of the central separation length when compared to 2D simulations.
基金supported by the National Natural Science Foundation of China(Nos.92252201 and 11721202)。
文摘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.
文摘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.
基金Supported by the National Natural Science Foundation of China(Nos.11972064,92052104)the Key Laboratory of Computational Aerodynamics,AVIC Aerodynamics Research Institute,China(No.YL2022XFX0405)the Fundamental Research Funds for the Central Universities,China.
文摘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.
文摘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.
基金support of the National Natural Science Foundation of China(Nos.U20A2069,U21B6003,12302389 and 12472337)the Advanced Aero-Power Innovation Workstation,China(No.HKCX2024-01-017)。
文摘Numerical simulations and theoretical models are developed in this paper for the Detonation-Wave/Boundary-Layer Interactions(DWBLIs)under reflections.Transient flow fields demonstrate the highly non-stationarity of the DWBLIs when Mach Reflection(MR)occur,and subsequent analyses show that the subsonic region introduced by the boundary layer exacerbates the instability.Further quantitative analyses show that viscosity has little effect on propulsive performance and the separation wave can be considered as an oblique detonation wave.Influence parameters to DWBLIs such as combustion chamber height,incoming Mach number,equivalence ratio,and inlet channel length are categorized and studied.Besides simulations,theoretical analytical modeling is established for Regular Reflection(RR)and MR of DWBLIs.Multiple formulas for the separation zone length are obtained according to the mass conservation under different transformation type between inviscid and viscid reflections.Comparison with the numerical simulations verifies the validity of the model and it can be further generalized to the curved DWBLIs.The developed model makes the theoretical solution process of DWBLIs possible and provides the key foundation for further analysis and solution.
