Come calcolare il numero di Prandtl (Pr): numero adimensionale fondamentale meccanica fluidi e trasmissione calore, rapporto tra diffusività di momento (viscosità cinematica) e diffusività termica del fluido. Formula: Pr = nu / alpha = (mu * c_p) / k, dove mu viscosità dinamica (Pa·s), c_p calore specifico a pressione costante (J/kg·K), k conducibilità termica (W/m·K), nu = mu/rho viscosità cinematica (m²/s), alpha = k/(rho*c_p) diffusività termica (m²/s). Il numero Prandtl è proprietà del fluido stesso (non del flusso), dipende solo dalla composizione e temperatura. Pr basso (Pr momento, strato limite termico più spesso del cinematico, calore si propaga rapidamente (gas, metalli liquidi). Pr alto (Pr > 1): momento si propaga più rapidamente, strato limite cinematico più spesso (liquidi viscosi, oli). Storia: Ludwig Prandtl (Germania 1875-1953) Göttingen Universität 1904 boundary layer theory paper III International Mathematical Congress Heidelberg 'Über Flüssigkeitsbewegung bei sehr kleiner Reibung' fondamentale paper introduzione concetto strato limite. Prandtl ricerca aerodinamica AVA Aerodynamische Versuchsanstalt Göttingen 1907 fondata. Strumento essenziale per: trasferimento calore convettivo ingegneria heat exchangers shell-and-tube counter-current flow;

HVAC progettazione heating ventilation air conditioning case impianti; aerodinamica industria automotive Italia FCA Stellantis Torino + Maserati + Ferrari + Lamborghini Bologna; industria aerospace Italia Avio Colleferro Vega rocket + Leonardo Roma 7000 employees + Piaggio Aero Genova; petrolchimica raffinerie Eni Milano Saipem; industria alimentare processi termici sterilizzazione Parma food valley; centrali nucleari Italia ENEA Frascati ITER fusion contributi; fluidi geotermici Italia Larderello Pisa centrali geotermiche prima centrale geotermoelettrica mondiale 1904 Piero Ginori Conti. Italia ricerca termofluidodinamica: Politecnico Milano DEIB Bioengineering DEAL Aerospaziale 1863 fondato; Politecnico Torino DET Energia 1859 fondato; Sapienza Roma DIETI Ingegneria Industriale 1303 fondata; Federico II Napoli DIETI Ingegneria Industriale 1224 fondata; Padova DEI Ingegneria Industriale 1222 Galileo professor; Bologna alma mater 1088 prima Europa Faculty Engineering DIN; Pisa SNS Scuola Normale Superiore 1810 Mathematics. Italia istituzioni ricerca: CNR Consiglio Nazionale Ricerche Roma 1923 8000 ricercatori 200 institutes ITC Costruzioni Bologna;

ENEA Roma + Frascati + Bologna energia ambiente fusione nucleare; ASI Agenzia Spaziale Italiana Roma 1988; INFN Frascati 1951 5000 ricercatori; Avio Colleferro propulsione spaziale Vega rocket; Leonardo Roma 1948 Defense Aerospace 50000 employees; Piaggio Aerospace Genova 1884. Italia industria automotive: FCA Stellantis Torino Mirafiori 2021 fusion FCA + PSA Peugeot 300000 employees worldwide; Ferrari Maranello 1947 4500 employees luxury sport cars; Lamborghini Sant Agata Bolognese 1963 Volkswagen subsidiary 2000 employees; Maserati Modena 1914 Stellantis luxury 2000 employees; Pagani Modena 2002 hypercar; Ducati Bologna 1926 motorcycles Volkswagen subsidiary 1500 employees; Aprilia Noale 1945 Piaggio Group; MV Agusta Schiranna Varese 1945; Iveco Trucks Torino 1975 CNH Industrial. Aria Pr 25°C = 0.713; acqua 25°C = 6.13; mercurio 25°C = 0.0248; olio motore 100°C = 1500; olio idraulico = 100; glicerolo = 12000; etanolo = 16; benzene = 7.

Formula numero di Prandtl

Definizione fondamentale: Pr = nu / alpha = (mu * c_p) / k. Numero adimensionale (no unità). nu viscosità cinematica m²/s = mu/rho dove mu Pa·s (Newton·secondo/m² = kg/m/s) viscosità dinamica e rho kg/m³ densità; alpha diffusività termica m²/s = k/(rho*c_p) dove k W/m/K conducibilità termica e c_p J/kg/K calore specifico isobaro.

Significato fisico: Pr rapporto tra velocità diffusione momento (transporto quantità moto) e velocità diffusione calore (transporto energia termica). Pr=1 entrambi processi diffondono ugualmente (Reynolds analogy). Pr<>1 momento diffonde rapidamente vs calore (oli viscosi).

Strato limite cinematico vs termico: relazione delta_thermal/delta_velocity ~ Pr^(-1/3) per flussi laminari. Pr delta_velocity strato termico più spesso; Pr=1 uguali; Pr>1 strato cinematico più spesso. Importante per design heat exchangers (water Pr=7) vs air conditioners (air Pr=0.71).

Correlazioni heat transfer: Nusselt number Nu = h*L/k dipende Re e Pr. Convezione laminare placca: Nu = 0.332 * Re^0.5 * Pr^(1/3); turbolenta tubo: Dittus-Boelter Nu = 0.023 * Re^0.8 * Pr^n (n=0.4 heating, 0.3 cooling). Sieder-Tate Nu = 0.023 * Re^0.8 * Pr^(1/3) * (mu/mu_w)^0.14 viscosity correction.

Variazione con temperatura: Pr dipende fortemente da T. Acqua: Pr=13.6 a 0°C; 7.0 a 20°C; 1.75 a 100°C; 0.86 a 300°C. Aria: Pr=0.71 stabile 0-1000°C. Olio motore: Pr=10000 a 0°C; 1500 a 100°C; 220 a 200°C. Mercurio: Pr=0.025 stabile metalli liquidi.

Storia e contesto

Ludwig Prandtl (Germania 1875-1953) Freising Bavaria nato, TU München mathematics 1893-1897, Hannover Universität engineering doctorate 1900 'Kipperscheinungen' lateral buckling thin-walled beams. Hannover Hochschule Technische 1900-1904 lecturer; Göttingen Universität full professor mechanics 1904-1947 chair Applied Mathematics Mechanics 43 anni carriera Göttingen direttore Institut für Angewandte Mathematik und Mechanik.

1904 paper III International Mathematical Congress Heidelberg August 1904 'Über Flüssigkeitsbewegung bei sehr kleiner Reibung' (Fluid Motion at Very Small Viscosity) 8 pagine fondamentale boundary layer theory rivoluzione meccanica fluidi separazione flusso D'Alembert paradox solved attaching wall flow vs free flow. AVA Aerodynamische Versuchsanstalt Göttingen 1907 fondata Prandtl direttore wind tunnel ricerca aerodinamica + Kaiser-Wilhelm-Institut für Strömungsforschung 1925 Göttingen. Prandtl studenti famosi: Theodore von Kármán (Hungary-USA 1881-1963) Caltech aerodinamica supersonic; Adolf Busemann (Germania 1901-1986) NASA Ames swept wing concept Mach cone; Albert Betz (Germania 1885-1968) Göttingen Betz limit turbine wind 59.3% maximum efficiency 1919 paper; Heinrich Blasius (Germania 1883-1970) TH Hamburg Blasius solution boundary layer 1908 PhD Göttingen Prandtl supervisor; Hermann Schlichting (Germania 1907-1982) Braunschweig Boundary Layer Theory textbook 1951 classic 8 editions; Karl Pohlhausen (Germania 1892-1980) Nuremberg; Walter Tollmien (Germania 1900-1968) Tollmien-Schlichting waves transition turbulence; Hermann Schultz-Grunow (Germania 1898-1976); Lothar Prandtl figlio Wilhelm Prandtl chemist.

Prandtl boundary layer textbook 1922 'Führer durch die Strömungslehre' classic German engineering textbook.

Concetti Prandtl: boundary layer 1904 separation transition turbulence; mixing length theory 1925 turbulent shear flow Prandtl mixing length; Prandtl number naming Schmidt 1925; Prandtl-Glauert compressibility correction 1927 subsonic compressibility; Lifting line theory 1918 finite wing aerodynamics induced drag; Prandtl-Meyer expansion fan supersonic flow.

Prandtl number nomenclatura: Ernst Schmidt (Germania 1892-1975) Munich 1925 paper introdusse termine 'Prandtl-Zahl' onorando Prandtl. Schmidt anche definì Schmidt number Sc = nu/D mass diffusivity. Numero adimensionale fluidodinamica heat transfer onore Ludwig Prandtl maestro.

Sviluppi successivi: Reynolds analogy 1874 Osborne Reynolds (UK 1842-1912) Manchester analogie momentum heat mass transfer Pr=Sc=1 ideale; Chilton-Colburn analogy 1934 j-factor for non-unity Pr Sc; Dittus-Boelter 1930 Berkeley Boelter empirical correlation tube flow; Sieder-Tate 1936 viscosity correction; Petukhov 1970 Russia improved correlation; Gnielinski 1976 Germany.

Italia ricerca termofluidodinamica: Politecnico Milano DEIB Bioengineering 1863 fondato 4500 employees Faculty Engineering Industrial; DEAL Dipartimento Energia Aerospaziale Politecnico Milano; Politecnico Torino DET Dipartimento Energia 1859 fondato Mirafiori automotive research center; Sapienza Roma DIETI Ingegneria Industriale 1303 fondata oldest Roma Faculty Engineering Industrial; Federico II Napoli DIETI 1224 fondata Federico II Hohenstaufen emperor Faculty Engineering; Padova DEI Ingegneria Industriale 1222 fondata Galileo professor 1592-1610; Bologna alma mater 1088 prima Europa Faculty Engineering DIN; Pisa SNS Scuola Normale Superiore 1810 Mathematics. CNR Consiglio Nazionale Ricerche Roma 1923 fondato 8000 ricercatori 200 institutes ITC Istituto Tecnologie Costruzione Bologna + INM Istituto Marina Genova ship hydrodynamics; ENEA Roma + Frascati + Bologna 1982 energia ambiente fusione nucleare DTT Divertor Tokamak Test Frascati; ASI Agenzia Spaziale Italiana Roma 1988; INFN Frascati 1951 5000 ricercatori; Avio Colleferro propulsione spaziale Vega rocket Solid Rocket Motor; Leonardo Roma 1948 Defense Aerospace 50000 employees worldwide F-35 partner Italy; Piaggio Aerospace Genova 1884 P-180 Avanti regional aircraft.

Italia geotermia storia: Larderello Pisa Toscana 1904 Piero Ginori Conti prima centrale geotermoelettrica mondiale 1904 5 lampadine accese, ENEL Larderello geotermica 800 MW capacity 30% energia geotermica mondiale Toscana centro mondo geotermia. Italia volcanologia INGV Istituto Nazionale Geofisica Vulcanologia Roma Etna Vesuvio Stromboli.

USA heat transfer research: MIT Cambridge MA mechanical engineering; Stanford Mechanical Palo Alto CA; Caltech Pasadena CA; UC Berkeley CA Lawrence Berkeley Lab; Princeton NJ; Yale CT. USA companies: GE General Electric Boston MA 1892 turbines aerospace; Westinghouse Pittsburgh PA 1886 nuclear; Bechtel San Francisco CA 1898 engineering; Fluor Irving TX 1912 engineering; Honeywell Charlotte NC 1906 industrial. Software heat transfer: ANSYS Canonsburg PA 1970 finite element CFD Fluent; COMSOL Stockholm 1986 multiphysics; OpenFOAM ESI Group open source CFD; Star-CCM Plus Siemens Digital Industries; Abaqus Dassault Systèmes Vélizy France; Aspen Plus chemical engineering simulation Cambridge MA 1981.

Tabella valori Pr

• Aria 0°C Pr: 0.713 standard atmosfera
• Aria 25°C Pr: 0.711 ambiente
• Aria 100°C Pr: 0.703
• Aria 500°C Pr: 0.683
• Aria 1000°C Pr: 0.731
• Acqua 0°C Pr: 13.6 alto
• Acqua 25°C Pr: 6.13
• Acqua 50°C Pr: 3.55
• Acqua 100°C Pr: 1.75
• Acqua 300°C Pr: 0.86
• Mercurio 25°C Pr: 0.0248 metallo liquido
• Sodio 100°C Pr: 0.011 nuclear coolant
• Piombo 400°C Pr: 0.024
• Olio motore 0°C Pr: 47000 estremo
• Olio motore 100°C Pr: 1500
• Olio idraulico 25°C Pr: 100
• Glicerolo 25°C Pr: 12500
• Etanolo 25°C Pr: 16
• Benzene 25°C Pr: 7
• Elio 25°C Pr: 0.68
• Idrogeno 25°C Pr: 0.71
• CO2 25°C Pr: 0.77
• R134a refrigerante Pr: 4-5 vapore
• Pr=nu/alpha: rapporto diffusività
• Pr=mu*cp/k: forma equivalente
• Reynolds analogy Pr=1: ideal heat-momentum
• Dittus-Boelter Nu: 0.023*Re^0.8*Pr^n
• Sieder-Tate: viscosity correction (mu/muw)^0.14
• Prandtl 1904 paper: Heidelberg congress 8 pagine
• Schmidt 1925: nomenclatura Prandtl-Zahl
• Italia Politecnico Milano 1863: DEIB DEAL
• Italia Politecnico Torino 1859: DET energia
• Italia Larderello 1904: prima centrale geotermoelettrica
• Italia ENEA Frascati 1982: energia ambiente
• Italia Leonardo 1948: 50000 employees aerospace
• Italia Avio Colleferro: Vega rocket
• Italia Ferrari Maranello 1947: F1 luxury
• Italia FCA Stellantis 2021: 300000 worldwide

Esempi numerici dettagliati

Esempio 1 - Aria 25°C ambiente. mu=1.81e-5 Pa·s, cp=1005 J/kg·K, k=0.0257 W/m·K. Pr = (1.81e-5 * 1005) / 0.0257 = 0.708. Vicino unità: heat e momentum diffondono ugualmente. Strato limite termico ~ cinematico. Italia: HVAC progettazione FCA Mirafiori Torino + Maserati Modena 1914 wind tunnel galleria vento test aerodinamici.

Esempio 2 - Acqua 25°C scambiatore calore. mu=8.9e-4 Pa·s, cp=4182 J/kg·K, k=0.598 W/m·K. Pr = (8.9e-4 * 4182) / 0.598 = 6.22. Pr>1: momento diffonde più rapido di calore. Strato limite cinematico più spesso. Importante design heat exchangers shell-and-tube. Italia: Politecnico Milano DEIB; Riello Verona caldaie 1922 6500 employees; Fondital Brescia caldaie.

Esempio 3 - Mercurio 25°C metallo liquido. mu=1.526e-3 Pa·s, cp=139.5 J/kg·K, k=8.34 W/m·K. Pr = (1.526e-3 * 139.5) / 8.34 = 0.0255. Pr<<1: calore diffonde molto più rapido di momento. Metalli liquidi nuclear reactor coolant fast breeder reactors. Italia ricerca nucleare: ENEA Frascati DTT Divertor Tokamak Test; INFN Frascati 1951.

Esempio 4 - Olio motore SAE 30 100°C. mu=0.486 Pa·s (alta viscosità), cp=1909 J/kg·K, k=0.144 W/m·K. Pr = (0.486 * 1909) / 0.144 = 6443. Pr>>1: momento domina. Strato limite cinematico molto spesso. Heat transfer lento diffusione termica in olio motore. Italia automotive: FCA Stellantis Torino Mirafiori; Ferrari Maranello 1947; Ducati Bologna 1926 motori.

Esempio 5 - Elio He gas leggero. mu=1.99e-5 Pa·s, cp=5193 J/kg·K, k=0.152 W/m·K. Pr = (1.99e-5 * 5193) / 0.152 = 0.680. Pr<1 gas tipico. He cryogenic cooling MRI superconductors. Italia: SAES Getters Lainate 1940 1500 employees criogenica He ultra-puro; INFN LNF Frascati cryogenic superconducting accelerator.

Limiti e considerazioni

1. Proprietà fluido isotermo: Pr = mu*cp/k a singola T, varia significativamente con temperatura (acqua 13.6 a 0°C, 1.75 a 100°C).
2. Fluidi newtoniani: Pr definito assume mu costante con shear rate. Non-newtoniani (polimeri, sangue) viscosità apparente varia.
3. Pressione effetti gas: Pr gas relativamente costante 0.7 ampio range pressione/temperatura. Liquidi p più rilevante.
4. Solo single phase: Pr applicabile fluidi monofase. Cambiamento fase (boiling, condensation) Nusselt correlazioni separate.
5. Compressibilità: gas ad alta velocità (Mach>0.3) Pr subisce correzione compressibilità Prandtl-Glauert.
6. Turbolenza vs laminare: Pr proprietà fluido. Turbolent Pr_t (eddy thermal/eddy momentum) ~ 0.85-0.95 vario.
7. Multicomponent: miscele liquidi/gas Pr media pesata + interaction effects. Modelli Wilke per gas.
8. Temperature estreme: criogenia Pr He liquido Pr~0.5 a 4K; alte T plasmi Pr varia notevolmente.
9. Reattività chimica: combustione reazioni Pr varia con composizione changing reactant→product.
10. Soluzioni concentrate: glicole etilenico antifreeze Pr varia 50-200 con concentrazione vs acqua pura.

Errori comuni

Confondere mu (dinamica) e nu (cinematica): mu Pa·s = kg/m·s viscosità dinamica; nu = mu/rho m²/s viscosità cinematica. Pr = mu*cp/k usa dinamica; Pr = nu/alpha usa cinematica equivalente.

Calore specifico cv vs cp: cp pressione costante (J/kg·K) corretto per Pr; cv volume costante. Per gas ideali cp = cv + R/M. Liquids cp ≈ cv (incompressibili).

Conducibilità k vs diffusività alpha: k W/m·K conduzione termica; alpha = k/(rho*cp) m²/s diffusività termica. Confusione frequente. Pr = nu/alpha relazione più diretta.

Pr a temperatura specifica: Pr varia molto con T per liquidi. Sempre indicare T riferimento (es. acqua 25°C Pr=6, non specificare implica ambiente).

Unità SI vs CGS: SI Pa·s (kg/m·s) vs cP centipoise (CGS) = 10^-3 Pa·s. Conversione importante: acqua 25°C mu = 0.89 cP = 8.9e-4 Pa·s.

Applicazioni pratiche

1. Heat exchangers shell-and-tube: Dittus-Boelter Nu=0.023 Re^0.8 Pr^n. Italia ricerca CNR ITC Bologna; Politecnico Milano DEAL.
2. HVAC progettazione: aria Pr=0.71 standard climatizzazione; FCA Stellantis Torino + Riello Verona caldaie.
3. Aerodinamica veicoli: galleria vento Politecnico Torino DET + Mirafiori; Ferrari Maranello + Lamborghini Sant'Agata + Ducati Bologna.
4. Aerospace: leading edge heating reentry; Avio Colleferro Vega rocket; Leonardo Roma 50000 employees.
5. Petrolchimica raffinerie: cracking heat transfer Pr oils. Eni Milano + Saipem; ENEA Roma.
6. Industria alimentare: sterilizzazione UHT pastorizzazione Italia Parma food valley; Tetra Pak Modena 1965; SIG Combibloc.
7. Nuclear reactor coolant: liquid metals Hg Na Pr<<1; ENEA Frascati DTT Divertor Tokamak; INFN Frascati 1951.
8. Fluidi geotermici: Larderello Pisa 1904 prima geotermoelettrica mondo Piero Ginori Conti; ENEL geotermia Toscana 800 MW.
9. Cryogenic engineering: He liquido Pr 0.5 superconductors MRI; SAES Getters Lainate 1940; Politecnico Milano.
10. Boundary layer aerospace: Prandtl 1904 Göttingen + Italia Leonardo F-35 partner; Piaggio Aerospace Genova 1884.

Aggiornato a luglio 2026