Equazione di stato: differenze tra le versioni

In chimica fisica e termodinamica, una equazione di stato è una equazione costitutiva che descrive lo stato della materia sotto un dato insieme di condizioni fisiche. Fornisce una relazione matematica tra due o più funzioni di stato associate alla materia, come temperatura, pressione, volume o energia interna. Le equazioni di stato sono utili nella la descrizione delle proprietà dei fluidi (e delle loro miscele), dei solidi e persino per descrivere l'interno delle stelle.

Il principale uso di una equazione di stato consiste nella determinazione dei parametri di stato relativi a gas e liquidi. Una di queste equazioni, tra le più semplici e utilizzate, è la legge dei gas perfetti valida per gas in condizioni di bassa pressione e alta temperatura (comportamento ideale). Tale trattazione matematica diviene però poco accurata in condizioni lontane dall'idealità e non è in grado di predire la liquefazione di un gas. Sono state quindi sviluppate diverse equazioni di stato più accurate per gas e liquidi. Attualmente non esiste una singola equazione di stato che descriva accuratamente le proprietà di tutte le sostanze nelle diverse condizioni possibili.

Oltre a equazioni che predicono il comportamento di gas e liquidi ci sono anche equazioni di stato per predire il volume dei solidi, ivi compresa la transizione tra stati cristallini differenti. Ci sono equazioni che modellano l'interno delle stelle, comprese le stelle di neutroni. Un concetto correlato è quello di fluido perfetto.

La legge di Boyle è stata forse la prima espressione di una equazione di stato. Nel 1662 Robert Boyle, un irlandese, compì una serie di esperimenti impiegando un tubo di vetro a J sigillato da un lato. Mise mercurio nel tubo, intrappolando così una piccola quantità di aria nel lato sigillato del tubo, dopodiché misurò accuratamente il volume del gas man mano che aggiungeva mercurio nel tubo. La pressione del gas poteva essere determinata dalla differenza di livello del mercurio nei due lati del tubo. Attraverso questi esperimenti Boyle notò che il volume del gas variava in maniera inversa rispetto alla pressione. In termini matematici questo può essere espresso come:

PV = costante

Questa relazione è stata attribuita anche a Edme Mariotte e ci si riferisce ad essa anche come Legge di Mariotte. L'opera di Mariotte però non venne pubblicata che nel 1676.

Nel 1787 il fisico francese Jacques Charles scoprì che ossigeno, azoto, idrogeno, anidride carbonica e aria si espandono allo stesso modo lungo il medesimo intervallo di 80K. Successivamente, nel 1802, Joseph Louis Gay-Lussac pubblicò risultati di esperienze simili, che indicavano una relazione lineare tra volume e temperatura

V₁/T₁ = V₂/T₂

Nel 1801 l'inglese John Dalton scoprì che la pressione totale esercitata da una miscela di n gas può essere definita come la somma

${\displaystyle P=p_{1}+p_{2}+...+p_{n}}$

dove ${\displaystyle p_{1}}$, ${\displaystyle p_{2}}$, ${\displaystyle p_{n}}$ rappresentano la pressione parziale di ogni componente.

Nel 1834 Émile Clapeyron combinò Le leggi di Boyle e di Charles nella prima formulazione della legge dei gas ideali. Inizialmente l'equazione era scritta come PV_m=R(T_C+267) (con la temperatura espressa in gradi Celsius). Indagini successive però rivelarono che il numero avrebbe dovuto essere 273.2, e in seguito la scala Celsius venne definita come 0 °C = 273.15 K, dando quindi:

PV_m=R(T_C+273.15)

La legge di Amagat (1880) afferma che il volune di una miscela gassosa a una data temperatura e pressione è la somma del volume che avrebbero i singoli componenti alle medesime condizioni. Questa è la prova sperimentale che il volume è una grandezza estensiva. In termini matematici:

${\displaystyle V=v_{1}+v_{2}+...+v_{n}}$

Nelle seguenti equazioni le variabili sono definite come segue. Si può utilizzare qualsiasi insieme coerente di unità di misura, anche se sono da preferire le unità SI

${\displaystyle P}$ pressione

${\displaystyle V}$ volume

${\displaystyle \rho }$ densità

${\displaystyle n}$ numero di moli di sostanza

${\displaystyle {\tilde {V}}={\frac {V}{n}}}$ volume molare, il volume di una mole di liquido o gas

${\displaystyle {\tilde {\rho }}={\frac {1}{\tilde {V}}}}$ densità molare

${\displaystyle T}$ temperatura (K)

${\displaystyle R}$ costante dei gas (8.314472 J/(mol·K))

L'equazione dei gas perfetti può essere scritta come:

${\displaystyle PV=nRT\,}$

L'equazione di stato dei gas ideali può anche essere espressa come:

${\displaystyle P=\rho (\gamma -1)e\,}$

dove ${\displaystyle \rho }$ è la densità, ${\displaystyle \gamma }$ l'indice adiabatico espresso dal rapporto tra il calore specifico a pressione costante e quello a volume costante (vedi anche trasformazione adiabatica), ed e l'energia interna. Questa forma è scritta puramente in termini di variabili intensive ed è molto utilizzata in simulazioni che implicano l'uso dell'equazioni di Eulero in quanto esprime la relazione tra l'energia interna e altre forme di energia (come quella cinetica), permettendo in tal modo il rispetto della conservazione dell'energia totale.

${\displaystyle \left(P+{\frac {a}{{\tilde {V}}^{2}}}\right)\left({\tilde {V}}-b\right)=RT}$

I parametri a, b e R sono costanti che dipendono dal fluido specifico. Possono essere così calcolati dalle proprietà critiche:

${\displaystyle a=3P_{c}V_{c}^{2}}$

${\displaystyle b={\frac {V_{c}}{3}}}$

${\displaystyle R={\frac {8P_{c}V_{c}}{3T_{c}}}}$

Proposta nel 1873, l'equazione di stato di Van der Waals è stata una delle prime equazioni a fornire risultati spiccatamente piu precisi dell'equazione di stato dei gas perfetti. In questa fondamentale equazione a viene chiamato parametro attrattivo e b parametro repulsivo o covolume. Pur essendo ben superiore all'equazione dei gas perfetti e nonostante predica correttamente la formazione di una fase liquida, l'accordo coi dati sperimentali è limitato nelle condizioni in cui il liquido si forma. Benché sia comunemente citata nei libri di testo e nelle pubblicazioni scientifiche per ragioni storiche e didattiche è ormai caduta in disuso: altre equazioni più moderne, di complessità solo lievemente maggiore, sono molto più accurate.

L'equazione di Van der Waals può essere considerata un "miglioramento" della legge dei gas ideali, sotto due aspetti:

1. Le molecole vengono considerate come particelle con un volume, non come punti materiali. Quindi V non può essere troppo piccolo, minore di una data costante caratteristica in funzione del tipo di gas. Abbiamo quindi l'aggiunta del termine correttivo (V - b) che sostituisce il semplice V.
2. Le molecole del gas ideale non interagiscono reciprocamente, in questo modo invece vengono considerate le forze attrattive intermolecolari, forze avvertibili entro una distanza di parecchi raggi molecolari. L'effetto finale risulta in una diminuzione della pressione totale esercitata dalle molecole di gas, per cui il termine P deve essere corretto addizionandogli un fattore che renda conto di tale diminuzione. Essendo la forza che agisce su ogni superficie di molecola proporzionale alla densità (~${\displaystyle \rho }$), la forza che agisce sull'insieme di molecole del gas è ~${\displaystyle \rho ^{2}}$~${\displaystyle {\frac {1}{V_{m}^{2}}}}$.

${\displaystyle \ P(V-b)=RTe^{-a/RTV}}$

Dove a è un parametro associato all'interazione tra le molecole e b tiene conto delle dimensioni delle molecole, similmente all'equazione di Van der Waals.

I parametri critici valgono:

${\displaystyle \ T_{c}={\frac {a}{4Rb}},P_{c}={\frac {a}{4e^{2}b^{2}}},V_{c}=2b}$

${\displaystyle {\frac {P{\tilde {V}}}{RT}}=1+{\frac {B}{\tilde {V}}}+{\frac {C}{{\tilde {V}}^{2}}}+{\frac {D}{{\tilde {V}}^{3}}}+\dots }$

${\displaystyle B=-V_{c}\,}$

${\displaystyle C={\frac {V_{c}^{2}}{3}}}$

${\displaystyle R={\frac {P_{c}V_{c}}{T_{c}}}}$

Benché non sia l'equazione di stato più conveniente da utilizzare, l'equazione del viriale è importante perché può essere derivata direttamente dalla meccanica statistica. Con appropriate assunzioni sulla forma matematica delle forze intermolecolari, si possono derivare espressioni analitiche per ognuno dei coefficienti. In questo caso B corrisponde a interazioni tra coppie di molecole, C a interazioni tra triplette, e così via. L'accuratezza raggiunta è maggiore con termini superiori.

${\displaystyle P={\frac {RT}{{\tilde {V}}-b}}-{\frac {a}{{\sqrt {T}}{\tilde {V}}\left({\tilde {V}}+b\right)}}}$

${\displaystyle a={\frac {0.42748R^{2}T_{c}^{2.5}}{P_{c}}}}$

${\displaystyle b={\frac {0.08664RT_{c}}{P_{c}}}}$

Introdotta nel 1949, l'equazione di stato Redlich-Kwong è stata un considerevole miglioramento rispetto ad altre equazioni di quel tempo. È ancora di un certo interesse, soprattutto per la sua forma relativamente semplice. Benché superiore all'equazione di Van der Waals, non è molto precisa riguardo alla fase liquida e non può quindi essere impiegata per un calcolo accurato degli equilibri liquido-vapore. Comunque può essere impiegata per questo scopo con l'ausilio di correlazioni separate per la fase liquida.

l'equazione di stato Redlich-Kwong è adeguata per il calcolo delle proprietà dei gas in condizioni in cui il rapporto tra la pressione e la pressione critica è minore della metà del rapporto tra la temperatura e la temperatura critica.

${\displaystyle P={\frac {RT}{{\tilde {V}}-b}}-{\frac {a\alpha (T)}{{\tilde {V}}\left({\tilde {V}}+b\right)}}}$

${\displaystyle a={\frac {0.42747R^{2}T_{c}^{2}}{P_{c}}}}$

${\displaystyle b={\frac {0.08664RT_{c}}{P_{c}}}}$

${\displaystyle \alpha (T)=\left(1+\left(0.48508+1.55171\omega -0.15613\omega ^{2}\right)\left(1-T_{r}^{0.5}\right)\right)^{2}}$

${\displaystyle T_{r}={\frac {T}{T_{c}}}}$

Dove ω è il fattore acentrico per la specie.

Per l'idrogeno:

${\displaystyle \alpha =1.202\exp \left(-0.30288T_{r}\right)}$

Nel 1972 Soave rimpiazzò il termine a/√(T) dell'equazione Redlich-Kwong con una funzione ${\displaystyle \alpha (T)}$ comprendente la temperatura e il fattore acentrico. La funzione α venne sviluppata per calcolare la pressione di vapore degli idrocarburi e l'equazione si comporta molto bene riguardo questi composti.

${\displaystyle P={\frac {RT}{{\tilde {V}}-b}}-{\frac {a\alpha (T)}{{\tilde {V}}^{2}+2b{\tilde {V}}-b^{2}}}}$

${\displaystyle a={\frac {0.45724R^{2}T_{c}^{2}}{P_{c}}}}$

${\displaystyle b={\frac {0.07780RT_{c}}{P_{c}}}}$

${\displaystyle \alpha (T)=\left(1+\left(0.37464+1.54226\omega -0.26992\omega ^{2}\right)\left(1-T_{r}^{0.5}\right)\right)^{2}}$

${\displaystyle T_{r}={\frac {T}{T_{c}}}}$

Dove ω è il fattore acentrico per la specie.

L'equazione Peng-Robinson venne sviluppata nel 1976 per soddisfare i seguenti requisiti:

1. I parametri avrebbero dovuto essere esprimibili in termini di proprietà critiche e fattore acentrico
2. Il modello avrebbe dovuto essere ragionevolmente accurato in prossimità del punto critico, particolarmente per il calcolo del fattore di comprimibilità e della densità del liquido
3. Le regole di miscelazione non avrebbero dovuto impiegare più di un singolo parametro di interazione binaria, parametro indipendente da temperatura, pressione e composizione
4. L'equazione avrebbe dovuto essere applicabile a tutti i calcoli delle proprietà dei fluidi utilizzati nell'ambito dei processi che utilizzavano gas naturale

Nel suo complesso, l'equazione di Peng e Robinson ha un'accuratezza simile a quella dell'equazione Redlich-Kwong-Soave, anche se è generalmente superiore nel predire la densità liquida di molti fluidi, specialmente quelli non polari.

${\displaystyle P={\tilde {\rho }}RT+\left(B_{0}RT-A_{0}-{\frac {C_{0}}{T^{2}}}\right){\tilde {\rho }}^{2}+\left(bRT-a-{\frac {d}{T}}\right){\tilde {\rho }}^{3}+\alpha a{\tilde {\rho }}^{6}+{\frac {c{\tilde {\rho }}^{3}}{T^{2}}}\left(1+\gamma {\tilde {\rho }}^{2}\right)\exp \left(-\gamma {\tilde {\rho }}^{2}\right)}$

dove ${\displaystyle \rho }$ è la densità molare (ottenibile come inverso del volume molare). L'equazione è stata scritta in termini di densità molare per motivi di praticità.

Le equazioni di stato BWR venne proposta all'inizio degli anni quaranta per l'utilizzo nell'industria petrolchimica (si ricordi che all'epoca l'equazione RK non era ancora stata introdotta), ed ebbe un grande successo. Il suo difetto sta nel richiedere parecchi parametri, peraltro di significato fisico non ben definito.

${\displaystyle P={\tilde {\rho }}RT+\left(B_{0}RT-A_{0}-{\frac {C_{0}}{T^{2}}}+{\frac {D_{0}}{T^{3}}}-{\frac {E_{0}}{T^{4}}}\right){\tilde {\rho }}^{2}+\left(bRT-a-{\frac {d}{T}}\right){\tilde {\rho }}^{3}+\alpha \left(a+{\frac {d}{T}}\right){\tilde {\rho }}^{6}+{\frac {c{\tilde {\rho }}^{3}}{T^{2}}}\left(1+\gamma {\tilde {\rho }}^{2}\right)\exp \left(-\gamma {\tilde {\rho }}^{2}\right)}$

È una versione più precisa dell'equazione precedente.

Valori dei vari parametri per 15 sostanze si trovano in: K.E. Starling, Fluid Properties for Light Petroleum Systems. Gulf Publishing Company (1973).

L'equazione di stato Elliott, Suresh e Donohue (ESD) è stata proposta nel 1990. L'equazione cerca di correggere una scarsa accuratezza dell'equazione Peng-Robinson per quanto riguarda il termine repulsivo di Van der Waals. L'equazione tiene conto degli effetti di forma di una molecola non polare e può essere estesa ai polimeri tramite l'aggiunta di un ulteriore termine (qui non mostrato). L'equazione è stata sviluppata attraverso simulazioni al computer e dovrebbe essere in grado di descrivere le proprietà fisiche essenziale delle dimensioni e della forma delle molecole, nonché le proprietà del legame idrogeno.

${\displaystyle {\frac {P{\tilde {V}}}{RT}}=Z=1+{\frac {4\left\langle c\eta \right\rangle }{1-1.9\eta }}-{\frac {9.5\left\langle qY\eta \right\rangle }{1+1.7745\left\langle Y\eta \right\rangle }}}$

Dove:

c = fattore di forma (shape factor)

${\displaystyle \eta =b\rho }$

${\displaystyle q=1+1.90476(c-1)}$

${\displaystyle Y=\exp({\frac {\epsilon }{kT}})-1.0617}$

Considerando l'acqua sottoposta e pressioni molto alte (condizione che si verifica tipicamente ad esempio nelle esplosioni nucleari sottacqua, nella litotripsia ultrasonica e nella sonoluminescenza) è utilizzabile l'equazione di stato irrigidita:

${\displaystyle p=\rho (\gamma -1)u-\gamma p^{0}\,}$

dove ${\displaystyle u}$ è l'energia interna per unità di massa, ${\displaystyle \gamma }$ è una costante empirica che assume solitamente valore 6.1 e ${\displaystyle p^{0}}$ è un'altra costante che rende conto dell'attrazione molecolare tra molecole di acqua. Il valore di quest'ultima è di circa 20000 atmosfere.

L'equazione assume questa forma perché la velocità del suono in acqua è data da ${\displaystyle c^{2}=\gamma (p+p^{0})/\rho }$.

Quindi l'acqua si può assumere che si comporti come un gas ideale che si trovi già sotto la pressione di 20000 atm, e ciò spiega il perché l'acqua è comunemente ritenuta incomprimibile: quando la pressione esterna varia da 1 a 2 atmosfere (da 100 kPa a 200 kPa), l'acqua si comporta analogamente a un gas ideale che passa da 20001 a 20002 atm (da 200.01 a 201.02 MPa).

Questa equazione non è in grado di predirre il calore specifico dell'acqua ma esistono poche alternative da utilizzare per processi che implicano condizioni tanto estreme.

Un fluido ultrarelativistico possiede equazione di stato

${\displaystyle p=c_{s}^{2}\mu }$

dove ${\displaystyle p}$ è la pressione, ${\displaystyle \mu }$ la densità di energia e ${\displaystyle c_{s}}$ una costante che si identifica con la velocità del suono.

L'equazione di stato per un gas ideale di Bose è

${\displaystyle PV_{m}=RT~{\frac {{\textrm {Li}}_{\alpha +1}(z)}{\zeta (\alpha )}}\left({\frac {T}{T_{c}}}\right)^{\alpha }}$

dove α è un esponente specifico del sistema (per esempio, in assenza di un campo potenziale α=3/2), z vale exp(μ/kT) dove μ è il potenziale chimico, Li è il polilogaritmo, ζ è la funzione zeta di Riemann e T_c è la temperatura critica alla quale inizia a formarsi il condensato di Bose-Einstein.

• Benedict, M., Webb, G. B., and Rubin, L. C., "An Empirical Equation for Thermodynamic Properties of Light Hydrocarbons and Their Mixtures: I. Methane, Ethane, Propane, and n-Butane", J. Chem. Phys., Vol. 8, No.4, pp. 334-345 (1940).
• ibid, "II. Mixtures of Methane, Ethane, Propane and n-Butane", J. Chem. Phys., Vol. 10, No. 12, pp. 747 - 758 (1942).
• ibid, "III. Constants for Twelve Hydrocarbons", Chem. Eng. Progr., Vol. 47, No. 8, pp. 419-422 (1951).
• ibid, "IV. Fugacities and Liquid-Vapor Equilibria", Chem. Engr. Progr., Vol. 47, No. 9, pp. 449-454 (1951).
• G. Soave, Equilibrium Constants from a Modified Redlich-Kwong. Equation of State, Chem. Eng. Sci., Vol. 27 pp. 1197-1203 (1972).
• Ding-Yu Peng,D. B. Robinson; Ind. Eng. Chem. Fundam., Vol. 15(1), pp. 59. [10] (1976).
• J.R. Elliott, Jr., S.J. Suresh, and M.D. Donohue, A simple equation of state for non-spherical and associating molecules, Ind. Eng. Chem. Res., Vol.29:1476 (1990).
• Elliott & Lira, "Introductory Chemical Engineering Thermodynamics", Prentice Hall (1999).

Portale Chimica

Portale Fisica