In this paper,we develop a fourth-order conservative wavelet-based shock-capturing scheme.The scheme is constructed by combining a wavelet collocation upwind method with the monotonic tangent of hyperbola for interfac...In this paper,we develop a fourth-order conservative wavelet-based shock-capturing scheme.The scheme is constructed by combining a wavelet collocation upwind method with the monotonic tangent of hyperbola for interface capturing(THINC)technique.We employ boundary variation diminishing(BVD)reconstruction to enhance the scheme’s effectiveness in handling shocks.First,we prove that wavelet collocation upwind schemes based on interpolating wavelets can be reformulated into a conservative form within the framework of wavelet theory,forming the foundation of the proposed scheme.The new fourthorder accurate scheme possesses significantly better spectral resolution than the fifth-and even seventh-order WENO-Z(weighted essentially non-oscillatory)schemes over the entire wave-number range.Moreover,the inherent low-pass filtering property of the wavelet bases allows them to filter high-frequency numerical oscillations,endowing the wavelet upwind scheme with robustness and accuracy in solving problems under extreme conditions.Notably,due to the wavelet multiresolution approximation,the proposed scheme possesses a distinctive shape-preserving property absent in the WENO-Z schemes and the fifth-order schemes with BVD reconstruction based on polynomials.Furthermore,compared to the fifth-order scheme with BVD reconstruction based on polynomials—which is significantly superior to the WENO schemes—the proposed scheme further enhances the ability to capture discontinuities.展开更多
Multiphase flows widely exist in various scientific and engineering fields,and strongly compressible multiphase flows commonly occur in practical applications,which makes them an important part of computational fluid ...Multiphase flows widely exist in various scientific and engineering fields,and strongly compressible multiphase flows commonly occur in practical applications,which makes them an important part of computational fluid dynamics.In this study,an axisymmetric adaptive multiresolution smooth particle hydrodynamics(SPH)model is proposed to solve various strongly compressible multiphase flow problems.In the present model,the governing equations are discretized in cylindrical polar coordinates,and an improved volume adaptive scheme is developed to better solve the problem of excessive volume change in strongly compressible multiphase flows.On this basis,combined with the adaptive particle refinement technique,an adaptive multiresolution scheme is proposed in this study.In addition,the high-order differential operator and diffusion correction term are utilized to improve the accuracy and stability.The effectiveness of the model is verified by testing four typical strongly compressible multiphase flow problems.By comparing the results of adaptive multiresolution SPH with other numerical results or experimental data,we can conclude that the present SPH method effectively models strongly compressible multiphase flows.展开更多
This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-laye...This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.展开更多
In this paper,we investigate the linear stability/instability of the planar Couette flow to the two-dimensional compressible Euler-Euler system for(ρ,u)and(n,v)with the sound speeds c_(1)and c_(2)respectively,coupled...In this paper,we investigate the linear stability/instability of the planar Couette flow to the two-dimensional compressible Euler-Euler system for(ρ,u)and(n,v)with the sound speeds c_(1)and c_(2)respectively,coupled each other through the drag force on T×R.It is shown in general for the different sound speeds c_(1)≠c_(2)that the perturbations of the densities(ρ,n)and the velocities(u,v)demonstrate the stability in any fixed finite time interval(0,T],besides,excluding the zero mode,the densities(ρ,n)and the irrotational components of the velocities(u,v)obey the algebraic time-growth rates,while the rotational components of the velocities(u,v)exhibit the algebraic time-decay rates due to the inviscid damping.For the case c_(1)=c_(2)(same sound speeds),it is proved that the relative velocity u-v decays faster than those of the velocities u,v caused by the inviscid damping,but the linear instability of the densities(ρ,n)and the irrotational components of the velocities(u,v)is also shown for some large time in spite of the inviscid damping.展开更多
Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have on...Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.展开更多
We study an initial boundary value problem for the Navier-Stokes equations of compressible viscous heat-conductive fluids in a 2-D periodic domain or the unit square domain. We establish a blow-up criterion for the lo...We study an initial boundary value problem for the Navier-Stokes equations of compressible viscous heat-conductive fluids in a 2-D periodic domain or the unit square domain. We establish a blow-up criterion for the local strong solutions in terms of the gradient of the velocity only, which coincides with the famous Beale-Kato-Majda criterion for ideal incompressible flows.展开更多
In this work, incompressible and compressible flows of background gas are characterized in argon inductively coupled plasma by using a fluid model, and the respective influence of the two flows on the plasma propertie...In this work, incompressible and compressible flows of background gas are characterized in argon inductively coupled plasma by using a fluid model, and the respective influence of the two flows on the plasma properties is specified. In the incompressible flow, only the velocity variable is calculated, while in the compressible flow, both the velocity and density variables are calculated. The compressible flow is more realistic; nevertheless, a comparison of the two types of flow is convenient for people to investigate the respective role of velocity and density variables. The peripheral symmetric profile of metastable density near the chamber sidewall is broken in the incompressible flow. At the compressible flow, the electron density increases and the electron temperature decreases. Meanwhile, the metastable density peak shifts to the dielectric window from the discharge center, besides for the peripheral density profile distortion, similar to the incompressible flow.The velocity profile at incompressible flow is not altered when changing the inlet velocity, whereas clear peak shift of velocity profile from the inlet to the outlet at compressible flow is observed as increasing the gas flow rate. The shift of velocity peak is more obvious at low pressures for it is easy to compress the rarefied gas. The velocity profile variations at compressible flow show people the concrete residing processes of background molecule and plasma species in the chamber at different flow rates. Of more significance is it implied that in the usual linear method that people use to calculate the residence time, one important parameter in the gas flow dynamics, needs to be rectified. The spatial profile of pressure simulated exhibits obvious spatial gradient. This is helpful for experimentalists to understand their gas pressure measurements that are always taken at the chamber outlet. At the end, the work specification and limitations are listed.展开更多
In this paper,the finite difference weighted essentially non-oscillatory (WENO) scheme is incorporated into the recently developed four kinds of lattice Boltzmann flux solver (LBFS) to simulate compressible flows,incl...In this paper,the finite difference weighted essentially non-oscillatory (WENO) scheme is incorporated into the recently developed four kinds of lattice Boltzmann flux solver (LBFS) to simulate compressible flows,including inviscid LBFS Ⅰ,viscous LBFS Ⅱ,hybrid LBFS Ⅲ and hybrid LBFS Ⅳ.Hybrid LBFS can automatically realize the switch between inviscid LBFS Ⅰ and viscous LBFS Ⅱ through introducing a switch function.The resultant hybrid WENO-LBFS scheme absorbs the advantages of WENO scheme and hybrid LBFS.We investigate the performance of WENO scheme based on four kinds of LBFS systematically.Numerical results indicate that the devopled hybrid WENO-LBFS scheme has high accuracy,high resolution and no oscillations.It can not only accurately calculate smooth solutions,but also can effectively capture contact discontinuities and strong shock waves.展开更多
This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiplerelaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, where...This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiplerelaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, whereas the modified single-relaxation-time （SRT） lattice Boltzmann scheme is applied for the evolution of potential energy distribution functions. The governing equations are discretized with the third-order Monotone Upwind Schemes for scalar conservation laws finite volume scheme. The choice of relaxation coefficients is discussed simply. Through the numerical simulations, it is found that compressible flows with strong shocks can be well simulated by present model. The numerical results agree well with the reference results and are better than that of the SRT version.展开更多
This work attempts to extend the fundamental theory for classic gas dynamics to viscous compressible flow,of which aeroacoustics will naturally be a special branch.As a continuation of Part I.Unbounded fluid(Mao et al...This work attempts to extend the fundamental theory for classic gas dynamics to viscous compressible flow,of which aeroacoustics will naturally be a special branch.As a continuation of Part I.Unbounded fluid(Mao et al.,2022),this paper studies the source of longitudinal field at solid boundary,caused by the on-wall kinematic and viscous dynamic coupling of longitudinal and transverse processes.We find that at this situation the easiest choice for the two independent thermodynamic variables is the dimensionless pressure P and temperature T.The two-level structure of boundary dynamics of longitudinal field is obtained by applying the continuity equation and its normal derivative to the surface.We show that the boundary dilatation flux represents faithfully the boundary production of vortex sound and entropy sound,and the mutual generation mechanism of the longitudinal and transverse fields on the boundary does not occur symmetrically"at the samc level,but appears along a zigzag route.At the first level,it is the pressure gradient that generates vorticity unidirectionally;while at the second level,it is the vorticity that generates dilatation unidirectionally.展开更多
The presence of a cavity changes the mean and fluctuating pressure distributions. Similarities are observed between a cylindrical cavity and a rectangular cavity for a compressible flow.The type of cavity flow field d...The presence of a cavity changes the mean and fluctuating pressure distributions. Similarities are observed between a cylindrical cavity and a rectangular cavity for a compressible flow.The type of cavity flow field depends on the diameter-to-depth ratio and the length-to-depth ratio.The feedback loop is responsible for the generation of discrete acoustic tones. In this study, the selfsustained oscillation for a compressible cylindrical cavity flow was investigated experimentally. For open-type cavities, the power spectra show that the strength of resonance depends on the diameterto-depth ratio（4.43–43.0） and the incoming boundary layer thickness-to-depth ratio（0.72–7.0）. The effective streamwise length is used as the characteristic length to estimate the Strouhal number. At higher modes, there is a large deviation from Rossiter＇s formula for rectangular cavities. The gradient-based searching method was used to evaluate the values of the empirical parameters. Less phase lag and a lower convection velocity are observed.展开更多
A modified SIMPLEC method which can solve compressible flows at low Mach number is introduced and used to study thermoacoustic waves induced by a rapid change of temperature at a solid wall and alternating- direction ...A modified SIMPLEC method which can solve compressible flows at low Mach number is introduced and used to study thermoacoustic waves induced by a rapid change of temperature at a solid wall and alternating- direction flows generated by thermoacoustic effects in a ta- pered resonator. The results indicate that the algorithm adopted in this paper can be used for calculating com- pressible flows and thermoacoustic waves. It is found that the pressure and velocity in the resonator behave as stand- ing waves, and the tapered resonator can suppress high- frequency harmonic waves as observed in a cylindrical res- onator.展开更多
The paper presents the implementation of non-Newtonian fluid properties for compressible multiphase solver in the open source framework OpenFOAM. The transport models for Power Law, Cross Power Law, Casson, Bird-Carre...The paper presents the implementation of non-Newtonian fluid properties for compressible multiphase solver in the open source framework OpenFOAM. The transport models for Power Law, Cross Power Law, Casson, Bird-Carreau and Herschel-Bulkley fluids were included in the thermophysical model library. Appropriate non-Newtonian liquids have been chosen from literature, and pressure driven test simulations are carried out. Therefore, the solver compressibleInterFoam is used to compute air-liquid mixture flows over a backward facing step. A validation of the novel models has been performed by means of a sample-based comparison of the strain rate viscosity relation. The theoretical rheological properties of the selected liquids agree well with the results of the simulated data.展开更多
In this paper,several exact expressions for the mean heat flux at the wall(qw)for the compressible turbulent channel flows are derived by using the internal energy equation or the total energy equation.Two different r...In this paper,several exact expressions for the mean heat flux at the wall(qw)for the compressible turbulent channel flows are derived by using the internal energy equation or the total energy equation.Two different routes,including the FIK method and the RD method,can be applied.The direct numerical simulation data of compressible channel flows at different Reynolds and Mach numbers verify the correctness of the derived formulas.Discussions related to the different energy equations,and different routes are carried out,and we may arrive at the conclusion that most of the formulas derived in the present work are just mathematical ones and that they generally are lacking in clear physical interpretation in our opinion.They can be used to estimate qw,but might not be suitable for exploring the underlying physics.展开更多
A highly efficient three-dimensional （31）） Lattice Boltzmann （LB） model for high-speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara [Phys. Rev. E 69 （200...A highly efficient three-dimensional （31）） Lattice Boltzmann （LB） model for high-speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara [Phys. Rev. E 69 （2004） 056702]. The convection term is discretized by the Non-oscillatory, containing No free parameters and Dissipative （NND） scheme, which effectively damps oscillations at discontinuities. To be more consistent with the kinetic theory of viscosity and to further improve the numerical stability, an additional dissipation term is introduced. Model parameters are chosen in such a way that the von Neumann stability criterion is satisfied. The new model is validated by well-known benchmarks, （i） Riemann problems, including the problem with Lax shock tube and a newly designed shock tube problem with high Mach number; （ii） reaction of shock wave on droplet or bubble. Good agreements are obtained between LB results and exact ones or previously reported solutions. The model is capable of simulating flows from subsonic to supersonic and capturing jumps resulted from shock waves.展开更多
This paper presents a deep reflection on the advective wave equations for velocity vector and dilatation discovered in the past decade.We show that these equations can form the theoretical basis of modern gas dynamics...This paper presents a deep reflection on the advective wave equations for velocity vector and dilatation discovered in the past decade.We show that these equations can form the theoretical basis of modern gas dynamics,because they dominate not only various complex viscous and heat-conducting gas flows but also their associated longitudinal waves,including aero-generated sound.Current aeroacoustics theory has been developing in a manner quite independently of gas dynamics;it is based on the advective wave equations for thermodynamic variables,say the exact Phillips equation of relative disturbance pressure as a representative one.However,these equations do not cover the fluid flow that generates and propagates sound waves.In using them,one has to assume simplified base-flow models,which we argue is the main theoretical obstacle to identifying sound source and achieving effective noise control.Instead,we show that the Phillips equation and alike is nothing but the first integral of the dilatation equation that also governs the longitudinal part of the flow field.Therefore,we conclude that modern aeroacoustics should merge back into the general unsteady gas dynamics as a special branch of it,with dilatation of multiple sources being a new additional and sharper sound variable.展开更多
The mechanism of shocklets is studied theoretically and numerically for the stationary fluid, uniform compressible flow, and boundary layer flow. The conditions that trigger shock waves for sound wave, weak discontinu...The mechanism of shocklets is studied theoretically and numerically for the stationary fluid, uniform compressible flow, and boundary layer flow. The conditions that trigger shock waves for sound wave, weak discontinuity, and Tollmien-Schlichting （T-S） wave in compressible flows are investigated. The relations between the three types of waves and shocklets are further analyzed and discussed. Different stages of the shocklet formation process are simulated. The results show that the three waves in compressible flows will transfer to shocklets only when the initial disturbance amplitudes are greater than the certain threshold values. In compressible boundary layers, the shocklets evolved from T-S wave exist only in a finite region near the surface instead of the whole wavefront.展开更多
The short-range property of interactions between scales in the compressible turbulent flow was examined. An estimation of the short-range scale scope and some formulae for the short-range eddy stress and heat transfe...The short-range property of interactions between scales in the compressible turbulent flow was examined. An estimation of the short-range scale scope and some formulae for the short-range eddy stress and heat transfer etc. were given. A concept of resonant-range interactions between extremely contiguous scales was introduced and some formulae for the resonant-range eddy stress and heat transfer etc. were also given. Multi-scale equations for the compressible turbulent flows were presented. The multi-scale equations are approximately closed and do not contain any empirical constants. The compressibility effects on turbulence are determined by the Farve averaged variables and the nonlinear relationships between the Farve- and physical-averaged variables.展开更多
We further develop the lattice Boltzmann （LB） model [Physica A 382 （2007） 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook （BGK） collision term in the LB equation, which...We further develop the lattice Boltzmann （LB） model [Physica A 382 （2007） 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook （BGK） collision term in the LB equation, which makes the model suitable for simulating flows with different Prandtl numbers. Secondly, the flux limiter finite difference （FLFD） scheme is employed to calculate the convection term of the LB equation, which makes the unphysical oscillations at discontinuities be effectively suppressed and the numerical dissipations be significantly diminished. The proposed model is validated by recovering results of some well-known benchmarks, including （i） The thermal Couette flow; （ii） One- and two-dlmenslonal FLiemann problems. Good agreements are obtained between LB results and the exact ones or previously reported solutions. The flexibility, together with the high accuracy of the new model, endows the proposed model considerable potential for tracking some long-standing problems and for investigating nonlinear nonequilibrium complex systems.展开更多
In this paper, a new flux limiter scheme with the splitting technique is successfully incorporated into a multiple-relaxation-time lattice Boltzmann （LB） model for shacked compressible flows. The proposed flux limit...In this paper, a new flux limiter scheme with the splitting technique is successfully incorporated into a multiple-relaxation-time lattice Boltzmann （LB） model for shacked compressible flows. The proposed flux limiter scheme is efficient in decreasing the artificial oscillations and numerical diffusion around the interface. Due to the kinetic nature, some interface problems being difficult to handle at the macroscopic level can be modeled more naturally through the LB method. Numerical simulations for the Richtmyer-Meshkov instability show that with the new model the computed interfaces are smoother and more consistent with physical analysis. The growth rates of bubble and spike present a satisfying agreement with the theoretical predictions and other numerical simulations.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.11925204).
文摘In this paper,we develop a fourth-order conservative wavelet-based shock-capturing scheme.The scheme is constructed by combining a wavelet collocation upwind method with the monotonic tangent of hyperbola for interface capturing(THINC)technique.We employ boundary variation diminishing(BVD)reconstruction to enhance the scheme’s effectiveness in handling shocks.First,we prove that wavelet collocation upwind schemes based on interpolating wavelets can be reformulated into a conservative form within the framework of wavelet theory,forming the foundation of the proposed scheme.The new fourthorder accurate scheme possesses significantly better spectral resolution than the fifth-and even seventh-order WENO-Z(weighted essentially non-oscillatory)schemes over the entire wave-number range.Moreover,the inherent low-pass filtering property of the wavelet bases allows them to filter high-frequency numerical oscillations,endowing the wavelet upwind scheme with robustness and accuracy in solving problems under extreme conditions.Notably,due to the wavelet multiresolution approximation,the proposed scheme possesses a distinctive shape-preserving property absent in the WENO-Z schemes and the fifth-order schemes with BVD reconstruction based on polynomials.Furthermore,compared to the fifth-order scheme with BVD reconstruction based on polynomials—which is significantly superior to the WENO schemes—the proposed scheme further enhances the ability to capture discontinuities.
基金Supported by the Guangxi Natural Science Foundation(Grant No.2021GXNSFBA196008)the Guangxi Science and Technology Development Program(Grant No.GuikeAD22035189).
文摘Multiphase flows widely exist in various scientific and engineering fields,and strongly compressible multiphase flows commonly occur in practical applications,which makes them an important part of computational fluid dynamics.In this study,an axisymmetric adaptive multiresolution smooth particle hydrodynamics(SPH)model is proposed to solve various strongly compressible multiphase flow problems.In the present model,the governing equations are discretized in cylindrical polar coordinates,and an improved volume adaptive scheme is developed to better solve the problem of excessive volume change in strongly compressible multiphase flows.On this basis,combined with the adaptive particle refinement technique,an adaptive multiresolution scheme is proposed in this study.In addition,the high-order differential operator and diffusion correction term are utilized to improve the accuracy and stability.The effectiveness of the model is verified by testing four typical strongly compressible multiphase flow problems.By comparing the results of adaptive multiresolution SPH with other numerical results or experimental data,we can conclude that the present SPH method effectively models strongly compressible multiphase flows.
基金supported by the National Basic Research Program of China(2009CB724100)the National Natural Science Foundation of China(10632050,10872205,11072248).
文摘This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.
基金supported by the National Natural Science Foundation of China(11931010,12326613,12331007)the Beijing Scholar Foundation of Beijing Municipal Committeethe key research project of Academy for Multidisciplinary Studies,Capital Normal University。
文摘In this paper,we investigate the linear stability/instability of the planar Couette flow to the two-dimensional compressible Euler-Euler system for(ρ,u)and(n,v)with the sound speeds c_(1)and c_(2)respectively,coupled each other through the drag force on T×R.It is shown in general for the different sound speeds c_(1)≠c_(2)that the perturbations of the densities(ρ,n)and the velocities(u,v)demonstrate the stability in any fixed finite time interval(0,T],besides,excluding the zero mode,the densities(ρ,n)and the irrotational components of the velocities(u,v)obey the algebraic time-growth rates,while the rotational components of the velocities(u,v)exhibit the algebraic time-decay rates due to the inviscid damping.For the case c_(1)=c_(2)(same sound speeds),it is proved that the relative velocity u-v decays faster than those of the velocities u,v caused by the inviscid damping,but the linear instability of the densities(ρ,n)and the irrotational components of the velocities(u,v)is also shown for some large time in spite of the inviscid damping.
基金supported by National Science Foundation of China (10576015)
文摘Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.
基金supported by the China Postdoctoral Science Foundation (20090450333)supported by the National Basic Research Program (2005CB321700)NSFC (40890154)
文摘We study an initial boundary value problem for the Navier-Stokes equations of compressible viscous heat-conductive fluids in a 2-D periodic domain or the unit square domain. We establish a blow-up criterion for the local strong solutions in terms of the gradient of the velocity only, which coincides with the famous Beale-Kato-Majda criterion for ideal incompressible flows.
基金Project supported by the National Natural Science Foundations of China(Grant No.11305023)
文摘In this work, incompressible and compressible flows of background gas are characterized in argon inductively coupled plasma by using a fluid model, and the respective influence of the two flows on the plasma properties is specified. In the incompressible flow, only the velocity variable is calculated, while in the compressible flow, both the velocity and density variables are calculated. The compressible flow is more realistic; nevertheless, a comparison of the two types of flow is convenient for people to investigate the respective role of velocity and density variables. The peripheral symmetric profile of metastable density near the chamber sidewall is broken in the incompressible flow. At the compressible flow, the electron density increases and the electron temperature decreases. Meanwhile, the metastable density peak shifts to the dielectric window from the discharge center, besides for the peripheral density profile distortion, similar to the incompressible flow.The velocity profile at incompressible flow is not altered when changing the inlet velocity, whereas clear peak shift of velocity profile from the inlet to the outlet at compressible flow is observed as increasing the gas flow rate. The shift of velocity peak is more obvious at low pressures for it is easy to compress the rarefied gas. The velocity profile variations at compressible flow show people the concrete residing processes of background molecule and plasma species in the chamber at different flow rates. Of more significance is it implied that in the usual linear method that people use to calculate the residence time, one important parameter in the gas flow dynamics, needs to be rectified. The spatial profile of pressure simulated exhibits obvious spatial gradient. This is helpful for experimentalists to understand their gas pressure measurements that are always taken at the chamber outlet. At the end, the work specification and limitations are listed.
基金This study was supported by the National Natural Science Foundation of China(Grants 11372168,11772179).
文摘In this paper,the finite difference weighted essentially non-oscillatory (WENO) scheme is incorporated into the recently developed four kinds of lattice Boltzmann flux solver (LBFS) to simulate compressible flows,including inviscid LBFS Ⅰ,viscous LBFS Ⅱ,hybrid LBFS Ⅲ and hybrid LBFS Ⅳ.Hybrid LBFS can automatically realize the switch between inviscid LBFS Ⅰ and viscous LBFS Ⅱ through introducing a switch function.The resultant hybrid WENO-LBFS scheme absorbs the advantages of WENO scheme and hybrid LBFS.We investigate the performance of WENO scheme based on four kinds of LBFS systematically.Numerical results indicate that the devopled hybrid WENO-LBFS scheme has high accuracy,high resolution and no oscillations.It can not only accurately calculate smooth solutions,but also can effectively capture contact discontinuities and strong shock waves.
基金supported by the Innovation Fund for Aerospace Science and Technology of China(Grant No.2009200066)the Aeronautical Science Fund of China(Grant No.20111453012)
文摘This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiplerelaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, whereas the modified single-relaxation-time （SRT） lattice Boltzmann scheme is applied for the evolution of potential energy distribution functions. The governing equations are discretized with the third-order Monotone Upwind Schemes for scalar conservation laws finite volume scheme. The choice of relaxation coefficients is discussed simply. Through the numerical simulations, it is found that compressible flows with strong shocks can be well simulated by present model. The numerical results agree well with the reference results and are better than that of the SRT version.
基金supported by the National Natural Science Foundation of China(Grant Nos.12102365,91752202,11472016,11621202,and 12272371).
文摘This work attempts to extend the fundamental theory for classic gas dynamics to viscous compressible flow,of which aeroacoustics will naturally be a special branch.As a continuation of Part I.Unbounded fluid(Mao et al.,2022),this paper studies the source of longitudinal field at solid boundary,caused by the on-wall kinematic and viscous dynamic coupling of longitudinal and transverse processes.We find that at this situation the easiest choice for the two independent thermodynamic variables is the dimensionless pressure P and temperature T.The two-level structure of boundary dynamics of longitudinal field is obtained by applying the continuity equation and its normal derivative to the surface.We show that the boundary dilatation flux represents faithfully the boundary production of vortex sound and entropy sound,and the mutual generation mechanism of the longitudinal and transverse fields on the boundary does not occur symmetrically"at the samc level,but appears along a zigzag route.At the first level,it is the pressure gradient that generates vorticity unidirectionally;while at the second level,it is the vorticity that generates dilatation unidirectionally.
基金support of the Ministry of Science and Technology (No. MOST 103-2923-E-006MY3)
文摘The presence of a cavity changes the mean and fluctuating pressure distributions. Similarities are observed between a cylindrical cavity and a rectangular cavity for a compressible flow.The type of cavity flow field depends on the diameter-to-depth ratio and the length-to-depth ratio.The feedback loop is responsible for the generation of discrete acoustic tones. In this study, the selfsustained oscillation for a compressible cylindrical cavity flow was investigated experimentally. For open-type cavities, the power spectra show that the strength of resonance depends on the diameterto-depth ratio（4.43–43.0） and the incoming boundary layer thickness-to-depth ratio（0.72–7.0）. The effective streamwise length is used as the characteristic length to estimate the Strouhal number. At higher modes, there is a large deviation from Rossiter＇s formula for rectangular cavities. The gradient-based searching method was used to evaluate the values of the empirical parameters. Less phase lag and a lower convection velocity are observed.
基金supported by the National Natural Science Foundation of China (50890182,10972226)
文摘A modified SIMPLEC method which can solve compressible flows at low Mach number is introduced and used to study thermoacoustic waves induced by a rapid change of temperature at a solid wall and alternating- direction flows generated by thermoacoustic effects in a ta- pered resonator. The results indicate that the algorithm adopted in this paper can be used for calculating com- pressible flows and thermoacoustic waves. It is found that the pressure and velocity in the resonator behave as stand- ing waves, and the tapered resonator can suppress high- frequency harmonic waves as observed in a cylindrical res- onator.
文摘The paper presents the implementation of non-Newtonian fluid properties for compressible multiphase solver in the open source framework OpenFOAM. The transport models for Power Law, Cross Power Law, Casson, Bird-Carreau and Herschel-Bulkley fluids were included in the thermophysical model library. Appropriate non-Newtonian liquids have been chosen from literature, and pressure driven test simulations are carried out. Therefore, the solver compressibleInterFoam is used to compute air-liquid mixture flows over a backward facing step. A validation of the novel models has been performed by means of a sample-based comparison of the strain rate viscosity relation. The theoretical rheological properties of the selected liquids agree well with the results of the simulated data.
基金the National Natural Science Foundation of China(Grant Nos.11822208,11772297,and 91852205)Guangdong provincial Key Laboratory(Grant No.2019B20203001).
文摘In this paper,several exact expressions for the mean heat flux at the wall(qw)for the compressible turbulent channel flows are derived by using the internal energy equation or the total energy equation.Two different routes,including the FIK method and the RD method,can be applied.The direct numerical simulation data of compressible channel flows at different Reynolds and Mach numbers verify the correctness of the derived formulas.Discussions related to the different energy equations,and different routes are carried out,and we may arrive at the conclusion that most of the formulas derived in the present work are just mathematical ones and that they generally are lacking in clear physical interpretation in our opinion.They can be used to estimate qw,but might not be suitable for exploring the underlying physics.
基金Supported by the Science Foundations of Laboratory of Computational PhysicalScience Foundation of China Academy of Engineering Physics under Grant Nos. 2009A0102005, 2009B0101012National Natural Science Foundation under Grant Nos. 10775018, 11074300, and 1107521 of China
文摘A highly efficient three-dimensional （31）） Lattice Boltzmann （LB） model for high-speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara [Phys. Rev. E 69 （2004） 056702]. The convection term is discretized by the Non-oscillatory, containing No free parameters and Dissipative （NND） scheme, which effectively damps oscillations at discontinuities. To be more consistent with the kinetic theory of viscosity and to further improve the numerical stability, an additional dissipation term is introduced. Model parameters are chosen in such a way that the von Neumann stability criterion is satisfied. The new model is validated by well-known benchmarks, （i） Riemann problems, including the problem with Lax shock tube and a newly designed shock tube problem with high Mach number; （ii） reaction of shock wave on droplet or bubble. Good agreements are obtained between LB results and exact ones or previously reported solutions. The model is capable of simulating flows from subsonic to supersonic and capturing jumps resulted from shock waves.
基金supported by the National Natural Science Foundation of China(Grant Nos.12102365,91752202 and 11472016)Luoqin Liu was supported by the Hundred Talents Program of the Chinese Academy of Sciences(CAS).
文摘This paper presents a deep reflection on the advective wave equations for velocity vector and dilatation discovered in the past decade.We show that these equations can form the theoretical basis of modern gas dynamics,because they dominate not only various complex viscous and heat-conducting gas flows but also their associated longitudinal waves,including aero-generated sound.Current aeroacoustics theory has been developing in a manner quite independently of gas dynamics;it is based on the advective wave equations for thermodynamic variables,say the exact Phillips equation of relative disturbance pressure as a representative one.However,these equations do not cover the fluid flow that generates and propagates sound waves.In using them,one has to assume simplified base-flow models,which we argue is the main theoretical obstacle to identifying sound source and achieving effective noise control.Instead,we show that the Phillips equation and alike is nothing but the first integral of the dilatation equation that also governs the longitudinal part of the flow field.Therefore,we conclude that modern aeroacoustics should merge back into the general unsteady gas dynamics as a special branch of it,with dilatation of multiple sources being a new additional and sharper sound variable.
基金supported by the National Natural Science Foundation of China(No.10872018)
文摘The mechanism of shocklets is studied theoretically and numerically for the stationary fluid, uniform compressible flow, and boundary layer flow. The conditions that trigger shock waves for sound wave, weak discontinuity, and Tollmien-Schlichting （T-S） wave in compressible flows are investigated. The relations between the three types of waves and shocklets are further analyzed and discussed. Different stages of the shocklet formation process are simulated. The results show that the three waves in compressible flows will transfer to shocklets only when the initial disturbance amplitudes are greater than the certain threshold values. In compressible boundary layers, the shocklets evolved from T-S wave exist only in a finite region near the surface instead of the whole wavefront.
文摘The short-range property of interactions between scales in the compressible turbulent flow was examined. An estimation of the short-range scale scope and some formulae for the short-range eddy stress and heat transfer etc. were given. A concept of resonant-range interactions between extremely contiguous scales was introduced and some formulae for the resonant-range eddy stress and heat transfer etc. were also given. Multi-scale equations for the compressible turbulent flows were presented. The multi-scale equations are approximately closed and do not contain any empirical constants. The compressibility effects on turbulence are determined by the Farve averaged variables and the nonlinear relationships between the Farve- and physical-averaged variables.
基金Supported by the Science Foundations of LCP and CAEP under Grant Nos. 2009A0102005 and 2009B0101012National Natural Science Foundation of China under Grant Nos. 11075021, 11074300, and 11074303+3 种基金National Basic Research Program (973 Program) under Grant No. 2007CB815105Fundamental Research Funds for the Central University under Grant No. 2010YS03Technology Support Program of LangFang under Grant Nos. 2010011029/30/31Science Foundation of NCIAE under Grant No. 2008-ky-13
文摘We further develop the lattice Boltzmann （LB） model [Physica A 382 （2007） 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook （BGK） collision term in the LB equation, which makes the model suitable for simulating flows with different Prandtl numbers. Secondly, the flux limiter finite difference （FLFD） scheme is employed to calculate the convection term of the LB equation, which makes the unphysical oscillations at discontinuities be effectively suppressed and the numerical dissipations be significantly diminished. The proposed model is validated by recovering results of some well-known benchmarks, including （i） The thermal Couette flow; （ii） One- and two-dlmenslonal FLiemann problems. Good agreements are obtained between LB results and the exact ones or previously reported solutions. The flexibility, together with the high accuracy of the new model, endows the proposed model considerable potential for tracking some long-standing problems and for investigating nonlinear nonequilibrium complex systems.
基金Supported by the Science Foundation of Laboratory of Computational Physics, Science Foundation of China Academy of Engineering Physics under Grant Nos. 2009A0102005, 2009B0101012National Basic Research Program of China under Grant No. 2007CB815105+1 种基金National Natural Science Foundation of China under Grant Nos. 11074300, 11075021, and 11074303the Fundamental Research Funds for the Central Universities under Grant No. 2010YS03
文摘In this paper, a new flux limiter scheme with the splitting technique is successfully incorporated into a multiple-relaxation-time lattice Boltzmann （LB） model for shacked compressible flows. The proposed flux limiter scheme is efficient in decreasing the artificial oscillations and numerical diffusion around the interface. Due to the kinetic nature, some interface problems being difficult to handle at the macroscopic level can be modeled more naturally through the LB method. Numerical simulations for the Richtmyer-Meshkov instability show that with the new model the computed interfaces are smoother and more consistent with physical analysis. The growth rates of bubble and spike present a satisfying agreement with the theoretical predictions and other numerical simulations.
