Mécanique des fluides

Exercice 54 ⭐️⭐️ Un écoulement visqueux, PC/L2

Un fluide newtonien de masse volumique ρ\displaystyle \rho constante, de viscosité dynamique η\displaystyle \eta se trouve entre deux plaques parallèles considérées comme d’extension infinie.

Une plaque, en y=0\displaystyle y=0 est fixe, et l’autre, en y=h\displaystyle y=h, est mobile de vitesse constante Vp=Vpux\displaystyle \overrightarrow{\mathrm{}V_{p}}=V_{p}\overrightarrow{u_{x}}. On se place en régime stationnaire. L’écoulement est unidirectionnel. On néglige la pesanteur.

  1. Montrer que la vitesse du fluide est de la forme V=V(y)ux\displaystyle \mathrm{}\overrightarrow{V}=V(y)\overrightarrow{u_{x}}.

  2. Montrer que V(y)=Vp(yhh22ηVpdPdxyh(1yh))\displaystyle \mathbb{}V(y)=V_{p}(\frac{y}{h}-\frac{h^{2}}{2\eta V_{p}}\frac{dP}{dx}\frac{y}{h}(1-\frac{y}{h})) et exprimer la loi de variation du champ de pression.

  3. Soit P\displaystyle \mathbb{}P' tel que dPdx\displaystyle \frac{dP'}{dx}=h22ηVpdPdx.\displaystyle \frac{h^{2}}{2\eta V_{p}}\frac{dP}{dx}. Montrer que dPdx\displaystyle \frac{dP'}{dx} est sans dimension. Que peut-on dire du profil des vitesses selon le signe de dPdx\displaystyle \frac{dP'}{dx} ?

  4. Déterminer le débit volumique Dv\displaystyle \mathrm{}D_{v} à travers une section d’abscisse x\displaystyle \mathrm{}x, de largeur \displaystyle \mathrm{}\ell selon l’axe(z)\displaystyle (z). Trouver la valeur de dPdx\displaystyle \mathrm{}\frac{dP}{dx} tel que Dv=0\displaystyle \mathrm{}D_{v}=0. Donner dans ce cas l’expression de V(y)\displaystyle \mathrm{}V(y) et représenter l’allure du champ des vitesses.

On rappelle l’équation de Navier-Stockes
ρ[Vt+(V.grad)V]=ηΔ(V)grad(P)+ρg\displaystyle \rho\left[\frac{\partial\overrightarrow{V}}{\partial t}+(\overrightarrow{V}.\overrightarrow{\textrm{grad}})\overrightarrow{V}\right]=\eta\Delta(\overrightarrow{V})-\overrightarrow{\textrm{grad}}(P)+\overrightarrow{\rho g}.

  1. Le fluide est mis en mouvement par la plaque dans la direction (Ox)\displaystyle \left(Ox\right), V=Vux\displaystyle \mathrm{}\overrightarrow{V}=V\overrightarrow{u_{x}}. L’écoulement est stationnaire, le champ des vitesses ne dépend pas de t\displaystyle t. Il y a invariance du problème par translation selon (Oz)\displaystyle \left(Oz\right) et (Ox)\displaystyle \left(Ox\right), donc V\displaystyle V ne dépend que de y\displaystyle y. On vérifie que divV=0\displaystyle \textrm{div}\overrightarrow{V}=0 pour cet écoulement incompressible.

  2. Dans l’équation de Navier-Stockes, les deux termes de gauche sont nuls car l’écoulement est stationnaire d’une part, et
    (V.grad)V=V(y)xV(y)ux=0\displaystyle (\overrightarrow{V}.\overrightarrow{\textrm{grad}})\overrightarrow{V}=V(y)\frac{\partial}{\partial x}V(y)\overrightarrow{u_{x}}=\overrightarrow{0}.
    Le poids est négligé, d’où ηΔ(V)grad(P)=0\displaystyle \eta\Delta(\overrightarrow{V})-\overrightarrow{\textrm{grad}}(P)=\overrightarrow{0}, soit
    {Py=Pz=0Px=dPdx=ηd2Vdy2.\displaystyle \begin{cases} \frac{\partial P}{\partial y}=\frac{\partial P}{\partial z}=0\\ \frac{\partial P}{\partial x}=\frac{dP}{dx}=\eta\frac{d^{2}V}{dy^{2}}. \end{cases}
    La dernière équation montre que dPdx=ηd2Vdy2\displaystyle \frac{dP}{dx}=\eta\frac{d^{2}V}{dy^{2}} est une constante puisque la pression ne dépend que de x\displaystyle x, et la vitesse ne dépend que y\displaystyle y.
    V(y)=12ηdPdxy2+ay+b\displaystyle \Rightarrow V(y)=\frac{1}{2\eta}\frac{dP}{dx}y^{2}+ay+b.

Les conditions aux limites V(y=0)=0\displaystyle V(y=0)=0 et V(y=h)=VP\displaystyle V(y=h)=V_{P} imposent b=0\displaystyle b=0 et
12ηdPdxh2+ah=VP\displaystyle \frac{1}{2\eta}\frac{dP}{dx}h^{2}+ah=V_{P}, soit a=VPh12ηdPdxh\displaystyle a=\frac{V_{P}}{h}-\frac{1}{2\eta}\frac{dP}{dx}h, qui donne bien l’expression donnée.
dPdx=ηd2Vdy2\displaystyle \frac{dP}{dx}=\eta\frac{d^{2}V}{dy^{2}} s’intègre en P(x)=P(0)+dPdxx\displaystyle P(x)=P(0)+\frac{dP}{dx}x. La pression varie de manière affine.

  1. [h22ηVp]L2(M.L1.T1)(L.T1)=M1.L2.T2\displaystyle \left[\frac{h^{2}}{2\eta V_{p}}\right]\equiv\frac{L^{2}}{\left(M.L^{-1}.T^{-1}\right)\left(L.T^{-1}\right)}=M^{-1}.L^{2}.T^{2}.
    [dPdx]M.L1.T2L=M.L2.T2\displaystyle \left[\frac{dP}{dx}\right]\equiv\frac{M.L^{-1}.T^{-2}}{L}=M.L^{-2}.T^{-2}. Le produit des 2 est bien adimensionné.

Posons α=dPdx\displaystyle \alpha=\frac{dP'}{dx} et u=y/h\displaystyle u=y/h.

V=VP[uαu(1u)]\displaystyle V=V_{P}\left[u-\alpha u\left(1-u\right)\right].

V(u)=VP(1α+2αu)\displaystyle V'(u)=V_{P}\left(1-\alpha+2\alpha u\right).

Si un extremum existe sur ]0,1[\displaystyle \left]0,1\right[, on l’obtient pour um=α12α\displaystyle u_{m}=\frac{\alpha-1}{2\alpha}. La double condition 0<um<1\displaystyle 0<u_{m}<1 impose α>1\displaystyle \alpha>1 ou α<1\displaystyle \alpha<-1.

Dans tous les cas, la plaque tend à entraîner le fluide dans la direction ux\displaystyle \overrightarrow{u_{x}}, le gradient de pression tend à déplacer le fluide dans le sens des pressions décroissantes.
Pour α<0\displaystyle \alpha<0, le gradient de pression et la vitesse de la plaque contribuent à déplacer le fluide dans le sens ux\displaystyle \overrightarrow{u_{x}}.
Pour α>0\displaystyle \alpha>0, le gradient de pression tend à entraîner le fluide selon ux\displaystyle -\overrightarrow{u_{x}}, alors que le mouvement de la plaque entraîne le fluide selon ux\displaystyle \overrightarrow{u_{x}}.

Pour α<1,um\displaystyle \alpha<-1, u_{m} correspond à une composante maximale de la vitesse selon ux\displaystyle \overrightarrow{u_{x}}. La plaque limite la vitesse du fluide à son voisinage.

Pour 1α1\displaystyle -1\leq\alpha\leq1, il n’y a pas d’extremum sur ]0,1[\displaystyle \left]0,1\right[.

Pour α>1,um\displaystyle \alpha>1, u_{m} correspond à une composante maximale de la vitesse selon ux\displaystyle -\overrightarrow{u_{x}}.

  1. Le débit volumique
    Dv=y=0hz=0V.dydzux=y=0hV(y)dy=u=01VP[uαu(1u)]hdu=hVP[(1α)12+α13]=hVP3α6.\displaystyle \begin{aligned}D_{v}&=\int_{y=0}^{h}\int_{z=0}^{\ell}\overrightarrow{V}.dydz\overrightarrow{u_{x}}=\ell\int_{y=0}^{h}V(y)dy\\&=\ell\int_{u=0}^{1}V_{P}\left[u-\alpha u\left(1-u\right)\right]hdu\\&=h\ell V_{P}\left[\left(1-\alpha\right)\frac{1}{2}+\alpha\frac{1}{3}\right]\\&=h\ell V_{P}\frac{3-\alpha}{6}.\end{aligned}

Le débit est nul pour α=3\displaystyle \alpha=3 (on est bien dans le cas α>1\displaystyle \alpha>1 où les effets du gradient de pression et de la vitesse de la plaque sont contraires), soit dPdx=6ηVPh2.\frac{dP}{dx}=\frac{6\eta V_{P}}{h^{2}}.

Dans ce cas V(y)=VPyh(3yh2)\displaystyle V(y)=V_{P}\frac{y}{h}\left(3\frac{y}{h}-2\right). La vitesse s’annule pour y=2h3\displaystyle y=\frac{2h}{3}, elle est extrêmale pour y=h3,Vm=VP3\displaystyle y=\frac{h}{3},V_{m}=-\frac{V_{P}}{3}.

Exercice 64 ⭐️⭐️⭐️ Vase de Tantale, PC/PSI/L2

Hypothèses: l’eau est un fluide parfait incompressible de masse volumique constante ρ\displaystyle \rho, et Ss\displaystyle S\gg s, H\displaystyle H et h0\displaystyle h_{0} sont inférieurs à 10 cm.

  1. La hauteur d’eau dans le récipient cylindrique est repérée par h\displaystyle h. On suppose que le siphon est amorcé. Déterminer le débit volumique D\displaystyle D dans le siphon en faisant des hypothèses adaptées.

  2. Donner l’équation différentielle vérifiée par h(t)\displaystyle h(t) et la résoudre pour h=h0\displaystyle h=h_{0} à t=0\displaystyle t=0. Quand la vidange s’arrête-t-elle?

  3. Un robinet est placé au-dessus de la cuve avec un débit d\displaystyle d. Initialement le récipient est vide. Discuter les différents régimes possibles. Dans le cas d’un régime périodique, déterminer la période.

Ecoulement parfait, homogène, incompressible 👉 Il manque l’hypothèse stationnaire pour utiliser la relation de Bernoulli sur une ligne de courant. Une approximation quasi stationnaire est-elle possible?

Comme la masse volumique est constante, l’écoulement est parfait, homogène, incompressible.
L’hypothèse d’un régime quasi-stationnaire dans le récipient est justifiée par le fait que Ss\displaystyle S\gg s, ce qui impose h˙=sSVsVs\displaystyle \left|\dot{h}\right|=\frac{s}{S}V_{s}\ll V_{s}, vitesse de sortie de l’eau (voir Remarque 1 en fin de corrigé).

Le théorème de Bernoulli s’applique sur une ligne de champ depuis un point A de la surface jusqu’à un point B situé au niveau de la sortie :

12ρvA2+p0+ρgh=12ρvB2+p0ρgH.\frac{1}{2}\rho v_{A}^{2}+p_{0}+\rho gh=\frac{1}{2}\rho v_{B}^{2}+p_{0}-\rho gH.

Or vA=h˙=sSvBvB\displaystyle v_{A}=-\dot{h}=\frac{s}{S}v_{B}\ll v_{B}, d’où v=2g(h+H).\displaystyle v=\sqrt{2g(h+H)}.

D=s2g(h+H).D=s\sqrt{2g(h+H)}.

  1. h˙=sS2g(h+H)\displaystyle \dot{h}=-\frac{s}{S}\sqrt{2g(h+H)}.

En séparant les variables dhh+H=sS2gdt\displaystyle \frac{dh}{\sqrt{h+H}}=-\frac{s}{S}\sqrt{2g}dt, qui s’intègre en

h+Hh0+H=sSg2t.\sqrt{h+H}-\sqrt{h_{0}+H}=-\frac{s}{S}\sqrt{\frac{g}{2}}t.

La vidange s’arrête quand le siphon se désamorce, pour h=0\displaystyle h=0.

  1. Au départ, le récipient se remplit jusqu’à ce que l’eau atteigne h=h0\displaystyle h=h_{0}. Alors le siphon s’amorce. Le récipient se remplit avec le débit d\displaystyle d et se vide avec le débit D(h)=s2g(h+H)\displaystyle D(h)=s\sqrt{2g(h+H)}.
  • Si d>D(h1)\displaystyle d>D(h_{1}), le niveau continue à monter et le récipient déborde.

  • Si D(h0)<dD(h1)\displaystyle D(h_{0})<d\leq D(h_{1}), le niveau de l’eau continue à monter et se stabilise à h>h_{0} tel que d=D(h).

  • Si D(0+)=s2gH<d<D(h0)\displaystyle D(0^{+})=s\sqrt{2gH}<d<D(h_{0}), le niveau, après être monté à h0\displaystyle h_{0} pour l’amorçage du siphon, va diminuer jusqu’à se stabiliser à 0<hh1\displaystyle 0<h\leq h_{1} tel que d=D(h)\displaystyle d=D(h).

  • Si d<s2gH\displaystyle d<s\sqrt{2gH}, le niveau , après être monté à h0\displaystyle h_{0} pour l’amorçage du siphon, va diminuer jusqu’à h=0\displaystyle h=0, le siphon va alors se désamorcer. Le débit de sortie s’annule alors que le recipient va recommencer à se remplir. Le niveau va à nouveau monter jusqu’à h0\displaystyle h_{0}, avant de diminuer. Un régime périodique s’établit et le niveau va osciller (oscillations de relaxation).
    Dans ce dernier cas, on parle du vase de Tantale par référence au héros mythologique Tantale, placé au milieu d’un fleuve qui s’assèche à chaque fois que Tantale se penche pour boire.
    Calculons la période des oscillations.
    La durée T1\displaystyle T_{1} du remplissage, de h=0\displaystyle h=0 à h0\displaystyle h_{0} est T1=Sh0d.T_{1}=\frac{Sh_{0}}{d}.
    La durée T2\displaystyle T_{2} de la vidange s’obtient avec Sh˙=s2g(h+H)+d\displaystyle S\dot{h}=-s\sqrt{2g(h+H)}+d. Posons u=h+H\displaystyle u=\sqrt{h+H}, d’où 2udu=dh\displaystyle 2udu=dh.
    En séparant les variables
    S2udu=dt(s2gu+d)=s2g(u0u)dt\displaystyle S2udu=dt\left(-s\sqrt{2g}u+d\right)=s\sqrt{2g}\left(u_{0}-u\right)dtu0=ds2g\displaystyle u_{0}=\frac{d}{s\sqrt{2g}},
    uu0udu=(1+u0u0u)du=s2S2gdt\displaystyle \frac{u}{u_{0}-u}du=(-1+\frac{u_{0}}{u_{0}-u})du=\frac{s}{2S}\sqrt{2g}dt, qui s’intégre en [uu0lnuu0]H+h0H=s2S2gT2\displaystyle \left[-u-u_{0}\ln\left|u-u_{0}\right|\right]_{\sqrt{H+h_{0}}}^{\sqrt{H}}=\frac{s}{2S}\sqrt{2g}T_{2}, soit
    T2=Ss2g[H+h0H+ds2gln(H+h0ds2gHds2g)],T_{2}=\frac{S}{s}\sqrt{\frac{2}{g}}\left[\sqrt{H+h_{0}}-\sqrt{H}+\frac{d}{s\sqrt{2g}}\ln\left(\frac{\sqrt{H+h_{0}}-\frac{d}{s\sqrt{2g}}}{\sqrt{H}-\frac{d}{s\sqrt{2g}}}\right)\right], homogène à défaut d’être joli 😏.
    La période T=T1+T2T=T_1+T_2

Remarque 1 — L’écoulement n’est pas stationnaire a priori. Cependant, soit V un ordre de grandeur du champ des vitesses dans le récipient à proximité de la conduite, et τ\displaystyle \tau une durée caractéristique de variation temporelle du champ des vitesses (et de h\displaystyle h), vtVτ\displaystyle \left\Vert \frac{\partial\overrightarrow{v}}{\partial t}\right\Vert \sim\frac{V}{\tau} et (v.grad)vV2h0\displaystyle \left\Vert \left(\overrightarrow{v.}\overrightarrow{\textrm{grad}}\right)\overrightarrow{v}\right\Vert \sim\frac{V^{2}}{h_{0}}.

Par conservation du débit volumique h˙=sSVV,\displaystyle \left|\dot{h}\right|=\frac{s}{S}V\ll V, soit h0τsSV.\displaystyle \frac{h_{0}}{\tau}\sim\frac{s}{S}V.

vt/(v.grad)vh0τVsS1\displaystyle \left\Vert \frac{\partial\overrightarrow{v}}{\partial t}\right\Vert /\left\Vert \left(\overrightarrow{v.}\overrightarrow{\textrm{grad}}\right)\overrightarrow{v}\right\Vert \sim\frac{h_{0}}{\tau V}\sim\frac{s}{S}\ll1, ce qui justifie l’approximation quasi stationnaire.

Remarque 2 — La pression minimale dans le tube, où la vitesse est constante, s’obtient en h=h0\displaystyle h=h_{0}. Etant données les valeurs de H\displaystyle H et h0\displaystyle h_{0}, P(h0)=P0ρg(h0+H)\displaystyle P(h_{0})=P_{0}-\rho g(h_{0}+H) reste suffisament élevée pour éviter cavitation et désamorçage à cause du tube.

Exercice 68 ⭐️⭐️ Oscillations dans un tube en U, Spé/L2

Un tube en U (section S\displaystyle S identique sur les différentes parties du tube), dont les deux bras sont verticaux, est rempli d’un volume SL\displaystyle SL d’eau. L’écoulement est parfait, le fluide incompressible de masse volumique ρ\displaystyle \rho, et on considère que la vitesse a même norme en tout point du fluide.

Au repos, la surface libre du liquide est à l’altitude commune z=0\displaystyle z=0. À la suite d’une perturbation initiale, il apparait une dénivellation entre les deux branches du tube h(t)=zA(t)zB(t)\displaystyle h(t)=z_{A}(t)-z_{B}(t). En faisant un bilan d’énergie mécanique, calculer la période T\displaystyle T des oscillations du liquide dans le tube.

La vitesse a même norme en tout point du fluide : v=z˙A=z˙B\displaystyle v=\left|\dot{z}_{A}\right|=\left|\dot{z}_{B}\right|.
L’énergie cinétique du fluide Ec=12ρSLz˙A2.E_{c}=\frac{1}{2}\rho SL\dot{z}_{A}^{2}.

Exprimons l’énergie potentielle du fluide en prenant Ep=0\displaystyle E_{p}=0 quand le fluide est à l’équilibre.
Sur le schéma, l’énergie potentielle de pesanteur supplémentaire de la partie de droite située au-dessus de z=0\displaystyle z=0 est (ρSzA)gzA2\displaystyle (\rho Sz_{A})g\frac{z_{A}}{2}, le centre de gravité de cette partie étant en zA/2\displaystyle z_{A}/2.

De même, il faut soustraire à l’énergie potentielle totale, l’énergie potentielle du fluide qui occuperait la zone situé en-dessous de z=0\displaystyle z=0, soit ρSzAgzB2\displaystyle \rho Sz_{A}g\frac{z_{B}}{2} (même volume car le fluide est incompressible).

Ep=ρgSzAzA2ρgSzAzB2=ρgSzA2.E_{p}=\rho gSz_{A}\frac{z_{A}}{2}-\rho gSz_{A}\frac{z_{B}}{2}=\rho gSz_{A}^{2}.

Le travail des forces de pression exercées sur les deux interfaces se compensent puisque les vitesses des deux interfaces sont opposées.

En négligeant tout phénomène dissipatif, l’énergie mécanique
Em=12ρSLz˙A2+ρgSzA2E_{m}=\frac{1}{2}\rho SL\dot{z}_{A}^{2}+\rho gSz_{A}^{2}
est constante.

En dérivant par rapport au temps, et pour zA˙0\displaystyle \dot{z_{A}}\neq0, Lz¨A+2gzA=0.L\ddot{z}_{A}+2gz_{A}=0.

C’est l’équation d’un oscillateur harmonique de période T=2πL2g.T=2\pi\sqrt{\frac{L}{2g}}.
En réalité, les forces de viscosité ramènent le fluide à l’équilibre.

Exercice 108 ⭐️⭐️⭐️ Amincissement d’un film vertical, PC/L2

On considère un film mince constitué par un liquide newtonien, de masse volumique constante ρ\displaystyle \rho, de viscosité cinématique ν=ηρ\displaystyle \nu=\frac{\eta}{\rho} en écoulement sur un substrat solide vertical et plan, et en contact avec l’atmosphère. Soit z=0\displaystyle z=0 la position du sommet du film, (Oz)\displaystyle \left(Oz\right) est l’axe vertical descendant, e0\displaystyle e_{0} l’épaisseur initiale uniforme du film et L\displaystyle L sa largeur invariante selon (Oy)\displaystyle (Oy). On suppose l’écoulement invariant par translation selon (Oy)\displaystyle (Oy). A cause de la gravité, ce film s’écoule et s’amincit lentement avec le temps, l’épaisseur du film est décrite par e(z,t).\displaystyle \mathbb{}e(z,t). La pression est considérée uniforme dans le film et on néglige tout phénomène lié à la tension superficielle. La viscosité dynamique de l’air est négligeable.

On rappelle l’équation de Navier-Stockes
ρ[vt+(v.grad)v]=ηΔ(v)grad(p)+ρg\displaystyle \rho\left[\frac{\partial\overrightarrow{v}}{\partial t}+(\overrightarrow{v}.\overrightarrow{\textrm{grad}})\overrightarrow{v}\right]=\eta\Delta(\overrightarrow{v})-\overrightarrow{\textrm{grad}}(p)+\overrightarrow{\rho g}.

On cherche à déterminer le champ des vitesses v(M)=v(x,z,t)ez\displaystyle \overrightarrow{v}(M)=v(x,z,t)\overrightarrow{e_{z}} correspondant à une épaisseur e(z,t)\displaystyle e(z,t) du film.

  1. Exprimer l’incompressibilité de l’écoulement. Dans le cas stationnaire, déterminer le champ des vitesses dans le fluide et le débit volumique Dv\displaystyle D_{v} en fonction de l’épaisseur du film e(z), \displaystyle \mathbb{}e(z),\textrm{ }ν,\displaystyle \nu, g.

  2. Dans l’approximation quasi-stationnaire, on considère que l’expression de Dv\displaystyle D_{v} obtenue au 1) reste valide bien que l’épaisseur soit maintenant une fonction de z\displaystyle z et du temps. En exprimant la conservation de la matière, déterminer, dans le cadre de cette approximation, l’équation vérifiée par l’épaisseur e(z,t).\displaystyle \mathrm{}e(z,t).

  3. Chercher la solution sous la forme e(z,t)=f(z)g(t)\displaystyle \mathrm{}e(z,t)=f(z)g(t).
    Le film ne peut que s’amincir. En déduire sur quel intervalle [0,h(t)]\displaystyle \left[0,h(t)\right] en z\displaystyle z, l’expression l’expression obtenue pour e(z,t)\displaystyle e(z,t) peut être valide. Représenter l’allure du film à un instant t\displaystyle t donné.

Pour t>0\displaystyle t>0, la vitesse étant non nulle pour z0\displaystyle z\neq0, e(z=0,t)0\displaystyle e(z=0,t)\rightarrow0.
Pour t0\displaystyle t\rightarrow0, e(z=0,t=0)=e0\displaystyle e(z=0,t=0)=e_{0}.

  1. L’écoulement est incompressible, divv=0vz=0\displaystyle \textrm{div}\overrightarrow{v}=0\Longrightarrow\frac{\partial v}{\partial z}=0.
    En régime stationnaire, vt=0\displaystyle \frac{\partial\overrightarrow{v}}{\partial t}=\overrightarrow{0}, et (v.grad)v=vvzez=0\displaystyle (\overrightarrow{v}.\overrightarrow{\textrm{grad}})\overrightarrow{v}=v\frac{\partial v}{\partial z}\overrightarrow{e_{z}}=\overrightarrow{0}. D’où, avec grad(p)=0\displaystyle \overrightarrow{\textrm{grad}}(p)=\overrightarrow{0}, ηv"(x)ez=ρg=ρgez\displaystyle \eta v"(x)\overrightarrow{e_{z}}=-\rho\overrightarrow{g}=-\rho g\overrightarrow{e_{z}}. Il s’en déduit
    v(x)=ρgηx2+Ax+B,v(x)=-\frac{\rho g}{\eta}x^{2}+Ax+B,
    A,B\displaystyle A,B sont des constantes.
    -En x=0\displaystyle x=0 au contact du solide immobile, v(0)=0B=0\displaystyle v(0)=0\Rightarrow B=0.
    -En x=e\displaystyle x=e , la force de viscosité au contact de l’air est négligeable
    v(e)=00=2ρgηe+A\displaystyle v'(e)=0\Rightarrow0=-\frac{2\rho g}{\eta}e+A.
    D’où v(x)=ρgηx2+2ρgηex\displaystyle v(x)=-\frac{\rho g}{\eta}x^{2}+\frac{2\rho g}{\eta}ex,
    v(x)=ρg2ηx(2ex).v(x)=\frac{\rho g}{2\eta}x\left(2e-x\right).
    Le débit volumique est Dv=L0ev(x)dx=Lρg2η[e313e3]\displaystyle D_{v}=L\int_{0}^{e}v(x)dx=L\frac{\rho g}{2\eta}\left[e^{3}-\frac{1}{3}e^{3}\right], soit Dv=gLνe33.D_{v}=\frac{gL}{\nu}\frac{e^{3}}{3}.

  2. Effectuons un bilan de matière entre t\displaystyle t et t+dt\displaystyle t+dt pour le système
    ouvert entre z\displaystyle z et z+dz\displaystyle z+dz :
    -masse entrant en z entre t\displaystyle t et t+dt\displaystyle t+dt : ρDv(z,t).dt\displaystyle \rho D_{v}(z,t).dt
    -masse sortant en z entre t\displaystyle t et t+dt\displaystyle t+dt : ρDv(z+dz,t).dt\displaystyle \rho D_{v}(z+dz,t).dt
    -variation de masse du système entre t\displaystyle t et t+dt\displaystyle t+dt : ρLdz.etdt\displaystyle \rho Ldz.\frac{\partial e}{\partial t}dt.
    D’où
    ρLdz.etdt=ρDv(z,t).dtρDv(z+dz,t).dt=ρDvz.dzdt\displaystyle \begin{aligned}\rho Ldz.\frac{\partial e}{\partial t}dt&=\rho D_{v}(z,t).dt-\rho D_{v}(z+dz,t).dt\\&=-\rho\frac{\partial D_{v}}{\partial z}.dzdt\end{aligned},

Avec Dvz=Let\displaystyle \frac{\partial D_{v}}{\partial z}=L\frac{\partial e}{\partial t}\Rightarrow
et=gνe2ez.\frac{\partial e}{\partial t}=-\frac{g}{\nu}e^{2}\frac{\partial e}{\partial z}.

  1. En remplaçant dans l’équation du 2), f(z)g(t)=gν[f(z)]2[g(t)]3f(z)\displaystyle f(z)g'(t)=-\frac{g}{\nu}\left[f(z)\right]^{2}\left[g(t)\right]^{3}f'(z).
    Supposons qu’il existe une valeur z0\displaystyle z_{0} de z\displaystyle z telle que f(z0)f(z0)0\displaystyle f'(z_{0})f(z_{0})\neq0,
    alors g(t)=gν[g(t)]3f(z0)f(z0)\displaystyle g'(t)=-\frac{g}{\nu}\left[g(t)\right]^{3}f'(z_{0})f(z_{0}).
    g(t)[g(t)]3=gνf(z0)f(z0)=K1 constante\displaystyle \Rightarrow\frac{g'(t)}{\left[g(t)\right]^{3}}=-\frac{g}{\nu}f'(z_{0})f(z_{0})=K_{1}\textrm{ constante},
    soit ddt(1g2(t))=2K1\displaystyle \frac{d}{dt}\left(\frac{1}{g^{2}(t)}\right)=-2K_{1} ou 1g2(t)=2K1t+K2\displaystyle \frac{1}{g^{2}(t)}=-2K_{1}t+K_{2},
    soit g(t)=±1K22K1t\displaystyle g(t)=\pm\frac{1}{\sqrt{K_{2}-2K_{1}t}}.
    Alors K1f(z)=gν[f(z)]2f(z)\displaystyle K_{1}f(z)=-\frac{g}{\nu}\left[f(z)\right]^{2}f'(z). Pour f(z)0\displaystyle f(z)\neq0, f(z)f(z)=K1νg\displaystyle f'(z)f(z)=-\frac{K_{1}\nu}{g}, soit ddx[f2(z)]=2K1νg\displaystyle \frac{d}{dx}\left[f^{2}(z)\right]=-2\frac{K_{1}\nu}{g} et f(z)=±2K1νgz+K3\displaystyle f(z)=\pm\sqrt{-2\frac{K_{1}\nu}{g}z+K_{3}}.
    Finalement, e(z,t)=2K1νgz+K3K22K1t\displaystyle e(z,t)=\sqrt{\frac{-2\frac{K_{1}\nu}{g}z+K_{3}}{K_{2}-2K_{1}t}}.
  • Pour t>0\displaystyle t>0, e(0,t)0\displaystyle e(0,t)\rightarrow0, d’où K3=0\displaystyle K_{3}=0 et e(z,t)=2νgK1zK22K1t\displaystyle e(z,t)=\sqrt{2\frac{\nu}{g}}\sqrt{\frac{-K_{1}z}{K_{2}-2K_{1}t}}

  • Pour t0\displaystyle t\rightarrow0, e(0,0)=e0K2=0\displaystyle e(0,0)=e_{0}\Rightarrow K_{2}=0 et e(z,t)=νgzt\displaystyle e(z,t)=\sqrt{\frac{\nu}{g}}\sqrt{\frac{z}{t}}.

Le film ne peut que s’amincir, donc e(z,t)e0zgνe02t\displaystyle e(z,t)\leq e_{0}\Rightarrow z\leq\frac{g}{\nu}e_{0}^{2}t.
L’expression (homogène avec ν\displaystyle \nu en m2\displaystyle ^2.s1\displaystyle ^{-1}), e(z,t)=νgzte(z,t)=\sqrt{\frac{\nu}{g}}\sqrt{\frac{z}{t}} n’est valide que pour zh(t)=gνe02t.z\leq h(t)=\frac{g}{\nu}e_{0}^{2}t.

Pour z>h(t),\displaystyle z>h(t), la solution e(z,t)=e0\displaystyle e(z,t)=e_{0} convient et vérifie l’équation
du 2) avec g(t)=0\displaystyle g'(t)=0 et f(z)=0\displaystyle f'(z)=0.

Exercice 149 ⭐️⭐️⭐️ Vents dominants, PC/L2

D’après Centrale Supelec
La Terre est assimilée à une sphère de centre O\displaystyle O et le référentiel géocentrique est considéré comme galiléen. Dans l’exercice on se place dans le référentiel terrestre RT\displaystyle \mathcal{R}_{T}, en rotation uniforme autour de l’axe OZ avec Ω=ΩuZ\displaystyle \overrightarrow{\Omega}=\Omega\overrightarrow{u}_{Z}.
À un point M de la surface terrestre, on associe un repère local (M,x,y,z)\displaystyle \left(M,x,y,z\right) avec (Mx)\displaystyle \left(Mx\right) orienté vers le sud, (My)\displaystyle \left(My\right) vers l’Est, et on considère que l’axe (Mz)\displaystyle \left(Mz\right) peut être confondu avec (Oz)\displaystyle \left(Oz\right).

On considère que l’écoulement de l’air est parfait (viscosité négligée) et on s’intéresse aux vents dominants au niveau de la mer, c’est-à-dire des écoulements à grande échelle et réguliers dans le temps.

  1. Le champ eulérien des vitesses est noté v(x,y,z).\displaystyle \overrightarrow{v}(x,y,z). Montrer qu’il obéit à l’équation ρa=ρggradP2ρΩv\rho\overrightarrow{a}=\rho\overrightarrow{g}-\overrightarrow{\textrm{grad}}P-2\rho\overrightarrow{\Omega}\land\overrightarrow{v}
    où on rappelle l’accélération a=DvDt=vt+(v.grad)v\displaystyle \overrightarrow{a}=\frac{D\overrightarrow{v}}{Dt}=\frac{\partial\overrightarrow{v}}{\partial t}+\left(\overrightarrow{v}.\overrightarrow{\textrm{grad}}\right)\overrightarrow{v}, P\displaystyle P la pression et g\displaystyle \overrightarrow{g} est le champ de pesanteur. Pourquoi la force d’inertie d’entraînement n’apparait pas ?

  2. Si l’air est en équilibre en tout point dans RT\displaystyle \mathcal{R}_{T}, quelle est l’équation vérifiée par la pression Peˊq(x,y,z)\displaystyle P_{\acute{e}q}(x,y,z)?
    Soit p(x,y,z)=P(x,y,z)Peˊq(x,y,z)\displaystyle p(x,y,z)=P(x,y,z)-P_{\acute{e}q}(x,y,z). Quelle est l’équation (1) reliant ρ,v,p\displaystyle \rho,\overrightarrow{v},p et Ω\displaystyle \overrightarrow{\Omega} ?

  3. On se limite dans la suite à l’étude d’écoulement stationnaires dans RT\displaystyle \mathcal{R}_{T}.
    Soit U\displaystyle U une grandeur caractéristique des vitesses d’un écoulement, L\displaystyle L une distance horizontale caractéristique de celui-ci. Le nombre de Rossby R0\displaystyle R_{0} est le rapport sans dimension entre le terme d’accélération et le terme dû à la force de Coriolis. Evaluer R0\displaystyle R_{0} pour un écoulement atmosphérique typique avec U=10 m.s1\displaystyle U=10\textrm{ m.s}^{-1}, L=1000 km\displaystyle L=1000\textrm{ km}. Commenter.

  4. On considère des écoulements géostrophiques, c’est-à-dire des écoulements atmosphériques pour lesquels R01\displaystyle R_{0}\ll1, et on ne garde que les termes dominants dans (1). La faible épaisseur de l’atmosphère (devant L\displaystyle L) permet de négliger la composante verticale vz\displaystyle v_{z}. Etablir l’expression de v\displaystyle \overrightarrow{v} en fonction de ez,grad p,ρ,Ω\displaystyle \overrightarrow{e_{z}},\overrightarrow{\textrm{grad }}p,\rho,\Omega et de la latitude λ\displaystyle \lambda. Montrer que la composante horizontale de grad p\displaystyle \overrightarrow{\textrm{grad }}p n’intervient pas.
    La valeur moyenne de la pression atmosphérique au sol a été étudiée sur l’ensemble du globe terrestre : voir ici. Il s’en déduit que la pression au sol ne dépend essentiellement que de la latitude λ\displaystyle \lambda, selon l’allure ci-dessous (on considère que la variation est schématiquement une fonction paire de λ\displaystyle \lambda).

    En déduire la direction des vents dominants (les alizés) entre les tropiques (entre λ=23\displaystyle \lambda=-23{^\circ}et λ=23\displaystyle \lambda=23{^\circ}). De même quelle est la direction des vents dominants en Europe de l’ouest ?

Référentiel terrestre non galiléen 👉 tenir compte des forces d’inertie. Le champ de pesanteur est la somme du champ gravitationnel et de la force d’inertie d’entraînement par unité de masse.
Evaluation des ordres de grandeur 👉 Pour la dérivation par rapport à l’espace, diviser par la longueur caractistique de variation de la fonction ; par exemple l’accélération convective (v.grad)vU2L\displaystyle \left\Vert \left(\overrightarrow{v}.\overrightarrow{\textrm{grad}}\right)\overrightarrow{v}\right\Vert \sim\frac{U^{2}}{L}

  1. Il s’agit de l’équation d’Euler dans le référentiel RT\displaystyle \mathcal{R}_{T} non galiléen. Il faut tenir compte des forces d’inertie de Coriolis et d’entraînement. La force de Coriolis volumique correspond au terme 2ρΩv\displaystyle -2\rho\overrightarrow{\varOmega}\land\overrightarrow{v}. La force volumique d’inertie d’entraînement est comprise dans le champ de pesanteur, somme du champ gravitationnel et du terme d’inertie d’entraînement.

  2. Pour v=0\displaystyle \overrightarrow{v}=\overrightarrow{0}, ρggradPeˊq=0.\rho\overrightarrow{g}-\overrightarrow{\textrm{grad}}P_{\acute{e}q}=\overrightarrow{0}.
    La pression ne varie que verticalement. Par différence
    ρ[vt+(v.grad)v]=grad p2ρΩv.\rho\left[\frac{\partial\overrightarrow{v}}{\partial t}+\left(\overrightarrow{v}.\overrightarrow{\textrm{grad}}\right)\overrightarrow{v}\right]=-\overrightarrow{\textrm{grad }}p-2\rho\overrightarrow{\Omega}\land\overrightarrow{v}.

  3. L’écoulement est stationnaire vt=0\displaystyle \frac{\partial\overrightarrow{v}}{\partial t}=\overrightarrow{0}.
    R0(v.grad)v2ΩvU2L2ΩUU2ΩL\displaystyle R_{0}\sim\frac{\left\Vert \left(\overrightarrow{v}.\overrightarrow{\textrm{grad}}\right)\overrightarrow{v}\right\Vert }{\left\Vert 2\overrightarrow{\varOmega}\land\overrightarrow{v}\right\Vert }\sim\frac{\frac{U^{2}}{L}}{2\Omega U}\sim\frac{U}{2\Omega L}.
    La période de rotation de la Terre est 24 heures, soit Ω=2π24.3600=7.105rad.s1\displaystyle \Omega=\frac{2\pi}{24.3600}=7.10^{-5}\textrm{rad.s}^{-1}, et R01014.10=7.102\displaystyle R_{0}\simeq\frac{10}{14.10}=7.10^{-2}. Le terme de Coriolis domine pour cet écoulement.

  4. Pour R01\displaystyle R_{0}\ll1, (1)\displaystyle \left(1\right) devient grad p=2ρΩv\displaystyle \overrightarrow{\textrm{grad }}p=-2\rho\overrightarrow{\Omega}\land\overrightarrow{v}.
    En multipliant vectoriellement par ez\displaystyle \overrightarrow{e_{z}},
    ezgrad p=2ρez[Ωv]=2ρ[ΩvzΩzv]=2ρΩzv\displaystyle \begin{aligned}\overrightarrow{e_{z}}\wedge\overrightarrow{\textrm{grad }}p&=-2\rho\overrightarrow{e_{z}}\wedge\left[\overrightarrow{\Omega}\land\overrightarrow{v}\right]\\&=-2\rho\left[\overrightarrow{\Omega}v_{z}-\Omega_{z}\overrightarrow{v}\right]\\&=2\rho\Omega_{z}\overrightarrow{v}\end{aligned}.
    Avec Ωz=Ωsinλ\displaystyle \Omega_{z}=\Omega\sin\lambda, v=12ρΩsinλezgrad p.\overrightarrow{v}=\frac{1}{2\rho\Omega\sin\lambda}\overrightarrow{e_{z}}\wedge\overrightarrow{\textrm{grad }}p. Seule la composante horizontale de p\displaystyle p intervient en raison du produit vectoriel.
    Dans l’hémisphère nord, le gradient horizontal de pression entre 0 et 23° de latitude est orienté du sud vers le nord (selon ex\displaystyle -\overrightarrow{e_{x}}), et sinλ>0\displaystyle \sin\lambda>0, donc v\displaystyle \overrightarrow{v} est orienté selon ey\displaystyle -\overrightarrow{e_{y}}, les alizés soufflent de l’est vers l’ouest.
    Dans l’hémisphère sud, le gradient de pression entre 0 et -23° de latitude est orienté du nord vers le sud (selon ex\displaystyle \overrightarrow{e_{x}}), et sinλ<0\displaystyle \sin\lambda<0, donc les alizés sont aussi des vents d’est.
    En Europe (35<λ<60)\displaystyle \left(35{^\circ}<\lambda<60{^\circ}\right), le gradient horizontal de pression est dirigé du nord vers le sud, les vents dominants sont des vents d’ouest.