Dynamique des fluides parfaits/Équations de Bernoulli

On reprend l'équation d'Euler locale (PFD local pour un fluide parfait) :

$-\overrightarrow{\mathrm{\nabla }p}+\rho \overrightarrow{g}=\rho .\overrightarrow{a}$

$\rho [\frac{\partial \overrightarrow{v}}{\partial t}+0.5\mathrm{\nabla }\stackrel{2}{\overrightarrow{v}}+(\overrightarrow{\mathrm{\nabla }}\wedge \overrightarrow{v})\wedge \overrightarrow{v}]=\rho \overrightarrow{g}-\overrightarrow{\mathrm{\nabla }p}$

Nous allons soumettre cette équation à quelques cas particuliers...

Pour une particule

Soit une particule de position $\overrightarrow{r}$ de paramétrage : $\overrightarrow{r}=\left(\begin{array}{c}x\\ y\\ z\end{array}\right)$. Ainsi $\overrightarrow{d_r}=\left(\begin{array}{c}dx\\ dy\\ dz\end{array}\right)$

Dans le cas d'un écoulement irrotationnel, le déplacement de n’importe quelles particules du fluide est en permanence parallèle à la direction globale de ce fluide : le produit vectoriel entre le déplacement et la vitesse est donc nul, ce qui implique :

$(\overrightarrow{\mathrm{\nabla }}\wedge \overrightarrow{v})\wedge \overrightarrow{v}=(\overrightarrow{\mathrm{\nabla }}\overrightarrow{v}\wedge \overrightarrow{v}).(\overrightarrow{v}\wedge \overrightarrow{dr})=0$

Il restera donc de l'équation d'Euler locale :

$\rho [\frac{\partial \overrightarrow{v}}{\partial t}+0.5\mathrm{\nabla }\stackrel{2}{\overrightarrow{v}}]=\rho \overrightarrow{g}-\overrightarrow{\mathrm{\nabla }p}$

Soit :

$\frac{\partial \overrightarrow{v}}{\partial t}.\overrightarrow{d_r}+0.5\mathrm{\nabla }\stackrel{2}{\overrightarrow{v}}.\overrightarrow{d_r}=\overrightarrow{g}.\overrightarrow{d_r}-\frac{1}{\rho }\overrightarrow{\mathrm{\nabla }p}.\overrightarrow{d_r}$

On peut aussi écrire cette expression sous cette forme :

$(\rho \frac{\partial \overrightarrow{v}}{\partial t}).\overrightarrow{d_r}+\mathrm{\nabla }(0.5\rho \stackrel{2}{\overrightarrow{v}}+p)\overrightarrow{d_r}-\rho \overrightarrow{g}.\overrightarrow{d_r}=0$

Pour un fluide parfait compressible ou incompressible soumis à l'accélération de pesanteur g

On se donne un axe vertical orienté vers le haut tel que :

$\overrightarrow{g}=\left(\begin{array}{c}0\\ 0\\ -g\end{array}\right)$

Mathématiquement, on peut écrire également :

$\overrightarrow{g}.\overrightarrow{dr}=-\overrightarrow{\mathrm{\nabla }}(gz).\overrightarrow{dr}$

Il vient finalement le théorème de Bernoulli généralisé valable pour un fluide parfait compressible ou incompressible lorsque l'axe vertical est orienté vers le haut :

$(\rho \frac{\partial \overrightarrow{v}}{\partial t}).\overrightarrow{d_r}+\mathrm{\nabla }(0.5\rho \stackrel{2}{\overrightarrow{v}}+p+\rho gz)\overrightarrow{d_r}=0$

On intègre l'équation le long d'une ligne de courant entre le point 1 et le point 2 :

$\underset{1}{\overset{2}{\int }}(\rho \frac{\partial \overrightarrow{v}}{\partial t}).\overrightarrow{d_r}+\underset{1}{\overset{2}{\int }}\mathrm{\nabla }(0.5\rho \stackrel{2}{\overrightarrow{v}}+p+\rho gz)\overrightarrow{d_r}=0$

Un écoulement est toujours irrotationnel le long d'une ligne de courant. Il faut retenir que l'équation de Bernoulli est valable pour tout écoulement irrotationnel, même si on ne suit pas une ligne de courant.

Pour un fluide parfait incompressible soumis à l'accélération de pesanteur g

La masse volumique du fluide est alors constante. L'expression précédente devient :

$\underset{1}{\overset{2}{\int }}(\rho \frac{\partial \overrightarrow{v}}{\partial t}).\overrightarrow{d_r}+[0.5\rho \stackrel{2}{\overrightarrow{v}}+p+\rho gz)\overrightarrow{d_r}{]}_1^2=0$

C'est le théorème de Bernoulli pour un fluide parfait incompressible le long d'une ligne de courant.

Pour l'écoulement permanent d'un fluide parfait

$\frac{\partial v}{\partial t}=0\Rightarrow [0.5\rho v^2+p+\rho gz)\overrightarrow{d_r}{]}_1^2=0\Rightarrow 0.5\rho ({v_2}^2-{v_1}^2)+(p_2-p_1)+\rho g(z_2-z_1)=0$

Cette dernière formule est valable pour un fluide parfait incompressible dans un écoulement permanent le long d'une ligne de courant.

$(0.5\rho v^2)$ est la pression dynamique. $(p+\rho gz)$ est la pression motrice. $(0.5\rho v^2+p+\rho gz)$ est la pression totale d'arrêt.

Modèle:Bas de page