Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-materia...Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.展开更多
In this paper, we are concerned with the local structural stability of one-dimensional shock waves in radiation hydrodynamics described by the isentropic Euler-Boltzmann equations. Even though in this radiation hydrod...In this paper, we are concerned with the local structural stability of one-dimensional shock waves in radiation hydrodynamics described by the isentropic Euler-Boltzmann equations. Even though in this radiation hydrodynamics model, the radiative effects can be understood as source terms to the isentropic Euler equations of hydrodynamics, in general the radiation field has singularities propagated in an angular domain issuing from the initial point across which the density is discontinuous. This is the major difficulty in the stability analysis of shocks. Under certain assumptions on the radiation parameters, we show there exists a local weak solution to the initial value problem of the one dimensional Euler-Boltzmann equations, in which the radiation intensity is continuous, while the density and velocity are piecewise Lipschitz continuous with a strong discontinuity representing the shock-front. The existence of such a solution indicates that shock waves are structurally stable, at least local in time, in radiation hydrodynamics.展开更多
基金Supported by the National Natural Science Foundation of China(10476011)~~
文摘Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.
基金supported by NNSF of China(10971134,11031001,91230102,11371250)
文摘In this paper, we are concerned with the local structural stability of one-dimensional shock waves in radiation hydrodynamics described by the isentropic Euler-Boltzmann equations. Even though in this radiation hydrodynamics model, the radiative effects can be understood as source terms to the isentropic Euler equations of hydrodynamics, in general the radiation field has singularities propagated in an angular domain issuing from the initial point across which the density is discontinuous. This is the major difficulty in the stability analysis of shocks. Under certain assumptions on the radiation parameters, we show there exists a local weak solution to the initial value problem of the one dimensional Euler-Boltzmann equations, in which the radiation intensity is continuous, while the density and velocity are piecewise Lipschitz continuous with a strong discontinuity representing the shock-front. The existence of such a solution indicates that shock waves are structurally stable, at least local in time, in radiation hydrodynamics.
