20,798 research outputs found
Adaptive multiresolution computations applied to detonations
A space-time adaptive method is presented for the reactive Euler equations
describing chemically reacting gas flow where a two species model is used for
the chemistry. The governing equations are discretized with a finite volume
method and dynamic space adaptivity is introduced using multiresolution
analysis. A time splitting method of Strang is applied to be able to consider
stiff problems while keeping the method explicit. For time adaptivity an
improved Runge--Kutta--Fehlberg scheme is used. Applications deal with
detonation problems in one and two space dimensions. A comparison of the
adaptive scheme with reference computations on a regular grid allow to assess
the accuracy and the computational efficiency, in terms of CPU time and memory
requirements.Comment: Zeitschrift f\"ur Physicalische Chemie, accepte
To Split or Not to Split, That Is the Question in Some Shallow Water Equations
In this paper we analyze the use of time splitting techniques for solving
shallow water equation. We discuss some properties that these schemes should
satisfy so that interactions between the source term and the shock waves are
controlled. This paper shows that these schemes must be well balanced in the
meaning expressed by Greenberg and Leroux [5]. More specifically, we analyze in
what cases it is enough to verify an Approximate C-property and in which cases
it is required to verify an Exact C-property (see [1], [2]). We also include
some numerical tests in order to justify our reasoning
An unconditionally stable algorithm for generalized thermoelasticity based on operator-splitting and time-discontinuous Galerkin finite element methods
An efficient time-stepping algorithm is proposed based on operator-splitting and the space–time discontinuous Galerkin finite element method for problems in the non-classical theory of thermoelasticity. The non-classical theory incorporates three models: the classical theory based on Fourier’s law of heat conduction resulting in a hyperbolic–parabolic coupled system, a non-classical theory of a fully-hyperbolic extension, and a combination of the two. The general problem is split into two contractive sub-problems, namely the mechanical phase and the thermal phase. Each sub-problem is discretized using the space–time discontinuous Galerkin finite element method. The sub-problems are stable which then leads to unconditional stability of the global product algorithm. A number of numerical examples are presented to demonstrate the performance and capability of the method
An Unstaggered Constrained Transport Method for the 3D Ideal Magnetohydrodynamic Equations
Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations
in more than one space dimension must either confront the challenge of
controlling errors in the discrete divergence of the magnetic field, or else be
faced with nonlinear numerical instabilities. One approach for controlling the
discrete divergence is through a so-called constrained transport method, which
is based on first predicting a magnetic field through a standard finite volume
solver, and then correcting this field through the appropriate use of a
magnetic vector potential. In this work we develop a constrained transport
method for the 3D ideal MHD equations that is based on a high-resolution wave
propagation scheme. Our proposed scheme is the 3D extension of the 2D scheme
developed by Rossmanith [SIAM J. Sci. Comp. 28, 1766 (2006)], and is based on
the high-resolution wave propagation method of Langseth and LeVeque [J. Comp.
Phys. 165, 126 (2000)]. In particular, in our extension we take great care to
maintain the three most important properties of the 2D scheme: (1) all
quantities, including all components of the magnetic field and magnetic
potential, are treated as cell-centered; (2) we develop a high-resolution wave
propagation scheme for evolving the magnetic potential; and (3) we develop a
wave limiting approach that is applied during the vector potential evolution,
which controls unphysical oscillations in the magnetic field. One of the key
numerical difficulties that is novel to 3D is that the transport equation that
must be solved for the magnetic vector potential is only weakly hyperbolic. In
presenting our numerical algorithm we describe how to numerically handle this
problem of weak hyperbolicity, as well as how to choose an appropriate gauge
condition. The resulting scheme is applied to several numerical test cases.Comment: 46 pages, 12 figure
- …