Bibliographie & annexes
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’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 ».
L'expertise technique et scientifique de référence
DOI (Digital Object Identifier)
Cet article fait partie de l’offre
Mathématiques
(202 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
Intégrateurs numériques d'ordre supérieur1. Reconnaître une équation différentielle raide
Les équations différentielles raides, écrites sous forme autonome, sont des systèmes d'équations différentielles :
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 :
où 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
![]()
pour h → 0. L'erreur accumulée sur un intervalle fixe de longueur t, c.-à-d., après
![]()
pas de calcul, se comporte donc comme
![]()
, et nous disons que la méthode est d'ordre 1.
pour h → 0. L'erreur accumulée sur un intervalle fixe de longueur t, c.-à-d., après
pas de calcul, se comporte donc comme
, et nous disons que la méthode est d'ordre 1.-
Méthode du trapèze
Une...
L'expertise technique et scientifique de référence
Cet article fait partie de l’offre
Mathématiques
(202 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
Reconnaître une équation différentielle raide
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...
L'expertise technique et scientifique de référence
Cet article fait partie de l’offre
Mathématiques
(202 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
