Come calcolare la portata di un tubo Venturi: inventato da Giovanni Battista Venturi (Italia 1746-1822 Universita' Modena Reggio Emilia) nel 1797 sulla base dell'equazione di Daniel Bernoulli (Olanda-Svizzera 1700-1782) "Hydrodynamica" 1738. Strumento essenziale per la misura della portata di fluidi in condotte chiuse: acqua negli acquedotti urbani (Acea Roma, A2A Milano, Hera Bologna gestori italiani), gas naturale nei gasdotti (Snam Rete Gas San Donato Milanese), oleodotti petroliferi (Eni Saipem), HVAC ventilazione industriale (Carrier USA, Daikin Giappone), processo alimentare-farmaceutico (GEA Group Germania, Tetra Pak Lund Svezia), medicale ventilatori polmonari (Philips Healthcare Eindhoven, Hamilton Medical Bonaduz Svizzera), centrali idroelettriche (Enel Green Power, Edison Foro Buonaparte Milano). Formula fondamentale: Q = Cd * A2 * sqrt(2*dP/(rho*(1-(A2/A1)^2))), dove Q portata volumetrica (m^3/s), Cd coefficient of discharge (tipico 0,98 turbolento), A1 area sezione tubo (m^2), A2 area gola Venturi (m^2), dP differenza pressione (Pa), rho densità' fluido (kg/m^3). Standard internazionale ISO 5167-4:2003 specifica geometria tubo Venturi con angolo convergente 21° e divergente 5-15°, coefficient Cd 0,984-0,995 in funzione Reynolds Re=rho*v*D/mu e diameter ratio beta=D2/D1 (raccomandato 0,3-0,75).

Formule tubo Venturi e Bernoulli

Equazione di Bernoulli (1738): P1 + 0,5*rho*v1^2 + rho*g*h1 = P2 + 0,5*rho*v2^2 + rho*g*h2, conservazione energia per fluido ideale incompressibile in moto stazionario lungo linea di flusso. Daniel Bernoulli (Groningen Olanda-Basilea Svizzera 1700-1782) "Hydrodynamica" pubblicato Strasburgo 1738. Pressione totale P_t = P_statica + 0,5*rho*v^2 + rho*g*h conservata. Per condotta orizzontale (h1=h2): P1 + 0,5*rho*v1^2 = P2 + 0,5*rho*v2^2, riduzione pressione = aumento velocità'.

Equazione continuità' massa: A1*v1 = A2*v2 (per fluido incompressibile, conservazione massa portata costante). v2/v1 = A1/A2, restringimento sezione produce aumento velocità'. Esempio Venturi D1=100 mm, D2=50 mm: A1/A2=4, v2=4*v1. Combinato con Bernoulli: dP = P1-P2 = 0,5*rho*(v2^2-v1^2) = 0,5*rho*v2^2*(1-(A2/A1)^2).

Formula Venturi base: v2 = sqrt(2*dP/(rho*(1-(A2/A1)^2))). Per beta=D2/D1 ratio: v2 = sqrt(2*dP/(rho*(1-beta^4))). Portata teorica Q_th = A2*v2. Portata reale Q = Cd * Q_th = Cd * A2 * sqrt(2*dP/(rho*(1-beta^4))). Cd coefficient discharge accounts attrito viscoso, contrazione vena contracta, deviazioni profilo velocità' parabolico vs uniforme idealizzato. Cd standard ISO 5167-4 = 0,984-0,995 turbolento, scende 0,75-0,90 laminare.

Reynolds number e regime flusso: Re = rho*v*D/mu = v*D/nu. Soglia: Re4000 turbolento (Cd stabile 0,98). Acqua 20°C: rho=998,2 kg/m^3, mu=1,002e-3 Pa*s. Aria 20°C 1 atm: rho=1,204 kg/m^3, mu=1,825e-5 Pa*s. Olio motore SAE 30 a 40°C: rho=870 kg/m^3, mu=0,04 Pa*s. Sangue 37°C: rho=1060 kg/m^3, mu=4e-3 Pa*s (non-Newtoniano shear-thinning).

Beta ratio e perdite carico: beta = D2/D1 (typically 0,3-0,75 ISO 5167). Beta basso 0,3: perdite carico permanenti basse (10% dP), alta precisione, range portata limitato. Beta alto 0,75: perdite carico fino 50% dP, maggior range portata. Recovery cone divergente 7-15° riduce perdite permanenti vs orifizio brusco. Comparazione: orifizio plate Cd~0,6 perdite 65%, nozzle Cd~0,95 perdite 30%, Venturi Cd~0,98 perdite solo 10-15% (vantaggio principale).

Storia e contesto

Giovanni Battista Venturi (Italia 1746-1822) Bibbiano Reggio Emilia fisico professore Modena Universita' nel 1797 pubblica "Recherches experimentales sur le principe de la communication laterale du mouvement dans les fluides" Parigi descrivendo sperimentalmente effetto restringimento condotta, riduzione pressione associata aumento velocità' fluido. Daniel Bernoulli (Groningen Olanda 1700, Basilea Svizzera 1782) figlio Johann Bernoulli matematico, "Hydrodynamica" 1738 Strasburgo equazione conservazione energia fluidodinamica preliminare alla teoria moderna. Leonhard Euler (Basilea-San Pietroburgo 1707-1783) generalizzo' equazioni Eulero per fluido ideale 1755 conservazione momentum. Claude-Louis Navier (Francia 1785-1836) e George Stokes (Gran Bretagna 1819-1903) formularono equazioni Navier-Stokes 1822-1845 viscosità' fluido reale. Clemens Herschel (Germania-USA 1842-1930) ingegnere Holyoke Water Power Co. Massachusetts 1887 brevettò "Herschel Venturi tube" applicazione ingegneristica industriale: misura portata acqua canali industria tessile Lowell Massachusetts. Osborne Reynolds (Gran Bretagna 1842-1912) Manchester University 1883 introdusse numero Reynolds adimensionale Re=rho*v*D/mu transizione laminare-turbolento. Ludwig Prandtl (Germania 1875-1953) Gottingen University 1904 boundary layer theory perdite attrito parete. Theodore von Karman (Ungheria-USA 1881-1963) Caltech 1930s aerodinamica supersonica e turbolenza.

Italia industria misura portata: Itron Liberec gruppo USA contatori acqua-gas-elettricità' presenza Italia, Endress+Hauser Reinach Svizzera misura processo industriale Italia, Yokogawa Electric Tokio Giappone field instruments, Emerson Process Daniel Houston USA fornitore oil gas. ARERA Autorita' Regolazione Energia Reti Ambiente Milano regolatore italiano gestione gas naturale, ENI San Donato Milanese maggior compagnia petrolifera, Snam Rete Gas San Donato Milanese trasporto gas, Acea Roma servizio idrico Lazio, A2A Milano servizio idrico Lombardia, Hera Bologna servizio idrico Emilia. Carrier UTC Farmington Connecticut USA HVAC commerciale, Daikin Osaka Giappone HVAC residenziale, Trane Ingersoll Rand Davidson NC USA HVAC industriale. Centrali idroelettriche italiane: Enel Green Power Roma 17 GW idroelettrico nazionale, Edison Foro Buonaparte Milano principale operatore privato. Industria alimentare-farmaceutica: GEA Group Dusseldorf Germania, Tetra Pak Lund Svezia, IMA Industria Macchine Automatiche Bologna automazione farma. Medicale ventilatori polmonari: Philips Healthcare Eindhoven Olanda, Hamilton Medical Bonaduz Svizzera, Vyaire Medical Mettawa USA, Drager Lubecca Germania. INFN Frascati esperimenti dinamica fluidi acceleratori. Centro Italiano Ricerche Aerospaziali CIRA Capua Caserta gallerie vento. CNR-INM Istituto Marina Roma idrodinamica navale-marina.

Politecnico Milano Bovisa Dipartimento Ingegneria Aerospaziale, Politecnico Torino Ingegneria Meccanica e Aerospaziale corso Aerodinamica Industriale, Sapienza Roma DIMA, Universita' Federico II Napoli Ingegneria Industriale.

Tabella tipici parametri Venturi e fluidi

• Acqua 20°C: rho=998,2 kg/m^3, mu=1,002e-3 Pa*s, nu=1,004e-6 m^2/s
• Aria 20°C 1 atm: rho=1,204 kg/m^3, mu=1,825e-5 Pa*s, nu=1,516e-5 m^2/s
• Olio motore SAE 30 40°C: rho=870 kg/m^3, mu=0,04 Pa*s, nu=4,6e-5 m^2/s
• Diesel 20°C: rho=832 kg/m^3, mu=2,5e-3 Pa*s, nu=3e-6 m^2/s
• Benzina 20°C: rho=750 kg/m^3, mu=6e-4 Pa*s, nu=8e-7 m^2/s
• Gas naturale (metano) 0°C 1 atm: rho=0,668 kg/m^3, mu=1,1e-5 Pa*s
• Glicerolo 20°C: rho=1261 kg/m^3, mu=1,412 Pa*s (alta viscosità')
• Mercurio 20°C: rho=13546 kg/m^3, mu=1,526e-3 Pa*s
• Etanolo 20°C: rho=789 kg/m^3, mu=1,2e-3 Pa*s
• Latte 20°C: rho=1030 kg/m^3, mu=2,5e-3 Pa*s
• Sangue 37°C: rho=1060 kg/m^3, mu=4e-3 Pa*s (non-Newtoniano)
• Beta ratio raccomandato: 0,3-0,75 (ISO 5167-4)
• Cd Venturi turbolento: 0,984-0,995 (ISO 5167-4)
• Reynolds soglia laminare-turbolento: Re=2300 transizione
• Angolo convergente: 21° standard ISO 5167-4
• Angolo divergente: 5-15° recovery cone
• Perdite permanenti Venturi: 10-15% dP (vs 65% orifizio)
• Range portata tipico: 4:1 (rapporto Qmax/Qmin)
• Precisione misura: ±0,5-1% (calibrazione ISO/IEC 17025)
• Temperatura operativa: -200°C a +600°C (materiali speciali)

Esempi numerici dettagliati

Esempio 1 - acqua acquedotto urbano D1=100 mm. Tubo PE100 acqua potabile 100 mm, Venturi D2=50 mm, beta=0,5. Differenza pressione misurata manometro digitale dP=10 kPa=10000 Pa. Acqua 20°C rho=998,2 kg/m^3, mu=1,002e-3 Pa*s. A1=78,54 cm^2, A2=19,63 cm^2. v2=sqrt(2*10000/(998,2*(1-0,5^4)))=sqrt(20000/935,8)=4,62 m/s. v1=v2*0,25=1,16 m/s. Q=Cd*A2*v2=0,98*19,63e-4*4,62=8,89e-3 m^3/s=32 m^3/h=534 L/min. Re=998*1,16*0,1/1e-3=115500 turbolento alto Cd stabile. Applicazione contatori acqua condominiali Itron-Sensus.

Esempio 2 - aria HVAC ventilazione D1=200 mm. Condotto ventilazione zincato D1=200 mm, Venturi D2=80 mm, beta=0,4. dP=500 Pa misurato Pitot tube digitale. Aria 20°C 1 atm rho=1,204 kg/m^3. A1=314,16 cm^2, A2=50,27 cm^2. v2=sqrt(2*500/(1,204*(1-0,4^4)))=sqrt(831/0,9692)=29,3 m/s. v1=4,69 m/s. Q=0,98*50,27e-4*29,3=0,144 m^3/s=8,67 m^3/min=520 m^3/h. Re=1,204*4,69*0,2/1,825e-5=61870 turbolento. Applicazione AHU Carrier Daikin uffici 20 ricambi aria/ora.

Esempio 3 - olio motore industriale D1=300 mm. Oleodotto trasferimento Eni D1=300 mm, Venturi D2=150 mm, beta=0,5. dP=50 kPa=50000 Pa. Olio SAE 30 40°C rho=870 kg/m^3, mu=0,04 Pa*s. A1=706,86 cm^2, A2=176,71 cm^2. v2=sqrt(100000/(870*0,9375))=sqrt(122697)=11,07 m/s. v1=2,77 m/s. Q=0,98*176,71e-4*11,07=0,1917 m^3/s=690 m^3/h=11500 L/min. Re=870*2,77*0,3/0,04=18083 turbolento. Applicazione raffineria Eni Sannazzaro Pavia distribuzione interna.

Esempio 4 - gas naturale gasdotto D1=500 mm. Gasdotto Snam Rete Gas D1=500 mm, Venturi D2=200 mm, beta=0,4. dP=30 mbar=3000 Pa (basse perdite gas naturale). Gas naturale (metano) 0°C 1 atm rho=0,668 kg/m^3 (compressibile, qui flusso lento incompressibile). A1=1963,5 cm^2, A2=314,16 cm^2. v2=sqrt(2*3000/(0,668*0,9744))=sqrt(9214)=96 m/s (alto). v1=15,4 m/s. Q=0,98*314,16e-4*96=2,96 m^3/s. Re=0,668*15,4*0,5/1,1e-5=4,67e5 turbolento alto. Per gas richiede correzione comprimibilita' Y=1-(0,41+0,35*beta^4)*dP/(kappa*P1) per dP/P1>0,01.

Esempio 5 - sangue micro-flusso medicale D1=6 mm. Cateter vascolare D1=6 mm, Venturi flow meter D2=3 mm, beta=0,5. dP=100 Pa (basso, bambino pediatrico). Sangue 37°C rho=1060 kg/m^3, mu=4e-3 Pa*s. A1=0,2827 cm^2, A2=0,07069 cm^2. v2=sqrt(200/(1060*0,9375))=sqrt(0,2014)=0,449 m/s. v1=0,112 m/s. Q=0,98*0,07069e-4*0,449=3,11e-6 m^3/s=11,2 mL/min=672 mL/h. Re=1060*0,112*0,006/4e-3=178 LAMINARE (Cd ridotto ~0,75). Applicazione monitor cardiopolmonare Philips IntelliVue ICU pediatrica Bambin Gesu Roma.

Limiti e considerazioni

1. Formula Venturi valida solo per fluido incompressibile (acqua, olio liquidi). Per gas richiede correzione Y comprimibilita' se dP/P1>0,01.
2. Cd coefficient discharge varia con Reynolds: 0,75-0,90 laminare, 0,95-0,99 turbolento alto. Calibrazione ISO 5167-4 raccomandata.
3. Beta ratio raccomandato 0,3-0,75: beta0,75 non standard.
4. Profilo velocità' parabolico laminare riduce Cd: assumere flusso pienamente sviluppato 30*D distanza monte.
5. Cavitation rischio se P_gola < pressione vapor saturo: riduce portata e danneggia gola Venturi.
6. Misura dP richiede manometro differenziale precisione ±0,1% (Endress+Hauser, Yokogawa, Emerson).
7. Fluido non-Newtoniano (sangue, polimeri): viscosità' apparente shear-rate dipendente, Cd diversi.
8. Tubo Venturi vs orifizio plate: Venturi perdite 10-15% vs orifizio 65%, ma costo 5-10x.
9. Standard ISO 5167-4:2003 specifica geometria classica Herschel: convergente 21°, divergente 5-15°.
10. Bidirezionalita': Venturi standard misura flusso unidirezionale, bidirezionale richiede design speciale.

Errori comuni

1. Confondere D1 (tubo) con D2 (gola): formula richiede beta=D2/D1<1, mai invertito.
2. Trascurare correzione comprimibilita' Y per gas con dP/P1>0,01: errore portata fino 20%.
3. Assumere Cd=1 ideale: realistico 0,98 turbolento, 0,75-0,90 laminare. Errore 2-25%.
4. Confondere Pa con kPa o mbar: conversione 1 kPa=1000 Pa, 1 mbar=100 Pa.
5. Trascurare densità' fluido a temperatura operativa: acqua 4°C 1000 kg/m^3 vs 80°C 971 kg/m^3 (3% errore).
6. Misurare pressione monte/valle Venturi distanza errata: ISO 5167 specifica tap location precisa.
7. Trascurare regime laminare per Re<2300: Cd diverso, formula precisa diversa (Hagen-Poiseuille).
8. Sottostimare cavitation per liquidi: P_gola = P1 - 0,5*rho*(v2^2-v1^2) deve > P_vapor saturo.
9. Assumere fluido Newtoniano per sangue/polimeri: viscosità' shear-rate dipendente cambia Cd.
10. Confondere portata massica m_dot=rho*Q (kg/s) con volumetrica Q (m^3/s).

Applicazioni numerate

1. Acquedotti urbani: Acea Roma servizio idrico Lazio, A2A Milano Lombardia, Hera Bologna Emilia distribuzione acqua potabile.
2. Gasdotti gas naturale: Snam Rete Gas San Donato Milanese 32000 km gasdotti Italia, distribuzione metano nazionale.
3. Oleodotti petroliferi: Eni Saipem San Donato Milanese estrazione e trasporto greggio raffinerie Sannazzaro Taranto.
4. HVAC commerciale-residenziale: Carrier UTC, Daikin, Trane Ingersoll Rand sistemi ventilazione AHU edifici.
5. Industria alimentare-farma: GEA Group Dusseldorf, Tetra Pak Lund, IMA Bologna Industria Macchine Automatiche.
6. Medicale ventilatori polmonari: Philips IntelliVue, Hamilton Medical Bonaduz, Drager Lubecca ICU ospedalieri.
7. Centrali idroelettriche: Enel Green Power 17 GW Italia, Edison Foro Buonaparte Milano operatore privato.
8. Aerospazio gallerie vento: CIRA Capua Caserta wind tunnel, Politecnico Milano Bovisa dipartimento aerospaziale.
9. Strumentazione misura: Endress+Hauser Reinach Svizzera, Yokogawa Tokio, Emerson Daniel Houston field instruments.
10. Educational didattica: Politecnico Milano Torino Sapienza Federico II corsi Fluidodinamica Ingegneria Industriale.

Per supporto

Per domande tecniche o segnalazioni: info@calcoi.com. Strumento didattico per fluidodinamica preliminare. Per progettazione professionale ISO 5167-4:2003 'Measurement of fluid flow by means of pressure differential devices - Part 4: Venturi tubes', Bernoulli D 1738 'Hydrodynamica' Strasburgo, Venturi GB 1797 'Recherches experimentales' Parigi, White FM 'Fluid Mechanics' 2016 McGraw-Hill, Munson 'Fundamentals of Fluid Mechanics' 2017 Wiley.

Aggiornato a luglio 2026