Validazione del calcolo

Un software di calcolo si giudica su una cosa sola: i numeri. In questa pagina pubblichiamo casi canonici svolti per esteso con le formule della pratica progettuale, replicabili con una calcolatrice, e li confrontiamo cifra per cifra con i risultati del motore di Gasnetics.

Il metodo: nessuna scatola nera

Trasparenza verso la norma e verso il calcolo manuale.

Il motore di Gasnetics implementa la correlazione di Renouard in forma quadratica, la formula di riferimento della pratica progettuale italiana per le reti di distribuzione in bassa e media pressione (il quadro normativo di riferimento è il DM 16-17/04/2008 con la UNI 9165):

P₁² − P₂² = 48,6 · s · L · Q1,82 / D4,82   [bar assoluti]

con s densità del gas relativa all'aria, L lunghezza in metri, Q portata in Sm³/h e D diametro interno in mm. Sulle reti magliate il motore risolve il bilancio nodale con il metodo di Newton-Raphson; sui casi di questa pagina — scelti apposta perché risolvibili anche a mano — la soluzione è ricalcolabile con una calcolatrice scientifica.

Tre regole di onestà:

  1. I numeri di questa pagina sono generati dal motore di produzione, lo stesso che esegue i calcoli nell'applicazione: non sono trascritti né aggiustati.
  2. Il confronto è con la formula enunciata e con il calcolo manuale, non con altri software: la fonte di verità è la norma e la letteratura, che chiunque può consultare.
  3. Ogni relazione di calcolo generata da Gasnetics include la stessa verifica: un'appendice «Sviluppo del calcolo» svolge per esteso un tratto del percorso critico della tua rete e ricalcola a mano tutti i tratti del percorso, mostrando lo scarto rispetto ai valori della relazione.

Gas usato in tutti i casi: gas naturale tipico di rete (93% CH₄, 4% C₂H₆, 2% N₂, 1% CO₂ molare), densità relativa s = 0,591215, densità ρ = 0,724179 kg/Sm³ (aria: 1,224900 kg/Sm³). Le pressioni relative assumono l'atmosfera standard 1,01325 bar. Gli sviluppi sono mostrati con valori arrotondati, ma il calcolo sottostante è in piena precisione: replicando con le cifre stampate si ritrova il risultato entro il decimillesimo di mbar.

Caso 1 — Tratto singolo in bassa pressione

Il mattone elementare: una sorgente, un'utenza, un tubo. Tutto verificabile riga per riga.

Pressione sorgente25 mbar rel
Lunghezza L100 m
Diametro interno D50,0 mm
Portata Q30 Sm³/h

Sviluppo del calcolo a mano

  1. Pressione assoluta di monte: P₁ = 1,01325 + 0,025 = 1,03825 bar.
  2. Resistenza del tratto: R = 48,6 · 0,591215 · 100 / 504,82 = 1,859285·10⁻⁵ bar²/(Sm³/h)1,82.
  3. Termine di portata: Q1,82 = 301,82 = 487,9336.
  4. Caduta quadratica: P₁² − P₂² = R · Q1,82 = 9,072075·10⁻³ bar².
  5. Pressione di valle: P₂ = √(1,03825² − 9,072075·10⁻³) = 1,0338718 bar = 20,6218 mbar rel.

Confronto con il motore

GrandezzaCalcolo a manoMotore GasneticsScarto
P₂ utenza [mbar rel]20,621820,6218< 10⁻¹⁰ mbar
Caduta di pressione [mbar]4,37824,3782< 10⁻¹⁰ mbar

Il motore risolve lo stesso caso come sistema nodale (7 iterazioni di Newton-Raphson fino a convergenza) e riporta anche la velocità del gas nel tratto: 4,16 m/s, entro il limite di bassa pressione. Rifacendo i conti con la densità relativa arrotondata a 4 decimali (s = 0,5912) lo scarto resta di circa 0,0001 mbar: il caso è replicabile davvero con una calcolatrice.

Caso 2 — Rete ramificata in bassa pressione

Un albero a tre utenze: le portate dei rami sono note per bilancio, quindi ogni tratto si verifica a mano in sequenza.

Sorgente S a 30 mbar rel (1,04325 bar ass). Da S parte una dorsale verso il nodo N1 (che consegna 25 Sm³/h), e da N1 due rami verso N2 (45 Sm³/h) e N3 (30 Sm³/h). Per bilancio ai nodi la dorsale trasporta 25+45+30 = 100 Sm³/h — e infatti il motore calcola esattamente 100,000000 / 45,000000 / 30,000000 Sm³/h sui tre tratti.

Sviluppo per tratto

TrattoL [m]D [mm]Q [Sm³/h]R [bar²/(Sm³/h)1,82]Q1,82P₁²−P₂² [bar²]
S → N1 (dorsale)1501001009,873534·10⁻⁷4365,15834,309954·10⁻³
N1 → N220060451,544267·10⁻⁵1020,57931,576047·10⁻²
N1 → N312060309,265601·10⁻⁶487,93364,520998·10⁻³

Applicando la caduta di ciascun tratto in cascata (P del nodo di valle = √(P²monte − caduta)):

Confronto con il motore

NodoCalcolo a mano [mbar rel]Motore Gasnetics [mbar rel]Scarto
N127,932327,9323< 10⁻¹⁰ mbar
N220,336120,3361< 10⁻¹⁰ mbar
N325,759025,7590< 10⁻¹⁰ mbar

Su una rete magliata (con anelli) la ripartizione delle portate non è più calcolabile a mano in sequenza — è esattamente il problema che il metodo nodale risolve — ma la fisica di ogni singolo tratto resta questa, ed è quella che l'appendice della relazione ricalcola e riconcilia tratto per tratto sul percorso critico.

Caso 3 — Correzione altimetrica

Un tratto in salita di 60 m: il gas naturale, più leggero dell'aria, in salita guadagna pressione relativa.

Pressione sorgente25 mbar rel
Lunghezza L300 m
Diametro interno D50,0 mm
Portata Q20 Sm³/h
Dislivello Δz+60 m

La correzione altimetrica usa la variabile ŷ = (P − c·(z − z medio))², con c = (ρaria − ρgas)·g/10⁵ espresso in bar/m — la spinta di galleggiamento della colonna di gas rispetto alla colonna d'aria. Le quote sono riferite alla quota media della rete perché solo i dislivelli hanno effetto fisico.

Sviluppo del calcolo a mano

  1. Coefficiente altimetrico: c = (1,224900 − 0,724179) · 9,80665 / 10⁵ = 4,910399·10⁻⁵ bar/m ≈ 0,0491 mbar per metro di salita.
  2. Quota media z medio = 30 m; correzione al nodo di monte: c·(0 − 30) = −1,473120·10⁻³ bar; al nodo di valle: +1,473120·10⁻³ bar.
  3. ŷ₁ = (1,03825 + 1,473120·10⁻³)² = 1,0810242 bar².
  4. Attrito: R·Q1,82 = 5,577855·10⁻⁵ · 233,2784 = 1,301193·10⁻² bar².
  5. ŷ₂ = 1,0810242 − 1,301193·10⁻² = 1,0680122 bar²; √ŷ₂ = 1,0334468 bar.
  6. Pressione di valle: P₂ = 1,0334468 + 1,473120·10⁻³ = 1,0349199 bar = 21,6699 mbar rel.

Confronto con il motore

GrandezzaCalcolo a manoMotore GasneticsScarto
P₂ utenza [mbar rel]21,669921,6699< 10⁻¹⁰ mbar

Il controllo di coerenza fisica: lo stesso tratto in piano consegnerebbe 18,7147 mbar. La salita di 60 m «regala» 2,9552 mbar (≈ 60 × 0,0491), perché il gas naturale è più leggero dell'aria — l'effetto, controintuitivo ma ben noto ai progettisti, che rende obbligatoria la correzione altimetrica sulle reti con dislivelli. Con un gas più denso dell'aria (GPL) il segno si inverte, e il motore lo gestisce con la stessa formula.

La stessa verifica, nella tua relazione

Questa pagina non chiede un atto di fiducia: ogni relazione di calcolo generata da Gasnetics contiene un'appendice «Sviluppo del calcolo e verifica di coerenza» che svolge per esteso, con i numeri della tua rete, un tratto del percorso critico, e ricalcola a mano la pressione di valle di tutti i tratti del percorso mostrando lo scarto rispetto ai valori della relazione. La coincidenza si controlla su ogni progetto, non solo sui casi di questa pagina.

I risultati del calcolo sono un supporto alla progettazione: la scelta progettuale e la verifica finale restano responsabilità del professionista.

Prova Gasnetics