424 research outputs found
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
Nonlinear Evolution Equations: Analysis and Numerics
The qualitative theory of nonlinear evolution equations is an
important tool for studying the dynamical behavior of systems in
science and technology. A thorough understanding of the complex
behavior of such systems requires detailed analytical and numerical
investigations of the underlying partial differential equations
Characteristic Evolution and Matching
I review the development of numerical evolution codes for general relativity
based upon the characteristic initial value problem. Progress in characteristic
evolution is traced from the early stage of 1D feasibility studies to 2D
axisymmetric codes that accurately simulate the oscillations and gravitational
collapse of relativistic stars and to current 3D codes that provide pieces of a
binary black hole spacetime. Cauchy codes have now been successful at
simulating all aspects of the binary black hole problem inside an artificially
constructed outer boundary. A prime application of characteristic evolution is
to extend such simulations to null infinity where the waveform from the binary
inspiral and merger can be unambiguously computed. This has now been
accomplished by Cauchy-characteristic extraction, where data for the
characteristic evolution is supplied by Cauchy data on an extraction worldtube
inside the artificial outer boundary. The ultimate application of
characteristic evolution is to eliminate the role of this outer boundary by
constructing a global solution via Cauchy-characteristic matching. Progress in
this direction is discussed.Comment: New version to appear in Living Reviews 2012. arXiv admin note:
updated version of arXiv:gr-qc/050809
Hyperbolic Techniques in Modelling, Analysis and Numerics
Several research areas are flourishing on the roots of the breakthroughs in conservation laws that took place in the last two decades. The meeting played a key role in providing contacts among the different branches that are currently developing. All the invitees shared the same common background that consists of the analytical and numerical techniques for nonlinear hyperbolic balance laws. However, their fields of applications and their levels of abstraction are very diverse. The workshop was the unique opportunity to share ideas about analytical issues like the fine-structure of singular solutions or the validity of entropy solution concepts. It turned out that generalized hyperbolic techniques are able to handle the challenges posed by new applications. The design of efficient structure preserving methods turned out to be the major line of development in numerical analysis
Recommended from our members
The numerical solution of stefan problems with front-tracking and smoothing methods
- …