ABSTRACTS
Michele Bartuccelli (Surrey) |
On a Class of Integrable Time-Dependent Dynamical Systems
|
We present some integrable time-dependent systems of classical dynamics,
and we apply the
results to the equation \ddot x + f(t) V_x(x) = 0,
with f(t) a positive differentiable function of time, and V(x) is any
smooth potential.
(Joint work with Guido Gentile and Kyriakos Georgiou.)
|
Larry Bates (Calgary) |
The Odd Symplectic Group in Geometry
|
We show how the odd symplectic group arises in contact geometry and the
study of periodic geodesics.
|
Anthony Bloch (Michigan) |
Dissipative Dynamics and Instabilities in Coupled Hamiltonian Systems
|
In this talk I will discuss dissipative behavior and stability
in infinite-dimensional mechanical systems which preserve energy.
In particular I will consider systems of oscillators
interacting with fields in both the classical and quantum settings.
|
Stefanella Boatto (Paris) |
Vortices in the plane and on a sphere: the non-linear
stability of a polygonal ring
|
Since the ninetheen century point-vortex systems have captured the
interest of physicists, celestial mechanicists, fluid mechanicists and
mathematicians.
In 1883 J.J. Thomson used point vortices for one of the first atomic
models (Thomson ``A Treatise on the Motion of Vortex Rings'', Macmillian
(1883)). In particular his interest was in configurations of identical
vortices equally spaced along the circumpherence of a circle, i.e. located
at the vertices of a regular polygon.
In this talk we will discuss the study of the linear and non-linear
stability of polygonal configuration of vortices. The non-linear stability
analysis is carried out by applying a sufficient criterion due to
Dirichlet. Linear and non-linear stability ranges coincide both in
the plane and on the sphere. Furthermore, in both geometries
the addition of a central vortex of strength K strongly modifies
the stability ranges (Joint work with Hildeberto Cabral).
|
Tom Bridges (Surrey) |
Instability and Breakup of Standing Waves (Spatially Periodic Breathers)
of Nonlinear Wave
Equations: Geometric Properties of O(2)-Equivariant Hamiltonian PDEs on
the
Real Line
|
A fundamental class of spatially periodic states of nonlinear wave
equations
on the real line is spatially periodic travelling waves, and
locally (small amplitude) the generic way that these states become unstable
is through a sideband (modulational, Benjamin-Feir, ...) instability.
If the PDE also has a reflection symmetry, it can support spatially
periodic
breathers (SPBs). At the linear level, these states consist of a
superposition of two travelling waves of the same amplitude travelling in
opposite directions. At the nonlinear level, they are spatially periodic
states which are also time-periodic. What is the appropriate
generalisation
of the sideband instability for these states?
Since the NLS equation is the appropriate model for the one-dimensional
case,
one might guess that the appropriate model for analysing the sideband
instability of SPBs is coupled NLS equations. This has been tried and
gives
incorrect results. Edgar Knobloch and collaborators have proposed coupled
NLS (or cGL in the non-conservative case) equations with additional
integral terms, resulting in a partial differential delay equation.
These equations predict the correct form of the sideband instability,
but are for weakly nonlinear SPBs, and do not give a geometric picture
of the nature of the instability. Surprisingly, these are the
only results in the literature about instability of SPBs!
We approach this problem in a new way which gives more complete
information
about the local properties, and extends to finite-amplitude SPBs.
The backbone of the argument is to use the multi-symplectic structures
framework. In this setting such states can be characterised abstractly.
Even in the MSS framework, a straightforward approach fails. Curiously,
the crucial idea is to embed the SPBs in a family of a general
quasiperiodic
two-wave interactions, and then take the limit as the two-waves
synchronise
as an SPB. This might appear to be a dubious limit! But, the two-wave
interaction has a toral (T^2) structure, and the limiting SPB also has a
toral (T^2) structure (because it is time periodic, and there is an S^1
orbit of such waves). This last point becomes to sound plausible when you
think about the analogous situation in finite-dimensions.
Although this analogy can not be stretched too far, an analogue in finite
dimensions is O(2) equivariant Hamiltonian systems, with the spherical
pendulum as prototype. The travelling waves correspond to conical pendulum
solutions and the standing waves correspond to planar pendulum solutions.
There is a torus of planar pendulum solutions obtained by acting with the
circle group on a single planar pendulum solution. Note however, that
there
is no analogue of the sideband instability in finite dimensions.
Back to the PDE: the decomposition of SPBs also sheds light on the nature
of the instability: it turns out that the usual sideband instability
carries
over, but it is doubled. The instability can be associated with a
Hamiltonian Hopf bifurcation with double eigenvalues (i.e. 8 dimensions
associated with the instability rather than 4), and the instability is
in some sense associated with a breakup of the two synchronised waves.
The theory is not rigorous at present, although one is hopeful that parts
of
it can be fully confirmed, but other parts may require special
consideration
due to small divisors. It will also be shown how the theory can be
applied to the instability of standing waves of the water-wave problem.
The above results can be described in some detail, and give a geometrical
picture of spatially-extended O(2)-equivariant Hamiltonian PDEs.
However, one of the most interesting outcomes is that the theory points
the
way towards a framework for more general wave interactions. For example,
how do generalized SPBs, such as spatially periodic patterns on a
hexagonal lattice, become unstable and/or breakup? Some preliminary
results on this problem and their application will be described.
|
Richard Cushman (Utrecht) |
The Rotation Number
and the Herpolhode Angle in the Euler Top
|
This talk will describe a the relation of the rotation
number of a motion of the Euler top on a two torus of constant energy
and length of angular momentum and the angle traced out by
the point of contact of the moment of inertia ellipsoid on
the invariant plane.
|
Gianne Derks (Surrey) |
Travelling Waves in a Singularly Perturbed Sine-Gordon Equation
|
The phase difference of the wave functions of the electrons in the two
superconductors of an ideal long Josephson junction can be modelled by
the sine-Gordon equation. This is a Hamiltonian PDE, which has a family
of stable travelling wave solutions, parametrised by their wave speed. If
a small bias current is applied to Josephson junction, the model has to
be modified by adding a perturbation consisting of both forcing and
dissipative terms, one of which is singular. This perturbation selects a
wave speed in the family of travelling wave solutions and this travelling
wave solution will persist in the perturbed system.
The persisting travelling wave is stable too. This can be shown by
analysing the spectral problem associated with the linearisation about
this wave. The solutions of the spectral problem of the unperturbed
problem are known in the literature. With a socalled Evans function, the
position eigenvalues of the perturbed spectral problem can be determined,
even in the singular case, and it follows that the wave is spectrally
stable. Nonlinear stability can be concluded for the regular case too.
|
Yuri Fedorov (Moscow) |
Integrable Nonholonomic Systems on Lie Groups
|
Integrable systems are usually associated with Hamiltonian ones.
Following V.Kozlov, A. Bloch, A.Veselov and others,
we consider a class of systems on unimodular Lie groups
with nonholonomic constraints, which (at least a priori) are not
Hamiltonian, but possess various tensor invariants, in particular, an
invariant
measure.
We show that for a wide class of metrics such systems are integrable
by quadratures and discuss their relations with various integrable
Hamiltonian systems. In partucular, we prove integrability of the motion of
n-dimensional Euler top under the action of some right-invariant constraints
on the group SO(n).
|
Katrin Gelfert (Dresden) |
All Volume Expanding Dynamical Systems Have Positive Topological Entropy.
|
One major significance of the topological entropy is its strong relation to
other dynamical invariants such as Lyapunov exponents, topological
pressure,
fractal dimension, and Hausdorff dimension, which provides our primary
motivation. Almost all previous investigations of the topological entropy
have been concerned with upper bounds. Exact formulas have been derived
under strong smoothness assumptions only. In this talk we will give lower
bounds of the topological entropy of $C^1$-smooth dynamical systems on
Riemannian manifolds which are sharp in some cases. They are formulated in
terms of the phase space dimension and of the exponential growth rates of a
singular value function of the tangent map. These rates correspond to the
deformation of $k$-volumes and can for instance be estimated in terms of
Lyapunov exponents.
Examples address Henon maps, the Lorenz system, the geodesic flow on a (not
necessarily compact) Riemannian manifold without conjugate points, and
skew
product systems.
|
Mark Groves (Loughborough) |
Nonlinear Water Waves and Spatial Dynamics
|
In this talk I will present a survey of the state of the art
in existence theories for small-amplitude three-dimensional
gravity-capillary water waves which are based upon the
`spatial dynamics' method. The existence theories relate to
travelling waves (which are stationary in a moving frame
of reference) or modulating pulses (which
consist of a permanent pulse-like envelope steadily advancing
in the laboratory frame and modulating an underlying wave-train).
The hydrodynamic problem
is formulated as an infinite-dimensional Hamiltonian
system in which an arbitrary horizontal spatial direction is the
time-like variable; the wave motions are supposed to
be periodic in a second, different horizontal direction.
Travelling waves can be detected using a centre-manifold
reduction technique by which the problem is reduced to a
locally equivalent Hamiltonian system with a finite number
of degrees of freedom. I will show how
virtually any bifurcation or resonance known in Hamiltonian systems theory
is included in the bifurcation catalogue for the reduced Hamiltonian system
and carry out nonlinear bifurcation theory for a representative
selection of bifurcation scenarios. In particular, I will present
existence results for doubly periodic travelling waves, multi-pulse
line solitary waves and further, rather exotic three-dimensional waves.
I will also present some examples in which a complete reduction to
a finite-dimensional reduced Hamiltonian system is not possible.
In particular, existence theories for modulating pulses yield a
reduction of the water-wave problem to a nonlinear wave equation
and existence theories for `periodically modulated' solitary
waves (which have a solitary-wave profile in the direction of
propagation and are periodic in the transverse direction) require
an infinite-dimensional version of the Lyapunov centre theorem.
The question of the existence of a
`fully localised solitary wave' (which decays to zero in all spatial
directions) will also be addressed.
|
Darryl Holm (Los Alamos) |
Nonlinear Balance and Exchange of
Stability in Dynamics of Solitons, Peakons, Ramps/Cliffs and Leftons
in a 1+1 Nonlinear Evolutionary PDE
|
Weak solitary wave solutions resulting from nonlinear
balances in fluid dynamics will be discussed by
dropping linear dispersion terms in an equation that arises in
shallow water waves at the next order beyond KdV.
Specifically, we study
m_t + u m_x + b u_x m = \nu m_{xx}
with momentum $m = u - \alpha^2 u_{xx}$ for fluid
velocity $u(x,t)$. This equation contains three constant
parameters: lengthscale $\alpha$, dimensionless balance parameter
$b$ and Navier-Stokes viscosity $\nu$. As we vary these
parameters, we find numerically and analytically that:
* The solutions of this 1+1 evolutionary PDE that emerge
in
different parameter regimes are pulsons,
ramps-and-cliffs,
peakons, compactons, leftons and solitons that act
like billiard balls, with no internal degrees of freedom
below the length scale $\alpha$.
* These solutions change stability as we vary the
parameter $b$.
They change their behavior at the special values
$b=0,\pm1,\pm2,\pm3$.
* This family of equations admits the classic Burgers
``ramps-and-cliffs'' solutions which are stable for
$-1 < b < 1$ with small viscosity.
For $b=0$ additional special solutions are available.
* For $b<-1$, the Burgers ramps-and-cliffs are unstable.
The stable pulse solutions that emerge for $b<-1$ move
leftward instead of rightward and they tend to a
stationary
pulse train whose individual pulses are each given
when
$\nu=0$ by
$u(x)=u_0sech^2(x/{2\alpha})$ for the case $b=-2$ and
by $u(x)=u_0sech(x/\alpha)$ for the case $b=-3$.
* For $b>1$, the Burgers ramps-and-cliffs are again
unstable. The stable solitary traveling wave that
develops for $b>1$ and $\nu=0$ is the ``peakon''
solution $u(x,t)=ce^{-|x-ct|/\alpha}$. Nonlinear
interactions
among these pulsons or peakons are governed by the
superposition of solutions for $b>1$ and $\nu=0$,
u(x,t) = \sum_{i=1}^N p_i(t) e^{-|x-q_i(t)|/\alpha}
These peakon solutions obey a finite dimensional
dynamical
system for the time-dependent speeds $p_i(t)$ and
positions $q_i(t)$. We study the peakon interactions
analytically, and we determine their fate numerically
under adding viscosity.
|
Andy Hone (Kent) |
Peakon Systems and Poisson Brackets
|
We consider a family of N-body dynamical systems which describe
the motion of peakons (peaked solitons) in a class of partial differential
differential equations which includes the (dispersionless) Camassa-Holm
equation and a new integrable equation isolated by Degasperis and Procesi.
We consider a further generalization to a one-parameter ``pulson'' family
of equations admitting pulses of arbitary symmetric shape. By requiring
that the equations admit a suitable Poisson structure, we obtain a
functional equation which has the peakon kernel as the only non-trivial
solution.
|
Peter Hydon (Surrey) |
The Discrete Variational Complex
|
The variational complex for differential equations is an
extension of the de Rham complex. On topologically trivial domains, this
complex can be used to construct conservation laws systematically
(without using symmetries and Noether's Theorem), and to reconstruct the
Lagrangian from a given system of Euler-Lagrange equations. This talk
describes a new variational complex for difference equations. The complex
is locally exact; in particular, an analogue of the Poincare lemma for
exact differential forms can be stated. Within this framework, we present
homotopy maps that allow one to calculate conservation laws of partial
difference equations and discrete Lagrangians for discrete Euler-Lagrange
systems.
|
Delia Ionescu (Munchen) |
Analysis of the Electrogravitational Kepler Problem
|
In the framework of Classical Mechanics, of General Relativity Theory
and of Relativistic theory of Gravity, the equations governing the trajectories
of charged mass point are given. A comparative analysis of the trajectories in
these theories is presented. It is used the notion of Birkhoffian of a manifold,
from the viewpoint of differential geometry, with the formalism of jets.
|
Toshihiro Iwai (Kyoto) |
Singularity of Many-Body Hamiltonians at Singular Configurations
|
The center-of-mass system of many bodies admits a natural action of the
rotation group SO(3). Singular configurations in the title of this talk
means
those configurations at which the rotation group SO(3) has non-trivial
isotropy
subgroups. Practically, the singular configurations under consideration
are
collinear configurations and a multiple collision.
Many-body Hamiltonians have such singularity at singular configurations
as two-body Hamiltonians do at the origin (ie, two-body collision), when
expressed
in terms spherical polar coordinates. However, the singularity is
not essential, so that the energy integral has finite value. This can be
proved by using the boundary conditions of wave functions at singular
configurations.
|
Bozidar Jovanovic (SANU, Belgrade) |
Reduction and Integrability
|
Suppose we are given a compact Riemannian manifold (Q,g) with
a completely integrable geodesic flow. Let G be a compact connected Lie
group acting freely on Q by isometries. The natural question arises:
will the geodesic flow on Q/G equipped with the submersion metric be
integrable? Under one natural assumption, we prove that the answer is
affirmative. New examples of manifolds with completely integrable geodesic
flows are obtained. [math-ph/0204048, to appear in Lett. Math. Phys.]
|
Melvin Leok (Caltech) |
A Discrete Theory of Connections on Principal Bundles
|
Motivated by applications to Discrete Lagrangian Reduction, we consider the
discrete analogue of the Atiyah sequence of a principal bundle, and relate a
splitting of the discrete Atiyah sequence with discrete horizontal lifts and
discrete connection forms. Continuous connections can be obtained by taking
the limit of discrete connections in a natural way.
This yields an isomorphism between (Q\times Q)/G and \tilde{G}\oplus(S\times
S). Both the discrete connection and the associated continuous connection
are necessary to express the Discrete Lagrange-Poincare operator in
coordinates.
|
Anna Litvak Hinenzon (Warwick) |
On Energy Surfaces and Instabilities
|
Geometrical study of energy surfaces of integrable Hamiltonian systems
lead
to upper bounds on the possible instability rate in the system after a
small
perturbation, which brakes the integrability, is applied.
This will be demonstrated using an example Hamiltonian with
degeneracies.
|
Karsten Matthies (Berlin/Warwick) |
Exponential Averaging for Hamiltonian Systems
|
We derive estimates on the
magnitude of non-adiabatic interaction between
a Hamiltonian partial differential equation and a high-frequency
nonlinear oscillator. Assuming spatial analyticity of the initial
conditions, we show that the dynamics can be transformed to
the uncoupled dynamics of an infinite-dimensional Hamiltonian system
and an anharmonic oscillator, up to coupling terms which are
exponentially small in a certain power of the frequency of the
oscillator. The result is derived from an abstract averaging theorem
for infinite-dimensional analytic evolution equations in Gevrey spaces.
An application is provided by a system of Nonlinear
Schr\"odinger Equations, coupled to a rapidly forcing single mode,
representing small-scale oscillations. (Joint work with Arnd Scheel).
|
Jerry Marsden (Caltech) |
Controlled Lagrangian and Hamiltonian Systems with Symmetry
|
The first part of the lecture will give the general
methodology of controlled Lagrangians methodology (work with Bloch and
Leonard) along with some examples of how it works in concrete
examples. We will then generalize this to a more general notion of
equivalence of controlled Lagrangian and Hamiltonian systems and prove
the equivalence of the Lagrangian and the Hamiltonian approaches (work
with Bloch, Chang, Leonard and Woolsey). The last part of the talk
(work with Chang, Bloch and Leonard) will give a general reduction
theory for controlled Lagrangian and controlled Hamiltonian systems
with symmetry. Reduction theory for these systems is needed in many
of the fundamental examples, such as a spacecraft with rotors and
underwater vehicle dynamics. We also show the equivalence of the
method of reduced controlled Lagrangian systems and that of reduced
controlled Hamiltonian systems for simple mechanical systems with
symmetry.
|
Chris McCord (Cincinatti) |
Integral Manifolds of the N-Body Problem
|
This talk will survey the efforts of the last decade to
understand the topology of the integral manifolds, and how that
topology constrains the dynamics of the N-body problem. We use
a homological approach to describe the manifolds, detect bifurcations
as the energy and angular momentum vary, and identify constraints on
the dynamics imposed by the homology. We will look at what has been
done, work in progress, and future prospects for this approach.
|
Marcel Oliver (Tubingen) |
Variational Asymptotics for Rotating Fluids near Geostrophy
|
Large scale flow in the atmosphere and in the ocean at mid-latitudes
can often be considered as nearly geostrophic, i.e. the pressure
gradient force is almost balanced by the Coriolis force, while inertia
is a lower order effect. It is therefore desirable to find reduced
models that represent large scale flows near geostrophy while being
computationally less expensive than the parent model. One important
consideration is the need to remove inertia-gravity waves from the
model, because such waves do not interact strongly with large scale
circulation, yet force severe step-size constraints upon direct
numerical simulations.
There are a number of well known models for geostrophic flow which
differ in the assumed parent model, in the precise scaling
assumptions, and the choice of coordinate system. It is often not
well understood how these different models relate to each other.
In this talk I will demonstrate an ansatz based on variational
asymptotics---involving an infinitesimal change of coordinates in the
low-Rossby number expansion of the variational principle---which can
reproduce a number of known models as well as new ones. This
formulation makes it easy distinguish models based on their analytical
properties (their well-posedness and regularity) and choose those most
consistent with the model setting: Vortical motion at large scales.
|
Maxim Pavlov (Moscow) |
Hydrodynamic Integrable Chains
|
The history of Benney systems has set of splendid steps obtained by namely
D.Benney, Y.Manin, B.Kupershmidt, V.Zakharov and J.Gibbons. Later was
showed that Benney system (in Zakharov reduction) can be interpreted as
dispersionless limit of KP (Kadomtzev-Petviashvili) hierarchy. Thus, most
interesting question is description of all possible reductions, and
concequently
corresponding Hamiltonian structures, symmetries and conservation laws.
We suggest approach to construct new integrable chains as Benney momentum
chain starting from arbitrary given integrable hydrodynamic type system.
|
Milena Radnovic (SANU, Belgrade) |
Caustics of Elliptical Billiard Trajectories
|
Cayley's type conditions for periodical trajectories of the billiard
within an ellipsoid in the Eucledean or Lobachevsky space of arbitrary
dimension are derived. Singular cases and their caustics are additionaly
discussed. We also show that Lobachevsky and Eucledean elliptic billiards
can be naturally considered as a part of a hierarchy of integrable
elliptical billiards. The results are joint with V.Dragovic and
B.Jovanovic.
|
Renzo Ricca (UCL) |
On Kelvin's Vortex Knots
|
In a memorable paper Lord Kelvin (1875) addressed the question of
evolution and stability of thin core vortex torus knots under
Euler's equations. The question remained open and has only recently
attracted new interest. We briefly outline some of our new
results (in the Localised Induction Approximation context) and
report on current work (done in collaboration with Xinyu He) to
tackle the original problem. In particular we show how geometric
and topological properties measured by the winding number, and
symmetry aspects, influence the evolution of these knots.
|
Jedrzej Sniatycki (Calgary) |
Poisson Reduction
|
I shall discuss successes and shortcomings of the program of Poisson reduction
of symmetries of Hamiltonian systems.
|
Britta Sommer (Aachen) |
Not So Proper Degeneracy of Systems Living According to Three
Different Time Scales - A KAM Theorem for the Spatial Lunar Problem
|
The classical KAM Theorem applies to perturbed non-degenerate
Hamiltonian systems. The frequency map of the unperturbed
Hamiltonian is a local diffeomorphism and all frequencies are of
the same order of magnitude. The persistence of invariant
diophantine tori can also be proven, if the unperturbed system is
a so called properly degenerate one. In this case, the
perturbation "removes the degeneracy". A new "unperturbed
Hamiltonian" is formed whose frequencies are of order 1 and of
order \varepsilon, where \varepsilon is the perturbation
parameter.
When we look at the the movement of a small body in the
gravitational potential of a planet and the sun, we encounter a
system that is not even properly degenerate. Usually, this system is
approximated by a perturbed Keplerian system and denoted as the
spatial Lunar problem. It turns out to be a three degree of
freedom system. Because of the superintegrability of the Kepler
Hamiltonian its unperturbed part does depend on one of three
actions only. The first order term of its perturbation does not
depend on all actions either. Thus, motion takes place according
to three different time scales.
In this talk we will show that there exist crucial similarities
between this kind of system and properly degenerate ones. This
leads to the formulation of an appropriate KAM Theorem. A short
description of the problems that have to be overcome when applying
the Theorem to the lunar problem finishes this talk.
|
Tadashi Tokieda (Montreal) |
Standing Vortices, Fidgeting Vortices, Dancing Vortices
|
''I would like to explain several
methods for finding (relative) equilibria and (relative) periodic
orbits of vortices''.
|
Alex Tsygvintsev (Loughborough) |
Integral invariants of H.Poincare
|
In this talk we remind the basic theory of integral invariants
created by H.Poncare. We describe the four principal types of integral
invariants, the relation with variational equations and sketch shortly the
application to three-body problem.
As a new result we give a non-academic example of hierarchy of integral
invariants arising in one problem of mechanics.
|
Jacques Vanneste (Edinburgh) |
Adiabatic invariance
and geometric angle for fluids in deforming domains
|
The response of a two-dimensional, incompressible and inviscid fluid
to a prescribed deformation of its boundary is investigated
under the assumption that the boundary deformation is slow compared to
the fluid motion. The evolution of both the Eulerian flow fields
and the Lagrangian particle positions is studied. The analogy of
the problem with finite-dimensional, near-integrable Hamiltonian
systems with slowly varying parameters is emphasised.
|
Cornelia Vizman (Timisoara) |
Central Extensions of Lie algebras of Symplectic Vector Fields
|
Using the result of Roger that the second cohomology group
of the Lie algebra of Hamiltonian vector fields on a closed symplectic
manifold is isomorphic to the first deRham cohomology group of the
manifold, we determine with the help of the Hochschild-Serre spectral
sequence the second cohomology group of the Lie algebra of symplectic
vector fields.
|
Matt West (Caltech) |
Geometric Collision Integrators
|
Mechanical systems with collisions pose a number of problems for standard
theoretical and numerical techniques due to their extremely nonsmooth
behavior. In this talk we formulate variational collision problems in both
continuous and discrete time. This yields an interesting class of
integrators which preserve well-defined geometric structures from the
continuous time problem. We use these as a basis to construct a highly
efficient explicit collision algorithm, which we demonstrate on a thin-shell
collision problem.
|
Claudia Wulff (Berlin/Warwick/Surrey) |
Approximate Momentum Conservation for Spatial Semidiscretizations
of Nonlinear Wave Equations
|
We prove that for the standard second order finite difference uniform
space discretization of the nonlinear wave equation with analytic
nonlinearity and analytic initial values momentum is conserved up to an
error which is exponentially small in the stepsize. We provide rigorous
estimates which are valid as long as the trajectories remain
bounded in a Gevrey space. The method of proof is that of
backward error analysis. We construct a modified equation which
interpolates the discretized system. For the standard finite difference
discretization scheme this modified system has translation symmetry and
is Hamiltonian and therefore conserves momentum. We show that both
conditions are essential for approximate momentum conservation of the
discretized system.
Joint work with Marcel Oliver (Tubingen/Bremen) and Matt West (CalTech).
|