Multisymplectic integrator explained

In mathematics, a multisymplectic integrator is a numerical method for the solution of a certain class of partial differential equations, that are said to be multisymplectic. Multisymplectic integrators are geometric integrators, meaning that they preserve the geometry of the problems; in particular, the numerical method preserves energy and momentum in some sense, similar to the partial differential equation itself. Examples of multisymplectic integrators include the Euler box scheme and the Preissman box scheme.

Multisymplectic equations

A partial differential equation (PDE) is said to be a multisymplectic equation if it can be written in the form

Kz _{t} + Lz _{x} = \nabla S(z), where z(t,x) is the unknown, K and L are (constant) skew-symmetric matrices and \nabla S denotes the gradient of S .^[1] This is a natural generalization of Jz _{t} = \nabla H(z) , the form of a Hamiltonian ODE.^[2] u _{tt} - u _{xx} = V'(u) , or more generally the nonlinear wave equation u _{tt} = \partial _{x} \sigma'(u _{x)} - f'(u) ,^[3] and the KdV equation u _{t} + uu _{x} + u _{xxx} = 0 .^[4]

Define the 2-forms

\omega and \kappa by \omega(u,v) = \langle Ku, v \rangle and \kappa(u,v) = \langle Lu, v \rangle where \langle ⋅ , ⋅ \rangle denotes the dot product. The differential equation preserves symplecticity in the sense that\partial _{t} \omega + \partial _{x} \kappa = 0 . ^[5] Taking the dot product of the PDE with u _{t} yields the local conservation law for energy:\partial _{t} E(u) + \partial _{x} F(u) = 0 where E(u) = S(u) - \tfrac 12 \kappa(u _{x,u)} , F(u) = \tfrac 12 \kappa(u _{t,u).} ^[6] The local conservation law for momentum is derived similarly:\partial _{t} I(u) + \partial _{x} G(u) = 0 where I(u) = \tfrac 12 \omega(u _{x,u)} , G(u) = S(u) - \tfrac 12 \omega(u _{t,u).}

The Euler box scheme

A multisymplectic integrator is a numerical method for solving multisymplectic PDEs whose numerical solution conserves a discrete form of symplecticity.^[7] One example is the Euler box scheme, which is derived by applying the symplectic Euler method to each independent variable.^[8]

The Euler box scheme uses a splitting of the skewsymmetric matrices

K and L of the form:\begin{align} K & = K _{+} + K _{-} with K _{-} = - K ^{T,} _{+} \\ L & = L _{+} + L _{-} with L _{-} = - L ^{T.} _{+} \end{align} For instance, one can take K _{+} and L _{+} to be the upper triangular part of K and L , respectively.^[9]

Now introduce a uniform grid and let

u _{n,i} denote the approximation to u(n\Delta{t}, i\Delta{x}) where \Delta{t} and \Delta{x} are the grid spacing in the time- and space-direction. Then the Euler box scheme is K _{+} \partial ^{+} _{t} u _{n,i} + K _{-} \partial ^{-} _{t} u _{n,i} + L _{+} \partial ^{+} _{x} u _{n,i} + L _{-} \partial ^{-} _{x} u _{n,i} = \nabla{S}(u _{n,i} ) where the finite difference operators are defined by\begin{align} \partial ^{+} _{t} u _{n,i} & = (u _{n+1,i} - u _{n,i})/(\Delta{t) }, & \partial_x^+ u_ &= \frac, \\[1ex]\partial_t^- u_ &= \frac, & \partial_x^- u_ &= \frac.\end ^[10] The Euler box scheme is a first-order method, which satisfies the discrete conservation law\partial ^{+} _{t} \omega _{n,i} + \partial ^{+} _{x} \kappa _{n,i} = 0 where \omega _{n,i} = d u _{n,i-1} \wedge K _{+} d u _{n,i} and \kappa _{n,i} = d u _{n-1,i} \wedge L _{+} d u _{n,i} . ^[11]

Preissman box scheme

Another multisymplectic integrator is the Preissman box scheme, which was introduced by Preissman in the context of hyperbolic PDEs.^[12] It is also known as the centred cell scheme.^[13] The Preissman box scheme can be derived by applying the Implicit midpoint rule, which is a symplectic integrator, to each of the independent variables.^[14] This leads to the scheme

K \partial ^{+} _{t} u _{n,i+1/2} + L \partial ^{+} _{x} u _{n+1/2,i} = \nabla{S}(u _{n+1/2,i+1/2} ), where the finite difference operators \partial ^{+} _{t} and \partial ^{+} _{x} are defined as above and the values at the half-integers are defined by u _{n,i+1/2} = (u _{n,i} + u _{n,i+1})/(2) , u _{n+1/2,i} = (u _{n,i} + u _{n+1,i})/(2) , u _{n+1/2,i+1/2} = (u _{n,i} + u _{n,i+1} + u _{n+1,i} + u _{n+1,i+1})/(4) . The Preissman box scheme is a second-order multisymplectic integrator which satisfies the discrete conservation law\partial ^{+} _{t} \omega _{n,i} + \partial ^{+} _{x} \kappa _{n,i} = 0 where \omega _{n,i} = d u _{n,i+1/2} \wedge K d u _{n,i+1/2} and \kappa _{n,i} = d u _{n+1/2,i} \wedge L d u _{n+1/2,i} . ^[15]

References

• .
• .
• .
• .
• .
• .

Notes and References

1. , p. 1374;, p. 335–336.
2. , p. 186.
3. , p. 335.
4. , p. 339–340.
5. , p. 186;, p. 336.
6. , p. 187;, p. 337–338.

7. .

8. .

9. .

10. .

11. .

12. refers to for the work by Preissman.
13. .

14. .

15. .