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.展开更多
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 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.展开更多
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 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.展开更多
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.展开更多
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 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.展开更多
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 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.展开更多
Twin boundary(TB)is a special and fundamental internal interface that plays a key role in altering the mechanical and physical properties of materials.However,the atomistic deformation mechanism of TB re-mains under d...Twin boundary(TB)is a special and fundamental internal interface that plays a key role in altering the mechanical and physical properties of materials.However,the atomistic deformation mechanism of TB re-mains under debate,of which the most concerned aspect is how TB would affect the mechanical strength and plasticity of a material.Herein,we introduce our new discovery that the pseudoelastic strain of a TB can recover with decomposition and escape of pile-up dislocations,demonstrated by imposing a sponta-neous pseudoelastic deformation with recoverable plastic bending strain up to 5.1%on a TB.We found that the steps on the curved TB gradually annihilated during the migration of the TB,which was in-duced by the slip of decomposition dislocations on the TB.The TB not only provides local strain harden-ing through interaction with dislocations during the loading stage but also acts as a channel for the fast movement of decomposition dislocations during the recovery stage.Beside,the TB can maintain excellent pseudoelasticity under a multicycle bending test,which may play an important role in improving the fa-tigue resistance of materials.These findings could open up a new avenue for optimizing the mechanical properties of materials by manipulating their twin boundaries at the nanoscale.展开更多
A reined global-local approach based on the scaled boundary inite element method（SBFEM） is proposed to improve the accuracy of predicted singular stress ield. The proposed approach is carried out in conjunction with...A reined global-local approach based on the scaled boundary inite element method（SBFEM） is proposed to improve the accuracy of predicted singular stress ield. The proposed approach is carried out in conjunction with two steps. First, the entire structure is analyzed by employing an arbitrary numerical method. Then, the interested region, which contains stress singularity, is re-solved using the SBFEM by placing the scaling center right at the singular stress point with the boundary conditions evaluated from the irst step imposed along the whole boundary including the side-faces. Beneiting from the semi-analytical nature of the SBFEM, the singular stress ield can be predicted accurately without highly reined meshes. It provides the FEM or other numerical methods with a rather simple and convenient way to improve the accuracy of stress analysis. Numerical examples validate the effectiveness of the proposed approach in dealing with various kinds of problems.展开更多
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,t...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.展开更多
Shock wave/boundary layer interaction(SWBLI)continues to pose a significant chal-lenge in the field of aerospace engineering.This paper aims to address this issue by proposing a novel approach for predicting aerodynam...Shock wave/boundary layer interaction(SWBLI)continues to pose a significant chal-lenge in the field of aerospace engineering.This paper aims to address this issue by proposing a novel approach for predicting aerodynamic coefficients and heat trans-fer in viscous supersonic and hypersonic flows using a high-order flux reconstruction technique.Currently,finite volume methods are extensively employed for the compu-tation of skin aerodynamic coefficients and heat transfer.Nevertheless,these numerical methods exhibit considerable susceptibility to a range of factors,including the inviscid flux function and the computational mesh.The application of high-order flux recon-struction techniques offers promising potential in alleviating these challenges.In contrast to other high-order methods,the flux reconstruction is combined with the lat-tice Boltzmann flux solver in this study.The current method evaluates the common inviscid flux at the cell interface by locally reconstructing the lattice Boltzmann equa-tion solution from macroscopic flow variables at solution points.Consequently,this framework performs a positivity-preserving,entropy-based adaptive filtering method for shock capturing.The present approach is validated by simulating the double Mach reflection,and then simulating some typical viscous problems.The results demonstrate that the current method accurately predicts aerodynamic coefficients and heat trans-fer,providing valuable insights into the application of high-order methods for shock wave/boundary layer interaction.展开更多
Deformation twinning, i.e., twin nucleation and twin growth （or twin boundary migration, TBM） activated by impinged basal slip at a symmetrical tilt grain boundary in HCP Mg, was examined with molecular dynamics （M...Deformation twinning, i.e., twin nucleation and twin growth （or twin boundary migration, TBM） activated by impinged basal slip at a symmetrical tilt grain boundary in HCP Mg, was examined with molecular dynamics （MD） simulations. The results show that the {1^-1^-21}-type twinning acts as the most preferential mode of twinning. Once such twins are formed, they are almost ready to grow. The TBM of such twins is led by pure atomic shuffling events. A secondary mode of twinning can also occur in our simulations. The {112^-2} twinning is observed at 10 K as the secondary twin. This secondary mode of twinning shows different energy barriers for nucleation as well as for growth compared with the {1^-1^-21}-type twining. In particular, TBMs in this case is triggered intrinsically by pyramidal slip at its twin boundary.展开更多
Three Plasma Synthetic Jet Actuators(PSJA)under the high-frequency actuation are used to control the Shock Wave Boundary Layer Interaction(SWBLI),a high-speed schlieren image processing method based on spatial Fourier...Three Plasma Synthetic Jet Actuators(PSJA)under the high-frequency actuation are used to control the Shock Wave Boundary Layer Interaction(SWBLI),a high-speed schlieren image processing method based on spatial Fourier transform as well as snapshot proper orthogonal decomposition were used to study the control effect of high-frequency plasma synthetic jet on lowfrequency unsteadiness of SWBLI.The analysis of the base flow shows that the separated shock wave actually has both large-and small-amplitude vibrations at low frequency.And the results revealed that the PSJA with an operating frequency of 2 k Hz has the ability to reduce the energy of low-frequency component of shock wave motion,indicating that the 2 k Hz actuation can effectively suppress low-frequency unsteadiness of the separated wave.Compared with the actuation frequency of 2 k Hz,the energy of low-frequency component of the shock wave is enhanced under the8 k Hz actuation,which aggravates the low-frequency unsteady motion of the shock wave.It is likely that the actuation frequency is too high,thus the intensity of the precursor shock wave induced by PSJA becomes weaker.Additionally,as the 4 k Hz actuation is applied,the pulsation of the separation region was enhanced,it is speculated that the actuation frequency is coupled with the oscillation frequency of the separation region.展开更多
基金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(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.
基金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.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.
基金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.
基金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 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.
文摘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.
基金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.
文摘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.
基金Y.Y.is supported by the Natural Science Foundation of China(Nos.51922017 and 51972009).
文摘Twin boundary(TB)is a special and fundamental internal interface that plays a key role in altering the mechanical and physical properties of materials.However,the atomistic deformation mechanism of TB re-mains under debate,of which the most concerned aspect is how TB would affect the mechanical strength and plasticity of a material.Herein,we introduce our new discovery that the pseudoelastic strain of a TB can recover with decomposition and escape of pile-up dislocations,demonstrated by imposing a sponta-neous pseudoelastic deformation with recoverable plastic bending strain up to 5.1%on a TB.We found that the steps on the curved TB gradually annihilated during the migration of the TB,which was in-duced by the slip of decomposition dislocations on the TB.The TB not only provides local strain harden-ing through interaction with dislocations during the loading stage but also acts as a channel for the fast movement of decomposition dislocations during the recovery stage.Beside,the TB can maintain excellent pseudoelasticity under a multicycle bending test,which may play an important role in improving the fa-tigue resistance of materials.These findings could open up a new avenue for optimizing the mechanical properties of materials by manipulating their twin boundaries at the nanoscale.
基金supported by the National Key Research and Development plan (Grant No. 2016YFB0201001)the Science Fund for Creative Research Groups of the National Natural Science Foundation of China (Grant No. 51421064)the National Natural Science Foundation of China (Grant No. 51138001)
文摘A reined global-local approach based on the scaled boundary inite element method（SBFEM） is proposed to improve the accuracy of predicted singular stress ield. The proposed approach is carried out in conjunction with two steps. First, the entire structure is analyzed by employing an arbitrary numerical method. Then, the interested region, which contains stress singularity, is re-solved using the SBFEM by placing the scaling center right at the singular stress point with the boundary conditions evaluated from the irst step imposed along the whole boundary including the side-faces. Beneiting from the semi-analytical nature of the SBFEM, the singular stress ield can be predicted accurately without highly reined meshes. It provides the FEM or other numerical methods with a rather simple and convenient way to improve the accuracy of stress analysis. Numerical examples validate the effectiveness of the proposed approach in dealing with various kinds of problems.
基金the financial support provided by the National Science and Technology Major Project (2017-Ⅱ-0007-0021)。
文摘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.
基金This study was supported by the National Natural Science Foundation of China(Grant No.12072158)the Natural Science Foundation of Jiangsu Province(Grant No.BK20231437)the Research Fund of Key Laboratory of Computational Aerodynamics,AVIC Aerodynamics Research Institute(Grant No.YL2022XFX0402).
文摘Shock wave/boundary layer interaction(SWBLI)continues to pose a significant chal-lenge in the field of aerospace engineering.This paper aims to address this issue by proposing a novel approach for predicting aerodynamic coefficients and heat trans-fer in viscous supersonic and hypersonic flows using a high-order flux reconstruction technique.Currently,finite volume methods are extensively employed for the compu-tation of skin aerodynamic coefficients and heat transfer.Nevertheless,these numerical methods exhibit considerable susceptibility to a range of factors,including the inviscid flux function and the computational mesh.The application of high-order flux recon-struction techniques offers promising potential in alleviating these challenges.In contrast to other high-order methods,the flux reconstruction is combined with the lat-tice Boltzmann flux solver in this study.The current method evaluates the common inviscid flux at the cell interface by locally reconstructing the lattice Boltzmann equa-tion solution from macroscopic flow variables at solution points.Consequently,this framework performs a positivity-preserving,entropy-based adaptive filtering method for shock capturing.The present approach is validated by simulating the double Mach reflection,and then simulating some typical viscous problems.The results demonstrate that the current method accurately predicts aerodynamic coefficients and heat trans-fer,providing valuable insights into the application of high-order methods for shock wave/boundary layer interaction.
基金Project(2012CB932202)supported by the National Basic Research Program of ChinaProjects(50890174,50971088)supported by the National Natural Science Foundation of China
文摘Deformation twinning, i.e., twin nucleation and twin growth （or twin boundary migration, TBM） activated by impinged basal slip at a symmetrical tilt grain boundary in HCP Mg, was examined with molecular dynamics （MD） simulations. The results show that the {1^-1^-21}-type twinning acts as the most preferential mode of twinning. Once such twins are formed, they are almost ready to grow. The TBM of such twins is led by pure atomic shuffling events. A secondary mode of twinning can also occur in our simulations. The {112^-2} twinning is observed at 10 K as the secondary twin. This secondary mode of twinning shows different energy barriers for nucleation as well as for growth compared with the {1^-1^-21}-type twining. In particular, TBMs in this case is triggered intrinsically by pyramidal slip at its twin boundary.
基金supported by the National Key R&D Program of China(No.2019YFA0405300)the National Natural Science Foundation of China(Nos.51907205 and 12002363)the Foundation Strengthening Fund,China(Nos.2019077 and 2019-053)。
文摘Three Plasma Synthetic Jet Actuators(PSJA)under the high-frequency actuation are used to control the Shock Wave Boundary Layer Interaction(SWBLI),a high-speed schlieren image processing method based on spatial Fourier transform as well as snapshot proper orthogonal decomposition were used to study the control effect of high-frequency plasma synthetic jet on lowfrequency unsteadiness of SWBLI.The analysis of the base flow shows that the separated shock wave actually has both large-and small-amplitude vibrations at low frequency.And the results revealed that the PSJA with an operating frequency of 2 k Hz has the ability to reduce the energy of low-frequency component of shock wave motion,indicating that the 2 k Hz actuation can effectively suppress low-frequency unsteadiness of the separated wave.Compared with the actuation frequency of 2 k Hz,the energy of low-frequency component of the shock wave is enhanced under the8 k Hz actuation,which aggravates the low-frequency unsteady motion of the shock wave.It is likely that the actuation frequency is too high,thus the intensity of the precursor shock wave induced by PSJA becomes weaker.Additionally,as the 4 k Hz actuation is applied,the pulsation of the separation region was enhanced,it is speculated that the actuation frequency is coupled with the oscillation frequency of the separation region.
