Calcolatore flusso disco rotante von Karman

Il calcolatore flusso del disco rotante di von Karman determina lo spessore dello strato limite delta, il momento torcente M e la velocità assiale di suzione w_inf per un disco infinito che ruota con velocità angolare Omega attorno al proprio asse z immerso in un fluido inizialmente in quiete, secondo la soluzione esatta delle equazioni di Navier-Stokes 3D incomprimibili trovata dal fisico ungaro-americano Theodore von Karman 1921 sul ZAMM (Zeitschrift fur Angewandte Mathematik und Mechanik) vol. 1, p. 233-252 "Uber laminare und turbulente Reibung". La variabile di similarita: zeta = z sqrt(Omega/nu) riduce il sistema PDE 3D in ODE 1D in zeta. I profili di velocità: u_r = r Omega F(zeta) radiale, u_theta = r Omega G(zeta) tangenziale, u_z = sqrt(nu Omega) H(zeta) assiale, con F(0)=0, G(0)=1 (no-slip parete), F(inf)=0, G(inf)=0 (fluido fermo all infinito), H(0)=0. Le 3 formule chiave:

• Spessore strato limite delta = 5,5 sqrt(nu/Omega) [m] (definizione 99% velocità esterna)
• Momento torcente per faccia M = 0,616/2 rho sqrt(nu) Omega^(3/2) R^4 [N m] (coefficiente Cochran 1934 numerico)
• Velocità assiale all infinito w_inf = -0,886 sqrt(nu Omega) [m/s] (suzione verso il disco, segno negativo). La fisica: forza centrifuga sulla parete espelle il fluido radialmente verso l esterno (effetto pompa centrifuga) e per continuità richiama assialmente fluido dall infinito (effetto suzione w_inf). Lo strato limite ha spessore costante indipendente da r (caso unico in fluidodinamica) perché il bilancio centrifuga vs viscosa non dipende da r. Il numero di Reynolds locale Re = r^2 Omega/nu, con transizione a turbolenza per Re ~ 3 10^5 (dischi grandi/veloci). Le 4 categorie di applicazioni:

• Storage HDD/floppy Omega = 300-1000 rad/s aria, delta = 1-2 mm (ora deprecato a favore SSD STMicroelectronics Catania-Agrate)
• Centrifughe biotech-laboratorio Omega = 50-500 rad/s acqua/buffer, delta = 0,3-0,8 mm (IIT Genova, IRCCS, ENEA Casaccia)
• Turbomacchine assiali Omega = 200-3000 rad/s gas/vapore (Ansaldo Energia Genova turbine, Avio SpA Colleferro motori Vega)
• Aerodinamica rotori elicotteri-EV Omega = 30-2000 rad/s (Leonardo Helicopters Cascina Costa, Stellantis EV Mirafiori, Ferrari ePower Maranello). La validazione CFD canonica: il flusso di von Karman e caso benchmark per OpenFOAM, ANSYS Fluent, SU2 usato dal CINECA Marconi100 + Leonardo HPC pre-exascale 250 PFlops 2023 per validare codici turbolenti k-omega SST e LES presso Politecnico Milano 1863, Politecnico Torino 1859, UniNa Federico II Napoli ingegneria aerospaziale

Tabelle disco rotante italiana

(Fonte: von Karman 1921 ZAMM 1, 233 + Cochran 1934 Proc. Cambridge Philos. Soc. 30, 365 coefficiente 0,616 + dati Politecnico Milano 1863 + Leonardo Helicopters Cascina Costa Vergiate 1948).

Esempi calcolo italiano disco rotante

• HDD 7200 rpm STMicroelectronics SGS-Thomson Catania 1987 + Agrate Brianza pre-deprecation SSD: Omega = 7200 2 pi/60 = 754 rad/s, nu_aria = 1,5e-5 m^2/s a 25 C, R = 0,047 m (3,5"), rho_aria = 1,225 kg/m^3. delta = 5,5 sqrt(1,5e-5/754) = 5,5 sqrt(1,99e-8) = 7,76e-4 m = 0,78 mm. Momento torcente M = 0,616 1,225 sqrt(1,5e-5) (754)^1,5 (0,047)^4/2 = 1,9e-6 N m. w_inf = -0,886 sqrt(1,5e-5 754) = -0,886 sqrt(0,0113) = -0,094 m/s (suzione 9 cm/s). Hai bisogno solo 2 mm clearance head-disk (legacy MR head IBM 1991, ora deprecato GMR + SSD NAND STMicroelectronics Catania-Agrate)
• Centrifuga biotech IIT Genova 2003 microfluidica purificazione DNA: Omega = 100 rad/s (1000 rpm), nu_acqua = 1,0e-6 m^2/s, R = 0,1 m, rho = 1000 kg/m^3. delta = 5,5 sqrt(1e-6/100) = 5,5e-4 m = 0,55 mm. M = 0,616 1000 sqrt(1e-6) (100)^1,5 (0,1)^4/2 = 0,616 1000 0,001 1000 1e-4/2 = 0,031 N m. w_inf = -0,886 sqrt(1e-6 100) = -0,886 0,01 = -8,86e-3 m/s (suzione 9 mm/s pellet formation). Applicazione IRCCS ENEA Casaccia 1962 sequenziamento DNA
• Turbina idraulica pala disco Ansaldo Energia Genova 1853: Omega = 300 rad/s (3000 rpm), nu_acqua = 1,0e-6 m^2/s, R = 0,5 m, rho = 1000 kg/m^3. delta = 5,5 sqrt(1e-6/300) = 3,17e-4 m = 0,32 mm. M = 0,616 1000 sqrt(1e-6) (300)^1,5 (0,5)^4/2 = 0,616 1000 0,001 5196 0,0625/2 = 10,0 N m. w_inf = -0,886 sqrt(1e-6 300) = -0,015 m/s. Significativo per dimensionamento boundary layer pala disco impulse Pelton turbines centrali idroelettriche Enel Edison
• Compressore aria assiale Avio SpA Colleferro 1908 motore Vega launcher ESA: Omega = 2000 rad/s, nu_aria = 1,5e-5, R = 0,2 m, rho = 1,225. delta = 5,5 sqrt(1,5e-5/2000) = 4,77e-4 m = 0,48 mm. M = 0,616 1,225 sqrt(1,5e-5) (2000)^1,5 (0,2)^4/2 = 1,04e-3 N m. w_inf = -0,886 sqrt(1,5e-5 2000) = -0,153 m/s. Avio Vega P80 booster solid + AVUM upper stage liquid LOX/RP1
• Rotore elicottero Leonardo Helicopters AW139 Cascina Costa Vergiate 1948: Omega = 30 rad/s (286 rpm), nu = 1,5e-5, R = 5,5 m, rho = 1,225. delta = 5,5 sqrt(1,5e-5/30) = 3,89e-3 m = 3,9 mm (boundary layer pala). M = 0,616 1,225 sqrt(1,5e-5) (30)^1,5 (5,5)^4/2 = 11,9 N m (per faccia disco virtuale; momento totale rotore 4 pale e ben superiore). w_inf = -0,886 sqrt(1,5e-5 30) = -1,88e-2 m/s. Modello disk actuator Glauert 1935 base downwash teoria momentum elicottero
• Motore brushless EV Stellantis Mirafiori Torino 1899 EV ex FIAT: Omega = 1500 rad/s (14.000 rpm), nu_aria = 1,5e-5, R = 0,08 m, rho = 1,225. delta = 5,5 sqrt(1,5e-5/1500) = 5,5e-4 m = 0,55 mm. M_drag = 3,1e-5 N m (perdita ventilazione air-gap statore-rotore). Applicazione Stellantis EV 500 + Ferrari ePower Maranello 1947 motori F1 ibridi MGU-K MGU-H
• Centrifuga laboratorio analisi Carlo Rubbia Nobel 1984 INFN Frascati 1955 + ENEA Casaccia 1962: Omega = 500 rad/s, nu_acqua = 1e-6, R = 0,15 m, rho = 1000. delta = 0,25 mm. M = 0,79 N m. w_inf = -0,020 m/s. Sample preparation isotopi medicale RT/PET radiofarmaci

Ricerca italiana fluidodinamica disco rotante

La ricerca italiana sulla fluidodinamica del disco rotante e turbomacchine:

• Politecnico Milano 1863 Dipartimento Energia Bovisa CFD turbomacchine k-omega SST + LES validazione benchmark von Karman; Politecnico Torino 1859 DIMEAS aerospaziale + Avio SpA Colleferro 1908 motori Vega P80 + AVUM Zefiro 23/9; UniNa Federico II Napoli 1224 ingegneria aerospaziale rotori elicotteri Leonardo
• Leonardo SpA ex Finmeccanica Roma 1948 + Leonardo Helicopters Cascina Costa Vergiate Varese 1948 ex Agusta AW139 (medium twin), AW169, AW109, NH90, AW609 tiltrotor: applicazione disco actuator + boundary layer pala rotore
• Ansaldo Energia Genova Sampierdarena 1853 turbine a gas GT26 GT36 H-Class 538 MW, turbine vapore impianti Enel Edison + nucleare ex CISE 1946 (programma Italia 1968-1987 Caorso Trino Latina Garigliano deprecato post-Chernobyl 1987 referendum); Avio SpA Colleferro 1908 ex Fiat Avio motori aeronautici turbofan (GE-Avio collaborazione GE9X Boeing 777X) + propulsione spaziale ESA Vega launcher
• STMicroelectronics ex SGS-Thomson Catania 1987 + Agrate Brianza Milano semiconduttori MEMS accelerometri smartphone + microcontrollori automotive STM32 + memoria SSD NAND deprecata HDD legacy + sensori Hall brushless EV motors
• IIT Istituto Italiano di Tecnologia Genova 2003 microfluidica lab-on-chip + biotecnologie centrifuge bioreactor; IRCCS Ospedali ricerca + ENEA Casaccia Roma 1962 + Frascati 1955 INFN + DTT Divertor Tokamak Test 2028 Frascati fusione nucleare boundary layer plasma divertor; SISSA Trieste 1978 + ICTP Trieste 1964 + Cineca Casalecchio Bologna 1969 Marconi100 30 PFlops + Leonardo HPC pre-exascale 250 PFlops 2023 (5 TOP500) simulazioni CFD turbolenza dischi rotanti DNS/LES validation case canonico OpenFOAM ANSYS Fluent SU2
• Stellantis ex FCA-PSA Mirafiori Torino 1899 ex FIAT motori EV + endotermici diesel benzina; Ferrari Maranello 1947 F1 motori ibridi MGU-K kinetic + MGU-H heat ePower 1.6L V6 turbo + Ducati Bologna 1926 motori MotoGP Desmosedici
• Eurofighter Typhoon Italia partnership 36% (insieme UK 33% Germania 21% Spagna 10%) Aeronautica Militare AMI + Leonardo + GE-Avio motori EJ200 turbofan; CAMM-X meteo radar metereologia AMI Pratica di Mare
• INRIM Torino 1923 metrologia Istituto Nazionale Ricerca Metrologica + INFN Frascati 1955 acceleratori DAFNE + LNGS Gran Sasso 1989 Borexino XENONnT dark matter + AGCM Antitrust 1990
• Premi Nobel italiani fisica: Marconi 1909, Fermi 1938, Segre 1959, Carlo Rubbia 1984 CERN W/Z bosons, Giacconi 2002 raggi X astrofisica, Parisi 2021 sistemi complessi disordinati
• PNRR Piano Nazionale Ripresa Resilienza Missione 2 transizione ecologica idrogeno verde + Missione 4 ricerca innovazione 6,9 mld euro. Bonus fiscali italiani ricerca-industria 2026: Credito d Imposta R&S art. 1 c. 200 L. 160/2019 fino 25% spese personale + macchinari + brevetti; Patent Box L. 178/2020 (DL 146/2021) deduzione 110% costi qualificati IP brevetti software know-how; Credito d Imposta Industria/Transizione 5.0 art. 38 DL 39/2024 fino 35% beni 4.0 + risparmio energetico 6%. Il flusso di von Karman resta il caso canonico di soluzione esatta Navier-Stokes 3D con applicazioni industriali italiane diversificate dal calcolo HPC alle turbine ai rotori

Aggiornato a luglio 2026