Frédéric Rousset

Alexander Ostermann, Frédéric Rousset, Katharina Schratz

We present a new filtered low-regularity Fourier integrator for the cubic nonlinear Schrödinger equation based on recent time discretization and filtering techniques. For this new scheme, we perform a rigorous error analysis and establish better convergence rates at low regularity than known for classical schemes in the literature so far. In our error estimates, we combine the better local error properties of the new scheme with a stability analysis based on general discrete Strichartz-type estimates. The latter allow us to handle a much rougher class of solutions as the error analysis can be carried out directly at the level of $L^2$ compared to classical results \black in dimension $d$, \black which are limited to higher-order (sufficiently smooth) Sobolev spaces $H^s$ with $s>d/2$. In particular, we are able to establish a global error estimate in $L^2$ for $H^1$ solutions which is roughly of order $τ^{ {1\over 2} + { 5-d \over 12} }$ in dimension $d \leq 3$ ($τ$ denoting the time discretization parameter). This breaks the "natural order barrier" of $τ^{1/2}$ for $H^1$ solutions which holds for classical numerical schemes (even in combination with suitable filter functions).

Clémentine Courtès, Frédéric Lagoutière, Frédéric Rousset

This article deals with the numerical analysis of the Cauchy problem for the Korteweg-de Vries equation with a finite difference scheme. We consider the Rusanov scheme for the hyperbolic flux term and a 4-points $θ$-scheme for the dispersive term. We prove the convergence under a hyperbolic Courant-Friedrichs-Lewy condition when $θ\geq \frac{1}{2}$ and under an "Airy" Courant-Friedrichs-Lewy condition when $θ<\frac{1}{2}$. More precisely, we get the first order convergence rate for strong solutions in the Sobolev space $H^s(\mathbb{R})$, $s \geq 6$ and extend this result to the non-smooth case for initial data in $H^s(\mathbb{R})$, with $s\geq \frac{3}{4}$ , to the price of a loss in the convergence order. Numerical simulations indicate that the orders of convergence may be optimal when $s\geq3$.

Katie Marsden, Frédéric Rousset, Katharina Schratz

We prove convergence of a filtered Lie splitting scheme for the wave maps equation with low regularity initial data in dimension 3. The convergence analysis is performed in discrete Bourgain spaces, as has proved fruitful for the low regularity analysis of the equation in the continuous setting. An important difficulty here is that the analysis of wave maps at low regularity requires the use of the null structure of the system, this structure thus has to be preserved at the discrete level to get an effective stable low regularity scheme. Since the null structure involves time derivatives, the scheme has to be designed carefully. The presence of time derivatives in the nonlinearity then constitutes the most significant source of numerical error. Nonetheless, we are able to prove convergence of the scheme for all subcritical initial data in $H^s$, $s>d/2$.

Ludwig Gauckler, Jianfeng Lu, Jeremy L. Marzuola et al.

Trigonometric time integrators are introduced as a class of explicit numerical methods for quasilinear wave equations. Second-order convergence for the semi-discretization in time with these integrators is shown for a sufficiently regular exact solution. The time integrators are also combined with a Fourier spectral method into a fully discrete scheme, for which error bounds are provided without requiring any CFL-type coupling of the discretization parameters. The proofs of the error bounds are based on energy techniques and on the semiclassical Gårding inequality.

Erwan Faou, Romain Horsin, Frédéric Rousset

We consider time discretizations of the Vlasov-HMF (Hamiltonian Mean-Field) equation based on splitting methods between the linear and non-linear parts. We consider solutions starting in a small Sobolev neighborhood of a spatially homogeneous state satisfying a linearized stability criterion (Penrose criterion). We prove that the numerical solutions exhibit a scattering behavior to a modified state, which implies a nonlinear Landau damping effect with polynomial rate of damping. Moreover, we prove that the modified state is close to the continuous one and provide error estimates with respect to the time stepsize.