Overall Objectives
Scientific Foundations
Application Domains
New Results
Other Grants and Activities

Section: Scientific Foundations

Structure-preserving numerical schemes for solving ordinary differential equations

Participants : François Castella, Philippe Chartier, Erwan Faou, Gilles Vilmart.

In many physical situations, the time-evolution of certain quantities may be written as a Cauchy problem for a differential equation of the form

Im1 $\mtable{...}$(1)

For a given y0 , the solution y(t) at time t is denoted Im2 ${\#981 _t{(y_0)}}$ . For fixed t , Im3 $\#981 _t$ becomes a function of y0 called the flow of (1 ). From this point of view, a numerical scheme with step size h for solving (1 ) may be regarded as an approximation $ \upper_phi$h of Im4 $\#981 _h$ . One of the main questions of geometric integration is whether intrinsic properties of Im3 $\#981 _t$ may be passed on to $ \upper_phi$h .

This question can be more specifically addressed in the following situations:

Reversible ODEs

The system (1 ) is said to be $ \rho$ -reversible if there exists an involutive linear map $ \rho$ such that

Im5 $\mtable{...}$(2)

It is then natural to require that $ \upper_phi$h satisfies the same relation. If this is so, $ \upper_phi$h is said to be symmetric . Symmetric methods for reversible systems of ODEs are just as much important as symplectic methods for Hamiltonian systems and offer an interesting alternative to symplectic methods.

ODEs with an invariant manifold

The system (1 ) is said to have an invariant manifold g whenever

Im6 $\mtable{...}$(3)

is kept globally invariant by Im3 $\#981 _t$ . In terms of derivatives and for sufficiently differentiable functions f and g , this means that

Im7 $\mtable{...}$

As an example, we mention Lie-group equations, for which the manifold has an additional group structure. This could possibly be exploited for the space-discretisation. Numerical methods amenable to this sort of problems have been reviewed in a recent paper [43] and divided into two classes, according to whether they use g explicitly or through a projection step. In both cases, the numerical solution is forced to live on the manifold at the expense of some Newton's iterations.

Hamiltonian systems

Hamiltonian problems are ordinary differential equations of the form:

Im8 $\mtable{...}$(4)

with some prescribed initial values (p(0), q(0)) = (p0, q0) and for some scalar function H , called the Hamiltonian. In this situation, H is an invariant of the problem. The evolution equation (4 ) can thus be regarded as a differential equation on the manifold

Im9 $\mtable{...}$

Besides the Hamiltonian function, there might exist other invariants for such systems: when there exist d invariants in involution, the system (4 ) is said to be integrable . Consider now the parallelogram P originating from the point Im10 ${{(p,q)}\#8712 \#8477 ^{2d}}$ and spanned by the two vectors Im11 ${\#958 \#8712 \#8477 ^{2d}}$ and Im12 ${\#951 \#8712 \#8477 ^{2d}}$ , and let $ \omega$($ \xi$, $ \eta$) be the sum of the oriented areas of the projections over the planes (pi, qi) of P ,

Im13 $\mtable{...}$

where J is the canonical symplectic matrix

Im14 $\mtable{...}$

A continuously differentiable map g from Im15 $\#8477 ^{2d}$ to itself is called symplectic if it preserves $ \omega$ , i.e. if

Im16 $\mtable{...}$

A fundamental property of Hamiltonian systems is that their exact flow is symplectic. Integrable Hamiltonian systems behave in a very remarkable way: as a matter of fact, their invariants persist under small perturbations, as shown in the celebrated theory of Kolmogorov, Arnold and Moser. This behavior motivates the introduction of symplectic numerical flows that share most of the properties of the exact flow. For practical simulations of Hamiltonian systems, symplectic methods possess an important advantage: the error-growth as a function of time is indeed linear, whereas it would typically be quadratic for non-symplectic methods.

Differential-algebraic equations

Whenever the number of differential equations is insufficient to determine the solution of the system, it may become necessary to solve the differential part and the constraint part altogether. Systems of this sort are called differential-algebraic systems. They can be classified according to their index, yet for the purpose of this expository section, it is enough to present the so-called index-2 systems

Im17 $\mtable{...}$(5)

where initial values (y(0), z(0)) = (y0, z0) are given and assumed to be consistent with the constraint manifold. By constraint manifold, we imply the intersection of the manifold

Im18 $\mtable{...}$

and of the so-called hidden manifold

Im19 $\mtable{...}$

This manifold Im20 ${\#8499 =\#8499 _1\#8898 \#8499 _2}$ is the manifold on which the exact solution (y(t), z(t)) of (5 ) lives.

There exists a whole set of schemes which provide a numerical approximation lying on Im21 $\#8499 _1$ . Furthermore, this solution can be projected on the manifold Im22 $\#8499 $ by standard projection techniques. However, it it worth mentioning that a projection destroys the symmetry of the underlying scheme, so that the construction of a symmetric numerical scheme preserving Im22 $\#8499 $ requires a more sophisticated approach.


Logo Inria