Benoit Grebert

Erwan Faou, Benoit Grebert, Eric Paturel

We consider {\em discretized} Hamiltonian PDEs associated with a Hamiltonian function that can be split into a linear unbounded operator and a regular nonlinear part. We consider splitting methods associated with this decomposition. Using a finite dimensional Birkhoff normal form result, we show the almost preservation of the {\em actions} of the numerical solution associated with the splitting method over arbitrary long time, provided the Sobolev norms of the initial data is small enough, and for asymptotically large level of space approximation. This result holds under {\em generic} non resonance conditions on the frequencies of the linear operator and on the step size. We apply this results to nonlinear Schrödinger equations as well as the nonlinear wave equation.}

Erwan Faou, Benoit Grebert, Eric Paturel

We consider Hamiltonian PDEs that can be split into a linear unbounded operator and a regular non linear part. We consider abstract splitting methods associated with this decomposition where no discretization in space is made. We prove a normal form result for the corresponding discrete flow under generic non resonance conditions on the frequencies of the linear operator and on the step size. This result implies the conservation of the regularity of the numerical solution associated with the splitting method over arbitrary long time, provided the initial data is small enough. This result holds for numerical schemes controlling the round-off error at each step to avoid possible high frequency energy drift. We apply this results to nonlinear Schrödinger equations as well as the nonlinear wave equation.}

Dario Bambusi, Erwan Faou, Benoit Grebert

In this paper we study the long time behavior of a discrete approximation in time and space of the cubic nonlinear Schrödinger equation on the real line. More precisely, we consider a symplectic time splitting integrator applied to a discrete nonlinear Schrödinger equation with additional Dirichlet boundary conditions on a large interval. We give conditions ensuring the existence of a numerical soliton which is close in energy norm to the continuous soliton. Such result is valid under a CFL condition between the time and space stepsizes. Furthermore we prove that if the initial datum is symmetric and close to the continuous soliton, then the associated numerical solution remains close to the orbit of the continuous soliton for very long times.

Erwan Faou, Benoit Grebert

We consider a wide class of semi linear Hamiltonian partial differential equa- tions and their approximation by time splitting methods. We assume that the nonlinearity is polynomial, and that the numerical tra jectory remains at least uni- formly integrable with respect to an eigenbasis of the linear operator (typically the Fourier basis). We show the existence of a modified interpolated Hamiltonian equation whose exact solution coincides with the discrete flow at each time step over a long time depending on a non resonance condition satisfied by the stepsize. We introduce a class of modified splitting schemes fulfilling this condition at a high order and prove for them that the numerical flow and the continuous flow remain close over exponentially long time with respect to the step size. For stan- dard splitting or implicit-explicit scheme, such a backward error analysis result holds true on a time depending on a cut-off condition in the high frequencies (CFL condition). This analysis is valid in the case where the linear operator has a discrete (bounded domain) or continuous (the whole space) spectrum.

Erwan Faou, Benoit Grebert

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.