Flot (mathématiques)

Le flot, coulée ou encore courant est, en mathématiques, un concept fondamental utilisé en géométrie différentielle.

La notion de flot permet notamment de modéliser le déplacement dans le temps des éléments d'un fluide. Pour ce faire, on crée une application α qui, à chaque point x de l'espace concerné par l'écoulement, associe un autre point α(x,t), correspondant à la position qu'aurait une particule du fluide à l'instant t, si elle avait été située en x à l'instant 0.

Flot associé à l'équation différentielle d'un pendule. Les abscisses représentent la position et les ordonnées la vitesse.

Une courbe intégrale d'un champ de vecteurs.

Le flot est associé à la notion de champ de vecteurs, c'est-à-dire à une application f, qui, à un point x d'un ouvert Ω d'un espace de Banach E, associe un vecteur de E. Un tel champ définit une équation différentielle du type α'(t) = f(α(t)), (c'est-à-dire, dans l'exemple du fluide, que la fonction f associe au point x à l'instant t, ses direction et vitesse de déplacement). Si la fonction f est localement lipschitzienne, pour chaque point x de Ω, il existe une solution maximale αModèle:Ind du problème de Cauchy constitué de cette équation différentielle et de la condition dite de Cauchy αModèle:Ind(0) = x. Vue comme une fonction de deux variables, t et x, l'application α est appelée le flot du champ f de vecteurs. Cette définition se généralise dans le cas d'un champ de vecteurs temporel (c'est-à-dire dépendant d'une variable t qui prend ses valeurs dans R) et dépendant d'un paramètre λ. Le flot et le champ de vecteurs deviennent des fonctions de trois variables : t, x et λ.

Si le champ de vecteurs f est régulier, le flot est le support de plusieurs théorèmes, piliers de la théorie des équations différentielles. Si la fonction f est de classe CModèle:Exp, le flot l'est aussi. Ce résultat est parfois considéré comme une forme avancée du théorème de Cauchy-Lipschitz. Si la fonction f ne dépend pas du temps, le théorème du redressement du flot indique que, localement, le champ de vecteurs est équivalent à un champ constant et cette équivalence transforme le flot en une fonction qui à (x, t) associe x + tv, où v est l'unique image du champ constant.

Le flot est utilisé dans diverses branches des mathématiques. En analyse qualitative des équations différentielles, il est le cadre d'expression de théorèmes, comme celui de Poincaré-Bendixson. On trouve la notion de flot de manière générale dans l'étude d'un système dynamique continu. En topologie algébrique, il est utilisé pour démontrer le théorème de la boule chevelue ou encore celui du point fixe de Brouwer ; des applications plus avancées définissent des flots caractéristiques de la géométrie des objets étudiés, tels que le flot de Ricci, outil de base utilisé par Grigori Perelman pour démontrer la conjecture de Poincaré. L'usage du flot dépasse le cadre strict des mathématiques ; ainsi, le flot de Ricci est à l'origine d'un des modes d'expression des équations de la relativité générale en physique.

Définitions

Introduction

Dans la suite de l'article, E désigne un espace de Banach. Le symbole (1) désigne l'équation différentielle x' = f(x), où f est fonction définie sur un ouvert Ω de E et à valeurs dans E et (2) l'équation x' = f(t, x). Dans le cas de l'équation (2), Ω désigne un ouvert de R×E et f est toujours une fonction définie sur Ω et à valeurs dans E. La condition de Cauchy appelée C dans l'article, est celle que vérifie une solution s de l'équation. Dans le cas de (1), elle signifie que s(0) = x₀, dans le cas (2) les notations choisies sont : s(t₀) = x₀. Il existe un troisième cas, plus général et qui permet en particulier l'étude des singularités de certaines équations différentielles. L'équation (3) est : x' = f(t, λ, x). Dans ce cas, Ω est un ouvert de R×F×E, où F est encore un espace de Banach. Dans les trois cas, la fonction f est supposée au moins continue, et localement lipschitzienne par rapport à x.

Le flot permet de formaliser un vocabulaire adapté pour l'étude de deux questions concernant les équations différentielles : la sensibilité à la condition initiale et le comportement asymptotique. Ces questions sont au cœur d'une branche mathématique appelée système dynamique.

Si la condition initiale, dite condition de Cauchy, est un peu modifiée, il se peut que la solution de l'équation soit de plus en plus éloignée de la courbe intégrale originelle, si la variable t augmente. Tel est le cas, par exemple, pour les systèmes chaotiques. L'étude de la régularité du flot apporte des premiers éléments de réponse.

Dans certaines conditions et pour une condition initiale donnée, la solution se prolonge à l'infini et « se stabilise » plus ou moins. Elle peut converger vers une valeur, s'approcher de plus en plus d'un comportement périodique, diverger ou encore adopter un comportement différent et appelé chaotique. L'étude de ces différents comportements est l'objet de la seconde question.

Sensibilité à la condition initiale

Pour l'étude de la sensibilité initiale, on étudie les flots locaux ou globaux, qui permettent d'énoncer une version plus forte du théorème de Cauchy-Lipschitz.

L'étude de la sensibilité à la condition initiale impose un vocabulaire et une représentation géométrique un peu différents de ceux utilisés dans les approches plus élémentaires. Pour comprendre l'origine de cette différence, le plus simple est de considérer le cas d'une équation différentielle autonome, c'est-à-dire de type (1), d'imaginer que Ω est un plan d'eau et que R représente le temps. le plan d'eau est agité par un courant, représenté par la fonction f, appelée champ de vecteurs. En dimension 2, on représente ce champ de vecteurs en associant à certains points x de Ω une représentation graphique du vecteur f(x), à l'image de la figure de droite. Une courbe intégrale satisfaisant à la condition de Cauchy C peut s'imaginer comme la trajectoire d'un bouchon placé dans l'eau à l'instant initial 0 et à la position x₀. Pour connaître d'un seul coup toutes les solutions de l'équation différentielle, il suffit de connaître le mouvement de la surface de l'eau, appelé flot, coulée ou encore courant. Dans le cas général de l'équation (2), c'est-à-dire celui où l'équation n'est pas nécessairement autonome, on dispose de la définition suivante :

Soit U un voisinage de x₀ tel que {t₀}×U soit inclus dans Ω et J un intervalle ouvert contenant t₀. Un flot local est la donnée d'une application β_t0 de J×U dans E tel que pour tout x de U, l'application de J dans E, qui à t associe β_t0(t, x) soit une courbe intégrale de condition de Cauchy (t₀, x) (c'est-à-dire que β_t0(t₀, x) = x).

Le théorème de Cauchy-Lipschitz garantit l'existence de courbes intégrales maximales. C'est-à-dire que, pour une condition de Cauchy donnée, il existe un intervalle maximal, une unique solution de l'équation définie sur cet intervalle et vérifiant la condition de Cauchy. Ce résultat permet d'enrichir les définitions associées au flot.

Soit x un point de E et JModèle:Ind l'intervalle égal au domaine de définition de la solution maximale satisfaisant la condition de Cauchy (t₀, x), qui peut être vide si, par exemple, (t₀, x) n'est pas élément de Ω. On note D(f) et Modèle:Refnec.

Modèle:Refconf de D(f) dans E, telle que l'application qui à t associe α_t0(t, x) soit la solution maximale de (2) associée à la condition de Cauchy (t₀, x).

Dans toute la suite de l'article α désigne le flot global du champ vectoriel f. Si α ne comporte pas d'indice, cela signifie que la fonction qui à t associe α(t, x) est la courbe intégrale maximale s vérifiant la condition de Cauchy : s(0) = x. Sinon, α_t0 désigne le flot global vérifiant α_t0(t₀, x) = x.

Comportement asymptotique

La deuxième question concerne le comportement asymptotique du flot, autrement dit ce qui se passe une fois le système stabilisé, si jamais il se stabilise. Pour la traiter, un vocabulaire spécifique existe ; il suppose en général que l'équation différentielle considérée est autonome, c'est-à-dire de type (1) Modèle:Supra.

L'orbite ou la trajectoire d'un point x de Ω est l'image d'une courbe intégrale maximale passant par x.

Le théorème de Cauchy-Lipschitz montre que les orbites forment une partition de Ω.

Un ensemble invariant par le flot est une réunion d'orbites.

Pour reprendre la métaphore de l'introduction, un ensemble A est invariant par le flot si un bouchon qui est dans A y a toujours été et y restera toujours. Modèle:Citation : on peut limiter l'étude du comportement asymptotique à une telle zone.

Un point de x de Ω est dit point d'équilibre, point fixe ou point stationnaire si son orbite est réduite au singleton {x} ou encore si x est un zéro de la fonction f. S'il existe un voisinage de x tel que toute courbe intégrale maximale qui rencontre ce voisinage converge vers x, le point est dit attractif.

Un point d'équilibre xModèle:Ind est un point sans courant, ou encore un point tel que, si l'on y place le bouchon, il reste indéfiniment immobile. Deux comportements différents peuvent se produire, le point peut être attractif ou non. Il peut exister une zone pas trop petite (d'intérieur non vide) tel que tout point dans cette zone finit par s'immobiliser en xModèle:Ind. Ce point semble « attirer » les trajectoires vers lui, pour cette raison, on parle de point attractif. Sinon, tout point se rapprochant de xModèle:Ind finit par s'éloigner, ce point d'équilibre est alors instable et il faut positionner le bouchon exactement sur le point pour qu'il y reste immobile. Dans le cas général, d'autres comportements que la convergence sont possibles :

Soit s une courbe intégrale définie sur R.

• Lorsque s admet une limite en Modèle:Math et en Modèle:Math, l'orbite est dite hétérocline si les deux limites sont distinctes et homocline sinon.
• Si s est périodique, l'orbite est dite périodique.
• L'ensemble ω-limite de l'orbite est l'ensemble des valeurs d'adhérence de la fonction s en Modèle:Math. On le note ω(x), où x est un point quelconque de l'orbite. On définit de même l'ensemble α-limite α(x) comme l'ensemble des valeurs d'adhérence de s en Modèle:Math.

Ces deux ensembles sont invariants par le flot. Si l'orbite est périodique, ses deux ensembles limites lui sont égaux.

Une notion connexe est celle d'attracteur. L'attracteur futur est le plus petit ensemble contenant tous les ensembles ω(x) si x décrit Ω, à l'exception, peut-être d'un ensemble de mesure nulle. L'attracteur passé correspond à la même définition, mais cette fois-ci avec les ensembles α-limites.

Sensibilité à la condition initiale

Résumé

Modèle:Article détaillé Dans sa version élémentaire, le théorème de Cauchy-Lipschitz montre l'existence et l'unicité d'une courbe intégrale maximale. On peut alors définir le flot global α comme la fonction qui, à un triplet (tModèle:Ind, tModèle:Ind, xModèle:Ind) avec (tModèle:Ind, xModèle:Ind) élément de Ω, associe (lorsqu'elle est définie) la valeur en tModèle:Ind de la solution maximale du problème de Cauchy xModèle:' = f(t, x), x(tModèle:Ind) = xModèle:Ind. On sait déjà que par rapport à tModèle:Ind, α est de classe CModèle:Exp si f est de classe CModèle:Exp, mais on va préciser son comportement par rapport au triplet (tModèle:Ind, tModèle:Ind, xModèle:Ind) :

• Le flot global est défini sur un ouvert et localement lipschitzien.
• Si la fonction f est de classe CModèle:Exp, il en est de même du flot global.

Pour certaines études, comme celles analysant des singularités de courbes intégrales, il est utile d'ajouter un paramètre à l'équation, qui prend la forme : x' = f(t, x, λ). Le flot global dépend alors du paramètre λ, choisi dans un Banach F. La question est celle de la régularité du flot global en fonction du paramètre. Le théorème décrivant cette situation est parfois dénommé Modèle:Citation. Si la fonction f est de classe CModèle:Exp, le flot, considéré comme une application sur un ouvert de R×R×E×F, est aussi de classe CModèle:Exp.

Continuité du flot

Dans un premier temps, on étudie le cas de l'équation (2), c'est-à-dire celui d'une équation non nécessairement autonome, mais sans paramètre. L'objectif est d'établir la continuité locale d'un flot, c'est-à-dire d'analyser ce qui se passe si la condition de Cauchy est peu modifiée.

On commence par préciser le comportement de α(t, tModèle:Ind, xModèle:Ind) par rapport à (tModèle:Ind, xModèle:Ind) dans le cas où t est proche de tModèle:Ind.

En affinant la technique de point fixe utilisée pour démontrer le théorème de Cauchy-Lipschitz classique, on parvient au résultat suivant :

Un flot discontinu ne peut pas se produire si le champ de vecteurs le définissant est continu et défini sur un ouvert.

Pour rendre global le résultat précédent de continuité « au voisinage de la diagonale », il reste encore un peu de travail. Il s'agit surtout de montrer que la configuration de droite ne peut pas se produire. En reprenant la métaphore de l'introduction, on suppose que l'étang contient un rocher. Sur la zone rouge, le flot contourne le rocher, sur la zone bleue il passe tout droit. Le flot possède alors des points de discontinuité. Pour s'en rendre compte, on considère l'évolution d'un voisinage d'un point symbolisé par une zone grise circulaire à un instant 0, sur la figure de droite. Ce point est situé à la frontière des deux zones. À l'instant t₀, ce voisinage est découpé en deux parties éloignées l'une de l'autre. Cette configuration peut se produire, par exemple si Ω n'est pas un ouvert. En revanche, dans les hypothèses étudiées, cette configuration est impossible :

Une remarque est utile pour le démontrer. Elle est formulée ici uniquement pour les équations autonomes. Le point β(s, x) désigne la position, à l'instant s, du point qui était en x à l'instant 0. Le point β(t, β(s, x)) désigne la position à l'instant t, du point qui était en β(s, x), à l'instant 0. Il désigne aussi la position à l'instant s + t du point qui était en position x à l'instant s. Ce qui signifie, si les différentes valeurs sont dans les bons domaines de définitions :

$\beta (t,\beta (s,x))=\beta (t+s,x)$

Ainsi, si les domaines de définitions s'y prêtent, il est possible de combiner les flots locaux, pour en construire un plus vaste. Cette technique permet de munir les flots locaux d'une structure de semi-groupe. Si le flot global est toujours défini sur R, les solutions intégrales maximales partageant une portion de domaine de définition forment un groupe, image de R par un morphisme de groupes.

Équation différentielle avec paramètre

On souhaite maintenant étudier le cas d'un flot local associé à l'équation :

$(3)\frac{\mathrm{d}x}{\mathrm{d}t}=f(t,x,\lambda )$.

On suppose que λ est élément d'un espace de Banach F, Ω est un ouvert de R×E×F et f est une fonction continue de Ω dans E, localement lipschitzienne par rapport à (λ, x).

Cette situation, en apparence plus générale que celle de l'équation (2), en est en fait un cas particulier, en remplaçant simplement, dans cette dernière :

• l'espace de Banach E (dans lequel chaque courbe intégrale prend ses valeurs) par E×F ;
• le champ de vecteurs (d'un ouvert de R×E dans E) par la fonction :$g:\mathrm{\Omega }\to E\times F,y=(t,x,\lambda )\mapsto (f(y),0)$ ;
• l'équation xModèle:' = f(t, x) par (xModèle:', λModèle:') = g(t, x, λ) ;
• la « condition initiale » (tModèle:Ind, xModèle:Ind) par (tModèle:Ind, xModèle:Ind, λModèle:Ind).

Le flot α dépend maintenant de trois variables : t, tModèle:Ind, xModèle:Ind et la « variable » supplémentaire λ, mais compte tenu de la définition de g, la valeur de λ le long d'une courbe intégrale reste constamment égale à sa valeur initiale λModèle:Ind. Cet artifice permet d'étendre directement aux équations à paramètres le résultat du paragraphe précédent : soit f une fonction continue, d'un ouvert de R×F×E dans E, localement lipschitzienne par rapport à ses variables dans F×E. Le flot global associé à l'équation (3) est défini sur un ouvert de R×R×F×E et localement lipschitzien.

Mais ce corollaire n'est d'aucune utilité pour la suite, contrairement au théorème suivant, qui généralise le lemme du paragraphe précédent et se démontre (directement) de la même manière : Modèle:Théorème

Le jeu d'écriture lui-même sera utile en revanche, pour étendre aux équations de type (3) un lemme (ci-dessous) sur la régularité CModèle:1.

Régularité du flot au voisinage de la diagonale

Les démonstrations associées à la régularité du flot ont pendant longtemps été un peu complexes. De manière indépendante, Modèle:Refsou et RobbinModèle:,Modèle:, ont trouvé une démonstration élémentaire utilisant le théorème des fonctions implicites. Cette fois, l'équation considérée est celle notée (1), c'est-à-dire qu'on ne traite que d'une équation différentielle autonome. On en déduit une partie du résultat analogue pour les équations non autonomes.

Soient x₀ un point de Ω et b un réel strictement positif tel que la courbe intégrale βModèle:Ind qui vaut xModèle:Ind en 0 soit définie au moins sur [–b, b]. Choisissons un ouvert Ω' contenant xModèle:Ind et un ouvert Ω" contenant Im(βModèle:Ind) – xModèle:Ind, tels que Ω' + Ω" soit inclus dans Ω. Notons enfin G l'espace de Banach des fonctions de classe CModèle:1 de [–b, b] dans E et nulles en 0, muni de la norme de la convergence uniforme de la dérivée, V l'ouvert de celles à valeurs dans Ω", et H l'espace de Banach des fonctions continues de [–b, b] dans E, muni de la norme de la convergence uniforme.

On considère l'application T de Ω'×V dans H définie par :

$\mathrm{\forall }x\in {\mathrm{\Omega }}^{\prime }\mathrm{\forall }\sigma \in V\mathrm{\forall }t\in [-b,b]T(x,\sigma )(t)={\sigma }^{\prime }(t)-f(x+\sigma (t))$.

Dire qu'un couple (x, σ) est un zéro de la fonction T revient à dire que x + σ est une solution sur [–b, b] de l'équation (1) avec la condition de Cauchy x + σ(0) = x. En particulier, (xModèle:Ind, βModèle:Ind – xModèle:Ind) est un zéro de T.

On démontre que si f est de classe CModèle:1 alors T aussi et pour b suffisamment petit, T satisfait les hypothèses du théorème des fonctions implicites au point (xModèle:Ind, βModèle:Ind – xModèle:Ind). Cela signifie que la fonction implicite donnée par l'équation T(x, σ) = 0 est de classe CModèle:1 au voisinage de xModèle:Ind. Or cette fonction est celle qui associe à x la fonction β(∙, x) – x définie sur [–b, b], ce qui permet de démontrer un premier résultat :

Ce lemme s'étend aux équations différentielles de type (3), par le jeu d'écriture ci-dessus pour le paramètre et un autre, analogue, pour la variable temporelle.

Or dès que f est de classe CModèle:1, la deuxième différentielle partielle de β vérifie une équation différentielle de type (3) :

$D_1{D}_2\beta (t,x)=D_{2}f(t,\beta (t,x))\circ D_2\beta (t,x)$.

Ceci permet de montrer par récurrence que pour toute équation de type (3), si f est de classe CModèle:Exp alors le flot α(t, tModèle:Ind, λ, xModèle:Ind) l'est aussi, par rapport à (t, λ, xModèle:Ind). Nous admettrons la régularité par rapport à tModèle:Ind, qui utilise entre autres, pour l'initialisation de la récurrence, le théorème de continuité ci-dessus. Moyennant quoi, on peut énoncer :

Régularité globale du flot

Le théorème de Cauchy-Lipschitz, sous sa forme élémentaire, garantit à la fois l'existence et l'unicité du flot global α, et sa régularité par rapport à la première variable. Le paragraphe précédent apporte aussi des informations. Il établit que si t est proche de 0 et si t et x varient peu, le flot est de classe CModèle:Exp, si f l'est. En revanche, la régularité globale du flot α suppose aussi que l'expression α(t, x) soit de classe CModèle:Exp en x, même si t n'est pas proche de 0. On a supposé ici que (t, x) est un point du domaine de définition de α.

La même technique que dans le paragraphe « Continuité du flot » Modèle:Supra permet de rendre global le résultat précédent de régularité du flot « au voisinage de la diagonale » pour une équation de type (2) — donc aussi, par le « jeu d'écriture » Modèle:Supra, pour une équation de type (3) :

• Théorème : Pour une équation différentielle non autonome à paramètre, si le champ de vecteurs f est de classe CModèle:Exp, le flot global l'est aussi.

Comportement asymptotique

Comportement à la limite du domaine

L'équation différentielle x' = –xModèle:2 n'admet pas de solution sur R tout entier si l'image de 0 n'est pas nulle.

Pour une valeur de x prise dans E, la courbe intégrale maximale t ↦ s(t) = α_t0(t, x) est définie sur un intervalle ouvert. Cet intervalle n'est pas toujours égal à R. Deux « obstacles » peuvent se présenter. Un premier cas est celui où la courbe a pour extrémité un point frontière de l'ouvert Ω. Si tel n'est pas le cas, il n'existe qu'un seul obstacle possible : la dérivée de s n'est pas bornée. Cela résulte de la proposition suivante :

• Soit b la borne supérieure du domaine de définition d'une courbe intégrale maximale s. Si b est fini et si sModèle:' est bornée au voisinage de b, alors s admet au point b une limite c et le couple (b, c) appartient à la frontière de Ω.

Un exempleModèle:, est donné par l'équation x' = –xModèle:2. La solution s'écrit :

$\alpha (t,x)=\frac{x}{xt+1}$.

Les courbes intégrales dont l'ordonnée à l'origine est différente de 0 ne se prolongent pas au-delà de –1/x.

Ensemble ω-limite

Le système dynamique de Lorenz montre un exemple où les ensembles ω-limites sont en général de structure complexe.

On suppose ici que l'équation différentielle étudiée est autonome. Une méthode pour étudier le comportement asymptotique est d'analyser l'ensemble ω-limite d'une orbite. Cet ensemble possède toujours les propriétés suivantes :

• Soit x un élément de Ω, l'ensemble ω-limite ω(x) est un fermé invariant par le flot.

En particulier, si y est un point de la trajectoire passant par x, alors ω(x) = ω(y). Si l'orbite de x n'est pas à valeurs dans un ensemble borné, la proposition peut s'avérer moins forte qu'il n'y parait : l'ensemble ω(x) peut être vide. Tel est le cas si la courbe intégrale s ayant pour condition de Cauchy s(0) = x est définie par s(t) = tv où v est un vecteur non nul de E. Si l'orbite est compacte, on a par contre :

• Avec les notations précédentes, si la trajectoire de x est à valeurs dans un compact K, ω(x) est un connexe compact non vide.

La trajectoire peut néanmoins être complexe, le système dynamique de Lorenz est un exemple exhibant des ensembles ω-limites compacts dont la dimension de Hausdorff n'est pas égale à 1.

Théorème de Poincaré-Bendixson

Modèle:Article détaillé Il existe un cas où le flot dispose d'un comportement asymptotique relativement simple, il se produit si l'équation différentielle associée est autonome, si l'espace E est le plan et si la courbe intégrale considérée est compacte. On dispose alors du théorème suivant, si α est le flot et x un point de Ω.

• Dans le cadre des hypothèses du paragraphe, la fonction qui à t associe α(t, x) est définie sur R. Si cette fonction n'est pas convergente, l'ensemble ω-limite est l'image d'une orbite cyclique.

Ce théorème se généralise mal. Si E est de dimension strictement supérieure à 2, le système dynamique de Lorenz montre que le résultat du théorème n'est plus généralement vrai.

Notes et références

Ouvrages cités

• Modèle:BergerGostiaux1 Modèle:Commentaire biblio
• Modèle:Article
• Modèle:Ouvrage
• Paul Malliavin, Géométrie différentielle intrinsèque, Hermann, 1972 Modèle:ISBN Modèle:Commentaire biblio
• Modèle:Ouvrage Modèle:Commentaire biblio
• Modèle:Lien web
• Modèle:Lien web, publié dans Géométries et dynamiques, Hermann, 2008 Modèle:ISBN, Modèle:P.