title: Equazione di Ergun
text: # Equazione di Ergun
  
  L'equazione di Ergun, ricavata dall'ingegnere chimico turco Sabri Ergun nel 1952, descrive le perdite di carico lungo un reattore a letto fisso.
  
  Tale equazione può essere scritta come:
  
      ${\frac {dP}{dz}}=-{\frac {G}{\rho g_{c}D_{p}}}\left({\frac {1-\phi }{\phi ^{3}}}\right)\left[{\frac {150\left(1-\phi \right)\mu }{D_{p}}}+1,75G\right]$
  
  dove:
  
  * $P=$ pressione, espressa in [Pa] (SI) oppure [lbf/ft²] (US)
  * $z=$ lunghezza del letto, espressa in [m] (SI) oppure [ft] (US)
  * $G=\rho u=$ velocità superficiale massica, espressa in [kg/m²·s] (SI) oppure [lbm/ft²·h] (US)
  * $\rho =$ densità del fluido, espressa in [kg/m³] (SI) oppure [lbm/ft³] (US)
  * $u=$ velocità superficiale, espressa in [m/s] (SI) oppure [ft/h] (US)
  * $g_{c}=$ fattore di conversione, che vale 1 (SI) oppure 32,174 lbm·ft/lbf·s² (US)
  * $D_{p}=$ diametro delle particelle di catalizzatore nel letto, espresso in [m] (SI) oppure [ft] (US)
  * $\phi =$ grado di vuoto del letto
  * $\mu =$ viscosità dinamica del fluido, espressa in [kg/m·s] (SI) oppure [lbm/ft²·h] (US).
  
  ## Caso gas
  
  Se il fluido che attraversa il letto è un gas, l'unico parametro che cambia con la pressione lungo il letto è la densità del gas stesso. Operando in condizioni stazionarie, la portata massica entrante Q0 nel letto eguaglia la portata massica uscente Q dal letto, per cui:
  
      $Q_{0}=Q$
      $\rho _{0}v_{0}=\rho v$
  
  con $v_{0}$ portata volumetrica entrante e $v$ portata volumetrica uscente. Per la legge dei gas perfetti si ha che:
  
      $v=v_{0}\left({\frac {P_{0}}{P}}\right)\left({\frac {T}{T_{0}}}\right)\left({\frac {F_{T}}{F_{T0}}}\right)$
  
  con $T$ temperatura e $F_{T}$ portata molare totale, da cui
  
      $\rho =\rho _{0}{\frac {v_{0}}{v}}=\rho _{0}\left({\frac {P}{P_{0}}}\right)\left({\frac {T_{0}}{T}}\right)\left({\frac {F_{T0}}{F_{T}}}\right)$
  
  Andando a sostituire tale relazione nell'Equazione di Ergun si ottiene:
  
      ${\frac {dP}{dz}}=-{\frac {G}{\rho _{0}g_{c}D_{p}}}\left({\frac {1-\phi }{\phi ^{3}}}\right)\left[{\frac {150\left(1-\phi \right)\mu }{D_{p}}}+1,75G\right]\left({\frac {P_{0}}{P}}\right)\left({\frac {T}{T_{0}}}\right)\left({\frac {F_{T}}{F_{T0}}}\right)$
  
  che si può semplificare, tenendo conto che i primi tre fattori del membro di destra sono costanti lungo il letto per condizioni di ingresso fissate, in:
  
      ${\frac {dP}{dz}}=-\beta _{0}\left({\frac {P_{0}}{P}}\right)\left({\frac {T}{T_{0}}}\right)\left({\frac {F_{T}}{F_{T0}}}\right)=-\beta _{0}\left({\frac {1}{y}}\right)\left({\frac {T}{T_{0}}}\right)\left({\frac {F_{T}}{F_{T0}}}\right)$
  
  Dividendo per la pressione a monte del letto:
  
      $\left({\frac {1}{P_{0}}}\right){\frac {dP}{dz}}=-\left({\frac {1}{P_{0}}}\right)\beta _{0}\left({\frac {1}{y}}\right)\left({\frac {T}{T_{0}}}\right)\left({\frac {F_{T}}{F_{T0}}}\right)$
      ${\frac {d\left(P/P_{0}\right)}{dz}}=-\beta _{0}\left({\frac {1}{y\cdot P_{0}}}\right)\left({\frac {T}{T_{0}}}\right)\left({\frac {F_{T}}{F_{T0}}}\right)$
      ${\frac {dy}{dz}}=-\beta _{0}\left({\frac {1}{y\cdot P_{0}}}\right)\left({\frac {T}{T_{0}}}\right)\left({\frac {F_{T}}{F_{T0}}}\right)$
  
  Ipotizzando di essere in condizioni isoterme, gli ultimi due fattori del membro di destra si semplificano:
  
      ${\frac {dy}{dz}}=-\beta _{0}\left({\frac {1}{y\cdot P_{0}}}\right)$
  
  dalla quale infine è possibile ricavare il rapporto tra pressione finale e iniziale lungo il letto,
  
      $y={\frac {P}{P_{0}}}={\sqrt {1-{\frac {2\beta _{0}z}{P_{0}}}}}$
  