Présentation
En anglaisRÉ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’articleABSTRACT
This article addresses the numerical integration of stiff differential equations and explains the phenomena that appear by using examples resulting from chemical reactions and discretized partial differential equations in space. A method for recognizing a stiff differential equation (standard methods, numerical phenomena, etc.) is proposed. Higher order numerical integrators are then introduced: implicit Runge-Kutta methods, multistep methods, extrapolation methods, and numerical integrator programming. In conclusion, methods for specific problems are proposed.
Auteur(s)
-
Gerhard WANNER : Professeurs - Département de mathématiques - Université de Genève
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 Intégration numérique des équations différentielles raides [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 ».
DOI (Digital Object Identifier)
Cet article fait partie de l’offre
Mathématiques
(166 articles en ce moment)
Cette offre vous donne accès à :
Une base complète d’articles
Actualisée et enrichie d’articles validés par nos comités scientifiques
Des services
Un ensemble d'outils exclusifs en complément des ressources
Un Parcours Pratique
Opérationnel et didactique, pour garantir l'acquisition des compétences transverses
Doc & Quiz
Des articles interactifs avec des quiz, pour une lecture constructive
Présentation
2. Intégrateurs numériques d'ordre supérieur
Les méthodes numériques comme la règle du trapèze, la règle du point-milieu et, surtout, la méthode d'Euler implicite (cf. § 1.4) sont certes adaptées à la résolution des problèmes raides, mais leur grand défaut est leur manque de précision. En pratique, on a besoin de méthodes d'ordre supérieur qui, toutefois, restent efficaces pour des problèmes raides. Nous allons discuter des méthodes de Runge-Kutta implicites, des BDF (méthodes multipas) ainsi que des méthodes d'extrapolation. Plus de détails et d'autres classes d'intégrateurs peuvent être trouvés dans [5].
2.1 Méthodes de Runge-Kutta implicites
Parmi les méthodes à un pas, les méthodes de Runge-Kutta implicites sont les plus importantes. Nous présentons deux cas particuliers.
HAUT DE PAGE
Considérons une équation différentielle et supposons connaître une approximation y n de la solution y (t n). Pour des nombres donnés, nous cherchons un polynôme u(t) de degré s qui vérifie :
L'approximation...
Cet article fait partie de l’offre
Mathématiques
(166 articles en ce moment)
Cette offre vous donne accès à :
Une base complète d’articles
Actualisée et enrichie d’articles validés par nos comités scientifiques
Des services
Un ensemble d'outils exclusifs en complément des ressources
Un Parcours Pratique
Opérationnel et didactique, pour garantir l'acquisition des compétences transverses
Doc & Quiz
Des articles interactifs avec des quiz, pour une lecture constructive
Intégrateurs numériques d'ordre supérieur
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...
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.
n 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).
n 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.
n https://www.l3harrisgeospatial.com/docs/lsode.html
LSODE est une programmation des méthodes multipas BDF.
n 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.
n http://www.ma.ic.ac.uk/~jcash/IVP_software/readme.html
MEBDF est synonyme de « Modified...
Cet article fait partie de l’offre
Mathématiques
(166 articles en ce moment)
Cette offre vous donne accès à :
Une base complète d’articles
Actualisée et enrichie d’articles validés par nos comités scientifiques
Des services
Un ensemble d'outils exclusifs en complément des ressources
Un Parcours Pratique
Opérationnel et didactique, pour garantir l'acquisition des compétences transverses
Doc & Quiz
Des articles interactifs avec des quiz, pour une lecture constructive