22 research outputs found
Exponentially accurate Hamiltonian embeddings of symplectic A-stable Runge--Kutta methods for Hamiltonian semilinear evolution equations
We prove that a class of A-stable symplectic Runge--Kutta time
semidiscretizations (including the Gauss--Legendre methods) applied to a class
of semilinear Hamiltonian PDEs which are well-posed on spaces of analytic
functions with analytic initial data can be embedded into a modified
Hamiltonian flow up to an exponentially small error. As a consequence, such
time-semidiscretizations conserve the modified Hamiltonian up to an
exponentially small error. The modified Hamiltonian is -close to the
original energy where is the order of the method and the time
step-size. Examples of such systems are the semilinear wave equation or the
nonlinear Schr\"odinger equation with analytic nonlinearity and periodic
boundary conditions. Standard Hamiltonian interpolation results do not apply
here because of the occurrence of unbounded operators in the construction of
the modified vector field. This loss of regularity in the construction can be
taken care of by projecting the PDE to a subspace where the operators occurring
in the evolution equation are bounded and by coupling the number of excited
modes as well as the number of terms in the expansion of the modified vector
field with the step size. This way we obtain exponential estimates of the form
with and ; for the semilinear wave
equation, , and for the nonlinear Schr\"odinger equation, . We give
an example which shows that analyticity of the initial data is necessary to
obtain exponential estimates
Resonances in long time integration of semi linear Hamiltonian PDEs
We consider a class of Hamiltonian PDEs that can be split into a linear
unbounded operator and a regular non linear part, and we analyze their
numerical discretizations by symplectic methods when the initial value is small
in Sobolev norms. The goal of this work is twofold: First we show how standard
approximation methods cannot in general avoid resonances issues, and we give
numerical examples of pathological behavior for the midpoint rule and
implicit-explicit integrators. Such phenomena can be avoided by suitable
truncations of the linear unbounded operator combined with classical splitting
methods. We then give a sharp bound for the cut-off depending on the time step.
Using a new normal form result, we show the long time preservation of the
actions for such schemes for all values of the time step, provided the initial
continuous system does not exhibit resonant frequencies
Energy-conserving methods for the nonlinear Schrödinger equation
In this paper, we further develop recent results in the numerical solution of Hamiltonian partial differential equations (PDEs) (Brugnano et al., 2015), by means of energy-conserving methods in the class of Line Integral Methods, in particular, the RungeâKutta methods named Hamiltonian Boundary Value Methods (HBVMs). We shall use HBVMs for solving the nonlinear Schrödinger equation (NLSE), of interest in many applications. We show that the use of energy-conserving methods, able to conserve a discrete counterpart of the Hamiltonian functional, confers more robustness on the numerical solution of such a problem
Metastable energy strata in numerical discretizations of weakly nonlinear wave equations
The quadratic nonlinear wave equation on a one-dimensional torus with small
initial values located in a single Fourier mode is considered. In this
situation, the formation of metastable energy strata has recently been
described and their long-time stability has been shown. The topic of the
present paper is the correct reproduction of these metastable energy strata by
a numerical method. For symplectic trigonometric integrators applied to the
equation, it is shown that these energy strata are reproduced even on long time
intervals in a qualitatively correct way.Comment: 28 pages, 9 figure