Physique
Agro-VétoBCPST2026
Cherchez d'abord seul. En cas de blocage, ouvrez une indication, réessayez, puis seulement la solution.
# Agro-Véto BCPST — Physique 2026 · Corrigé guidé
*Corrigé rédigé par IA (La Belle Copie), non officiel. Signalez toute erreur.*
## Présentation du sujet
**Thème.** « Les cheveux de Pélé » : ces fibres de verre volcanique naissent quand une goutte de lave projetée est étirée par le vent, puis figée par le refroidissement, avant d'être transportée parfois à des kilomètres. Épreuve de 3 heures, calculatrice autorisée, 11 pages, 33 questions. Les données numériques sont regroupées en page 3 ; trois documents (photographies, simulations, profil topographique) sont à exploiter.
**Découpage.**
| Partie | Questions | Contenu | Notions du programme |
|---|---|---|---|
| I.A | Q1–Q2 | Taille d'une goutte arrachée par le vent | analyse dimensionnelle, volume d'une sphère |
| I.B | Q3–Q5 | Quantité de mouvement d'un cylindre qui s'étire | conservation de la masse, intégrale sur un système continu |
| I.C | Q6–Q9 | Équation (1) : temps caractéristiques et longueur d'équilibre | lois d'échelle, équilibre mécanique, esprit critique |
| I.D | Q10 | Rôle de la viscosité | exploitation de photographies, tension de surface, viscosité |
| II.A–B | Q11–Q18 | Refroidissement puis solidification d'un cylindre | loi de Newton, premier principe, équation du premier ordre, changement d'état |
| II.C | Q19–Q21 | Discussion du modèle thermique | comparaison de temps, diffusion thermique, critique |
| III.A | Q22–Q25 | Tir balistique sans vent | PFD, trajectoire, portée, temps de vol |
| III.B | Q26–Q33 | Tir dans le vent | nombre de Reynolds, frottement quadratique, méthode d'Euler, régime permanent, lecture de document |
**Difficulté et longueur.** Le sujet n'a pas de difficulté calculatoire : il teste surtout la modélisation et la discussion. Les questions de cours à ne pas laisser passer sont nombreuses : Q1, Q2, Q3, Q11 à Q17, Q22 à Q25, Q26–Q27 et Q28. Les questions de commentaire (Q7, Q9, Q10, Q19, Q21, Q25, Q29, Q30) sont celles où l'on se distingue : un commentaire chiffré et comparatif vaut plus qu'une impression.
**Gestion du temps.** L'énoncé conseille de répartir le temps de manière équilibrée : une heure par partie environ. Les valeurs nécessaires à la suite sont toujours redonnées ($V_0 = 10$ mm³, $\tau_a \simeq 4$ ms, $t_{tot} \simeq 11{,}2$ s, $k \simeq 1$ m⁻¹, $Re_c \simeq 100$) : une erreur en amont ne bloque jamais l'aval.
---
## Partie I — Formation d'un cheveu
On construit les ordres de grandeur de la formation d'un filament : la taille d'une goutte (analyse dimensionnelle), la dynamique de l'étirement (quantité de mouvement d'un cylindre déformable), puis les temps caractéristiques lus sur une équation non linéaire qu'on ne résout pas. La partie se termine par une question qualitative sur la viscosité.
### I.A. Volume des gouttes
### Question 1
> [!INDICATION]
> Exprimer les dimensions de $\gamma$, $\rho_a$ et $u$ dans les unités de base (masse, longueur, temps).
> [!INDICATION]
> Chercher $r_0 = \gamma^\alpha\rho_a^\beta u^\delta$ et identifier les exposants de M, L et T.
#### Solution
$[\gamma] = \mathrm{N\cdot m^{-1}} = \mathrm{kg\cdot s^{-2}}$, $[\rho_a] = \mathrm{kg\cdot m^{-3}}$, $[u] = \mathrm{m\cdot s^{-1}}$. On cherche $r_0 = K\gamma^\alpha\rho_a^\beta u^\delta$ ($K$ sans dimension) :
$$\mathrm{kg}^{\alpha+\beta}\,\mathrm{m}^{-3\beta+\delta}\,\mathrm{s}^{-2\alpha-\delta} = \mathrm{m}
\iff
\begin{cases} \alpha+\beta = 0\\ -3\beta+\delta = 1\\ -2\alpha-\delta = 0 \end{cases}
\iff \alpha = 1,\ \beta = -1,\ \delta = -2 .$$
La solution est unique : à un facteur numérique près,
$$\boxed{r_0 = \frac{\gamma}{\rho_a u^2}}$$
Lecture physique : $\gamma/r_0$ est la surpression capillaire qui maintient la goutte, $\rho_a u^2$ la pression dynamique du vent qui la déchire ; $r_0$ est la taille pour laquelle elles s'équilibrent.
> [!METHODE]
> Pour « vérifier » une formule par analyse dimensionnelle, il suffit de calculer la dimension du second membre. Chercher les exposants est un peu plus long mais prouve en plus que c'est la seule combinaison possible.
> [!PIEGE]
> $\rho_a$ est la masse volumique de l'**air**, pas celle de la lave : c'est le vent qui fragmente.
### Question 2
> [!INDICATION]
> Calculer $r_0$, puis le volume d'une sphère de ce rayon.
> [!INDICATION]
> Convertir $r_0$ en millimètres avant de l'élever au cube limite les erreurs de puissances de 10.
#### Solution
$$r_0 = \frac{0{,}5}{1{,}3\times 20^2} = \frac{0{,}5}{520} = 9{,}6\times10^{-4}\ \mathrm{m} \approx 0{,}96\ \mathrm{mm}$$
$$V_0 = \frac43\pi r_0^3 = \frac43\pi\,(0{,}96\ \mathrm{mm})^3 \qquad \boxed{V_0 \approx 3{,}7\ \mathrm{mm^3} = 3{,}7\times10^{-9}\ \mathrm{m^3}}$$
C'est bien « légèrement trop faible » devant les $10$ mm³ retenus ensuite : facteur 2,7 sur le volume, soit seulement 1,4 sur le rayon.
### I.B. Accélération d'un cylindre qui s'étire
On passe de la goutte au cylindre. Les trois questions établissent le membre de gauche de l'équation (1).
### Question 3
> [!INDICATION]
> Deux hypothèses de l'énoncé suffisent : l'une porte sur la masse, l'autre sur la masse volumique.
#### Solution
Il n'y a pas d'échange de matière avec l'extérieur : le cylindre est un système fermé, de masse $m$ constante. La masse volumique $\rho$ de la lave est supposée constante (liquide incompressible). Le volume $V = m/\rho$ est donc constant, égal au volume initial :
$$\boxed{S(t)\,L(t) = V_0}$$
> [!PIEGE]
> La conservation de la masse seule ne suffit pas : un gaz enfermé conserve sa masse mais pas son volume. Il faut citer aussi l'incompressibilité.
### Question 4
> [!INDICATION]
> Exprimer le volume de la tranche entre $x$ et $x+\mathrm{d}x$, puis utiliser la question 3 pour éliminer $S$.
> [!INDICATION]
> À instant fixé, $L$ ne dépend pas de $x$ : on peut le sortir de l'intégrale.
#### Solution
La tranche a pour volume $S\,\mathrm{d}x$ avec $S = V_0/L$ :
$$\boxed{\mathrm{d}m = \frac{\rho V_0}{L}\,\mathrm{d}x}$$
$$p = \int_0^L \frac{x}{L}\,v(L)\,\frac{\rho V_0}{L}\,\mathrm{d}x = \frac{\rho V_0\,v(L)}{L^2}\int_0^L x\,\mathrm{d}x = \frac{\rho V_0\,v(L)}{L^2}\cdot\frac{L^2}{2}
\qquad \boxed{p = \frac{\rho V_0}{2}\,v(L)}$$
Contrôle : avec $m = \rho V_0$, $p = m\,v(L)/2 = m\,v_G$, où $v_G = v(L/2)$ est la vitesse du centre de masse, situé au milieu du cylindre.
> [!METHODE]
> Pour un système continu, on découpe en éléments de masse $\mathrm{d}m$, on exprime $\mathrm{d}m$ avec la variable d'intégration, puis on intègre. Le contrôle $p = m\,v_G$ permet de vérifier le résultat.
### Question 5
> [!INDICATION]
> Que représente $v(L)$ pour la longueur $L(t)$ du cylindre, sachant que la face gauche est fixe ?
#### Solution
La face gauche est fixe en $x = 0$ et la face droite, matérielle, se trouve en $x = L(t)$ : sa vitesse est $v(L) = \mathrm{d}L/\mathrm{d}t$. Avec $\rho$ et $V_0$ constants :
$$m\,a(t) = \frac{\mathrm{d}p}{\mathrm{d}t} = \frac{\rho V_0}{2}\,\frac{\mathrm{d}}{\mathrm{d}t}\left(\frac{\mathrm{d}L}{\mathrm{d}t}\right)
\qquad \boxed{m\,a(t) = \frac{\rho V_0}{2}\,\frac{\mathrm{d}^2L}{\mathrm{d}t^2}(t)}$$
Comme $m = \rho V_0$, cela s'écrit $a = \frac12\ddot L$ : le centre de masse, au milieu, accélère deux fois moins que l'extrémité libre.
> [!PIEGE]
> L'identification $v(L) = \mathrm{d}L/\mathrm{d}t$ n'est pas donnée par l'énoncé : elle se justifie par le fait que la face gauche est fixe.
### I.C. Équation différentielle et temps caractéristiques
L'équation (1) n'a pas de solution analytique simple. On en extrait des ordres de grandeur en comparant ses termes deux à deux.
### Question 6
> [!INDICATION]
> Ne garder dans (1) que le terme d'inertie et le terme de vent.
> [!INDICATION]
> Loi d'échelle : remplacer chaque grandeur par son ordre de grandeur, $L$ par une longueur typique et $\mathrm{d}^2L/\mathrm{d}t^2$ par une longueur divisée par le carré d'un temps.
#### Solution
Sans le terme capillaire, (1) devient
$$\frac{\rho V_0}{2}\frac{\mathrm{d}^2L}{\mathrm{d}t^2} \simeq \frac12\rho_a C_D \frac{V_0}{L}u^2 .$$
À temps courts $L \sim 2r_0$, et pour le temps caractéristique $\tau_a$ on estime $\dfrac{\mathrm{d}^2L}{\mathrm{d}t^2} \sim \dfrac{2r_0}{\tau_a^2}$. Avec $C_D \simeq 1$ :
$$\frac{\rho V_0}{2}\cdot\frac{2r_0}{\tau_a^2} \sim \frac12\rho_a\frac{V_0}{2r_0}u^2
\quad\Longrightarrow\quad \tau_a^2 \sim \frac{4r_0^2}{u^2}\,\frac{\rho}{\rho_a}
\quad\Longrightarrow\quad \boxed{\tau_a = \frac{2r_0}{u}\sqrt{\frac{\rho}{\rho_a}}}$$
Contrôle : $\tau_a = \dfrac{2\times9{,}6\times10^{-4}}{20}\times\sqrt{\dfrac{2300}{1{,}3}} \approx 9{,}6\times10^{-5}\times42 \approx 4{,}0\times10^{-3}$ s, la valeur annoncée en Q7.
Interprétation : $2r_0/u$ est le temps que met l'air pour passer devant la goutte ; le facteur $\sqrt{\rho/\rho_a} \approx 42$ traduit l'inertie de la lave, bien plus dense que l'air.
> [!METHODE]
> Raisonner par loi d'échelle, c'est écrire explicitement les substitutions ($L \to 2r_0$, $\ddot L \to 2r_0/\tau^2$) avant de calculer. On obtient ainsi une relation entre ordres de grandeur, sans résoudre l'équation.
> [!PIEGE]
> Il faut remplacer **les deux** occurrences de $L$ (dans $\ddot L$ et dans $V_0/L$) par $2r_0$ pour retrouver exactement la formule imposée.
### Question 7
> [!INDICATION]
> Attention à la conversion des mm³ en m³.
> [!INDICATION]
> Pour commenter, comparer les deux temps et dire ce que leur rapport implique pour l'étirement.
#### Solution
$V_0 = 10\ \mathrm{mm^3} = 10\times(10^{-3}\ \mathrm{m})^3 = 1{,}0\times10^{-8}\ \mathrm{m^3}$.
$$\tau_c = \sqrt{\frac{\rho V_0}{\gamma}} = \sqrt{\frac{2{,}3\times10^3\times1{,}0\times10^{-8}}{0{,}5}} = \sqrt{4{,}6\times10^{-5}}\ \mathrm{s} \qquad \boxed{\tau_c \approx 6{,}8\times10^{-3}\ \mathrm{s}}$$
Commentaire : $\tau_c/\tau_a \approx 1{,}7$. Les deux temps sont du même ordre de grandeur, quelques millisecondes : l'étirement par le vent et le rappel capillaire sont en compétition. Le vent agit un peu plus vite, ce qui permet à l'étirement de s'amorcer avant que la tension de surface ne ramène la goutte vers la sphère. Si l'on avait $\tau_a \gg \tau_c$, la goutte resterait sphérique ; si l'on avait $\tau_a \ll \tau_c$, le vent l'étirerait sans résistance.
> [!PIEGE]
> $1\ \mathrm{mm^3} = 10^{-9}\ \mathrm{m^3}$, pas $10^{-3}$ ni $10^{-6}$.
### Question 8
> [!INDICATION]
> Que vaut l'accélération à l'équilibre ?
> [!INDICATION]
> Poser $X = V_0/L_{eq}$ : l'équation devient $aX = bX^{1/2}$, qu'on simplifie par $X^{1/2} > 0$.
#### Solution
À l'équilibre, $\ddot L = 0$ :
$$\frac12\rho_a C_D \frac{V_0}{L_{eq}}u^2 = \gamma\left(\frac{V_0}{L_{eq}}\right)^{1/2}
\ \Longrightarrow\ \left(\frac{V_0}{L_{eq}}\right)^{1/2} = \frac{2\gamma}{\rho_a C_D u^2}
\ \Longrightarrow\ \boxed{L_{eq} = V_0\left(\frac{\rho_a C_D u^2}{2\gamma}\right)^2 = \frac{\rho_a^2 C_D^2 u^4 V_0}{4\gamma^2}}$$
Application numérique ($C_D = 1$) : $\dfrac{\rho_a u^2}{2\gamma} = \dfrac{1{,}3\times400}{2\times0{,}5} = 520\ \mathrm{m^{-1}}$, donc
$$L_{eq} = 1{,}0\times10^{-8}\times520^2 \qquad \boxed{L_{eq} \approx 2{,}7\times10^{-3}\ \mathrm{m} = 2{,}7\ \mathrm{mm}}$$
Avec $\rho_a u^2 = \gamma/r_0$, on obtient la forme compacte $L_{eq} = V_0/(4r_0^2)$ : la section d'équilibre est $(2r_0)^2$.
> [!PIEGE]
> $u^4 = 20^4 = 1{,}6\times10^5$. L'erreur de puissance est fréquente.
### Question 9
> [!INDICATION]
> Comparer $L_{eq}$ à la taille de la goutte de départ et à la longueur des cheveux donnée en introduction.
> [!INDICATION]
> Regarder la structure de l'équation (1) : y a-t-il un terme qui dissipe l'énergie ? quelles données de l'énoncé n'y figurent pas ?
#### Solution
- **Comparaison aux observations.** Une sphère de $10$ mm³ a un diamètre de $2{,}7$ mm : $L_{eq}$ est à peine supérieur à la taille de la goutte d'origine. Le modèle ne prévoit presque pas d'étirement, alors que les cheveux mesurent plusieurs centimètres, au moins 10 fois plus. **Le modèle n'est pas crédible quantitativement.** Il donne seulement l'échelle millimétrique et le bon sens de variation ($L_{eq}$ croît avec $u$, décroît avec $\gamma$).
- **Pas de dissipation.** L'équilibre est stable (si $L < L_{eq}$ le vent l'emporte et allonge, si $L > L_{eq}$ la capillarité l'emporte et raccourcit), mais (1) ne contient aucun terme d'amortissement : le système oscille autour de $L_{eq}$ sans s'y arrêter. $L_{eq}$ n'est pas une longueur finale.
- **Forces incomplètes.** La force appelée « force de viscosité », $F_c = -\gamma S^{1/2}$, ne contient que la tension de surface : c'est une force capillaire. La viscosité $\eta$, donnée et très élevée, n'apparaît nulle part ; la solidification, qui fige la forme, non plus.
- **Sensibilité.** $L_{eq} \propto u^4$ : $u = 40$ m·s⁻¹ donnerait $4{,}3$ cm. Une rafale change l'ordre de grandeur, ce qui rend la prédiction peu robuste.
- **Traînée.** Elle est écrite avec la section droite $S$ ; pour un cylindre en travers du vent, l'aire exposée serait plutôt $2RL$.
> [!METHODE]
> Discuter la crédibilité d'un modèle : (1) comparer chiffres à l'appui à une observation donnée dans l'énoncé ; (2) repérer ce qui manque dans l'équation (dissipation, données non utilisées) ; (3) dire dans quel sens va l'erreur.
### I.D. Viscosité
### Question 10
> [!INDICATION]
> Les deux expériences ne diffèrent que par un paramètre : lequel ? Décrire ce qui change sur les photographies.
> [!INDICATION]
> Distinguer le phénomène moteur de la rupture du filament et celui qui la freine.
> [!INDICATION]
> Transposer à la lave : sa viscosité vaut 2500 fois celle de l'eau, et elle augmente quand la lave refroidit.
#### Solution
- **Eau (figure 3).** Le pont liquide s'étrangle et s'amincit en un filament qui se rompt en quelques dizaines de millisecondes (images espacées de $8{,}3$ ms, rupture entre la 4ᵉ et la 5ᵉ). Le filament se fragmente ensuite en gouttelettes : la tension de surface tend à minimiser l'aire de l'interface, et une sphère a moins de surface qu'un long cylindre de même volume.
- **Huile siliconée (figure 4).** Même tension de surface, viscosité 100 fois plus grande : le filament s'amincit beaucoup plus lentement et reste continu, long et régulier. Le temps de rupture ($\sim 1$ s) dépasse largement la durée de l'expérience.
- **Rôle de la viscosité.** La tension de surface est le moteur de la rupture ; la viscosité s'oppose aux déformations rapides et ralentit l'écoulement du liquide hors de l'étranglement. À tension de surface égale, un fluide plus visqueux peut donc être étiré en un long filament sans se fragmenter.
- **Application à la lave.** La lave est 2500 fois plus visqueuse que l'eau : elle se comporte comme l'huile, pas comme l'eau. Sa viscosité augmente en outre très fortement quand elle refroidit, ce qui finit par figer le filament étiré par le vent. Là où la capillarité l'emporte, la matière se rassemble en gouttes qui se solidifient : ce sont les larmes de Pélé (figures 1 et 5), souvent attachées à la base des cheveux.
> [!PIEGE]
> La viscosité n'étire pas le filament : c'est le vent qui l'étire. La viscosité l'empêche de se fragmenter pendant l'étirement.
> [!JURY]
> Une question qualitative sur documents se traite en trois temps : ce que l'on observe (description précise des images, avec les échelles de temps données dans les légendes), le mécanisme physique, puis la transposition au cas étudié.
---
## Partie II — Solidification des cheveux de Pélé
Modèle thermique d'un cylindre de rayon $R = 0{,}56$ mm et de longueur $L = 1$ cm, à température uniforme : refroidissement par conducto-convection (équation du premier ordre), puis palier de solidification. On discute ensuite les hypothèses.
### II.A. Étape de refroidissement
### Question 11
> [!INDICATION]
> Loi de Newton : le flux est proportionnel à la surface et à l'écart de température. Choisir le signe pour un flux **reçu**.
#### Solution
$$\boxed{\phi = h\,S_{tot}\,(T_a - T)}$$
Contrôle du signe : pour $T > T_a$, $\phi < 0$, le cylindre perd de l'énergie, ce qui est bien le cas.
> [!PIEGE]
> Écrire $\phi = hS_{tot}(T - T_a)$ pour un flux reçu conduirait en Q13 à une température qui croît indéfiniment.
### Question 12
> [!INDICATION]
> Exprimer $S_{lat}$ et $S$ en fonction de $R$ et $L$, puis calculer le rapport avec les valeurs de l'énoncé.
#### Solution
$$\frac{S_{lat}}{S} = \frac{2\pi RL}{\pi R^2} = \frac{2L}{R} = \frac{2\times1{,}0\times10^{-2}}{5{,}6\times10^{-4}} \approx 36$$
$S_{tot} = S_{lat} + 2S = S_{lat}(1 + R/L)$ : les deux bases ne comptent que pour $5{,}6$ %. On les néglige :
$$\boxed{S_{tot} \simeq 2\pi RL \qquad \phi \simeq 2\pi RL\,h\,(T_a - T)}$$
### Question 13
> [!INDICATION]
> Définir le système, puis appliquer le premier principe entre $t$ et $t+\mathrm{d}t$ à une phase condensée.
> [!INDICATION]
> Exprimer la masse en fonction de $\rho$, $R$ et $L$, puis simplifier.
#### Solution
Système : le cylindre de lave, fermé, de masse $m = \rho\pi R^2 L$, phase condensée. Entre $t$ et $t+\mathrm{d}t$, sa variation d'énergie interne est $m\,c\,\mathrm{d}T$ et il reçoit le transfert thermique $\phi\,\mathrm{d}t$ :
$$\rho\pi R^2L\,c\,\frac{\mathrm{d}T}{\mathrm{d}t} = 2\pi RL\,h\,(T_a - T).$$
En divisant par $2\pi RLh$ :
$$\boxed{\frac{\rho R c}{2h}\,\frac{\mathrm{d}T}{\mathrm{d}t}(t) + T(t) = T_a}$$
$L$ disparaît : seul compte le rapport volume/surface, $R/2$. Un cheveu plus fin refroidit plus vite.
> [!METHODE]
> Bilan thermique : système, système fermé, $\mathrm{d}U = mc\,\mathrm{d}T$ pour une phase condensée, transfert $\phi\,\mathrm{d}t$, puis mise sous forme canonique.
### Question 14
> [!INDICATION]
> Comparer à la forme canonique $\tau\,\mathrm{d}T/\mathrm{d}t + T = T_{\infty}$.
#### Solution
$$\boxed{\tau = \frac{\rho R c}{2h}} \qquad \tau = \frac{2{,}3\times10^3\times5{,}6\times10^{-4}\times1{,}5\times10^3}{2\times100} \approx 9{,}7\ \mathrm{s}$$
Unité : $\dfrac{\mathrm{kg\,m^{-3}\cdot m\cdot J\,kg^{-1}K^{-1}}}{\mathrm{W\,m^{-2}K^{-1}}} = \dfrac{\mathrm{J}}{\mathrm{W}} = \mathrm{s}$.
### Question 15
> [!INDICATION]
> Solution générale = solution particulière constante + solution de l'équation homogène ; la constante se fixe avec $T(0)$.
#### Solution
Solution particulière : $T = T_a$. Solution homogène : $A\,e^{-t/\tau}$. Avec $T(0) = T_i$, $A = T_i - T_a$ :
$$\boxed{T(t) = T_a + (T_i - T_a)\,e^{-t/\tau}}$$
Contrôles : $T(0) = T_i$, $T \to T_a$ quand $t \to \infty$, décroissance puisque $T_i > T_a$.
> [!PIEGE]
> L'exponentielle porte sur l'écart $T - T_a$, pas sur $T$ : $T(t) = T_i\,e^{-t/\tau}$ est faux.
### Question 16
> [!INDICATION]
> Isoler l'exponentielle dans $T(t_1) = T_s$, puis prendre le logarithme.
#### Solution
$$e^{-t_1/\tau} = \frac{T_s - T_a}{T_i - T_a} \qquad \boxed{t_1 = \tau\,\ln\frac{T_i - T_a}{T_s - T_a}}$$
Seuls des écarts de température interviennent, on peut rester en °C : $t_1 = 9{,}66\times\ln\dfrac{1410}{660} = 9{,}66\times0{,}759$, soit $\boxed{t_1 \approx 7{,}3\ \mathrm{s}}$.
> [!PIEGE]
> Un $t_1$ négatif signale une fraction retournée dans le logarithme.
### II.B. Étape de changement d'état
### Question 17
> [!INDICATION]
> Pendant le changement d'état, la température est fixée : que dire du flux échangé ?
> [!INDICATION]
> Égaler l'énergie libérée par la solidification complète à l'énergie évacuée pendant $t_2$.
#### Solution
À $T = T_s$ constante, le flux reçu est constant : $\phi = 2\pi RL\,h\,(T_a - T_s) < 0$. La solidification de la masse $m = \rho\pi R^2L$ libère l'énergie $m\,\ell_{fus}$ (opposé de l'enthalpie de fusion), évacuée vers l'air pendant $t_2$ :
$$-\rho\pi R^2L\,\ell_{fus} = \phi\,t_2 = -2\pi RL\,h\,(T_s - T_a)\,t_2
\qquad \boxed{t_2 = \frac{\rho R\,\ell_{fus}}{2h\,(T_s - T_a)}}$$
$$t_2 = \frac{2{,}3\times10^3\times5{,}6\times10^{-4}\times4{,}0\times10^5}{2\times100\times660} \qquad \boxed{t_2 \approx 3{,}9\ \mathrm{s}}$$
On retrouve $t_{tot} = t_1 + t_2 \approx 11{,}2$ s, la valeur de Q19. Forme équivalente : $t_2 = \tau\,\dfrac{\ell_{fus}}{c\,(T_s - T_a)} = 9{,}7\times0{,}40$.
> [!PIEGE]
> L'équation (2) ne s'applique plus pendant le palier (l'énoncé précise « avant que celui-ci ne change d'état ») : le bilan est ici un simple produit puissance × durée.
### Question 18
> [!INDICATION]
> Deux morceaux : une portion d'exponentielle, puis un palier.
> [!INDICATION]
> Pour faire apparaître $\tau$ : la tangente à l'origine d'une exponentielle coupe son asymptote au bout d'une durée $\tau$.
#### Solution

- Repères : $T_i = 1450$ °C, $T_s = 700$ °C, $T_a = 40$ °C ; $t_1 \approx 7{,}3$ s, $\tau \approx 9{,}7$ s, $t_{tot} \approx 11{,}2$ s.
- Sur $[0\,;t_1]$ : exponentielle décroissante et convexe, qui tendrait vers l'asymptote $T = T_a$.
- La tangente à l'origine, de pente $-(T_i - T_a)/\tau$, coupe l'asymptote en $t = \tau$ ; comme $\tau > t_1$, elle se prolonge au-delà de la fin du refroidissement.
- En $t_1$, la pente vaut environ $-68$ °C·s⁻¹, puis le palier est horizontal jusqu'à $t_{tot}$ : point anguleux.
| $t$ (s) | 0 | 2 | 4 | 6 | 7,3 | 7,3 → 11,2 |
|---|---|---|---|---|---|---|
| $T$ (°C) | 1450 | 1186 | 972 | 798 | 700 | 700 (palier) |
> [!PIEGE]
> Le palier est rigoureusement horizontal, et $\tau$ ne se lit pas en $t_1$ : il faut prolonger la tangente jusqu'à l'asymptote.
### II.C. Résultats et discussion sur le modèle
### Question 19
> [!INDICATION]
> Calculer les rapports de $t_{tot}$ aux deux temps de la partie I.
> [!INDICATION]
> Que signifie une séparation de plusieurs ordres de grandeur entre la mise en forme et la solidification ? Qu'est-ce que la capillarité a le temps de faire ?
#### Solution
$$\frac{t_{tot}}{\tau_a} \approx \frac{11{,}2}{4\times10^{-3}} \approx 3\times10^3, \qquad \frac{t_{tot}}{\tau_c} \approx \frac{11{,}2}{6{,}8\times10^{-3}} \approx 2\times10^3 .$$
- La solidification est environ mille fois plus lente que la mise en forme. L'étirement a lieu alors que la lave est entièrement liquide ; la forme est déterminée par la mécanique avant que la lave ne fige. On peut traiter séparément formation et solidification et, comme dans la partie II, considérer la géométrie fixée pendant le refroidissement.
- En contrepartie, la capillarité dispose de plusieurs milliers de temps $\tau_c$ pour fragmenter le filament avant qu'il ne fige. Si les cheveux existent, c'est grâce à la viscosité de la lave et à sa forte augmentation au refroidissement (Q10).
> [!METHODE]
> Comparer deux temps caractéristiques très différents permet de découpler les phénomènes : pendant le processus lent, le rapide est terminé (ou en équilibre).
### Question 20
> [!INDICATION]
> La diffusion relie une distance et un temps par une loi en $\ell^2 \sim D\,t$.
> [!INDICATION]
> Quelle distance la chaleur doit-elle parcourir pour atteindre la surface d'échange principale ?
#### Solution
La chaleur sort par la surface latérale (Q12) : la distance à parcourir du cœur à la surface est le rayon $R$.
$$\boxed{\tau_{diff} \sim \frac{R^2}{D_{th}}} \qquad \tau_{diff} = \frac{(5{,}6\times10^{-4})^2}{3\times10^{-6}} \approx 0{,}10\ \mathrm{s}$$
$\tau_{diff}/\tau \approx 0{,}01$ : la température s'homogénéise environ cent fois plus vite qu'elle ne varie par les échanges avec l'air. **L'hypothèse de température uniforme est justifiée.**
Remarque (hors exigence) : on retrouve ce résultat avec le nombre de Biot $hR/\lambda$, où $\lambda = \rho c D_{th} \approx 10$ W·m⁻¹·K⁻¹ ; $Bi \approx 5\times10^{-3} \ll 1$.
> [!PIEGE]
> Avec $L = 1$ cm comme distance, on trouverait $33$ s et la conclusion inverse. Le choix de la distance doit être justifié par la géométrie du flux.
### Question 21
> [!INDICATION]
> Reprendre une à une les hypothèses du modèle : modes de transfert, géométrie, coefficients, nature du changement d'état.
> [!INDICATION]
> À 1450 °C, la lave est incandescente. Quel mode de transfert le modèle oublie-t-il ?
#### Solution
- **Rayonnement négligé.** C'est la critique la plus forte : la lave incandescente rayonne. Ordre de grandeur avec la loi de Stefan (émissivité $\approx 0{,}9$) : $\varepsilon\sigma T_i^4 \approx 4\times10^5$ W·m⁻² à 1450 °C, contre $h(T_i - T_a) \approx 1{,}4\times10^5$ W·m⁻² pour la conducto-convection. Le rayonnement diminue fortement en refroidissant ($T^4$), mais en le prenant en compte $t_{tot}$ serait environ deux fois plus court (de l'ordre de 5 s).
- **Géométrie fixe.** Le cylindre continue de s'étirer et de s'affiner ; comme $\tau \propto R$, un vrai cheveu (largeur d'une fraction de millimètre) se solidifie plus vite.
- **$h$ constant.** Il dépend de la vitesse relative de l'air et du rayon, qui varient pendant le vol.
- **Changement d'état idéalisé.** La lave est un mélange de silicates : elle se solidifie sur un intervalle de températures. Les cheveux sont des verres (solide amorphe) : pas de palier net ni de chaleur latente bien définie.
- **Propriétés constantes.** $c$, $\rho$ et surtout la viscosité varient sur 750 °C ; la solidification commence en surface avant le cœur.
- **Air à 40 °C.** Près de l'éruption, l'air et les gaz volcaniques sont bien plus chauds.
> [!JURY]
> Une critique de modèle vaut d'autant plus qu'elle est hiérarchisée, chiffrée quand c'est possible et accompagnée du sens de l'erreur (ici : le modèle surestime $t_{tot}$). Il est inutile de critiquer l'hypothèse de température uniforme, qui vient d'être validée en Q20.
---
## Partie III — Mécanique : transport des cheveux de Pélé
D'abord le tir balistique classique (sans vent, sans frottement), qui montre qu'on ne peut pas expliquer le transport. Puis le vent : nombre de Reynolds, frottement quadratique, résolution numérique par la méthode d'Euler, régime permanent, et confrontation au relief du Piton de la Fournaise.
### III.A. Vent et frottements négligés
### Question 22
> [!INDICATION]
> Système, référentiel, bilan des forces, PFD, projection sur les deux axes.
> [!INDICATION]
> Intégrer deux fois en utilisant les conditions initiales sur la vitesse puis sur la position.
#### Solution
Système : le cheveu, point matériel de masse $m$ ; référentiel terrestre galiléen ; seule force : le poids $-mg\,\vec e_y$. Le PFD donne
$$\frac{\mathrm{d}^2x}{\mathrm{d}t^2} = 0, \qquad \frac{\mathrm{d}^2y}{\mathrm{d}t^2} = -g .$$
Avec $\dot x(0) = v_0\cos\alpha$, $\dot y(0) = v_0\sin\alpha$, puis $x(0) = y(0) = 0$ :
$$\dot x = v_0\cos\alpha,\quad \dot y = -gt + v_0\sin\alpha
\qquad \boxed{x(t) = v_0\cos\alpha\;t \qquad y(t) = -\frac12 g t^2 + v_0\sin\alpha\;t}$$
### Question 23
> [!INDICATION]
> Éliminer le temps entre $x(t)$ et $y(t)$.
#### Solution
Pour $\cos\alpha \neq 0$, $t = x/(v_0\cos\alpha)$ :
$$\boxed{y(x) = x\tan\alpha - \frac{g\,x^2}{2v_0^2\cos^2\alpha}}$$
Parabole tournée vers le bas. Pour $\alpha = \pi/4$ et les valeurs de l'énoncé : $y = x - 0{,}1\,x^2$ (en m).
### Question 24
> [!INDICATION]
> La portée est la racine non nulle de $y(x) = 0$.
> [!INDICATION]
> Faire apparaître $\sin 2\alpha$.
#### Solution
$y(D) = 0$, $D \neq 0$ : $D = \dfrac{2v_0^2\cos^2\alpha\tan\alpha}{g}$, soit
$$\boxed{D = \frac{v_0^2\sin 2\alpha}{g}}$$
$D$ est maximale pour $\sin 2\alpha = 1$ : $\boxed{\alpha = \pi/4}$ et $\boxed{D_{max} = \dfrac{v_0^2}{g} = 10\ \mathrm{m}}$.
**Pas de risque sanitaire** : les retombées restent à une dizaine de mètres du cratère, dans une zone déjà dangereuse et sans cultures, alors que l'introduction parle de kilomètres. Le modèle sans vent ne peut pas expliquer la catastrophe décrite.
> [!PIEGE]
> La question contient quatre demandes : l'expression de $D$, l'angle optimal, la valeur de $D_{max}$ et la discussion sanitaire. La dernière est facile à oublier.
### Question 25
> [!INDICATION]
> Le temps de vol est la racine non nulle de $y(t) = 0$.
> [!INDICATION]
> Pour commenter, comparer à un autre temps du sujet, par exemple celui de la partie II.
#### Solution
$$\boxed{\tau_{vol} = \frac{2v_0\sin\alpha}{g}} \qquad \alpha = \frac{\pi}{4} : \ \tau_{vol} = \frac{2\times10\times\frac{\sqrt2}{2}}{10} = \sqrt2 \approx 1{,}4\ \mathrm{s}$$
Commentaire : le vol est très court, bien plus court que la durée de solidification $t_{tot} \approx 11$ s. Sans vent, la particule retomberait encore liquide, ce qui contredit l'hypothèse « déjà solidifié ». Un vol d'une seconde ne peut pas transporter un cheveu à des kilomètres. De plus, avec $k \simeq 1$ m⁻¹ et $v \sim 10$ m·s⁻¹, le temps caractéristique du freinage par l'air est $1/(kv) \sim 0{,}1$ s, bien inférieur à $\tau_{vol}$ : négliger l'air n'est pas justifié.
### III.B. En présence de vent
### Question 26
> [!INDICATION]
> Le nombre de Reynolds se construit avec les propriétés du fluide qui s'écoule, une vitesse et une longueur.
> [!INDICATION]
> L'air contourne le cylindre perpendiculairement à son axe : quelle dimension « voit »-il ?
#### Solution
$$Re = \frac{\rho_a\,V\,d}{\eta_a}$$
$\rho_a$, $\eta_a$ : masse volumique et viscosité de l'air ; $V \sim u$ : vitesse de l'air par rapport au cylindre ; $d = 2R$ : diamètre, dimension transverse à l'écoulement.
$$Re = \frac{1{,}3\times20\times1{,}12\times10^{-3}}{2\times10^{-5}} \qquad \boxed{Re \approx 1{,}5\times10^3}$$
> [!PIEGE]
> Utiliser $\rho$ et $\eta$ de la lave n'a pas de sens : le Reynolds caractérise l'écoulement de l'air autour du solide.
### Question 27
> [!INDICATION]
> Comparer $Re$ au Reynolds critique et utiliser le rappel de l'énoncé.
#### Solution
$Re \approx 1{,}5\times10^3 \gg Re_c \approx 100$ : la convection l'emporte sur la viscosité, **les frottements sont quadratiques** en la vitesse relative :
$$\vec F = -\frac12\rho_a C_D S_{ref}\,\|\vec v - \vec u\|\,(\vec v - \vec u),$$
ce qui est bien la forme des équations (3). Avec $R$ ou $L$ comme longueur, $Re$ vaut $7\times10^2$ ou $1{,}3\times10^4$ : la conclusion ne change pas.
> [!PIEGE]
> Un grand nombre de Reynolds signifie que la viscosité est négligeable devant l'inertie, et non que le fluide est très visqueux.
### Question 28
> [!INDICATION]
> Méthode d'Euler explicite : nouvelle valeur = ancienne valeur + dérivée × pas. Quelle est la dérivée de $y$ ? celle de $v_y$ ?
> [!INDICATION]
> S'inspirer des lignes 7 et 9, déjà écrites.
#### Solution
$\dot y = v_y$ et $\dot v_y = -g - k\|\vec v - \vec u\|\,v_y$, avec $\|\vec v - \vec u\| = \sqrt{(v_x - u)^2 + v_y^2}$ :
```python
y_new = y+vy*dt
vy_new = vy-g*dt-k*np.sqrt((vx-u)**2+vy**2)*vy*dt
```
> [!PIEGE]
> Le pas `dt` multiplie **tous** les termes de la dérivée, pesanteur comprise. Dans la ligne 10, le facteur final est `vy` (vent horizontal), alors que la norme contient `vx-u`.
### Question 29
> [!INDICATION]
> Où se trouve l'origine du repère ? À quelle surface correspond $y = 0$ ?
> [!INDICATION]
> Le terrain autour d'un volcan est-il horizontal ?
#### Solution
La boucle se poursuit tant que $y \geq 0$, c'est-à-dire tant que le cheveu est au-dessus de l'altitude de son point de lancement, le sommet. Elle s'arrête au premier point calculé d'altitude négative : le sol est modélisé par le plan horizontal passant par le sommet.
Ce choix est discutable : autour du sommet, le terrain descend. Un cheveu qui passe sous $y = 0$ n'a pas atterri, il survole encore le flanc. La simulation s'arrête trop tôt et sous-estime fortement la portée. Il faudrait comparer $y$ à l'altitude du terrain sous le cheveu (profil du document 3), par exemple avec une condition `while y >= profil(x)`. C'est ce que préparent les questions suivantes, en prolongeant la résolution au-delà de $y = 0$.
> [!JURY]
> « Interpréter puis commenter » appelle deux réponses distinctes : traduire la condition en langage physique, puis juger sa pertinence par rapport à la situation réelle.
### Question 30
> [!INDICATION]
> Décrire les deux phases de chaque courbe et l'effet de $u$ sur la seconde.
> [!INDICATION]
> En fin de trajectoire, à quelle vitesse le cheveu descend-il, et à quelle vitesse avance-t-il ? Comparer à l'altitude du sommet.
> [!INDICATION]
> Pour la dispersion, lister ce qui varie d'un cheveu à l'autre ou d'un instant à l'autre.
#### Solution
- **Allure.** Une courte phase de montée (altitude maximale inférieure à 1 m, contre $2{,}5$ m sans frottement), étendue sur quelques mètres, puis une **droite** descendante. La pente de cette droite diminue quand $u$ augmente : pour une même descente de $2{,}35$ m, le cheveu parcourt environ 20 m pour $u = 15$ m·s⁻¹ et 37 m pour $u = 30$ m·s⁻¹. Les conditions initiales sont vite oubliées.
- **Transport kilométrique.** En régime établi, le cheveu avance à la vitesse du vent et tombe lentement, à vitesse constante (environ 3 m·s⁻¹) : il ne perd que 160 m d'altitude par kilomètre parcouru pour $u = 20$ m·s⁻¹. Lancé du sommet à plus de 2,6 km d'altitude, au-dessus d'un relief qui descend lui aussi, il peut parcourir des kilomètres.
- **Dispersion.** La distance parcourue pour une dénivelée donnée vaut $\Delta y\cdot u/\sqrt{g/k} = \Delta y\cdot u\sqrt{k/g}$ : elle est proportionnelle à $u$, qui fluctue (rafales), et croît comme $\sqrt{k}$ avec $k \propto 1/R$, donc les cheveux les plus fins vont plus loin. La direction du vent varie aussi (étalement en azimut : une distance devient une surface), ainsi que les points d'émission et les conditions initiales. Les retombées couvrent donc de grandes surfaces.
> [!PIEGE]
> Les trajectoires ne sont pas des paraboles plus ou moins étalées : leur partie finale est rectiligne, ce que les questions 31 et 32 vont expliquer.
### Question 31
> [!INDICATION]
> Si les vitesses tendent vers des limites, que deviennent leurs dérivées ?
> [!INDICATION]
> Exploiter d'abord l'équation verticale pour montrer que la vitesse relative ne s'annule pas, puis l'équation horizontale.
#### Solution
On admet que $v_x$ et $v_y$ tendent vers des limites finies $v_{x,\infty}$ et $v_{y,\infty}$ ; les dérivées tendent vers $0$. Avec $N = \sqrt{(v_{x,\infty} - u)^2 + v_{y,\infty}^2}$, le système (3) donne
$$0 = -kN\,(v_{x,\infty} - u) \qquad 0 = -g - kN\,v_{y,\infty}.$$
La seconde impose $kN\,v_{y,\infty} = -g \neq 0$ : donc $N \neq 0$ et $v_{y,\infty} < 0$. La première donne alors $v_{x,\infty} = u$. Par suite $N = |v_{y,\infty}| = -v_{y,\infty}$ et la seconde s'écrit $k\,v_{y,\infty}^2 = g$ :
$$\boxed{v_{x,\infty} = u \qquad v_{y,\infty} = -\sqrt{\frac{g}{k}}}$$
C'est un **régime permanent** (de vitesse limite) : le cheveu suit le vent horizontalement et tombe à vitesse constante, son poids étant compensé par la traînée. Numériquement $v_{y,\infty} \approx -3{,}2$ m·s⁻¹ ; le temps d'établissement est de l'ordre de $1/\sqrt{gk} \approx 0{,}3$ s, ce qui explique la brièveté du transitoire du document 2.
> [!PIEGE]
> Oublier la valeur absolue dans $N = |v_{y,\infty}|$ conduit à un mauvais signe : un cheveu qui monterait indéfiniment.
### Question 32
> [!INDICATION]
> La pente d'une courbe paramétrée par le temps est le rapport des dérivées $\dot y/\dot x$.
#### Solution
$$\frac{\mathrm{d}y}{\mathrm{d}x} = \frac{v_y}{v_x}\quad(v_x \to u \neq 0) \qquad \boxed{\left(\frac{\mathrm{d}y}{\mathrm{d}x}\right)_\infty = -\frac{1}{u}\sqrt{\frac{g}{k}}}$$
La pente est constante : la trajectoire tend vers une droite, comme sur le document 2.
| $u$ (m·s⁻¹) | 15 | 20 | 25 | 30 |
|---|---|---|---|---|
| pente $-\sqrt{g/k}/u$ | −0,21 | −0,16 | −0,13 | −0,11 |
Ces valeurs coïncident avec les parties rectilignes du document 2.
### Question 33
> [!INDICATION]
> Lire sur la figure 8 la dénivelée totale et la distance horizontale jusqu'au niveau de la mer.
> [!INDICATION]
> Le cheveu reste en l'air tant qu'il est au-dessus du terrain : comparer deux pentes.
#### Solution
**Pente du flanc.** Sur la figure 8, l'altitude passe de $2{,}6$ km au sommet ($x = 0$) à $0$ pour $x \approx 11{,}8$ km :
$$p \approx \frac{2{,}6}{11{,}8} \approx 0{,}22 \quad(\text{environ } 12^\circ).$$
**Condition.** Partant du sommet, le cheveu s'éloigne durablement du sol si sa trajectoire descend moins vite que le flanc :
$$\frac{1}{u}\sqrt{\frac{g}{k}} < p \iff \boxed{u > u_{min} = \frac{1}{p}\sqrt{\frac{g}{k}}} \qquad u_{min} \approx \frac{3{,}16}{0{,}22} \approx 14\ \mathrm{m\cdot s^{-1}} \approx 50\ \mathrm{km\cdot h^{-1}}.$$

Pour $u = 20$ m·s⁻¹ (pente 0,16 inférieure à 0,22), le cheveu survole le flanc jusqu'à la côte, à plus de 12 km du cratère : il atteint les zones agricoles. Sous 14 m·s⁻¹ environ, il retombe sur les pentes du volcan. Le profil n'est pas de pente uniforme (replat sommital jusqu'à 1,4 km, tronçon plus raide entre 6 et 8 km, pente douce près de la côte) : ce seuil n'est qu'un ordre de grandeur, à 10 ou 20 % près.
> [!PIEGE]
> Le sens de l'inégalité : le cheveu va loin si sa trajectoire est **moins** pentue que le sol, donc si $u$ est **grand**. Trouver $u < u_{min}$ contredirait le document 2.
---
## Bilan
**Méthodes à retenir.**
- Analyse dimensionnelle et lois d'échelle (Q1, Q6, Q20) : construire une longueur ou un temps avec les grandeurs pertinentes, et estimer une dérivée seconde par $L/\tau^2$.
- Système continu (Q4) : découper en $\mathrm{d}m$, intégrer, contrôler par $p = m\,v_G$.
- Bilan thermique d'une phase condensée (Q13, Q17) : premier principe entre $t$ et $t+\mathrm{d}t$, puis bilan global pendant un changement d'état.
- Équation linéaire du premier ordre (Q14 à Q16) : forme canonique, temps caractéristique, solution, durée pour atteindre une valeur.
- Comparaison de temps caractéristiques (Q7, Q19, Q20, Q25) : c'est le fil conducteur du sujet. Dès qu'un temps est calculé, le comparer à un autre et en tirer une conséquence physique.
- Méthode d'Euler (Q28) et régime permanent d'un système non linéaire (Q31) : dérivées nulles, puis résolution dans le bon ordre.
**Prérequis à revoir.** Loi de Newton et signe d'un flux reçu ; solution de $\tau\dot T + T = T_\infty$ ; chaleur latente ; tir balistique ; nombre de Reynolds et types de frottement ; écriture d'un pas d'Euler en Python.
**Questions transposables.** Le refroidissement d'un objet par conducto-convection suivi d'un changement d'état (Q11–Q18) et la validation de l'uniformité de la température par un temps de diffusion (Q20) se retrouvent dans de nombreux sujets. La chute dans un fluide avec vitesse limite (Q31) est un classique. La lecture de la condition d'arrêt d'une boucle de simulation (Q29) est une compétence de plus en plus demandée.
Aide sur le corrigé
Connectez-vous pour demander une explication.