# Rotor (météorologie)

Un rotor est un mouvement tourbillonnaire de l'air associé à des ondes orographiques. Il est localisé sous la partie laminaire et peut donner lieu à des turbulences extrêmement violentes qui peuvent briser un aéronef. Il est souvent matérialisé par des cumulus fractus (ou fractocumulus).

## Cause

Un rotor est un phénomène météorologique associé à des ondes de gravité. L'air va osciller verticalement en aval de l'obstacle et en application directe des lois de conservation, des fluctuations du vent horizontal vont se produire. Des rouleaux vont donc se former et sont généralement le siège de turbulences qui sont amplifiées en fin de journée lorsque le sol est surchauffé.
Les ondes orographiques produisent des ascendances (et descendances) très puissantes où la vitesse verticale peut dépasser 10 m/s. Toutefois, près du sol l'air est instable car l'air en redescendant le long de la montagne se réchauffe suivant le gradient adiabatique sec (g/Cp). De plus, par temps clair, la sous-couche va être réchauffée à partir du sol et donc des ascendances thermiques vont se former. La combinaison de ces ascendances thermiques et des mouvements tourbillonnaires associé aux ondes orographiques va rendre l'air turbulent. Des sous-tourbillons vont se former.
La transition entre une rafale ascendante et une rafale descendante à l'intérieur du rotor peut s'effectuer sur une cinquantaine de mètres environ. Des sous-tourbillons sont également formés à la limite amont entre le rotor principal et l'onde orographique proprement dite$^,$. Ceci est dû à un effet de cisaillement entre le tourbillon (rotor) et l'onde qui génère une instabilité de Kelvin-Helmholtz.

## Types de rotor

Rolf Hertenstein a défini deux classes de rotors:
* les rotors associés à des ondes piégées (type I)
* et les rotors associés à un ressaut hydraulique (type II).

### Rotor associé à des ondes piégées

Ce type de rotor est de loin le plus courant. Les ondes sont quasi stationnaires et le mouvement de l'air dans le rotor lui-même est quasi circulaire. En général la turbulence à l'intérieur du rotor est de faible à sévère mais en pratique jamais extrême. Ces rotors sont souvent matérialisés par une ligne de nuages en rouleau constituée de cumulus ou cumulus fractus. Ces nuages, peu épais, se situent en général au niveau de la ligne de crête.
Une explication simple et élégante de la formation de rotors associés à des ondes piégées a été donnée par Paul Queney$^,$ et Richard Scorer.
Cette théorie est basée sur une version simplifiée de l'équation de Scorer. Cette théorie est appelée modèle de l'œil de chat car les lignes de courant en aval d'une chaîne de montagnes ont la forme d'un œil de félin (figure 5 dans la boîte déroulante, elles sont aussi montrées de manière approximative dans la figure 8). Cette théorie est basée sur la loi de conservation de la masse. En particulier, cette théorie explique pourquoi le vent au sol sous un rotor est souvent inversé et que le rotor a la forme d'un tambour horizontal. Ce modèle est un modèle bi-dimensionnel en (x,z) et donc s'appuie sur des hypothèses simplificatrices.
On considère un modèle en 2 dimensions où la chaîne de montagnes est rectiligne et de hauteur constante. On suppose qu'à l'amont de la montagne, le champ de vitesses est axisymétrique et est de la forme :
${\vec {U}}(z)=U{\vec {i}}$
où $U=\gamma z$ pour $z\leq z_{M}$ et que $U=\gamma z_{M}=U_{M}$ sinon.
On suppose qu'à l'aval de la montagne se forment des ondes stationnaires que l'on modèle par un champ de vitesses de type ${\vec {U}}+{\vec {u}}$ et l'on cherche à exprimer ${\vec {u}}$.
On écrit : ${\vec {u}}=u{\vec {i}}+w{\vec {k}}$
Comme la vitesse du vent est très nettement inférieure à la vitesse du son, on peut considérer le fluide comme étant incompressible et donc l'on a :
$\nabla {\vec {u}}=0$
et donc il existe une fonction de courant $\psi$ telle que :
$u=-{\partial \psi \over \partial z}\qquad w={\partial \psi \over \partial x}$
On suppose que le flot est barotropique et que la perturbation du flot est irrotationnelle. On a donc :
${\partial w \over \partial x}-{\partial u \over \partial z}=0$
On vérifie alors aisément que
$\Delta \psi ={\partial ^{2}\psi \over \partial x^{2}}+{\partial ^{2}\psi \over \partial z^{2}}=0$
On suppose que l'extension horizontale du phénomène par rapport à l'extension verticale est très grande et donc que :
$\left|{\partial ^{2}\psi \over \partial x^{2}}\right|\ll \left|{\partial ^{2}\psi \over \partial z^{2}}\right|$
On obtient alors :
${\partial ^{2}\psi \over \partial z^{2}}\approx {\partial ^{2}\psi \over \partial x^{2}}+{\partial ^{2}\psi \over \partial z^{2}}=0$
La solution générale de cette équation laplacienne est donc :
$\psi (x,z)=f(x)+zg(x)$
On définit la fonction de courant globale $\Psi$ comme étant :
$\Psi (x,z)={1 \over 2}\gamma z^{2}+f(x)+zg(x)$
On vérifie aisément que :
$w(x,z)={\partial \Psi \over \partial x}=f'(x)+zg'(x)$
Il est clair qu'à z = 0, on a w(x,0) = 0. On obtient donc $\forall x\quad f'(x)=0$ et on peut donc supposer que $f(x)=0$
On obtient donc :
$\Psi (x,z)={1 \over 2}\gamma z^{2}+zg(x)$
On suppose que la fonction $\Psi$ est périodique et que g(x) peut être représentée par une fonction sinusoïdale. On écrit alors :
$\Psi (x,z)={1 \over 2}\gamma z^{2}+az\sin(kx)$
Soit N la fréquence de Brunt-Väisälä. On suppose que le nombre d'onde k vaut :
$k={N \over U_{M}}$
La vitesse du vent vertical est donnée par :
    $w(x,z)={\partial \Psi \over \partial x}={\partial \left({1 \over 2}\gamma z^{2}+az\sin(kx)\right) \over \partial x}$
Et donc :
    $w(x,z)=akz\cos(kx)$
k est le nombre d'onde défini par :
    $k={N \over U_{M}}$
où N est la fréquence de Brunt-Väisälä.
On déterminera a infra.
On rappelle que $u(x,z)=-{\partial \psi \over \partial z}$ et donc :
$u(x,z)=-a\sin(kx)$
Le vent horizontal en aval sera donc :
$U_{h}=U+u=\gamma z-a\sin(kx)$
On constate alors que près du sol, le vent sera inversé lorsque $kx=\left(n+{1 \over 2}\right)\pi$ c'est-à-dire:
$x={n\pi \over k}$
On remarque que
${d \over dt}\Psi (x(t),z(t))={\partial \Psi \over \partial x}{dx \over dt}+{\partial \Psi \over \partial z}{dz \over dt}={\partial \Psi \over \partial x}(U+u)+{\partial \Psi \over \partial z}w=-w(U+u)+(U+u)w=0$
La trajectoire d'une parcelle d'air suivra donc les lignes de niveau définies par $\Psi (x,z)=C$.
On va démontrer que la circulation d'air se décompose en une boucle fermée appelée rotor et un courant synoptique.
Une ligne de courant est donc définie par :
$C=\Psi (x,z)={1 \over 2}\gamma z^{2}+az\sin(kx)$
On pose $t=\tan {kx \over 2}$ On rappelle que $\sin(kx)={2t \over 1+t^{2}}$
On substitue et donc :
${1 \over 2}\gamma z^{2}+az{2t \over 1+t^{2}}-C=0$
Une telle courbe en (z,t) est une courbe elliptique et donc est une courbe de genre 1. Une telle courbe est donc constituée d'une boucle et d'une courbe partant à l'infini.
Dans le cas C = 0 la courbe est polynomiale (et se factorise) et l'on a :
$z={4at \over \gamma (1+t^{2})}={2a \over \gamma }\sin(kx)$
Les solutions explicites en ce qui concerne les lignes de courant sont les suivantes :
$z={-a\sin(kx)\pm {\sqrt {a^{2}\sin ^{2}(kx)+\gamma C}} \over \gamma }$
Dans le cas C = 0, on a :
$z={-2a\sin(kx) \over \gamma }$
Le sommet du rotor est donc situé à $z_{r}={2a \over \gamma }$. Si l'on fait l'hypothèse supplémentaire que le rotor se situe au niveau de la ligne de crête, connaissant γ (donné par les conditions synoptiques), on peut en déduire a. Soit h la hauteur de la montagne, on a alors :
    $h=z_{r}={2a \over \gamma }$
et donc :
    $a={h\gamma \over 2}$
On obtient donc la vitesse verticale suivante :
    $w(x,z)={h\gamma Nz \over 2U_{M}}\cos(kx)$
Donc, si l'on prend h = 2000 m, γ = 0.02 s¯¹, N= 0.01 Hz, U$_M$ = 40 m/s, à 1000 mètres d'altitude, le courant ascendant maximum sera :
    $w_{M}={2\times 10^{3}\times 2\times 10^{-2}\times 10^{-2}\times 1\times 10^{3} \over 2\times 4\times 10^{1}}=2.5$ m/s
À 2000 mètres, l'ascendance sera de 5 m/s. Ces valeurs numériques correspondent grossièrement aux courants ascendants habituellement relevés.
La figure ci-dessous représente les lignes de courant pour k = 0.001 m¯¹, a = 20 m/s, γ = 0.02 s¯¹.
La turbulence est maximale en amont et au-dessus du rotor. De plus, les vitesses ascensionnelles les plus élevées se trouvent juste en amont du rotor comme les courbes de température isentropique le montrent dans la référence (Grubišić, 2004). De plus, lorsque le rotor est matérialisé par un nuage-mur, l'ascendance devient encore plus forte à cause de la libération de chaleur latente. Dans les cas extrêmes, on peut avoir affaire à un rotor de type II qui est décrit ci-dessous.

### Rotor de ressaut hydraulique

Ce type de rotor est rare et est semble-t-il associé à une inversion marquée de température au sommet de la montagne. Les simulations numériques effectuées semblent confirmer le phénomène. En outre ces rotors ont tendance à se produire en fin d'après midi lorsque le sol est surchauffé. Ces rotors sont extrêmement violents et peuvent atteindre des altitudes élevées (jusqu'à 7 km au-dessus du sol). Ces rotors peuvent être matérialisés par des nuages dont le flanc au vent est quasi vertical et dont la hauteur peut atteindre 8 000 à 9 000 mètres. Ceci peut s'expliquer par le fait que l'air est en forte ascension dans le nuage et donc empêche les gouttelettes d'eau ou cristaux de glace de redescendre vers le sol. À titre d'illustration, l'image ci-contre montre un ressaut hydraulique où la turbulence est semble-t-il intense. La figure 5 de la référence Hertenstein, 2005 montre un nimbostratus de grande extension verticale associé à un ressaut hydraulique.
Au cours du Sierra Wave Project, un pilote de planeur (Larry Edgar) avait rencontré le 25 avril 1955 une turbulence extrême en entrant dans un nuage de rotor qui grandit considérablement en l'espace de quelques secondes. Le planeur s'était engouffré dans le nuage et le planeur eut subi des accélérations gigantesques (16 g)$^,$$^,$$^,$ qui brisèrent l'aéronef, blessèrent le pilote qui fut aussi victime d'un voile noir. Il semblerait que le pilote de planeur eût rencontré un rotor de type II (ressaut hydraulique). Le pilote bien que blessé put s'éjecter en parachute.
Les figures 6 et 7 montrent un phénomène de ressaut hydraulique en aval de chutes naturelle et artificielle. En particulier, la figure 7 montre un double ressaut hydraulique.

### Rotors migratoires

Ces rotors sont des tourbillons se déplaçant dans le sens du vent à une vitesse inférieure à celle du vent synoptique. Ils se forment en général en présence d'une activité thermique ou après le passage d'un front froid. Ils sont souvent matérialisés par des petits cumulus fractus et peuvent souvent être confondus avec des ascendances thermiques. Les fractocumulus marquant ces rotors ont un mouvement en apparence circulaire. Ces rotors ont une durée de vie relativement courte qui est de l'ordre de quelques minutes.

### Rotors associés à des cyclones tropicaux

Sous l'ouragan Fran, Joshua Wurman a identifié un mouvement tourbillonnaire d'axe horizontal s'apparentant à des rotors associés à des ondes orographiques. La distance entre 2 rouleaux était de l'ordre de 600 m. La vitesse du vent synoptique à 1 000 mètres était de 50 à 60 m/s (100 à 120 nœuds) alors qu'en certains endroits au sol, le vent n'était que de 15 à 35 m/s (30 à 70 nœuds). Ces rotors étaient parallèles à la direction générale du vent. La partie descendante du rotor engendrait des vents au sol de 120 nœuds tandis que la partie ascendante du rotor ralentissait les vents au sol de l'ordre de 30 à 70 nœuds.

## Vol à voile
### Turbulence

Soit d la distance entre une rafale montante et descendante, soit ±w la vitesse de la rafale verticale et soit v la vitesse horizontale de l'aéronef. L'accélération que l'aéronef va subir sera :
$a={2w\times v \over d}$
Donc, plus l'aéronef volera vite, plus les accélérations seront importantes.
La différence de vitesse verticale entre les deux rafales peut être de l'ordre de 20 m/s. Un aéronef (avion de ligne) qui traverserait cette zone à 300 nœuds (150 m/s) va donc subir une accélération verticale de 6 g (60 m/s$^2$) et va donc probablement se briser. Toutefois, les planeurs sont beaucoup plus lents. En supposant que le planeur traverse le rotor à 40 nœuds (20 m/s), l'accélération verticale sera seulement de 0,8 g ce qui est nettement en deçà des limites de vol.
La turbulence peut être définie de manière plus rigoureuse en introduisant la notion de tourbillon (parfois appelée « vorticité » par calque de l'anglais) qui est définie de la manière suivante :
    ${\vec {\eta }}={\vec {\nabla }}\times {\vec {v}}$
où v est le champ de vecteurs vitesse. Le tourbillon est exprimé en s¯¹. Ainsi, le facteur de charge subi par un aéronef sera :
    $a=g+{\eta \times v \over 2}$
Au cours de l'expérience T-REX, des tourbillons de 0,4 s¯¹ ont été estimés. Ainsi, un avion de ligne volant à 200 m/s subira un facteur de charge de 10 + 0.4 × 200/2 = 50 m/s$^2$. Ce qui correspondrait à 5 fois l'accélération de la gravité ce qui est suffisant pour détruire l'aéronef. Lors de l'accident du Boeing 737 à Colorado Springs, il a été estimé que le tourbillon a pu atteindre 0.6 s¯¹.
A contrario, un tourbillon de 0,4 Hz appliqué à un planeur volant à 20 m/s engendre un facteur de charge de seulement 8 m/s$^2$, ce qui est moins que l'accélération de la pesanteur et donc ne détruira en aucun cas ce type d'aéronef.
On peut aussi remarquer que le tourbillon est maximal à l'avant du rotor et donc la turbulence sera maximale dans la zone de transition entre le rotor et la partie laminaire de l'onde traduisant une instabilité de Kelvin-Helmholtz à la jonction des deux masses d'air$^,$.
La turbulence est « acceptable » du côté au vent (de faible à modérée) (partie ascendante du rotor), tandis qu'elle peut être sévère à extrême du côté sous le vent (partie descendante du rotor). Ainsi, un planeur peut exploiter la partie ascendante du rotor sans prendre trop de risques alors que la partie descendante accumule 2 problèmes : le planeur se retrouve rapidement au sol et les turbulences peuvent être dangereuses.

### Exploitation du rotor

De plus, les rotors sont relativement bien organisés et il est donc possible pour un bon pilote de planeur d'exploiter les ascendances thermiques et les ascendances liées au rotor lui-même. Ainsi, on notera que la dimension d'un rotor est de l'ordre de 1 à 2 kilomètres. Toutefois, ce rotor, qui est relativement étendu, est composé de sous-tourbillons dont l'étendue est de l'ordre de quelques dizaines de mètres. Les rotors de type I sont peu dangereux pour le vélivole si celui-ci se déplace suffisamment lentement et reste à l'extérieur du nuage. Au cours de l'ascension à l'intérieur d'un rotor, le pilote devra ouvrir le virage du côté au vent et resserrer le virage du côté sous le vent.
Lorsque le planeur a atteint une altitude suffisante à l'intérieur du rotor, le pilote va se diriger en direction de la montagne pour contacter l'onde proprement dite. Juste avant de contacter l'onde elle-même, le planeur va être soumis aux très forts sous-tourbillons mentionnés plus haut. Le pilote va subir pendant un temps très court (de l'ordre de quelques secondes) une très forte turbulence puis le pilote aura l'impression de voler dans de l'huile avec le variomètre souvent bloqué en montée.
Comme le rotor correspond à une séparation du flot en altitude avec le sol (théorie de l'œil de chat de Paul Queney), le long de la ligne de séparation entre les deux systèmes (onde laminaire et rotor), la température potentielle sera uniforme le long de cette ligne de séparation. Ainsi si l'air est suffisamment humide, un nuage (cumulus fractus) va se former dans la partie supérieure du rotor. Comme il y a libération de chaleur latente dans le rotor, l'air va devenir instable et donc des turbulences extrêmes vont se former à l'intérieur de ce cumulus à l'air faussement inoffensif. Ceci explique pourquoi plusieurs pilotes ont brisé leur machines en volant dans ces soi-disant nuages inoffensifs$^,$$^,$.
L'exploitation d'un rotor migratoire s'apparente à un travail de Sisyphe. En effet ces rotors dérivent lentement dans le sens du vent et disparaissent au bout de quelques minutes. Le pilote aura donc été poussé vers l'arrière et devra remonter face au vent pour contacter un nouveau rotor. Cette stratégie ne sera valable que si la perte d'altitude en remontant face au vent est inférieure au gain d'altitude à l'intérieur du rotor. Comme la vitesse du vent est supérieure à la vitesse de déplacement du rotor, à chaque cycle, le pilote spiralant dans le rotor devra ouvrir le virage face au vent.

## Mesures à l'intérieur des rotors

50 ans après le Sierra Wave Project un projet similaire de mesures à l'intérieur d'un système ondulatoire a été effectué en mars-avril 2004 ayant pour nom « Sierra Rotors Project ». Les résultats des mesures ont été comparés avec des simulations numériques et ces mesures ont montré qu'au premier ressaut, les rotors pouvaient atteindre 6 000 m soit 20 000 pieds et peuvent affecter des aéronefs volant nettement au-dessus de la ligne de crêtes$^,$.

## Effets terrestres et maritimes

Comme dans le cas de l'effet de foehn, les rotors peuvent se rendre au sol lorsque les conditions sont propices. Ils peuvent alors donner des rafales importantes en aval des montages qui causent des dégâts aux infrastructures, comme les bâtiments, ou donnent des vents de travers qui vont renverser des véhicules sur les routes ou des aéronefs dans les aéroports.
Lorsqu'ils se propagent sur un plan d'eau ils peuvent également affecter la navigation. Par exemple, les voiliers de type Gréement bermudien sont vulnérables aux rotors. En aval d'îles de taille modeste mais d'altitude élevée, des ondes orographiques se forment comme au-dessus la Mer Adriatique en Dalmatie ou au-dessus la Mer Égée. Ces ondes engendrent des rotors dans la couche sous-laminaire et donc de la turbulence horizontale et verticale. La maîtrise du voilier devient alors problématique car la poussée du vent ordinaire dans la voile disparaît lorsque la voile devient horizontale. Le voilier peut alors chavirer. Les marins nomment cette condition rafale descendante ce qui est incorrect car ces rafales ne sont pas associées à des orages, ce sont en fait des tourbillons horizontaux.