Calcolatore Perdita di Carico Tubazioni (Darcy-Weisbach)

Calcola la perdita di carico distribuita ΔH_f in una tubazione con la formula di Darcy-Weisbach. Fattore di attrito da Colebrook-White/Swamee-Jain, numero di Reynolds, regime laminare/turbolento. Gratis, in 5 lingue.

Parametri di calcolo
Questo calcolatore ti è stato utile?

Che cos'è la perdita di carico in una tubazione

Quando un fluido scorre in un tubo, l'attrito con la parete interna e la turbolenza del flusso convertono progressivamente energia meccanica (pressione + velocità + quota) in calore, che si disperde nell'ambiente. Questa dissipazione si manifesta come una caduta di pressione tra ingresso e uscita: la perdita di carico, indicata con ΔH_f (in metri di colonna del fluido) o ΔP_f (in pascal). È il parametro chiave per dimensionare tubazioni di acquedotti, impianti di riscaldamento, oleodotti, gasdotti, circuiti di raffreddamento e per scegliere la pompa che deve vincere queste perdite fornendo la prevalenza necessaria. Sottostimarla comporta portate inferiori al richiesto; sovrastimarla porta a pompe sovradimensionate con maggiore consumo energetico. La formula di riferimento è quella di Darcy-Weisbach (1857), universalmente valida per fluidi newtoniani incomprimibili in condotti a sezione circolare.

La formula di Darcy-Weisbach

La perdita di carico distribuita (lungo la lunghezza della tubazione) è espressa da:

ΔH_f = f · (L/D) · v² / (2g)  [m di colonna fluido]

oppure in termini di pressione: ΔP_f = ρ·g·ΔH_f = f·(L/D)·(ρv²/2) [Pa]. I termini sono: f fattore di attrito di Darcy (adimensionale, dipende da regime e rugosità), L lunghezza della tubazione (m), D diametro interno (m), v velocità media del fluido (m/s, v = Q/A con A = π·D²/4), g = 9,81 m/s² accelerazione di gravità. Il termine v²/(2g) è la quota cinetica equivalente della velocità (in metri di fluido), e il fattore f·L/D è il numero di quote cinetiche di energia dissipate. Esempio: acqua in tubo DN100 (D=100 mm), L=200 m, Q=10 L/s → v = 1,27 m/s, con f ≈ 0,020 → ΔH_f ≈ 3,3 m di colonna d'acqua (32 kPa).

Il numero di Reynolds e i regimi di moto

Il numero di Reynolds è il parametro adimensionale che distingue il tipo di flusso:

Re = ρ·v·D/μ = v·D/ν

con μ viscosità dinamica (Pa·s) e ν = μ/ρ viscosità cinematica (m²/s). Per acqua a 20 °C: ν = 1,004·10⁻⁶ m²/s. Per olio motore SAE 30: ν ≈ 100·10⁻⁶ m²/s (100 volte più viscoso). I tre regimi:

  • Laminare (Re < 2300): il fluido scorre in strati paralleli, senza mescolamento. Profilo parabolico di velocità (max al centro, zero alle pareti). Fattore di attrito noto analiticamente: f = 64/Re (Poiseuille). Tipico per oli, tubi capillari, portate ridotte.
  • Transizione (2300 < Re < 4000): flusso instabile, alterna caratteristiche laminari e turbolente. Fattore di attrito indeterminato, si preferisce evitarlo nel dimensionamento.
  • Turbolento (Re > 4000): presenza di vortici, mescolamento tridimensionale, profilo di velocità appiattito. Fattore di attrito f dipende sia da Re sia dalla rugosità relativa ε/D. È il regime tipico di acquedotti, impianti idrici, gasdotti.

Esempio: acqua v = 1,27 m/s in tubo D = 100 mm → Re = 126 800 (turbolento). Olio v = 0,3 m/s in tubo D = 20 mm → Re = 64 (laminare).

Il fattore di attrito: da Colebrook a Swamee-Jain

In regime turbolento il fattore di attrito è dato dall'equazione implicita di Colebrook-White (1939):

1/√f = −2·log₁₀[ε/(3,7·D) + 2,51/(Re·√f)]

Questa equazione è trascendente (f compare sia a destra che a sinistra) e richiede iterazioni per essere risolta. Nel 1976 Swamee e Jain proposero un'approssimazione esplicita con precisione dell'1-2% rispetto a Colebrook:

f = 0,25 / [log₁₀(ε/(3,7·D) + 5,74/Re^0,9)]²

valida per 5·10³ < Re < 10⁸ e ε/D < 0,05. Il calcolatore usa questa formulazione. La rugosità assoluta ε (in mm o m) dipende dal materiale del tubo:

  • PVC nuovo: 0,007 mm.
  • Rame nuovo: 0,0015 mm.
  • Acciaio nuovo: 0,046 mm.
  • Acciaio zincato: 0,15 mm.
  • Acciaio con incrostazioni: 0,5-3 mm.
  • Ghisa nuova: 0,25 mm; ghisa vecchia: 1-3 mm.
  • Calcestruzzo: 0,3-3 mm.
  • Vetro: 0,001 mm (praticamente liscio).

La rugosità relativa ε/D è ciò che conta nel calcolo: un tubo grande con ε = 0,5 mm può essere "idraulicamente liscio" se D è grande (ε/D piccolo).

Il diagramma di Moody: il grafico universale

Il diagramma di Moody (Lewis Ferry Moody, 1944) è una rappresentazione grafica di f in funzione di Re, con curve parametrizzate per diversi valori di ε/D. È lo strumento tradizionale (pre-calcolatrici scientifiche) per leggere f. Sul diagramma si distinguono quattro zone:

  • Zona laminare (Re < 2300): retta di pendenza -1 in scala log-log, corrispondente a f = 64/Re. Le curve di rugosità coincidono qui perché il flusso laminare non "sente" la parete oltre lo strato limite viscoso.
  • Zona di transizione: gap tra Re 2300 e ~4000, tratteggiato per indicare l'indeterminazione.
  • Zona turbolento liscia: curve di rugosità piccolissima (ε/D < 10⁻⁵) tendono asintoticamente a quella di Prandtl-Kármán, dove il fattore di attrito dipende solo da Re: f = 0,316/Re^0,25 (Blasius, valida fino Re = 10⁵).
  • Zona pienamente turbolenta (Re grande e ε/D non trascurabile): le curve diventano orizzontali → f dipende solo da ε/D, non da Re. In questa zona f = [2·log₁₀(3,7·D/ε)]⁻² (Kármán-Nikuradse).

Il calcolatore automatizza la lettura del diagramma di Moody con la formula di Swamee-Jain.

Velocità consigliate e criteri di dimensionamento

Le velocità consigliate variano per applicazione (fonti UNI-EN 806-3 per acqua sanitaria, ASHRAE per HVAC):

  • Acqua potabile in acquedotto: 0,5-1,5 m/s (economia energetica); tolleranza max 2 m/s per brevi tratti.
  • Aspirazione di pompa: 0,5-1,5 m/s (evita cavitazione).
  • Mandata di pompa: 1,5-3 m/s.
  • Acqua calda sanitaria: 0,5-1,5 m/s (limitazione contro erosione e rumore).
  • Riscaldamento circolazione forzata: 0,5-1,2 m/s (silenziosità).
  • Vapore saturo: 15-40 m/s.
  • Gas naturale distribuzione: 5-20 m/s.
  • Aria in canali di ventilazione: 2-8 m/s (silenziosità).

Il criterio pratico "perdita di carico unitaria" è j = ΔH_f/L: valore raccomandato 0,005-0,03 m/m per acqua fredda in acquedotti (Bonicelli, Marchi-Rubatta). Sopra 0,05 m/m la pompa consuma troppo; sopra 0,15 m/m ci sono rischi di erosione e cavitazione. Un buon compromesso economico è la velocità ottima economica, dove la somma tra costo del tubo (grande → costoso) e costo pompaggio (piccolo → tanti kWh persi) è minima. Formula empirica: v_opt ≈ (Q/K)^0,2 con K = 3·10⁻⁴ m per acqua.

Perdite concentrate e prevalenza pompa

Oltre alle perdite distribuite (calcolate qui con Darcy-Weisbach) esistono le perdite concentrate in curve, raccordi, valvole, filtri, riduzioni. Si esprimono come:

ΔH_conc = ζ · v²/(2g)

con ζ coefficiente di perdita concentrata (adimensionale) tabellato per ogni pezzo: gomito 90° a raggio corto ζ ≈ 0,9; curva 90° a raggio largo ζ ≈ 0,3; valvola a sfera aperta ζ ≈ 0,05; valvola di regolazione parzialmente aperta ζ = 1-100. In un impianto residenziale la somma delle perdite concentrate può raggiungere il 30-50% delle distribuite. La prevalenza totale richiesta alla pompa è:

H_pompa = ΔH_geo + Σ ΔH_f + Σ ΔH_conc + ΔH_utenza

dove ΔH_geo è il dislivello geometrico da vincere, Σ ΔH_f la somma di tutte le perdite distribuite nei tratti, Σ ΔH_conc la somma delle perdite concentrate, ΔH_utenza la pressione residua richiesta al punto di consegna (tipicamente 10-15 m per rubinetti civili). La potenza elettrica della pompa è P_ele = ρ·g·Q·H_pompa/η_pompa, con η_pompa 60-85% per pompe centrifughe residenziali.

Come usare il calcolatore

Inserisci il diametro interno del tubo D in mm (non il diametro nominale DN, che è quello esterno: es. tubo Fe DN100 ha D_interno ≈ 105 mm). Inserisci la lunghezza L del tubo in metri (dalla sorgente al punto di consegna). Inserisci la portata Q in L/s (per convertire: 1 m³/h = 0,278 L/s; 60 L/min = 1 L/s). Inserisci la rugosità assoluta ε del tubo in mm (vedi tabella: rame 0,0015; PVC 0,007; acciaio nuovo 0,046; ghisa 0,25; acciaio con incrostazioni 0,5-3). Inserisci la densità del fluido ρ (kg/m³, acqua fredda 998; acqua calda 60°C 983; olio motore 850; glicole 40% 1060). Inserisci la viscosità cinematica ν (m²/s: acqua 20°C 1,004·10⁻⁶; acqua 60°C 0,467·10⁻⁶; olio motore 100·10⁻⁶; glicole 40% a 20°C 3,5·10⁻⁶). Il calcolatore restituisce la velocità v, il numero di Reynolds Re, il fattore di attrito f (regime laminare o turbolento con Swamee-Jain), la perdita di carico ΔH_f in metri di colonna, la caduta di pressione ΔP_f in kPa, la potenza dissipata in kW (utile per dimensionare la pompa) e una valutazione qualitativa.

Domande frequenti

Che differenza c'è tra il fattore di attrito di Darcy e quello di Fanning?

Il fattore di Fanning f_F è definito come metà di quello di Darcy: f_F = f_D/4. Le due formulazioni coesistono nella letteratura, generando spesso confusione. Il Fanning è preferito in ingegneria chimica e petrolifera USA; il Darcy è preferito in ingegneria civile e nella comunicazione europea/italiana. Il diagramma di Moody usa Darcy. Verifica sempre la definizione: per un flusso laminare f_D = 64/Re mentre f_F = 16/Re. Il calcolatore qui usa esclusivamente Darcy.

Perché l'equazione di Colebrook è implicita e serve iterare?

Perché il fattore di attrito f appare sia a destra (dentro il logaritmo) sia a sinistra (√f), quindi non si può isolare in forma chiusa. Storicamente si risolveva con: iterazione fixed-point (partendo da un f iniziale, si ricalcola fino a convergenza in 3-5 iterazioni), tramite abachi (nomogrammi tabellari), o grafico di Moody. Oggi si usano approssimazioni esplicite con precisione dell'1-3%: Swamee-Jain 1976 (quella usata qui), Haaland 1983 (formulazione alternativa più semplice), Chen 1979 (più precisa ma complessa). Per calcoli critici (ottimizzazione economica di grandi reti) si preferisce risolvere Colebrook per iterazione con precisione a 10⁻⁶.

Il calcolatore vale anche per gas (aria, metano)?

Sì, se il flusso è a bassa velocità (Mach < 0,3, cioè circa v < 100 m/s per aria a 20 °C), il gas può essere considerato incomprimibile e la formula Darcy-Weisbach si applica direttamente. Basta usare la densità ρ e la viscosità cinematica ν del gas alle condizioni di esercizio (per aria 20 °C, 1 atm: ρ = 1,20 kg/m³, ν = 1,5·10⁻⁵ m²/s). Per flussi ad alta velocità o comprimibili (aria in gasdotti a alta pressione con ΔP/P > 0,1) va usata la formula di Weymouth o Panhandle che tengono conto della variazione di densità lungo il tubo.

Come confrontare Darcy-Weisbach con formule empiriche (Hazen-Williams, Manning)?

Hazen-Williams è una formula empirica del 1902 per acqua a 20 °C: ΔH = 10,67·L·(Q/C)^1,852 / D^4,87 con C coefficiente di scabrezza (100-140 per tubi comuni). Facile da usare a mente/tabella ma valida solo per acqua a T ambiente e tubi non lisci (C > 100). Sottostima f a Re bassi (< 10⁵). Molto usata in USA per reti idriche. Manning è per canali a superficie libera. Darcy-Weisbach è universale (qualsiasi fluido, qualsiasi Re, qualsiasi rugosità) ed è la formula raccomandata dalle norme ISO e da IWA (International Water Association) per calcoli professionali. Le differenze pratiche sono ±5-15% tra Hazen-Williams e Darcy-Weisbach per acqua a portate normali.

Come tengo conto delle perdite localizzate (raccordi, valvole)?

Con il metodo delle lunghezze equivalenti o dei coefficienti di perdita concentrata ζ. Ogni raccordo, curva, valvola aggiunge una perdita ΔH_conc = ζ·v²/(2g). Valori tipici: gomito 90° raggio corto ζ = 0,9; gomito 90° raggio largo 0,3; T ramo dritto 0,3; T ramo con derivazione 0,9-1,5; valvola sfera aperta 0,05; valvola globo aperta 6-10; riduzione graduale 0,05-0,2; allargamento improvviso 1,0 (perdita di Borda-Carnot); filtro Y-strainer 2-3. In alternativa, si esprime ogni pezzo in metri equivalenti di tubo (L_eq = ζ·D/f): un gomito 90° in DN50 (D=50mm) con f=0,025 e ζ=0,9 equivale a L_eq = 0,9·0,050/0,025 = 1,8 m di tubo dritto. La somma dei L_eq va aggiunta alla lunghezza reale nel calcolatore.

Cosa succede se il flusso è nel regime di transizione (2300 < Re < 4000)?

Il regime di transizione è instabile: il flusso oscilla erraticamente tra laminare e turbolento con perdite di carico variabili nel tempo. È una zona da evitare nel dimensionamento perché non si può prevedere con precisione ΔH. Il calcolatore fa un'interpolazione lineare tra f_laminare a Re=2300 e f_turbolento a Re=4000 come approssimazione, ma il valore reale può oscillare del ±20-30%. Se il tuo calcolo cade in questa zona, meglio ridimensionare il tubo (diametro minore → v maggiore → Re maggiore → sicuramente turbolento) o accettare che la portata reale sarà instabile. Alcuni dispositivi (impianti pieno-vuoto, sistemi con dispositivi a controllo pulsato) generano naturalmente flussi di transizione, ma il dimensionamento va fatto con margine.

Come confrontare i risultati con software professionale (EPANET, WaterCAD)?

Software come EPANET (US EPA, gratuito, open source) e WaterCAD/Bentley (commerciale) risolvono reti idriche complesse con centinaia di nodi e tubi, applicando iterativamente le equazioni di Colebrook o Hazen-Williams a ogni tratto e bilanciando la conservazione di massa a ogni nodo. Il singolo tubo dà lo stesso risultato di questo calcolatore (±1% per Colebrook vs Swamee-Jain). Le differenze diventano significative solo se: (a) il fluido è comprimibile (gas ad alta pressione, che i software gestiscono con formule apposite); (b) c'è transitorio idraulico (colpo d'ariete, che richiede analisi non stazionaria con metodo delle caratteristiche); (c) rete magliata con più anelli, dove i software risolvono il sistema di equazioni non lineari a tutti i nodi contemporaneamente.

Formula

{"latex":"\Delta H_{f}=f\,\dfrac{L}{D}\,\dfrac{v^{2}}{2g};\;\;Re=\dfrac{v\,D}{\nu};\;\;f_{lam}=\dfrac{64}{Re};\;\;f_{turb}=\dfrac{0{,}25}{\left[\log_{10}\!\left(\tfrac{\varepsilon}{3{,}7 D}+\tfrac{5{,}74}{Re^{0{,}9}}\right)\right]^{2}}","variables":{"it":"ΔH_f perdita di carico distribuita (m di colonna fluido); f fattore di attrito di Darcy (adimensionale); L lunghezza del tubo (m); D diametro interno (m); v velocità media del fluido (m/s), v = Q/A con A = πD²/4; g = 9,81 m/s² accelerazione di gravità; Re numero di Reynolds; ν viscosità cinematica del fluido (m²/s); ε rugosità assoluta della parete interna (mm o m); ρ densità (kg/m³) usata per convertire ΔH in ΔP = ρgΔH. Regimi: laminare Re < 2300 (Poiseuille); turbolento Re > 4000 (Colebrook-White approssimata con Swamee-Jain 1976 esplicita). Transizione 2300-4000: instabile, si usa interpolazione lineare.","en":"ΔH_f distributed head loss (m of fluid column); f Darcy friction factor (dimensionless); L pipe length (m); D internal diameter (m); v mean fluid velocity (m/s), v = Q/A with A = πD²/4; g = 9.81 m/s² gravity; Re Reynolds number; ν fluid kinematic viscosity (m²/s); ε absolute wall roughness (mm or m); ρ density (kg/m³) used to convert ΔH into ΔP = ρgΔH. Regimes: laminar Re < 2300 (Poiseuille); turbulent Re > 4000 (Colebrook-White approximated by explicit Swamee-Jain 1976). Transition 2300-4000: unstable, linear interpolation.","es":"ΔH_f pérdida de carga distribuida (m columna fluido); f factor de fricción de Darcy; L longitud del tubo (m); D diámetro interno (m); v velocidad media (m/s), v = Q/A con A = πD²/4; g = 9,81 m/s²; Re número de Reynolds; ν viscosidad cinemática (m²/s); ε rugosidad absoluta (mm); ρ densidad (kg/m³) para ΔP = ρgΔH. Regímenes: laminar Re < 2300 (Poiseuille); turbulento Re > 4000 (Colebrook-White aproximada por Swamee-Jain 1976). Transición 2300-4000: interpolación lineal.","fr":"ΔH_f perte de charge répartie (m colonne fluide) ; f coefficient de Darcy ; L longueur du tuyau (m) ; D diamètre intérieur (m) ; v vitesse moyenne (m/s), v = Q/A avec A = πD²/4 ; g = 9,81 m/s² ; Re nombre de Reynolds ; ν viscosité cinématique (m²/s) ; ε rugosité absolue (mm) ; ρ densité (kg/m³) pour ΔP = ρgΔH. Régimes : laminaire Re < 2300 (Poiseuille) ; turbulent Re > 4000 (Colebrook-White approché par Swamee-Jain 1976 explicite). Transition 2300-4000 : interpolation linéaire.","de":"ΔH_f verteilter Druckhöhenverlust (m Fluidsäule); f Darcy-Reibungsbeiwert; L Rohrlänge (m); D Innendurchmesser (m); v mittlere Geschwindigkeit (m/s), v = Q/A mit A = πD²/4; g = 9,81 m/s²; Re Reynolds-Zahl; ν kinematische Viskosität (m²/s); ε absolute Wandrauhigkeit (mm); ρ Dichte (kg/m³) für ΔP = ρgΔH. Strömungsformen: laminar Re < 2300 (Poiseuille); turbulent Re > 4000 (Colebrook-White per Swamee-Jain 1976 approximiert). Übergang 2300-4000: lineare Interpolation."}}

Commenti (0)

Ancora nessun commento. Sii il primo!