195 research outputs found
Backward error analysis for multisymplectic discretizations of Hamiltonian PDEs
Several recently developed multisymplectic schemes for Hamiltonian PDEs have
been shown to preserve associated local conservation laws and constraints very
well in long time numerical simulations. Backward error analysis for PDEs, or
the method of modified equations, is a useful technique for studying the
qualitative behavior of a discretization and provides insight into the
preservation properties of the scheme. In this paper we initiate a backward
error analysis for PDE discretizations, in particular of multisymplectic box
schemes for the nonlinear Schrodinger equation. We show that the associated
modified differential equations are also multisymplectic and derive the
modified conservation laws which are satisfied to higher order by the numerical
solution. Higher order preservation of the modified local conservation laws is
verified numerically.Comment: 12 pages, 6 figures, accepted Math. and Comp. Simul., May 200
Geometric numerical schemes for the KdV equation
Geometric discretizations that preserve certain Hamiltonian structures at the
discrete level has been proven to enhance the accuracy of numerical schemes. In
particular, numerous symplectic and multi-symplectic schemes have been proposed
to solve numerically the celebrated Korteweg-de Vries (KdV) equation. In this
work, we show that geometrical schemes are as much robust and accurate as
Fourier-type pseudo-spectral methods for computing the long-time KdV dynamics,
and thus more suitable to model complex nonlinear wave phenomena.Comment: 22 pages, 14 figures, 74 references. Other author's papers can be
downloaded at http://www.lama.univ-savoie.fr/~dutykh
Simple bespoke preservation of two conservation laws
Conservation laws are among the most fundamental geometric properties of a partial differential equation (PDE), but few known finite difference methods preserve more than one conservation law. All conservation laws belong to the kernel of the Euler operator, an observation that was first used recently to construct approximations symbolically that preserve two conservation laws of a given PDE. However, the complexity of the symbolic computations has limited the effectiveness of this approach. The current paper introduces some key simplifications that make the symbolic-numeric approach feasible. To illustrate the simplified approach, we derive bespoke finite difference schemes that preserve two discrete conservation laws for the Korteweg-de Vries (KdV) equation and for a nonlinear heat equation. Numerical tests show that these schemes are robust and highly accurate compared to others in the literature
Multi-Symplectic Simulation on Soliton-Collision for Nonlinear Perturbed Schrödinger Equation
Seeking solitary wave solutions and revealing their interactional characteristics for nonlinear evolution equations help us lot to comprehend the motion laws of the microparticles. As a local nonlinear dynamic behavior, the soliton-collision is difficult to be reproduced numerically. In this paper, the soliton-collision process in the nonlinear perturbed Schrödinger equation is simulated employing the multi-symplectic method. The multi-symplectic formulations are derived including the multi-symplectic form and three local conservation laws of the nonlinear perturbed Schrödinger equation. Employing the implicit midpoint rule, we construct a multi-symplectic scheme, which is equivalent to the Preissmann box scheme, for the nonlinear perturbed Schrödinger equation. The elegant structure-preserving properties of the multi-symplectic scheme are illustrated by the tiny maximum absolute residual of the discrete multi-symplectic structure at each time step in the numerical simulations. The effects of the perturbation strength on the soliton-collision in the nonlinear perturbed Schrödinger equation are reported in the numerical results in detail
Geometric Numerical Integration (hybrid meeting)
The topics of the workshop
included interactions between geometric numerical integration and numerical partial differential equations;
geometric aspects of stochastic differential equations;
interaction with optimisation and machine learning;
new applications of geometric integration in physics;
problems of discrete geometry, integrability, and algebraic aspects
Small Collaboration: Numerical Analysis of Electromagnetic Problems (hybrid meeting)
The classical theory of electromagnetism describes the interaction of electrically charged particles through electromagnetic forces, which are carried by the electric and magnetic fields. The propagation of the electromagnetic fields can be described by Maxwell's equations. Solving Maxwell's equations numerically is a challenging problem which appears in many different technical applications. Difficulties arise for instance from material interfaces or if the geometrical features are much larger than or much smaller than a typical wavelength. The spatial discretization needs to combine good geometrical flexibility with a relatively high order of accuracy.
The aim of this small-scale, week-long interactive mini-workshop jointly organized by the University of Duisburg-Essen and the University of Twente, and kindly hosted at the MFO, is to bring together experts in non-standard and mixed finite elements methods with experts in the field of electromagnetism
High-Resolution Mathematical and Numerical Analysis of Involution-Constrained PDEs
Partial differential equations constrained by involutions provide the highest fidelity mathematical models for a large number of complex physical systems of fundamental interest in critical scientific and technological disciplines. The applications described by these models include electromagnetics, continuum dynamics of solid media, and general relativity. This workshop brought together pure and applied mathematicians to discuss current research that cuts across these various disciplines’ boundaries. The presented material illuminated fundamental issues as well as evolving theoretical and algorithmic approaches for PDEs with involutions. The scope of the material covered was broad, and the discussions conducted during the workshop were lively and far-reaching
- …