Symmetry-preserving numerical schemes

Alexander Bihlo,Francis Valiquette
DOI: https://doi.org/10.48550/arXiv.1608.02557
2016-12-07
Abstract:In these lectures we review two procedures for constructing finite difference numerical schemes that preserve symmetries of differential equations. The first approach is based on Lie's infinitesimal symmetry generators, while the second method uses the novel theory of equivariant moving frames. The advantages of both techniques are discussed and illustrated with the Schwarzian differential equation, the Korteweg-de Vries equation and Burgers' equation. Numerical simulations are presented and innovative techniques for obtaining better invariant numerical schemes are introduced. New research directions and open problems are indicated at the end of these notes.
Numerical Analysis,Mathematical Physics
What problem does this paper attempt to address?
The core problem that this paper attempts to solve is to develop numerical integration methods that can preserve the symmetries of differential equations. Specifically, the authors focus on how to construct finite - difference numerical schemes that can maintain the symmetric properties of the original differential equations during the discretization process. The importance of this goal lies in the fact that by preserving the geometric properties of the system, the long - term accuracy of numerical results can be improved, especially when dealing with physical systems with specific symmetries. The paper discusses two main methods to achieve this goal: 1. **The method based on infinitesimal generators of Lie groups**: This method uses the infinitesimal generators of Lie groups to construct difference invariants and then build symmetry - preserving numerical schemes. In this way, it can be ensured that the numerical scheme converges to the original differential equation in the continuous limit and preserves its symmetry. 2. **The equivariant moving frame method**: This is a relatively new method that constructs difference invariants through the equivariant moving frame. Compared with the traditional infinitesimal generator method, this method provides a more systematic way to find suitable combinations of difference invariants, thereby constructing symmetry - preserving numerical schemes. The paper also explores the advantages of these two methods and conducts numerical simulations through specific examples (such as the Schwarzian differential equation, the Korteweg - de Vries equation, and the Burgers equation) to show the performance of these symmetry - preserving numerical schemes in practical applications. In addition, the authors also propose some innovative techniques to further improve the symmetry - preserving numerical schemes, especially for the grid - winding problem encountered when dealing with partial differential equations. In summary, this paper aims to show, through theoretical analysis and numerical experiments, how to effectively construct and apply symmetry - preserving numerical schemes, thereby improving the accuracy and stability of numerical simulations.