The computation of compressible flows becomesmore challengingwhen the Mach number has different orders of magnitude.When the Mach number is of order one,modern shock capturing methods are able to capture shocks and ot...The computation of compressible flows becomesmore challengingwhen the Mach number has different orders of magnitude.When the Mach number is of order one,modern shock capturing methods are able to capture shocks and other complex structures with high numerical resolutions.However,if theMach number is small,the acoustic waves lead to stiffness in time and excessively large numerical viscosity,thus demanding much smaller time step and mesh size than normally needed for incompressible flow simulation.In this paper,we develop an all-speed asymptotic preserving(AP)numerical scheme for the compressible isentropic Euler and Navier-Stokes equations that is uniformly stable and accurate for all Mach numbers.Our idea is to split the system into two parts:one involves a slow,nonlinear and conservative hyperbolic system adequate for the use of modern shock capturing methods and the other a linear hyperbolic system which contains the stiff acoustic dynamics,to be solved implicitly.This implicit part is reformulated into a standard pressure Poisson projection system and thus possesses sufficient structure for efficient fast Fourier transform solution techniques.In the zero Mach number limit,the scheme automatically becomes a projection method-like incompressible solver.We present numerical results in one and two dimensions in both compressible and incompressible regimes.展开更多
The quasi-neutral limit of the multi-dimensional non-isentropic bipolar Euler-Poisson system is considered in the present paper. It is shown that for well-prepared initial data the smooth solution of the non-isentropi...The quasi-neutral limit of the multi-dimensional non-isentropic bipolar Euler-Poisson system is considered in the present paper. It is shown that for well-prepared initial data the smooth solution of the non-isentropic bipolar Euler-Poisson system converges strongly to the compressible non-isentropic Euler equations as the Debye length goes to zero.展开更多
This paper is concerned with the dissipation of solutions of the isentropic Navier-Stokes equations in even and bigger than two multi-dimensions. Pointwise estimates of the time-asymptotic shape of the solutions are o...This paper is concerned with the dissipation of solutions of the isentropic Navier-Stokes equations in even and bigger than two multi-dimensions. Pointwise estimates of the time-asymptotic shape of the solutions are obtained and the generalized Huygan's principle is exhibited. The approch of the paper is based on the detailed analysis of the Green function of Iinearized system. This is used to study the coupling of nonlinear diffesion waves.展开更多
The possibility of pressure control with the structural change of a safety valve is investigated Safety valve is commonly used as safety devices for numerous applications which include boilers,ships,industrial plant...The possibility of pressure control with the structural change of a safety valve is investigated Safety valve is commonly used as safety devices for numerous applications which include boilers,ships,industrial plants,and piping Setting and stopping pressures of a safety valve, p set and p sto ,are traditionally adjusted with a fine tuning of seat ring and valve ring heights, h sr and h vr However, it is not easy to achieve the proper setting and stopping pressures of a safety valve in practice The depth of inside and outside grooves in a valve, d i and d o are modified and their effects on setting and stopping pressures of a safety vlave are tested The most appropriate values appear 1 0 mm in d i and 0 5～1 0 mm in d o,respectively The valve ring height, h vr ,shows that the best results can be achieved at 2 3 mm for setting pressures of 0 1～0 4 MPa and 1 0 mm for setting pressures of 0 5～1 0 MPa The stopping pressures increases with the increase of seat ring height, h sr , upto certain h sr value and then becomes independent to the seat ring height This implies that there exists the optimum h sr ,which provides the largest flow rate and the proper stopping pressure Stopping pressures of a safety valve are adjusted with the seat ring and valve ring heights This study,however,demonstrated that the modification of value grooves also changes setting and stopping pressures of a safety valve Therefore,the proper selection in dimensions of the inside and outside grooves should be considered for the safety valve design展开更多
The complete form of the vertical vorticity tendency equation （the complete-form vorticity equation） is derived from the Ertel potential vorticity equation to contain thermodynamic factors. In this study, a new comp...The complete form of the vertical vorticity tendency equation （the complete-form vorticity equation） is derived from the Ertel potential vorticity equation to contain thermodynamic factors. In this study, a new complete-form vorticity equation, which has the same form as the original complete-form vorticity equation, is deduced from the absolute vorticity vector equation combined with the continuity equation and the expression of three-dimensional （3D） entropy gradient. By comparing the complete-form vorticity equation with the classical vertical vorticity equation, it is found that regardless of whether or not the isentropic surface is tilting, the two vorticity equations are in essence the same. The ＂baroclinic term＂ of the complete-form vorticity equation is exactly equal to the solenoidal term of the classical one, and there is a significant amount of cancellation between the two baroclinic items （the ＂slantwise term＂ and the horizontal vorticity change term） in the complete-form vorticity equation. In operational weather analysis, the tilt of the isentropic surface can be diagnosed according to the density of the isotherm on the upper-level isobaric map. For synoptic-scale motion, the vertical vorticity produced by the tilt of the isentropic surface is due to the contribution of atmospheric baroclinicity, which is measured by the solenoid. The 3D solenoid is parallel to the isentropic surface, so the more tilted the isentropic surface, the bigger the projection of the 3D solenoid in the vertical direction. The baroclinic contribution can be interpreted based on the PV thinking theory, but the relationship between the vorticity field and the potential vorticity field is not immediate.展开更多
基金J.-G.Liu was supported by NSF grant DMS 10-11738.J.Haack and S.Jin were supported by NSF grant DMS-0608720the NSF FRG grant”Collaborative research on Kinetic Description of Multiscale Phenomena:Modeling,Theory and Computation”(NSF DMS-0757285).
文摘The computation of compressible flows becomesmore challengingwhen the Mach number has different orders of magnitude.When the Mach number is of order one,modern shock capturing methods are able to capture shocks and other complex structures with high numerical resolutions.However,if theMach number is small,the acoustic waves lead to stiffness in time and excessively large numerical viscosity,thus demanding much smaller time step and mesh size than normally needed for incompressible flow simulation.In this paper,we develop an all-speed asymptotic preserving(AP)numerical scheme for the compressible isentropic Euler and Navier-Stokes equations that is uniformly stable and accurate for all Mach numbers.Our idea is to split the system into two parts:one involves a slow,nonlinear and conservative hyperbolic system adequate for the use of modern shock capturing methods and the other a linear hyperbolic system which contains the stiff acoustic dynamics,to be solved implicitly.This implicit part is reformulated into a standard pressure Poisson projection system and thus possesses sufficient structure for efficient fast Fourier transform solution techniques.In the zero Mach number limit,the scheme automatically becomes a projection method-like incompressible solver.We present numerical results in one and two dimensions in both compressible and incompressible regimes.
基金supported by National Natural Science Foundation of China(Grant No.40890154)National Basic Research Program(Grant No.2005CB321700)+5 种基金supported by National Natural Science Foundation of China(Grant No.10701011)supported by National Natural Science Foundation of China(Grant No.10431060)Beijing Nova Program,Program for New Century Excellent Talentsin University,Huo Ying Dong Foundation(Grant No.111033)supported by National Natural Science Foundation of China(Grant No.10901011)Beijing Municipal Natural Science Foundation(Grant No.1102009)Foundation for Talents of Beijing(Grant No.20081D0501500171)
文摘The quasi-neutral limit of the multi-dimensional non-isentropic bipolar Euler-Poisson system is considered in the present paper. It is shown that for well-prepared initial data the smooth solution of the non-isentropic bipolar Euler-Poisson system converges strongly to the compressible non-isentropic Euler equations as the Debye length goes to zero.
基金Supported in part by National Natural Science Foundationof China (19871065) Hua-Cheng Grant
文摘This paper is concerned with the dissipation of solutions of the isentropic Navier-Stokes equations in even and bigger than two multi-dimensions. Pointwise estimates of the time-asymptotic shape of the solutions are obtained and the generalized Huygan's principle is exhibited. The approch of the paper is based on the detailed analysis of the Green function of Iinearized system. This is used to study the coupling of nonlinear diffesion waves.
文摘The possibility of pressure control with the structural change of a safety valve is investigated Safety valve is commonly used as safety devices for numerous applications which include boilers,ships,industrial plants,and piping Setting and stopping pressures of a safety valve, p set and p sto ,are traditionally adjusted with a fine tuning of seat ring and valve ring heights, h sr and h vr However, it is not easy to achieve the proper setting and stopping pressures of a safety valve in practice The depth of inside and outside grooves in a valve, d i and d o are modified and their effects on setting and stopping pressures of a safety vlave are tested The most appropriate values appear 1 0 mm in d i and 0 5～1 0 mm in d o,respectively The valve ring height, h vr ,shows that the best results can be achieved at 2 3 mm for setting pressures of 0 1～0 4 MPa and 1 0 mm for setting pressures of 0 5～1 0 MPa The stopping pressures increases with the increase of seat ring height, h sr , upto certain h sr value and then becomes independent to the seat ring height This implies that there exists the optimum h sr ,which provides the largest flow rate and the proper stopping pressure Stopping pressures of a safety valve are adjusted with the seat ring and valve ring heights This study,however,demonstrated that the modification of value grooves also changes setting and stopping pressures of a safety valve Therefore,the proper selection in dimensions of the inside and outside grooves should be considered for the safety valve design
基金Supported by the National Natural Science Foundation of China(41475042 and 41175043)China Meteorological Administration Special Public Welfare Research Fund(GYHY201406002)
文摘The complete form of the vertical vorticity tendency equation （the complete-form vorticity equation） is derived from the Ertel potential vorticity equation to contain thermodynamic factors. In this study, a new complete-form vorticity equation, which has the same form as the original complete-form vorticity equation, is deduced from the absolute vorticity vector equation combined with the continuity equation and the expression of three-dimensional （3D） entropy gradient. By comparing the complete-form vorticity equation with the classical vertical vorticity equation, it is found that regardless of whether or not the isentropic surface is tilting, the two vorticity equations are in essence the same. The ＂baroclinic term＂ of the complete-form vorticity equation is exactly equal to the solenoidal term of the classical one, and there is a significant amount of cancellation between the two baroclinic items （the ＂slantwise term＂ and the horizontal vorticity change term） in the complete-form vorticity equation. In operational weather analysis, the tilt of the isentropic surface can be diagnosed according to the density of the isotherm on the upper-level isobaric map. For synoptic-scale motion, the vertical vorticity produced by the tilt of the isentropic surface is due to the contribution of atmospheric baroclinicity, which is measured by the solenoid. The 3D solenoid is parallel to the isentropic surface, so the more tilted the isentropic surface, the bigger the projection of the 3D solenoid in the vertical direction. The baroclinic contribution can be interpreted based on the PV thinking theory, but the relationship between the vorticity field and the potential vorticity field is not immediate.
