439 research outputs found
Variational Integrators for Reduced Magnetohydrodynamics
Reduced magnetohydrodynamics is a simplified set of magnetohydrodynamics
equations with applications to both fusion and astrophysical plasmas,
possessing a noncanonical Hamiltonian structure and consequently a number of
conserved functionals. We propose a new discretisation strategy for these
equations based on a discrete variational principle applied to a formal
Lagrangian. The resulting integrator preserves important quantities like the
total energy, magnetic helicity and cross helicity exactly (up to machine
precision). As the integrator is free of numerical resistivity, spurious
reconnection along current sheets is absent in the ideal case. If effects of
electron inertia are added, reconnection of magnetic field lines is allowed,
although the resulting model still possesses a noncanonical Hamiltonian
structure. After reviewing the conservation laws of the model equations, the
adopted variational principle with the related conservation laws are described
both at the continuous and discrete level. We verify the favourable properties
of the variational integrator in particular with respect to the preservation of
the invariants of the models under consideration and compare with results from
the literature and those of a pseudo-spectral code.Comment: 35 page
Stable Unitary Integrators for the Numerical Implementation of Continuous Unitary Transformations
The technique of continuous unitary transformations has recently been used to
provide physical insight into a diverse array of quantum mechanical systems.
However, the question of how to best numerically implement the flow equations
has received little attention. The most immediately apparent approach, using
standard Runge-Kutta numerical integration algorithms, suffers from both severe
inefficiency due to stiffness and the loss of unitarity. After reviewing the
formalism of continuous unitary transformations and Wegner's original choice
for the infinitesimal generator of the flow, we present a number of approaches
to resolving these issues including a choice of generator which induces what we
call the "uniform tangent decay flow" and three numerical integrators
specifically designed to perform continuous unitary transformations efficiently
while preserving the unitarity of flow. We conclude by applying one of the flow
algorithms to a simple calculation that visually demonstrates the many-body
localization transition.Comment: 13 pages, 4 figures, Comments welcom
Splitting and composition methods in the numerical integration of differential equations
We provide a comprehensive survey of splitting and composition methods for
the numerical integration of ordinary differential equations (ODEs). Splitting
methods constitute an appropriate choice when the vector field associated with
the ODE can be decomposed into several pieces and each of them is integrable.
This class of integrators are explicit, simple to implement and preserve
structural properties of the system. In consequence, they are specially useful
in geometric numerical integration. In addition, the numerical solution
obtained by splitting schemes can be seen as the exact solution to a perturbed
system of ODEs possessing the same geometric properties as the original system.
This backward error interpretation has direct implications for the qualitative
behavior of the numerical solution as well as for the error propagation along
time. Closely connected with splitting integrators are composition methods. We
analyze the order conditions required by a method to achieve a given order and
summarize the different families of schemes one can find in the literature.
Finally, we illustrate the main features of splitting and composition methods
on several numerical examples arising from applications.Comment: Review paper; 56 pages, 6 figures, 8 table
Adaptive Geometric Numerical Integration of Mechanical Systems
This thesis is about structure preserving numerical integration of initial value problems, i.e., so called geometric numerical integrators. In particular, we are interested in how time-step adaptivity can be achieved in conjunction with structure preserving properties without destroying the good long time integration properties which are typical for geometric integration methods. As a specific application we consider dynamic simulations of rolling bearings and rotor dynamical problems. The work is part of a research collaboration between SKF (www.skf.com) and the Centre of Mathematical Sciences at Lund University
Multisymplectic geometry, variational integrators, and nonlinear PDEs
This paper presents a geometric-variational approach to continuous and
discrete mechanics and field theories. Using multisymplectic geometry, we show
that the existence of the fundamental geometric structures as well as their
preservation along solutions can be obtained directly from the variational
principle. In particular, we prove that a unique multisymplectic structure is
obtained by taking the derivative of an action function, and use this structure
to prove covariant generalizations of conservation of symplecticity and
Noether's theorem. Natural discretization schemes for PDEs, which have these
important preservation properties, then follow by choosing a discrete action
functional. In the case of mechanics, we recover the variational symplectic
integrators of Veselov type, while for PDEs we obtain covariant spacetime
integrators which conserve the corresponding discrete multisymplectic form as
well as the discrete momentum mappings corresponding to symmetries. We show
that the usual notion of symplecticity along an infinite-dimensional space of
fields can be naturally obtained by making a spacetime split. All of the
aspects of our method are demonstrated with a nonlinear sine-Gordon equation,
including computational results and a comparison with other discretization
schemes.Comment: LaTeX2E, 52 pages, 11 figures, to appear in Comm. Math. Phy
Invariant Discretization Schemes Using Evolution-Projection Techniques
Finite difference discretization schemes preserving a subgroup of the maximal
Lie invariance group of the one-dimensional linear heat equation are
determined. These invariant schemes are constructed using the invariantization
procedure for non-invariant schemes of the heat equation in computational
coordinates. We propose a new methodology for handling moving discretization
grids which are generally indispensable for invariant numerical schemes. The
idea is to use the invariant grid equation, which determines the locations of
the grid point at the next time level only for a single integration step and
then to project the obtained solution to the regular grid using invariant
interpolation schemes. This guarantees that the scheme is invariant and allows
one to work on the simpler stationary grids. The discretization errors of the
invariant schemes are established and their convergence rates are estimated.
Numerical tests are carried out to shed some light on the numerical properties
of invariant discretization schemes using the proposed evolution-projection
strategy
- …