Conservation of energy, momentum and actions in numerical discretizations of non-linear wave equations

For classes of symplectic and symmetric time-stepping methods— trigonometric integrators and the Störmer-Verlet or leapfrog method—applied to spectral semi-discretizations of semilinear wave equations in a weakly non-linear setting, it is shown that energy, momentum, and all harmonic actions are approximately preserved over long times. For the case of interest where the CFL number is not a small parameter, such results are outside the reach of standard backward error analysis. Here, they are instead obtained via a modulated Fourier expansion in tim

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 $O(h^p)$-close to the original energy where $p$ is the order of the method and $h$ 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 $O(\exp(-c/h^{1/(1+q)}))$ with $c>0$ and $q \geq 0$; for the semilinear wave equation, $q=1$, and for the nonlinear Schr\"odinger equation, $q=2$. We give an example which shows that analyticity of the initial data is necessary to obtain exponential estimates

Asynchronous Variational Integrators

We describe a new class of asynchronous variational integrators (AVI) for nonlinear elastodynamics. The AVIs are distinguished by the following attributes: (i) The algorithms permit the selection of independent time steps in each element, and the local time steps need not bear an integral relation to each other; (ii) the algorithms derive from a spacetime form of a discrete version of Hamilton’s variational principle. As a consequence of this variational structure, the algorithms conserve local momenta and a local discrete multisymplectic structure exactly. To guide the development of the discretizations, a spacetime multisymplectic formulation of elastodynamics is presented. The variational principle used incorporates both configuration and spacetime reference variations. This allows a unified treatment of all the conservation properties of the system.A discrete version of reference configuration is also considered, providing a natural definition of a discrete energy. The possibilities for discrete energy conservation are evaluated. Numerical tests reveal that, even when local energy balance is not enforced exactly, the global and local energy behavior of the AVIs is quite remarkable, a property which can probably be traced to the symplectic nature of the algorith

Multi-Symplectic Integrators for Nonlinear Wave Equations

Symplectic (area-preserving) integrators for Hamiltonian ordinary differential equations have shown to be robust, efficient and accurate in long-term calculations. In this thesis, we show how symplectic integrators have a natural generalization to Hamiltonian PDEs by introducing the concept of multi-symplectic partial differential equations (PDEs). In particular, we show that multi-symplectic PDEs have an underlying spatio-temporal multi-symplectic structure characterized by a multi-symplectic conservation law MSCL). Then multi-symplectic integrators (MSIs) are numerical schemes that preserve exactly the MSCL. Remarkably, we demonstrate that, although not designed to do so, MSIs preserve very well other associated local conservation laws and global invariants, such as the energy and the momentum, for very long periods of time. We develop two types of MSIs, based on finite differences and Fourier spectral approximations, and illustrate their superior performance over traditional integrators by deriving new numerical schemes to the well known 1D nonlinear Schrödinger and sine-Gordon equations and the 2D Gross-Pitaevskii equation. In sensitive regimes, the spectral MSIs are not only more accurate but are better at capturing the spatial features of the solutions. In particular, for the sine-Gordon equation we show that its phase space, as measured by the nonlinear spectrum associated with it, is better preserved by spectral MSIs than by spectral non-symplectic Runge-Kutta integrators. Finally, to further understand the improved performance of MSIs, we develop a backward error analysis of the multi-symplectic centered-cell discretization for the nonlinear Schrödinger equation. We verify that the numerical solution satisfies to higher order a nearby modified multi-symplectic PDE and its modified multi-symplectic energy conservation law. This implies, that although the numerical solution is an approximation, it retains the key feature of the original PDE, namely its multi-symplectic structure

Variational Structures in Cochain Projection Based Variational Discretizations of Lagrangian PDEs

Compatible discretizations, such as finite element exterior calculus, provide a discretization framework that respect the cohomological structure of the de Rham complex, which can be used to systematically construct stable mixed finite element methods. Multisymplectic variational integrators are a class of geometric numerical integrators for Lagrangian and Hamiltonian field theories, and they yield methods that preserve the multisymplectic structure and momentum-conservation properties of the continuous system. In this paper, we investigate the synthesis of these two approaches, by constructing discretization of the variational principle for Lagrangian field theories utilizing structure-preserving finite element projections. In our investigation, compatible discretization by cochain projections plays a pivotal role in the preservation of the variational structure at the discrete level, allowing the discrete variational structure to essentially be the restriction of the continuum variational structure to a finite-dimensional subspace. The preservation of the variational structure at the discrete level will allow us to construct a discrete Cartan form, which encodes the variational structure of the discrete theory, and subsequently, we utilize the discrete Cartan form to naturally state discrete analogues of Noether's theorem and multisymplecticity, which generalize those introduced in the discrete Lagrangian variational framework by Marsden et al. [29]. We will study both covariant spacetime discretization and canonical spatial semi-discretization, and subsequently relate the two in the case of spacetime tensor product finite element spaces.Comment: 44 pages, 1 figur