Inthermodynamics, a departure function is defined for any thermodynamic property as the difference between the property as computed for an ideal gas and the property of the species as it exists in the real world, for a specified temperature T and pressure P. Common departure functions include those for enthalpy, entropy, and internal energy.
Departure functions are used to calculate real fluid extensive properties (i.e. properties which are computed as a difference between two states). A departure function gives the difference between the real state, at a finite volume or non-zero pressure and temperature, and the ideal state, usually at zero pressure or infinite volume and temperature.
For example, to evaluate enthalpy change between two points h(v1,T1) and h(v2,T2) we first compute the enthalpy departure function between volume v1 and infinite volume at T = T1, then add to that the ideal gas enthalpy change due to the temperature change from T1toT2, then subtract the departure function value between v2 and infinite volume.
Departure functions are computed by integrating a function which depends on an equation of state and its derivative.
Departure functions for Peng–Robinson equation of state[edit]
The Peng–Robinson equation of state relates the three interdependent state properties pressure P, temperature T, and molar volume Vm. From the state properties (P, Vm, T), one may compute the departure function for enthalpy per mole (denoted h) and entropy per mole (s):[2]
Typically, one knows two of the three state properties (P, Vm, T), and must compute the third directly from the equation of state under consideration. To calculate the third state property, it is necessary to know three constants for the species at hand: the critical temperatureTc, critical pressurePc, and the acentric factorω. But once these constants are known, it is possible to evaluate all of the above expressions and hence determine the enthalpy and entropy departures.
![{\displaystyle {\begin{aligned}{\frac {H^{\mathrm {ig} }-H}{RT}}&=\int _{V}^{\infty }\left[T\left({\frac {\partial Z}{\partial T}}\right)_{V}\right]{\frac {dV}{V}}+1-Z\\[2ex]{\frac {S^{\mathrm {ig} }-S}{R}}&=\int _{V}^{\infty }\left[T\left({\frac {\partial Z}{\partial T}}\right)_{V}-1+Z\right]{\frac {dV}{V}}-\ln Z\\[2ex]{\frac {G^{\mathrm {ig} }-G}{RT}}&=\int _{V}^{\infty }(1-Z){\frac {dV}{V}}+\ln Z+1-Z\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/36c1b804e6f2a2e03c65b788e9efdb0679801fa0)
![{\displaystyle {\begin{aligned}h_{T,P}-h_{T,P}^{\mathrm {ideal} }&=RT_{C}\left[T_{r}(Z-1)-2.078(1+\kappa ){\sqrt {\alpha }}\ln \left({\frac {Z+2.414B}{Z-0.414B}}\right)\right]\\[1.5ex]s_{T,P}-s_{T,P}^{\mathrm {ideal} }&=R\left[\ln(Z-B)-2.078\kappa \left({\frac {1+\kappa }{\sqrt {T_{r}}}}-\kappa \right)\ln \left({\frac {Z+2.414B}{Z-0.414B}}\right)\right]\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7e57561471341070dd3e7b2ac32e915248e1b11b)
