Contactez-nous
Reconnaître une équation différentielle raide
Intégration numérique des équations différentielles raides
AF653 v1 Article de référence

Reconnaître une équation différentielle raide
Intégration numérique des équations différentielles raides

Auteur(s) : Ernst HAIRER, Gerhard WANNER

Date de publication : 10 oct. 2007 | Read in English

Logo Techniques de l'Ingenieur Cet article est réservé aux abonnés
Pour explorer cet article plus en profondeur Consulter l'extrait gratuit

Déjà abonné ?

Présentation

1 - Reconnaître une équation différentielle raide

2 - Intégrateurs numériques d'ordre supérieur

3 - Méthodes pour des problèmes particuliers

Sommaire

Présentation

RÉSUMÉ

Cet article traite de l’intégration numérique des équations différentielles raides, et explique les phénomènes qui apparaissent, en s'appuyant sur des exemples issus des réactions chimiques et des équations aux dérivées partielles discrétisées en espace. Une méthode de reconnaissance d’une équation différentielle raide (méthodes standards, phénomène numériques, etc.) est proposée. Sont ensuite introduits les intégrateurs numériques d’ordre supérieur : méthodes de Runge-Kutta implicites, méthodes multipas, méthodes d’extrapolation, programmation des intégrateurs numériques. Pour terminer, des méthodes pour des problèmes particuliers sont proposés.

Lire cet article issu d'une ressource documentaire complète, actualisée et validée par des comités scientifiques.

Lire l’article

Auteur(s)

INTRODUCTION

Dans de nombreuses applications, la dynamique d'un système peut être modélisée par des équations différentielles. L'étude de systèmes mécaniques (par exemple en astronomie ou en dynamique moléculaire), l'analyse des circuits électriques ou la théorie du contrôle (robotique) nous fournissent de tels problèmes. Souvent, pour les problèmes dits raides, les méthodes standards ne fournissent pas une solution correcte en un temps de calcul acceptable.

Ce dossier récapitulatif explique les phénomènes qui apparaissent dans les équations différentielles raides, en s'appuyant sur des exemples issus des réactions chimiques ainsi que des équations aux dérivées partielles discrétisées en espace. Les propriétés essentielles des intégrateurs numériques pour la résolution des équations raides sont discutées (A-stabilité, domaine de stabilité). Pour des problèmes généraux, les méthodes de Runge-Kutta implicites, les méthodes multipas (BDF) et les méthodes d'extrapolation sont traitées. Pour des problèmes raides particuliers de grande dimension sont également abordées les méthodes explicites avec grande région de stabilité, les méthodes de séparation et les méthodes implicites-explicites. Une liste de programmes informatiques du domaine public est donnée en Documentation [Doc. AF 653].

Comme références sur la résolution numérique des équations différentielles raides, le lecteur pourra consulter les ouvrages généraux suivants [1] [2] [3] [4] [5] [6], mentionnées en « Pour en savoir plus ».

Logo Techniques de l'Ingenieur

Cet article est réservé aux abonnés.
Il vous reste 92 % à découvrir.

Pour explorer cet article Consulter l'extrait gratuit

Déjà abonné ?


DOI (Digital Object Identifier)

https://doi.org/10.51257/a-v1-af653

Article inclus dans l'offre

"Mathématiques"

(170 articles)

Une base complète d’articles

Actualisée et enrichie d’articles validés par nos comités scientifiques.

Des contenus enrichis

Quiz, médias, tableaux, formules, vidéos, etc.

Des modules pratiques

Opérationnels et didactiques, pour garantir l'acquisition des compétences transverses.

Des avantages inclus

Un ensemble de services exclusifs en complément des ressources.

Voir l'offre

1. Reconnaître une équation différentielle raide

Les équations différentielles raides, écrites sous forme autonome, sont des systèmes d'équations différentielles :

y˙=f(y)y(t0)=y0

où les méthodes standards (explicites) se révèlent inefficaces. Il est cependant utile de commencer par présenter ces méthodes, car elles permettent de comprendre les phénomènes et les difficultés qui interviennent dans la résolution des problèmes raides. Pour plus de détails, on pourra consulter le dossier [AF 1 220] des Techniques de l'Ingénieur, référence [42].

1.1 Méthodes standards pour équations différentielles

  • Méthode d'Euler

    L'algorithme le plus simple pour calculer des approximations numériques y n de la solution y(t) aux points t 1, t 2, t 3, … est donné par :

    yn +1 = yn + hf  (yn )

    h = t n +1 − t n est appelé longueur de pas.

    Cette formule représente le début de la série de Taylor pour la solution, et possède donc une erreur locale de O(h2) pour h → 0. L'erreur accumulée sur un intervalle fixe de longueur t, c.-à-d., après ...

Logo Techniques de l'Ingenieur

Cet article est réservé aux abonnés.
Il vous reste 92 % à découvrir.

Pour explorer cet article Consulter l'extrait gratuit

Déjà abonné ?


Lecture en cours
Reconnaître une équation différentielle raide

Article inclus dans l'offre

"Mathématiques"

(170 articles)

Une base complète d’articles

Actualisée et enrichie d’articles validés par nos comités scientifiques.

Des contenus enrichis

Quiz, médias, tableaux, formules, vidéos, etc.

Des modules pratiques

Opérationnels et didactiques, pour garantir l'acquisition des compétences transverses.

Des avantages inclus

Un ensemble de services exclusifs en complément des ressources.

Voir l'offre

Sommaire
Sommaire

BIBLIOGRAPHIE

  • (1) - WILLOUGHBY (R.A.), éd -   Stiff differential systems  -  . Plenum Press, New York, The IBM Research Symposia Series (1974).

  • (2) - MIRANKER (W.L.) -   Numerical methods for stiff equations and singular perturbation problems ; volume 5 of Mathematics and its Applications  -  . D. Reidel Publishing Co., Dordrecht (1981).

  • (3) - DEKKER (K.), VERWER (J.G.) -   Stability of Runge-Kutta methods for stiff nonlinear differential equations, volume 2 of CWI Monographs  -  . North-Holland Publishing Co., Amsterdam (1984).

  • (4) - AIKEN (R.C.), éd -   Stiff computation  -  . Oxford University Press, New York (1985).

  • (5) - HAIRER (E.), WANNER (G.) -   Solving Ordinary Differential Equations II. Stiff and Differential-Algebraic Problems  -  . Springer Series in Computational Mathematics 14. Springer-Verlag, Berlin, 2e édition (1996).

  • (6) - HUNDSDORFER (W.), VERWER...

1 Programmes informatiques

Pour la résolution des équations différentielles raides, il existe d'excellents programmes du domaine public qui peuvent être téléchargés sur internet.

http://www.unige.ch/~hairer/software.html

RADAU5 est basé sur une méthode de Runge-Kutta implicite (Radau IIA, § 2.1.1) d'ordre 5.

RADAU est une version avec ordre variable (ordres 5, 9, et 13).

SEULEX est la méthode d'extrapolation basée sur la méthode d'Euler linéairement implicite, discutée au paragraphe 2.3.

ROCK est une méthode de Runge-Kutta explicite avec un grand intervalle de stabilité (§ 3.1.2).

http://www.netlib.org/ode/index.html

DASSL est un programme célèbre de L. Petzold basé sur les méthodes BDF .

RKC est la méthode explicite discutée dans le paragraphe 3.1.1.

https://www.l3harrisgeospatial.com/docs/lsode.html

LSODE est une programmation des méthodes multipas BDF.

http://www.zib.de/Numerik/numsoft/CodeLib/ivpode.en.html

EULSIM est la méthode d'extrapolation basée sur la méthode d'Euler linéairement implicite, discutée au paragraphe 2.3.

LIMEX est une extension aux problèmes différentiels-algébriques.

http://www.ma.ic.ac.uk/~jcash/IVP_software/readme.html

MEBDF est synonyme de « Modified...

Logo Techniques de l'Ingenieur

Cet article est réservé aux abonnés.
Il vous reste 95 % à découvrir.

Pour explorer cet article Consulter l'extrait gratuit

Déjà abonné ?


Article inclus dans l'offre

"Mathématiques"

(170 articles)

Une base complète d’articles

Actualisée et enrichie d’articles validés par nos comités scientifiques.

Des contenus enrichis

Quiz, médias, tableaux, formules, vidéos, etc.

Des modules pratiques

Opérationnels et didactiques, pour garantir l'acquisition des compétences transverses.

Des avantages inclus

Un ensemble de services exclusifs en complément des ressources.

Voir l'offre