Équation de Helmholtz (mécanique des fluides)

En mécanique des fluides l'équation de Helmholtz est l'équation de transport de la vorticité (ou tourbillon) pour l'écoulement d'un fluide incompressible barotrope. Elle a été établie par Hermann von Helmholtz en 1858.

Cette équation est utilisée pour l'étude des mécanismes de turbulence et, dans le cas non visqueux, pour l'étude de l'atmosphère.

Équations de Navier-Stokes pour un milieu incompressible

Les équations de Navier-Stokes pour un milieu incompressible s'écrivent

• Équation d'incompressibilité

$\mathrm{\nabla }\cdot 𝐕=0$

• Équation de bilan de la quantité de mouvement

$\frac{\partial 𝐕}{\partial t}+\mathrm{\nabla }\cdot \left(𝐕𝐕\right)=-\frac{1}{\rho }\mathrm{\nabla }p+\nu {\mathrm{\nabla }}^2𝐕+𝐠$

où ρ est la masse volumique supposée variable, p la pression, V la vitesse et ν la viscosité cinématique du fluide. g est une force extérieure.

Équation de Helmholtz

L'équation de quantité de mouvement ci-dessus peut s'écrire

$\frac{\partial 𝐕}{\partial t}+\frac{1}{2}\mathrm{\nabla }{V}^2+\mathrm{\Omega }\times 𝐕=-\frac{1}{\rho }\mathrm{\nabla }p+\nu {\mathrm{\nabla }}^2𝐕+𝐠$

On a introduit dans cette équation la vorticité

$\mathrm{\Omega }=\mathrm{\nabla }\times 𝐕$

en prenant le rotationnel de cette nouvelle équation il vient (on suppose g irrotationnel)

$\frac{\partial \mathrm{\Omega }}{\partial t}+\mathrm{\nabla }\times (\mathrm{\Omega }\times 𝐕)=\frac{1}{{\rho }^2}\mathrm{\nabla }\rho \times \mathrm{\nabla }p+\nu {\mathrm{\nabla }}^2\mathrm{\Omega }$

Le second terme se développe de la façon suivante

$\mathrm{\nabla }\times (\mathrm{\Omega }\times 𝐕)=(𝐕\cdot \mathrm{\nabla })\mathrm{\Omega }-𝐕\underset{\equiv 0}{\underbrace{\mathrm{\nabla }\cdot \mathrm{\Omega }}}-(\mathrm{\Omega }\cdot \mathrm{\nabla })𝐕+\mathrm{\Omega }\underset{=0}{\underbrace{\mathrm{\nabla }\cdot 𝐕}}$

Finalement l'équation de conservation de la vorticité s'écrit

$\frac{\partial \mathrm{\Omega }}{\partial t}+(𝐕\cdot \mathrm{\nabla })\mathrm{\Omega }=(\mathrm{\Omega }\cdot \mathrm{\nabla })𝐕+\nu {\mathrm{\nabla }}^2\mathrm{\Omega }+\mathrm{\Pi }$

où Π est le vecteur barocline

$\mathrm{\Pi }=\frac{1}{{\rho }^2}\mathrm{\nabla }\rho \times \mathrm{\nabla }p$

Ce terme est nul pour un fluide barotrope. On obtient alors l'équation de Helmholtz

$\frac{\partial \mathrm{\Omega }}{\partial t}+(𝐕\cdot \mathrm{\nabla })\mathrm{\Omega }=(\mathrm{\Omega }\cdot \mathrm{\nabla })𝐕+\nu {\mathrm{\nabla }}^2\mathrm{\Omega }$