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
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.
Lois de Newton et conservation
- Premier cycle (DUT/BTS)
- Licence (L1–L3)
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
- Deuxième loi sous forme de quantité de mouvement.
- Troisième loi : les forces internes s’annulent, la quantité de mouvement totale est conservée.
- Théorème de l’énergie cinétique.
- 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
- v₁ = (2 − 1)·3 / 3 = 1,00 m/s
- v₂ = 2·2·3 / 3 = 4,00 m/s
- p: 2·3 = 6 → 2·1 + 1·4 = 6 kg·m/s
- 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
- m/s
- v₂ = 2m₁u₁ / (m₁ + m₂).
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
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
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
- Traînée quadratique, opposée à la vitesse.
- Couplées par |v| : pas de solution élémentaire.
- Vitesse limite : la traînée équilibre le poids.
- 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
- k = ρC_DA/(2m) = 0,00636 1/m
- vacuum: R = v₀² sin 70° / g = 153,3 m, h = 26,8 m, t = 4,68 s
- RK4 (Δt = 1 ms): R = 93,9 m, h = 20,2 m, t = 4,05 s
- RK4 (Δt = 0.5 ms): R = 93,9 m
- 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
- m
- Changez seulement l’angle de lancement dans le code.
- Attendez-vous à environ 1 m de plus qu’à 35° : avec traînée, le meilleur angle est d’environ 41,0°.
Solution complète
R(45°) = 94,8 m, contre 93,9 m à 35° et 163,2 m dans le vide à 45°.
Exercices
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
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
- Loi de Newton avec force de rappel et amortissement visqueux.
- Équation caractéristique.
- Sous-amorti : paire complexe ; partie réelle = taux de décroissance, partie imaginaire = pulsation amortie.
- 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
- ω_n = √16 = 4,000 rad/s
- ζ = 2 / (2√16) = 0,250
- ω_d = 4√(1 − 0.0625) = 3,8730 rad/s
- s = -1,0 ± 3,8730 i
- 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
-
- ζ = c / (2√(km)).
- √(8 · 0,5) = 2.
Solution complète
ζ = 1,2 / (2·2) = 0,30 ; ω_n = 4 rad/s.
Exercices
À 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
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
- Le principe de Hamilton donne l’équation d’Euler–Lagrange.
- Pendule, avec l’angle comme coordonnée.
- Période exacte par conservation de l’énergie ; la moyenne arithmético-géométrique calcule rapidement l’intégrale elliptique K.
- É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
- T₀ = 2π√(1/9.80665) = 2,0064 s
- cos 30° = 0.866025
- AGM(1, 0.866025) = 0.931808
- 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
- s
- T₀ = 2π√(2/g).
- Diviser par AGM(1, cos 45°) ≈ 0,8472.
Solution complète
T₀ = 2,8375 s ; T = 3,3492 s (18 % de plus).
Exercices
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
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
- Équations de Maxwell (SI).
- Dans le vide, on prend le rotationnel de la loi de Faraday : équation d’onde de vitesse c.
- Force de Lorentz : giration circulaire perpendiculaire à B.
- 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
- r_g = (1.6726e-27 · 4e5) / (1.6022e-19 · 5e-9) = 835 km
- ω_g = eB/m = 0,4789 rad/s
- T_g = 2π/ω_g = 13,12 s
- f_pe = 20,08 kHz
- 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
- m
- r_g = m_e v / (eB) avec m_e = 9,109×10⁻³¹ kg.
Solution complète
r_g = 9,109e-31 · 1e6 / (1,602e-19 · 3e-5) = 0,190 m.
Exercices
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
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
- Distribution des vitesses de Maxwell–Boltzmann.
- Moments de la distribution (équipartition).
- Premier principe (énergie) et deuxième principe (entropie).
- 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
- R = k_B N_A = 8.314463 J/(mol·K)
- v_rms = √(3 · 8.314463 · 300 / 0.0280134) = 516,8 m/s
- 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
- m/s
- M = 0,0319988 kg/mol.
- Les molécules plus lourdes sont plus lentes : v ∝ 1/√M.
Solution complète
v_rms = √(3RT/M) = 483,6 m/s (contre 516,8 m/s pour N₂).
Exercices
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
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
- Navier–Stokes incompressible et continuité.
- L’adimensionnement ne laisse qu’un paramètre : le nombre de Reynolds.
- Écoulement permanent en conduite : solution exacte (Hagen–Poiseuille), valable tant que Re reste sous environ 2300.
- 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.458e-6 · 216.65^1.5 / (216.65 + 110.4) = 0,00001422 Pa·s
- Re = 0.3639 · 250 · 5 / μ
- 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
-
- μ(288,15 K) par Sutherland ≈ 1,79×10⁻⁵ Pa·s.
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
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
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
- Rythme d’horloge en champ faible : terme de potentiel (RG) plus terme de vitesse (RR).
- Satellite contre sol, orbite circulaire v² = GM/a (mouvement de l’horloge au sol négligé ici).
- 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
- GR: 45,65 μs/day
- SR: -7,21 μs/day
- net: 38,44 μs/day
- 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
- ″/century
- Par orbite : 6πGM☉ / (c²a(1 − e²)), avec GM☉ = 1,32712440018×10²⁰ m³/s².
- Orbites par siècle = 36 525 / 224,701 ; convertir les radians en secondes d’arc.
Solution complète
Δϖ ≈ 8,62″ par siècle pour Vénus, contre 42,98″ pour Mercure (Will 2014 donne 42,98″ pour Mercure).
Exercices
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
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
- D. Halliday, R. Resnick, J. Walker, Fundamentals of Physics, 12th ed., Wiley, 2021.
- J. R. Taylor, Classical Mechanics, University Science Books, 2005.
- H. Goldstein, C. Poole, J. Safko, Classical Mechanics, 3rd ed., Addison-Wesley, 2002.
- 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.
- D. V. Schroeder, An Introduction to Thermal Physics, Addison-Wesley, 2000.
- F. M. White, Fluid Mechanics, 8th ed., McGraw-Hill, 2016.
- 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).
- W. H. Press et al., Numerical Recipes, 3rd ed., Cambridge University Press, 2007.
- MIT OpenCourseWare, 8.01 Classical Mechanics and 8.02 Electricity and Magnetism.
- J. D. Anderson Jr., Fundamentals of Aerodynamics, 6th ed., McGraw-Hill, 2017.