Fréquence de Brunt-Väisälä
La fréquence de Brunt-Väisälä (ou Brunt-Vaisala) est la fréquence d'oscillation d'une particule fluide déplacée verticalement dans un environnement stable autour de sa position initiale paramétrisée par David Brunt et Vilho Väisälä. Elle correspond à la fréquence d'une onde de gravité qui joue un rôle très important dans les échanges énergétiques des écoulements géophysiques, notamment en dynamique atmosphérique et pour l'océanographie physique. Par exemple, entre autres paramètres, la fréquence de Brunt-Väisälä contrôle la hauteur et l'espacement entre les rues de cumulus ou les altocumulus lenticularis en aval de montagnes, ainsi que celui entre les crêtes de houle en pleine mer.
![](http://up.wiki.x.io/wikipedia/commons/thumb/b/b1/Wave_cloud.jpg/220px-Wave_cloud.jpg)
Théorie
modifierLe concept de l'oscillation et de la fréquence de Brunt-Väisälä provient de l'application de la deuxième des trois lois de Newton dans un milieu stablement verticalement stratifié. On peut expliquer la nature oscillatoire des fluides stratifiés en pensant à une particule de fluide dont la densité augmente avec la profondeur. Lorsqu'elle se trouve déplacée verticalement en dehors de sa position d'équilibre, sa densité devient plus grande ou plus faible que le fluide environnant et une force de restitution excédentaire, la pesanteur ou la poussée d'Archimède respectivement, apparaît et tend à la ramener vers le point d'équilibre[1]. En général, la particule dépasse l'équilibre sur son chemin de retour, car la force a induit une accélération. Ce phénomène, entretenu, déclenche une oscillation dont la « fréquence » (ou, à proprement parler, la pulsation[2]) est[3],[4] :
où est l'accélération locale de la pesanteur, est le déplacement de la parcelle et est la densité potentielle définie comme la densité qu'aurait une parcelle de fluide déplacée adiabatiquement à une pression de référence (souvent choisie comme un bar dans le cas de l'atmosphère terrestre).
Pour un gaz parfait, on a l'égalité : où est la masse volumique, la pression et l'indice adiabatique de l'air. Avec les variables thermodynamiques usuelles, on peut donc écrire
Le calcul de cette fréquence est une façon de connaître la stabilité de l'air :
- une oscillation se produit si et seulement si , c'est-à-dire si le nombre sous la racine carrée est positif et la racine carrée est donc réelle. Cela impose que soit négatif, ce qui se traduit par le fait que le gradient de densité est négatif (la stratification du milieu doit être telle que la densité diminue avec l'altitude) ; dans le cas contraire, le nombre sous la racine est négatif et sa racine carrée est un nombre imaginaire pur. L'interprétation physique en est que l'oscillation se dissipe, comme c'est le cas dans un fluide dont la stratification n'est pas stable et où se produit de la convection : la parcelle déplacée devient par exemple moins dense que son environnement et accélère dans la même direction que le déplacement initial (pas d'oscillation) ;
- si , la stabilité est « neutre » car cela signifie qu'il n'y a pas de variation de densité. La parcelle déplacée demeura à sa nouvelle altitude (atmosphère) ou profondeur (océan) une fois le déplacement terminé.
Dans l'atmosphère
modifierLa densité est reliée directement à la température et au contenu en vapeur d'eau de la parcelle d'air. Soit la température potentielle. L'équation devient, dans ce milieu[5] :
- , où est l'altitude au-dessus du sol.
Dans l'atmosphère terrestre typique, la valeur de N[2] est de 0,012 s−1. La période de l'oscillation étant , elle est de l'ordre de huit minutes[5].
Modèle théorique
Il va être démontré que dans une masse d'air stable, une parcelle d'air à laquelle on a apporté une perturbation va osciller verticalement avec la fréquence N définie par :
où est la température potentielle à l'altitude z. Ensuite, on va définir le nombre de Froude qui est déduit de l'équation établie dans cette section. Ce nombre de Froude prédit l'existence (ou non existence) d'un phénomène de blocage. La démonstration détaillée qui suit a pour base la référence[6].
On considère un volume de contrôle de surface S compris entre la hauteur z et la hauteur z + δ z où δ z est une quantité infiniment petite. On suppose que la pression à l'altitude z est p(z) et à l'altitude z + δ z, la pression est p(z + δ z). Soit ρ la masse volumique de l'air. La masse de la parcelle d'air est donc m = ρ S δ z. Les forces qui s'appliquent à la parcelle d'air sont:
la pression sur la face inférieure : p(z) S
la pression sur la face supérieure : - p(z+ δ z) S
la gravité : - g ρ S δ z
La force exercée sur la parcelle d'air est donc :
L'accélération a de la parcelle d'air sera donc :
On obtient donc (on note que δ z est un infiniment petit) :
On peut simplifier et donc :
En outre, soit ρ₀ la densité de l'air extérieur. On a
Donc,
Finalement,
Et donc:
On utilise la loi des gaz parfaits. On a
où T est la température absolue de la parcelle d'air et Cp et Cv sont les chaleurs spécifiques à pression ou volume constants. Donc,
- .
La pression de la parcelle d'air p₀ est égale à la pression du milieu extérieur p. Soit T₀ la température extérieure. Donc :
- .
On obtient donc :
Donc,
L'air est un gaz di-atomique et donc :
- .
On définit la température potentielle comme suit :
La température potentielle est donc la température qu'aurait la parcelle d'air si elle était comprimée adiabatiquement à la pression standard au niveau de la mer. Comme la pression de la parcelle d'air est égale à la pression extérieure, si et sont les températures potentielles respectives de la parcelle d'air et de l'air extérieur, on a:
On obtient donc :
Plus spécifiquement, on écrit :
De même :
Donc,
Donc,
La parcelle d'air est en ascension adiabatique et donc : et donc :
On note que et donc:
On suppose maintenant que ne dépend pas de z et donc est uniforme. On définit la quantité N² comme suit :
- .
On a donc en première approximation :
On note que :
où ε est un nombre infiniment petit. On obtient donc :
L'ensemble des nombres hyperréels est un corps et donc, l'on peut diviser l'équation ci-dessus par δz et donc,
En éliminant ε qui est infiniment petit, on obtient donc :
En intégrant cette équation, on obtient donc :
On a de même
et donc :
Par définition :
On obtient alors l'équation différentielle linéaire suivante :
La solution générale de cette équation différentielle s'écrit :
où a et b sont 2 constantes dépendant des conditions initiales.
Supposons qu'à t = 0, la vitesse verticale de la parcelle d'air est w et que h = 0. La solution de l'équation différentielle ci-dessus s'écrit :
La déflexion maximale en hauteur de la parcelle d'air est w/N. Par conséquent, lorsqu'une masse d'air de vitesse horizontale u rencontre une montagne de hauteur h > u/N, la masse d'air ne pourra pas franchir la montagne et l'on sera en présence d'un phénomène de blocage en amont. Le critère d'existence ou de non existence de phénomène de blocage est déterminé par la valeur du nombre de Froude météorologique défini par :
Si le nombre de Froude est plus grand que 1, il n' y a pas de blocage et dans le cas contraire il y a blocage.
La quantité N définie par :
- .
est appelée fréquence de Brunt-Väisälä.
Application numérique
Dans l'atmosphère standard, le gradient adiabatique est g/C_p = 9.75 K/km et l'on a d T / d z = - 6.5 K/km. Donc, on a
On obtient donc :
Dans l'océan
modifierDans l'océan, la densité in-situ dépend de la température T , de la salinité S et de la pression p : .
La variation de densité n'est pas linéaire selon la température (la densité maximale de l'eau non salée est à 4 °C et la densité change soudainement dans la couche de glace de surface, entre autres facteurs). Lorsqu'une particule de fluide est déplacée verticalement de manière adiabatique (c'est-à-dire sans que T et S ne soient modifiées) la variation de densité due à une variation de niveau est[7] :
Où :
- est la compressibilité ;
- est la variation de la pression p avec la profondeur z orientée vers la surface ;
- est la densité moyenne (due à l'approximation de Boussinesq, généralement , de sorte que lorsque la particule se déplace vers la surface ( ) la densité diminue ( car puisque z est orienté vers la surface).
C'est cette densité modifiée par la pression qu'il faut comparer à la densité environnante pour obtenir la fréquence de Brünt-Väisälä :
Cette formule peut aussi s'écrire en termes de densité potentielle référencée localement :
Propriétés
modifierLes ondes de gravité ont plusieurs propriétés qui s'interprètent à partir de leur fréquence, parmi lesquelles on remarque :
- la direction de propagation de ces ondes dépend de la fréquence du forçage et aussi de la fréquence de Brunt-Väisälä locale (stratification de densité locale) ;
- la vitesse de phase (vitesse de propagation des fronts d'onde) et la vitesse de groupe (vitesse avec laquelle l'énergie des ondes est transmise) des ondes internes sont perpendiculaires.
En utilisant l'approximation de Boussinesq, on peut trouver la relation de dispersion des ondes générées par :
où est la fréquence d'excitation utilisée, est la fréquence de Brunt-Väisälä et est l'angle du vecteur de propagation par rapport à la horizontal.
Bibliographie
modifier- James R. Holton, An Introduction to Dynamic Meteorology, 4e édition, New York, Academic Press, 535 p. (ISBN 978-0-12-354015-7 et 0-12-354015-1, lire en ligne), p. 50-53
- Lighthill, J., Waves in Fluids, Cambridge University Press,
- Mowbray,D.E. et B.S.H.Rarity, « A theoretical and experimental investigation of the phase configuration of internal waves of small amplitude in a density stratified liquid », Journal of Fluid Mechanics, no 28, , p. 1-16
- Rogers, R. R. et Yau, M. K., Short Course in Cloud Physics, 3e édition, Butterworth-Heinemann, , 304 p. (ISBN 0-7506-3215-1), p. 30-35EAN 9780750632157
- Tritton, D.J., Physical Fluid Dynamics. 2e édition, Oxford University Press,
Notes et références
modifier- Occurrences en français de « fréquence de Brunt-Vaisala », sur TERMIUM Plus, la banque de données terminologiques et linguistiques du gouvernement du Canada..
- Bien qu'on emploie souvent le terme « fréquence », il s'agit en fait d'une vitesse angulaire fictive, ou pulsation ; ainsi le nombre d'oscillations par seconde n'est pas de mais de , et la période n'est pas de mais de secondes
- L. Gostiaux, T. Dauxois et E. Falcon, École normale supérieure de Lyon, « Réflexion critique d’ondes internes de gravité en fluides stratifiés », Compte-rendus [sic] de la 8e Rencontre du Non-Linéaire, Paris, Université Paris-Diderot, (lire en ligne [archive du ] [PDF], consulté le ).
- Radyadour Kh Zeytounian, Modélisation asymptotique en mécanique des fluides newtoniens, Springer, , 226 p. (ISBN 9783540578383, lire en ligne), p. 152.
- Rogers et Yau, p. 33-35
- (en) James R Holton, Introduction to dynamic meteorology (Fourth edition), vol. 88, Amsterdam, Elsevier Academic press, , 526 p. (ISBN 0-12-354015-1, lire en ligne), p. 52
- (en) G. K. Vallis, Atmospheric and Oceanic Fluid Dynamics, Cambridge,
Voir aussi
modifierArticles connexes
modifierLien externe
modifier- Laboratoire d'hydrodynamique, « Ondes de gravité internes », École polytechnique (France) (consulté le ).