8,454 research outputs found
Entropy Stable Finite Volume Approximations for Ideal Magnetohydrodynamics
This article serves as a summary outlining the mathematical entropy analysis
of the ideal magnetohydrodynamic (MHD) equations. We select the ideal MHD
equations as they are particularly useful for mathematically modeling a wide
variety of magnetized fluids. In order to be self-contained we first motivate
the physical properties of a magnetic fluid and how it should behave under the
laws of thermodynamics. Next, we introduce a mathematical model built from
hyperbolic partial differential equations (PDEs) that translate physical laws
into mathematical equations. After an overview of the continuous analysis, we
thoroughly describe the derivation of a numerical approximation of the ideal
MHD system that remains consistent to the continuous thermodynamic principles.
The derivation of the method and the theorems contained within serve as the
bulk of the review article. We demonstrate that the derived numerical
approximation retains the correct entropic properties of the continuous model
and show its applicability to a variety of standard numerical test cases for
MHD schemes. We close with our conclusions and a brief discussion on future
work in the area of entropy consistent numerical methods and the modeling of
plasmas
Numerical simulation of conservation laws with moving grid nodes: Application to tsunami wave modelling
In the present article we describe a few simple and efficient finite volume
type schemes on moving grids in one spatial dimension combined with appropriate
predictor-corrector method to achieve higher resolution. The underlying finite
volume scheme is conservative and it is accurate up to the second order in
space. The main novelty consists in the motion of the grid. This new dynamic
aspect can be used to resolve better the areas with large solution gradients or
any other special features. No interpolation procedure is employed, thus
unnecessary solution smearing is avoided, and therefore, our method enjoys
excellent conservation properties. The resulting grid is completely
redistributed according the choice of the so-called monitor function. Several
more or less universal choices of the monitor function are provided. Finally,
the performance of the proposed algorithm is illustrated on several examples
stemming from the simple linear advection to the simulation of complex shallow
water waves. The exact well-balanced property is proven. We believe that the
techniques described in our paper can be beneficially used to model tsunami
wave propagation and run-up.Comment: 46 pages, 7 figures, 7 tables, 94 references. Accepted to
Geosciences. Other author's papers can be downloaded at
http://www.denys-dutykh.com
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
Capturing nonclassical shocks in nonlinear elastodynamic with a conservative finite volume scheme
For a model of nonlinear elastodynamics, we construct a finite volume scheme
which is able to capture nonclassical shocks (also called undercompressive
shocks). Those shocks verify an entropy inequality but are not admissible in
the sense of Liu. They verify a kinetic relation which describes the jump, and
keeps an information on the equilibrium between a vanishing dispersion and a
vanishing diffusion. The scheme pre-sented here is by construction exact when
the initial data is an isolated nonclassical shock. In general, it does not
introduce any diffusion near shocks, and hence nonclas-sical solutions are
correctly approximated. The method is fully conservative and does not use any
shock-tracking mesh. This approach is tested and validated on several test
cases. In particular, as the nonclassical shocks are not diffused at all, it is
possible to obtain large time asymptotics
- …