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

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

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.

Provalo nel Physics Lab

Leggi di Newton e conservazione

  • Corso post-diploma
  • Laurea triennale

Percorso: fisica · Circa 25 min · Numeri verificati in fase di generazione:

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

  1. F=dpdt,p=mv\mathbf F=\frac{d\mathbf p}{dt},\qquad \mathbf p=m\mathbf v
    Seconda legge in forma di quantità di moto.
  2. F12=−F21  ⇒  ddt(p1+p2)=0\mathbf F_{12}=-\mathbf F_{21}\;\Rightarrow\;\frac{d}{dt}(\mathbf p_1+\mathbf p_2)=0
    Terza legge: le forze interne si annullano, la quantità di moto totale si conserva.
  3. W=∫F⋅dr=Δ(12mv2)W=\int\mathbf F\cdot d\mathbf r=\Delta\left(\tfrac12mv^2\right)
    Teorema dell’energia cinetica.
  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}
    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
  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
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

  1. m/s

    Soluzione completa

    v₂ = 2·1·4/4 = 2,0 m/s; il carrello leggero rimbalza a -2,0 m/s.

Esercizi

  1. 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

Percorso: fisica · Circa 40 min · Numeri verificati in fase di generazione:

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

  1. mv˙=mg−12ρCDA ∣v∣ vm\dot{\mathbf v}=m\mathbf g-\tfrac12\rho C_DA\,|\mathbf v|\,\mathbf v
    Resistenza quadratica, opposta alla velocità.
  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}
    Accoppiate tramite |v|: nessuna soluzione elementare.
  3. vt=2mgρCDAv_t=\sqrt{\frac{2mg}{\rho C_DA}}
    Velocità limite: la resistenza bilancia il peso.
  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’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
  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
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

  1. m

    Soluzione completa

    R(45°) = 94,8 m, contro 93,9 m a 35° e 163,2 m nel vuoto a 45°.

Esercizi

  1. 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

Percorso: fisica · Circa 35 min · Numeri verificati in fase di generazione:

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

  1. mx¨+cx˙+kx=0m\ddot x+c\dot x+kx=0
    Legge di Newton con forza elastica e smorzamento viscoso.
  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}}
    Equazione caratteristica.
  3. s=−ζωn±i ωn1−ζ2(ζ<1)s=-\zeta\omega_n\pm i\,\omega_n\sqrt{1-\zeta^2}\quad(\zeta<1)
    Sottosmorzato: coppia complessa; parte reale = tasso di decadimento, parte immaginaria = pulsazione smorzata.
  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}
    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
  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
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

    Soluzione completa

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

Esercizi

  1. 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

Percorso: fisica · Circa 45 min · Numeri verificati in fase di generazione:

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

  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
    Il principio di Hamilton dà l’equazione di Eulero–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
    Pendolo, con l’angolo come coordinata.
  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})}
    Periodo esatto dalla conservazione dell’energia; la media aritmetico-geometrica calcola in fretta l’integrale ellittico 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}
    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
  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
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

  1. s

    Soluzione completa

    T₀ = 2,8375 s; T = 3,3492 s (il 18 % in più).

Esercizi

  1. 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

Percorso: fisica · Circa 40 min · Numeri verificati in fase di generazione:

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

  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}
    Equazioni di 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}}
    Nel vuoto si prende il rotore della legge di Faraday: equazione d’onda con velocità 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}
    Forza di Lorentz: girazione circolare perpendicolare a B.
  4. fpe=12πne2ε0mef_{pe}=\frac{1}{2\pi}\sqrt{\frac{n e^2}{\varepsilon_0 m_e}}
    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
  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
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

  1. m

    Soluzione completa

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

Esercizi

  1. 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

Percorso: fisica · Circa 35 min · Numeri verificati in fase di generazione:

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

  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)}
    Distribuzione delle velocità di 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}}
    Momenti della distribuzione (equipartizione).
  3. dU=δQ−δW,dS≥δQTdU=\delta Q-\delta W,\qquad dS\ge\frac{\delta Q}{T}
    Primo principio (energia) e secondo principio (entropia).
  4. η≤ηC=1−TcTh\eta\le\eta_C=1-\frac{T_c}{T_h}
    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
  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
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

  1. m/s

    Soluzione completa

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

Esercizi

  1. 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

Percorso: fisica · Circa 40 min · Numeri verificati in fase di generazione:

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

  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 incomprimibile e 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’adimensionalizzazione lascia un solo parametro: il numero di 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}
    Flusso stazionario in condotta: soluzione esatta (Hagen–Poiseuille), valida finché Re resta sotto circa 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}
    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. μ = 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⁷
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

    Soluzione completa

    Re = 1,225 · 30 · 0,2 / 1,789×10⁻⁵ ≈ 410757: circa 80 volte più piccolo dell’ala in crociera.

Esercizi

  1. 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

Percorso: fisica · Circa 45 min · Numeri verificati in fase di generazione:

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

  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}
    Ritmo dell’orologio in campo debole: termine di potenziale (RG) più termine di velocità (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 contro suolo, orbita circolare v² = GM/a (moto dell’orologio a terra trascurato qui).
  3. Δϖ=6πGM⊙c2a(1−e2) per orbit\Delta\varpi=\frac{6\pi GM_\odot}{c^2a(1-e^2)}\ \text{per orbit}
    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
  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
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

  1. ″/century

    Soluzione completa

    Δϖ ≈ 8,62″ per secolo per Venere, contro 42,98″ per Mercurio (Will 2014 riporta 42,98″ per Mercurio).

Esercizi

  1. 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

Un oscillatore smorzato ha autovalori −0,5 ± 3i. È…

In che verso derivano in totale gli orologi GPS rispetto a quelli a terra?

Raddoppiare la temperatura assoluta di un gas moltiplica la velocità quadratica media per…

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

Aiuta / Dona