Corso di idraulica: 1 Introduzione all’Idraulica e all’Idrodinamica

1 – Introduzione all’idraulica e all’idrodinamica

Questa trattazione fornisce i concetti fondamentali dell’idraulica/ idrodinamica, le equazioni di base e numerosi esempi numerici svolti passo-passo. Unitá di misura: SI (m, s, kg, Pa, N).


1. Definizioni fondamentali

Fluido

Un fluido è una sostanza che si deforma continuamente sotto l’azione di una tensione tangenziale infinitesima; comprende liquidi e gas.

Densità ρ\rho

Massa per unità di volume:

ρ=mV[kgm3]\rho = \frac{m}{V}\quad [\mathrm{kg\,m^{-3}}]

Pressione pp

Forza per unità di superficie:

p=FA[Pa=Nm2]p = \frac{F}{A}\quad [\mathrm{Pa} = N\,m^{-2}]

Esempio: una forza F=100 NF=100\ \mathrm{N} su un’area A=0.02 m2A=0{.}02\ \mathrm{m^2}p=100/0.02=5000 Pap=100/0.02=5000\ \mathrm{Pa}.

Portata (volumetrica) QQ e portata di massa m˙\dot m

  • Portata volumetrica:

Q=dVdt[m3s1]Q = \frac{dV}{dt}\quad[\mathrm{m^3\,s^{-1}}]
  • Portata di massa:

m˙=ρQ[kgs1]\dot m = \rho\,Q \quad [\mathrm{kg\,s^{-1}}]

Velocità media vv e relazione con QQ

Per sezione di area AA:

Q=Av v=QAQ = A\,v \qquad \Rightarrow\ v = \frac{Q}{A}

Viscosità

  • Viscosità dinamica (o assoluta) μ\mu [Pa·s] (o N·s·m⁻²): misura della resistenza allo scorrimento.

  • Viscosità cinematica ν\nu [m²·s⁻¹]:

ν=μρ.\nu = \frac{\mu}{\rho}.

Esempio numerico (acqua a ~20°C): μ1.002×103 Pas\mu \approx 1{.}002\times10^{-3}\ \mathrm{Pa\cdot s}, ρ998 kgm3\rho \approx 998\ \mathrm{kg\,m^{-3}}

ν=1.002×1039981.004×106 m2/s.\nu=\frac{1{.}002\times10^{-3}}{998}\approx 1{.}004\times10^{-6}\ \mathrm{m^2/s}.

(Valore calcolato: ν1.0040106 m2/s\nu \approx 1{.}0040\cdot 10^{-6}\ \mathrm{m^2/s}.)


2. Fluidi ideali vs fluidi reali

  • Fluido ideale (inviscid, incompressibile): μ=0\mu = 0 e spesso ρ=\rho= costante. Modello utile per semplificare equazioni (es. equazione di Bernoulli).

  • Fluido reale: μ>0\mu>0; presenta attrito viscoso, transizione laminare/turbulenta, strati limite, perdite d’energia. I gas possono essere trattati come comprimibili quando le variazioni di densità non sono trascurabili.

Osservazioni pratiche:

  • Liquidi (acqua) spesso si trattano incomprimibili (buona approssimazione) se velocità ≪ velocità del suono.

  • Gas: assumere incomprimibilità solo se Mach Ma<0.3Ma < 0{.}3 (circa).


3. Equazioni dimensionali e numeri adimensionali

scopo

I numeri adimensionali condensano il rapporto tra effetti fisici concorrenti (inertia vs viscosità, gravità vs inerzia, compressibilità vs inerzia, ...). Sono fondamentali per similitudine e scala.

Reynolds ReRe — rapporto inerzia/viscosità

Re=ρvLμ=vLνRe = \frac{\rho\, v \, L}{\mu} = \frac{v\,L}{\nu}
  • LL: lunghezza caratteristica (es. diametro di una tubazione DD).

  • Interpretazione: alto ReRe → forze d’inerzia dominanti → tendenza alla turbolenza; basso ReRe → viscosità dominante → flusso laminare.

Regole pratiche (tubo circolare):

  • Re<2300Re < 2300 → laminare,

  • 2300Re40002300 \lesssim Re \lesssim 4000 → transizione,

  • Re>4000Re > 4000 → turbolento.

Esempio 1 (calcolo Re): acqua a 20°C: ρ=998 kg/m3\rho=998\ \mathrm{kg/m^3}, μ=1.002×103 Pas\mu=1.002\times10^{-3}\ \mathrm{Pa\cdot s}. Prendiamo tubo D=0.05 mD=0{.}05\ \mathrm{m} e v=1.5 m/sv=1{.}5\ \mathrm{m/s}.

Re=9981.50.051.002×10374700.Re=\frac{998\cdot 1{.}5\cdot 0{.}05}{1{.}002\times 10^{-3}}\approx 74\,700.

Risultato: Re7.47×104Re\approx 7{.}47\times10^{4} → regime turbulento.
(Valore calcolato: Re74700.6Re\approx 74\,700.6.)


Froude FrFr — rapporto inerzia/gravitazione (importante per superfici libere)

Due definizioni equivalenti:

Fr=vgLoFr2=v2gL.Fr = \frac{v}{\sqrt{g\,L}} \quad\text{o}\quad Fr^2 = \frac{v^2}{gL}.
  • LL = lunghezza caratteristica (es. profondità, lunghezza nave).

  • Fr1Fr \ll 1: effetti gravitazionali dominano (onde lente), Fr1Fr\gg1: inerzia domina.

Esempio: v=1.5 m/sv=1{.}5\ \mathrm{m/s}, L=0.5 mL=0{.}5\ \mathrm{m}:

Fr=1.59.810.50.677.Fr = \frac{1{.}5}{\sqrt{9{.}81\cdot 0{.}5}}\approx 0{.}677.

(Valore calcolato: Fr0.6773Fr\approx0{.}6773; Fr20.4587Fr^2\approx0{.}4587.)


Mach MaMa — rapporto velocità/c

Ma=vcMa = \frac{v}{c}
  • cc: velocità del suono nel mezzo (aria c343 m/sc\approx343\ \mathrm{m/s} a 20°C).

  • Regola pratica: se Ma<0.3Ma<0{.}3, effetti comprimibili trascurabili.

Esempio: v=100 m/sv=100\ \mathrm{m/s} in aria:

Ma=1003430.292<0.3Ma=\frac{100}{343}\approx0{.}292 < 0{.}3

→ approssimazione incompressibile ancora accettabile ma vicina al limite.


4. Principi generali: continuità, equilibrio, conservazione

4.1 Equazione di continuità (conservazione della massa)

Forma integrale (controllo di volume VV):

ddtVρdV+SρvndS=0\frac{d}{dt}\int_{V}\rho\, dV + \oint_{S}\rho\, \mathbf{v}\cdot \mathbf{n}\, dS = 0

Forma differenziale (generale):

ρt+(ρv)=0.\frac{\partial \rho}{\partial t} + \nabla\cdot(\rho\,\mathbf{v}) = 0.

Incompressibile (ρ\rho costante) si riduce a:

v=0.\nabla\cdot\mathbf{v} = 0.

Esempio pratico (1-D stazionario): in regime stazionario la massa che entra = massa che esce → Q1=Q2Q_1=Q_2.
Se Q=0.01 m3/sQ=0{.}01\ \mathrm{m^3/s}, area A1=πD12/4A_1=\pi D_1^2/4 con D1=0.05 mD_1=0{.}05\ \mathrm{m} e D2=0.02 mD_2=0{.}02\ \mathrm{m}. Calcoliamo le velocità:

  • A1=π(0.05)2/40.0019635 m2A_1 = \pi (0{.}05)^2/4 \approx 0{.}0019635\ \mathrm{m^2},

  • A2=π(0.02)2/40.00031416 m2A_2 = \pi (0{.}02)^2/4 \approx 0{.}00031416\ \mathrm{m^2},

v1=QA10.010.00196355.093 m/s,v2=QA20.010.0003141631.831 m/s.v_1=\frac{Q}{A_1}\approx \frac{0.01}{0.0019635}\approx 5{.}093\ \mathrm{m/s}, \quad v_2=\frac{Q}{A_2}\approx \frac{0.01}{0.00031416}\approx 31{.}831\ \mathrm{m/s}.

(Valori calcolati: v15.09296 m/sv_1\approx 5{.}09296\ \mathrm{m/s}, v231.83099 m/sv_2\approx 31{.}83099\ \mathrm{m/s}.)


4.2 Equilibrio e conservazione della quantità di moto (Navier–Stokes)

Forma vettoriale generale (per fluido viscoso, compressibile):

ρ(vt+vv)=p+μ2v+ρg+festerni\rho\left(\frac{\partial \mathbf{v}}{\partial t} + \mathbf{v}\cdot\nabla\mathbf{v}\right) = -\nabla p + \mu \nabla^2 \mathbf{v} + \rho \mathbf{g} + \mathbf{f}_{\text{esterni}}
  • Termine di sinistra: accelerazione (inerziale).

  • p-\nabla p: gradiente di pressione.

  • μ2v\mu\nabla^2\mathbf{v}: viscous diffusion (forze viscose).

  • ρg\rho\mathbf{g}: peso specifico.

Limite inviscido (fluido ideale): μ=0\mu=0equazione di Euler:

ρ(vt+vv)=p+ρg.\rho\left(\frac{\partial \mathbf{v}}{\partial t} + \mathbf{v}\cdot\nabla\mathbf{v}\right) = -\nabla p + \rho \mathbf{g}.

4.3 Conservazione dell’energia e equazione di Bernoulli

Per un fluido incompressibile, stazionario e inviscido lungo una stessa linea di corrente:

p+12ρv2+ρgz=costantep + \tfrac12 \rho v^2 + \rho g z = \text{costante}
  • pp: pressione statica.

  • 12ρv2\tfrac12\rho v^2: pressione dinamica (energia cinetica per unità volume).

  • ρgz\rho g z: energia potenziale per unità volume.

Forma in termini di carichi (head) dividendo per ρg\rho g:

pρg+v22g+z=costante (head).\frac{p}{\rho g} + \frac{v^2}{2g} + z = \text{costante (head)}.

Se si includono perdite (pipe con attrito), versione energetica:

p1ρg+v122g+z1=p2ρg+v222g+z2+hL\frac{p_1}{\rho g} + \frac{v_1^2}{2g} + z_1 = \frac{p_2}{\rho g} + \frac{v_2^2}{2g} + z_2 + h_L

dove hLh_L è la perdita di carico dovuta a attrito e turbolenza.


5. Perdite per attrito: Darcy–Weisbach

Per tubazione circolare:

hf=fLDv22gh_f = f \,\frac{L}{D}\,\frac{v^2}{2g}
  • ff: fattore di attrito (Darcy). Per flusso turbolento l’interpretazione e calcolo di ff richiede la carta di Moody o formule approssimate (es. Blasius).

Blasius (approssimazione per flussi turbolenti lisci, 4.000Re1054{.}000\lesssim Re\lesssim10^5):

f0.3164Re0.25f \approx 0{.}3164\,Re^{-0.25}

Esempio 3 (perdita di carico e pressione)
Usiamo il caso già calcolato: acqua, D=0.05 mD=0{.}05\ \mathrm{m}, v=1.5 m/sv=1{.}5\ \mathrm{m/s}, L=100 mL=100\ \mathrm{m}, Re74.7006Re\approx74{.}7006.

  1. Calcolare ff con Blasius:

f0.3164Re0.250.3164(7.47006×104)0.250.01914.f\approx 0.3164\cdot Re^{-0.25}\approx 0{.}3164\cdot (7.47006\times10^4)^{-0.25}\approx 0{.}01914.

(Valore calcolato: f0.019138f\approx 0.019138.)

  1. Calcolare perdita di carico hfh_f:

hf=fLDv22g=0.0191381000.051.5229.814.3895 m.h_f = f\frac{L}{D}\frac{v^2}{2g} =0.019138\cdot\frac{100}{0.05}\cdot\frac{1.5^2}{2\cdot 9.81} \approx 4{.}3895\ \mathrm{m}.

(Valore calcolato: hf4.38954 mh_f\approx 4{.}38954\ \mathrm{m}.)

  1. Pressione persa (Δp) corrispondente:

Δp=ρghf9989.814.3895442.975 kPa.\Delta p = \rho\,g\,h_f \approx 998\cdot 9.81\cdot 4.38954 \approx 42{.}975\ \mathrm{kPa}.

(Valore calcolato: Δp42.975 kPa\Delta p\approx 42{.}975\ \mathrm{kPa} ≃ 0.43 bar.)

  1. Potenza idraulica richiesta dal pompaggio (senza tenere conto di efficienza):

P=ΔpQP = \Delta p \cdot Q

con Q=Av=πD2/4vQ=A v = \pi D^2/4 \cdot v . Per D=0.05 mD=0{.}05\ \mathrm{m}, A0.0019635 m2A\approx 0{.}0019635\ \mathrm{m^2} e Q0.00294524 m3/sQ\approx 0{.}00294524\ \mathrm{m^3/s}. Quindi

P42.975103 Pa0.00294524 m3/s126.6 W.P \approx 42{.}975\cdot10^3\ \mathrm{Pa}\cdot 0{.}00294524\ \mathrm{m^3/s} \approx 126{.}6\ \mathrm{W}.

(Valore calcolato: P126.57 WP\approx126{.}57\ \mathrm{W}.) Se la pompa ha efficienza η\eta, potenza elettrica Pelec=P/ηP_\text{elec}=P/\eta.


6. Bernoulli: esempio di velocità da dislivello (vasca → ugello)

Problema: acqua scorre da una vasca libera con livello a quota superiore di Δz=10 m\Delta z = 10\ \mathrm{m}; trascurando perdite, calcolare la velocità di uscita vv.

Applicando Bernoulli tra livello della vasca (dove v0v\approx0) e uscita:

ρgΔz=12ρv2v=2gΔz.\rho g \Delta z = \tfrac12 \rho v^2 \quad\Rightarrow\quad v=\sqrt{2 g \Delta z}.

Con g=9.81 m/s2g=9{.}81\ \mathrm{m/s^2}:

v=29.811014.007 m/s.v=\sqrt{2\cdot 9{.}81\cdot 10}\approx 14{.}007\ \mathrm{m/s}.

(Valore calcolato: v14.0071 m/sv\approx 14{.}0071\ \mathrm{m/s}.)

Portata se l’ugello ha diametro D=0.05 mD=0{.}05\ \mathrm{m}:

Q=Av=πD24v0.001963514.00710.02750 m3/s (27.50 L/s).Q=A v=\frac{\pi D^2}{4}v\approx 0{.}0019635\cdot 14{.}0071 \approx 0{.}02750\ \mathrm{m^3/s}\ (27{.}50\ \mathrm{L/s}).

(Massa fluente: m˙=ρQ9980.0275027.45 kg/s\dot m = \rho Q \approx 998\cdot0{.}02750\approx 27{.}45\ \mathrm{kg/s}.)


7. Sommario pratico dei passaggi di calcolo (checklist)

  1. Identificare le grandezze caratteristiche: LL (diametro, profondità), vv (velocità media), ρ\rho, μ\mu.

  2. Calcolare QQ o vv con Q=AvQ=A v.

  3. Stimare ReRe per conoscere il regime.

  4. Se necessario calcolare ff (Moody/Blasius) e perdite hfh_f.

  5. Applicare Bernoulli con/o senza perdite per stime di velocità o pressione.

  6. Convertire head ↔ pressione: p=ρghp=\rho g h.

  7. Stimare potenza: P=ΔpQP=\Delta p\cdot Q e dividere per η\eta per potenza elettrica richiesta.


8. Riferimenti rapidi alle formule principali

  • Densità: ρ=m/V\rho=m/V.

  • Pressione: p=F/Ap=F/A.

  • Portata volumetrica: Q=AvQ=A v.

  • Portata massa: m˙=ρQ\dot m=\rho Q.

  • Viscosità cinematica: ν=μ/ρ\nu=\mu/\rho.

  • Reynolds: Re=ρvLμ=vLνRe=\dfrac{\rho v L}{\mu}=\dfrac{v L}{\nu}.

  • Froude: Fr=vgLFr=\dfrac{v}{\sqrt{gL}}.

  • Mach: Ma=vcMa=\dfrac{v}{c}.

  • Continuità (differenziale): ρt+(ρv)=0\dfrac{\partial \rho}{\partial t}+\nabla\cdot(\rho\mathbf{v})=0.

  • Navier–Stokes: ρ(tv+vv)=p+μ2v+ρg\rho(\partial_t\mathbf{v}+\mathbf{v}\cdot\nabla\mathbf{v})=-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}.

  • Bernoulli (invicid, incompressible): p+12ρv2+ρgz=costantep+\tfrac12\rho v^2+\rho g z=\text{costante}.

  • Darcy–Weisbach: hf=fLDv22gh_f=f\dfrac{L}{D}\dfrac{v^2}{2g}.


9. Note pratiche e limiti

  • Le equazioni semplificate (Bernoulli, Blasius) sono utilissime per stime veloci ma vanno applicate entro i loro limiti: flusso stazionario, incompressibile, senza fenomeni transitori violenti, ecc.

  • La turbolenza richiede modelli (RANS, LES) per predizione accurata; il fattore di attrito ff dipende anche dalla rugosità relativa ϵ/D\epsilon/D.

  • Per gas a velocità elevate considerare la comprimibilità (campo di validità MaMa).

  • Per applicazioni ingegneristiche reali si eseguono sempre prove sperimentali o simulazioni CFD per verificare i calcoli teorici.

Commenti

Post popolari in questo blog

Corso di chimica: Reazioni chimiche

Corso di Stampa 3D: 3 – Modelling 3D – Software

Corso di Fondamenti di Matematica e Logica: 3 Insiemi Relazioni e Funzioni