Come calcolare il numero di Grashof (Gr): numero adimensionale fondamentale convezione naturale (free convection) che misura il rapporto tra forze di galleggiamento (buoyancy, dovute alla differenza di densità per gradiente termico) e forze viscose (viscosity-induced friction). Formula: Gr = (g * β * ΔT * L³) / ν² dove g accelerazione gravitazionale (9.81 m/s²), β coefficiente di espansione termica volumetrica (1/K), ΔT differenza temperatura superficie-ambiente (K), L lunghezza caratteristica geometrica (m), ν viscosità cinematica del fluido (m²/s). Per gas ideali β = 1/T (T temperatura assoluta in Kelvin); aria 20°C β ~ 1/293 = 3.43e-3 K^-1; acqua varia con T, 20°C β = 2.07e-4 K^-1; alcoli organici 0.001-0.002 K^-1.

Numero Rayleigh Ra = Gr * Pr combina Grashof e Prandtl per analisi heat transfer convezione naturale. Critical Ra ~ 1e9 transizione laminare-turbolento parete verticale; Ra 1e9 turbulent. Storia: Franz Grashof (Germania 1826-1893) Karlsruhe Polytechnikum (oggi KIT) 1875 paper Theoretische Maschinenlehre similitude theory thermal convection introduzione gruppi adimensionali natural convection prima sistematica analisi. Edgar Buckingham (USA 1867-1940) NIST 1914 pi theorem dimensional analysis.

Strumento essenziale per: HVAC riscaldamento ventilazione radiatori convezione naturale Italia Riello Verona 1922 6500 + Aermec Verona 1961 + Galletti Bologna 1909; raffreddamento elettronica CPU GPU heat sinks senza ventola passive cooling; edilizia bioedilizia Italia thermal mass walls passive cooling case naturali; atmosfera meteorologia Italia ARPA + INGV Roma; geyser termici Italia Larderello Pisa Toscana 1904 prima centrale geotermoelettrica mondo; solar collector Italia Enel Green Power Roma + Archimede Solar Energy Massa Martana Perugia 2009; vetro fusione Italia industrie vetro Murano Venezia; industria alimentare sterilizzazione UHT Parma food valley.

Italia ricerca termofluidodinamica: Politecnico Milano DEIB Bioengineering 1863 + DEAL Aerospaziale; Politecnico Torino DET Energia 1859; Sapienza Roma DIETI 1303; Federico II Napoli DIETI 1224; Padova DEI 1222 Galileo professor; Bologna alma mater 1088 Faculty Engineering DIN; Pisa SNS Scuola Normale Superiore 1810; Pavia 1361. CNR Roma 1923 8000 ricercatori; ENEA Roma + Frascati + Bologna 1982; ASI Roma 1988; INFN Frascati 1951; CINECA Bologna 1969 supercomputing CFD natural convection simulations. Italia industria HVAC radiatori naturale convezione: Riello Verona 1922 6500 employees caldaie radiatori; Aermec Bevilacqua Verona 1961 climatizzazione + radiatori; Galletti Bologna 1909 ventilconvettori radiatori; Fondital Brescia 1971 caldaie radiatori; Beretta Vimercate Milano caldaie murali; Vaillant Italia subsidiary German parent; Bosch Termotecnica Italia subsidiary; Daikin Italia subsidiary Japan parent; Climaveneta Bassano Vicenza 1971 chiller air handlers; Lu-Ve Uboldo Varese 1985 condensatori.

Italia industria edilizia bioedilizia: Mapei Milano 1937 11000 chemicals adhesives; Italcementi Bergamo 1864 cement Heidelberg subsidiary; Buzzi Unicem Casale Monferrato 1907 cement; Knauf Italia subsidiary German parent gypsum; Saint-Gobain Italia subsidiary glass ceramics.

Italia geotermia: 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; Monte Amiata Toscana geotermia secondaria; Bagnore plants 30 MW.

Italia industria solare: Enel Green Power Roma 2008 5500 employees rinnovabili worldwide; Archimede Solar Energy Massa Martana Perugia 2009 first solar molten salt parabolic trough; Eni gas&luce Roma; Edison Energia Milano 1884.

Formula Grashof

Definizione fondamentale: Gr = (g * β * ΔT * L³) / ν². Numero adimensionale (units cancel: m/s² * 1/K * K * m³ / (m²/s)² = m^4/s² / m^4/s² = 1). g gravità 9.81 m/s² Earth; β coefficiente espansione termica volumetrica (1/K); ΔT |Ts - T∞| differenza temperatura assoluta superficie-ambiente; L lunghezza caratteristica (altezza parete verticale, diametro tubo orizzontale, lato cubo); ν viscosità cinematica = μ/ρ (m²/s).

Coefficiente espansione termica β: per gas ideali β = 1/T (T temperatura assoluta K), aria 20°C β = 1/293 = 3.41e-3 K^-1, aria 60°C β = 1/333 = 3.0e-3 K^-1; per liquidi tabulato: acqua 20°C β=2.07e-4, 60°C β=5.23e-4, 100°C β=7.82e-4 K^-1; alcohol etilico 1.10e-3; glicerolo 4.95e-4; mercurio 1.81e-4; oli vegetali 0.7e-3.

Numero Rayleigh Ra = Gr * Pr: combina Grashof + Prandtl Pr (rapporto diffusività momento/termica). Ra è il parametro principale heat transfer convezione naturale. Critical Ra: 1e9 transizione laminare-turbolento parete verticale; per layer fluido orizzontale (Bénard convection) Ra critico = 1708 onset cellule convettive Rayleigh-Bénard; Ra < 1708 conduzione pura.

Correlazioni Nusselt natural convection: Churchill-Chu vertical plate Nu = (0.825 + 0.387*Ra^(1/6) / (1 + (0.492/Pr)^(9/16))^(8/27))² valida 10^-1 < Ra < 10^12; orizzontale plate hot up Nu = 0.54*Ra^(1/4) (10^4-10^7), Nu = 0.15*Ra^(1/3) (10^7-10^11); cilindro orizzontale Churchill-Chu Nu = (0.6 + 0.387*Ra^(1/6) / (1+(0.559/Pr)^(9/16))^(8/27))²; sfera Nu = 2 + 0.589*Ra^(1/4)/(1+(0.469/Pr)^(9/16))^(4/9).

Numeri correlati: Reynolds Re = vL/ν forced convection (assente in natural). Prandtl Pr = ν/α diffusività. Nusselt Nu = h*L/k heat transfer coefficient. Bond Bo = ρgL²/σ gravity vs surface tension. Richardson Ri = Gr/Re² rapporto natural/forced convection. Ri >> 1 dominante natural; Ri << 1 dominante forced; Ri ~ 1 mixed convection.

Storia e contesto

Franz Grashof (Germania 1826-1893) Düsseldorf born, Karlsruhe Polytechnikum (TH Karlsruhe oggi KIT 1825 erstes Polytechnikum Welt) 1854-1893 long career, Polytechnischen Schule Karlsruhe rector 1872-1893. Grashof founder VDI Verein Deutscher Ingenieure 1856 Düsseldorf primo presidente. 1875 textbook 'Theoretische Maschinenlehre' Vol II Hydraulik introduzione gruppi adimensionali similitude theory natural convection. Grashof number nomenclatura attribuita posteriore a Wilhelm Nusselt 1915 dimensional analysis paper. Edgar Buckingham (USA 1867-1940) NIST Bureau Standards 1914 paper Physical Review 4:345 pi theorem dimensional analysis sistematizzò. Continuazione natural convection: Wilhelm Nusselt (Germania 1882-1957) Munich Karlsruhe Berlin 1915 paper Mitteilungen Forschungsarbeiten 38; Ernst Schmidt (Germania 1892-1975) Munich Schmidt-Zahl 1925 + heat transfer textbook 'Einführung in die technische Thermodynamik'. Lord Rayleigh (UK 1842-1919) Cambridge Cavendish Nobel Physics 1904 1916 paper Philosophical Magazine 32:529 'On Convection Currents in a Horizontal Layer of Fluid' Bénard convection layer fluido orizzontale critical Ra=1708 instability onset cellule convettive. Henri Bénard (France 1874-1939) Paris 1900 PhD thesis Sorbonne first observation Bénard cells Rayleigh-Bénard convection. Modern textbook references: Frank Incropera (USA 1939-) Notre Dame 'Fundamentals of Heat and Mass Transfer' 1981 textbook 7 editions Wiley classic;

Adrian Bejan (Romania-USA 1948-) Duke entropy generation minimization constructal theory 1996; Yunus Cengel (Turkey-USA 1953-) UNR Heat Transfer textbook; Stuart Churchill + Stuart Chu 1975 paper Industrial Engineering Chemistry Fundamentals 14:344 Churchill-Chu correlation vertical plate. Italia ricerca termofluidodinamica natural convection: Politecnico Milano DEIB Bioengineering 1863 fondato + DEAL Aerospaziale + DICEAM Industriale 4500 employees; Politecnico Torino DET Energia 1859 Mirafiori automotive research center; Sapienza Roma DIETI Ingegneria Industriale 1303 fondata oldest Roma; Federico II Napoli DIETI 1224 Federico II Hohenstaufen; Padova DEI 1222 fondata Galileo professor 1592-1610; Bologna alma mater 1088 prima Europa DIN. CNR Consiglio Nazionale Ricerche Roma 1923 fondato 8000 ricercatori 200 institutes ITC Costruzioni Bologna heat transfer; ENEA Roma + Frascati + Bologna 1982 energia ambiente fusione DTT Divertor Tokamak Test Frascati; ASI Agenzia Spaziale Italiana Roma 1988; INFN Frascati 1951; CINECA Bologna 1969 supercomputing HPC CFD natural convection simulations. Italia industria HVAC radiatori convezione: Riello Verona 1922 fondata Pilade Riello 6500 employees worldwide caldaie radiatori riscaldamento; Aermec Bevilacqua Verona 1961 fondata Costanzo Aldegheri 1300 employees climatizzazione + radiatori convezione naturale; Galletti Bologna 1909 fondata Carlo Galletti 1000 employees ventilconvettori radiatori;

Fondital Brescia 1971 caldaie radiatori 1500 employees; Beretta Vimercate Milano caldaie murali; Vaillant Italia Aprilia Roma subsidiary German parent Vaillant 17000 worldwide; Bosch Termotecnica Italia subsidiary Bosch Stuttgart 1886; Daikin Italia subsidiary Japan parent Daikin Industries 70000; Climaveneta Bassano Vicenza 1971 1500 chiller air handlers; Lu-Ve Uboldo Varese 1985 1100 condensatori; Roca Group Spagna subsidiary; Lasian Spagna; Vasco Belgio. Italia industria edilizia bioedilizia thermal mass: Mapei Milano 1937 fondata Rodolfo Squinzi 11000 worldwide chemicals adhesives building materials; Italcementi Bergamo 1864 fondata Cesare Pesenti cement Heidelberg subsidiary 2016 acquisition 11000; Buzzi Unicem Casale Monferrato 1907 fondata Federico Buzzi cement 9000 worldwide; Cementir Holding Roma 1947 cement; Knauf Italia subsidiary Iphofen Germania parent; Saint-Gobain Italia subsidiary Paris parent 1665 + Glass Vetri Murano Venezia. Italia geotermia storia: Larderello Pisa Toscana 1904 Piero Ginori Conti prima centrale geotermoelettrica mondiale 1904 5 lampadine accese ENEL acquired 1962 Larderello geotermica 800 MW capacity 30% energia geotermica mondiale Toscana centro mondo geotermia + Monte Amiata Toscana geotermia secondaria + Bagnore plants 30 MW; INGV Istituto Nazionale Geofisica Vulcanologia Roma 1999 Etna Vesuvio Stromboli + idrogeologia.

Italia industria solare: Enel Green Power Roma 2008 fondata 5500 employees rinnovabili worldwide solar wind hydro; Archimede Solar Energy Massa Martana Perugia 2009 first solar molten salt parabolic trough technology; Eni gas&luce Roma 2017 retail energy; Edison Energia Milano 1884; Iren Reggio Emilia utility multi-utility 9000; A2A Milano utility multi-servizio 13000. USA natural convection research: MIT Cambridge MA 1861 mechanical engineering; Stanford Mechanical Palo Alto CA 1925; Caltech Pasadena CA 1891; UC Berkeley CA 1868 Lawrence Berkeley Lab; Princeton NJ 1746; Cornell NY 1865; Carnegie Mellon Pittsburgh PA 1900; Georgia Tech Atlanta GA 1885; Purdue West Lafayette IN 1869; Notre Dame Indiana 1842 Frank Incropera textbook author. USA companies HVAC: Carrier Palm Beach Gardens FL 1915 53000 HVAC; Trane Davidson NC 1885 38000; Johnson Controls Milwaukee WI 1885 100000 + York 1874; Lennox International Richardson TX 1895; Rheem Atlanta GA 1925 water heaters; Goodman Houston TX 1975. USA software CFD natural convection: ANSYS Canonsburg PA 1970 Fluent finite element 5500; COMSOL Stockholm 1986 multiphysics; OpenFOAM open source; Star-CCM+ Siemens Plano TX.

Tabella valori Grashof

• Gr = (g*β*ΔT*L³)/ν²: definizione
• Ra = Gr*Pr: Rayleigh combinato
• Critical Ra ~ 1e9: laminar-turbulent vertical plate
• Critical Ra = 1708: Rayleigh-Bénard onset
• β = 1/T gas ideali (T in K)
• β acqua 20°C: 2.07e-4 K^-1
• β acqua 60°C: 5.23e-4 K^-1
• β alcohol etilico: 1.10e-3 K^-1
• β glicerolo: 4.95e-4 K^-1
• β aria 20°C: 3.41e-3 K^-1
• ν aria 25°C: 1.6e-5 m²/s
• ν acqua 25°C: 8.9e-7 m²/s
• g Earth: 9.81 m/s²
• g Mars: 3.71 m/s²
• g Moon: 1.62 m/s²
• Churchill-Chu vertical: Nu correlation
• Lord Rayleigh 1916 paper: Bénard convection
• Henri Bénard 1900: PhD Sorbonne
• Critical Ra Bénard: 1708 onset cells
• Grashof 1875 textbook: Theoretische Maschinenlehre
• Grashof Karlsruhe Polytechnikum: 1854-1893
• Grashof VDI 1856: founder + presidente
• Italia Larderello 1904: prima geotermoelettrica
• Italia ENEL Larderello: 800 MW geotermia
• Italia Riello Verona 1922: 6500 caldaie
• Italia Aermec Verona 1961: climatizzazione
• Italia Galletti Bologna 1909: ventilconvettori
• Italia Fondital Brescia 1971: 1500 caldaie
• Italia Climaveneta Vicenza 1971: chillers
• Italia Politecnico Milano DEIB 1863: 4500 + DEAL
• Italia Politecnico Torino DET 1859
• Italia Sapienza Roma DIETI 1303
• Italia Mapei Milano 1937: 11000 building
• Italia Italcementi Bergamo 1864: 11000
• Italia Enel Green Power Roma 2008: 5500 rinnovabili
• Italia Archimede Solar Energy 2009: solar parabolic

Esempi numerici dettagliati

Esempio 1 - Parete verticale aria 70°C in stanza 20°C. Aria, g=9.81, β=3.43e-3 K^-1 (1/293), ΔT=50 K, L=1 m altezza parete, ν=1.6e-5 m²/s. Gr = (9.81*3.43e-3*50*1)/(1.6e-5)² = 1.68/2.56e-10 = 6.57e9 turbulent. Ra = Gr*Pr = 6.57e9*0.71 = 4.66e9 turbulent vertical plate. Italia: Riello Verona 1922 6500 caldaie + radiatori; Aermec Verona 1961 climatizzazione.

Esempio 2 - Pentola acqua bollente. Acqua, g=9.81, β=5.23e-4 K^-1 (60°C average), ΔT=70 K (Ts=363K=90°C boiling, T∞=293K=20°C ambient), L=0.1 m diametro pentola, ν=1.0e-6. Gr = (9.81*5.23e-4*70*0.001)/(1e-6)² = 3.59e-4/1e-12 = 3.59e8 transizione. Ra = Gr*Pr = 3.59e8*6 = 2.15e9 turbulent boiling. Italia: industria food Italia Lavazza Torino 1895 + Tetra Pak Modena 1965.

Esempio 3 - Stanza riscaldata aria. Aria, g=9.81, β=3.43e-3 K^-1, ΔT=10 K (radiatore 30°C, ambiente 20°C), L=3 m altezza stanza, ν=1.6e-5. Gr = (9.81*3.43e-3*10*27)/(2.56e-10) = 9.085/2.56e-10 = 3.55e10 turbulent indoor convection. Italia: Galletti Bologna 1909 ventilconvettori; Fondital Brescia 1971; Riello Verona 1922.

Esempio 4 - Radiatore convezione naturale. Aria intorno radiatore 80°C in stanza 20°C, β=3.43e-3 K^-1, ΔT=60 K, L=0.5 m altezza radiatore, ν=1.6e-5. Gr = (9.81*3.43e-3*60*0.125)/2.56e-10 = 0.252/2.56e-10 = 9.86e8 transizione. Italia HVAC: Riello Verona 1922 6500 + Aermec Verona 1961 + Beretta Vimercate Milano + Bosch Termotecnica Italia subsidiary.

Esempio 5 - Microchip raffreddamento passivo. CPU 80°C aria ambiente 25°C, β=3.43e-3, ΔT=55 K, L=0.01 m heat sink fin height, ν=1.6e-5. Gr = (9.81*3.43e-3*55*1e-6)/2.56e-10 = 1.85e-6/2.56e-10 = 7.23e3 low laminar weak natural convection. Pr=0.71: Ra = 5.13e3 weak. Italia: STMicroelectronics Catania Agrate microfabbricazione MEMS; Magneti Marelli Corbetta 1919 41000 automotive electronics.

Limiti e considerazioni

1. Solo flussi naturali no-forced: Grashof per natural convection only. Forced + natural mixed: Richardson Ri = Gr/Re² > 1 natural domina.
2. Boussinesq approximation: assume β*ΔT < 50 K aria errori >5%, considerare densità variabile.
3. Newtoniani fluidi: assume ν costante. Non-newtoniani (polimeri, gel) richiedono trattazione separata.
4. Fluido omogeneo: assume single-phase fluid. Phase change boiling/condensation richiedono Stefan number + altri.
5. Gravità verticale: assume g costante verticale. Spazio (microgravità) g=0 no convection naturale.
6. Steady-state: applicabile flussi steady-state. Transienti richiedono trattazione tempo-dipendente.
7. Geometria specifica: correlazioni Nu specifiche geometria (parete verticale, cilindro orizzontale, sfera). Forme irregolari FEA.
8. Boundary conditions: T_wall costante vs q_wall costante danno correlazioni diverse.
9. Atmospheric stability: in atmosphere stratification stabile/instabile influisce convezione, ricerca atmospheric science specifica.
10. Surface tension effects: scala micro/nano-fluidica Bond number Bo competing.

Errori comuni

L³ vs L² o L: Gr dipende da L³ (L cubed). Errore comune usare L² o L. Triplicare L aumenta Gr 27x.

β unità: K^-1 (per Kelvin). Per gas ideali β = 1/T (T Kelvin). Errore usare °C: T(K) = T(°C) + 273.15.

ν cinematica vs μ dinamica: ν = μ/ρ m²/s (kinematic); μ Pa·s (dynamic). Grashof usa ν kinematic. Errore confondere.

ΔT positiva: ΔT = |Ts - T∞| valore assoluto. Ts > T∞ heated wall buoyancy upward; Ts < T∞ cooled wall buoyancy downward. Errore segno.

Critical Ra confusion: 1708 Bénard layer orizzontale; 1e9 transizione vertical plate. Errore confondere applicazioni.

Applicazioni pratiche

1. HVAC radiatori convezione naturale: Italia Riello Verona 1922 6500 + Aermec Verona 1961 + Galletti Bologna 1909 + Fondital Brescia 1971 + Beretta Vimercate.
2. Raffreddamento elettronica passivo: heat sinks senza ventola. Italia: STMicroelectronics Catania Agrate; Magneti Marelli Corbetta 1919.
3. Edilizia bioedilizia: thermal mass walls case naturali. Italia: Mapei Milano 1937 + Italcementi Bergamo 1864 + Buzzi Unicem 1907.
4. Atmosfera meteorologia: thunderstorms convection, cumulus clouds. Italia: ARPA + INGV Roma 1999.
5. Geotermia: Italia Larderello Pisa Toscana 1904 prima centrale geotermoelettrica mondo + Monte Amiata + Bagnore plants 30 MW; ENEL Larderello 800 MW.
6. Solar collector flat plate: Italia Enel Green Power Roma 2008 + Archimede Solar Energy Massa Martana Perugia 2009.
7. Vetro fusione: convezione naturale fluido. Italia: vetreria Murano Venezia historical glassmaking.
8. Industria alimentare: cottura forni convezione naturale. Italia: SACMI Imola 1919; SMEG Guastalla Reggio Emilia 1948 cucine elettrodomestici.
9. Reattori chimici: stirred tank reactors. Italia: Eni Milano 1953 32000; Saipem Milano 1957 32000.
10. Ricerca aerospazio microgravità: g=0 spazio no natural convection; ASI Roma 1988 + esperimenti ISS Italian astronauts Cristoforetti Parmitano Nespoli.

Aggiornato a luglio 2026