La physique pas à pas

Mécanique, traînée, oscillations, lagrangiens, électromagnétisme et plasma, thermodynamique, fluides et la relativité dont le GPS a besoin : démonstrations, exemples résolus et exercices, du premier cycle au doctorat.

Constantes utilisées sur cette page

c=299 792 458 m/sc=299\,792\,458\ \text{m/s}
e=1.602 176 634×10−19 Ce=1.602\,176\,634\times10^{-19}\ \text{C}
kB=1.380 649×10−23 J/Kk_B=1.380\,649\times10^{-23}\ \text{J/K}
μ0=1.256 637 062 12×10−6 N/A2\mu_0=1.256\,637\,062\,12\times10^{-6}\ \text{N/A}^2
GM⊕=3.986 004 418×1014 m3/s2GM_\oplus=3.986\,004\,418\times10^{14}\ \text{m}^3/\text{s}^2
g0=9.80665 m/s2g_0=9.80665\ \text{m/s}^2

Constantes : CODATA 2018 (c, e, m_p, m_e, ε₀, μ₀, k_B et N_A tels que publiés). Chaque nombre de cette page est calculé par le script de génération du site.

Essayez dans le Physics Lab

Lois de Newton et conservation

  • Premier cycle (DUT/BTS)
  • Licence (L1–L3)

Parcours : physique · Environ 25 min · Nombres vérifiés à la génération:

Vous saurez

Énoncer les lois de Newton ; établir la conservation de la quantité de mouvement et de l’énergie ; résoudre un choc élastique à 1 D.

Prérequis

Algèbre, vecteurs, dérivées. La physique expliquée

D’abord l’intuition

Les forces changent le mouvement. Si aucune force extérieure n’agit sur un groupe d’objets, leur quantité de mouvement totale ne peut pas changer ; et si aucune énergie n’est perdue en chaleur ou en son, leur énergie cinétique totale non plus.

🧒 Explique-moi comme si j’avais 10 ans

Pousse un chariot vide : il accélère ; pousse un chariot plein : il accélère moins. Et quand deux autos tamponneuses se cognent, l’élan que l’une perd, l’autre le gagne : le total ne change pas.

Pas à pas

  1. F=dpdt,p=mv\mathbf F=\frac{d\mathbf p}{dt},\qquad \mathbf p=m\mathbf v
    Deuxième loi sous forme de quantité de mouvement.
  2. F12=−F21  ⇒  ddt(p1+p2)=0\mathbf F_{12}=-\mathbf F_{21}\;\Rightarrow\;\frac{d}{dt}(\mathbf p_1+\mathbf p_2)=0
    Troisième loi : les forces internes s’annulent, la quantité de mouvement totale est conservée.
  3. W=∫F⋅dr=Δ(12mv2)W=\int\mathbf F\cdot d\mathbf r=\Delta\left(\tfrac12mv^2\right)
    Théorème de l’énergie cinétique.
  4. v1=(m1−m2)u1+2m2u2m1+m2,v2=(m2−m1)u2+2m1u1m1+m2v_1=\frac{(m_1-m_2)u_1+2m_2u_2}{m_1+m_2},\qquad v_2=\frac{(m_2-m_1)u_2+2m_1u_1}{m_1+m_2}
    On résout ensemble la conservation de la quantité de mouvement et de l’énergie cinétique pour un choc élastique.

Exemple résolu

Données
Chariot 1 : m₁ = 2 kg à u₁ = 3 m/s heurte le chariot 2 : m₂ = 1 kg au repos. Choc élastique sur rail horizontal sans frottement.
Inconnue
Vitesses finales v₁ et v₂.
Stratégie
Appliquer les formules du choc élastique, puis vérifier quantité de mouvement et énergie.
Étapes
  1. v₁ = (2 − 1)·3 / 3 = 1,00 m/s
  2. v₂ = 2·2·3 / 3 = 4,00 m/s
  3. p: 2·3 = 6 → 2·1 + 1·4 = 6 kg·m/s
  4. KE: 9,0 J → 9,0 J
Réponse
v₁ = 1,0 m/s, v₂ = 4,0 m/s (tous deux vers l’avant).
Vérification
Quantité de mouvement (6 kg·m/s) et énergie cinétique (9 J) inchangées, et la vitesse relative s’inverse (3 m/s avant, 3 m/s après), comme il se doit pour un choc élastique.
Portée
La même algèbre décrit une assistance gravitationnelle d’une sonde, vue depuis le référentiel de la planète.
Voir le code Python
m1, u1, m2, u2 = 2.0, 3.0, 1.0, 0.0
v1 = ((m1-m2)*u1 + 2*m2*u2) / (m1+m2)
v2 = ((m2-m1)*u2 + 2*m1*u1) / (m1+m2)
print(v1, v2, m1*v1 + m2*v2, 0.5*m1*v1**2 + 0.5*m2*v2**2)

Exercices corrigés automatiquement

  1. m/s

    Solution complète

    v₂ = 2·1·4/4 = 2,0 m/s ; le chariot léger repart en arrière à -2,0 m/s.

Exercices

  1. Une balle tombe de 20 m sans vitesse initiale (sans air). Quelle est sa vitesse en bas ?

    Voir la solution

    mgh = ½mv², donc v = √(2gh) = 19,81 m/s, indépendamment de la masse.

Erreur fréquente

Traiter la quantité de mouvement comme un nombre et non un vecteur : un chariot qui recule a une quantité de mouvement négative.

Références

  • D. Halliday, R. Resnick, J. Walker, Fundamentals of Physics, 12th ed., Wiley, 2021.
  • J. R. Taylor, Classical Mechanics, University Science Books, 2005.
  • MIT OpenCourseWare, 8.01 Classical Mechanics and 8.02 Electricity and Magnetism.

Projectiles avec traînée (EDO numérique)

  • Licence (L1–L3)
  • Master

Parcours : physique · Environ 40 min · Nombres vérifiés à la génération:

Vous saurez

Écrire les équations du mouvement avec traînée quadratique, voir pourquoi elles n’ont pas de solution analytique, et les intégrer par Runge–Kutta d’ordre 4.

Prérequis

Deuxième loi de Newton, vecteurs, Python de base. La physique expliquée

D’abord l’intuition

La traînée croît comme le carré de la vitesse et s’oppose toujours à elle. Cela couple les mouvements horizontal et vertical : on avance la solution par petits pas de temps sur ordinateur.

🧒 Explique-moi comme si j’avais 10 ans

Lance une balle : l’air la freine sans arrêt. Elle va moins loin et moins haut que dans l’espace, et elle redescend plus raide qu’elle n’est montée.

Pas à pas

  1. mv˙=mg−12ρCDA ∣v∣ vm\dot{\mathbf v}=m\mathbf g-\tfrac12\rho C_DA\,|\mathbf v|\,\mathbf v
    Traînée quadratique, opposée à la vitesse.
  2. v˙x=−k∣v∣vx,v˙y=−g−k∣v∣vy,k=ρCDA2m\dot v_x=-k|\mathbf v|v_x,\quad \dot v_y=-g-k|\mathbf v|v_y,\qquad k=\frac{\rho C_DA}{2m}
    Couplées par |v| : pas de solution élémentaire.
  3. vt=2mgρCDAv_t=\sqrt{\frac{2mg}{\rho C_DA}}
    Vitesse limite : la traînée équilibre le poids.
  4. sn+1=sn+Δt6(k1+2k2+2k3+k4)\mathbf s_{n+1}=\mathbf s_n+\tfrac{\Delta t}{6}(\mathbf k_1+2\mathbf k_2+2\mathbf k_3+\mathbf k_4)
    RK4 : l’erreur globale décroît comme Δt⁴.

Exemple résolu

Données
Balle générique : m = 0,145 kg, diamètre 0,074 m, C_D = 0,35, ρ = 1,225 kg/m³ ; lancée à 40 m/s (131 ft/s), à 35° au-dessus d’un sol plat.
Inconnue
Portée, hauteur maximale et durée de vol, avec et sans air.
Stratégie
Formule exacte dans le vide ; RK4 avec Δt = 1 ms avec traînée, puis Δt divisé par deux pour vérifier la convergence.
Étapes
  1. k = ρC_DA/(2m) = 0,00636 1/m
  2. vacuum: R = v₀² sin 70° / g = 153,3 m, h = 26,8 m, t = 4,68 s
  3. RK4 (Δt = 1 ms): R = 93,9 m, h = 20,2 m, t = 4,05 s
  4. RK4 (Δt = 0.5 ms): R = 93,9 m
  5. v_t = 39,3 m/s
Réponse
Avec traînée : portée 93,9 m (contre 153,3 m), apogée 20,2 m (contre 26,8 m), durée 4,05 s (contre 4,68 s).
Vérification
Diviser le pas par deux change la portée de moins de 1 mm : l’intégration a convergé. La vitesse de lancement vaut à peu près la vitesse limite (39,3 m/s) : au départ la traînée est aussi grande que le poids, une forte réduction est donc attendue.
Portée
Le même intégrateur, avec poussée et densité variable en plus, sert à calculer les trajectoires des fusées et des capsules en rentrée.
Voir le code Python
import math
m, d, cd, rho, g = 0.145, 0.074, 0.35, 1.225, 9.80665
k = 0.5*rho*cd*math.pi*d*d/4/m
def f(s):
    x, y, vx, vy = s; v = math.hypot(vx, vy)
    return [vx, vy, -k*v*vx, -g - k*v*vy]
s, dt = [0, 0, 40*math.cos(math.radians(35)), 40*math.sin(math.radians(35))], 1e-3
while s[1] >= 0:
    a = f(s); b = f([s[i]+dt/2*a[i] for i in range(4)])
    c = f([s[i]+dt/2*b[i] for i in range(4)]); e = f([s[i]+dt*c[i] for i in range(4)])
    s = [s[i]+dt/6*(a[i]+2*b[i]+2*c[i]+e[i]) for i in range(4)]
print(s[0])

Exercices corrigés automatiquement

  1. m

    Solution complète

    R(45°) = 94,8 m, contre 93,9 m à 35° et 163,2 m dans le vide à 45°.

Exercices

  1. Pourquoi la balle redescend-elle plus raide qu’elle n’est montée ?

    Voir la solution

    La traînée retire de la vitesse horizontale tout le temps, tandis que la gravité redonne de la vitesse verticale à la descente : l’angle de descente est plus raide.

Erreur fréquente

Appliquer la traînée séparément en −k v_x² et −k v_y². La traînée dépend de la vitesse totale |v|, pas de chaque composante seule.

Références

  • J. R. Taylor, Classical Mechanics, University Science Books, 2005.
  • W. H. Press et al., Numerical Recipes, 3rd ed., Cambridge University Press, 2007.
  • F. M. White, Fluid Mechanics, 8th ed., McGraw-Hill, 2016.

Oscillateurs : simple, amorti et valeurs propres

  • Licence (L1–L3)
  • Master

Parcours : physique · Environ 35 min · Nombres vérifiés à la génération:

Vous saurez

Résoudre l’oscillateur harmonique amorti par son équation caractéristique ; relier les valeurs propres à la pulsation propre, à l’amortissement et au facteur de qualité ; faire le lien avec la stabilité des avions.

Prérequis

EDO linéaires à coefficients constants, nombres complexes. La physique expliquée

D’abord l’intuition

Presque tout système près de l’équilibre se comporte comme une masse sur un ressort avec un peu de frottement : ponts, circuits, molécules, avions. Deux nombres le résument : la pulsation propre et l’amortissement.

🧒 Explique-moi comme si j’avais 10 ans

Une balançoire va et vient toute seule. Si tu arrêtes de pousser, l’air et la chaîne qui grince lui prennent peu à peu son énergie et elle s’arrête. Les ingénieurs mesurent à quelle vitesse elle oscille et à quelle vitesse elle se calme.

Pas à pas

  1. mx¨+cx˙+kx=0m\ddot x+c\dot x+kx=0
    Loi de Newton avec force de rappel et amortissement visqueux.
  2. x=est  ⇒  s2+2ζωns+ωn2=0,ωn=k/m, ζ=c2kmx=e^{st}\;\Rightarrow\;s^2+2\zeta\omega_ns+\omega_n^2=0,\quad \omega_n=\sqrt{k/m},\ \zeta=\frac{c}{2\sqrt{km}}
    Équation caractéristique.
  3. s=−ζωn±i ωn1−ζ2(ζ<1)s=-\zeta\omega_n\pm i\,\omega_n\sqrt{1-\zeta^2}\quad(\zeta<1)
    Sous-amorti : paire complexe ; partie réelle = taux de décroissance, partie imaginaire = pulsation amortie.
  4. x˙=(01−ωn2−2ζωn)x,Q=12ζ\dot{\mathbf x}=\begin{pmatrix}0&1\\-\omega_n^2&-2\zeta\omega_n\end{pmatrix}\mathbf x,\qquad Q=\frac{1}{2\zeta}
    Même chose sous forme de système du premier ordre : les valeurs propres de la matrice sont les racines s. C’est exactement l’analyse du mode d’incidence d’un avion.

Exemple résolu

Données
m = 1 kg, k = 16 N/m, c = 2 N·s/m.
Inconnue
ω_n, ζ, ω_d, valeurs propres et Q.
Stratégie
Appliquer les définitions, puis la formule des racines.
Étapes
  1. ω_n = √16 = 4,000 rad/s
  2. ζ = 2 / (2√16) = 0,250
  3. ω_d = 4√(1 − 0.0625) = 3,8730 rad/s
  4. s = -1,0 ± 3,8730 i
  5. Q = 1/(2ζ) = 2,0
Réponse
ω_n = 4 rad/s, ζ = 0,25, s = -1 ± 3,873i s⁻¹, Q = 2.
Vérification
|s|² = 1 + 15 = 16 = ω_n², comme il se doit ; avec ζ < 1 le mouvement oscille, et l’amplitude est divisée par e chaque seconde.
Portée
Lire la stabilité directement sur les valeurs propres est l’outil central de l’automatique et de la mécanique du vol (voir le mode d’incidence du parcours aéro).
Voir le code Python
import numpy as np
m, c, k = 1.0, 2.0, 16.0
A = np.array([[0, 1], [-k/m, -c/m]])
s = np.linalg.eigvals(A)
wn = np.sqrt(k/m); z = c/(2*np.sqrt(k*m))
print(s, wn, z, 1/(2*z))

Exercices corrigés automatiquement

    Solution complète

    ζ = 1,2 / (2·2) = 0,30 ; ω_n = 4 rad/s.

Exercices

  1. À quoi ressemble l’amortissement critique (ζ = 1), et où l’utilise-t-on ?

    Voir la solution

    Une racine réelle double : le retour au repos le plus rapide sans dépassement. Les ferme-portes et les aiguilles d’instruments sont conçus près de cette valeur.

Erreur fréquente

Confondre ω_n et ω_d, ou rad/s et hertz : f = ω / (2π).

Références

  • J. R. Taylor, Classical Mechanics, University Science Books, 2005.
  • D. Halliday, R. Resnick, J. Walker, Fundamentals of Physics, 12th ed., Wiley, 2021.

Mécanique lagrangienne et hamiltonienne

  • Master
  • Doctorat / recherche

Parcours : physique · Environ 45 min · Nombres vérifiés à la génération:

Vous saurez

Établir les équations du mouvement par le principe de Hamilton ; traiter le pendule exactement ; voir comment les symétries donnent des lois de conservation (Noether).

Prérequis

Bases du calcul des variations, dérivées partielles, le thème des oscillateurs. La physique expliquée

D’abord l’intuition

On écrit une seule fonction scalaire, L = énergie cinétique − énergie potentielle, dans les coordonnées naturelles (un angle, un rayon). Les équations du mouvement et les quantités conservées en sortent automatiquement, sans chercher les forces de liaison.

🧒 Explique-moi comme si j’avais 10 ans

Au lieu de suivre chaque poussée et chaque traction, on peut décrire une balançoire seulement par son énergie de mouvement et son énergie de hauteur. La nature choisit le chemin qui garde un équilibre spécial entre les deux, et tout le mouvement découle de cette seule règle.

Pas à pas

  1. S=∫L(q,q˙,t) dt,δS=0  ⇒  ddt∂L∂q˙−∂L∂q=0S=\int L(q,\dot q,t)\,dt,\quad \delta S=0\;\Rightarrow\;\frac{d}{dt}\frac{\partial L}{\partial\dot q}-\frac{\partial L}{\partial q}=0
    Le principe de Hamilton donne l’équation d’Euler–Lagrange.
  2. L=12mℓ2θ˙2+mgℓcos⁡θ  ⇒  θ¨=−gℓsin⁡θL=\tfrac12m\ell^2\dot\theta^2+mg\ell\cos\theta\;\Rightarrow\;\ddot\theta=-\frac{g}{\ell}\sin\theta
    Pendule, avec l’angle comme coordonnée.
  3. T=4ℓg K(k),k=sin⁡θ02  ⟺  T=2πℓ/gAGM⁡(1,cos⁡θ02)T=4\sqrt{\frac{\ell}{g}}\,K(k),\quad k=\sin\frac{\theta_0}{2}\;\Longleftrightarrow\;T=\frac{2\pi\sqrt{\ell/g}}{\operatorname{AGM}(1,\cos\frac{\theta_0}{2})}
    Période exacte par conservation de l’énergie ; la moyenne arithmético-géométrique calcule rapidement l’intégrale elliptique K.
  4. H=∑ipiq˙i−L,q˙=∂H∂p, p˙=−∂H∂q;∂L∂q=0⇒p=constH=\sum_i p_i\dot q_i-L,\quad \dot q=\frac{\partial H}{\partial p},\ \dot p=-\frac{\partial H}{\partial q};\qquad \frac{\partial L}{\partial q}=0\Rightarrow p=\text{const}
    Équations de Hamilton. Noether : une coordonnée cyclique implique un moment conservé (par ex. le moment cinétique dans un champ central).

Exemple résolu

Données
Pendule ℓ = 1 m, lâché sans vitesse à θ₀ = 60°, g = 9,80665 m/s².
Inconnue
Période aux petits angles et période exacte.
Stratégie
T₀ = 2π√(ℓ/g), puis diviser par AGM(1, cos 30°).
Étapes
  1. T₀ = 2π√(1/9.80665) = 2,0064 s
  2. cos 30° = 0.866025
  3. AGM(1, 0.866025) = 0.931808
  4. T = 2,0064 / 0.931808 = 2,1532 s
Réponse
T₀ = 2,006 s ; T exacte = 2,153 s, soit 7,3 % de plus.
Vérification
La série T ≈ T₀(1 + θ₀²/16) donne environ 6,9 % à 60° ; les termes suivants ajoutent le reste. Les grandes oscillations sont toujours plus lentes que la formule aux petits angles.
Portée
L’approche lagrangienne s’étend aux bras robotiques, à l’attitude des engins spatiaux et à la théorie des perturbations orbitales, où les schémas de forces deviennent impraticables.
Voir le code Python
import math
def agm(a, b):
    while abs(a-b) > 1e-15: a, b = (a+b)/2, math.sqrt(a*b)
    return a
L, g, th0 = 1.0, 9.80665, math.radians(60)
T0 = 2*math.pi*math.sqrt(L/g)
print(T0, T0/agm(1, math.cos(th0/2)))

Exercices corrigés automatiquement

  1. s

    Solution complète

    T₀ = 2,8375 s ; T = 3,3492 s (18 % de plus).

Exercices

  1. Dans un problème à force centrale en coordonnées polaires, quel moment est conservé et pourquoi ?

    Voir la solution

    L ne dépend pas de l’angle φ (symétrie de rotation), donc p_φ = mr²φ̇, le moment cinétique, est constant : c’est la deuxième loi de Kepler.

Erreur fréquente

Utiliser sin θ ≈ θ à grande amplitude. Au-delà d’environ 20°, la période aux petits angles est fausse de plus de 0,7 %.

Références

  • H. Goldstein, C. Poole, J. Safko, Classical Mechanics, 3rd ed., Addison-Wesley, 2002.
  • J. R. Taylor, Classical Mechanics, University Science Books, 2005.

Électromagnétisme, équations de Maxwell et plasma spatial

  • Licence (L1–L3)
  • Master

Parcours : physique · Environ 40 min · Nombres vérifiés à la génération:

Vous saurez

Énoncer les équations de Maxwell, en déduire la vitesse des ondes c, et calculer le rayon de giration et la fréquence plasma du vent solaire.

Prérequis

Analyse vectorielle (div, rot), deuxième loi de Newton. La physique expliquée

D’abord l’intuition

Quatre équations décrivent tous les effets électriques et magnétiques, des moteurs à la radio et à la lumière elle-même. Dans l’espace, où le gaz est ionisé, les mêmes lois font tourner les particules chargées autour des lignes de champ.

🧒 Explique-moi comme si j’avais 10 ans

Les aimants et les charges électriques poussent et tirent sans se toucher. Un aimant fait tourner en spirale les particules chargées venues du Soleil autour de ses lignes, comme des perles sur un fil : c’est ainsi que le champ magnétique terrestre les guide vers les pôles et crée les aurores.

Pas à pas

  1. ∇ ⁣⋅E=ρε0,∇ ⁣⋅B=0,∇ ⁣×E=−∂B∂t,∇ ⁣×B=μ0J+μ0ε0∂E∂t\nabla\!\cdot\mathbf E=\frac{\rho}{\varepsilon_0},\quad \nabla\!\cdot\mathbf B=0,\quad \nabla\!\times\mathbf E=-\frac{\partial\mathbf B}{\partial t},\quad \nabla\!\times\mathbf B=\mu_0\mathbf J+\mu_0\varepsilon_0\frac{\partial\mathbf E}{\partial t}
    Équations de Maxwell (SI).
  2. ∇2E=μ0ε0∂2E∂t2  ⇒  c=1μ0ε0\nabla^2\mathbf E=\mu_0\varepsilon_0\frac{\partial^2\mathbf E}{\partial t^2}\;\Rightarrow\;c=\frac{1}{\sqrt{\mu_0\varepsilon_0}}
    Dans le vide, on prend le rotationnel de la loi de Faraday : équation d’onde de vitesse c.
  3. mv˙=q(E+v×B)  ⇒  rg=mv⊥∣q∣B,ωg=∣q∣Bmm\dot{\mathbf v}=q(\mathbf E+\mathbf v\times\mathbf B)\;\Rightarrow\;r_g=\frac{mv_\perp}{|q|B},\quad \omega_g=\frac{|q|B}{m}
    Force de Lorentz : giration circulaire perpendiculaire à B.
  4. fpe=12πne2ε0mef_{pe}=\frac{1}{2\pi}\sqrt{\frac{n e^2}{\varepsilon_0 m_e}}
    Fréquence plasma électronique : les ondes radio en dessous ne traversent pas le plasma.

Exemple résolu

Données
Proton du vent solaire près de la Terre : v⊥ = 400 km/s, champ interplanétaire B = 5 nT ; densité électronique n = 5 cm⁻³ (valeurs typiques).
Inconnue
Rayon de giration, période de giration et fréquence plasma électronique.
Stratégie
Mettre les constantes CODATA dans r_g, ω_g et f_pe ; tout en SI (n = 5×10⁶ m⁻³).
Étapes
  1. r_g = (1.6726e-27 · 4e5) / (1.6022e-19 · 5e-9) = 835 km
  2. ω_g = eB/m = 0,4789 rad/s
  3. T_g = 2π/ω_g = 13,12 s
  4. f_pe = 20,08 kHz
  5. c = 1/√(μ₀ε₀) = 299792458 m/s
Réponse
r_g ≈ 835 km, période de giration 13,1 s, f_pe ≈ 20,1 kHz.
Vérification
Le rayon de giration est bien plus petit que la magnétosphère terrestre (des dizaines de milliers de km) : le plasma se comporte comme un fluide lié au champ. Le c obtenu à partir de μ₀ et ε₀ redonne la valeur définie 299 792 458 m/s.
Portée
Ces échelles déterminent comment la météo spatiale atteint satellites et réseaux électriques, et pourquoi certaines fréquences radio sont réfléchies par l’ionosphère.
Voir le code Python
import math
mp, e, me, eps0 = 1.67262192369e-27, 1.602176634e-19, 9.1093837015e-31, 8.8541878128e-12
v, B, n = 4.0e5, 5e-9, 5e6
rg = mp*v/(e*B); wg = e*B/mp
fpe = math.sqrt(n*e*e/(eps0*me))/(2*math.pi)
print(rg/1e3, 2*math.pi/wg, fpe/1e3)

Exercices corrigés automatiquement

  1. m

    Solution complète

    r_g = 9,109e-31 · 1e6 / (1,602e-19 · 3e-5) = 0,190 m.

Exercices

  1. Pourquoi ∇·B = 0 est-il important physiquement ?

    Voir la solution

    Il n’y a pas de monopôles magnétiques : les lignes de champ ne commencent ni ne finissent, elles se referment. Les codes numériques de plasma doivent le respecter, sinon ils créent de fausses forces.

Erreur fréquente

Oublier les préfixes : 5 nT = 5×10⁻⁹ T, et 5 cm⁻³ = 5×10⁶ m⁻³.

Références

  • D. J. Griffiths, Introduction to Electrodynamics, 4th ed., Cambridge University Press, 2017.
  • F. F. Chen, Introduction to Plasma Physics and Controlled Fusion, 3rd ed., Springer, 2016.
  • E. Tiesinga et al., “CODATA Recommended Values of the Fundamental Physical Constants: 2018”, Rev. Mod. Phys. 93, 025010 (2021).

Thermodynamique et physique statistique

  • Premier cycle (DUT/BTS)
  • Licence (L1–L3)
  • Master

Parcours : physique · Environ 35 min · Nombres vérifiés à la génération:

Vous saurez

Relier la température à la vitesse moléculaire par la distribution de Maxwell–Boltzmann ; appliquer les premier et deuxième principes ; calculer le rendement de Carnot.

Prérequis

Loi des gaz parfaits, intégrales gaussiennes. La physique expliquée

D’abord l’intuition

La température est l’énergie cinétique moyenne par molécule. Compter les façons dont les molécules peuvent se partager l’énergie (l’entropie) explique pourquoi un moteur thermique ne peut jamais transformer toute sa chaleur en travail.

🧒 Explique-moi comme si j’avais 10 ans

L’air est fait de minuscules molécules qui foncent plus vite qu’un avion de ligne. « Chaud » veut juste dire qu’elles vont plus vite. Et la chaleur va toujours du chaud vers le froid, jamais l’inverse toute seule.

Pas à pas

  1. f(v)=4π(m2πkBT)3/2v2e−mv2/(2kBT)f(v)=4\pi\left(\frac{m}{2\pi k_BT}\right)^{3/2}v^2e^{-mv^2/(2k_BT)}
    Distribution des vitesses de Maxwell–Boltzmann.
  2. 12m⟨v2⟩=32kBT  ⇒  vrms=3kBTm,vˉ=8kBTπm\tfrac12m\langle v^2\rangle=\tfrac32k_BT\;\Rightarrow\;v_{rms}=\sqrt{\frac{3k_BT}{m}},\qquad \bar v=\sqrt{\frac{8k_BT}{\pi m}}
    Moments de la distribution (équipartition).
  3. dU=δQ−δW,dS≥δQTdU=\delta Q-\delta W,\qquad dS\ge\frac{\delta Q}{T}
    Premier principe (énergie) et deuxième principe (entropie).
  4. η≤ηC=1−TcTh\eta\le\eta_C=1-\frac{T_c}{T_h}
    Aucun moteur entre deux sources ne dépasse Carnot.

Exemple résolu

Données
Azote N₂ (masse molaire 28,0134 g/mol) à 300 K (27 °C).
Inconnue
Vitesses quadratique moyenne et moyenne.
Stratégie
Utiliser R = k_B N_A et M en kg/mol.
Étapes
  1. R = k_B N_A = 8.314463 J/(mol·K)
  2. v_rms = √(3 · 8.314463 · 300 / 0.0280134) = 516,8 m/s
  3. v̄ = √(8 · 8.314463 · 300 / (π · 0.0280134)) = 476,2 m/s
Réponse
v_rms ≈ 517 m/s, v̄ ≈ 476 m/s.
Vérification
Le rapport v_rms/v̄ = √(3π/8) ≈ 1,085, comme pour tout gaz de Maxwell–Boltzmann ; les deux dépassent un peu la vitesse du son dans l’air (environ 347 m/s à 300 K), comme attendu.
Portée
Les vitesses moléculaires fixent les vitesses de réaction, la rapidité des fuites et de la diffusion des gaz, et pourquoi les gaz légers s’échappent des atmosphères planétaires.
Voir le code Python
import math
kB, NA = 1.380649e-23, 6.02214076e23
R, M, T = kB*NA, 0.0280134, 300.0
print(math.sqrt(3*R*T/M), math.sqrt(8*R*T/(math.pi*M)))

Exercices corrigés automatiquement

  1. m/s

    Solution complète

    v_rms = √(3RT/M) = 483,6 m/s (contre 516,8 m/s pour N₂).

Exercices

  1. Un moteur fonctionne entre 600 K et 300 K. Quel est le meilleur rendement possible ?

    Voir la solution

    η_C = 1 − 300/600 = 0,50, soit 50 %. Les moteurs réels font nettement moins à cause des frottements et des transferts de chaleur à vitesse finie.

Erreur fréquente

Utiliser des grammes par mole au lieu de kilogrammes par mole, ce qui rend les vitesses √1000 ≈ 32 fois trop petites.

Références

  • D. V. Schroeder, An Introduction to Thermal Physics, Addison-Wesley, 2000.
  • D. Halliday, R. Resnick, J. Walker, Fundamentals of Physics, 12th ed., Wiley, 2021.
  • E. Tiesinga et al., “CODATA Recommended Values of the Fundamental Physical Constants: 2018”, Rev. Mod. Phys. 93, 025010 (2021).

Fluides : Navier–Stokes et nombre de Reynolds

  • Licence (L1–L3)
  • Master
  • Doctorat / recherche

Parcours : physique · Environ 40 min · Nombres vérifiés à la génération:

Vous saurez

Écrire les équations de Navier–Stokes incompressibles, les adimensionner pour trouver le nombre de Reynolds, et résoudre exactement l’écoulement en conduite (Poiseuille).

Prérequis

Analyse vectorielle, deuxième loi de Newton, le thème atmosphère du parcours aéro. La physique expliquée

D’abord l’intuition

Navier–Stokes, c’est la deuxième loi de Newton pour une petite parcelle de fluide : pression, viscosité et inertie s’affrontent. Que l’écoulement soit régulier ou turbulent dépend du rapport inertie/viscosité.

🧒 Explique-moi comme si j’avais 10 ans

Le miel coule tout doucement, mais l’eau d’un robinet ouvert à fond éclabousse et tourbillonne. Un seul nombre, le nombre de Reynolds, dit ce qui va se passer : grand, ça tourbillonne ; petit, ça coule en douceur.

Pas à pas

  1. ρ(∂u∂t+u⋅∇u)=−∇p+μ∇2u,∇ ⁣⋅u=0\rho\left(\frac{\partial\mathbf u}{\partial t}+\mathbf u\cdot\nabla\mathbf u\right)=-\nabla p+\mu\nabla^2\mathbf u,\qquad \nabla\!\cdot\mathbf u=0
    Navier–Stokes incompressible et continuité.
  2. u∗=uU, x∗=xL  ⇒  …=−∇∗p∗+1Re∇∗2u∗,Re=ρULμ\mathbf u^*=\frac{\mathbf u}{U},\ \mathbf x^*=\frac{\mathbf x}{L}\;\Rightarrow\;\ldots=-\nabla^*p^*+\frac{1}{Re}\nabla^{*2}\mathbf u^*,\quad Re=\frac{\rho UL}{\mu}
    L’adimensionnement ne laisse qu’un paramètre : le nombre de Reynolds.
  3. u(r)=Δp4μL(R2−r2)  ⇒  Q=πR4Δp8μLu(r)=\frac{\Delta p}{4\mu L}(R^2-r^2)\;\Rightarrow\;Q=\frac{\pi R^4\Delta p}{8\mu L}
    Écoulement permanent en conduite : solution exacte (Hagen–Poiseuille), valable tant que Re reste sous environ 2300.
  4. μ(T)=1.458×10−6 T3/2T+110.4  Pa⋅s\mu(T)=\frac{1.458\times10^{-6}\,T^{3/2}}{T+110.4}\ \ \text{Pa·s}
    Loi de Sutherland pour la viscosité de l’air (celle de l’atmosphère standard de 1976).

Exemple résolu

Données
Aile d’avion de ligne générique en croisière : corde L = 5 m, U = 250 m/s, à 11 km : ρ = 0,3639 kg/m³, T = 216,65 K.
Inconnue
Viscosité μ et nombre de Reynolds Re.
Stratégie
Sutherland pour μ, puis Re = ρUL/μ.
Étapes
  1. μ = 1.458e-6 · 216.65^1.5 / (216.65 + 110.4) = 0,00001422 Pa·s
  2. Re = 0.3639 · 250 · 5 / μ
  3. Re ≈ 3.2 × 10⁷
Réponse
μ ≈ 1,42×10⁻⁵ Pa·s, Re ≈ 3,2×10⁷ : la couche limite est surtout turbulente.
Vérification
L’atmosphère standard de 1976 donne μ = 1,4216×10⁻⁵ Pa·s à 11 km. Un Re de plusieurs dizaines de millions est typique des avions réels, bien au-dessus des maquettes en soufflerie : d’où l’importance des effets d’échelle.
Portée
Re détermine la traînée, la transition et le décollement ; le reproduire est le problème central des essais en soufflerie et de la modélisation de la turbulence en CFD.
Voir le code Python
rho, U, L, T = 0.363918, 250.0, 5.0, 216.65
mu = 1.458e-6*T**1.5/(T + 110.4)
print(mu, rho*U*L/mu)

Exercices corrigés automatiquement

    Solution complète

    Re = 1,225 · 30 · 0,2 / 1,789×10⁻⁵ ≈ 410757 : environ 80 fois plus petit que l’aile en croisière.

Exercices

  1. Doubler le rayon d’une conduite à perte de charge égale multiplie le débit laminaire par combien ?

    Voir la solution

    Q ∝ R⁴, donc par 16. Exemple : R = 0,5 mm, L = 10 cm, Δp = 1 kPa, eau μ = 1,0 mPa·s donne Q = 0,245 mL/s.

Erreur fréquente

Utiliser la viscosité cinématique ν (m²/s) là où il faut la viscosité dynamique μ (Pa·s), ou l’inverse : Re = ρUL/μ = UL/ν.

Références

  • F. M. White, Fluid Mechanics, 8th ed., McGraw-Hill, 2016.
  • J. D. Anderson Jr., Fundamentals of Aerodynamics, 6th ed., McGraw-Hill, 2017.

La relativité pour les orbites : horloges GPS et Mercure

  • Master
  • Doctorat / recherche

Parcours : physique · Environ 45 min · Nombres vérifiés à la génération:

Vous saurez

Estimer les dérives d’horloge relativistes restreinte et générale d’un satellite GPS, et l’avance relativiste du périhélie de Mercure.

Prérequis

Mécanique orbitale (vitesse circulaire, vis-viva), développements en série. La physique expliquée

D’abord l’intuition

Les horloges en mouvement retardent (relativité restreinte) et celles plus bas dans la gravité retardent aussi (relativité générale). Pour le GPS, l’effet de gravité l’emporte, et l’ensemble est assez grand pour que le système échoue en quelques minutes sans correction.

🧒 Explique-moi comme si j’avais 10 ans

Les horloges des satellites GPS avancent un tout petit peu plus vite que celles au sol, parce que la gravité est plus faible là-haut. Sans correction, ta carte se décalerait d’environ 12 km par jour.

Pas à pas

  1. dτdt≈1+Φc2−v22c2,Φ=−GMr\frac{d\tau}{dt}\approx1+\frac{\Phi}{c^2}-\frac{v^2}{2c^2},\qquad \Phi=-\frac{GM}{r}
    Rythme d’horloge en champ faible : terme de potentiel (RG) plus terme de vitesse (RR).
  2. Δff∣GR=GMc2(1R⊕−1a),Δff∣SR=−GM2ac2\left.\frac{\Delta f}{f}\right|_{GR}=\frac{GM}{c^2}\left(\frac1{R_\oplus}-\frac1a\right),\qquad \left.\frac{\Delta f}{f}\right|_{SR}=-\frac{GM}{2ac^2}
    Satellite contre sol, orbite circulaire v² = GM/a (mouvement de l’horloge au sol négligé ici).
  3. Δϖ=6πGM⊙c2a(1−e2) per orbit\Delta\varpi=\frac{6\pi GM_\odot}{c^2a(1-e^2)}\ \text{per orbit}
    Avance du périhélie en RG, à partir de l’équation de l’orbite de Schwarzschild.

Exemple résolu

Données
Orbite GPS a = 26 561,75 km (circulaire), R⊕ = 6378,137 km, GM⊕ = 3,986004418×10¹⁴ m³/s².
Inconnue
Gain d’horloge par jour dû à la RG et à la RR, le total, et l’erreur de distance équivalente.
Stratégie
Calculer les deux décalages relatifs, multiplier par 86 400 s, puis le temps net par c.
Étapes
  1. GR: 45,65 μs/day
  2. SR: -7,21 μs/day
  3. net: 38,44 μs/day
  4. c · 38,44 μs = 11,5 km/day
Réponse
Environ +45,7 μs/jour (gravité) et -7,2 μs/jour (vitesse) : total ≈ +38,4 μs/jour, soit environ 12 km d’erreur de distance par jour.
Vérification
Ashby (2003), qui inclut la rotation et l’aplatissement de la Terre, donne environ +45,7 et −7,1 μs/jour, soit 38,6 μs/jour : notre modèle simple concorde à 0,2 μs/jour près.
Portée
Les horloges des satellites GPS sont volontairement réglées pour retarder avant le lancement, afin de rester synchrones avec le sol une fois en orbite. La relativité est une correction d’ingénierie quotidienne.
Voir le code Python
GM, c, R, a = 3.986004418e14, 299792458.0, 6378137.0, 26561750.0
gr = GM*(1/R - 1/a)/c**2*86400e6
sr = -GM/(2*a*c**2)*86400e6
print(gr, sr, gr + sr, (gr + sr)*1e-6*c/1e3)

Exercices corrigés automatiquement

  1. ″/century

    Solution complète

    Δϖ ≈ 8,62″ par siècle pour Vénus, contre 42,98″ pour Mercure (Will 2014 donne 42,98″ pour Mercure).

Exercices

  1. Pourquoi la RG modifie-t-elle plus l’orbite de Mercure que celle de la Terre ?

    Voir la solution

    L’avance par orbite varie comme 1/(a(1−e²)), et Mercure fait environ 4 fois plus d’orbites par siècle : l’effet par siècle est donc bien plus grand près du Soleil.

Erreur fréquente

Ajouter les termes RR et RG avec le même signe. La vitesse du satellite ralentit son horloge (négatif), son altitude l’accélère (positif).

Références

  • N. Ashby, “Relativity in the Global Positioning System”, Living Reviews in Relativity 6, 1 (2003).
  • C. M. Will, “The Confrontation between General Relativity and Experiment”, Living Reviews in Relativity 17, 4 (2014).
  • E. Tiesinga et al., “CODATA Recommended Values of the Fundamental Physical Constants: 2018”, Rev. Mod. Phys. 93, 025010 (2021).

Vérifiez votre compréhension

Un oscillateur amorti a pour valeurs propres −0,5 ± 3i. Il est…

Dans quel sens les horloges des satellites GPS dérivent-elles, au total, par rapport au sol ?

Doubler la température absolue d’un gaz multiplie la vitesse quadratique moyenne par…

Pour aller plus loin

Ensuite : Taylor et Goldstein pour la mécanique, Griffiths et Chen pour les champs et le plasma, Schroeder pour la physique thermique, et la revue en accès libre d’Ashby pour la relativité dans le GPS.

Pour les enseignants

Usage conseillé : une séance de 50 minutes par section. Parcours premier cycle/licence : Newton et conservation → projectiles avec traînée (les étudiants modifient le Python RK4 et comparent avec le Physics Lab) → oscillateurs → thermodynamique. Parcours master : mécanique lagrangienne → électromagnétisme et plasma spatial → Navier–Stokes et nombre de Reynolds → relativité pour le GPS. Chaque réponse numérique est calculée par le script de génération à partir des constantes CODATA 2018 et sert de corrigé. Aligné sur les thèmes MIT OCW 8.01/8.02.

Références

Aider / Faire un don