When better fuel-air mixing in the combustion chamber or a reduction in base drag are required in vehicles,rockets,and aeroplanes,the base pressure control is activated.Controlling the base pressure and drag is necess...When better fuel-air mixing in the combustion chamber or a reduction in base drag are required in vehicles,rockets,and aeroplanes,the base pressure control is activated.Controlling the base pressure and drag is necessary in both scenarios.In this work,semi-circular ribs with varying diameters(2,4,and 6 mm)positioned at six distinct positions(0.5D,1D,1.5D,2D,3D,and 4D)inside a square duct with a side of 15 mm are proposed as an efficient way to apply the passive control technique.In-depth research is done on optimising rib size for various rib sites.According to this study,the base pressure rises as rib height increases.Furthermore,the optimal location for the semi-circular ribs with a diameter of 2 mm is at 0.5D.The 1D location appears to be optimal for the 4 mm size as well.For the 6 mm size,however,the 4D position fills this function.展开更多
The purpose of this work is to shed light on the effect of the pivot position on the surface pressure distribution over a 3D wing in different flight conditions.The study is intended to support the design and developm...The purpose of this work is to shed light on the effect of the pivot position on the surface pressure distribution over a 3D wing in different flight conditions.The study is intended to support the design and development of aerospace vehicles where stability analysis,performance optimization,and aircraft design are of primary importance.The following parameters are considered:Mach numbers(M)of 1.3,1.8,2.3,2.8,3.3,and 3.8,angle of incidence(θ)in the range from 5°to 25°,pivot position from h=0.2 to 1.The results of the CFD numerical simulations match available analytical data,thereby providing evidence for the reliability of the used approach.The findings provide valuable insights into the relationship between the surface pressure distribution,the Mach number and the angle of incidence.展开更多
By analyzing the characteristics of low Mach number perfect gas flows, a novel Slightly Compressible Model (SCM) for low Mach number perect gas flows is derived. In view of numerical calculations, this model is proved...By analyzing the characteristics of low Mach number perfect gas flows, a novel Slightly Compressible Model (SCM) for low Mach number perect gas flows is derived. In view of numerical calculations, this model is proved very efficient, for it is kept within thep-v frame but does not have to satisfy the time consuming divergence-free condition in order to get the incompressible Navier-Stokes equation solution. Writing the equations in the form of conservation laws, we have derived the characteristic systems which are necessary for numerical calculations. A cell-centered finite-volume method with flux difference upwind-biased schemes is used for the equation solutions and a new Exact Newton Relaxation (ENR) implicit method is developed. Various computed results are presented to validate the present model. Laminar flow solutions over a circular cylinder with wake developing and vortex shedding are presented. Results for inviscid flow over a sphere are compared in excellent agreement with the exact analytic incompressible solution. Three-dimensional viscous flow solutions over sphere and prolate spheroid are also calculated and compared well with experiments and other incompressible solutions. Finally, good convergent performances are shown for sphere viscous flows.展开更多
In the present study,the base pressure variations induced by the presence of a cavity,known to have a strong influence of the behaviour of supersonic projectiles,are investigated through numerical solution of the bala...In the present study,the base pressure variations induced by the presence of a cavity,known to have a strong influence of the behaviour of supersonic projectiles,are investigated through numerical solution of the balance equations for mass,momentum,and energy.An area ratio of four is considered and numerical simulations are carried out at Mach M=1.2,1.4,1.6,and 1.8 assuming no cavity or cavity locations 0.5D,1D,1.5D,and 2D.The inlet pressure of the nozzle is considered as a flow variable.The Taguchi method is also used,and the considered cases are then analyzed using a full factorial experimental design.The results show that the cavity is effective in increasing the base pressure for the conditions examined.For other nozzle pressure ratios,cavities do not lead to passive control due the change in the reattachment length.The distribution of wall pressure reveals that,in general,a cavity used to implement passive control of the base pressure does not adversely influence the flow pattern in the domain.展开更多
This paper is concerned with the low Mach number limit for the compressible Navier-Stokes equations in an exterior domain. We present here an approach based on Strichartz estimate defined on a non trapping exterior do...This paper is concerned with the low Mach number limit for the compressible Navier-Stokes equations in an exterior domain. We present here an approach based on Strichartz estimate defined on a non trapping exterior domain and we will be able to show the compactness and strong convergence of the velocity vector field.展开更多
This work aims to compute stability derivatives in the Newtonian limit in pitch when the Mach number tends to infinity.In such conditions,these stability derivatives depend on the Ogive’s shape and not the Mach numbe...This work aims to compute stability derivatives in the Newtonian limit in pitch when the Mach number tends to infinity.In such conditions,these stability derivatives depend on the Ogive’s shape and not the Mach number.Generally,the Mach number independence principle becomes effective from M=10 and above.The Ogive nose is obtained through a circular arc on the cone surface.Accordingly,the following arc slopes are consideredλ=5,10,15,−5,−10,and−15.It is found that the stability derivatives decrease due to the growth inλfrom 5 to 15 and vice versa.Forλ=5 and 10,the damping derivative declines with an increase inλfrom 5 to 10.Yet,for the damping derivatives,the minimum location remains at a pivot position,h=0.75 for large values ofλ.Hence,whenλ=−15,the damping derivatives are independent of the cone angles for most pivot positions except in the early twenty percent of the leading edge.展开更多
On super-sonic or trans-sonic planar cascade wind tunnel of free jet intermittent type, wind blowing experiments were performed on the typical sections of stator and rotor blades in the last stage of ultra-ultra-criti...On super-sonic or trans-sonic planar cascade wind tunnel of free jet intermittent type, wind blowing experiments were performed on the typical sections of stator and rotor blades in the last stage of ultra-ultra-critical steam turbine with extra-long blade of 1200mm. The influences of attack angle and Mach number on the aerodynamic performances of these sections of the blade profiles were verified, and their operating ranges were also specified.展开更多
We investigate the low Mach number limit for the isentropic compressible NavierStokes equations with a revised Maxwell's law(with Galilean invariance) in R^(3). By applying the uniform estimates of the error syste...We investigate the low Mach number limit for the isentropic compressible NavierStokes equations with a revised Maxwell's law(with Galilean invariance) in R^(3). By applying the uniform estimates of the error system, it is proven that the solutions of the isentropic Navier-Stokes equations with a revised Maxwell's law converge to that of the incompressible Navier-Stokes equations as the Mach number tends to zero. Moreover, the convergence rates are also obtained.展开更多
Considering droplet phenomena at low Mach numbers,large differences in the magnitude of the occurring characteristic waves are presented.As acoustic phenomena often play a minor role in such applications,classical exp...Considering droplet phenomena at low Mach numbers,large differences in the magnitude of the occurring characteristic waves are presented.As acoustic phenomena often play a minor role in such applications,classical explicit schemes which resolve these waves suffer from a very restrictive timestep restriction.In this work,a novel scheme based on a specific level set ghost fluid method and an implicit-explicit(IMEX)flux splitting is proposed to overcome this timestep restriction.A fully implicit narrow band around the sharp phase interface is combined with a splitting of the convective and acoustic phenomena away from the interface.In this part of the domain,the IMEX Runge-Kutta time discretization and the high order discontinuous Galerkin spectral element method are applied to achieve high accuracies in the bulk phases.It is shown that for low Mach numbers a significant gain in computational time can be achieved compared to a fully explicit method.Applica-tions to typical droplet dynamic phenomena validate the proposed method and illustrate its capabilities.展开更多
The interaction of a shock wave with a spherical helium bubble is investigated numerically by using the high- resolution piecewise parabolic method (PPM), in which the viscous and turbulence effects are both conside...The interaction of a shock wave with a spherical helium bubble is investigated numerically by using the high- resolution piecewise parabolic method (PPM), in which the viscous and turbulence effects are both considered. The bubble is of the same size and is accelerated by a planar shock of different Mach numbers (Ma). The re- suits of low Ma cases agree quantitatively with those of experiments [G. Layes, O. Le M4tayer. Phys. Fluids 19 (2007) 042105]. With the increase of Ma, the final geometry of the bubble becomes quite different, the com- pression ratio is highly raised, and the time-dependent mean bubble velocity is also influenced. The compression ratios measured can be well normalized when Ma is low, while less agreement has been achieved for high Ma cases. In addition, the mixedness between two fluids is enhanced greatly as Ma increases. Some existed scaling laws of these quantities for the shock wave strength cannot be directly applied to high Ma cases.展开更多
Mach number is a key metric in the evaluation of wind tunnel flow field performance.This complex process of wind tunnel test mainly has the problems of nonlinearity and time lag.In order to overcome the problems and c...Mach number is a key metric in the evaluation of wind tunnel flow field performance.This complex process of wind tunnel test mainly has the problems of nonlinearity and time lag.In order to overcome the problems and control the Mach number stability,this paper proposes a new method of Mach number prediction based on a nonlinear autoregressive exogenous-genetic algorithm-Elman(NARX-GA-Elman)model,which adopts NARX as the basic framework,determines the order of the input variables by using the false nearest neighbor(FNN),and uses the dynamic nonlinear network Elman to fit the model,and finally uses the global optimization algorithm GA to optimize the weight thresholds in the model to establish the Mach number prediction model with optimal performance under single working condition.By comparing with the traditional algorithm,the prediction accuracy of the model is improved by 61.5%,and the control accuracy is improved by 55.7%,which demonstrates that the model has very high prediction accuracy and good stability performance.展开更多
The test section’s Mach number in wind tunnel testing is a significant metric for evaluating system performance.The quality of the flow field in the wind tunnel is contingent upon the system's capacity to maintai...The test section’s Mach number in wind tunnel testing is a significant metric for evaluating system performance.The quality of the flow field in the wind tunnel is contingent upon the system's capacity to maintain stability across various working conditions.The process flow in wind tunnel testing is inherently complex,resulting in a system characterized by nonlinearity,time lag,and multiple working conditions.To implement the predictive control algorithm,a precise Mach number prediction model must be created.Therefore,this report studies the method for Mach number prediction modelling in wind tunnel flow fields with various working conditions.Firstly,this paper introduces a continuous transonic wind tunnel.The key physical quantities affecting the flow field of the wind tunnel are determined by analyzing its structure and blowing process.Secondly,considering the nonlinear and time-lag characteristics of the wind tunnel system,a CNN-LSTM model is employed to establish the Mach number prediction model by combining the 1D-CNN algorithm with the LSTM model,which has long and short-term memory functions.Then,the attention mechanism is incorporated into the CNN-LSTM prediction model to enable the model to focus more on data with greater information importance,thereby enhancing the model's training effectiveness.The application results ultimately demonstrate the efficacy of the proposed approach.展开更多
In this paper, we study the low Mach number limit of a compressible nonisothermal model for nematic liquid crystals in a bounded domain. We establish the uniform estimates with respect to the Mach number, and thus pro...In this paper, we study the low Mach number limit of a compressible nonisothermal model for nematic liquid crystals in a bounded domain. We establish the uniform estimates with respect to the Mach number, and thus prove the convergence to the solution of the incompressible model for nematic liquid crystals.展开更多
A novel extension to SMAC scheme is proposed for variable density flows under low Mach number approximation. The algorithm is based on a predictor—corrector time integration scheme that employs a projection method fo...A novel extension to SMAC scheme is proposed for variable density flows under low Mach number approximation. The algorithm is based on a predictor—corrector time integration scheme that employs a projection method for the momentum equation. A constant-coefficient Poisson equation is solved for the pressure following both the predictor and corrector steps to satisfy the continuity equation at each time step. The proposed algorithm has second order centrally differenced convective fluxes with upwinding based on Cell Peclet number while diffusive flux are viscous fourth order accurate. Spatial discretization is performed on a collocated grid system that offers computational simplicity and straight forward extension to curvilinear coordinate systems. The algorithm is kinetic energy preserving. Further in this paper robustness and accuracy are demonstrated by performing test on channel flow with non-Boussinesq condition on different temperature ratios.展开更多
Convective heating of the rocket base caused by high-temperature reverse flow has long been a focus of thermal protection research.With distinctive structural characteristics,the base thermal environment of a twin-noz...Convective heating of the rocket base caused by high-temperature reverse flow has long been a focus of thermal protection research.With distinctive structural characteristics,the base thermal environment of a twin-nozzle engine proves more susceptible to the recirculation region than its multi-nozzle counterparts.During the transonic stage,significant alterations in the flow field structure at the rocket base strongly influence the recirculation region.This study investigated the thermal environment of the rocket base with a twin-nozzle configuration in freestream at Mach numbers of 0.6 to 3.0.Results indicate that the freestream Mach number significantly affects the thermal environment at the rocket base during the transonic stage.The increase of Mach number from 0.6 to 1.0 causes the convective heating of the rocket base to increase by 7.7 times.This phenomenon arises due to the plume-induced shock wave caused by the impact of the supersonic free shear layer and plume shear layer while the flight speed exceeds the sound speed.The interaction between the shock wave and the shear layer amplifies turbulence in the recirculation region and at the inflection point,resulting in a stronger high-temperature reverse flow.In addition,the cause of low-altitude base heating was analyzed,and it was found that the mechanism is different from the high-temperature countercurrent effect caused by plume interaction.展开更多
This paper verifies the low Mach number limit of the non-isentropic compressible magnetohydrodynamic(MHD)equations with or without the magnetic diffusion in a three-dimensional bounded domain when the temperature vari...This paper verifies the low Mach number limit of the non-isentropic compressible magnetohydrodynamic(MHD)equations with or without the magnetic diffusion in a three-dimensional bounded domain when the temperature variation is large but finite.The uniform estimates of strong solutions are established in a short time interval independent of the Mach number,provided that the slip boundary condition for the velocity and the Neumann boundary condition for the temperature are imposed and the initial data is well-prepared.展开更多
This paper develops a low-diffusion robust flux splitting scheme termed TVAP to achieve the simulation of wide-ranging Mach number flows.Based on Toro-Vazquez splitting approach,the new scheme splits inviscid flux int...This paper develops a low-diffusion robust flux splitting scheme termed TVAP to achieve the simulation of wide-ranging Mach number flows.Based on Toro-Vazquez splitting approach,the new scheme splits inviscid flux into convective system and pressure system.This method introduces Mach number splitting function and numerical sound speed to evaluate advection system.Meanwhile,pressure diffusion term,pressure momentum flux,interface pressure and interface velocity are constructed to measure pressure system.Then,typical test problems are utilized to systematically assess the robustness and accuracy of the resulting scheme.Matrix stability analysis and a series of numerical cases,such as double shear-layer problem and hypersonic viscous flow over blunt cone,demonstrate that TVAP scheme achieves excellent low diffusion,shock stability,contact discontinuity and low-speed resolution,and is potentially a good candidate for wide-ranging Mach number flows.展开更多
Based on modified Leishman-Beddoes (L-B) state space model at low Mach number (lower than 0.3), the airfoil aeroelastic system is presented in this paper. The main modifications for L-B model include a new dynamic...Based on modified Leishman-Beddoes (L-B) state space model at low Mach number (lower than 0.3), the airfoil aeroelastic system is presented in this paper. The main modifications for L-B model include a new dynamic stall criterion and revisions of normal force and pitching moment coefficient. The bifurcation diagrams, the limit cycle oscillation (LCO) phase plane plots and the time domain response figures are applied to investigating the stall flutter bifurcation behavior of airfoil aeroelastic systems with symmetry or asymmetry. It is shown that the symmetric periodical oscillation happens after subcritical bifurcation caused by dynamic stall, and the asymmetric periodical oscillation, which is caused by the interaction of dynamic stall and static divergence, only happens in the airfoil aeroelastic system with asymmetry. Validations of the modified L-B model and the airfoil aeroelastic system are presented with the experimental airload data of NACA0012 and OA207 and experimental stall flutter data of NACA0012 respectively. Results demonstrate that the airfoil aeroelastic system presented in this paper is effective and accurate, which can be applied to the investigation of airfoil stall flutter at low Mach number.展开更多
We propose an accurate and robust Roe-type scheme applied to the compressible Euler system at all Mach numbers.To study the occurrence of unstable modes during the shock wave computation,a shock instability analysis o...We propose an accurate and robust Roe-type scheme applied to the compressible Euler system at all Mach numbers.To study the occurrence of unstable modes during the shock wave computation,a shock instability analysis of several Roe-type schemes is carried out.This analysis approach allows to propose a simple and effective modification to eliminate shock instability of the Roe method for hypersonic flows.A desirable feature of this modification is that it does not resort to any additional numerical dissipation on linear degenerate waves to suppress the shock instability.With an all Mach correction strategy,the modified Roe-type scheme is further extended to solve flow problems at all Mach numbers.Numerical results that are obtained for various test cases indicate that the new scheme has a good performance in terms of accuracy and robustness.展开更多
An all speed scheme for the Isentropic Euler equations is presented in thispaper. When the Mach number tends to zero, the compressible Euler equations converge to their incompressible counterpart, in which the density...An all speed scheme for the Isentropic Euler equations is presented in thispaper. When the Mach number tends to zero, the compressible Euler equations converge to their incompressible counterpart, in which the density becomes a constant. Increasing approximation errors and severe stability constraints are the main difficultyin the low Mach regime. The key idea of our all speed scheme is the special semiimplicit time discretization, in which the low Mach number stiff term is divided intotwo parts, one being treated explicitly and the other one implicitly. Moreover, the fluxof the density equation is also treated implicitly and an elliptic type equation is derivedto obtain the density. In this way, the correct limit can be captured without requesting the mesh size and time step to be smaller than the Mach number. Compared withprevious semi-implicit methods [11,13,29], firstly, nonphysical oscillations can be suppressed by choosing proper parameter, besides, only a linear elliptic equation needs tobe solved implicitly which reduces much computational cost. We develop this semiimplicit time discretization in the framework of a first order Local Lax-Friedrichs (orRusanov) scheme and numerical tests are displayed to demonstrate its performances.展开更多
基金supported by the Structures and Materials(S&M)Research Lab of Prince Sultan Universitysupport of Prince Sultan University in paying the article processing charges(APC)for this publication.
文摘When better fuel-air mixing in the combustion chamber or a reduction in base drag are required in vehicles,rockets,and aeroplanes,the base pressure control is activated.Controlling the base pressure and drag is necessary in both scenarios.In this work,semi-circular ribs with varying diameters(2,4,and 6 mm)positioned at six distinct positions(0.5D,1D,1.5D,2D,3D,and 4D)inside a square duct with a side of 15 mm are proposed as an efficient way to apply the passive control technique.In-depth research is done on optimising rib size for various rib sites.According to this study,the base pressure rises as rib height increases.Furthermore,the optimal location for the semi-circular ribs with a diameter of 2 mm is at 0.5D.The 1D location appears to be optimal for the 4 mm size as well.For the 6 mm size,however,the 4D position fills this function.
文摘The purpose of this work is to shed light on the effect of the pivot position on the surface pressure distribution over a 3D wing in different flight conditions.The study is intended to support the design and development of aerospace vehicles where stability analysis,performance optimization,and aircraft design are of primary importance.The following parameters are considered:Mach numbers(M)of 1.3,1.8,2.3,2.8,3.3,and 3.8,angle of incidence(θ)in the range from 5°to 25°,pivot position from h=0.2 to 1.The results of the CFD numerical simulations match available analytical data,thereby providing evidence for the reliability of the used approach.The findings provide valuable insights into the relationship between the surface pressure distribution,the Mach number and the angle of incidence.
基金The project supported by the Basic Research on Frontier Problems in Fluid and Aerodynamics in Chinathe National Natural Science Foundation of China (19772069)
文摘By analyzing the characteristics of low Mach number perfect gas flows, a novel Slightly Compressible Model (SCM) for low Mach number perect gas flows is derived. In view of numerical calculations, this model is proved very efficient, for it is kept within thep-v frame but does not have to satisfy the time consuming divergence-free condition in order to get the incompressible Navier-Stokes equation solution. Writing the equations in the form of conservation laws, we have derived the characteristic systems which are necessary for numerical calculations. A cell-centered finite-volume method with flux difference upwind-biased schemes is used for the equation solutions and a new Exact Newton Relaxation (ENR) implicit method is developed. Various computed results are presented to validate the present model. Laminar flow solutions over a circular cylinder with wake developing and vortex shedding are presented. Results for inviscid flow over a sphere are compared in excellent agreement with the exact analytic incompressible solution. Three-dimensional viscous flow solutions over sphere and prolate spheroid are also calculated and compared well with experiments and other incompressible solutions. Finally, good convergent performances are shown for sphere viscous flows.
文摘In the present study,the base pressure variations induced by the presence of a cavity,known to have a strong influence of the behaviour of supersonic projectiles,are investigated through numerical solution of the balance equations for mass,momentum,and energy.An area ratio of four is considered and numerical simulations are carried out at Mach M=1.2,1.4,1.6,and 1.8 assuming no cavity or cavity locations 0.5D,1D,1.5D,and 2D.The inlet pressure of the nozzle is considered as a flow variable.The Taguchi method is also used,and the considered cases are then analyzed using a full factorial experimental design.The results show that the cavity is effective in increasing the base pressure for the conditions examined.For other nozzle pressure ratios,cavities do not lead to passive control due the change in the reattachment length.The distribution of wall pressure reveals that,in general,a cavity used to implement passive control of the base pressure does not adversely influence the flow pattern in the domain.
文摘This paper is concerned with the low Mach number limit for the compressible Navier-Stokes equations in an exterior domain. We present here an approach based on Strichartz estimate defined on a non trapping exterior domain and we will be able to show the compactness and strong convergence of the velocity vector field.
文摘This work aims to compute stability derivatives in the Newtonian limit in pitch when the Mach number tends to infinity.In such conditions,these stability derivatives depend on the Ogive’s shape and not the Mach number.Generally,the Mach number independence principle becomes effective from M=10 and above.The Ogive nose is obtained through a circular arc on the cone surface.Accordingly,the following arc slopes are consideredλ=5,10,15,−5,−10,and−15.It is found that the stability derivatives decrease due to the growth inλfrom 5 to 15 and vice versa.Forλ=5 and 10,the damping derivative declines with an increase inλfrom 5 to 10.Yet,for the damping derivatives,the minimum location remains at a pivot position,h=0.75 for large values ofλ.Hence,whenλ=−15,the damping derivatives are independent of the cone angles for most pivot positions except in the early twenty percent of the leading edge.
文摘On super-sonic or trans-sonic planar cascade wind tunnel of free jet intermittent type, wind blowing experiments were performed on the typical sections of stator and rotor blades in the last stage of ultra-ultra-critical steam turbine with extra-long blade of 1200mm. The influences of attack angle and Mach number on the aerodynamic performances of these sections of the blade profiles were verified, and their operating ranges were also specified.
基金Yuxi HU was supported by the NNSFC (11701556)the Yue Qi Young Scholar ProjectChina University of Mining and Technology (Beijing)。
文摘We investigate the low Mach number limit for the isentropic compressible NavierStokes equations with a revised Maxwell's law(with Galilean invariance) in R^(3). By applying the uniform estimates of the error system, it is proven that the solutions of the isentropic Navier-Stokes equations with a revised Maxwell's law converge to that of the incompressible Navier-Stokes equations as the Mach number tends to zero. Moreover, the convergence rates are also obtained.
基金support provided by the Deutsche Forschun-gsgemeinschaft(DFG,German Research Foundation)through the project GRK 2160/1“Droplet Interaction Technologies”and through the project no.457811052
文摘Considering droplet phenomena at low Mach numbers,large differences in the magnitude of the occurring characteristic waves are presented.As acoustic phenomena often play a minor role in such applications,classical explicit schemes which resolve these waves suffer from a very restrictive timestep restriction.In this work,a novel scheme based on a specific level set ghost fluid method and an implicit-explicit(IMEX)flux splitting is proposed to overcome this timestep restriction.A fully implicit narrow band around the sharp phase interface is combined with a splitting of the convective and acoustic phenomena away from the interface.In this part of the domain,the IMEX Runge-Kutta time discretization and the high order discontinuous Galerkin spectral element method are applied to achieve high accuracies in the bulk phases.It is shown that for low Mach numbers a significant gain in computational time can be achieved compared to a fully explicit method.Applica-tions to typical droplet dynamic phenomena validate the proposed method and illustrate its capabilities.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11232011 and 11402262the 111 Project under Grant No B07033the China Postdoctoral Science Foundation Funded Project under Grant No 2014M561833
文摘The interaction of a shock wave with a spherical helium bubble is investigated numerically by using the high- resolution piecewise parabolic method (PPM), in which the viscous and turbulence effects are both considered. The bubble is of the same size and is accelerated by a planar shock of different Mach numbers (Ma). The re- suits of low Ma cases agree quantitatively with those of experiments [G. Layes, O. Le M4tayer. Phys. Fluids 19 (2007) 042105]. With the increase of Ma, the final geometry of the bubble becomes quite different, the com- pression ratio is highly raised, and the time-dependent mean bubble velocity is also influenced. The compression ratios measured can be well normalized when Ma is low, while less agreement has been achieved for high Ma cases. In addition, the mixedness between two fluids is enhanced greatly as Ma increases. Some existed scaling laws of these quantities for the shock wave strength cannot be directly applied to high Ma cases.
基金funded by the National Natural Science Foundation of China(No.61503069)the Fundamental Research Funds for the Central Universities(N150404020).
文摘Mach number is a key metric in the evaluation of wind tunnel flow field performance.This complex process of wind tunnel test mainly has the problems of nonlinearity and time lag.In order to overcome the problems and control the Mach number stability,this paper proposes a new method of Mach number prediction based on a nonlinear autoregressive exogenous-genetic algorithm-Elman(NARX-GA-Elman)model,which adopts NARX as the basic framework,determines the order of the input variables by using the false nearest neighbor(FNN),and uses the dynamic nonlinear network Elman to fit the model,and finally uses the global optimization algorithm GA to optimize the weight thresholds in the model to establish the Mach number prediction model with optimal performance under single working condition.By comparing with the traditional algorithm,the prediction accuracy of the model is improved by 61.5%,and the control accuracy is improved by 55.7%,which demonstrates that the model has very high prediction accuracy and good stability performance.
基金funded by the National Natural Science Foundation of China(No.61503069)the Fundamental Research Funds for the Central Universities(N150404020).
文摘The test section’s Mach number in wind tunnel testing is a significant metric for evaluating system performance.The quality of the flow field in the wind tunnel is contingent upon the system's capacity to maintain stability across various working conditions.The process flow in wind tunnel testing is inherently complex,resulting in a system characterized by nonlinearity,time lag,and multiple working conditions.To implement the predictive control algorithm,a precise Mach number prediction model must be created.Therefore,this report studies the method for Mach number prediction modelling in wind tunnel flow fields with various working conditions.Firstly,this paper introduces a continuous transonic wind tunnel.The key physical quantities affecting the flow field of the wind tunnel are determined by analyzing its structure and blowing process.Secondly,considering the nonlinear and time-lag characteristics of the wind tunnel system,a CNN-LSTM model is employed to establish the Mach number prediction model by combining the 1D-CNN algorithm with the LSTM model,which has long and short-term memory functions.Then,the attention mechanism is incorporated into the CNN-LSTM prediction model to enable the model to focus more on data with greater information importance,thereby enhancing the model's training effectiveness.The application results ultimately demonstrate the efficacy of the proposed approach.
基金supported by NSFC(11171154)supported in part by by NSFC(11671193)A Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions
文摘In this paper, we study the low Mach number limit of a compressible nonisothermal model for nematic liquid crystals in a bounded domain. We establish the uniform estimates with respect to the Mach number, and thus prove the convergence to the solution of the incompressible model for nematic liquid crystals.
文摘A novel extension to SMAC scheme is proposed for variable density flows under low Mach number approximation. The algorithm is based on a predictor—corrector time integration scheme that employs a projection method for the momentum equation. A constant-coefficient Poisson equation is solved for the pressure following both the predictor and corrector steps to satisfy the continuity equation at each time step. The proposed algorithm has second order centrally differenced convective fluxes with upwinding based on Cell Peclet number while diffusive flux are viscous fourth order accurate. Spatial discretization is performed on a collocated grid system that offers computational simplicity and straight forward extension to curvilinear coordinate systems. The algorithm is kinetic energy preserving. Further in this paper robustness and accuracy are demonstrated by performing test on channel flow with non-Boussinesq condition on different temperature ratios.
基金supported by the Outstanding Research Project of Shen Yuan Honors CollegeAcademic Excellence Foundation of BUAA for PhD StudentsNational Key Laboratory of aerospace liquid propulsion
文摘Convective heating of the rocket base caused by high-temperature reverse flow has long been a focus of thermal protection research.With distinctive structural characteristics,the base thermal environment of a twin-nozzle engine proves more susceptible to the recirculation region than its multi-nozzle counterparts.During the transonic stage,significant alterations in the flow field structure at the rocket base strongly influence the recirculation region.This study investigated the thermal environment of the rocket base with a twin-nozzle configuration in freestream at Mach numbers of 0.6 to 3.0.Results indicate that the freestream Mach number significantly affects the thermal environment at the rocket base during the transonic stage.The increase of Mach number from 0.6 to 1.0 causes the convective heating of the rocket base to increase by 7.7 times.This phenomenon arises due to the plume-induced shock wave caused by the impact of the supersonic free shear layer and plume shear layer while the flight speed exceeds the sound speed.The interaction between the shock wave and the shear layer amplifies turbulence in the recirculation region and at the inflection point,resulting in a stronger high-temperature reverse flow.In addition,the cause of low-altitude base heating was analyzed,and it was found that the mechanism is different from the high-temperature countercurrent effect caused by plume interaction.
基金supported by National Natural Science Foundation of China(Grant Nos.11971477,12131007 and 11761141008)the Fundamental Research Funds for the Central Universitiesthe Research Funds of Renmin University of China(Grant No.18XNLG30)。
文摘This paper verifies the low Mach number limit of the non-isentropic compressible magnetohydrodynamic(MHD)equations with or without the magnetic diffusion in a three-dimensional bounded domain when the temperature variation is large but finite.The uniform estimates of strong solutions are established in a short time interval independent of the Mach number,provided that the slip boundary condition for the velocity and the Neumann boundary condition for the temperature are imposed and the initial data is well-prepared.
基金supported by the Space Science and Technology Fund Project of China(No.2020-HT-XG-14)。
文摘This paper develops a low-diffusion robust flux splitting scheme termed TVAP to achieve the simulation of wide-ranging Mach number flows.Based on Toro-Vazquez splitting approach,the new scheme splits inviscid flux into convective system and pressure system.This method introduces Mach number splitting function and numerical sound speed to evaluate advection system.Meanwhile,pressure diffusion term,pressure momentum flux,interface pressure and interface velocity are constructed to measure pressure system.Then,typical test problems are utilized to systematically assess the robustness and accuracy of the resulting scheme.Matrix stability analysis and a series of numerical cases,such as double shear-layer problem and hypersonic viscous flow over blunt cone,demonstrate that TVAP scheme achieves excellent low diffusion,shock stability,contact discontinuity and low-speed resolution,and is potentially a good candidate for wide-ranging Mach number flows.
基金Aeronautical Science Foundation of China (08A52003)Science and Technology Foundation of Rotorcraft Aeromechanics Laboratory (9140C4001010901)
文摘Based on modified Leishman-Beddoes (L-B) state space model at low Mach number (lower than 0.3), the airfoil aeroelastic system is presented in this paper. The main modifications for L-B model include a new dynamic stall criterion and revisions of normal force and pitching moment coefficient. The bifurcation diagrams, the limit cycle oscillation (LCO) phase plane plots and the time domain response figures are applied to investigating the stall flutter bifurcation behavior of airfoil aeroelastic systems with symmetry or asymmetry. It is shown that the symmetric periodical oscillation happens after subcritical bifurcation caused by dynamic stall, and the asymmetric periodical oscillation, which is caused by the interaction of dynamic stall and static divergence, only happens in the airfoil aeroelastic system with asymmetry. Validations of the modified L-B model and the airfoil aeroelastic system are presented with the experimental airload data of NACA0012 and OA207 and experimental stall flutter data of NACA0012 respectively. Results demonstrate that the airfoil aeroelastic system presented in this paper is effective and accurate, which can be applied to the investigation of airfoil stall flutter at low Mach number.
基金This work was supported by the National Natural Science Foundation of China(No.11472004)the Foundation of Innovation of NUDT(No.B150106).
文摘We propose an accurate and robust Roe-type scheme applied to the compressible Euler system at all Mach numbers.To study the occurrence of unstable modes during the shock wave computation,a shock instability analysis of several Roe-type schemes is carried out.This analysis approach allows to propose a simple and effective modification to eliminate shock instability of the Roe method for hypersonic flows.A desirable feature of this modification is that it does not resort to any additional numerical dissipation on linear degenerate waves to suppress the shock instability.With an all Mach correction strategy,the modified Roe-type scheme is further extended to solve flow problems at all Mach numbers.Numerical results that are obtained for various test cases indicate that the new scheme has a good performance in terms of accuracy and robustness.
基金supported by the French"Commissariat´a l’Energie Atomique(CEA)"(Centre de Saclay)in the frame of the contract"ASTRE",#SAV 34160the Marie Curie Actions of the European Commission in the frame of the DEASE project(MESTCT-2005-021122)。
文摘An all speed scheme for the Isentropic Euler equations is presented in thispaper. When the Mach number tends to zero, the compressible Euler equations converge to their incompressible counterpart, in which the density becomes a constant. Increasing approximation errors and severe stability constraints are the main difficultyin the low Mach regime. The key idea of our all speed scheme is the special semiimplicit time discretization, in which the low Mach number stiff term is divided intotwo parts, one being treated explicitly and the other one implicitly. Moreover, the fluxof the density equation is also treated implicitly and an elliptic type equation is derivedto obtain the density. In this way, the correct limit can be captured without requesting the mesh size and time step to be smaller than the Mach number. Compared withprevious semi-implicit methods [11,13,29], firstly, nonphysical oscillations can be suppressed by choosing proper parameter, besides, only a linear elliptic equation needs tobe solved implicitly which reduces much computational cost. We develop this semiimplicit time discretization in the framework of a first order Local Lax-Friedrichs (orRusanov) scheme and numerical tests are displayed to demonstrate its performances.