La fisica passo dopo passo
Meccanica, resistenza dell’aria, oscillazioni, lagrangiane, elettromagnetismo e plasma, termodinamica, fluidi e la relatività che serve al GPS: dimostrazioni, esempi svolti ed esercizi, dal post-diploma al dottorato.
Costanti usate in questa pagina
Costanti: CODATA 2018 (c, e, m_p, m_e, ε₀, μ₀, k_B e N_A come pubblicate). Ogni numero di questa pagina è calcolato dallo script di build del sito.
Leggi di Newton e conservazione
- Corso post-diploma
- Laurea triennale
Saprai
Enunciare le leggi di Newton; ricavare la conservazione della quantità di moto e dell’energia; risolvere un urto elastico in 1 D.
Prerequisiti
Algebra, vettori, derivate. La fisica spiegata
Prima l’intuizione
Le forze cambiano il moto. Se nessuna forza esterna agisce su un gruppo di oggetti, la loro quantità di moto totale non può cambiare; e se non si perde energia in calore o suono, nemmeno l’energia cinetica totale.
🧒 Spiegamelo come se avessi 10 anni
Spingi un carrello vuoto e accelera; spingine uno pieno e accelera meno. E quando due autoscontri si urtano, la spinta che uno perde l’altro la guadagna: il totale resta uguale.
Passo dopo passo
- Seconda legge in forma di quantità di moto.
- Terza legge: le forze interne si annullano, la quantità di moto totale si conserva.
- Teorema dell’energia cinetica.
- Si risolvono insieme conservazione della quantità di moto e dell’energia cinetica per un urto elastico.
Esempio svolto
- Dati
- Carrello 1: m₁ = 2 kg a u₁ = 3 m/s urta il carrello 2: m₂ = 1 kg fermo. Urto elastico su binario orizzontale senza attrito.
- Incognita
- Velocità finali v₁ e v₂.
- Strategia
- Applicare le formule dell’urto elastico, poi verificare quantità di moto ed energia.
- Passaggi
- 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
- Risposta
- v₁ = 1,0 m/s, v₂ = 4,0 m/s (entrambi in avanti).
- Controllo di plausibilità
- Quantità di moto (6 kg·m/s) ed energia cinetica (9 J) invariate, e la velocità relativa si inverte (3 m/s prima, 3 m/s dopo), come deve in un urto elastico.
- Significato
- La stessa algebra descrive una fionda gravitazionale di una sonda, vista dal sistema del pianeta.
Mostra il codice 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)Esercizi con correzione automatica
- m/s
- v₂ = 2m₁u₁ / (m₁ + m₂).
Soluzione completa
v₂ = 2·1·4/4 = 2,0 m/s; il carrello leggero rimbalza a -2,0 m/s.
Esercizi
Una palla cade da 20 m da ferma (senza aria). Con che velocità arriva in fondo?
Mostra la soluzione
mgh = ½mv², quindi v = √(2gh) = 19,81 m/s, indipendentemente dalla massa.
Errore frequente
Trattare la quantità di moto come un numero e non come un vettore: un carrello che arretra ha quantità di moto negativa.
Riferimenti
- 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.
Proiettili con resistenza dell’aria (EDO numerica)
- Laurea triennale
- Laurea magistrale
Saprai
Scrivere le equazioni del moto con resistenza quadratica, capire perché non hanno soluzione in forma chiusa e integrarle con Runge–Kutta del quarto ordine.
Prerequisiti
Seconda legge di Newton, vettori, Python di base. La fisica spiegata
Prima l’intuizione
La resistenza cresce col quadrato della velocità e si oppone sempre ad essa. Questo accoppia moto orizzontale e verticale, quindi facciamo avanzare la soluzione a piccoli passi di tempo al computer.
🧒 Spiegamelo come se avessi 10 anni
Lancia una palla e l’aria continua a frenarla. Non va lontano né in alto come nello spazio, e scende più ripida di come è salita.
Passo dopo passo
- Resistenza quadratica, opposta alla velocità.
- Accoppiate tramite |v|: nessuna soluzione elementare.
- Velocità limite: la resistenza bilancia il peso.
- RK4: l’errore globale cala come Δt⁴.
Esempio svolto
- Dati
- Palla generica: m = 0,145 kg, diametro 0,074 m, C_D = 0,35, ρ = 1,225 kg/m³; lanciata a 40 m/s (131 ft/s), a 35° sopra un terreno piano.
- Incognita
- Gittata, altezza massima e tempo di volo, con e senza aria.
- Strategia
- Formula chiusa nel vuoto; RK4 con Δt = 1 ms con resistenza, poi Δt dimezzato per verificare la convergenza.
- Passaggi
- 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
- Risposta
- Con resistenza: gittata 93,9 m (contro 153,3 m), apice 20,2 m (contro 26,8 m), tempo 4,05 s (contro 4,68 s).
- Controllo di plausibilità
- Dimezzare il passo cambia la gittata di meno di 1 mm: l’integrazione è convergente. La velocità di lancio è circa quella limite (39,3 m/s), quindi al lancio la resistenza è grande quanto il peso: una forte riduzione è attesa.
- Significato
- Lo stesso integratore, con spinta e densità variabile, serve a calcolare le traiettorie di razzi e capsule in rientro.
Mostra il codice 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])Esercizi con correzione automatica
- m
- Cambia solo l’angolo di lancio nel codice.
- Aspettati solo circa 1 m in più che a 35°: con resistenza l’angolo migliore è circa 41,0°.
Soluzione completa
R(45°) = 94,8 m, contro 93,9 m a 35° e 163,2 m nel vuoto a 45°.
Esercizi
Perché la palla scende più ripida di come è salita?
Mostra la soluzione
La resistenza toglie velocità orizzontale per tutto il tempo, mentre la gravità aggiunge velocità verticale in discesa: l’angolo di discesa è più ripido.
Errore frequente
Applicare la resistenza separatamente come −k v_x² e −k v_y². Il modulo dipende dalla velocità totale |v|, non da ogni componente da sola.
Riferimenti
- 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.
Oscillatori: semplice, smorzato e autovalori
- Laurea triennale
- Laurea magistrale
Saprai
Risolvere l’oscillatore armonico smorzato con l’equazione caratteristica; legare gli autovalori a pulsazione naturale, smorzamento e fattore di qualità; collegare alla stabilità degli aerei.
Prerequisiti
EDO lineari a coefficienti costanti, numeri complessi. La fisica spiegata
Prima l’intuizione
Quasi ogni sistema vicino all’equilibrio si comporta come una massa su una molla con un po’ di attrito: ponti, circuiti, molecole e aerei. Due numeri lo riassumono: pulsazione naturale e smorzamento.
🧒 Spiegamelo come se avessi 10 anni
Un’altalena va avanti e indietro da sola. Se smetti di spingere, l’aria e la catena che cigola le tolgono pian piano energia e si ferma. Gli ingegneri misurano quanto veloce oscilla e quanto in fretta si calma.
Passo dopo passo
- Legge di Newton con forza elastica e smorzamento viscoso.
- Equazione caratteristica.
- Sottosmorzato: coppia complessa; parte reale = tasso di decadimento, parte immaginaria = pulsazione smorzata.
- La stessa cosa come sistema del primo ordine: gli autovalori della matrice sono le radici s. È proprio l’analisi di corto periodo di un aereo.
Esempio svolto
- Dati
- m = 1 kg, k = 16 N/m, c = 2 N·s/m.
- Incognita
- ω_n, ζ, ω_d, autovalori e Q.
- Strategia
- Sostituire nelle definizioni, poi nella formula delle radici.
- Passaggi
- ω_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
- Risposta
- ω_n = 4 rad/s, ζ = 0,25, s = -1 ± 3,873i s⁻¹, Q = 2.
- Controllo di plausibilità
- |s|² = 1 + 15 = 16 = ω_n², come deve; con ζ < 1 il moto oscilla e l’ampiezza cala di un fattore e ogni secondo.
- Significato
- Leggere la stabilità direttamente dagli autovalori è lo strumento centrale di controlli e dinamica del volo (vedi il corto periodo nel percorso aero).
Mostra il codice 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))Esercizi con correzione automatica
-
- ζ = c / (2√(km)).
- √(8 · 0,5) = 2.
Soluzione completa
ζ = 1,2 / (2·2) = 0,30; ω_n = 4 rad/s.
Esercizi
Com’è lo smorzamento critico (ζ = 1) e dove si usa?
Mostra la soluzione
Una radice reale doppia: il ritorno a riposo più rapido senza sovraelongazione. Chiudiporta e lancette degli strumenti sono progettati vicino a questo valore.
Errore frequente
Confondere ω_n con ω_d, o rad/s con hertz: f = ω / (2π).
Riferimenti
- J. R. Taylor, Classical Mechanics, University Science Books, 2005.
- D. Halliday, R. Resnick, J. Walker, Fundamentals of Physics, 12th ed., Wiley, 2021.
Meccanica lagrangiana e hamiltoniana
- Laurea magistrale
- Dottorato / ricerca
Saprai
Ricavare le equazioni del moto dal principio di Hamilton; trattare il pendolo in modo esatto; vedere come le simmetrie danno leggi di conservazione (Noether).
Prerequisiti
Basi del calcolo delle variazioni, derivate parziali, il tema degli oscillatori. La fisica spiegata
Prima l’intuizione
Si scrive una sola funzione scalare, L = energia cinetica − energia potenziale, nelle coordinate naturali (un angolo, un raggio). Equazioni del moto e grandezze conservate escono da sole, senza cercare le forze vincolari.
🧒 Spiegamelo come se avessi 10 anni
Invece di seguire ogni spinta e trazione, puoi descrivere un’altalena solo con la sua energia di moto e la sua energia di altezza. La natura sceglie il percorso che tiene un equilibrio speciale tra le due, e da quell’unica regola segue tutto il moto.
Passo dopo passo
- Il principio di Hamilton dà l’equazione di Eulero–Lagrange.
- Pendolo, con l’angolo come coordinata.
- Periodo esatto dalla conservazione dell’energia; la media aritmetico-geometrica calcola in fretta l’integrale ellittico K.
- Equazioni di Hamilton. Noether: una coordinata ciclica implica un momento conservato (ad es. il momento angolare in una forza centrale).
Esempio svolto
- Dati
- Pendolo ℓ = 1 m, lasciato da fermo a θ₀ = 60°, g = 9,80665 m/s².
- Incognita
- Periodo per piccoli angoli e periodo esatto.
- Strategia
- T₀ = 2π√(ℓ/g), poi dividere per AGM(1, cos 30°).
- Passaggi
- 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
- Risposta
- T₀ = 2,006 s; T esatto = 2,153 s, il 7,3 % in più.
- Controllo di plausibilità
- La serie T ≈ T₀(1 + θ₀²/16) dà circa il 6,9 % a 60°; i termini successivi aggiungono il resto. Le grandi oscillazioni sono sempre più lente della formula per piccoli angoli.
- Significato
- L’approccio lagrangiano si estende a bracci robotici, assetto dei veicoli spaziali e teoria delle perturbazioni orbitali, dove i diagrammi delle forze diventano impraticabili.
Mostra il codice 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)))Esercizi con correzione automatica
- s
- T₀ = 2π√(2/g).
- Dividi per AGM(1, cos 45°) ≈ 0,8472.
Soluzione completa
T₀ = 2,8375 s; T = 3,3492 s (il 18 % in più).
Esercizi
In un problema a forza centrale in coordinate polari, quale momento si conserva e perché?
Mostra la soluzione
L non dipende dall’angolo φ (simmetria di rotazione), quindi p_φ = mr²φ̇, il momento angolare, è costante: è la seconda legge di Keplero.
Errore frequente
Usare sin θ ≈ θ a grande ampiezza. Oltre circa 20° il periodo per piccoli angoli sbaglia di più dello 0,7 %.
Riferimenti
- H. Goldstein, C. Poole, J. Safko, Classical Mechanics, 3rd ed., Addison-Wesley, 2002.
- J. R. Taylor, Classical Mechanics, University Science Books, 2005.
Elettromagnetismo, equazioni di Maxwell e plasma spaziale
- Laurea triennale
- Laurea magistrale
Saprai
Enunciare le equazioni di Maxwell, ricavare la velocità delle onde c e calcolare raggio di girazione e frequenza di plasma del vento solare.
Prerequisiti
Analisi vettoriale (div, rot), seconda legge di Newton. La fisica spiegata
Prima l’intuizione
Quattro equazioni descrivono ogni effetto elettrico e magnetico, dai motori alla radio fino alla luce stessa. Nello spazio, dove il gas è ionizzato, le stesse leggi fanno girare le particelle cariche attorno alle linee di campo.
🧒 Spiegamelo come se avessi 10 anni
Magneti e cariche elettriche spingono e tirano senza toccarsi. Un magnete fa girare a spirale le particelle cariche del Sole attorno alle sue linee, come perline su un filo: così il campo magnetico terrestre le guida verso i poli e crea le aurore.
Passo dopo passo
- Equazioni di Maxwell (SI).
- Nel vuoto si prende il rotore della legge di Faraday: equazione d’onda con velocità c.
- Forza di Lorentz: girazione circolare perpendicolare a B.
- Frequenza di plasma elettronica: le onde radio al di sotto non attraversano il plasma.
Esempio svolto
- Dati
- Protone del vento solare vicino alla Terra: v⊥ = 400 km/s, campo interplanetario B = 5 nT; densità elettronica n = 5 cm⁻³ (valori tipici).
- Incognita
- Raggio di girazione, periodo di girazione e frequenza di plasma elettronica.
- Strategia
- Inserire le costanti CODATA in r_g, ω_g e f_pe; tutto in SI (n = 5×10⁶ m⁻³).
- Passaggi
- 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
- Risposta
- r_g ≈ 835 km, periodo di girazione 13,1 s, f_pe ≈ 20,1 kHz.
- Controllo di plausibilità
- Il raggio di girazione è molto più piccolo della magnetosfera terrestre (decine di migliaia di km), quindi il plasma si comporta come un fluido legato al campo. La c ottenuta da μ₀ ed ε₀ riproduce il valore definito 299 792 458 m/s.
- Significato
- Queste scale decidono come il meteo spaziale raggiunge satelliti e reti elettriche, e perché la ionosfera riflette alcune frequenze radio.
Mostra il codice 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)Esercizi con correzione automatica
- m
- r_g = m_e v / (eB) con m_e = 9,109×10⁻³¹ kg.
Soluzione completa
r_g = 9,109e-31 · 1e6 / (1,602e-19 · 3e-5) = 0,190 m.
Esercizi
Perché ∇·B = 0 è importante fisicamente?
Mostra la soluzione
Non esistono monopoli magnetici: le linee di campo non iniziano né finiscono, si richiudono. I codici numerici di plasma devono rispettarlo, altrimenti producono forze false.
Errore frequente
Dimenticare i prefissi: 5 nT = 5×10⁻⁹ T, e 5 cm⁻³ = 5×10⁶ m⁻³.
Riferimenti
- 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).
Termodinamica e meccanica statistica
- Corso post-diploma
- Laurea triennale
- Laurea magistrale
Saprai
Legare la temperatura alla velocità molecolare con la distribuzione di Maxwell–Boltzmann; applicare primo e secondo principio; calcolare il rendimento di Carnot.
Prerequisiti
Legge dei gas perfetti, integrali gaussiani. La fisica spiegata
Prima l’intuizione
La temperatura è l’energia cinetica media per molecola. Contare i modi in cui le molecole si spartiscono l’energia (l’entropia) spiega perché un motore termico non può mai trasformare tutto il calore in lavoro.
🧒 Spiegamelo come se avessi 10 anni
L’aria è fatta di minuscole molecole che sfrecciano più veloci di un aereo di linea. «Caldo» vuol dire solo che vanno più veloci. E il calore va sempre dal caldo al freddo, mai al contrario da solo.
Passo dopo passo
- Distribuzione delle velocità di Maxwell–Boltzmann.
- Momenti della distribuzione (equipartizione).
- Primo principio (energia) e secondo principio (entropia).
- Nessun motore tra due sorgenti supera Carnot.
Esempio svolto
- Dati
- Azoto N₂ (massa molare 28,0134 g/mol) a 300 K (27 °C).
- Incognita
- Velocità quadratica media e media.
- Strategia
- Usare R = k_B N_A e M in kg/mol.
- Passaggi
- 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
- Risposta
- v_rms ≈ 517 m/s, v̄ ≈ 476 m/s.
- Controllo di plausibilità
- Il rapporto v_rms/v̄ = √(3π/8) ≈ 1,085, come per ogni gas di Maxwell–Boltzmann; entrambe superano un po’ la velocità del suono nell’aria (circa 347 m/s a 300 K), come previsto.
- Significato
- Le velocità molecolari fissano le velocità di reazione, quanto in fretta i gas fuoriescono e diffondono, e perché i gas leggeri sfuggono dalle atmosfere planetarie.
Mostra il codice 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)))Esercizi con correzione automatica
- m/s
- M = 0,0319988 kg/mol.
- Le molecole più pesanti sono più lente: v ∝ 1/√M.
Soluzione completa
v_rms = √(3RT/M) = 483,6 m/s (contro 516,8 m/s per N₂).
Esercizi
Un motore lavora tra 600 K e 300 K. Qual è il massimo rendimento possibile?
Mostra la soluzione
η_C = 1 − 300/600 = 0,50, cioè 50 %. I motori reali ottengono sensibilmente meno per attrito e scambi di calore a velocità finita.
Errore frequente
Usare grammi per mole invece di chilogrammi per mole, rendendo le velocità √1000 ≈ 32 volte troppo piccole.
Riferimenti
- 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).
Fluidi: Navier–Stokes e numero di Reynolds
- Laurea triennale
- Laurea magistrale
- Dottorato / ricerca
Saprai
Scrivere le equazioni di Navier–Stokes incomprimibili, adimensionalizzarle per trovare il numero di Reynolds e risolvere esattamente il flusso in condotta (Poiseuille).
Prerequisiti
Analisi vettoriale, seconda legge di Newton, il tema atmosfera del percorso aero. La fisica spiegata
Prima l’intuizione
Navier–Stokes è la seconda legge di Newton per una piccola porzione di fluido: pressione, viscosità e inerzia si sfidano. Che il flusso sia regolare o turbolento dipende dal rapporto tra inerzia e viscosità.
🧒 Spiegamelo come se avessi 10 anni
Il miele scende liscio, ma l’acqua da un rubinetto aperto al massimo schizza e fa vortici. Un solo numero, il numero di Reynolds, ti dice cosa succede: grande, vortici; piccolo, liscio.
Passo dopo passo
- Navier–Stokes incomprimibile e continuità.
- L’adimensionalizzazione lascia un solo parametro: il numero di Reynolds.
- Flusso stazionario in condotta: soluzione esatta (Hagen–Poiseuille), valida finché Re resta sotto circa 2300.
- Legge di Sutherland per la viscosità dell’aria (quella dell’atmosfera standard 1976).
Esempio svolto
- Dati
- Ala generica di aereo di linea in crociera: corda L = 5 m, U = 250 m/s, a 11 km: ρ = 0,3639 kg/m³, T = 216,65 K.
- Incognita
- Viscosità μ e numero di Reynolds Re.
- Strategia
- Sutherland per μ, poi Re = ρUL/μ.
- Passaggi
- μ = 1.458e-6 · 216.65^1.5 / (216.65 + 110.4) = 0,00001422 Pa·s
- Re = 0.3639 · 250 · 5 / μ
- Re ≈ 3.2 × 10⁷
- Risposta
- μ ≈ 1,42×10⁻⁵ Pa·s, Re ≈ 3,2×10⁷: lo strato limite è per lo più turbolento.
- Controllo di plausibilità
- L’atmosfera standard 1976 riporta μ = 1,4216×10⁻⁵ Pa·s a 11 km. Un Re di decine di milioni è tipico degli aerei reali, molto sopra i modelli in galleria del vento: ecco perché contano gli effetti di scala.
- Significato
- Re decide resistenza, transizione e separazione; riprodurlo è il problema centrale delle prove in galleria e della modellazione della turbolenza in CFD.
Mostra il codice 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)Esercizi con correzione automatica
-
- μ(288,15 K) da Sutherland ≈ 1,79×10⁻⁵ Pa·s.
Soluzione completa
Re = 1,225 · 30 · 0,2 / 1,789×10⁻⁵ ≈ 410757: circa 80 volte più piccolo dell’ala in crociera.
Esercizi
Raddoppiare il raggio di un tubo a pari caduta di pressione moltiplica la portata laminare per quanto?
Mostra la soluzione
Q ∝ R⁴, quindi per 16. Esempio: R = 0,5 mm, L = 10 cm, Δp = 1 kPa, acqua μ = 1,0 mPa·s dà Q = 0,245 mL/s.
Errore frequente
Usare la viscosità cinematica ν (m²/s) dove serve quella dinamica μ (Pa·s), o viceversa: Re = ρUL/μ = UL/ν.
Riferimenti
- F. M. White, Fluid Mechanics, 8th ed., McGraw-Hill, 2016.
- J. D. Anderson Jr., Fundamentals of Aerodynamics, 6th ed., McGraw-Hill, 2017.
Relatività per le orbite: orologi GPS e Mercurio
- Laurea magistrale
- Dottorato / ricerca
Saprai
Stimare le derive d’orologio relativistiche ristretta e generale di un satellite GPS e l’avanzamento relativistico del perielio di Mercurio.
Prerequisiti
Meccanica orbitale (velocità circolare, vis-viva), sviluppi in serie. La fisica spiegata
Prima l’intuizione
Gli orologi in moto rallentano (relatività ristretta) e quelli più in basso nella gravità rallentano (relatività generale). Per il GPS vince l’effetto della gravità, e i due insieme sono abbastanza grandi da far fallire il sistema in pochi minuti senza correzione.
🧒 Spiegamelo come se avessi 10 anni
Gli orologi dei satelliti GPS vanno un pochino più veloci di quelli a terra, perché lassù la gravità è più debole. Senza correggerlo, la tua mappa sbaglierebbe di circa 12 km al giorno.
Passo dopo passo
- Ritmo dell’orologio in campo debole: termine di potenziale (RG) più termine di velocità (RR).
- Satellite contro suolo, orbita circolare v² = GM/a (moto dell’orologio a terra trascurato qui).
- Avanzamento del perielio in RG dall’equazione dell’orbita di Schwarzschild.
Esempio svolto
- Dati
- Orbita GPS a = 26 561,75 km (circolare), R⊕ = 6378,137 km, GM⊕ = 3,986004418×10¹⁴ m³/s².
- Incognita
- Guadagno dell’orologio al giorno da RG e RR, il totale e l’errore di distanza equivalente.
- Strategia
- Calcolare i due spostamenti relativi, moltiplicare per 86 400 s, poi il tempo netto per c.
- Passaggi
- GR: 45,65 μs/day
- SR: -7,21 μs/day
- net: 38,44 μs/day
- c · 38,44 μs = 11,5 km/day
- Risposta
- Circa +45,7 μs/giorno (gravità) e -7,2 μs/giorno (velocità): netto ≈ +38,4 μs/giorno, cioè circa 12 km di errore di distanza al giorno.
- Controllo di plausibilità
- Ashby (2003), che include rotazione e schiacciamento terrestri, dà circa +45,7 e −7,1 μs/giorno, netto 38,6 μs/giorno: il nostro modello semplice concorda entro 0,2 μs/giorno.
- Significato
- Gli orologi dei satelliti GPS vengono regolati di proposito per andare più lenti prima del lancio, così in orbita restano allineati con quelli a terra. La relatività è una correzione ingegneristica quotidiana.
Mostra il codice 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)Esercizi con correzione automatica
- ″/century
- Per orbita: 6πGM☉ / (c²a(1 − e²)), con GM☉ = 1,32712440018×10²⁰ m³/s².
- Orbite per secolo = 36 525 / 224,701; converti i radianti in secondi d’arco.
Soluzione completa
Δϖ ≈ 8,62″ per secolo per Venere, contro 42,98″ per Mercurio (Will 2014 riporta 42,98″ per Mercurio).
Esercizi
Perché la RG cambia l’orbita di Mercurio più di quella della Terra?
Mostra la soluzione
L’avanzamento per orbita va come 1/(a(1−e²)), e Mercurio compie circa 4 volte più orbite per secolo, quindi l’effetto per secolo è molto maggiore vicino al Sole.
Errore frequente
Sommare i termini RR e RG con lo stesso segno. La velocità del satellite rallenta l’orologio (negativo), la quota lo accelera (positivo).
Riferimenti
- 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).
Verifica la comprensione
Approfondisci
Poi: Taylor e Goldstein per la meccanica, Griffiths e Chen per campi e plasma, Schroeder per la fisica termica e la rassegna open access di Ashby sulla relatività nel GPS.
Per insegnanti e docenti
Uso suggerito: una lezione di 50 minuti per sezione. Percorso post-diploma/triennale: Newton e conservazione → proiettili con resistenza (gli studenti modificano il Python RK4 e confrontano col Physics Lab) → oscillatori → termodinamica. Percorso magistrale: meccanica lagrangiana → elettromagnetismo e plasma spaziale → Navier–Stokes e numero di Reynolds → relatività per il GPS. Ogni risposta numerica è calcolata dallo script a partire dalle costanti CODATA 2018 e funge da chiave. Allineato ai temi MIT OCW 8.01/8.02.
Riferimenti
- 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.