Engineering
Volo a vela nelle termiche
Vola con un aliante bianco da 15 m sopra un’ampia valle alpina in un pomeriggio estivo vigoroso: spirala in una termica sotto il suo cumulo, sali verso la base della nube, poi plana fino alla nube successiva e sali di nuovo. La polare delle velocità a V³ + b/V dà le velocità di minima discesa e di massima efficienza e un’efficienza massima di 40; in una termica gaussiana l’inclinazione che sale meglio ha una forma chiusa tramite la funzione W di Lambert: troppo poco inclinato e il cerchio resta fuori dal nucleo, troppo inclinato e la discesa in virata si mangia l’ascendenza. Tra una termica e l’altra la velocità di MacCready è la tangente alla polare tracciata dalla salita attesa, e la velocità di trasferimento misurata segue V c/(c + w_s). Disattiva la centratura e un cerchio fuori centro scende fuori dall’ascendenza; indebolisci la termica e l’aliante atterra in un campo. Ogni joule di quota e di velocità è seguito dall’aria ascendente fino alla resistenza.

Cosa succede
Un aliante da 15 m e 400 kg (10,5 m², allungamento 21,4, fattore di Oswald 0,92, C_D0 = 0.0095, C_Lmax = 1.6) è una massa puntiforme in aria calma di densità 1,1 kg/m³. La portanza ½ρV²S·C_L è perpendicolare alla velocità rispetto all’aria, inclinata di φ attorno a essa; la resistenza ½ρV²S·(C_D0 + k·C_L²), k = 1/(π·e·AR). In regime stazionario questa è la polare w_s = a·V³ + b/V, a = ρS·C_D0/2W, b = 2kW/ρS: minima discesa 0,65 m/s a 81 km/h, efficienza massima 40,4 a 107 km/h; in virata il raggio è V²/(g·tan φ) e la discesa a·V³ + b/(V·cos²φ), con l’elica esatta (sin γ = D/W, L·cos φ = W·cos γ) risolta accanto. Tre termiche gaussiane w = W0·exp(−r²/R²), da 3,5 m/s e 150 m, si attenuano come 1 − ((z − z1)/300 m)² sotto la base delle nubi a 2000 m. L’inclinazione migliore ha forma chiusa tramite la funzione di Lambert: u = −2·W(−√(b·q/(V·W0))/2), q = V^4/(gR)², tan φ = √(q/u), 47,8° e +1,76 m/s a 94 km/h; una salita concentrica si ferma a z1 + 300 m·√(1 − s/w(r)) = 1932 m. La velocità di MacCready risolve 2a·V^4 − c·V − 2b = 0 con il metodo di Ferrari, 157 km/h con 1,8 m/s, e la velocità di trasferimento è V·c/(c + w_s). L’autopilota imposta C_L per un g alla velocità attuale più un PI sull’errore rispetto a una velocità di riferimento che varia a 4 m/s², rolla fino a 34°/s, tiene il C_L della virata livellata sotto il 97 % di C_Lmax, centra riducendo l’inclinazione man mano che l’ascendenza netta cresce lungo la traiettoria e punta verso la nube successiva. Sottopassi di 20 ms tengono costanti C_L e l’inclinazione, prendono la forza nel punto medio e avanzano con la velocità media, così che m·g·z + m·v²/2 varia esattamente del lavoro dell’aria meno il lavoro della resistenza, a meno dell’arrotondamento; la quota di uscita, l’ingresso in termica, il traguardo e un atterraggio sono localizzati dentro il passo per bisezione.
Da provare
Guarda la prova di distanza: la salita all’inclinazione ottimale, la planata alla velocità di MacCready e la velocità della tratta rispetto alla previsione. Passa a Termica e imposta l’inclinazione a mano da 25° a 60°; disattiva la centratura con uno scostamento; vola in Planata a diverse velocità di crociera e leggi l’efficienza; aggiungi zavorra e guarda le velocità crescere con √m. Tieni la quota di uscita sotto il tetto.
Il modello e i suoi limiti
Una massa puntiforme senza dinamica di rotazione: l’autopilota imposta l’inclinazione tramite un limite di velocità di rollio e C_L direttamente; una polare parabolica senza l’aumento di resistenza ad alto C_L; densità dell’aria costante e niente vento, così che le termiche non derivano (rimandato a una fase successiva); termiche come colonne gaussiane stazionarie di intensità fissa senza discendenze attorno né turbolenza; un variometro netto ideale. Le montagne sono scenografia che il volo non incontra mai. La valle è disegnata in scala 1:300, l’aliante quattro volte più grande. Un modello didattico, non un simulatore di volo.
Equazioni del modulo
Equazioni w_s(V) = a·V³ + b/V, a = ρ·S·C_D0/(2W), b = 2k·W/(ρ·S), V_ms = (b/3a)^(1/4), V_bg = (b/a)^(1/4), (L/D)max = 1/(2√(C_D0·k)) c(φ) = W0·exp(−r²/R²) − a·V³ − b/(V·cos²φ), r = V²/(g·tan φ), tan φ* = √(q/u), u = −2·W(−√(b·q/(V·W0))/2), q = V^4/(g·R)² m·dv/dt = L + D + m·g, |L| = ½ρ·V_a²·S·C_L across v_a banked by φ, |D| = ½ρ·V_a²·S·(C_D0 + k·C_L²), v_a = v − w(x)·e_z 2a·V^4 − c·V − 2b = 0 (MacCready), V_avg = V·c/(c + w_s(V)), Δ(m·g·z + ½m·v²) = ∫F·w dt − ∫F·(v − w) dt


