Come calcolare il numero di Biot (Bi): numero adimensionale fondamentale heat transfer transitorio (transient conduction) che misura il rapporto tra resistenza convezione esterna superficie solido e resistenza conduzione interna del solido stesso. Formula: Bi = h * Lc / k_solid dove h coefficiente di convezione tra fluido e superficie solido (W/m²·K), Lc lunghezza caratteristica del solido (m, tipicamente V/A_surface volume/area superficie), k_solid conducibilità termica del solido (W/m·K) NB: solido non fluido. Bi rappresenta come la temperatura del solido evolve durante riscaldamento/raffreddamento transitorio.

Criterio lumped capacitance method (sistema concentrato): Bi < 0.1 temperatura del solido uniforme spazialmente, può essere trattata come singola temperatura T(t) variabile nel tempo. Soluzione lumped: T(t) - T_inf = (T_0 - T_inf) * exp(-t/tau), dove tau = rho*c_p*V / (h*A) constante tempo (s). Errore <5% Bi<0.1, <2% Bi<0.01.

Bi > 0.1: gradienti significativi temperatura nel solido, soluzione richiede analisi spaziale completa: equazione del calore Fourier con condizioni iniziali e contorno. Soluzione analitica serie Fourier o Heisler charts grafici (Heisler 1947) o approssimazione one-term.

Storia: Jean-Baptiste Biot (Francia 1774-1862) Sorbonne Paris 1804 paper Bibliothèque Britannique 27:310 conduction calore solido prima trattazione matematica della propagazione del calore in solidi. Joseph Fourier (Francia 1768-1830) 1822 paper 'Théorie analytique de la chaleur' equazione del calore Fourier law q = -k*grad(T). Wilhelm Nusselt 1915 dimensional analysis Buckingham. Heisler charts 1947 Massachusetts Institute Technology MIT solutions transient heat conduction shapes simple slab cylinder sphere.

Strumento essenziale per: raffreddamento componenti elettronici CPU GPU heat sinks; industria alimentare sterilizzazione UHT pastorizzazione cottura time-temperature treatment Italia Parma food valley; metallurgia trattamenti termici tempra acciai quenching cooling; industria automotive Italia FCA Stellantis Ferrari brake disc cooling; aerospace Italia Avio Colleferro Vega rocket Leonardo Roma 50000 employees re-entry shock heating; edilizia bioedilizia Italia thermal mass walls; meteorologia raffreddamento corpo umano hypothermia survival.

Italia ricerca termofluidodinamica: Politecnico Milano DEIB Bioengineering 1863 fondato 4500 employees + DEAL Aerospaziale + DICEAM Industriale + DICEAM Civile; Politecnico Torino DET Energia 1859 fondato Mirafiori automotive research; Sapienza Roma DIETI Ingegneria Industriale 1303 fondata; Federico II Napoli DIETI 1224 fondata; Padova DEI 1222 Galileo professor; Bologna alma mater 1088 Faculty Engineering DIN; Pisa SNS Scuola Normale Superiore 1810; Pavia 1361; Genova DIME 1471. Italia istituzioni ricerca: CNR Roma 1923 8000 ricercatori; ENEA Roma + Frascati + Bologna 1982 energia ambiente fusione; ASI Roma 1988; INFN Frascati 1951; Avio Colleferro Vega rocket; Leonardo Roma 1948 50000 employees aerospace defense F-35 partner; Maire Tecnimont Roma 2005 6500 engineering oil&gas chemicals. Italia industria alimentare termica: Tetra Pak Italia Modena 1965 + Lund Sweden parent packaging UHT sterilization; Barilla Parma 1877 8500 employees pasta produzione termica; Galbani Milano 1882 caseificio termici; Parmalat Parma 1961 latte UHT pastorizzazione;

Ferrero Alba 1946 30000 employees cioccolato Nutella temper chocolate; Lavazza Torino 1895 4000 caffè torrefazione termici; illy Trieste 1933 caffè; Mutti Parma 1899 conserve pomodoro sterilizzazione; Bauli Castel d'Azzano Verona 1922 dolciario lievitazione termiche. Italia metallurgia termiche tempra: Tenova Castellanza Varese 1990 impianti metallurgici metallurgical equipment 2300 employees; Danieli Buttrio Udine 1914 9000 employees impianti acciaio worldwide largest; Sider Energy Brescia 1958; SMS Group Mönchengladbach Germania subsidiary Italia; Marcegaglia Mantova 1959 acciaio profilati; Riva Group Genova 1954; ArcelorMittal Italia Taranto 2017 8 milioni tonnellate/anno largest steel plant Europe.

Formula numero Biot

Definizione fondamentale: Bi = h * Lc / k_solid. Numero adimensionale (no unità). h coefficiente convezione esterno (W/m²·K), Lc lunghezza caratteristica del solido (V/A_surface volume su area superficie), k_solid conducibilità termica del solido (W/m·K). NB: critico distinguere k_solido vs k_fluido. Nu = h*L/k_fluido convezione; Bi = h*Lc/k_solido conduzione interna.

Significato fisico: Bi rapporto resistenza convezione esterna (1/h) su resistenza conduzione interna (Lc/k_solid). Bi<>1: convezione veloce, conduzione lenta, gradienti spaziali importanti. Lumped capacitance method approssimazione utilizzata Bi<0.1 errore <5%.

Lumped capacitance method: Bi<0.1 valido. Soluzione T(t) = T_inf + (T_0 - T_inf) * exp(-t/tau) dove tau = rho * c_p * V / (h * A) = rho * c_p * Lc / h time constant. Tempo per raffreddare/riscaldare al 63.2% (1-1/e) iniziale-finale. 5*tau ~ 99% completamento.

Heisler charts (1947): Bi>0.1 usano grafici trial dimensionless time tau* = alpha*t/Lc² (Fourier number Fo) e dimensionless temperature theta = (T(t) - T_inf)/(T_0 - T_inf). Three classic shapes: slab piano, cilindro, sfera. One-term approximation: theta_0 = C * exp(-zeta_1² * Fo) per Fo>0.2 (zeta_1 root function tabulato).

Lunghezza caratteristica: per slab piano spessore 2L: Lc = L (mezzo spessore); cilindro infinito raggio R: Lc = R/2; sfera raggio R: Lc = R/3; geometria irregolare Lc = V/A_surface. Distinguere Lc transient conduction (V/A) vs L convection correlations (geometric).

Storia e contesto

Jean-Baptiste Biot (Francia 1774-1862) Paris born, École Polytechnique 1794 mathematics 1797 graduate primo classe, Sorbonne professor mathematics 1800-1809 + Collège de France professor 1800-1862 long carrière 62 anni, Académie des Sciences 1803 elected. Biot polymath: astronomy (Biot crater Moon), geomagnetism (Biot-Savart law magnetism Felix Savart 1820), polarization light 1812 (Biot-Fresnel law optics), spectroscopy, conductione del calore 1804 paper Bibliothèque Britannique 27:310 'Mémoire sur la propagation de la chaleur' first mathematical treatment heat propagation in solids preceded Fourier 1822. Joseph Fourier (Francia 1768-1830) Auxerre Paris École Polytechnique 1794 mathematics + Egypt expedition 1798 Napoleon scientific commission 1801-1802 prefect Isère Grenoble, 1822 magnum opus 'Théorie analytique de la chaleur' Paris classica heat equation Fourier law q = -k*grad(T) + Fourier series + Fourier transform. Académie Sciences Paris secrétaire perpétuel 1822-1830. M. P. Heisler (USA) 1947 paper Transactions ASME 69:227 'Temperature Charts for Induction and Constant Temperature Heating' grafici universali dimensionless Heisler charts trans heat conduction shapes simple. 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;

Yunus Cengel (Turkey-USA 1953-) UNR; Theodore Bergman (USA) John Lavine 7th-8th edition Incropera. Italia ricerca termofluidodinamica: Politecnico Milano DEIB Bioengineering 1863 fondato + DEAL Aerospaziale + DICEAM 4500 employees Faculty Engineering Mechanical; Politecnico Torino DET Dipartimento Energia 1859 fondato + DENERG + DIATI Mirafiori automotive research center FIAT FCA Stellantis collaboration; Sapienza Roma DIETI Ingegneria Industriale 1303 fondata oldest Roma; Federico II Napoli DIETI 1224 Federico II Hohenstaufen emperor; Padova DEI 1222 fondata Galileo professor 1592-1610; Bologna alma mater 1088 prima Europa Faculty Engineering DIN Industrial Engineering; Pisa SNS Scuola Normale Superiore 1810 Mathematics. CNR Consiglio Nazionale Ricerche Roma 1923 fondato 8000 ricercatori 200 institutes ITC Costruzioni Bologna; ENEA Roma + Frascati + Bologna 1982 energia ambiente fusione DTT Divertor Tokamak Test Frascati; ASI Agenzia Spaziale Italiana Roma 1988; INFN Frascati 1951 5000 ricercatori; CINECA Bologna 1969 supercomputing CFD heat transfer simulations; Avio Colleferro propulsione spaziale Vega rocket Solid Rocket Motor Z23 Z40 Z9; Leonardo Roma 1948 Defense Aerospace 50000 employees worldwide F-35 partner; Piaggio Aerospace Genova 1884; Maire Tecnimont Roma 2005 6500 employees engineering oil&gas chemicals worldwide.

Italia industria alimentare termica processi: Tetra Pak Italia Modena 1965 (subsidiary Tetra Pak International AB Lund Sweden 1951 founded Ruben Rausing) packaging UHT sterilization aseptic process; Barilla Parma 1877 fondata Pietro Barilla 8500 employees worldwide pasta produzione termica essiccazione; Galbani Milano 1882 fondata Galbani Egidio caseificio termici stagionatura formaggi Reggio; Parmalat Parma 1961 fondata Calisto Tanzi 16000 worldwide latte UHT pastorizzazione 60-100°C; Ferrero Alba 1946 fondata Pietro Ferrero 30000 worldwide cioccolato Nutella tempering chocolate 27-31°C; Lavazza Torino 1895 fondata Luigi Lavazza 4000 caffè torrefazione termica 200-220°C; illy Trieste 1933 fondata Francesco Illy caffè; Mutti Parma 1899 conserve pomodoro sterilizzazione 90-100°C; Bauli Castel d'Azzano Verona 1922 dolciario; Granarolo Bologna 1957 latte; Conserve Italia Bologna 2002 fusion conserve. Italia metallurgia trattamenti termici: Tenova Castellanza Varese 1990 fondata 2300 employees impianti metallurgici metallurgical equipment 60+ countries; Danieli Buttrio Udine 1914 fondata 9000 employees impianti acciaio worldwide largest equipment supplier; Sider Energy Brescia 1958 forni industriali; Marcegaglia Mantova 1959 acciaio profilati; Riva Group Genova 1954 acciaio worldwide; ArcelorMittal Italia Taranto 2017 ex-ILVA 8 milioni tonnellate/anno crude steel Europe largest plant; Acciaierie Beltrame Vicenza 1896; Pittini Group Udine 1960; Feralpi Lonato Brescia 1968;

Cogne Acciai Speciali Aosta 1916; Acciaierie Valbruna Vicenza 1925.

Tabella valori Biot

• Bi = h*Lc/k_solid: definizione
• Bi < 0.01: very low, lumped highly accurate <2% error
• Bi < 0.1: low, lumped capacitance valid <5%
• Bi 0.1-1: moderate, Heisler charts
• Bi 1-10: high, conduction-dominated
• Bi > 100: extreme, T_surface = T_fluid
• Lumped tau: rho*cp*V/(h*A)
• Time constant: 5*tau ~ 99% completamento
• Slab Lc: spessore mezzo L
• Cilindro Lc: R/2
• Sfera Lc: R/3
• Acciaio k: 50 W/m·K
• Rame k: 400 W/m·K
• Alluminio k: 240 W/m·K
• Mattone k: 0.7 W/m·K
• Cibo k: 0.4-0.6 W/m·K
• Acqua liquida k: 0.6 W/m·K
• Aria k: 0.026 W/m·K
• h aria libera: 5-25 W/m²K
• h aria forced: 25-250 W/m²K
• h acqua libera: 100-1000 W/m²K
• h acqua forced: 500-15000 W/m²K
• h boiling: 2500-50000 W/m²K
• h condensation: 5000-100000 W/m²K
• Biot 1804 paper: Bibliothèque Britannique 27:310
• Fourier 1822: Théorie analytique chaleur
• Heisler 1947: ASME 69:227 charts
• Italia Politecnico Milano 1863: DEIB DEAL
• Italia Politecnico Torino 1859: DET energia
• Italia Tetra Pak Modena 1965: UHT packaging
• Italia Parmalat Parma 1961: 16000 latte UHT
• Italia Barilla Parma 1877: 8500 pasta
• Italia Ferrero Alba 1946: 30000 Nutella
• Italia Tenova Varese 1990: 2300 metallurgical
• Italia Danieli Udine 1914: 9000 acciaio worldwide
• Italia ArcelorMittal Taranto 2017: 8M tonnellate/anno

Esempi numerici dettagliati

Esempio 1 - Sfera acciaio raffreddamento aria. Sfera acciaio R=0.05 m raffreddata aria. h=25 W/m²K convezione naturale, k_steel=50 W/m·K. Lc = R/3 = 0.017 m. Bi = 25*0.017/50 = 0.0083 < 0.01 very low. Lumped highly accurate <2% error. Italia: trattamenti termici Tenova Varese 1990 + Danieli Udine 1914.

Esempio 2 - Mattone forno acqua. Mattone Lc=0.1 m immerso acqua bollente convezione forced h=500 W/m²K, k_brick=0.7 W/m·K. Bi = 500*0.1/0.7 = 71.4 high. Gradienti grandi nel solido, T_surface ~ T_water 100°C, interno raffreddato lentamente. Heisler charts richieste. Italia: industria laterizi mattoni; calce calcestruzzo CTI Comitato Termotecnico Italiano Milano 1933.

Esempio 3 - Componente rame olio. Componente rame Lc=0.01 m raffreddato olio h=200 W/m²K, k_copper=400 W/m·K. Bi = 200*0.01/400 = 0.005 very low. Lumped <2% error. Time const tau = rho*cp*Lc/h. Italia: cooling componenti elettronici Magneti Marelli Corbetta + Italtel Milano 1921 telecomunicazioni.

Esempio 4 - Cibo sterilizzazione vapore. Conserva alimentare cilindrica Lc=0.02 m vapore saturo h=5000 W/m²K, k_food=0.5 W/m·K. Bi = 5000*0.02/0.5 = 200 extreme high. T_surface = T_steam 121°C subito; interno raffreddato lento serie Fourier. Italia: Mutti Parma 1899 conserve pomodoro sterilizzazione; Tetra Pak Modena 1965; Parmalat Parma 1961.

Esempio 5 - Plate alluminio acqua time constant. Lastra alluminio Lc=0.005 m raffreddata acqua h=100, k=200, rho*cp=3.5e6 J/m³K. Bi = 100*0.005/200 = 0.0025 very low. tau = 3.5e6*0.005/100 = 175 s = 2.9 min costante tempo. 5*tau = 875 s = 14.6 min completamento 99%. Italia: heat sinks elettronica STMicroelectronics Catania Agrate; Magneti Marelli Corbetta Milano 1919.

Limiti e considerazioni

1. Lumped capacitance Bi<0.1: applicabile solo Bi<0.1 (errore <5%); oltre richiede Heisler charts o serie Fourier.
2. Geometria semplice: Heisler charts solo per slab piano, cilindro infinito, sfera. Forme irregolari richiedono FEA finite element.
3. Single-phase solido: applicabile solido omogeneo isotropo. Composti, multi-fase richiedono trattazione separata.
4. h costante: assume h coefficiente convezione costante nel tempo e spazio. Variazioni T richiedono iterativo.
5. Phase change esclusa: melting solidification freezing/thaw richiedono Stefan number + latent heat treatment.
6. Anisotropia: solidi anisotropi (legno, compositi) richiedono k tensoriale diverso x-y-z directions.
7. Convezione naturale variabile: in geometrie complesse h varia con posizione (Rayleigh number locali).
8. Radiation effects: ad alta temperatura T>500°C radiation contribution h_rad importante combinarlo con h_conv.
9. Steady state vs transient: Bi applicabile transient only. Steady state usa thermal resistance R = L/(k*A).
10. Internal heat generation: se solido genera calore (Joule heating, reazioni nucleari) trattazione separata Q_gen + Bi.

Errori comuni

k_solid vs k_fluid: Bi usa k del solido (acciaio=50, rame=400, alluminio=240, mattone=0.7); Nu usa k del fluido (acqua=0.6, aria=0.026). Errore comune confondere.

Lunghezza caratteristica: per slab Lc = mezzo spessore L; cilindro Lc = R/2; sfera Lc = R/3; irregolare Lc = V/A. Errore usare R intero invece R/2 o R/3.

Lumped Bi<0.1: criterio lumped Bi<0.1 errore <5%; oltre Heisler charts. Errore applicare lumped Bi=0.5 dà errori 30-50%.

Tempo costante tau: tau = rho*cp*Lc/h = rho*cp*V/(h*A). Tempo per raggiungere 63.2% finale. 5*tau ~ 99%. Errore confondere tau con tempo totale.

T_inf vs T_initial: T(t) = T_inf + (T_0 - T_inf)*exp(-t/tau). T_inf temperatura ambiente fluido (asintotico); T_0 temperatura iniziale solido. Errore inverte i due.

Applicazioni pratiche

1. Cooling componenti elettronici: CPU GPU heat sinks Bi<0.1 lumped. Italia: Magneti Marelli Corbetta + STMicroelectronics Catania.
2. Industria alimentare: sterilizzazione UHT pastorizzazione cottura time-temperature; Italia Parma food valley Mutti 1899 + Parmalat 1961 + Barilla 1877 + Ferrero Alba 1946; Tetra Pak Modena 1965.
3. Metallurgia tempra: cooling rate quenching steels martensite; Italia Tenova Varese 1990 + Danieli Udine 1914 + Marcegaglia + ArcelorMittal Taranto.
4. Industria automotive brake: brake disc cooling F1 racing; Italia FCA Stellantis Torino + Ferrari Maranello + Lamborghini + Ducati Bologna.
5. Aerospace re-entry: heat shield ablative tile thermal protection systems; Italia Avio Colleferro Vega rocket; Leonardo Roma 50000 employees F-35.
6. Edilizia bioedilizia: thermal mass walls passive cooling buildings; Italia case naturali Mapei Milano 1937 11000 employees adhesives chemicals.
7. Body cooling hypothermia: Newton's law of cooling lumped capacitance human body Bi<0.1 immersione acqua fredda survival time.
8. Quench glass tempering: tempered glass quenching automatically Bi>0.1 series solutions; Italia Saint-Gobain Italia subsidiary; vetreria Murano Venezia.
9. Fuel rod nuclear reactor: UO2 ceramic fuel internal heat generation + cladding cooling; Italia ENEA Frascati DTT Divertor Tokamak; INFN Frascati 1951.
10. Heat treatment dental implants: titanium implants heat treatment biocompatibility; Italia ricerca biomedical Politecnico Milano DEIB Bioengineering.

Aggiornato a luglio 2026