Calcolo Scambiatori Di Calore A Fascio Tubiero Excel

Calcola le prestazioni termiche del tuo scambiatore a fascio tubiero con precisione professionale. Inserisci i parametri tecnici per ottenere risultati dettagliati e grafici interattivi.

Guida Completa al Calcolo degli Scambiatori di Calore a Fascio Tubiero con Excel

Gli scambiatori di calore a fascio tubiero (shell and tube) rappresentano la soluzione più diffusa nell’industria per il trasferimento di calore tra due fluidi. Questo articolo fornisce una guida tecnica dettagliata per il dimensionamento e la verifica di questi apparati utilizzando fogli di calcolo Excel, con particolare attenzione agli aspetti termodinamici, fluidodinamici e costruttivi.

Principi Fondamentali degli Scambiatori a Fascio Tubiero

Gli scambiatori shell and tube consistono in:

• Fascio tubiero: serie di tubi attraverso cui scorre uno dei fluidi (generalmente quello più pulito o a pressione più elevata)
• Mantello (shell): contenitore esterno dove scorre il secondo fluido
• Testate (header): elementi che distribuiscono/raccolgono il fluido nei tubi
• Deflettori (baffles): piastre che dirigono il fluido nel mantello per migliorare lo scambio termico

Il funzionamento si basa sulla differenza di temperatura tra i due fluidi (ΔT) che determina il flusso termico secondo l’equazione fondamentale:

Q = U × A × ΔT_ml

Dove:

• Q: potenza termica scambiata (W)
• U: coefficiente globale di scambio termico (W/m²K)
• A: area di scambio (m²)
• ΔT_ml: differenza media logaritmica di temperatura (K)

Metodologia di Calcolo Step-by-Step

1. Definizione dei parametri operativi
   • Portate massiche dei fluidi (kg/s)
   • Temperature di ingresso/uscita (°C)
   • Proprietà termofisiche (calore specifico, conducibilità, viscosità)
2. Calcolo del bilancio termico

   Verifica che: Q_caldo = ṁ_h × c_p,h × (T_h,in – T_h,out) ≈ Q_freddo = ṁ_c × c_p,c × (T_c,out – T_c,in)

3. Determinazione della ΔT_ml

   Per scambiatori in equicorrente/controcorrente:

   ΔT_ml = [(T_h,in – T_c,out) – (T_h,out – T_c,in)] / ln[(T_h,in – T_c,out)/(T_h,out – T_c,in)]

   Per configurazioni più complesse (2-4 passaggi) si applica un fattore di correzione F (da grafici TEMA o calcolo numerico).

4. Stima del coefficiente globale U

   1/U = 1/h_i + (d_o/d_i) × (1/h_o) + (d_o × ln(d_o/d_i))/(2k) + R_f,i + R_f,o

   Dove h_i e h_o sono i coefficienti convettivi interno/esterno, calcolabili con correlazioni empiriche (Dittus-Boelter per turbolento, Sieder-Tate per laminare).

5. Calcolo dell’area richiesta

   A = Q / (U × F × ΔT_ml)

6. Verifica delle perdite di carico

   Lato tubi: ΔP = f × (L/d) × (ρv²/2) × N_pass

   Lato mantello: ΔP = 8 × J_f × (D_s/d_o) × (L/B) × (ρv²/2)

   Dove J_f è il fattore di attrito per fascio tubiero (da grafici o equazioni specifiche).

Implementazione in Excel: Struttura del Foglio di Calcolo

Un foglio Excel ben strutturato per questi calcoli dovrebbe includere:

Sezione	Contenuto	Formule Chiave
Input	Parametri operativi (temperature, portate, proprietà fluidi)	=SEERC.Viscosità(T; “Acqua”)
Bilancio Termico	Calcolo Q e verifica coerenza	=m_dot_h*Cp_h*(T_h_in-T_h_out)
ΔTml	Calcolo differenza media logaritmica	=((T1-T4)-(T2-T3))/LOG((T1-T4)/(T2-T3))
Coefficienti Convettivi	hi e ho con correlazioni	=0.023*Re^0.8*Pr^n*(k/d)
Coefficiente Globale	Calcolo 1/U	=1/hi+do/(di*ho)+…+Rf
Dimensionamento	Area richiesta e geometria	=Q/(U*F*ΔTml)
Perdite di Carico	ΔP tubi e mantello	=f*(L/d)*(rho*v^2/2)*Npass
Output	Risultati finali e grafici	=SEERC.Grafico(Temperatura; Posizione)

Correlazioni Empiriche per il Calcolo dei Coefficienti Convettivi

La precisione del calcolo dipende dalle correlazioni utilizzate per determinare i coefficienti di scambio termico convettivo:

Regime	Correlazione	Campo di Applicazione	Precisione
Turbolento (Re > 10,000)	Dittus-Boelter: Nu = 0.023 × Re^0.8 × Pr^n	0.7 < Pr < 160 L/d > 60	±25%
Transizione (2,300 < Re < 10,000)	Gnielinski: Nu = (f/8)(Re-1000)Pr/(1+12.7(f/8)^0.5(Pr^2/3-1))	0.5 < Pr < 2000 3000 < Re < 5×10^6	±20%
Laminare (Re < 2,300)	Sieder-Tate: Nu = 1.86 × (Re × Pr × d/L)^1/3 × (μ/μw)^0.14	0.48 < Pr < 16,700 0.0044 < (μ/μw) < 9.75	±30%
Lato Mantello (Kern)	ho = 0.36 × k/de × Re^0.55 × Pr^1/3 × (μ/μw)^0.14	2000 < Re < 1,000,000	±35%

Per il calcolo del fattore di attrito nei tubi, la correlazione di Colebrook-White (implicita) può essere approssimata con:

1/√f = -2.0 × log[(ε/d)/3.7 + 2.51/(Re × √f)]

Dove ε è la rugosità assolute del tubo (per acciaio commercialmente liscio, ε ≈ 0.045 mm).

Esempio Pratico di Calcolo

Consideriamo un caso reale di scambiatore per riscaldare 10 kg/s di acqua da 20°C a 60°C utilizzando olio termico che si raffredda da 150°C a 120°C. I parametri geometrici sono:

• Diametro tubi: 25.4 mm (1″)
• Spessore: 2 mm
• Materiale: Acciaio inox (k = 16 W/mK)
• Passaggi: 1-2 (1 passaggio mantello, 2 passaggi tubi)
• Numero tubi: 120
• Lunghezza: 4 m
• Diametro mantello: 400 mm

Passo 1: Bilancio termico

Q = ṁ_acqua × c_p,acqua × (60-20) = 10 × 4.18 × 40 = 1,672 kW

Verifica lato olio: Q = 10 × 2.2 × (150-120) = 1,650 kW (differenza < 2% → accettabile)

Passo 2: ΔT_ml

ΔT₁ = 150 – 60 = 90°C

ΔT₂ = 120 – 20 = 100°C

ΔT_ml = (90 – 100)/ln(90/100) = 94.87°C

Fattore correzione F per configurazione 1-2: ~0.95 (da grafici TEMA)

Passo 3: Coefficienti convettivi

Lato acqua (tubi):

• Re = (4 × ṁ)/(π × d × μ) ≈ 22,000 (turbolento)
• Nu = 0.023 × 22,000^0.8 × 5.8^0.4 ≈ 120
• h_i = Nu × k/d ≈ 120 × 0.6 / 0.0214 ≈ 3,360 W/m²K

Lato olio (mantello):

• Re ≈ 8,500 (transizione)
• Nu ≈ 80 (correlazione Kern)
• h_o ≈ 80 × 0.12 / 0.0254 ≈ 378 W/m²K

Passo 4: Coefficiente globale

1/U = 1/3,360 + (0.0254/0.0214)/378 + (0.0254 × ln(25.4/21.4))/(2×16) + 0.0002 + 0.0002

U ≈ 300 W/m²K

Passo 5: Area richiesta

A = 1,672,000 / (300 × 0.95 × 94.87) ≈ 61.5 m²

Area disponibile: π × d × L × N_tubi = π × 0.0214 × 4 × 120 ≈ 32.3 m² → Sottodimensionato!

Soluzioni possibili:

• Aumentare il numero di tubi a 230 (A ≈ 61 m²)
• Allungare i tubi a 7.5 m
• Utilizzare tubi con alette (finned tubes)

Ottimizzazione del Design

L’ottimizzazione di uno scambiatore shell and tube richiede un bilanciamento tra:

1. Prestazioni termiche: massimizzare U × A
2. Perdite di carico: limitare ΔP < 50-100 kPa
3. Costi: minimizzare materiali e complessità costruttiva
4. Manutenibilità: accesso per pulizia e ispezione

Strategie di ottimizzazione:

• Variazione del passo dei tubi: passo triangolare (30°) aumenta U del 20-30% rispetto a passo quadrato
• Spaziatura dei deflettori: B/L ≈ 0.2-0.5 (dove B = spaziatura, L = diametro mantello)
• Configurazione passaggi: 1-2 offre il miglior compromesso tra efficienza e ΔP
• Materiali: rame per alta conducibilità (k=400 W/mK), titanio per resistenza alla corrosione
• Fouling: prevedere fattori di incrostazione realistici (0.0001-0.0005 m²K/W per acqua trattata)

Un analisi costi-benefici tipica mostra che:

• Aumentare la velocità dei fluidi del 20% → +15% su U ma +60% su ΔP
• Ridurre il diametro dei tubi del 20% → +25% su U ma +100% su ΔP
• Usare alette → +300% su A esterna ma +20% costo e +30% ΔP

Validazione dei Risultati

La validazione del modello Excel può essere effettuata attraverso:

1. Confronti con software commerciali:
   • HTRI Xchanger Suite (standard industriale)
   • Aspen Exchanger Design & Rating
   • COMSOL Multiphysics (CFD)
2. Dati sperimentali:
   • Misure in impianti pilota
   • Test su prototipi (norma TEMA Class R)
3. Correlazioni pubblicate:
   • Grafici TEMA per F e J_f
   • Tabelle di proprietà termofisiche (NIST REFPROP)

Una discrepanza < 10% tra il modello Excel e i software professionali è generalmente accettabile per applicazioni industriali non critiche. Per progetti ad alta precisione (es. scambiatori per centrali nucleari), sono richieste analisi CFD 3D.

Errori Comuni da Evitare

Nella pratica ingegneristica, gli errori più frequenti includono:

1. Sottostima del fouling:
   • Utilizzare sempre fattori realistici (es. 0.00035 m²K/W per acqua di torre)
   • Prevedere sovradimensionamento del 10-20% per applicazioni con fluidi sporchi
2. Scelta errata delle correlazioni:
   • Verificare sempre il regime di flusso (Re)
   • Per fluidi non-newtoniani, usare correlazioni specifiche (es. Metzner-Reed)
3. Trascurare le perdite di carico:
   • ΔP eccessive possono richiedere pompe più costose
   • In applicazioni con vapore, ΔP > 10% della pressione assoluta causa flash
4. Ignorare gli effetti termici secondari:
   • Variazione delle proprietà con la temperatura (es. μ(T))
   • Condensazione/surriscaldamento in cambiamenti di fase
5. Errori di unità di misura:
   • Assicurarsi che tutte le equazioni usino unità coerenti (SI o imperiali)
   • Attenzione a conversioni come 1 BTU/h = 0.293 W

Applicazioni Industriali Tipiche

Gli scambiatori shell and tube trovano applicazione in:

Settore	Applicazione Specifiche	Range Tipico	Materiali Comuni
Oil & Gas	Raffreddamento gas compressi Riscaldamento greggio	1-50 MW 300-600°C	Acciaio al carbonio Acciaio inox 316
Energia	Condensatori turbine Surriscaldatori	10-500 MW 50-300°C	Rame (tubi) Acciaio (mantello)
Chimica/Farmaceutica	Reattori a doppio mantello Cristallizzatori	0.1-10 MW -20 a 200°C	Acciaio inox 316L Hastelloy
Alimentare	Pasteurizzatori Concentratori	0.05-5 MW 60-150°C	Acciaio inox 304 Titano
HVAC	Chiller ad assorbimento Recupero calore	0.01-2 MW 5-90°C	Rame Alluminio

Per ciascuna applicazione, la normativa di riferimento include:

• TEMA (Tubular Exchanger Manufacturers Association) – Standard per design e fabbricazione
• ASME Sec. VIII – Requisiti per pressione
• API 660 – Scambiatori per industria petrolifera
• PED 2014/68/UE – Direttiva europea pressione

Strumenti Avanzati per la Progettazione

Oltre ai fogli Excel, gli ingegneri professionisti utilizzano:

1. Software dedicati:
   • HTRI Xist: analisi termica e fluidodinamica dettagliata
   • Aspen EDR: integrazione con simulazioni di processo
   • COMSOL: modelli CFD 3D per geometrie complesse
2. Librerie Python:

   import CoolProp.CoolProp as CP
   import numpy as np
   
   # Calcolo proprietà acqua a 80°C
   T = 80 + 273.15  # K
   props = {
       'cp': CP.PropsSI('C', 'T', T, 'P', 101325, 'Water'),
       'k': CP.PropsSI('L', 'T', T, 'P', 101325, 'Water'),
       'mu': CP.PropsSI('V', 'T', T, 'P', 101325, 'Water'),
       'rho': CP.PropsSI('D', 'T', T, 'P', 101325, 'Water')
   }
                   

3. Standard di calcolo:
   • Bell-Delaware method per lato mantello
   • Kern’s method per stima rapida
   • Donohue equation per ΔP lato mantello

L’integrazione tra questi strumenti permette di:

• Ottimizzare il design con algoritmi genetici
• Valutare l’impatto di variazioni operative (es. riduzione portata del 10%)
• Generare documentazione automatica per la certificazione

Fonti Autorevoli:

U.S. Department of Energy – Heat Exchanger Design Handbook MIT – Heat Exchanger Fundamentals NIST REFPROP – Thermophysical Properties Database