Periodic Reporting for period 2 - LAHACODE (Low-regularity and high oscillations: numerical analysis and computation of dispersive evolution equations)
Reporting period: 2021-11-01 to 2023-04-30
We developed a general framework of low regularity integrators for an abstract evolution equation which involves not only dispersive type equations, but also hyperbolic problems and parabolic problems (such as the wave equation, nonlinear heat equation, the complex Ginzburg-Landau equation, the half wave and Klein--Gordon equations). The main idea is based on new oscillatory integrators which allow us to embed the underlying oscillations of the partial differential equation into the numerical discretisation. This result was published in SIAM J. Num. Analysis. Furthermore, we could develop a new type of low-regularity integrator for Navier-Stokes equations, coupled with a stabilized finite element method in space. Unlike the other low-regularity integrators for nonlinear dispersive equations, which are all fully explicit in time, the proposed method is semi-implicit in time in order to preserve the energy-decay structure of NS equations. First-order convergence of the proposed method is established independent of the viscosity coefficient μ, under weaker regularity conditions than other existing numerical methods, including the semi-implicit Euler method and classical exponential integrators.
We introduced discrete Bourgain spaces, which allow us to establish low regularity error estimates for the nonlinear Schrödinger equation under periodic boundary condition. This result will appear in the Journal of the European Mathematical Society (JEMS). We could furthermore establish sharp L^2 error estimates for the periodic Korteweg--de Vries equation and, with the aid of discrete Strichartz type estimates, for the nonlinear Schrödinger equation on the full space in up to three spatial dimensions. The latter results are published in Foundations of Computational Mathematics and Mathematics of Computations.
Recent years have seen an increasing amount of research devoted to the development of so-called resonance-based methods for dispersive nonlinear partial differential equations. In many situations, this new class of methods allows for approximations in a much more general setting (e.g. for rough data) than, for instance, classical splitting or exponential integrator methods. However, they lack one important property: the preservation of geometric structures. This is particularly drastic in the case of the Korteweg--de Vries (KdV) equation which is a fundamental model in the broad field of dispersive equations that is completely integrable, possessing infinitely many conserved quantities, an important property which we wish to capture -- at least up to some degree -- also on the discrete level. A revolutionary step in this direction was set by the theory of geometric numerical integration resulting in the development of a wide range of structure-preserving algorithms for Hamiltonian systems. However, in general, these methods rely heavily on highly regular solutions. State-of-the-art low-regularity integrators, on the other hand, poorly preserve the geometric structure of the underlying PDE. We made a first step towards bridging the gap between low regularity and structure preservation. We introduced a novel symplectic (in the Hamiltonian picture) resonance-based method on the example of the KdV equation that allows for low-regularity approximations to the solution while preserving the underlying geometric structure of the continuous problem on the discrete level. This is a joint work with Georg Maierhofer, see also arxiv https://arxiv.org/abs/2205.05024
We introduced a general framework of low regularity integrators which allows us to approximate the time dynamics of a large class of equations, including parabolic and hyperbolic problems, as well as dispersive equations, up to arbitrary high order on general domains. The structure of the local error of the new schemes is driven by nested commutators which in general require (much) lower regularity assumptions than classical methods do. Our main idea lieds in embedding the central oscillations of the nonlinear PDE into the numerical discretisation. The latter is achieved by a novel decorated tree formalism inspired by singular SPDEs with Regularity Structures and allows us to control the nonlinear interactions in the system up to arbitrary high order on the infinite dimensional (continuous) as well as finite dimensional (discrete) level. This is a joint work with Yvonne Alama Bronsard and Yvain Bruned, see also https://arxiv.org/abs/2202.01171
We introduced a new class of numerical schemes which allow for low regularity approximations to the expectation. This quantity plays an important role in physics, in particular in the study of wave turbulence where one needs to adopt a statistical approach in order to obtain deep insight into the generic long-time behaviour of solutions to dispersive equations. Our new class of schemes is based on Wick's theorem and Feynman diagrams together with a resonance based discretisation set in a more general context: we introduce a novel combinatorial structure called paired decorated forests which are two decorated trees whose decorations on the leaves come in pair. The character of the scheme draws its inspiration from the treatment of singular stochastic partial differential equations via Regularity Structures. In contrast to classical approaches, we do not discretize the PDE itself, but rather its expectation. This allows us to heavily exploit the optimal resonance structure and underlying gain in regularity on the finite dimensional (discrete) level. This is a joint work with Yvonne Alama Bronsard and Yvain Bruned, see also https://arxiv.org/abs/2205.02156.