For compressible two-phase displacement problem,the modified upwind finite difference fractionalsteps schemes are put forward.Some techniques,such as calculus of variations,commutative law of multiplicationof differen...For compressible two-phase displacement problem,the modified upwind finite difference fractionalsteps schemes are put forward.Some techniques,such as calculus of variations,commutative law of multiplicationof difference operators,decomposition of high order difference operators,the theory of prior estimates and tech-niques are used.Optimal order estimates in L^2 norm are derived for the error in the approximate solution.Thismethod has already been applied to the numerical simulation of seawater intrusion and migration-accumulationof oil resources.展开更多
In this paper,a new sharp-interface approach to simulate compressible multiphase flows is proposed.The new scheme consists of a high order WENO finite volume scheme for solving the Euler equations coupled with a high ...In this paper,a new sharp-interface approach to simulate compressible multiphase flows is proposed.The new scheme consists of a high order WENO finite volume scheme for solving the Euler equations coupled with a high order pathconservative discontinuous Galerkin finite element scheme to evolve an indicator function that tracks the material interface.At the interface our method applies ghost cells to compute the numerical flux,as the ghost fluid method.However,unlike the original ghost fluid scheme of Fedkiw et al.[15],the state of the ghost fluid is derived from an approximate-state Riemann solver,similar to the approach proposed in[25],but based on a much simpler formulation.Our formulation leads only to one single scalar nonlinear algebraic equation that has to be solved at the interface,instead of the system used in[25].Away from the interface,we use the new general Osher-type flux recently proposed by Dumbser and Toro[13],which is a simple but complete Riemann solver,applicable to general hyperbolic conservation laws.The time integration is performed using a fully-discrete one-step scheme,based on the approaches recently proposed in[5,7].This allows us to evolve the system also with time-accurate local time stepping.Due to the sub-cell resolution and the subsequent more restrictive time-step constraint of the DG scheme,a local evolution for the indicator function is applied,which is matched with the finite volume scheme for the solution of the Euler equations that runs with a larger time step.The use of a locally optimal time step avoids the introduction of excessive numerical diffusion in the finite volume scheme.Two different fluids have been used,namely an ideal gas and a weakly compressible fluid modeled by the Tait equation.Several tests have been computed to assess the accuracy and the performance of the new high order scheme.A verification of our algorithm has been carefully carried out using exact solutions as well as a comparison with other numerical reference solutions.The material interface is resolved sharply and accurately without spurious oscillations in the pressure field.展开更多
基金Supported by the Major State Basic Research Program of China (Grant No.1999032803)the National Natural Science Foundation of China (Grant No.10372052,10271066)the Decorate Foundation of the Ministry Education of China (Grant No.20030422047)
文摘For compressible two-phase displacement problem,the modified upwind finite difference fractionalsteps schemes are put forward.Some techniques,such as calculus of variations,commutative law of multiplicationof difference operators,decomposition of high order difference operators,the theory of prior estimates and tech-niques are used.Optimal order estimates in L^2 norm are derived for the error in the approximate solution.Thismethod has already been applied to the numerical simulation of seawater intrusion and migration-accumulationof oil resources.
文摘In this paper,a new sharp-interface approach to simulate compressible multiphase flows is proposed.The new scheme consists of a high order WENO finite volume scheme for solving the Euler equations coupled with a high order pathconservative discontinuous Galerkin finite element scheme to evolve an indicator function that tracks the material interface.At the interface our method applies ghost cells to compute the numerical flux,as the ghost fluid method.However,unlike the original ghost fluid scheme of Fedkiw et al.[15],the state of the ghost fluid is derived from an approximate-state Riemann solver,similar to the approach proposed in[25],but based on a much simpler formulation.Our formulation leads only to one single scalar nonlinear algebraic equation that has to be solved at the interface,instead of the system used in[25].Away from the interface,we use the new general Osher-type flux recently proposed by Dumbser and Toro[13],which is a simple but complete Riemann solver,applicable to general hyperbolic conservation laws.The time integration is performed using a fully-discrete one-step scheme,based on the approaches recently proposed in[5,7].This allows us to evolve the system also with time-accurate local time stepping.Due to the sub-cell resolution and the subsequent more restrictive time-step constraint of the DG scheme,a local evolution for the indicator function is applied,which is matched with the finite volume scheme for the solution of the Euler equations that runs with a larger time step.The use of a locally optimal time step avoids the introduction of excessive numerical diffusion in the finite volume scheme.Two different fluids have been used,namely an ideal gas and a weakly compressible fluid modeled by the Tait equation.Several tests have been computed to assess the accuracy and the performance of the new high order scheme.A verification of our algorithm has been carefully carried out using exact solutions as well as a comparison with other numerical reference solutions.The material interface is resolved sharply and accurately without spurious oscillations in the pressure field.
