An experimental study of cavity oscillating flow carried out on subsonic wall jet facilities in an anechoic room is summarized. The jet exit Mach number range is from 0.2 to 0.8. The effects of the flow Mach number (...An experimental study of cavity oscillating flow carried out on subsonic wall jet facilities in an anechoic room is summarized. The jet exit Mach number range is from 0.2 to 0.8. The effects of the flow Mach number ( Ma ) and the cavity depth ( D ) on the oscillation are studied. It is found that for L/D =4, (shallow cavity), the oscillation is mainly due to the self exciting of the free shear layer above the cavity opening, for L/D =2, the acoustic resonance is responsible for the oscillation. Preliminary tests are performed to study the suppression effect of the leading edge tone excitation on cavity flow, and considerable reduction of oscillation has been achieved when Ma ≤0.6.展开更多
Combining the detached eddy simulation(DES)method and Ffowcs Williams-Hawkings(FW-H)equation,the effect of bogie cavity end wall inclination on the flow field and aerodynamic noise in the bogie region is numerically s...Combining the detached eddy simulation(DES)method and Ffowcs Williams-Hawkings(FW-H)equation,the effect of bogie cavity end wall inclination on the flow field and aerodynamic noise in the bogie region is numerically studied.First,the simulation is conducted based on a simplified cavity-bogie model,including five cases with different inclination angles of the front and rear walls of the cavity.By comparing and analyzing the flow field and acoustic results of the five cases,the influence of the regularity and mechanism of the bogie cavity end wall inclination on the flow field and the aerodynamic noise of the bogie region are revealed.Then,the noise reduction strategy determined by the results of the simplified cavity-bogie model is applied to a three-car marshaling train model to verify its effectiveness when applied to the real train.The results reveal that the forward inclination of the cavity front wall enlarges the influence area of shear vortex structures formed at the leading edge of the cavity and intensifies the interaction between the vortex structures and the front wheelset,frontmotor,and front gearbox,resulting in the increase of the aerodynamic noise generated by the bogie itself.The backward inclination of the cavity rear wall is conducive to guiding the vortex structures flow out of the cavity and weakening the interaction between the shear vortex structures and the cavity rear wall,leading to the reduction of the aerodynamic noise generated by the bogie cavity.Inclining the rear end wall of the foremost bogie cavity of the head car is a feasible aerodynamic noise reduction measure for high-speed trains.展开更多
Prevention and control measures of spontaneous combustion of coal and gas accumulation in a goaf require an accurate description of its gas flow state.However,the commonly used fluid dynamics in porous media is not su...Prevention and control measures of spontaneous combustion of coal and gas accumulation in a goaf require an accurate description of its gas flow state.However,the commonly used fluid dynamics in porous media is not suitable for the new-born goaf with fracture cavity combination,multi-scale,and large blocks.In this study,we propose a cavity flow algorithm to accurately describe the gas flow state in the new-born goaf.The genetic algorithm(GA)is used to randomly generate the binary matrix of a goaf caving shape.The difference between the gas flow state calculated by the lattice Boltzmann method(LBM)and the measured data at the boundary or internal measuring points of the real goaf is taken as the GA fitness value,and the real goaf caving shape and the gas flow state are quickly addressed by GA.The experimental model of new-born goaf is established,and the laser Doppler anemometry(LDA)experiment is carried out.The results show that the Jaccard similarity coefficient between the reconstructed caving shape and the real caving shape is 0.7473,the mean square error between the calculated wind speed and the LDA-measured value is 0.0244,and the R2 coefficient is 0.8986,which verify the feasibility of the algorithm.展开更多
Heat and mass transfer of a circular-shaped porous moist object inside a two-dimensional triangle cavity is investigated by using finite element method.The porous object is considered to be a moist food sample,located...Heat and mass transfer of a circular-shaped porous moist object inside a two-dimensional triangle cavity is investigated by using finite element method.The porous object is considered to be a moist food sample,located in the middle of the cavity with inlet and outlet ports with different configurations of inlet/outlet ports.Convective drying performance is numerically assessed for different values of Reynolds numbers(between 50 and 250),dry air inlet temperature(between 40 and 80℃)and different locations of the port.It is observed that changing the port locations has significant impacts on the flow recirculaitons inside the triangular chamber while convective drying performance is highly affected.The moisture content reduces with longer time and for higher Reynolds number(Re)values.Case P4 where inlet and outlet ports are in the middle of the walls provides the most effective configuration in terms of convective drying performance while the worst case is seen for P1 case where inlet and outlet are closer to the corners of the chamber.There is a 192% difference between the moisture reduction of these two cases at Re=250,T=80℃ and t=120 min.展开更多
By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improv...By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number (Re) flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.展开更多
Supersonic flows past two-dimensional cavities with/without control are investigated by the direct numerical simulation (DNS). For an uncontrolled cavity, as the thickness of the boundary layer declines, transition ...Supersonic flows past two-dimensional cavities with/without control are investigated by the direct numerical simulation (DNS). For an uncontrolled cavity, as the thickness of the boundary layer declines, transition of the dominant mode from the steady mode to the Rossiter Ⅱ mode and then to the Rossiter III mode is observed due to the change of vortex-corner interactions. Meanwhile, a low frequency mode appears. However, the wake mode observed in a subsonic cavity flow is absent in the current simulation. The oscillation frequencies obtained from a global dynamic mode decomposition (DMD) approach are consistent with the local power spectral density (PSD) analysis. The dominant mode transition is clearly shown by the dynamic modes obtained from the DMD. A passive control technique of substituting the cavity trailing edge with a quarter-circle is studied. As the effective cavity length increases, the dominant mode transition from the Rossiter Ⅱ mode to the Rossiter Ⅲ mode occurs. With the control, the pressure oscillations are reduced significantly. The interaction of the shear layer and the recirculation zone is greatly weakened, combined with weaker shear layer instability, responsible for the suppression of pressure oscillations. Moreover, active control using steady subsonic mass injection upstream of a cavity leading edge can stabilize the flow.展开更多
The reduced-order model (ROM) for the two-dimensional supersonic cavity flow based on proper orthogonal decomposition (POD) and Galerkin projection is investigated. Presently, popular ROMs in cavity flows are base...The reduced-order model (ROM) for the two-dimensional supersonic cavity flow based on proper orthogonal decomposition (POD) and Galerkin projection is investigated. Presently, popular ROMs in cavity flows are based on an isentropic assumption, valid only for flows at low or moderate Mach numbers. A new ROM is constructed involving primitive variables of the fully compressible Navier-Stokes (N-S) equations, which is suitable for flows at high Mach numbers. Compared with the direct numerical simulation (DNS) results, the proposed model predicts flow dynamics (e.g., dominant frequency and amplitude) accurately for supersonic cavity flows, and is robust. The comparison between the present transient flow fields and those of the DNS shows that the proposed ROM can capture self-sustained oscillations of a shear layer. In addition, the present model reduction method can be easily extended to other supersonic flows.展开更多
Flow induced oscillation in a cavity and its suppression by means of acoustic excitation were studied both experimentally and numerically. In the experiment it was found that with the leading edge pure tone excitation...Flow induced oscillation in a cavity and its suppression by means of acoustic excitation were studied both experimentally and numerically. In the experiment it was found that with the leading edge pure tone excitation at some frequencies and intensities. the flow-induced oscillation in the cavity could be greatly suppressed. Cavity flows both with and without acoustic excitation were studied by solving the 2-D time-dependent Reynolds averaged Navier Stokes equations using explicit predictor-corrector difference algorithm of MacCormack. Effects of turbulence were simulated via Cebeci-Smith turbulence mode with relaxation modification. The computational and experimental results are compared. and good agreement is obtained.展开更多
To date,there are very few studies on the transition beyond second Hopf bifurcation in a lid-driven square cavity,due to the difficulties in theoretical analysis and numerical simulations.In this paper,we study the ch...To date,there are very few studies on the transition beyond second Hopf bifurcation in a lid-driven square cavity,due to the difficulties in theoretical analysis and numerical simulations.In this paper,we study the characteristics of the third Hopf bifurcation in a driven square cavity by applying a consistent fourth-order compact finite difference scheme rectently developed by us.We numerically identify the critical Reynolds number of the third Hopf bifurcation located in the interval of(13944.7021,13946.5333)by the method of bisection.Through Fourier analysis,it is discovered that the flow becomes chaotic with a characteristic of period-doubling bifurcation when the Reynolds number is beyond the third bifurcation critical interval.Nonlinear time series analysis further ascertains the flow chaotic behaviors via the phase diagram,Kolmogorov entropy and maximal Lyapunov exponent.The phase diagram changes interestingly from a closed curve with self-intersection to an unclosed curve and the attractor eventually becomes strange when the flow becomes chaotic.展开更多
A streamline upwind/Petrov-Galerkin (SUPG) finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for-...A streamline upwind/Petrov-Galerkin (SUPG) finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for- mulation of momentum equations. Using the penalty function method, the continuity equation is simplified and the pres- sure of the momentum equations is eliminated. The lid-driven cavity flow problem is solved using the present model. It is shown that steady flow simulations are computable up to Re = 27500, and the present results agree well with previous solutions. Tabulated results for the properties of the primary vortex are also provided for benchmarking purposes.展开更多
Based on finite volume method, subsonic and transonic flow in 3-D cavity of different length-to-depth ratios are numerically investigated by using Navier-Stokes equations with k-ε SST two-equation turbulence model an...Based on finite volume method, subsonic and transonic flow in 3-D cavity of different length-to-depth ratios are numerically investigated by using Navier-Stokes equations with k-ε SST two-equation turbulence model and coupled implicit algorithm. The cavity streamline patterns and the static pressure coefficient distributions on the cavity floor are shown, and the flow characteristics of the cavity and the floor pressure coefficient distributions are analyzed. Numerical results demon- strate that the flow characteristics of the cavity at subsonic and transonic speeds are different from that of supersonic ones; length-to-depth ratio is the main factor that affects the flow characteristics of the cavity at subsonic and transonic number has a neglectable effect on the cavity floor pressure distributions. speeds and causes changes of the cavity flow type; Mach cavity flow fields at subsonic and transonic speeds and the展开更多
The fractional volumetric lattice Boltzmann method with much better stability was used to simulate two-chmensional cavity flows. Because the effective viscosity was reduced by the fraction factor, it is very effective...The fractional volumetric lattice Boltzmann method with much better stability was used to simulate two-chmensional cavity flows. Because the effective viscosity was reduced by the fraction factor, it is very effective for simulating high Reynolds number flows. Simulations were carried out on a uniform grids system. The stream lines and the velocity profiles obtained from the simulations agree well with the standard lattice Boltzmann method simulations. Comparisons of detailed flow patterns with other studies via location of vortex centers are also satisfactory.展开更多
Employing Rayleigh's method, the collapse of a vaporous bubble in an incompressible liquid with surface tension is analysed. The expressions of time versus radius, bubble-wall velocity and pressure developed at co...Employing Rayleigh's method, the collapse of a vaporous bubble in an incompressible liquid with surface tension is analysed. The expressions of time versus radius, bubble-wall velocity and pressure developed at collapse are thus introduced.Finally, the numerical solution of velocity and pressure field in the liquid surrounding the cavity is also given.展开更多
Numerical simulations are used to investigate the self-sustained oscillating flows past an open cavity. The two-dimensional incompressible Navier-Stokes equations are solved directly by using the finite difference met...Numerical simulations are used to investigate the self-sustained oscillating flows past an open cavity. The two-dimensional incompressible Navier-Stokes equations are solved directly by using the finite difference method for cavities with an upstream laminar boundary layer. A series of simulations are performed for a variety of cavity length-to-depth ratio. The results show the switching among some flow modes including non-oscillation mode, shear layer mode and wake mode. The variation of the Strouhal number is in favorable agreement with available experimental data. The results of flow fields in the cavity reveal the relationship between the cavity shear layer oscillation modes and recirculating vortices in the cavity.展开更多
Physics-informed neural networks(PINNs)are proved methods that are effective in solving some strongly nonlinear partial differential equations(PDEs),e.g.,Navier-Stokes equations,with a small amount of boundary or inte...Physics-informed neural networks(PINNs)are proved methods that are effective in solving some strongly nonlinear partial differential equations(PDEs),e.g.,Navier-Stokes equations,with a small amount of boundary or interior data.However,the feasibility of applying PINNs to the flow at moderate or high Reynolds numbers has rarely been reported.The present paper proposes an artificial viscosity(AV)-based PINN for solving the forward and inverse flow problems.Specifically,the AV used in PINNs is inspired by the entropy viscosity method developed in conventional computational fluid dynamics(CFD)to stabilize the simulation of flow at high Reynolds numbers.The newly developed PINN is used to solve the forward problem of the two-dimensional steady cavity flow at Re=1000 and the inverse problem derived from two-dimensional film boiling.The results show that the AV augmented PINN can solve both problems with good accuracy and substantially reduce the inference errors in the forward problem.展开更多
Nowadays,most researchers focus on the cavity shedding mechanisms of unsteady cavitating flows over different objects,such as 2D/3D hydrofoils,venturi-type section,axisymmetric bodies with different headforms,and so o...Nowadays,most researchers focus on the cavity shedding mechanisms of unsteady cavitating flows over different objects,such as 2D/3D hydrofoils,venturi-type section,axisymmetric bodies with different headforms,and so on.But few of them pay attention to the differences of cavity shedding modality under different cavitation numbers in unsteady cavitating flows over the same object.In the present study,two kinds of shedding patterns are investigated experimentally.A high speed camera system is used to observe the cavitating flows over an axisymmetric blunt body and the velocity fields are measured by a particle image velocimetry(PIV)technique in a water tunnel for different cavitation conditions.The U-type cavitating vortex shedding is observed in unsteady cavitating flows.When the cavitation number is 0.7,there is a large scale cavity rolling up and shedding,which cause the instability and dramatic fluctuation of the flows,while at cavitation number of 0.6,the detached cavities can be conjunct with the attached part to induce the break-off behavior again at the tail of the attached cavity,as a result,the final shedding is in the form of small scale cavity and keeps a relatively steady flow field.It is also found that the interaction between the re-entrant flow and the attached cavity plays an important role in the unsteady cavity shedding modality.When the attached cavity scale is insufficient to overcome the re-entrant flow,it deserves the large cavity rolling up and shedding just as that at cavitation number of 0.7.Otherwise,the re-entrant flow is defeated by large enough cavity to induce the cavity-combined process and small scale cavity vortexes shedding just as that of the cavitation number of0.6.This research shows the details of two different cavity shedding modalities which is worthful and meaningful for the further study of unsteady cavitation.展开更多
Multiple steady solutions and hysteresis phenomenon in the square cavity flows driven by the surface with antisymmetric velocity profile are investigated by numerical simulation and bifurcation analysis.A high order s...Multiple steady solutions and hysteresis phenomenon in the square cavity flows driven by the surface with antisymmetric velocity profile are investigated by numerical simulation and bifurcation analysis.A high order spectral element method with the matrix-free pseudo-arclength technique is used for the steady-state solution and numerical continuation.The complex flow patterns beyond the symmetry-breaking at Re≈320 are presented by a bifurcation diagram for Re<2500.The results of stable symmetric and asymmetric solutions are consistent with those reported in literature,and a new unstable asymmetric branch is obtained besides the stable branches.A novel hysteresis phenomenon is observed in the range of 2208<Re<2262,where two pairs of stable and two pairs of unstable asymmetric steady solutions beyond the stable symmetric state coexist.The vortices near the sidewall appear when the Reynolds number increases,which correspond to the bifurcation of topology structure,but not the bifurcation of Navier-Stokes equations.The hysteresis is proposed to be the result of the combined mechanisms of the competition and coalescence of secondary vortices.展开更多
When the variable geometry hypersonic inlet is sealed with ceramic wafers,the cavity flows inside the sealing chamber can be affected by the boundary layer near the side wall.To study the influence of the boundary lay...When the variable geometry hypersonic inlet is sealed with ceramic wafers,the cavity flows inside the sealing chamber can be affected by the boundary layer near the side wall.To study the influence of the boundary layer thickness near the side wall on the flow and leakage characteristics in sealing chamber,the numerical calculation of the cavity flow in the sealing chamber under different inflow boundary layer thicknesses is carried out.The results show that three-dimensional cavity flow structures are close to being asymmetric,and the entrance pressure of the leakage path can also be affected by asymmetry;with the increase of the thickness of the boundary layer,the pressure at the cavity floor and the seal entrance decreases.Finally,the existing leakage prediction model is modified according to the distribution rule of the cavity floor and the flow properties in the leakage path.展开更多
As an integral part of the internal air system of aero-engines,the axial throughflow of the cooling air can interact with the cavity flow between the rotating compressor disks,forming a threedimensional,unsteady,and u...As an integral part of the internal air system of aero-engines,the axial throughflow of the cooling air can interact with the cavity flow between the rotating compressor disks,forming a threedimensional,unsteady,and unstable flow field.The flow characteristics in an engine-like rotating multi-stage cavity with throughflow were investigated using particle image velocimetry,flow visualization technology and three-dimensional unsteady Reynolds-Averaged Navier-Stokes (RANS)simulations.The focus of current research was to understand the distribution of the mean swirl ratio and its variation with a wide range of non-dimensional parameters in the co-rotating cavity with high inlet pre-swirl axial throughflow.The maximum axial Reynolds number and rotational Reynolds numbers could reach 4.41×10^(4)and 1.24×10^(6),respectively.The velocity measurement results indicate that the mean swirl ratio is greater than 1 and decreases with an increase in the radial position.The flow structure is dominated by the Rossby number,and two different flow patterns (flow penetration and flow stratification) are identified and confirmed by flow visualization images.In the absence of buoyancy,the flow penetration caused by the precession of the throughflow makes it easier for the throughflow to reach a high radius region.Satisfactory consistency of results between measurements and numerical calculations is obtained.This study provides a theoretical basis and data support for toroidal vortex breakdown,which is of practical significance for the design of high-pressure compressor cavities.展开更多
Experimental and computational analysis has been carried out by many researchers on supersonic cavity flow,but detailed analysis based on Rossiter’s model still requires some insight.In the current study an open rect...Experimental and computational analysis has been carried out by many researchers on supersonic cavity flow,but detailed analysis based on Rossiter’s model still requires some insight.In the current study an open rectangular cavity with a length to depth ratio of 2(L/D=2)and Mach number at the inlet as 1.71,was considered as a baseline configuration for experimental analysis.Statistical techniques such as power spectral density(PSD),correlation,and overall sound pressure level(OASPL)were carried out on the unsteady pressure data,to analyze the aero-acoustic flow physics.High-speed schlieren images were processed to obtain spatially coherent modes by proper orthogonal decomposition(POD).The analysis was extended for different dimensions of subcavities on the aft,floor,and front wall.This detailed analysis of these configurations with different dimensions and combinations revealed the various flow features and mode frequencies in supersonic cavity.As the front wall subcavity act as a passive control device,reducing the overall sound pressure level inside the cavity whereas,the aft wall subcavity acts as a passive resonator with distinct harmonic fluidresonant modes,a similar phenomenon was observed for floor subcavity at different locations.A novel method was employed to analyze Rossiter’s model and its applicability in estimating experimental modes was verified,as it accurately predicted the dominant frequencies with a maximum of 2.74% uncertainty among all the configurations.展开更多
文摘An experimental study of cavity oscillating flow carried out on subsonic wall jet facilities in an anechoic room is summarized. The jet exit Mach number range is from 0.2 to 0.8. The effects of the flow Mach number ( Ma ) and the cavity depth ( D ) on the oscillation are studied. It is found that for L/D =4, (shallow cavity), the oscillation is mainly due to the self exciting of the free shear layer above the cavity opening, for L/D =2, the acoustic resonance is responsible for the oscillation. Preliminary tests are performed to study the suppression effect of the leading edge tone excitation on cavity flow, and considerable reduction of oscillation has been achieved when Ma ≤0.6.
基金supported by National Natural Science Foundation of China(12172308)National Key Research and Development Program of China(2020YFA0710902).
文摘Combining the detached eddy simulation(DES)method and Ffowcs Williams-Hawkings(FW-H)equation,the effect of bogie cavity end wall inclination on the flow field and aerodynamic noise in the bogie region is numerically studied.First,the simulation is conducted based on a simplified cavity-bogie model,including five cases with different inclination angles of the front and rear walls of the cavity.By comparing and analyzing the flow field and acoustic results of the five cases,the influence of the regularity and mechanism of the bogie cavity end wall inclination on the flow field and the aerodynamic noise of the bogie region are revealed.Then,the noise reduction strategy determined by the results of the simplified cavity-bogie model is applied to a three-car marshaling train model to verify its effectiveness when applied to the real train.The results reveal that the forward inclination of the cavity front wall enlarges the influence area of shear vortex structures formed at the leading edge of the cavity and intensifies the interaction between the vortex structures and the front wheelset,frontmotor,and front gearbox,resulting in the increase of the aerodynamic noise generated by the bogie itself.The backward inclination of the cavity rear wall is conducive to guiding the vortex structures flow out of the cavity and weakening the interaction between the shear vortex structures and the cavity rear wall,leading to the reduction of the aerodynamic noise generated by the bogie cavity.Inclining the rear end wall of the foremost bogie cavity of the head car is a feasible aerodynamic noise reduction measure for high-speed trains.
基金This work was supported by the Natural Science Foundation of China(Nos.51774169 and 51574142)the National Key Research and Development Program of China(No.2017YFC0804401).
文摘Prevention and control measures of spontaneous combustion of coal and gas accumulation in a goaf require an accurate description of its gas flow state.However,the commonly used fluid dynamics in porous media is not suitable for the new-born goaf with fracture cavity combination,multi-scale,and large blocks.In this study,we propose a cavity flow algorithm to accurately describe the gas flow state in the new-born goaf.The genetic algorithm(GA)is used to randomly generate the binary matrix of a goaf caving shape.The difference between the gas flow state calculated by the lattice Boltzmann method(LBM)and the measured data at the boundary or internal measuring points of the real goaf is taken as the GA fitness value,and the real goaf caving shape and the gas flow state are quickly addressed by GA.The experimental model of new-born goaf is established,and the laser Doppler anemometry(LDA)experiment is carried out.The results show that the Jaccard similarity coefficient between the reconstructed caving shape and the real caving shape is 0.7473,the mean square error between the calculated wind speed and the LDA-measured value is 0.0244,and the R2 coefficient is 0.8986,which verify the feasibility of the algorithm.
文摘Heat and mass transfer of a circular-shaped porous moist object inside a two-dimensional triangle cavity is investigated by using finite element method.The porous object is considered to be a moist food sample,located in the middle of the cavity with inlet and outlet ports with different configurations of inlet/outlet ports.Convective drying performance is numerically assessed for different values of Reynolds numbers(between 50 and 250),dry air inlet temperature(between 40 and 80℃)and different locations of the port.It is observed that changing the port locations has significant impacts on the flow recirculaitons inside the triangular chamber while convective drying performance is highly affected.The moisture content reduces with longer time and for higher Reynolds number(Re)values.Case P4 where inlet and outlet ports are in the middle of the walls provides the most effective configuration in terms of convective drying performance while the worst case is seen for P1 case where inlet and outlet are closer to the corners of the chamber.There is a 192% difference between the moisture reduction of these two cases at Re=250,T=80℃ and t=120 min.
基金Project supported by the National Natural Science Foundation of China (Grant No 70271069).
文摘By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number (Re) flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.
基金supported by the National Natural Science Foundation of China(Nos.11232011 and11402262)the 111 Project of China(No.B07033)+1 种基金the China Postdoctoral Science Foundation(No.2014M561833)the Fundamental Research Funds for the Central Universities
文摘Supersonic flows past two-dimensional cavities with/without control are investigated by the direct numerical simulation (DNS). For an uncontrolled cavity, as the thickness of the boundary layer declines, transition of the dominant mode from the steady mode to the Rossiter Ⅱ mode and then to the Rossiter III mode is observed due to the change of vortex-corner interactions. Meanwhile, a low frequency mode appears. However, the wake mode observed in a subsonic cavity flow is absent in the current simulation. The oscillation frequencies obtained from a global dynamic mode decomposition (DMD) approach are consistent with the local power spectral density (PSD) analysis. The dominant mode transition is clearly shown by the dynamic modes obtained from the DMD. A passive control technique of substituting the cavity trailing edge with a quarter-circle is studied. As the effective cavity length increases, the dominant mode transition from the Rossiter Ⅱ mode to the Rossiter Ⅲ mode occurs. With the control, the pressure oscillations are reduced significantly. The interaction of the shear layer and the recirculation zone is greatly weakened, combined with weaker shear layer instability, responsible for the suppression of pressure oscillations. Moreover, active control using steady subsonic mass injection upstream of a cavity leading edge can stabilize the flow.
基金Project supported by the National Natural Science Foundation of China(Nos.11232011,11402262,11572314,and 11621202)
文摘The reduced-order model (ROM) for the two-dimensional supersonic cavity flow based on proper orthogonal decomposition (POD) and Galerkin projection is investigated. Presently, popular ROMs in cavity flows are based on an isentropic assumption, valid only for flows at low or moderate Mach numbers. A new ROM is constructed involving primitive variables of the fully compressible Navier-Stokes (N-S) equations, which is suitable for flows at high Mach numbers. Compared with the direct numerical simulation (DNS) results, the proposed model predicts flow dynamics (e.g., dominant frequency and amplitude) accurately for supersonic cavity flows, and is robust. The comparison between the present transient flow fields and those of the DNS shows that the proposed ROM can capture self-sustained oscillations of a shear layer. In addition, the present model reduction method can be easily extended to other supersonic flows.
文摘Flow induced oscillation in a cavity and its suppression by means of acoustic excitation were studied both experimentally and numerically. In the experiment it was found that with the leading edge pure tone excitation at some frequencies and intensities. the flow-induced oscillation in the cavity could be greatly suppressed. Cavity flows both with and without acoustic excitation were studied by solving the 2-D time-dependent Reynolds averaged Navier Stokes equations using explicit predictor-corrector difference algorithm of MacCormack. Effects of turbulence were simulated via Cebeci-Smith turbulence mode with relaxation modification. The computational and experimental results are compared. and good agreement is obtained.
基金Project supported by the National Natural Science Foundation of China(Grant No.12162001)the Natural Science Foundation of Ningxia(Grant No.2019AAC03129)the Construction Project of First-Class Disciplines in Ningxia Higher Education(Grant No.NXYLXK2017B09)。
文摘To date,there are very few studies on the transition beyond second Hopf bifurcation in a lid-driven square cavity,due to the difficulties in theoretical analysis and numerical simulations.In this paper,we study the characteristics of the third Hopf bifurcation in a driven square cavity by applying a consistent fourth-order compact finite difference scheme rectently developed by us.We numerically identify the critical Reynolds number of the third Hopf bifurcation located in the interval of(13944.7021,13946.5333)by the method of bisection.Through Fourier analysis,it is discovered that the flow becomes chaotic with a characteristic of period-doubling bifurcation when the Reynolds number is beyond the third bifurcation critical interval.Nonlinear time series analysis further ascertains the flow chaotic behaviors via the phase diagram,Kolmogorov entropy and maximal Lyapunov exponent.The phase diagram changes interestingly from a closed curve with self-intersection to an unclosed curve and the attractor eventually becomes strange when the flow becomes chaotic.
基金the National Natural Science Foundation of China (Grants 41372301 and 51349011)the Preeminent Youth Talent Project of Southwest University of Science and Technology (Grant 13zx9109)
文摘A streamline upwind/Petrov-Galerkin (SUPG) finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for- mulation of momentum equations. Using the penalty function method, the continuity equation is simplified and the pres- sure of the momentum equations is eliminated. The lid-driven cavity flow problem is solved using the present model. It is shown that steady flow simulations are computable up to Re = 27500, and the present results agree well with previous solutions. Tabulated results for the properties of the primary vortex are also provided for benchmarking purposes.
文摘Based on finite volume method, subsonic and transonic flow in 3-D cavity of different length-to-depth ratios are numerically investigated by using Navier-Stokes equations with k-ε SST two-equation turbulence model and coupled implicit algorithm. The cavity streamline patterns and the static pressure coefficient distributions on the cavity floor are shown, and the flow characteristics of the cavity and the floor pressure coefficient distributions are analyzed. Numerical results demon- strate that the flow characteristics of the cavity at subsonic and transonic speeds are different from that of supersonic ones; length-to-depth ratio is the main factor that affects the flow characteristics of the cavity at subsonic and transonic number has a neglectable effect on the cavity floor pressure distributions. speeds and causes changes of the cavity flow type; Mach cavity flow fields at subsonic and transonic speeds and the
文摘The fractional volumetric lattice Boltzmann method with much better stability was used to simulate two-chmensional cavity flows. Because the effective viscosity was reduced by the fraction factor, it is very effective for simulating high Reynolds number flows. Simulations were carried out on a uniform grids system. The stream lines and the velocity profiles obtained from the simulations agree well with the standard lattice Boltzmann method simulations. Comparisons of detailed flow patterns with other studies via location of vortex centers are also satisfactory.
文摘Employing Rayleigh's method, the collapse of a vaporous bubble in an incompressible liquid with surface tension is analysed. The expressions of time versus radius, bubble-wall velocity and pressure developed at collapse are thus introduced.Finally, the numerical solution of velocity and pressure field in the liquid surrounding the cavity is also given.
文摘Numerical simulations are used to investigate the self-sustained oscillating flows past an open cavity. The two-dimensional incompressible Navier-Stokes equations are solved directly by using the finite difference method for cavities with an upstream laminar boundary layer. A series of simulations are performed for a variety of cavity length-to-depth ratio. The results show the switching among some flow modes including non-oscillation mode, shear layer mode and wake mode. The variation of the Strouhal number is in favorable agreement with available experimental data. The results of flow fields in the cavity reveal the relationship between the cavity shear layer oscillation modes and recirculating vortices in the cavity.
基金Project supported by the Fundamental Research Funds for the Central Universities of China(No.DUT21RC(3)063)the National Natural Science Foundation of China(No.51720105007)the Baidu Foundation(No.ghfund202202014542)。
文摘Physics-informed neural networks(PINNs)are proved methods that are effective in solving some strongly nonlinear partial differential equations(PDEs),e.g.,Navier-Stokes equations,with a small amount of boundary or interior data.However,the feasibility of applying PINNs to the flow at moderate or high Reynolds numbers has rarely been reported.The present paper proposes an artificial viscosity(AV)-based PINN for solving the forward and inverse flow problems.Specifically,the AV used in PINNs is inspired by the entropy viscosity method developed in conventional computational fluid dynamics(CFD)to stabilize the simulation of flow at high Reynolds numbers.The newly developed PINN is used to solve the forward problem of the two-dimensional steady cavity flow at Re=1000 and the inverse problem derived from two-dimensional film boiling.The results show that the AV augmented PINN can solve both problems with good accuracy and substantially reduce the inference errors in the forward problem.
基金Supported by National Natural Science Foundation of China(Grant Nos.51209004,11172040)Beijing Municipal Natural Science Foundation of China(Grant No.3144043)
文摘Nowadays,most researchers focus on the cavity shedding mechanisms of unsteady cavitating flows over different objects,such as 2D/3D hydrofoils,venturi-type section,axisymmetric bodies with different headforms,and so on.But few of them pay attention to the differences of cavity shedding modality under different cavitation numbers in unsteady cavitating flows over the same object.In the present study,two kinds of shedding patterns are investigated experimentally.A high speed camera system is used to observe the cavitating flows over an axisymmetric blunt body and the velocity fields are measured by a particle image velocimetry(PIV)technique in a water tunnel for different cavitation conditions.The U-type cavitating vortex shedding is observed in unsteady cavitating flows.When the cavitation number is 0.7,there is a large scale cavity rolling up and shedding,which cause the instability and dramatic fluctuation of the flows,while at cavitation number of 0.6,the detached cavities can be conjunct with the attached part to induce the break-off behavior again at the tail of the attached cavity,as a result,the final shedding is in the form of small scale cavity and keeps a relatively steady flow field.It is also found that the interaction between the re-entrant flow and the attached cavity plays an important role in the unsteady cavity shedding modality.When the attached cavity scale is insufficient to overcome the re-entrant flow,it deserves the large cavity rolling up and shedding just as that at cavitation number of 0.7.Otherwise,the re-entrant flow is defeated by large enough cavity to induce the cavity-combined process and small scale cavity vortexes shedding just as that of the cavitation number of0.6.This research shows the details of two different cavity shedding modalities which is worthful and meaningful for the further study of unsteady cavitation.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11902043 and 11772065)the Science Challenge Project(Grant No.TZ2016001).
文摘Multiple steady solutions and hysteresis phenomenon in the square cavity flows driven by the surface with antisymmetric velocity profile are investigated by numerical simulation and bifurcation analysis.A high order spectral element method with the matrix-free pseudo-arclength technique is used for the steady-state solution and numerical continuation.The complex flow patterns beyond the symmetry-breaking at Re≈320 are presented by a bifurcation diagram for Re<2500.The results of stable symmetric and asymmetric solutions are consistent with those reported in literature,and a new unstable asymmetric branch is obtained besides the stable branches.A novel hysteresis phenomenon is observed in the range of 2208<Re<2262,where two pairs of stable and two pairs of unstable asymmetric steady solutions beyond the stable symmetric state coexist.The vortices near the sidewall appear when the Reynolds number increases,which correspond to the bifurcation of topology structure,but not the bifurcation of Navier-Stokes equations.The hysteresis is proposed to be the result of the combined mechanisms of the competition and coalescence of secondary vortices.
基金supported by the Opening Foundation of National State Key Laboratory of High Temperature Gas Dynamics(No.2021KF07)。
文摘When the variable geometry hypersonic inlet is sealed with ceramic wafers,the cavity flows inside the sealing chamber can be affected by the boundary layer near the side wall.To study the influence of the boundary layer thickness near the side wall on the flow and leakage characteristics in sealing chamber,the numerical calculation of the cavity flow in the sealing chamber under different inflow boundary layer thicknesses is carried out.The results show that three-dimensional cavity flow structures are close to being asymmetric,and the entrance pressure of the leakage path can also be affected by asymmetry;with the increase of the thickness of the boundary layer,the pressure at the cavity floor and the seal entrance decreases.Finally,the existing leakage prediction model is modified according to the distribution rule of the cavity floor and the flow properties in the leakage path.
文摘As an integral part of the internal air system of aero-engines,the axial throughflow of the cooling air can interact with the cavity flow between the rotating compressor disks,forming a threedimensional,unsteady,and unstable flow field.The flow characteristics in an engine-like rotating multi-stage cavity with throughflow were investigated using particle image velocimetry,flow visualization technology and three-dimensional unsteady Reynolds-Averaged Navier-Stokes (RANS)simulations.The focus of current research was to understand the distribution of the mean swirl ratio and its variation with a wide range of non-dimensional parameters in the co-rotating cavity with high inlet pre-swirl axial throughflow.The maximum axial Reynolds number and rotational Reynolds numbers could reach 4.41×10^(4)and 1.24×10^(6),respectively.The velocity measurement results indicate that the mean swirl ratio is greater than 1 and decreases with an increase in the radial position.The flow structure is dominated by the Rossby number,and two different flow patterns (flow penetration and flow stratification) are identified and confirmed by flow visualization images.In the absence of buoyancy,the flow penetration caused by the precession of the throughflow makes it easier for the throughflow to reach a high radius region.Satisfactory consistency of results between measurements and numerical calculations is obtained.This study provides a theoretical basis and data support for toroidal vortex breakdown,which is of practical significance for the design of high-pressure compressor cavities.
基金the support provided by the Department of Aerospace Engineering,IIST,Trivandrum,India,in carrying out this research work.
文摘Experimental and computational analysis has been carried out by many researchers on supersonic cavity flow,but detailed analysis based on Rossiter’s model still requires some insight.In the current study an open rectangular cavity with a length to depth ratio of 2(L/D=2)and Mach number at the inlet as 1.71,was considered as a baseline configuration for experimental analysis.Statistical techniques such as power spectral density(PSD),correlation,and overall sound pressure level(OASPL)were carried out on the unsteady pressure data,to analyze the aero-acoustic flow physics.High-speed schlieren images were processed to obtain spatially coherent modes by proper orthogonal decomposition(POD).The analysis was extended for different dimensions of subcavities on the aft,floor,and front wall.This detailed analysis of these configurations with different dimensions and combinations revealed the various flow features and mode frequencies in supersonic cavity.As the front wall subcavity act as a passive control device,reducing the overall sound pressure level inside the cavity whereas,the aft wall subcavity acts as a passive resonator with distinct harmonic fluidresonant modes,a similar phenomenon was observed for floor subcavity at different locations.A novel method was employed to analyze Rossiter’s model and its applicability in estimating experimental modes was verified,as it accurately predicted the dominant frequencies with a maximum of 2.74% uncertainty among all the configurations.