Description lagrangienne

En représentation lagrangienne, la position M à l'instant t de la particule qui se trouvait en ${\displaystyle M_{0}}$ à l'instant 0 est donnée par une relation du type

${\displaystyle M=f(M_{0},t)\,}$.

Cela correspond à la description paramétrique de la trajectoire en coordonnées cartésiennes :

${\displaystyle x=f_{x}(x_{0},y_{0},z_{0},t)\,}$

${\displaystyle y=f_{y}(x_{0},y_{0},z_{0},t)\,}$

${\displaystyle z=f_{z}(x_{0},y_{0},z_{0},t)\,}$.

Cette méthode présente un inconvénient : le référentiel se déplace avec le fluide. Il est donc difficile de connaître l'état du fluide en un point donné de l'espace et du temps.

La représentation d'Euler définit à tout instant la valeur d'une grandeur (par exemple une composante de la vitesse) associée à un point fixe de l'écoulement. La variation de cette grandeur au cours du temps est alors décrite par une dérivée partielle parfois appelée dérivée eulérienne.

La représentation de Lagrange suit une particule dans son mouvement. La variation précédente est alors représentée par la dérivée particulaire ou dérivée totale ou dérivée lagrangienne. Elle tient compte non seulement de la variation locale du paramètre au cours du temps mais aussi de la variation de celui-ci liée au déplacement de la particule.

Le lien entre la description lagrangienne et la description eulérienne est démontré ici pour un mouvement à une dimension. Pendant l'intervalle de temps ${\displaystyle \mathrm {d} t}$, une particule située en ${\displaystyle x}$ à l'instant ${\displaystyle t}$ s'est déplacée en ${\displaystyle x+v\mathrm {d} t}$. La variation de la grandeur ${\displaystyle f}$ s'écrit donc :

${\displaystyle f(x+v\mathrm {d} t,t+\mathrm {d} t)-f(x,t)=\left({\frac {\partial f}{\partial t}}+v{\frac {\partial f}{\partial x}}\right)\mathrm {d} t=\mathrm {d} f}$

Remarquons que le résultat obtenu peut être établi à l'aide d'un développement de Taylor bi-dimensionnel à l'ordre 1. Ici, ${\displaystyle \mathrm {d} f}$ est en fait une approximation du premier ordre de ce développement de Taylor. En divisant par ${\displaystyle \mathrm {d} t}$ on obtient la dérivée particulaire qui s'écrit avec ${\displaystyle {\frac {\mathrm {D} f}{\mathrm {D} t}}}$ comme notation la plus utilisée :

${\displaystyle {\frac {\mathrm {D} f}{\mathrm {D} t}}={\frac {\partial f}{\partial t}}+v{\frac {\partial f}{\partial x}}}$.

La formule se généralise à trois dimensions en introduisant le gradient de la grandeur ${\displaystyle f}$ :

${\displaystyle {\frac {\mathrm {D} f}{\mathrm {D} t}}={\frac {\partial f}{\partial t}}+{\overrightarrow {V}}\cdot {\overrightarrow {\mathrm {grad} }}\,f}$.

L'opérateur ${\displaystyle {\overrightarrow {V}}\cdot {\overrightarrow {\mathrm {grad} }}}$ est dénommé opérateur d'advection.

Soient les notations ${\displaystyle f_{L}}$ et ${\displaystyle f_{E}}$ désignant une même propriété d'un fluide, exprimée dans une description lagrangienne (${\displaystyle f_{L}}$ ) ou eulérienne (${\displaystyle f_{E}}$). Considérons la particule ${\displaystyle p}$, ayant pour coordonnée spatiale ${\displaystyle x}$. On peut noter :

${\displaystyle f_{L}(p,t)=f_{E}(x,t)}$

Dans le référentiel de la particule ${\displaystyle p}$, la propriété ${\displaystyle f_{L}(p,t)}$ dépend uniquement du temps. Ainsi, on peut contracter l'écriture pour obtenir une expression de ${\displaystyle f_{L}(p,t)}$ qui soit spécifique à la particule ${\displaystyle p}$ :

${\displaystyle f_{L}(p,t)=fp_{L}(t)}$

Au temps ${\displaystyle t+dt}$, la particule ${\displaystyle p}$ possède la coordonnée spatiale ${\displaystyle x+dx}$. L'équivalence entre la description lagrangienne et eulérienne s'écrit alors de la manière qui suit:

${\displaystyle fp_{L}(t+dt)=f_{E}(x+dx,t+dt)}$

Naturellement, l'expression suivante peut alors être établie :

${\displaystyle fp_{L}(t+dt)-fp_{L}(t)=f_{E}(x+dx,t+dt)-f_{E}(x,t)}$

Il est alors clair qu'un développement de Taylor permet d'établir le lien entre la dérivée particulaire d'ordre ${\displaystyle 1}$ de ${\displaystyle f_{E}}$ et la dérivée temporelle de ${\displaystyle fp_{L}}$. Pour le membre de gauche, un développement limité donne :

${\displaystyle fp_{L}(t+dt)-fp_{L}(t)=dt{\frac {dfp_{L}}{dt}}+{\frac {dt^{2}}{2!}}{\frac {d^{2}fp_{L}}{dt^{2}}}+...+o(dt^{n})}$

Pour le membre de droite, un développement limité bidimensionnel donne :

${\displaystyle f_{E}(x+dx,t+dt)-f_{E}(x,t)=dt{\frac {\partial f_{E}}{\partial t}}+dx{\frac {\partial f_{E}}{\partial x}}+dxdt{\frac {\partial ^{2}f_{E}}{\partial x\partial t}}+{\frac {dt^{2}}{2!}}{\frac {\partial ^{2}f_{E}}{\partial ^{2}t}}+{\frac {dx^{2}}{2!}}{\frac {\partial ^{2}f_{E}}{\partial ^{2}x}}+...+o(dt^{n},dx^{n})}$

En regroupant les termes de même ordre, il vient d'abord pour l'ordre ${\displaystyle 1}$:

${\displaystyle dt{\frac {dfp_{L}}{dt}}=dt{\frac {\partial f_{E}}{\partial t}}+dx{\frac {\partial f_{E}}{\partial x}}}$

Et évidemment, avec ${\displaystyle D_{t}}$ désignant la dérivée totale d'ordre ${\displaystyle 1}$ :

${\displaystyle {\frac {dfp_{L}}{dt}}={\frac {\partial f_{E}}{\partial t}}+{\frac {dx}{dt}}{\frac {\partial f_{E}}{\partial x}}=D_{t}f_{E}}$

Il est intéressant de remarquer que la dérivée totale de ${\displaystyle f_{E}}$ correspond en fait à la simple dérivée temporelle de ${\displaystyle fp_{L}}$. Cela peut en effet faciliter l'intégration d'équations aux dérivées partielles faisant intervenir la dérivée totale. Ainsi, la dérivée totale de ${\displaystyle f_{E}}$ peut être remplacée par la dérivée temporelle de ${\displaystyle fp_{L}}$. Une éventuelle intégration suppose bien sûr d'être réalisée en remplaçant ${\displaystyle f_{E}}$ par ${\displaystyle fp_{L}}$. Le résultat peut ensuite être reconverti en remplaçant, à l'inverse, ${\displaystyle fp_{L}}$ par ${\displaystyle f_{E}}$ dans la solution éventuelle.

Pour le second ordre, on peut introduire la notion de dérivée totale d'ordre ${\displaystyle 2}$, notée ${\displaystyle D_{t}^{2}}$ et en posant ${\displaystyle u={\frac {dx}{dt}}}$:

${\displaystyle {\frac {dt^{2}}{2}}{\frac {d^{2}fp_{L}}{dt^{2}}}=dxdt{\frac {\partial ^{2}f_{E}}{\partial x\partial t}}+{\frac {dt^{2}}{2}}{\frac {\partial ^{2}f_{E}}{\partial ^{2}t}}+{\frac {dx^{2}}{2}}{\frac {\partial ^{2}f_{E}}{\partial ^{2}x}}}$

Vient alors l'expression qui suit :

${\displaystyle D_{t}^{2}f_{E}={\frac {d^{2}fp_{L}}{dt^{2}}}={\frac {\partial ^{2}f_{E}}{\partial ^{2}t}}+u^{2}{\frac {\partial ^{2}f_{E}}{\partial ^{2}x}}+2u{\frac {\partial ^{2}f_{E}}{\partial x\partial t}}}$

Dans le cadre de cette description, et ${\displaystyle \rho }$ désignant la densité du fluide, ${\displaystyle \rho ({\overrightarrow {x}},t)}$ désigne la densité du fluide qui, initialement (temps 0), était à la position ${\displaystyle {\overrightarrow {x}}}$ et se trouve désormais (temps t) en ${\displaystyle {\overrightarrow {X}}({\overrightarrow {x}},t)}$... mais cela peut bien sûr s'appliquer à n'importe quelle autre fonction décrivant une propriété locale du fluide.

Cette description donne une bonne idée de ce qui se passe dans le fluide, par exemple, si ${\displaystyle {\frac {\mathrm {d} \rho }{\mathrm {d} t}}<0}$, alors on peut affirmer que le fluide s'étend (la densité de la particule fluide baisse). En particulier, on peut faire un bilan des forces s'appliquant à la particule fluide que l'on suit, et appliquer la relation fondamentale de la dynamique en écrivant que la somme des forces vaut la masse de la particule fluide multipliée par son accélération, écrite comme ${\displaystyle {\frac {\mathrm {d} {\overrightarrow {u}}}{\mathrm {d} t}}}$ où ${\displaystyle {\overrightarrow {u}}={\frac {\mathrm {d} {\overrightarrow {X}}}{\mathrm {d} t}}}$ est la vitesse.

• Dynamique des fluides
• Mécanique des milieux continus
• (en) Lagrangian description dans l'article Continuum mechanics