Formula di Chézy

Nell'ambito dell'ingegneria idraulica, la formula di Chézy, o formula di Chézy-Tadini[senza fonte], sviluppata dall'ingegnere francese Antoine de Chézy, è una formula empirica utilizzata per calcolare la velocità di un fluido a pelo libero in condizioni di moto uniforme[1] e principalmente turbolento[2], il cui campo di applicazione è stato poi ampliato anche alle correnti in pressione.

Storia

I presupposti della formula possono essere ricondotti ad Albert Brahms, geometra di un principato tedesco, il quale scrisse che, a differenza di una sfera posta su un piano inclinato, l'acqua fluente in un canale inclinato non si muove di moto uniformemente accelerato, bensì di moto uniforme perché l'accelerazione è controbilanciata dall'attrito; affermò inoltre che le velocità sono proporzionali alla radice quadrata della pendenza del fondo e scrisse che:[3]

(EN)

«the values of friction at equal slopes of water surface are to each other, in case of open flowing waters, as the areas wetted by the water are to the quantities that flow over them.»

(IT)

«i valori di attrito a parità di pendenze della superficie dell'acqua sono tra loro, nel caso di correnti a pelo libero, come le aree bagnate dall'acqua sono rispetto alle quantità che fluiscono su di esse.»

(Albert Brahms, citato da Clemens Herschel in 115 experiments on the carrying capacity of large, riveted, metal conduits, up to six feet per second of velocity of flow[4])

Tuttavia Brahms non fornì alcuna formula.[3]

Ritratto di Antoine Chézy
Ritratto di Antoine de Chézy eseguito da Louis Jean Desprez nel 1772–76

Nel 1769, l'ingegnere Jean-Rodolphe Perronet e il suo assistente Antoine de Chézy, all'epoca ispettore generale dei ponti e delle strade, erano stati designati a relazionare sul progetto di un canale presso il fiume Yvette, al fine di fornire acqua a Parigi;[3] a partire da quell'anno, Chézy raccolse dati sperimentali sul canale in terra Courpalet e sul fiume Senna[1], che inviò insieme alle sue conclusioni a Perronet in un documento francese datato 1775 con il titolo Tesi sulla velocità di flusso in un dato canale.[5]

Nella relazione intitolata Formule pour Trouver la Vitesse de l'Eau Conduit dan une Rigole donnée[6][7], datata 1776, riportò la formula:[5]

v = χ ⋅ A / P ⋅ i f {\displaystyle v=\chi \cdot {\sqrt {A/P\cdot i_{f}\;}}}

dove χ {\displaystyle \chi } è un fattore di resistenza al flusso, i f {\displaystyle i_{f}} è la pendenza di fondo, A {\displaystyle A} e P {\displaystyle P} sono rispettivamente l'area bagnata e il perimetro bagnato della sezione trasversale.

Chézy fornì inoltre il valore del coefficiente χ {\displaystyle \chi } calcolato per il canale Courpalet, pari a 31 m1/2/s, e per il fiume Senna, pari a 44 m1/2/s.[8][5]

Perronet non diede risalto al lavoro svolto da Chézy, il quale non ricercò mai visibilità; fu solo nel 1797 che Louis-François Letourneur, membro del Direttorio francese, assegnò a Chézy l'incarico di direttore dell'École des Ponts et Chaussées (l'attuale École des Ponts ParisTech); nel 1804, successivamente alla morte di Chézy[9], Gaspard de Prony e Pierre-Charles Lesage, ne riconobbero il merito dei contributi scientifici.[10]

La formula

La formula è la seguente:

v = χ ⋅ R h ⋅ J {\displaystyle v=\chi \cdot {\sqrt {R_{h}\cdot J}}}

dove:

  • v {\displaystyle v} è la velocità media nella sezione idraulica trasversale in m/s;
  • χ {\displaystyle \chi } è il coefficiente di Chézy in m1/2/s;
  • R h = A / P {\displaystyle R_{h}=A/P} è il raggio idraulico della sezione trasversale in metri, pari al rapporto tra area bagnata e perimetro bagnato nella sezione trasversale;
  • J {\displaystyle J} è la cadente piezometrica, adimensionale.

Nella sua forma originariamente proposta[5], è valida per il solo moto uniforme, essendo J = i f {\displaystyle J=i_{f}} , dove i f {\displaystyle i_{f}} è la pendenza del fondo alveo.

Nel caso di correnti a pelo libero, viene utilizzata anche nella seguente forma:[11]

Q = χ ⋅ A ⋅ R h ⋅ J {\displaystyle Q=\chi \cdot A\cdot {\sqrt {R_{h}\cdot J}}}

dove:

  • Q {\displaystyle Q} è la portata in m3/s;
  • A {\displaystyle A} è l'area bagnata della sezione trasversale in m2.

Derivazione matematica

La formula di Chézy può essere derivata matematicamente da due assunzioni:[12]

  1. la forza F r {\displaystyle F_{r}} resistente al flusso del liquido, agente su un'area A f {\displaystyle A_{f}} di fondo alveo, è direttamente proporzionale alla velocità v {\displaystyle v} del flusso al quadrato, cioè F r A f = C ⋅ v 2 {\displaystyle {\frac {F_{r}}{A_{f}}}=C\cdot v^{2}} , essendo C {\displaystyle C} una costante di proporzionalità, ed esprimendo l'area in funzione del perimetro bagnato della sezione trasversale fluviale e della lunghezza del tratto considerato ( A f = P ⋅ L {\displaystyle A_{f}=P\cdot L} ) si ha F r = C ⋅ P ⋅ L ⋅ v 2 {\displaystyle F_{r}=C\cdot P\cdot L\cdot v^{2}} ;
  2. nel moto uniforme, la componente effettiva della forza di gravità che causa il moto, cioè la componente parallela al fondo alveo F g , p {\displaystyle F_{g,p}} , è uguale alla forza resistente al moto stesso F r {\displaystyle F_{r}} ; quindi sapendo che F g , p = γ w ⋅ A ⋅ L ⋅ sin ⁡ ( θ ) ≈ γ w ⋅ A ⋅ L ⋅ tan ⁡ ( θ ) = γ w ⋅ A ⋅ L ⋅ i f {\displaystyle F_{g,p}=\gamma _{w}\cdot A\cdot L\cdot \sin(\theta )\approx \gamma _{w}\cdot A\cdot L\cdot \tan(\theta )=\gamma _{w}\cdot A\cdot L\cdot i_{f}} [13], essendo γ w {\displaystyle \gamma _{w}} il peso specifico del liquido, A {\displaystyle A} l'area della sezione trasversale bagnata e i f {\displaystyle i_{f}} la pendenza del fondo alveo, segue che v 2 = γ w C ⋅ A P ⋅ i f {\displaystyle v^{2}={\frac {\gamma _{w}}{C}}\cdot {\frac {A}{P}}\cdot i_{f}} . Ponendo le quantità costanti pari al coefficiente χ 2 = γ w C {\displaystyle \chi ^{2}={\frac {\gamma _{w}}{C}}} e sostituendo con la definizione di raggio idraulico R h = A P {\displaystyle R_{h}={\frac {A}{P}}} , si ha infine: v = χ ⋅ R h ⋅ i f {\displaystyle v=\chi \cdot {\sqrt {R_{h}\cdot i_{f}}}} .

Estensione alle correnti in pressione

Nonostante la formula di Chezy abbia origine nel '700, è più recente la sua estensione alle correnti in pressione.[non chiaro]

Tenendo presente che R h {\displaystyle R_{h}} è il raggio idraulico, dato dal rapporto

R h = A P {\displaystyle R_{h}={\frac {A}{P}}}

in cui A rappresenta l'area della sezione del condotto e P rappresenta il contorno, o perimetro, del condotto.

R h = D / 4 {\displaystyle R_{h}=D/4} (per un condotto a sezione circolare)

e considerata J come la cadente dei carichi totali la formula di Chézy pur uguale nella scrittura è valida per le correnti in pressione.[non chiaro]

V = χ R h ⋅ J {\displaystyle V=\chi {\sqrt {R_{h}\cdot J}}}

Tuttavia è spesso utilizzata per calcolare la cadente dei carichi, nella forma:

J = V 2 χ 2   R h {\displaystyle J={\frac {V^{2}}{\chi ^{2}\ R_{h}}}}

Coefficiente di Chézy

L'equazione di Chézy, se confrontata con l'equazione di Darcy-Weisbach, permette di ottenere il coefficiente di Chézy in funzione dell'indice di resistenza λ {\displaystyle \lambda } , ottenendo:[14]

χ = 8 g λ {\displaystyle \chi ={\sqrt {\frac {8g}{\lambda }}}}

Tuttavia, data la poca praticità del calcolo di λ {\displaystyle \lambda } , per la determinazione del coefficiente di Chézy, diversi autori hanno fornito espressioni empiriche, cosiddette "pratiche", valide solo nel caso di moto puramente turbolento data l'assenza di dipendenza dal numero di Reynolds;[2] nel seguito si riportano le formule principali espresse nel S.I..

formula di Kutter semplificata

Inizialmente proposta da E. Ganguillet e W. R. Kutter nel 1869 in una forma più complessa, ottenuta da elaborazioni basate sulle misure di deflusso in canali di vario tipo, incluse quelle effettuate da Bazin pubblicate nel 1865[15] e quelle effettuate su fiumi europei e sul Mississippi[16], è stata successivamente semplificata nella seguente forma:[14]

χ = 100 1 + m R h {\displaystyle \chi ={\frac {100}{1+{\frac {m}{\sqrt {R_{h}}}}}}}

in cui m {\displaystyle m} è l'indice di scabrezza, ha le dimensioni della radice quadrata di una lunghezza ed è variabile tra 0,12 m1/2 per tubazioni in acciaio nuove e 0,45 m1/2 per tubazioni in ghisa con forti incrostazioni.[14]

formula di Bazin

Proposta da H. Bazin nel 1897, sulla base di dati ottenuti da piccoli canali sperimentali, nella sua seconda formulazione è equivalente alla formula di Kutter, ma considerata meno soddisfacente di quest'ultima[17] e ha la seguente forma:[14]

χ = 87 1 + γ R h {\displaystyle \chi ={\frac {87}{1+{\frac {\gamma }{\sqrt {R_{h}}}}}}}

in cui γ {\displaystyle \gamma } è l'indice di scabrezza, ha le dimensioni della radice quadrata di una lunghezza, di cui Bazin fornisce alcuni valori tabulati.[17]

Eguagliando l'espressione di Kutter con quella di Bazin è possibile ottenere la relazione tra i due coefficienti m {\displaystyle m} e γ {\displaystyle \gamma } :

γ = 0 , 87 ⋅ m − 0 , 13 ⋅ R h {\displaystyle \gamma =0,87\cdot m-0,13\cdot {\sqrt {R_{h}}}} .

formula di Gauckler-Strickler

Lo stesso argomento in dettaglio: Formula di Gauckler-Strickler.

Nel 1868 Philippe Gaspard Gauckler propose la seguente formula per il coefficiente di Chézy:

χ = k s ⋅ R h 1 / 6 {\displaystyle \chi =k_{s}\cdot R_{h}^{1/6}}

in cui k s {\displaystyle k_{s}} è chiamato "coefficiente di Gauckler-Strickler", ha le dimensioni di [m1/3/s], è una misura inversamente proporzionale alla scabrezza della parete ed è reperibile in tabelle, essendo variabile tra circa 140 m1/3/s per i tubi nuovi in acciaio a 65 m1/3/s per tubi in ghisa con forti incrostazioni.[18]

Strickler presentò nel 1923, in maniera indipendente da Gauckler, la stessa formula proposta[19], oltre che una scala di valori per l'indice di scabrezza[18].

Essendo una relazione di tipo monomio, risulta di semplice utilizzo per applicazioni analitiche.

formula di Manning

Inizialmente proposta nel 1889 da Robert Manning, sulla base di sette formule differenti basate sulle misure di Bazin[15] e di 170 osservazioni di taratura, è stata successivamente modificata nella forma seguente:[20]

v = 1 n ⋅ R h 2 / 3 ⋅ J {\displaystyle v={\frac {1}{n}}\cdot R_{h}^{2/3}\cdot {\sqrt {J}}}

in cui n {\displaystyle n} è il coefficiente di scabrezza, detto "numero di Manning", ha le dimensioni di [s/m1/3], e misura in maniera direttamente proporzionale la scabrezza della superficie.

Confrontata con la formula di Chézy fornisce l'espressione del coefficiente di Chézy in funzione del numero di Manning:[21]

χ = 1 n ⋅ R h 1 / 6 {\displaystyle \chi ={\frac {1}{n}}\cdot R_{h}^{1/6}}

Per la scelta del coefficiente di Manning è stata fornita una metodologia di calcolo sulla base delle caratteristiche dell'alveo[22], oltre che dei valori tabellati[23].

Confrontando tale espressione con quella di Gauckler-Strickler si determina che il numero di Manning è pari al reciproco del coefficiente di Gauckler-Strickler:

n = 1 k s {\displaystyle n={\frac {1}{k_{s}}}}

Note

  1. ^ a b Chow (1959), p. 93.
  2. ^ a b Çengel et al. (2007), p. 278.
  3. ^ a b c Clemens (1897), p. 73.
  4. ^ sulla base del testo originale Anfangs-Gründe der Deich - und Wasser-Baukunst, Brahms (1754)
  5. ^ a b c d Khoury (2004).
  6. ^ (EN) Nikolaos D. Katopodes, Free-Surface Flow:: Shallow Water Dynamics, Butterworth-Heinemann, 30 agosto 2018, p. 462, ISBN 978-0-12-815488-5. URL consultato il 27 aprile 2021.
  7. ^ che tradotto in italiano è "Formula per trovare la velocità uniforme che l'acqua avrà in un fosso o in un canale di cui sia nota la pendenza"
  8. ^ i valori dei coefficienti riportati sono stati convertiti nell'unità di misura del sistema internazionale, a partire dal sistema di misura anticamente utilizzato in Francia, presente nel documento originale.
  9. ^ avvenuta nel 1798
  10. ^ Clemens (1897), pp. 117-119.
  11. ^ ottenuta moltiplicando ambo i membri per A {\displaystyle A} , ricordando che Q = v ⋅ A {\displaystyle Q=v\cdot A}
  12. ^ Chow (1959), pp. 93-94.
  13. ^ per valori di θ {\displaystyle \theta } piccoli, come risultano essere nella variabilità delle pendenze fluviali, si può confondere il seno con la tangente ( sin ⁡ ( θ ) ≈ tan ⁡ ( θ ) {\displaystyle \sin(\theta )\approx \tan(\theta )} ); ad esempio, per un'inclinazione θ = 20 ∘ {\displaystyle \theta =20^{\circ }} , pari a una pendenza i f = tan ⁡ ( 20 ∘ ) = 0 , 3640 {\displaystyle i_{f}=\tan(20^{\circ })=0,3640} , l'errore che si commette è t a n ( 20 ∘ ) − s i n ( 20 ∘ ) = 0 , 022 {\displaystyle tan(20^{\circ })-sin(20^{\circ })=0,022}
  14. ^ a b c d Çengel et al. (2007), p. 279.
  15. ^ a b (FR) Henry Darcy e Henry Émile Bazin, Recherches hydrauliques enterprises, Imprimerie Nationale, 1865, OCLC 38698151.
  16. ^ Chow (1959), p. 94.
  17. ^ a b Chow (1959), p. 95.
  18. ^ a b Çengel et al. (2007), p. 280.
  19. ^ Chow (1959), p. 99.
  20. ^ Chow (1959), pp. 98-99.
  21. ^ Chow (1959), p. 100.
  22. ^ Cowan (1956), pp. 473-475.
  23. ^ Chow (1959), pp. 110-113.

Bibliografia

  • Gianfranco Becciu, Fondamenti di costruzioni idrauliche, UTET scienze tecniche, 2010, ISBN 978-88-598-0522-9, OCLC 875255671. URL consultato il 13 giugno 2021.
  • (DE) Albert Brahms, Anfangs-Gründe der Deich - und Wasser-Baukunst, Schuster, 1754, ISBN 3-7963-0273-4, OCLC 612058507. Ospitato su http://www.manuscriptorium.com/apps/index.php?direct=record&pid=AIPDIG-STK___B142SYSTEMNU0SMC7B5-cs.
  • Yunus A. Çengel e John M. Cimbala, Meccanica dei fluidi, a cura di Giuseppe Cozzo e Cinzia Santoro, McGraw-Hill, 2007, ISBN 978-88-386-6384-0, OCLC 799749775.
  • Duilio Citrini e Giorgio Noseda, Idraulica, 2ª ed., Milano, Casa Editrice Ambrosiana, 1987, ISBN 8808081044.
  • (EN) Ven Te Chow, Open-channel Hydraulics, New York, McGraw-Hill, 1959, ISBN 978-0-07-085906-7.
  • (EN) Herschel Clemens, 115 experiments on the carrying capacity of large, riveted, metal conduits, up to six feet per second of velocity of flow, New York, John Wiley & Sons, 1897.
  • (EN) Woody L. Cowan, Estimating hydraulic roughness coefficients, in Agricultural Engineering, vol. 37, n. 7, luglio 1956.
  • (EN) Fadi Khoury, History of the Chézy Formula, su chezy.sdsu.edu, giugno 2004. URL consultato il 23 marzo 2021.

Voci correlate

  Portale Fisica
  Portale Ingegneria