Équation de bilan de la quantité de mouvement/Forme globale

Début de la boite de navigation du chapitre

fin de la boite de navigation du chapitre

La forme globale de l'équation de bilan de la quantité de mouvement est un équation équivalente à la forme générale du principe fondamental de la dynamique. Pour un volume $V$ donné, on peut écrire que la dérivée totale de la quantité de mouvement est égale à la somme des forces qui s'exercent sur le fluide contenu dans le volume $V$.

$\frac{\mathrm{d}\overrightarrow{p}}{\mathrm{d}t}=\sum _V\overrightarrow{F}$.

Il existe deux familles de forces :

• les forces de contact, surfaciques (pression, viscosité),
• les forces volumiques, à distance (poids, force électromagnétique).

L'équation peut se décomposer ainsi :

$\frac{\mathrm{d}\overrightarrow{p}}{\mathrm{d}t}=\overrightarrow{F_S}+\overrightarrow{F_V}$.

Nous allons établir plusieurs expressions de la forme intégrale, ou forme globale, de l'équation de bilan de la quantité de mouvement.

Dérivée particulaire

La dérivée particulaire de la quantité de mouvement

$\frac{\mathrm{d}\overrightarrow{p}}{\mathrm{d}t}=\frac{\mathrm{d}}{\mathrm{d}t}\int \int {\int }_V\rho \overrightarrow{v}\mathrm{d}V$

peut s'exprimer de plusieurs façons :

$\frac{\mathrm{d}\overrightarrow{p}}{\mathrm{d}t}=\int \int {\int }_V(\frac{\partial (\rho \overrightarrow{v})}{\partial t}+\overrightarrow{\mathrm{d}\mathrm{i}\mathrm{v}}(\rho \overrightarrow{v}\otimes \overrightarrow{v}))\mathrm{d}V=\int \int {\int }_V(\frac{\mathrm{d}(\rho \overrightarrow{v})}{\mathrm{d}t}+\rho \overrightarrow{v}\text{div}\overrightarrow{v})\mathrm{d}V=\int \int {\int }_V\rho \frac{\mathrm{d}\overrightarrow{v}}{\mathrm{d}t}\mathrm{d}V$.

Ici la deuxième expression peut se simplifier : en partant de la dérivée particulaire de la densité volumique de quantité de mouvement

$\frac{\mathrm{d}(\rho \overrightarrow{v})}{\mathrm{d}t}=\rho \frac{\mathrm{d}\overrightarrow{v}}{\mathrm{d}t}+\frac{\mathrm{d}\rho }{\mathrm{d}t}\overrightarrow{v}$,

puis en ajoutant de part et d'autre le terme $\rho \overrightarrow{v}\text{div}\overrightarrow{v}$, il vient :

$\frac{\mathrm{d}(\rho \overrightarrow{v})}{\mathrm{d}t}+\rho \overrightarrow{v}\text{div}\overrightarrow{v}=\rho \frac{\mathrm{d}\overrightarrow{v}}{\mathrm{d}t}+\underset{=0}{\underbrace{(\frac{\mathrm{d}\rho }{\mathrm{d}t}+\rho \text{div}\overrightarrow{v})}}\overrightarrow{v}$,

où l'on reconnaît, au dessus de l'accolade, la forme non-conservative de l'équation de continuité.

La dérivée particulaire de la quantité de mouvement prend alors un expression très simple :

$\frac{\mathrm{d}\overrightarrow{p}}{\mathrm{d}t}=\int \int {\int }_V\rho \frac{\mathrm{d}\overrightarrow{v}}{\mathrm{d}t}\mathrm{d}V$.

Forces volumiques

Les forces de volume peuvent s'exprimer simplement :

$\overrightarrow{F_V}=\int \int {\int }_V\overrightarrow{f_V}\mathrm{d}V$.

Dans la grande majorité des cas en mécanique des fluides, la seule force à distance que subit le fluide est son poids. Dans ce cas :

$\overrightarrow{F_V}=\int \int {\int }_V\mathrm{d}m\overrightarrow{g}=\int \int {\int }_V\rho \overrightarrow{g}\mathrm{d}V$.

Cette force est elle-même souvent négligée dans le cas des gaz.

Forces surfaciques

Les forces de surface qui s'appliquent sur la surface $S$ fermée, frontière du volume $V$ sont un peu plus subtiles à détailler.

On peut noter :

$\overrightarrow{F_S}=\int \subset \supset {\int }_S\overrightarrow{\mathrm{d}{F}_S}$

Elles sont modélisées par une matrice nommée tenseur des contraintes.

$\overline{\stackrel{‾}{\mathrm{\tau }}}=\left(\begin{array}{ccc}{\sigma }_{xx}& {\tau }_{xy}& {\tau }_{xz}\\ {\tau }_{yx}& {\sigma }_{yy}& {\tau }_{yz}\\ {\tau }_{zx}& {\tau }_{zy}& {\sigma }_{zz}\end{array}\right)$.

de sorte que

$\overrightarrow{\mathrm{d}{F}_S}=\overline{\stackrel{‾}{\tau }}\times \overrightarrow{\mathrm{d}S}=\left(\begin{array}{ccc}{\sigma }_{xx}& {\tau }_{xy}& {\tau }_{xz}\\ {\tau }_{yx}& {\sigma }_{yy}& {\tau }_{yz}\\ {\tau }_{zx}& {\tau }_{zy}& {\sigma }_{zz}\end{array}\right)\times \left(\begin{array}{c}\mathrm{d}{S}_x\\ \mathrm{d}{S}_y\\ \mathrm{d}{S}_z\end{array}\right)=\left(\begin{array}{c}{\sigma }_{xx}\mathrm{d}{S}_x+{\tau }_{xy}\mathrm{d}{S}_y+{\tau }_{xz}\mathrm{d}{S}_z\\ {\tau }_{yx}\mathrm{d}{S}_x+{\sigma }_{yy}\mathrm{d}{S}_y+{\tau }_{yz}\mathrm{d}{S}_z\\ {\tau }_{zx}\mathrm{d}{S}_x+{\tau }_{zy}\mathrm{d}{S}_y+{\sigma }_{zz}\mathrm{d}{S}_z\end{array}\right)$.

On peut exprimer les forces de surface :

$\overrightarrow{F_S}=\int \subset \supset {\int }_S\overline{\stackrel{‾}{\tau }}\times \overrightarrow{\mathrm{d}S}=\int \int {\int }_V\overrightarrow{\text{div}}\overline{\stackrel{‾}{\mathrm{\tau }}}\mathrm{d}V$,

où $\overrightarrow{\text{div}}\overline{\stackrel{‾}{\mathrm{\tau }}}$ est la divergence de la matrice $\overline{\stackrel{‾}{\mathrm{\tau }}}$ qui est définie de la façon suivante :

$\overrightarrow{\text{div}}\overline{\stackrel{‾}{\mathrm{\tau }}}=\left(\begin{array}{c}\frac{\partial {\sigma }_{xx}}{\partial x}+\frac{\partial {\tau }_{xy}}{\partial y}+\frac{\partial {\tau }_{xz}}{\partial z}\\ \frac{\partial {\tau }_{yx}}{\partial x}+\frac{\partial {\sigma }_{yy}}{\partial y}+\frac{\partial {\tau }_{yz}}{\partial z}\\ \frac{\partial {\tau }_{zx}}{\partial x}+\frac{\partial {\tau }_{zy}}{\partial y}+\frac{\partial {\sigma }_{zz}}{\partial z}\end{array}\right)$.

Équation bilan

Le bilan de la quantité de mouvement peut donc s'écrire :

$\int \int {\int }_V(\frac{\partial (\rho \overrightarrow{v})}{\partial t}+\overrightarrow{\mathrm{d}\mathrm{i}\mathrm{v}}(\rho \overrightarrow{v}\otimes \overrightarrow{v}))\mathrm{d}V=\int \int {\int }_V\overrightarrow{f_V}\mathrm{d}V+\int \subset \supset {\int }_S\overline{\stackrel{‾}{\tau }}\times \overrightarrow{\mathrm{d}S}$,

ou encore plus simplement

$\int \int {\int }_V\rho \frac{\mathrm{d}\overrightarrow{v}}{\mathrm{d}t}\mathrm{d}V=\int \int {\int }_V\overrightarrow{f_V}\mathrm{d}V+\int \subset \supset {\int }_S\overline{\stackrel{‾}{\tau }}\times \overrightarrow{\mathrm{d}S}$.