Corrigé détaillé ENS Physique Chimie L PC 2022

Les forces exercées par le fluide environnant sur les parois délimitant le volume élémentaire \mathrm{d}^3\tau = \mathrm{d}r \times r\mathrm{d}\theta \times \mathrm{d}z sont des actions de contact surfaciques compressives. Pour une surface élémentaire d'aire \mathrm{d}S orientée par sa normale unitaire extérieure \vec{n}_{\mathrm{ext}}, la force pressante exercée sur \mathrm{d}^3\tau s'écrit :

Les flèches représentant l'action de la pression pointent donc vers l'intérieur du volume élémentaire.

En faisant abstraction des faces supérieure et inférieure de normales \pm\vec{e}_z, on dénombre quatre surfaces latérales :

1. Face intérieure en r (aire \mathrm{d}S_r = r\,\mathrm{d}\theta\,\mathrm{d}z) : de normale sortante -\vec{e}_r, la pression y vaut p(r, \theta, z) et l'action s'exerce selon +\vec{e}_r.
2. Face extérieure en r+\mathrm{d}r (aire \mathrm{d}S_{r+\mathrm{d}r} = (r+\mathrm{d}r)\,\mathrm{d}\theta\,\mathrm{d}z) : de normale sortante +\vec{e}_r, la pression y vaut p(r+\mathrm{d}r, \theta, z) et l'action s'exerce selon -\vec{e}_r.
3. Face latérale en \theta (aire \mathrm{d}S_\theta = \mathrm{d}r\,\mathrm{d}z) : de normale sortante -\vec{e}_\theta, la pression y vaut p(r, \theta, z) et l'action s'exerce selon +\vec{e}_\theta.
4. Face latérale en \theta+\mathrm{d}\theta (aire \mathrm{d}S_{\theta+\mathrm{d}\theta} = \mathrm{d}r\,\mathrm{d}z) : de normale sortante +\vec{e}_\theta, la pression y vaut p(r, \theta+\mathrm{d}\theta, z) et l'action s'exerce selon -\vec{e}_\theta.

Le système présentant une symétrie de révolution autour de l'axe (Oz), la pression est invariante par rotation autour de cet axe. On a donc :

Les arguments de la pression se réduisent ainsi à (r, z) et (r+\mathrm{d}r, z).

Considérons le système fermé constitué par la particule fluide occupant le volume élémentaire \mathrm{d}^3\tau = r\,\mathrm{d}r\,\mathrm{d}\theta\,\mathrm{d}z centré en \mathrm{M}(r,\theta,z) à l'instant t. Sa masse est :

L'étude est menée dans le référentiel du laboratoire \mathcal{R}_0, supposé galiléen.

Cinématique : Le liquide tourne en bloc à la vitesse angulaire uniforme \vec{\Omega} = \Omega\,\vec{e}_z. Chaque particule fluide décrit un cercle de rayon r à la vitesse angulaire \Omega constante ; son vecteur accélération dans \mathcal{R}_0 est purement centripète :

Bilan des forces extérieures :

• Le poids de l'élément de fluide :

  \mathrm{d}\vec{P} = \mathrm{d}m\,\vec{g} = \rho\,\vec{g}\,\mathrm{d}^3\tau

• La résultante des forces de pression exercées par le fluide environnant sur la frontière fermée \partial(\mathrm{d}^3\tau) de la particule :

  \mathrm{d}\vec{F}_p = \oiint_{\partial(\mathrm{d}^3\tau)} -p\,\mathrm{d}^2\vec{S} = -\vec{\nabla}p(\mathrm{M})\,\mathrm{d}^3\tau
• Les forces visqueuses : le fluide étant entraîné en rotation de bloc à la vitesse \vec{v} = r\Omega\,\vec{e}_\theta, il n'y a aucun cisaillement relatif entre couches fluides adjacentes (\partial v/\partial r - v/r = \Omega - \Omega = 0). Les forces de viscosité sont donc rigoureusement nulles.

Principe fondamental de la dynamique : Appliqué à la particule fluide dans le référentiel galiléen \mathcal{R}_0 :

En remplaçant par les expressions précédentes :

En simplifiant par le volume élémentaire \mathrm{d}^3\tau \neq 0 et en réorganisant les termes, on obtient :

D'après le résultat de la question Q2, le gradient de pression s'écrit :

Le champ de pression présentant une symétrie de révolution (p = p(r,z)), l'expression du gradient donnée dans le formulaire de l'énoncé conduit au système d'équations aux dérivées partielles :

En intégrant la première relation par rapport à r :

où f est une fonction de la seule variable z. En dérivant cette expression par rapport à z et en identifiant avec la seconde équation, on obtient :

où A \in \mathbb{R} est une constante d'intégration.

L'expression de la pression en tout point du fluide s'écrit donc :

D'après la question Q3, la pression au sein du liquide s'écrit :

En l'absence d'effets de tension superficielle, la surface libre est à la pression atmosphérique uniforme p_0. Tout point de cette surface de coordonnées (r, \theta, Z(r)) vérifie donc :

En particulier au centre de la cuve, pour r = 0, l'altitude vaut Z(0) = Z_0, ce qui donne :

L'expression de la pression sur la surface libre devient ainsi :

En simplifiant par la masse volumique \rho > 0, on obtient :

d'où l'équation de la surface libre :

Considérons le liquide contenu dans le récipient comme système fermé dans le référentiel galiléen \mathcal{R}_0. En régime permanent, le centre de masse du fluide est immobile selon l'axe vertical (Oz), si bien que la résultante des forces extérieures verticales s'annule :

Ces forces verticales comprennent :

• Le poids du liquide : \vec{P} = -mg\,\vec{e}_z = -\rho g \pi R^2 h\,\vec{e}_z (le volume total de liquide incompressible étant V = \pi R^2 h) ;
• L'action de la pression atmosphérique uniforme p_0 sur la surface libre : par projection sur la section horizontale du cylindre, sa résultante vaut \vec{F}_{\text{atm}} = -p_0 \pi R^2\,\vec{e}_z ;
• L'action des parois latérales, purement radiale par symétrie de révolution (composante selon \vec{e}_z nulle) ;

• L'action du fond du récipient (en z=0) sur le liquide :

  \vec{F}_{\text{fond}\to\text{liq}} = \left( \iint_{\text{fond}} p(r, 0)\,\mathrm{d}S \right) \vec{e}_z = \left( \int_0^R p(r, 0) \times 2\pi r\,\mathrm{d}r \right) \vec{e}_z

L'équilibre vertical impose donc :

D'après la question Q3 et la question Q4, la pression à l'altitude z s'écrit p(r, z) = p_0 + \rho g (Z(r) - z), d'où sur le fond en z=0 :

En injectant cette expression dans la relation d'équilibre :

En remplaçant Z(r) = Z_0 + \dfrac{\Omega^2}{2g}r^2 (établi à la question Q4) :

On en déduit l'altitude au centre :

L'altitude de la surface libre s'écrit alors :

Cette équation est bien de la forme Z(r) = h + k(2r^2 - R^2) avec la constante réelle positive :

Étude et représentation graphique de Q(s) :

En introduisant la variable adimensionnée s = r/R \in [0, 1], la fonction s'écrit :

Ses particularités sont :

• Q est une fonction parabolique croissante de s sur [0, 1] ;
• Tangente horizontale au centre : Q'(0) = 0 avec un minimum en s=0 valant Q(0) = -kR ;
• Annulation en s = 1/\sqrt{2} \simeq 0{,}71 : à cette distance de l'axe, la hauteur du liquide est inchangée par rapport au repos (Z = h) ;
• Maximum atteint sur la paroi en s=1, avec Q(1) = +kR ;
• Symétrie des écarts : l'abaissement au centre |Z(0) - h| = kR^2 est strictement égal à la surélévation au bord Z(R) - h = kR^2, ce qui traduit la conservation du volume de liquide incompressible.

D'après l'équation (3) établie à la question Q5, la hauteur de la surface libre est donnée par :

La fonction Z(r) est strictement croissante sur [0, R]. Elle admet donc :

• son minimum au centre (r = 0) :

  Z(0) = h - k R^2

• son maximum à la périphérie (r = R) :

  Z(R) = h + k R^2

L'hypothèse \mathcal{H}_\Omega impose deux conditions :

1. Le fond reste entièrement immergé : le niveau minimal du liquide ne doit pas atteindre le fond situé à l'altitude z = 0, soit :

   Z(0) \ge 0 \iff h - k R^2 \ge 0 \iff k R^2 \le h

2. Le liquide n'est pas éjecté hors du cylindre : le niveau maximal du liquide ne doit pas dépasser le bord supérieur du cylindre situé à l'altitude z = H, soit :

   Z(R) \le H \iff h + k R^2 \le H \iff k R^2 \le H - h

Dans le cas particulier où h = \dfrac{H}{2}, on a :

Les deux conditions sont donc rigoureusement identiques et s'écrivent :

En remplaçant k = \dfrac{\Omega^2}{4g}, il vient :

Puisque la vitesse angulaire \Omega est positive, la condition sur \Omega traduisant l'hypothèse \mathcal{H}_\Omega est :

Dans le référentiel tournant lié au récipient, l'élément de fluide de volume \mathrm{d}^3\tau et de masse \mathrm{d}m = \rho\,\mathrm{d}^3\tau est soumis à deux forces volumiques :

• le poids : \mathrm{d}\vec{P} = \mathrm{d}m_g\,\vec{g} = \rho\,\mathrm{d}^3\tau\,\vec{g} ;
• la force d'inertie d'entraînement (force centrifuge) : \mathrm{d}\vec{F}_{\text{ie}} = \mathrm{d}m_i\,r\Omega^2\vec{e}_r = \rho\,\mathrm{d}^3\tau\,r\Omega^2\vec{e}_r.

La condition d'équilibre hydrostatique impose que le gradient de pression équilibre ces forces volumiques :

La surface libre étant une isobare (p = p_0), elle est orthogonale à chaque point au gradient de pression, et donc au champ de pesanteur effective \vec{g}_{\text{eff}} = -g\vec{e}_z + r\Omega^2\vec{e}_r.

La raison fondamentale pour laquelle la constante k = \dfrac{\Omega^2}{4g}, et donc la géométrie de la surface libre, est indépendante de \rho repose sur le principe d'équivalence entre masse inerte et masse pesante (m_i = m_g) :

Plaçons-nous dans le référentiel tournant \mathcal{R}_1, en rotation uniforme à la vitesse angulaire \vec{\Omega} = \Omega\,\vec{e}_z par rapport au laboratoire.

Dans ce référentiel, la pesanteur apparente (ou effective) subie par un point matériel de masse m s'écrit :

D'après la question Q2, le gradient de pression dans le fluide en rotation s'écrit précisément :

La surface libre du fluide étant une isobare (p = p_0), elle est en tout point orthogonale au gradient de pression, donc orthogonale au vecteur \vec{g}_{\text{eff}}.

Considérons le flotteur sphérique de masse m_f, de volume total V et de volume immergé V_{\text{imm}}. Le gradient de pression étant supposé uniforme à l'échelle du flotteur, la résultante des forces de pression exercées par le liquide sur sa partie immergée (poussée d'Archimède généralisée) s'écrit :

Le flotteur est également soumis à son poids effectif (poids réel et force d'inertie d'entraînement centrifuge) :

La résultante des forces agissant sur le flotteur dans le référentiel tournant est donc :

1. Position d'équilibre privilégiée sur la surface :

Puisque le vecteur \vec{g}_{\text{eff}} est rigoureusement normal à la surface libre en tout point, la résultante \vec{F} est purement normale à la surface libre. La composante des forces tangentielle à la surface libre est donc identiquement nulle :

Le flotteur ne subit aucune force le sollicitant vers l'axe de rotation ou vers la paroi.

2. Modification de la situation de flottaison :

À l'équilibre selon la direction normale à la surface libre, on a :

• Volume immergé : Le volume immergé V_{\text{imm}} (et donc la fraction volumique immergée V_{\text{imm}}/V) ne dépend ni de \Omega ni de la position r. Il est strictement le même qu'au repos.
• Plan de flottaison : Le plan de flottaison reste tangent à la surface libre locale ; il est donc incliné d'un angle \alpha(r) par rapport à l'horizontale, avec \tan\alpha = \frac{r\Omega^2}{g}.
• Intensité des forces : L'intensité du poids effectif et de la poussée d'Archimède est multipliée par le facteur \sqrt{1 + \frac{r^2\Omega^4}{g^2}} \ge 1.

L'équation de la surface libre, avec pour origine le sommet \mathrm{S}, est donnée par :

Soit \alpha l'angle que fait la tangente au profil méridien du miroir avec l'horizontale en un point M(r, W(r)). La pente s'écrit :

Un rayon incident se propageant selon -\vec{e}_z fait un angle d'incidence i = \alpha avec la normale locale à la surface. D'après les lois de Snell-Descartes pour la réflexion, l'angle de réflexion vaut également i = \alpha. Le rayon réfléchi est ainsi dévié d'un angle :

par rapport à la verticale descendante, en direction de l'axe optique (\mathrm{S}z).

Suivant l'indication de l'énoncé, choisissons le point M particulier tel que :

Pour cette inclinaison, la déviation vaut D = 90^\circ : le rayon réfléchi émerge horizontalement (selon -\vec{e}_r). Il coupe donc l'axe (\mathrm{S}z) au foyer \mathrm{F} exactement à la même altitude que le point de réflexion M :

La condition \tan\alpha = 4kr = 1 donne la coordonnée radiale de ce point :

L'altitude correspondante vaut alors :

En remplaçant k = \dfrac{\Omega^2}{4g}, on obtient finalement :

D'après le résultat de la question 9, la distance focale f du télescope parabolique est reliée à la vitesse angulaire \Omega par :

On en déduit l'expression de la vitesse angulaire :

Avec g = 10\text{ m}\cdot\text{s}^{-2} et f = 10\text{ m}, on obtient :

Pour exprimer cette vitesse angulaire en tours par minute, notée usuellement N, on utilise la relation de conversion :

En prenant \pi \simeq 3{,}14 et \sqrt{2} \simeq 1{,}41, le dénominateur vaut \pi\sqrt{2} \simeq 4{,}44, d'où :

1. Nature des objets et formation des images

Les étoiles sont situées à des distances astronomiques très grandes devant les dimensions du télescope : elles constituent des objets à l'infini.

• Pour une étoile \mathrm{E}_1 située rigoureusement sur l'axe optique (Sz) du télescope (direction de visée au zénith), les rayons lumineux incidents sont parallèles à l'axe (Sz) (direction -\vec{e}_z). D'après la propriété du réflecteur parabolique démontrée à la question Q9, tous ces rayons se réfléchissent en convergeant exactement au foyer principal image F :

  \mathrm{E}'_1 = F

• Si l'étoile \mathrm{E}_2 est également située strictement sur le même axe de visée, ses rayons incidents sont également parallèles à l'axe optique. Son image se forme donc au même foyer principal :

  \mathrm{E}'_2 = \mathrm{E}'_1 = F

  Quelle que soit la distance de ces étoiles à la Terre, leurs images se superposent en F.

• Si l'on interprète « dans l'axe de visée » au sens de deux étoiles distinctes résolubles présentes dans le champ de visée du télescope : \mathrm{E}_1 est sur l'axe optique et converge en F, tandis que \mathrm{E}_2 présente une direction inclinée d'un angle \alpha par rapport à l'axe. Les rayons incidents issus de \mathrm{E}_2, parallèles entre eux et inclinés de \alpha, convergent après réflexion en un foyer secondaire \mathrm{E}'_2 situé dans le plan focal image.

2. Position du détecteur CCD

Pour capter des images nettes d'objets situés à l'infini, la surface sensible du détecteur CCD doit être placée dans le plan focal image du miroir :

L'utilisation d'un réflecteur parabolique présente un avantage déterminant par rapport à un miroir sphérique :

• Stigmatisme rigoureux pour un objet à l'infini sur l'axe : par propriété géométrique de la parabole, tout rayon incident parallèle à l'axe optique converge exactement au foyer image F, quelle que soit sa distance à l'axe.
• Absence d'aberration de sphéricité : un miroir sphérique ne réalise un stigmatisme qu'approché, valable uniquement dans les conditions de Gauss (rayons paraxiaux, proches du sommet et peu inclinés). Pour un miroir sphérique de grande dimension, les rayons marginaux (éloignés de l'axe) se croisent plus près du sommet que les rayons paraxiaux, formant une tache focale le long de l'axe appelée aberration géométrique de sphéricité.

Le réflecteur parabolique permet donc de s'affranchir de l'aberration de sphéricité et d'ouvrir largement le miroir (grand diamètre) afin de collecter un maximum de lumière tout en conservant une image parfaitement nette au foyer pour les objets situés sur l'axe.

L'utilisation d'un télescope de grand diamètre D présente deux intérêts physiques majeurs :

1. Augmentation du pouvoir collecteur (luminosité) : La surface collectrice du miroir primaire est proportionnelle au carré de son diamètre :

   S = \frac{\pi D^2}{4}

   Le flux lumineux collecté en provenance d'un objet céleste lointain est donc proportionnel à D^2. Un diamètre plus grand permet de capter des flux très faibles, rendant ainsi observables des objets très peu lumineux (étoiles de magnitude élevée, galaxies lointaines) ou de diminuer les temps d'exposition nécessaires pour un rapport signal sur bruit donné.

2. Amélioration du pouvoir de résolution angulaire (séparation) : L'image d'un point à l'infini à travers une pupille circulaire de diamètre D est limitée par la diffraction sous la forme d'une tache d'Airy dont le demi-angle d'ouverture est donné par le critère de Rayleigh :

   \theta_{\min} \simeq 1{,}22\,\frac{\lambda}{D}

   Augmenter le diamètre D diminue la diffraction et permet ainsi de séparer des sources angulairement très proches ou de distinguer des détails plus fins de l'objet observé.

L'hypothèse \mathcal{H}_v rassemble les hypothèses physiques suivantes :

1. Écoulement orthoradial (ou azimutal) pur : les composantes radiale et axiale de la vitesse sont rigoureusement nulles :

   v_r = 0 \quad \text{et} \quad v_z = 0

   Le mouvement des particules fluides s'effectue donc uniquement selon des trajectoires circulaires concentriques coaxiales à l'axe (\mathrm{O}z) (absence de recirculation méridienne ou de courant secondaire).

2. Symétrie de révolution (axisymétrie) : le champ de vitesse est invariant par rotation autour de l'axe (\mathrm{O}z), il ne dépend donc pas de la coordonnée angulaire \theta :

   \frac{\partial \vec{v}}{\partial \theta} = \vec{0}

3. Écoulement invariant selon l'axe vertical : la vitesse ne dépend pas de l'altitude z :

   \frac{\partial \vec{v}}{\partial z} = \vec{0}

   On néglige ainsi tout effet de couche limite et de frottement visqueux au fond du récipient (z = 0), ainsi qu'à la surface libre.

Considérons une particule fluide de volume élémentaire \mathrm{d}^3\tau = r\,\mathrm{d}r\,\mathrm{d}\theta\,\mathrm{d}z et de masse \mathrm{d}m = \rho\,\mathrm{d}^3\tau, animée de la vitesse \vec{v}(M, t) = v(r, t)\,\vec{e}_\theta(\theta). Son énergie cinétique élémentaire dans le référentiel galiléen \mathcal{R}_0 s'écrit :

L'énergie cinétique totale E_{\mathrm{c}}(t) du liquide s'obtient par sommation sur l'ensemble du volume \mathcal{V}(t) occupé par le fluide à l'instant t, délimité par le fond du récipient en z = 0, la surface libre en z = Z(r, t), la paroi en r = R et la symétrie axiale sur \theta \in [0, 2\pi] :

Comme l'intégrande est indépendant de \theta et de z, les intégrations selon ces deux variables sont immédiates :

On en déduit l'expression de l'énergie cinétique sous forme d'une intégrale simple sur le rayon r :

Dans l'état initial, le liquide tourne en bloc avec le récipient à la vitesse angulaire uniforme \vec{\Omega} = \Omega \vec{e}_z. Le champ de vitesse s'écrit donc :

1. Énergie cinétique E_{\mathrm{c}}^{\star} pour une surface plane et horizontale :

Si la surface libre restait plane et horizontale à l'altitude Z(r) = h, l'expression établie à la question Q15 donnerait :

La masse totale du liquide incompressible étant m = \rho \pi R^2 h, on en déduit :

2. Énergie cinétique réelle E_{\mathrm{c}0} :

En tenant compte du profil parabolique de la surface libre Z(r) = h + k(2r^2 - R^2) issu de la question Q5, l'énergie cinétique s'écrit :

En identifiant avec la forme E_{\mathrm{c}0} = E_{\mathrm{c}}^{\star}(1 + \alpha_{\mathrm{c}}), on obtient directement :

Soit un élément de volume de liquide \mathrm{d}^3\tau = r\,\mathrm{d}r\,\mathrm{d}\theta\,\mathrm{d}z de masse \mathrm{d}m = \rho\,\mathrm{d}^3\tau, situé à l'altitude z. Compte tenu du champ de pesanteur uniforme \vec{g} = -g\,\vec{e}_z et de la référence choisie pour l'énergie potentielle (E_{\mathrm{p}} = 0 en z = 0), l'énergie potentielle gravitationnelle élémentaire s'écrit :

L'énergie potentielle gravitationnelle totale du liquide à l'instant t s'obtient en intégrant sur le volume \mathcal{V}(t) délimité par r \in [0, R], \theta \in [0, 2\pi] et z \in [0, Z(r, t)] :

On obtient ainsi :

D'après la question Q17, l'énergie potentielle gravitationnelle du liquide s'écrit sous forme intégrale :

Énergie potentielle de référence E_{\mathrm{p}}^\star : Si la surface libre restait plane et horizontale, on aurait Z(r) = h pour tout r \in [0, R]. L'énergie potentielle associée s'écrirait alors :

La masse totale du liquide s'exprimant par m = \rho \pi R^2 h, on obtient immédiatement :

Énergie potentielle initiale réelle E_{\mathrm{p}0} : Dans l'état initial, la surface libre parabolique est décrite par l'équation (3) (question Q5) :

Développons son carré :

Calculons l'intégrale \int_0^R r Z^2(r) \, \mathrm{d}r terme à terme :

• Premier terme :

  \int_0^R r h^2 \, \mathrm{d}r = \frac{1}{2} h^2 R^2

• Terme croisé (nul en vertu de la conservation du volume du liquide incompressible) :

  \int_0^R 2 h k (2r^3 - r R^2) \, \mathrm{d}r = 2 h k \left[ \frac{r^4}{2} - \frac{r^2 R^2}{2} \right]_0^R = 2 h k \left( \frac{R^4}{2} - \frac{R^4}{2} \right) = 0

• Dernier terme :

  \begin{aligned} \int_0^R r k^2 (4r^4 - 4r^2 R^2 + R^4) \, \mathrm{d}r &= k^2 \int_0^R (4r^5 - 4r^3 R^2 + r R^4) \, \mathrm{d}r \\ &= k^2 \left[ \frac{2}{3} r^6 - r^4 R^2 + \frac{1}{2} r^2 R^4 \right]_0^R \\ &= k^2 R^6 \left( \frac{2}{3} - 1 + \frac{1}{2} \right) = \frac{1}{6} k^2 R^6 \end{aligned}

On en déduit l'expression de l'intégrale :

En multipliant par \pi \rho g, il vient :

avec :

Dans l'état final (t \to \infty), le liquide s'est entièrement immobilisé :

Par conservation de la masse d'un fluide incompressible, la surface libre retrouve sa forme plane et horizontale d'altitude uniforme h. D'après la définition de E_{\mathrm{p}}^\star donnée à la question Q18 :

L'énergie mécanique finale s'écrit donc :

Dans l'état initial (t = 0), d'après les résultats des questions Q16 et Q18 :

La variation d'énergie mécanique s'en déduit :

En identifiant avec la forme proposée \Delta E_{\mathrm{m}} = -E_{\mathrm{c}}^\star (1 + \alpha), on a :

Exprimons le rapport E_{\mathrm{p}}^\star / E_{\mathrm{c}}^\star à l'aide des définitions :

En utilisant \alpha_{\mathrm{c}} = \dfrac{k R^2}{3h} et \alpha_{\mathrm{p}} = \dfrac{k^2 R^4}{3h^2}, le second terme devient :

Rappelons la relation établie à la question Q5 : k = \dfrac{\Omega^2}{4g}, soit \dfrac{g}{\Omega^2} = \dfrac{1}{4k}. Il vient alors :

On en déduit l'expression de \alpha en fonction de k :

En remplaçant enfin k par son expression k = \dfrac{\Omega^2}{4g} :

On a bien la forme requise avec :

Considérons le système fermé constitué par le liquide de masse m entre l'état initial (à t = 0^+, liquide en mouvement) et l'état final (t \to \infty, liquide totalement immobilisé).

Le premier principe de la thermodynamique appliqué à un système en mouvement macroscopique dans le référentiel galiléen \mathcal{R}_0 s'écrit :

où \Delta E_{\mathrm{m}} = \Delta E_{\mathrm{c}} + \Delta E_{\mathrm{p}} est la variation de l'énergie mécanique macroscopique (analysée à la question Q19), \Delta U la variation de son énergie interne, Q le transfert thermique reçu et W le travail des forces extérieures autres que la pesanteur (cette dernière étant prise en compte dans E_{\mathrm{p}}).

Détaillons les échanges énergétiques avec l'extérieur au cours de la phase d'arrêt :

• Transfert thermique Q : l'énoncé stipule que tout échange de chaleur entre le liquide et le milieu extérieur (parois de la cuve et air ambiant) est négligeable, d'où :

  Q = 0

• Travail des forces de pression atmosphérique W_{\text{atm}} : la pression atmosphérique p_0 étant supposée uniforme et constante, le travail exercé sur la surface libre s'écrit :

  W_{\text{atm}} = -\int_{V_{\text{i}}}^{V_{\text{f}}} p_0 \,\mathrm{d}V = -p_0 \Delta V

  Le liquide étant incompressible, son volume global demeure rigoureusement constant (\Delta V = 0), d'où W_{\text{atm}} = 0.

• Travail des forces exercées par les parois de la cuve W_{\text{cuve}} : la cuve est maintenue immobile pour tout t > 0. Les points de la paroi ont une vitesse nulle dans \mathcal{R}_0. Par condition d'adhérence du fluide visqueux sur les parois fixes, la vitesse du fluide au contact des parois est nulle. La puissance des forces de contact (pression et frottement visqueux) exercées par la cuve sur le liquide est donc identiquement nulle à tout instant :

  \mathcal{P}_{\text{cuve}\to\text{liq}} = 0 \implies W_{\text{cuve}} = 0

Le travail total reçu par le liquide est donc nul :

On en déduit le bilan énergétique global :

En injectant le résultat de la question Q19, il vient :

Le liquide étant assimilé à une phase condensée indilatable et incompressible, sa capacité thermique massique à pression constante se confond avec celle à volume constant. La variation d'énergie interne entre l'état initial et l'état final s'écrit donc :

D'après le bilan énergétique établi à la question Q20, on a :

En utilisant l'expression de E_{\mathrm{c}}^{\star} = \frac{1}{4} m R^2 \Omega^2 (question Q16) et celle de \alpha = \frac{1}{8} \frac{R^2\Omega^2}{gh} (question Q19), il vient :

En simplifiant par la masse m, on obtient :

Application numérique :

• R = 4\text{ cm} = 0{,}04\text{ m} ;
• h = 4\text{ cm} = 0{,}04\text{ m} ;
• g = 10\text{ m}\cdot\text{s}^{-2} ;
• \Omega = 6\text{ rad}\cdot\text{s}^{-1} ;
• c_p = 4\times 10^3\text{ J}\cdot\text{kg}^{-1}\cdot\text{K}^{-1}.

On calcule le terme correctif \alpha :

Ainsi, 1 + \alpha \simeq 1{,}02 \simeq 1. La variation de température vaut alors :

Commentaire : La hausse de température est extrêmement faible (de l'ordre de quelques microkelvins) et totalement indétectable avec un thermomètre usuel. Cela s'explique par la très forte capacité thermique de l'eau (c_p = 4\times 10^3\text{ J}\cdot\text{kg}^{-1}\cdot\text{K}^{-1}) mise en regard de l'énergie cinétique massique macroscopique très modeste mise en jeu (\frac{1}{4}R^2\Omega^2 \sim 10^{-2}\text{ J}\cdot\text{kg}^{-1}). L'hypothèse d'une transformation isotherme est donc rigoureusement justifiée.

Considérons le liquide comme système fermé entre :

• l'état initial : rotation en bloc à la vitesse angulaire \Omega,
• l'état final : liquide totalement au repos dans \mathcal{R}_0.

Ces deux états d'équilibre mécanique étant rigoureusement identiques à ceux de la question précédente, la variation d'énergie mécanique \Delta E_{\mathrm{m}} < 0 est strictement la même dans les deux cas.

Appliquons le premier principe de la thermodynamique dans le référentiel du laboratoire \mathcal{R}_0 :

L'enceinte étant toujours supposée athermane (Q = 0), on a :

Distinguons les deux situations pour le travail mécanique W reçu par le liquide :

• Immobilisation brutale : les parois de la cuve sont immobiles dans \mathcal{R}_0 pour tout t > 0. Les forces surfaciques exercées par la paroi sur le fluide s'appliquent en des points de vitesse nulle dans \mathcal{R}_0, de sorte que la puissance reçue est nulle. Le travail reçu par le fluide est donc nul :

  W_{\text{brutal}} = 0 \implies \Delta U_{\text{brutal}} = -\Delta E_{\mathrm{m}}

  Toute l'énergie mécanique initiale perdue est alors convertie en énergie interne par dissipation visqueuse interne.

• Ralentissement progressif : le récipient est en rotation à une vitesse angulaire \Omega_c(t) > 0 qui décroît jusqu'à zéro. Le liquide tournant plus vite que le récipient en cours de décélération, la paroi exerce sur le liquide un couple résistant de freinage :

  \vec{\Gamma}_{\text{paroi}\to\text{liq}} = \Gamma_z \vec{e}_z \quad \text{avec} \quad \Gamma_z < 0

  La puissance mécanique reçue par le liquide dans \mathcal{R}_0 s'écrit :

  \mathcal{P}_{\text{reçue}} = \vec{\Gamma}_{\text{paroi}\to\text{liq}} \cdot \vec{\Omega}_c(t) = \Gamma_z \Omega_c(t) < 0

  Le travail mécanique reçu par le liquide est donc strictement négatif :

  W_{\text{progressif}} = \int \mathcal{P}_{\text{reçue}}\,\mathrm{d}t < 0

  Le fluide fournit un travail mécanique au récipient (qui le transmet à son système de freinage extérieur).

Il s'ensuit que :

Comme \Delta U = m c_p \Delta T, on en déduit que :

La variation de température du liquide aurait donc été plus faible.

Dans le cas d'une rotation en bloc du fluide à vitesse angulaire uniforme \Omega (mouvement de solide indéformable), le champ de vitesse s'écrit :

Lors d'un tel mouvement, la distance entre deux particules fluides quelconques reste invariante au cours du temps : il n'y a aucune déformation relative du milieu, et par conséquent les forces de frottement visqueux doivent être strictement nulles, soit \tau = 0.

• Si l'on omettait le terme en v/r, la contrainte s'écrirait \tau = \eta \frac{\partial v}{\partial r} = \eta \Omega \neq 0, ce qui prédirait une dissipation visqueuse lors d'un simple mouvement de corps rigide, ce qui est absurde.

• Avec l'expression (9), on constate que :

  \begin{aligned} \frac{\partial v}{\partial r} - \frac{v}{r} &= \Omega - \frac{\Omega\,r}{r} = 0 \implies \tau = 0 \end{aligned}

La présence du terme -v/r est donc indispensable pour garantir l'absence de contrainte de cisaillement visqueuse lors d'une rotation globale de solide indéformable.

Considérons le domaine annulaire délimité par les rayons r et r+\mathrm{d}r, et les cotes z et z+\mathrm{d}z. Il possède quatre surfaces frontières : deux surfaces horizontales en z et z+\mathrm{d}z, et deux surfaces cylindriques en r et r+\mathrm{d}r.

1. Inventaire des forces agissant sur le domaine annulaire :

• Action à distance : le poids de l'élément de volume,

  \mathrm{d}^2\vec{P} = -\rho\,\mathrm{d}^2V\,g\,\vec{e}_z = -2\pi\rho g r\,\mathrm{d}r\,\mathrm{d}z\,\vec{e}_z

• Forces de pression (actions normales) :

  • Sur les bases horizontales : forces dirigées selon \pm\vec{e}_z ;
  • Sur les faces cylindriques intérieure (r) et extérieure (r+\mathrm{d}r) : forces normales dirigées radialement selon \pm\vec{e}_r.

• Forces de viscosité (actions tangentielles de cisaillement) :

  • Sur les surfaces supérieure et inférieure : d'après les hypothèses du sujet, les contraintes de cisaillement de normale \pm\vec{e}_z sont nulles ;
  • Sur la surface cylindrique interne en r : force de cisaillement tangentielle exercée par le fluide intérieur (r' < r), colinéaire à -\vec{e}_\theta ;
  • Sur la surface cylindrique externe en r+\mathrm{d}r : force de cisaillement tangentielle exercée par le fluide extérieur (r' > r+\mathrm{d}r), colinéaire à +\vec{e}_\theta.

2. Identification des composantes de moment selon l'axe (\mathrm{O}z) :

Le moment élémentaire en \mathrm{O} d'une force \vec{F} appliquée en \mathrm{M} = r\vec{e}_r + z\vec{e}_z a pour projection sur l'axe (\mathrm{O}z) :

• Poids et forces de pression verticales : ces forces sont colinéaires à \vec{e}_z, d'où \overrightarrow{\mathrm{OM}} \wedge \vec{e}_z = (r\vec{e}_r + z\vec{e}_z) \wedge \vec{e}_z = -r\,\vec{e}_\theta, qui est orthogonal à \vec{e}_z. Leur moment selon (\mathrm{O}z) est strictement nul.
• Forces de pression radiales : ces forces sont colinéaires à \vec{e}_r et coupent l'axe (\mathrm{O}z). En effet, \overrightarrow{\mathrm{OM}} \wedge \vec{e}_r = (r\vec{e}_r + z\vec{e}_z) \wedge \vec{e}_r = -z\,\vec{e}_\theta, d'où une composante selon (\mathrm{O}z) nulle.

• Forces de cisaillement visqueux : elles sont orthoradiales (selon \vec{e}_\theta). Or :

  \overrightarrow{\mathrm{OM}} \wedge \vec{e}_\theta = (r\vec{e}_r + z\vec{e}_z) \wedge \vec{e}_\theta = r\,\vec{e}_z - z\,\vec{e}_r

  La projection sur \vec{e}_z vaut donc r \neq 0.

1. Moment cinétique du domaine élémentaire :

Le domaine annulaire élémentaire a pour volume \mathrm{d}^2V = 2\pi r\,\mathrm{d}r\,\mathrm{d}z et pour masse :

Tous les points de cet anneau sont situés à la même distance r de l'axe (\mathrm{O}z) et sont animés de la vitesse orthoradiale \vec{v}(r, t) = v(r, t)\,\vec{e}_\theta.

Le moment cinétique élémentaire en \mathrm{O} d'une particule de masse \mathrm{d}^3m située en \overrightarrow{\mathrm{OM}} = r\vec{e}_r + z\vec{e}_z vaut :

En intégrant sur \theta \in [0, 2\pi], le terme en \vec{e}_r(\theta) s'annule par symétrie. La composante selon l'axe (\mathrm{O}z) du moment cinétique du domaine annulaire s'écrit donc :

soit :

2. Théorème du moment cinétique et équation aux dérivées partielles :

D'après le bilan établi à la question précédente, seules les contraintes de cisaillement agissant sur les parois cylindriques de rayons r et r+\mathrm{d}r créent un moment non nul selon (\mathrm{O}z) :

• En r+\mathrm{d}r : le fluide extérieur (r' > r+\mathrm{d}r) exerce sur le domaine annulaire une contrainte \vec{\tau}(r+\mathrm{d}r, t) = \tau(r+\mathrm{d}r, t)\,\vec{e}_\theta à travers l'aire \mathrm{d}^2S_e = 2\pi(r+\mathrm{d}r)\,\mathrm{d}z. Le moment résultant en \mathrm{O} selon (\mathrm{O}z) est :

  \mathrm{d}\mathcal{M}_{z, e} = (r+\mathrm{d}r) \times \left[\tau(r+\mathrm{d}r, t)\,2\pi(r+\mathrm{d}r)\,\mathrm{d}z\right] = 2\pi (r+\mathrm{d}r)^2 \tau(r+\mathrm{d}r, t)\,\mathrm{d}z

• En r : par la troisième loi de Newton, le fluide intérieur (r' < r) exerce sur le domaine annulaire la force opposée à celle que le domaine exerce sur lui, soit une contrainte -\vec{\tau}(r, t) = -\tau(r, t)\,\vec{e}_\theta à travers l'aire \mathrm{d}^2S_i = 2\pi r\,\mathrm{d}z. Le moment résultant en \mathrm{O} selon (\mathrm{O}z) est :

  \mathrm{d}\mathcal{M}_{z, i} = -2\pi r^2 \tau(r, t)\,\mathrm{d}z

Le moment total selon (\mathrm{O}z) appliqué au système fermé de fluide contenu dans l'anneau vaut ainsi :

D'après le théorème du moment cinétique appliqué à ce système fermé dans le référentiel galiléen \mathcal{R}_0, et en tenant compte de la condition admise \frac{\mathrm{d}\vec{v}}{\mathrm{d}t} = \frac{\partial\vec{v}}{\partial t} :

ce qui donne :

En simplifiant par 2\pi\,\mathrm{d}r\,\mathrm{d}z, on obtient l'équation aux dérivées partielles liant v(r, t) et \tau(r, t) :

D'après la question Q25, le théorème du moment cinétique appliqué au domaine annulaire élémentaire fournit la relation :

En injectant l'expression de la contrainte de cisaillement issue de l'équation (9) :

il vient :

Dérivons cette expression par rapport à r :

En divisant par r^2, on obtient :

Or, d'après l'expression (1) du formulaire, le laplacien scalaire d'une fonction ne dépendant que de r (ici v(r, t) à t fixé) s'écrit :

On en déduit :

En introduisant la viscosité cinématique \nu = \dfrac{\eta}{\rho}, on parvient à l'équation aux dérivées partielles demandée :

1. Condition initiale (t = 0^+) :

Avant l'arrêt (t = 0^-), le liquide est animé d'un mouvement de rotation solide uniforme de vitesse angulaire \Omega, soit un champ de vitesse initial v(r, 0^-) = \Omega r.

Par inertie mécanique, le cœur du fluide ne s'arrête pas instantanément : pour tout point strictement intérieur au cylindre (0 \le r < R), on a :

2. Conditions aux limites spatiales :

• À la paroi latérale (r = R) :
  Le récipient est immobilisé pour tout t > 0. D'après la condition d'adhérence d'un fluide visqueux sur une paroi solide immobile, la vitesse relative du fluide par rapport à la paroi est nulle :

  \boxed{v(R, t) = 0 \quad \forall t > 0}

• Sur l'axe de symétrie (r = 0) :
  Le vecteur unitaire orthoradial \vec{e}_\theta n'étant pas défini de manière univoque en r = 0, la continuité et la régularité du champ de vitesse \vec{v}(0, t) = v(0, t)\vec{e}_\theta au centre imposent que la vitesse s'y annule :

  \boxed{v(0, t) = 0 \quad \forall t \ge 0}

Remarque asymptotique :
Lorsque t \to +\infty, les frottements visqueux dissipent toute l'énergie cinétique du liquide, conduisant à l'arrêt complet :

D'après la relation (11) de l'énoncé, la vorticité \omega(r, t) s'écrit :

En multipliant par r, il vient :

En intégrant cette relation par rapport à la variable radiale muette u entre 0 et r, à l'instant t fixé :

Le membre de gauche s'évalue immédiatement :

D'après la question Q27, la vitesse reste finie sur l'axe et s'annule en u = 0 (v(0, t) = 0). Le terme limite en u=0 est donc nul :

On en déduit :

Soit, pour r > 0 :

D'après la question Q28, le champ de vitesse s'exprime en fonction de la vorticité par :

1. Calcul des dérivées partielles de v(r, t) :

• Dérivée par rapport au temps t : la dérivation sous le signe intégral donne immédiatement :

  \frac{\partial v}{\partial t}(r, t) = \frac{1}{r}\int_0^r u\,\frac{\partial \omega}{\partial t}(u, t)\,\mathrm{d}u

• Dérivée première par rapport à r : la fonction I(r, t) = \int_0^r u\,\omega(u, t)\,\mathrm{d}u a pour dérivée spatiale \frac{\partial I}{\partial r} = r\,\omega(r, t). En appliquant la règle de Leibniz :

  \begin{aligned} \frac{\partial v}{\partial r}(r, t) &= -\frac{1}{r^2}\int_0^r u\,\omega(u, t)\,\mathrm{d}u + \frac{1}{r}\left[r\,\omega(r, t)\right] \\ &= -\frac{v(r, t)}{r} + \omega(r, t) \end{aligned}

  soit sous forme intégrale :

  \frac{\partial v}{\partial r}(r, t) = \omega(r, t) - \frac{1}{r^2}\int_0^r u\,\omega(u, t)\,\mathrm{d}u

• Dérivée seconde par rapport à r : en dérivant la forme \frac{\partial v}{\partial r} = \omega - \frac{v}{r} par rapport à r :

  \begin{aligned} \frac{\partial^2 v}{\partial r^2}(r, t) &= \frac{\partial \omega}{\partial r}(r, t) - \frac{1}{r}\frac{\partial v}{\partial r}(r, t) + \frac{v(r, t)}{r^2} \\ &= \frac{\partial \omega}{\partial r}(r, t) - \frac{1}{r}\left(\omega(r, t) - \frac{v(r, t)}{r}\right) + \frac{v(r, t)}{r^2} \\ &= \frac{\partial \omega}{\partial r}(r, t) - \frac{\omega(r, t)}{r} + \frac{2v(r, t)}{r^2} \end{aligned}

  soit sous forme intégrale :

  \frac{\partial^2 v}{\partial r^2}(r, t) = \frac{\partial \omega}{\partial r}(r, t) - \frac{\omega(r, t)}{r} + \frac{2}{r^3}\int_0^r u\,\omega(u, t)\,\mathrm{d}u

2. Équation aux dérivées partielles vérifiée par \omega(r, t) : Le champ de vitesse étant indépendant de \theta et z, le laplacien scalaire s'écrit, conformément au formulaire (1) de l'énoncé :

On calcule alors le terme intervenant dans l'équation (10) :

L'équation de Navier-Stokes projetée sur \vec{e}_\theta (équation (10) établie à la question Q26) s'écrit donc simplement :

En appliquant l'opérateur \frac{1}{r}\frac{\partial}{\partial r}(r\,\cdot) de part et d'autre :

Puisque l'ordre des dérivations spatiale et temporelle peut être interverti (théorème de Schwarz) :

De plus, la vorticité \omega ne dépendant ni de \theta ni de z, le laplacien de \omega est donné par le formulaire (1) :

On en déduit finalement :

1. Propriétés de l'équation (12) et phénomène décrit

L'équation vérifiée par la vorticité \omega(r, t) :

présente les propriétés remarquables suivantes :

• Linéarité : elle est linéaire et homogène, ce qui permet l'application du principe de superposition.
• Ordres de dérivation : elle est du premier ordre par rapport au temps t et du second ordre par rapport aux coordonnées d'espace.
• Irréversibilité : la dissymétrie entre l'ordre 1 en temps et l'ordre 2 en espace brise l'invariance par renversement du temps (t \to -t) ; il s'agit d'une équation décrivant un phénomène dissipatif.

Cette équation décrit un phénomène de diffusion moléculaire de la vorticité (ou diffusion visqueuse de la quantité de mouvement), gouverné par la viscosité cinématique \nu qui joue le rôle de coefficient de diffusion. Elle rend compte de l'atténuation et de l'étalement spatial de la rotation du fluide vers les parois immobiles sous l'effet des frottements visqueux.

2. Analogie avec un autre domaine de la physique

Cette équation est formellement identique à :

• l'équation de la chaleur régissant le champ de température T en conduction thermique pure (loi de Fourier) :

  \frac{\partial T}{\partial t} = D_{\mathrm{th}} \Delta T \quad \text{avec } D_{\mathrm{th}} = \frac{\lambda}{\rho c}

• ou encore l'équation de diffusion de matière (première et deuxième lois de Fick) :

  \frac{\partial c}{\partial t} = D \Delta c

La situation étudiée est ainsi directement analogue au refroidissement d'un cylindre solide initialement à température uniforme, dont la paroi latérale extérieure serait brutalement maintenue à température nulle.

3. Définition du flux surfacique de vorticité

En exprimant l'équation de diffusion sous la forme d'un bilan local de conservation de la vorticité :

et en constatant que \Delta \omega = \mathrm{div}(\vec{\nabla}\omega), on identifie le vecteur densité de flux de vorticité :

Le champ \omega ne dépendant que de r et t, le gradient s'écrit \vec{\nabla}\omega = \frac{\partial \omega}{\partial r}\vec{e}_r, de sorte que \vec{j}_\omega = \varphi_\omega(r, t)\,\vec{e}_r.

Le flux surfacique algébrique de vorticité à travers une surface cylindrique d'axe (Oz) orientée vers l'extérieur par \vec{e}_r est donc :

Vérification dimensionnelle :

d'où :

ce qui correspond bien à la dimension imposée par l'énoncé.

1. Solution stationnaire générale de l'équation de diffusion de la vorticité

En régime stationnaire, la dérivée temporelle de la vorticité s'annule : \dfrac{\partial \omega}{\partial t} = 0. L'équation (12) se réduit à l'équation de Laplace :

Compte tenu des invariances (symétrie de révolution selon \theta et indépendance en z d'après l'hypothèse \mathcal{H}_v), le laplacien scalaire donné par l'équation (1) de l'énoncé s'écrit :

Pour r > 0, cela conduit par intégrations successives à :

où C_1 est une constante réelle. On en déduit la forme générale de la solution stationnaire :

2. Solution à retenir dans le cadre de notre étude

Le domaine occupé par le liquide comprend l'axe de rotation r = 0. Or, le terme en \ln(r/R) diverge quand r \to 0. Pour des raisons de régularité physique, le champ de vorticité doit rester borné sur l'axe :

La vorticité stationnaire est donc spatialement uniforme :

3. Expression de la vitesse v(r) associée

En utilisant la relation intégrale établie à la question Q28 :

Il vient ainsi :

4. Situations physiques décrites

Le champ de vitesse associé est de la forme v(r) = \Omega^* r avec \Omega^* = \dfrac{\omega_0}{2}, ce qui correspond au profil de vitesse d'une rotation solide (ou rotation en bloc) à la vitesse angulaire constante \Omega^*. Cette solution stationnaire peut décrire :

• L'état initial (t < 0) : le liquide tourne en bloc avec le récipient à la vitesse angulaire \Omega, soit \Omega^* = \Omega. On a alors \omega_0 = 2\Omega et v(r) = \Omega r.
• L'état final (t \to \infty) : le récipient étant immobilisé, la condition d'adhérence à la paroi latérale v(R) = 0 impose \Omega^* = 0, soit \omega_0 = 0 et v(r) = 0. Cela décrit l'arrêt total du liquide au repos.

5. Commentaire

Pour un champ de vitesse linéaire v(r) = \Omega^* r, la contrainte de cisaillement visqueux s'annule identiquement :

Il n'y a donc aucun frottement visqueux interne ni dissipation d'énergie : la rotation solide et le repos complet sont les seuls régimes stationnaires possibles d'un fluide newtonien occupant un cylindre plein. Tout autre profil de vitesse induit un cisaillement qui se dissipe au cours du temps par diffusion visqueuse.

1. Profil radial de vitesse v(r, t)

Pour t = 0^-, le profil est rectiligne :

Pour t > 0, la condition aux limites à la paroi v(R, t) = 0 s'applique immédiatement. Une couche limite d'épaisseur caractéristique \delta(t) \sim \sqrt{\nu t} pénètre dans le fluide depuis la paroi r = R. Aux temps courts (t_1), le cœ ur du fluide conserve une rotation quasi-solide (v \approx \Omega r), puis pour t_2 et t_3, la diffusion de la quantité de mouvement ralentit l'ensemble du volume jusqu'à l'arrêt total (v \to 0 quand t \to \infty).

2. Profil radial de vorticité \omega(r, t) et sens du flux \varphi_\omega

À t = 0^-, le liquide étant en bloc :

Pour t > 0, la paroi freine le fluide, ce qui induit une perte de vorticité en périphérie. Le gradient spatial de vorticité est négatif (\frac{\partial \omega}{\partial r} \le 0) et la tangente en r = 0 est horizontale par parité/symétrie axiale. Le vecteur densité de flux de vorticité s'écrit :

Puisque \frac{\partial \omega}{\partial r} < 0, la composante algébrique est positive : \varphi_\omega > 0. Le flux surfacique de vorticité est donc dirigé vers l'extérieur (du cœ ur du fluide vers la paroi).

Considérons le liquide contenu dans le récipient, occupant un volume caractéristique V \approx \pi R^2 h. Dans le cadre du modèle simplifié proposé, la quantité globale de vorticité Q_\omega(t) s'écrit :

Pendant la durée \mathrm{d}t, la variation de cette quantité de vorticité est due au flux surfacique sortant à travers la frontière du domaine. Les hypothèses de l'énoncé précisent qu'aucune contrainte de cisaillement n'agit sur les surfaces d'altitudes extrêmes (normale \pm \vec{e}_z) ; la diffusion s'effectue donc uniquement à travers la surface latérale du cylindre, d'aire S_{\mathrm{lat}} = 2\pi R h.

Le bilan global de quantité de vorticité s'écrit ainsi :

où \varphi_\omega(R, t) est la composante radiale du flux surfacique de vorticité défini à la question Q30 :

Dans ce modèle très simple à une seule dimension spatiale, le gradient radial à la paroi peut être approché par l'écart caractéristique de vorticité entre la paroi (où le fluide est immobilisé, \omega(R, t) \approx 0) et le centre (\omega(0, t) = \omega_0(t)), sur la distance caractéristique R :

Il vient alors :

En injectant cette relation dans le bilan global :

Soit, en simplifiant par \pi h :

On obtient l'équation différentielle régissant l'évolution de la vorticité :

avec pour temps caractéristique d'immobilisation :

1. Estimation du temps caractéristique \tau

La viscosité cinématique de l'eau vaut :

Avec le rayon du cylindre R = 4\text{ cm} = 4 \times 10^{-2}\text{ m}, le temps caractéristique de diffusion visqueuse radiale déterminé à la question Q33 s'évalue selon :

(ou \tau \sim R^2/\nu = 1600\text{ s} \approx 27\text{ min} en ordre de grandeur dimensionnel).

2. Explications de l'écart avec la valeur expérimentale \tau_{\exp} \sim 1\text{ min}

Le modèle de diffusion purement radiale surestime considérablement la durée de ralentissement : le fluide s'arrête en réalité beaucoup plus rapidement (\tau_{\exp} \ll \tau). Cet écart s'explique par les hypothèses simplificatrices de \mathcal{H}_v (détaillées à la question Q14) :

• Absence de dépendance en z et omission du frottement au fond : le modèle néglige les contraintes de cisaillement sur les surfaces horizontales de normale \pm\vec{e}_z. Or, le fond du récipient est immobile en z = 0, imposant par adhérence \vec{v}(r, \theta, z=0, t) = \vec{0}. Le gradient vertical de vitesse \partial v/\partial z dissipe une part majeure de l'énergie cinétique.

• Écoulement purement orthoradial (v_r = 0, v_z = 0) : à l'arrêt brutal des parois, le gradient de pression radial dans le fluide entraîne le liquide du fond vers le centre car la vitesse orthoradiale y est fortement ralentie (la force centrifuge ne compense plus le gradient de pression). Cela crée un écoulement secondaire méridien (circulation d'Ekman) qui pompe le fluide du cœur vers la couche limite visqueuse du fond. Ce phénomène classique de « spin-down » conduit à un temps caractéristique de freinage :

  \tau_{\text{Ekman}} \sim \frac{h}{\sqrt{\nu \Omega}} \sim \frac{0{,}04}{\sqrt{10^{-6} \times 2\pi}} \sim 16\text{ s} \sim 1\text{ min}

  qui correspond parfaitement à la valeur observée expérimentalement.

• Régime d'écoulement initial : le nombre de Reynolds basé sur l'échelle globale vaut \mathrm{Re} = \frac{\Omega R^2}{\nu} \sim \frac{2\pi \times 16 \times 10^{-4}}{10^{-6}} \sim 10^4. L'écoulement initialement cisaillé peut être turbulent, ce qui augmente le transport effectif de quantité de mouvement par rapport à la seule viscosité moléculaire laminaire.

Le nom officiel est le N-[(4-hydroxy-3-méthoxyphényl)méthyl]nonanamide :

• La partie acide carboxylique donne la chaîne acyle linéaire à 9 carbones :

  \mathrm{CH}_3-(\mathrm{CH}_2)_7-\mathrm{C}(=\mathrm{O})-

  liée à l'azote de la fonction amide -\mathrm{NH}-.

• Le groupe fixé sur l'azote est le groupe (4-hydroxy-3-méthoxyphényl)méthyle, constitué d'un pont méthylène -\mathrm{CH}_2- relié en position 1' d'un cycle phényle substitué par :

  • un groupe hydroxyle -\mathrm{OH} en position 4' (position para) ;
  • un groupe méthoxyle -\mathrm{OCH}_3 en position 3' (position méta).

On en déduit la représentation topologique du nonivamide :

Le nonivamide (N-[(4-hydroxy-3-méthoxyphényl)méthyl]nonanamide) comporte :

• le fragment vanillyle (4-hydroxy-3-méthoxybenzyl) : 8 carbones, 1 cycle aromatique (4 insaturations), 2 atomes d'oxygène ;
• la fonction amide : 1 insaturation, 1 oxygène, 1 azote ;
• la chaîne acyle linéaire (nonanoyle) : 9 carbones, aucune insaturation supplémentaire.

Sa formule brute est donc :

Détaillons les formules brutes des capsaïcinoïdes du tableau 1 :

• Capsaïcine : chaîne acyle en \mathrm{C}_{10} avec une double liaison \mathrm{C=C}, soit \mathrm{C}_{18}\mathrm{H}_{27}\mathrm{NO}_3 ;
• Dihydrocapsaïcine : chaîne acyle en \mathrm{C}_{10} saturée, soit \mathrm{C}_{18}\mathrm{H}_{29}\mathrm{NO}_3 ;
• Nordihydrocapsaïcine : le préfixe nor- indique la perte d'un groupement \mathrm{CH}_2 par rapport à la dihydrocapsaïcine. La chaîne acyle saturée comporte 9 carbones (acide 7-méthyloctanoïque), ce qui conduit à la formule brute \mathrm{C}_{17}\mathrm{H}_{27}\mathrm{NO}_3 ;
• Homodihydrocapsaïcine et Homocapsaïcine : chaînes acyles comportant 11 carbones, soit respectivement \mathrm{C}_{19}\mathrm{H}_{31}\mathrm{NO}_3 et \mathrm{C}_{19}\mathrm{H}_{29}\mathrm{NO}_3.

Le seul composé présentant la même formule brute que le nonivamide est donc la nordihydrocapsaïcine.

Le nonivamide et la nordihydrocapsaïcine possèdent les mêmes fonctions chimiques (phénol, éther, amide secondaire) et le même groupe aromatique, mais diffèrent par la ramification de leur chaîne alkyle (linéaire en \mathrm{C}_9 pour le nonivamide, ramifiée à terminaison isopropyle pour la nordihydrocapsaïcine).

La molécule de nonivamide présente une nature amphiphile :

• une chaîne alkyle grasse en \mathrm{C_8} fortement apolaire et hydrophobe ;
• un groupement amide et une fonction phénol/méthoxy capables d'engager des interactions dipolaires fortes et des liaisons hydrogène (donneur et accepteur).

Une approche qualitative naïve des interactions soluté-solvant prédit qu'un solvant apolaire dissoudra la partie lipophile, tandis qu'un solvant polaire et protique solvatera efficacement le pôle polaire. Or, les données expérimentales du tableau 2 révèlent des incohérences majeures :

1. Solubilités infimes dans deux solvants aux antipodes : l'eau (solvant le plus polaire et fortement protique) et le n-hexane (solvant totalement apolaire, aprotique) présentent tous deux une solubilité extrêmement faible et du même ordre de grandeur :

   x_{\mathrm{sol}}(\text{eau}) = 1{,}93 \times 10^{-5} \quad \text{et} \quad x_{\mathrm{sol}}(n\text{-hexane}) = 7{,}00 \times 10^{-5}
2. Absence de corrélation avec la polarité : la solubilité ne varie nullement de manière monotone avec la polarité ou la permittivité diélectrique relative des solvants (\varepsilon_r(\text{hexane}) \approx 2 < \varepsilon_r(\text{acétone}) \approx 21 < \varepsilon_r(\text{méthanol}) \approx 33 < \varepsilon_r(\text{DMSO}) \approx 47 < \varepsilon_r(\text{eau}) \approx 78). L'eau, solvant le plus polaire, possède la solubilité la plus faible.
3. Solubilités élevées dans des solvants de natures diverses : l'acétone (polaire aprotique), le méthanol (polaire protique) et le DMSO (très polaire aprotique) présentent tous de très fortes solubilités (x_{\mathrm{sol}} \in [0{,}23 \,;\, 0{,}32]).

Une approche purement qualitative centrée sur les interactions soluté-solvant est donc insuffisante. La dissolution résulte d'une compétition thermodynamique complexe entre :

• la rupture des interactions au sein du cristal solide (cohésion par liaisons hydrogène intermoléculaires et empilement aromatique \pi-\pi) ;
• la rupture des interactions solvant-solvant (très coûteuse dans l'eau en raison du réseau dense de liaisons hydrogène, conduisant à l'effet hydrophobe) ;
• les effets entropiques associés à la mise en solution et à la solvatation.

L'équilibre de dissolution du nonivamide s'écrit :

À l'équilibre physico-chimique sous la température T et la pression standard P^\circ, l'égalité des potentiels chimiques entre le corps pur solide et le soluté en solution s'écrit :

• Le nonivamide solide est sous forme de corps pur dans sa propre phase, son activité vaut a_{\text{s}} = 1 :

  \mu(\text{nonivamide(s)}) = \mu_{\text{s}}^\circ(T)

• En solution, en adoptant une convention d'activité égale à la fraction molaire (a_{\text{solv}} = x_{\text{sol}}) pour la solution idéale :

  \mu(\text{nonivamide(solv)}) = \mu_{\text{solv}}^\circ(T) + R T \ln(x_{\text{sol}})

À l'équilibre, on a donc :

On en déduit l'expression de la solubilité x_{\text{sol}} :