483 research outputs found
Entropy stable DGSEM for nonlinear hyperbolic systems in nonconservative form with application to two-phase flows
In this work, we consider the discretization of nonlinear hyperbolic systems
in nonconservative form with the high-order discontinuous Galerkin spectral
element method (DGSEM) based on collocation of quadrature and interpolation
points (Kopriva and Gassner, J. Sci. Comput., 44 (2010), pp.136--155; Carpenter
et al., SIAM J. Sci. Comput., 36 (2014), pp.~B835-B867). We present a general
framework for the design of such schemes that satisfy a semi-discrete entropy
inequality for a given convex entropy function at any approximation order. The
framework is closely related to the one introduced for conservation laws by
Chen and Shu (J. Comput. Phys., 345 (2017), pp.~427--461) and relies on the
modification of the integral over discretization elements where we replace the
physical fluxes by entropy conservative numerical fluxes from Castro et al.
(SIAM J. Numer. Anal., 51 (2013), pp.~1371--1391), while entropy stable
numerical fluxes are used at element interfaces. Time discretization is
performed with strong-stability preserving Runge-Kutta schemes. We use this
framework for the discretization of two systems in one space-dimension: a
system with a nonconservative product associated to a
linearly-degenerate field for which the DGSEM fails to capture the physically
relevant solution, and the isentropic Baer-Nunziato model. For the latter, we
derive conditions on the numerical parameters of the discrete scheme to further
keep positivity of the partial densities and a maximum principle on the void
fractions. Numerical experiments support the conclusions of the present
analysis and highlight stability and robustness of the present schemes
Coupling techniques for nonlinear hyperbolic equations. III. The well-balanced approximation of thick interfaces
We continue our analysis of the coupling between nonlinear hyperbolic
problems across possibly resonant interfaces. In the first two parts of this
series, we introduced a new framework for coupling problems which is based on
the so-called thin interface model and uses an augmented formulation and an
additional unknown for the interface location; this framework has the advantage
of avoiding any explicit modeling of the interface structure. In the present
paper, we pursue our investigation of the augmented formulation and we
introduce a new coupling framework which is now based on the so-called thick
interface model. For scalar nonlinear hyperbolic equations in one space
variable, we observe that the Cauchy problem is well-posed. Then, our main
achievement in the present paper is the design of a new well-balanced finite
volume scheme which is adapted to the thick interface model, together with a
proof of its convergence toward the unique entropy solution (for a broad class
of nonlinear hyperbolic equations). Due to the presence of a possibly resonant
interface, the standard technique based on a total variation estimate does not
apply, and DiPerna's uniqueness theorem must be used. Following a method
proposed by Coquel and LeFloch, our proof relies on discrete entropy
inequalities for the coupling problem and an estimate of the discrete entropy
dissipation in the proposed scheme.Comment: 21 page
Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations
The application of a new implicit unconditionally stable high resolution total variation diminishing (TVD) scheme to steady state calculations. It is a member of a one parameter family of explicit and implicit second order accurate schemes developed by Harten for the computation of weak solutions of hyperbolic conservation laws. This scheme is guaranteed not to generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments show that this scheme not only has a rapid convergence rate, but also generates a highly resolved approximation to the steady state solution. A detailed implementation of the implicit scheme for the one and two dimensional compressible inviscid equations of gas dynamics is presented. Some numerical computations of one and two dimensional fluid flows containing shocks demonstrate the efficiency and accuracy of this new scheme
A radiation-hydrodynamics scheme valid from the transport to the diffusion limit
We present in this paper the numerical treatment of the coupling between
hydrodynamics and radiative transfer. The fluid is modeled by classical
conservation laws (mass, momentum and energy) and the radiation by the grey
moment system. The scheme introduced is able to compute accurate
numerical solution over a broad class of regimes from the transport to the
diffusive limits. We propose an asymptotic preserving modification of the HLLE
scheme in order to treat correctly the diffusion limit. Several numerical
results are presented, which show that this approach is robust and have the
correct behavior in both the diffusive and free-streaming limits. In the last
numerical example we test this approach on a complex physical case by
considering the collapse of a gas cloud leading to a proto-stellar structure
which, among other features, exhibits very steep opacity gradients.Comment: 29 pages, submitted to Journal of Computational physic
- …