WikiPrépaLivrets

Mines Physique 1 MP 2023Sujet, corrigé et rapport du jury

Fonctions spéciales en physique : W de Lambert, intégrales elliptiques et fonction d'erreur

Pas encore noté
Faisable en Sup

Faisable en sup ?

Faisable en partie en MPSI : 25 questions sur 34, parties I et II en entier. Autres questions faisables en MPSI

Détail par partie
  • Partie I : La fonction de W de LAMBERT15/15Faisable
  • Partie II : L'intégrale elliptique de première espèce7/7Faisable
  • Partie III : La fonction d'erreur de GAUSS : erf(χ)3/12En partie

Téléchargements

Présentation du sujet

Difficulté moyenne
Fonctions spéciales en physique : W de Lambert, intégrale elliptique et fonction erf
Afficher ou masquer la section

Le sujet montre comment des fonctions spéciales permettent de résoudre analytiquement des problèmes physiques. Il comporte trois parties indépendantes : un tir de projectile avec et sans frottements (fonction W de Lambert), le pendule simple aux grands angles (intégrale elliptique de première espèce), puis la diffusion thermique dans un demi-espace appliquée à la formation d'une croûte de lave (fonction erf).

  1. 1Partie I : la fonction W de Lambertpremière annéeChute libre avec vitesse initiale, puis ajout d'un frottement linéaire, schéma d'Euler en Python, et recherche de l'instant d'impact et de la portée à l'aide de la fonction W.
  2. 2Partie II : l'intégrale elliptique de première espècepremière annéePendule simple sans frottement et perte de l'isochronisme des oscillations lorsque l'amplitude n'est plus petite.
  3. 3Partie III : la fonction d'erreur de Gaussdeuxième annéeÉquation de la chaleur dans un demi-espace (problème de Stefan), solution en erf, puis croissance d'une croûte de lave comparée à des données expérimentales.

Difficulté moyenne. Le jury a trouvé la première partie très sélective alors qu'elle était classique, tandis que les questions mathématiques et l'équation de la chaleur ont été globalement bien traitées ; la fin du sujet a été très peu abordée.

L'épreuve en chiffres

Moyenne 11,68 / 20 · écart-type 3,68 · 5 091 présents · où vous situez-vous ?
Afficher ou masquer la section
Moyenne
11,68/ 20
Écart-type
3,68
Présents
5 091
Coefficient
3
Durée
3 h
moyenne 11,6805101520
Deux tiers des copies environ (moyenne ± écart-type)

Votre note sur 20 à ce sujet, en conditions de concours.

Source : document officiel du concours, épreuve du 4 mai 2023. Notes publiées par le concours (après harmonisation le cas échéant). Courbe : estimation par une loi normale.

Ce qu'a observé le jury

6 erreurs relevées
Conditions initiales mal exploitées · Manque de sens physique sur les trajectoires · Schéma d'Euler importé au lieu d'être écrit
Afficher ou masquer la section

Les questions à caractère mathématique ont été globalement bien traitées, et presque toutes les copies savent établir l'équation de la chaleur. En revanche, la résolution d'équations différentielles avec conditions initiales a fortement trié les candidats dans la partie I. Le jury regrette un manque de sens physique sur le tir de projectile.

Les erreurs les plus sanctionnées

  1. 1
    Conditions initiales mal exploitéesQ5

    Beaucoup de candidats n'ont pas su intégrer les équations du mouvement avec frottement ou exploiter les conditions initiales. Une position initiale nulle n'entraîne pas forcément une constante d'intégration nulle.

    « n’implique pas nécessairement une constante d’intégration nulle »
  2. 2
    Manque de sens physique sur les trajectoiresQ6, Q7, Q9

    Des résultats donnant une altitude ou une abscisse négative, ou des trajectoires qui s'arrêtent en plein vol, n'ont pas été remis en question.

    « Les copies manquaient cruellement de sens physique »
  3. 3
    Schéma d'Euler importé au lieu d'être écritQ11

    La question demandait de coder soi-même le schéma d'Euler, et non d'appeler une fonction toute faite.

    « il était attendu une implémentation explicite du schéma d’Euler, et non une importation de la fonction d’une bibliothèque Python »
  4. 4
    Équation de cours fausseQ15

    Une équation censée être parfaitement connue a souvent été écrite avec des erreurs d'homogénéité ou de signe, sanctionnées.

    « De nombreuses inhomogénéités ou des mauvais signes étaient implacablement sanctionnés pour une équation censée être parfaitement connue. »
  5. 5
    Calculs forcés vers le résultatQ17, Q27

    Les démonstrations qui aboutissent au résultat de l'énoncé sans justification ont été sanctionnées, notamment le signe de la dérivée angulaire lors de la prise de racine carrée.

    « Cette question a fait l’objet de nombreuses tentatives d’arnaques, systématiquement sanctionnées. »
  6. 6
    Dimension d'un coefficient de diffusionQ24

    Trouver la dimension de D, pourtant très simple, a posé problème à de nombreux candidats.

Ce qui a été bien réussi

  • L'établissement de l'équation de la chaleur est maîtrisé par une très grande majorité de copies (Q24).
  • Plusieurs questions de la partie I ont été très bien traitées (Q1, Q3, Q4, Q10).
  • Le tracé et la méthode des rectangles semblent maîtrisés (Q18).
  • Les rares copies qui ont repéré la coquille du programme Python de la Q19 ont été valorisées.

Conseils du jury

  • Vérifier le sens physique des résultats (signe, ordre de grandeur, allure de trajectoire) avant de poursuivre.
  • S'entraîner à résoudre des équations différentielles en exploitant soigneusement les conditions initiales.
  • Savoir écrire soi-même un schéma d'Euler et un programme Python élémentaire.
  • Donner chaque application numérique avec une unité et légender les axes des courbes : sinon la réponse vaut 0.
  • Présenter une copie soignée, sans abréviations ambiguës, avec des résultats encadrés.

Synthèse rédigée par WikiPrépa à partir du rapport officiel du jury (à télécharger en PDF). Les citations sont extraites du rapport.

Corrigé détaillé

34 questions
Indice, stratégie puis corrigé complet pour chaque question.
Afficher ou masquer la section

34 questions corrigées une à une : un indice, puis la stratégie, puis le corrigé complet avec le résultat encadré. Chaque question est vérifiée avant publication.

  1. Q1Rappeler la définition d'un référentiel galiléen. Dans quelle mesure le référentiel terrestre peut-il être…
  2. Q2Établir les équations horaires du mouvement. Montrer que le mouvement est plan.
  3. Q3Établir l'équation de la trajectoire. Quelle est la forme de la trajectoire ? Est-elle symétrique ?
  4. Q4Déterminer les coordonnées du sommet de la trajectoire. Définir la portée du tir et établir son expression…
  5. Q5Quelle est la dimension du coefficient ? Définir à partir de un temps caractéristique. Le mouvement…
  6. Q6Établir, en fonction,,, et, les nouvelles équations horaires du mouvement.
  7. Q7Dans la situation où, simplifier les équations horaires de la trajectoire et donner l'allure du mouvement.
  8. Q8Dans la situation où, simplifier les équations horaires du mouvement en faisant apparaitre une vitesse…
Voir les 26 autres questions
  1. Q9Déduire des résultats précédents, l'allure globale de la trajectoire dans une situation où le temps de vol…
  2. Q10Tracer l'allure de la trajectoire pour un temps de vol de l'ordre de.
  3. Q11Dresser le tableau de variation de la fonction et déterminer la valeur de son minimum global. La fonction W…
  4. Q12On peut montrer que:. Quelle est la valeur de ? On souhaite appliquer le schéma d'Euler explicite avec un…
  5. Q13En déduire à quel instant le projectile touche le sol. On posera.
  6. Q14On rappelle que par définition où est la fonction identité:. En déduire que la portée est donnée par.
  7. Q15À l'aide de la figure 2, déterminer la valeur numérique de l'angle assurant la portée maximale pour, et.
  8. Q16Établir l'équation différentielle du mouvement vérifiée par la fonction.
  9. Q17On fait l'approximation des petits angles tels que. Établir dans ces conditions la période des oscillations…
  10. Q18Déterminer l'expression générale de sans faire l'approximation des petits angles. En déduire que la période…
  11. Q19Après avoir tracé le graphe de la fonction pour, illustrer le principe de la méthode des rectangles médians…
  12. Q20Recopier et compléter le code suivant permettant de calculer par la méthode des rectangles médians.
  13. Q21En utilisant la figure 4, pour un pendule tel que, évaluer lorsque. Quel est le décalage temporel induit par…
  14. Q22Dans quelle situation courante rencontre-t-on la cycloïde ?
  15. Q23Comment explique-t-on de nos jours le résultat erroné obtenu par Lord KELVIN ?
  16. Q24Rappeler la définition du vecteur, densité de flux thermique. Quelle est sa dimension ? Rappeler la loi de…
  17. Q25Quelle est l'énergie thermique reçue par cette tranche entre et ? Pourquoi étudie-t-on une tranche «…
  18. Q26Quelle est l'équation vérifiée par ? Déterminer les valeurs de, et.
  19. Q27Montrer que
  20. Q28En utilisant la fonction, montrer que. On donne. En déduire une expression de faisant apparaître une…
  21. Q29En reprenant l'équation de la question 27, montrer que
  22. Q30Quelle est l'énergie libérée par la solidification à la température d'une tranche de lave de surface en…
  23. Q31Toute l'énergie libérée par la solidification doit être évacuée par diffusion dans la lave solide car la…
  24. Q32En déduire que
  25. Q33Quel algorithme peut on utiliser pour obtenir la constante numériquement ? Expliquer en quelques mots son…
  26. Q34À l'aide de la figure 7, estimer la valeur numérique de. En déduire l'épaisseur de la croûte de lave six…

Description

Le sujet montre comment des fonctions spéciales permettent de résoudre analytiquement des problèmes physiques. Il comporte trois parties indépendantes : un tir de projectile avec et sans frottements (fonction W de Lambert), le pendule simple aux grands angles (intégrale elliptique de première espèce), puis la diffusion thermique dans un demi-espace appliquée à la formation d'une croûte de lave (fonction erf).

Ces sujets peuvent vous intéresser

Lecture du sujet en ligne

L'énoncé complet, avec les formules et les figures, sans ouvrir le PDF.
Afficher ou masquer la section

ÉCOLE DES PONTS PARISTECH, ISAE-SUPAERO, ENSTA PARIS, TÉLÉCOM PARIS, MINES PARIS, MINES SAINT-ÉTIENNE, MINES NANCY, IMT ATLANTIQUE, ENSAE PARIS, CHIMIE PARISTECH - PSL.
Concours Mines-Télécom, Concours Centrale-Supélec (Cycle International).

CONCOURS 2023
PREMIÈRE ÉPREUVE DE PHYSIQUE
Durée de l'épreuve : 3 heures
L'usage de la calculatrice ou de tout dispositif électronique est interdit.
Les candidats sont priés de mentionner de façon apparente
sur la première page de la copie :
PHYSIQUE I - MP
L'énoncé de cette épreuve comporte 7 pages de texte.
Si, au cours de l'épreuve, un candidat repère ce qui lui semble être une erreur d'énoncé, il le signale sur sa copie et poursuit sa composition en expliquant les raisons des initiatives qu'il est amené à prendre.

Fonctions spéciales

Ce sujet comporte trois parties indépendantes.
Bon nombre de problèmes rencontrés en physique peuvent être résolus à l'aide de «fonctions spéciales ». Ces fonctions définies mathématiquement sont implémentées dans de nombreuses bibliothèques informatiques (comme scipy) et peuvent être utilisées aussi simplement qu'une fonction sinus ou racine carrée qui sont elles aussi d'une certaine manière des fonctions spéciales et tout aussi analytiques . . .
On rencontre bien souvent des résolutions numériques de problèmes physiques alors que l'utilisation de ces fonctions spéciales permet une résolution complète et analytique. Ce problème se propose d'illustrer l'intérêt de ces «fonctions spéciales».

I La fonction de W de Lambert

I.A Tir d'un projectile sans frottements

Un projectile assimilé à un point matériel de masse m est lancé à partir du sol en O avec une vitesse initiale v⃗_0 ∈ (O, u⃗_y, u⃗_z) et faisant un angle θ_0 avec l'horizontale dans le référentiel terrestre supposé galiléen.
    1. Rappeler la définition d'un référentiel galiléen. Dans quelle mesure le référentiel terrestre peut-il être supposé galiléen?
    1. Établir les équations horaires du mouvement.
Montrer que le mouvement est plan.
    1. Établir l'équation de la trajectoire. Quelle est la forme de la trajectoire? Est-elle symétrique?
Figure 1 - Tir d'un projectile
    1. Déterminer les coordonnées du sommet S de la trajectoire. Définir la portée ℓ du tir et établir son expression. Quel est l'angle θ_0 assurant un tir de portée maximale?

I.B Tir d'un projectile avec frottements

On considère maintenant que le projectile est soumis à une force de frottements proportionnelle à la vitesse: f⃗ = − αv⃗ avec α > 0.
    1. Quelle est la dimension du coefficient α ? Définir à partir de α un temps caractéristique τ. Le mouvement reste-t-il plan?
      ◻ − 6. Établir, en fonction g, τ, v_0 = ‖v⃗_0‖, θ_0 et t, les nouvelles équations horaires du mouvement.
    1. Dans la situation où t≪τ, simplifier les équations horaires de la trajectoire et donner l'allure du mouvement.
    1. Dans la situation où t≫τ, simplifier les équations horaires du mouvement en faisant apparaitre une vitesse limite v_∞.
      Où retombe le projectile?
    1. Déduire des résultats précédents, l'allure globale de la trajectoire dans une situation où le temps de vol est grand devant τ, en séparant la trajectoire en trois phases.
    1. Tracer l'allure de la trajectoire pour un temps de vol de l'ordre de τ.

I.C La portée maximale d'un tir avec frottement

    1. Dresser le tableau de variation de la fonction T : χ ↦ T(χ) = χe^χ et déterminer la valeur β de son minimum global.
      La fonction W de Lambert est définie comme étant la fonction réciproque de T sur [β, + ∞[. Reproduire le graphe de T représenté sur la partie gauche de la figure 2 et expliquer comment en déduire l'allure de W représenté sur la partie droite.
Figure 2 - Représentations graphiques de T(χ) (à gauche) et W(χ) (à droite)
    1. On peut montrer que : (χ + exp[W(χ)])W^′(χ) = 1. Quelle est la valeur de W(0) ?
On souhaite appliquer le schéma d'Euler explicite avec un pas h = 0.0001 pour résoudre cette équation différentielle. Donner le code python permettant d'obtenir une représentation graphique de W(χ) sur l'intervalle [0; 2, 5[.
La fonction W(χ) est implémentée dans scipy. On peut l'appeler avec : from scipy.special import lambertw.
On montre que si ad ≠ 0, la solution de l'équation at + b + ce^(dt) = 0 pour l'inconnue t est donnée par l'expression
t = − b/a − 1/d W((cd)/aexp(− (bd)/a)).
◻ − 13. En déduire à quel instant t^∗ > 0 le projectile touche le sol. On posera u = − (1 + (v_0 sinθ_0)/(gτ)).
◻ − 14. On rappelle que par définition Wexp(W) = Id où Id est la fonction identité : χ ↦ χ.
En déduire que la portée est donnée par ℓ = τv_0 cosθ_0(1 − W(ue^u)/u).
En posant γ = v_0/v_∞, on montre que l'angle initial donnant la portée maximale est :
θ_(max) = {arcsin((γW((γ^2 − 1)/e))/(γ^2 − 1 − W((γ^2 − 1)/e))), si γ ≠ 1; arcsin(1/(e − 1)) ≃ 35, 6^∘, si γ = 1
    1. À l'aide de la figure 2, déterminer la valeur numérique de l'angle assurant la portée maximale pour v_0 = 10 m ⋅ s^(− 1), g = 9, 8 m ⋅ s^(− 2) et τ = 0, 4 s.

II L'intégrale elliptique de première espèce

Dans toute cette partie on néglige les frottements de l'air.
On étudie un pendule simple constitué d'une masse ponctuelle m et d'une tige rigide de longueur ℓ et de masse négligeable, astreint à évoluer dans un plan vertical ( O, u⃗_x, u⃗_y ).
On repère sa position par l'angle θ(t). À t = 0 on lâche le pendule sans vitesse initiale avec θ(t = 0) = θ_0 ∈ ]0, π/2[.
    1. Établir l'équation différentielle du mouvement vérifiée par la fonction θ(t).
    1. On fait l'approximation des petits angles tels que sinθ ∼ θ.
Figure 3 - Pendule simple
Établir dans ces conditions la période T_0 des oscillations. Quelle est la propriété remarquable de la période dans le cadre de cette approximation?
    1. Déterminer l'expression générale de (dθ)/(dt) sans faire l'approximation des petits angles.
En déduire que la période T des oscillations du pendule est donnée par :
T = (2T_0)/π∫_0^(θ_0)(dθ)/(√(2(cosθ − cosθ_0))).
La propriété remarquable de la question précédente est-elle conservée?
En effectuant le changement de variable sinθ/2 = sinφsin(θ_0)/2, on montre que :
T = (2T_0)/πK(sin^2(θ_0)/2) avec K(χ) = ∫_0^(π/2)(dφ)/(√(1 − χsin^2 φ)).
On souhaite calculer l'intégrale K(χ) par la méthode des rectangles médians pour un angle θ_0 = π/3.
    1. Après avoir tracé le graphe de la fonction χ ↦ 1 + √χ pour χ ∈ [0; 9], illustrer le principe de la méthode des rectangles médians pour calculer le réel I = ∫_0^9(√χ + 1)dχ en utilisant 9 rectangles.
      Si on double le nombre de rectangles utilisés qu'en est-il de la différence entre la valeur exacte de I et la valeur approchée numériquement par la méthode des rectangles médians?
    1. Recopier et compléter le code suivant permettant de calculer K(χ) par la méthode des rectangles médians.
import math as m
def f(x,phi):
    return
S = 0.
N = 100
a = 0.
b = m.pi/2.
pas = ............
theta_0 = m.pi/3.
x = m.sin(theta_0)**2
for i in range(N):
    phi =
    S =
print(pas * S)
La fonction χ ↦ K(χ) est nommée intégrale elliptique complète de première espèce. Elle est implémentée dans scipy. On peut l'appeler directement avec : from scipy.special import ellipk.
◻ − 21. En utilisant la figure 4 , pour un pendule tel que T_0 = 1 s, évaluer T lorsque θ_0 = 50^∘. Quel est le décalage temporel induit par la prise en compte de l'approximation des petits angles si l'on envisage de mesurer une heure?
Au XVII ^e siècle les puissances maritimes désiraient posséder des instruments précis pour la mesure du temps afin de faciliter la navigation (notamment pour déterminer la longitude). Les rois de France et d'Angleterre avaient offert des prix importants à qui serait capable de réaliser un chronomètre précis, fiable et utilisable en mer.
Figure 4 − θ_0 ↦ T(θ_0)/T_0
◻ − 22. Dans quelle situation courante rencontre-t-on la cycloïde?

III La fonction d'erreur de GAUSS : erf(χ)

III.A Introduction au problème de Stefan

Un certain nombre de problèmes géologiques importants peuvent être modélisés par le chauffage ou le refroidissement instantané d'un demi-espace semi-infini. Au milieu du XIX ^e siècle Lord Kelvin a ainsi utilisé cette idée pour estimer l'âge de la Terre. Il supposa qu'à la surface le flux d'énergie thermique résultait du refroidissement d'un flux initialement chaud de la Terre et a conclu que l'âge de la Terre était environ 65 millions d'années. On retrouve ces phénomènes en étudiant le refroidissement de la lithosphère océanique ou l'évolution d'une coulée de magma.
    1. Comment explique-t-on de nos jours le résultat erroné obtenu par Lord Kelvin?
On étudie un milieu matériel semi-infini défini par y > 0 dont la surface subit un changement instantané de température. Initialement à t = 0^−, le demi-espace est à la température uniforme T_1; pour t > 0, la surface y = 0 est maintenue à une température constante T_0. Si T_1 > T_0, le milieu matériel se refroidit et sa température diminue. La situation est représentée à la figure 5 pour le cas T_1 > T_0.
Figure 5 - Évolution de la température
Le flux thermique élémentaire, défini comme la quantité d'énergie traversant une surface élémentaire dS pendant dt, est noté δφ_Q.
    1. Rappeler la définition du vecteur j⃗_Q, densité de flux thermique. Quelle est sa dimension?
Rappeler la loi de Fourier, ainsi que ses conditions d'application.
En déduire la dimension de la conductivité thermique κ.
On étudie une tranche mésoscopique de sol de masse m de masse volumique ρ et de capacité thermique massique c comprise entre y et y + dy de surface S.
    1. Quelle est l'énergie thermique δQ reçue par cette tranche entre t et t + dt ?
Pourquoi étudie-t-on une tranche «mésoscopique»?
Établir l'expression de sa variation d'énergie interne dU en fonction de (∂j_Q)/(∂y), S, dy et dt puis en fonction de ρ, c, S, (∂T)/(∂t), dy et dt.
En déduire l'équation de la chaleur à une dimension (∂T)/(∂t) = D(∂^2 T)/(∂y^2) dans laquelle on précisera l'expression et la dimension du coefficient D de diffusion thermique.
En déduire l'expression d'une longueur caractéristique L en fonction de D et du temps t.
On introduit la température adimensionnée
θ(y, t) = (T(y, t) − T_1)/(T_0 − T_1).
    1. Quelle est l'équation vérifiée par θ(y, t) ?
Déterminer les valeurs de θ(y > 0, t = 0), θ(y = 0, t > 0) et θ(y → + ∞, t > 0).
On introduit une variable de similarité sans dimension η = y/(2√(Dt)) et on suppose que θ n'est une fonction que de cette seule variable η.
◻ − 27. Montrer que
(d^2 θ(η))/(dη^2) + 2η(dθ(η))/(dη) = 0.
◻ − 28. En utilisant la fonction φ(η) = (dθ(η))/(dη), montrer que θ(η) = 1 − 2/(√π)∫_0^η e^(− z^2) dz.
On donne ∫_0^(+ ∞)e^(− z^2) dz = (√π)/2. En déduire une expression de T(y, t) faisant apparaître une intégrale.
La fonction χ ↦ 2/(√π)∫_0^χ e^(− z^2) dz est appelée fonction d'erreur de GAUSS, elle est implémentée dans scipy.
Elle est souvent notée erf(χ). On peut l'appeler directement en utilisant la commande : from scipy.special import erf.

III.B Formation d'une croûte de lave solide.

Dans cette dernière partie on s'intéresse à une coulée de lave en fusion et à la formation d'une croûte solide à sa surface. On étudie alors l'augmentation de l'épaisseur de cette croûte en fonction du temps.
À la surface extérieure, en y = 0, la lave est en contact avec l'air à la température constante T_0. La lave en fusion à la température T_f est donc soudainement portée à la température T_0 à t = 0. Dans ces conditions, la couche superficielle de la lave se solidifie, et on note y_s(t) l'épaisseur de la couche de lave solide.
Nous devons donc résoudre l'équation de la chaleur dans l'espace 0 ⩽ y ⩽ y_s(t) avec comme conditions aux limites T = T_0 en y = 0, et T = T_f en y = y_s(t), et comme condition initiale y_s = 0 à t = 0.
Figure 6 - Formation d'une croûte de lave solide
La position y_s(t) de l'interface de transition de phase est une fonction a priori inconnue du temps. Comme dans la situation précédente il n'y a pas d'échelle de longueur définie dans ce problème. Pour cette raison, on travaillera également avec la variable de similarité sans dimension η = y/(2√(Dt)).
On utilisera également la température adimensionnée
θ(y, t) = (T(y, t) − T_0)/(T_f − T_0)
La profondeur de l'interface de solidification y_s(t) doit enfin s'adapter à la longueur caractéristique de la diffusion thermique. Nous supposerons que celle-ci varie proportionnellement à la racine carrée du temps, de telle sorte que : η_s = (y_s(t))/(2√(Dt)) = cte = λ. Cette constante est inconnue et reste à déterminer.
◻ − 29. En reprenant l'équation de la question 27, montrer que
θ(η) = (erf(η))/(erf(λ)).
Afin d'obtenir l'expression puis la valeur de la constante λ, nous allons étudier la solidification d'une tranche de lave d'épaisseur d y_s entre les instants t et t + dt
    1. Quelle est l'énergie δQ libérée par la solidification à la température T_f d'une tranche dy_s de lave de surface S en fonction de la masse volumique ρ de la lave en fusion et l'enthalpie de fusion massique : Δh_(sol → liq).
    1. Toute l'énergie libérée par la solidification doit être évacuée par diffusion dans la lave solide car la lave en fusion reste à la température T_f. Montrer que :
ρΔh_(sol → liq)(T_f)(dy_s(t))/(dt) = κ((∂T)/(∂y))_(y = y_s)
Figure 7 - Graphe de λ ↦ (exp(− λ^2))/(λerf(λ))
    1. En déduire que
(exp(− λ^2))/(λerf(λ)) = (√π)/(c(T_f − T_0))Δh_(sol → liq)(T_f).
    1. Quel algorithme peut on utiliser pour obtenir la constante λ numériquement?
      Expliquer en quelques mots son fonctionnement.
On donne les valeurs numériques suivantes :
- Δh_(sol → liq)(T_f) = 400 kJ ⋅ kg^(− 1), - ρ = 2600 kg ⋅ m^(− 3); - c = 1 kJ ⋅ kg^(− 1) ⋅ K^(− 1), - D = 7 × 10^(− 7)SI; - T_f − T_0 = 1000 K, - √π ∼ 1, 77
    1. À l'aide de la figure 7, estimer la valeur numérique de λ.
En déduire l'épaisseur de la croûte de lave six mois après l'éruption.
Comparer votre résultat à ceux de la figure 8 tirés d'une expérience ^1.
Figure 8 - Épaisseurs des croûtes de lave solides à la surface des lacs de lave dans les trois cratères à fosse Kilauea lki (1959), Alae (1963) et Makaopuhi (1965) sur le volcan Kilauea, Hawaii (Wright et al., 1976), et résultat théorique.

FIN DE L'ÉPREUVE


    1. Wright, T. L., Peck, D. L., and Shaw, H. R. (1976). Kilauea lava lakes : Natural laboratories for study of cooling, crystallization, and differentiation of basaltic magma. In The Geophysics of the Pacific Ocean Basin and its Margin, eds. G. H. Sutton, M. H. Manghnani, R. Moberly, and E. U. McAfee, vol. 19 of Geophysical Monograph Series, Washington, D.C. : American Geophysical Union, pp. 375-90

Questions fréquentes

4 questions
Sur quels chapitres porte le sujet Mines physique 1 MP 2023 ?
Afficher ou masquer la section

Sur quels chapitres porte le sujet Mines physique 1 MP 2023 ?

Sur la mécanique du point (tir avec et sans frottements), le pendule simple aux grands angles et la diffusion thermique avec changement d'état. Le fil conducteur est l'usage de fonctions spéciales : W de Lambert, intégrale elliptique et erf.

Quelles erreurs le jury des Mines a-t-il relevées en physique 1 MP 2023 ?

Des conditions initiales mal exploitées (Q5), un manque de sens physique sur les trajectoires, un schéma d'Euler importé d'une bibliothèque au lieu d'être codé (Q11) et des démonstrations forcées vers le résultat (Q17, Q27).

Le sujet Mines physique 1 MP 2023 est-il faisable en première année ?

Les parties I et II portent sur la mécanique du point et le pendule simple, abordables dès la première année. La partie III repose sur la diffusion thermique, vue en deuxième année.

Le sujet Mines physique 1 MP 2023 est-il difficile ?

Le jury a trouvé la partie I très sélective, surtout pour la résolution des équations différentielles. Les questions mathématiques et l'équation de la chaleur ont en revanche été bien réussies.

Pas de description pour le moment