The derivatives of the general NRTL equations are very useful for VLE calculations, but are incredibly cumbersome to perform by hand in the general form.
General NRTL equations
The general equation for
for species
in a mixture of
components is[1]:
![{\displaystyle \ln(\gamma _{i})={\frac {\displaystyle \sum _{j=1}^{n}{x_{j}\tau _{ji}G_{ji}}}{\displaystyle \sum _{k=1}^{n}{x_{k}G_{ki}}}}+\sum _{j=1}^{n}{\frac {x_{j}G_{ij}}{\displaystyle \sum _{k=1}^{n}{x_{k}G_{kj}}}}{\left({\tau _{ij}-{\frac {\displaystyle \sum _{m=1}^{n}{x_{m}\tau _{mj}G_{mj}}}{\displaystyle \sum _{k=1}^{n}{x_{k}G_{kj}}}}}\right)}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d4330523a3603ef7bd15b498c4dbe9f1bbcc17f5) | | (1.1) |
with
![{\displaystyle G_{ij}={\text{exp}}\left({-\alpha _{ij}\tau _{ij}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d0f77ca288885b6041f6e647312ab06fc1d29bd5) | | (1.2) |
![{\displaystyle \alpha _{ij}=\alpha _{ij_{0}}+\alpha _{ij_{1}}T}](https://wikimedia.org/api/rest_v1/media/math/render/svg/324d37d0f837df04b11a7aebb0059e4352978d80) | | (1.3) |
![{\displaystyle \tau _{ij}=A_{ij}+{\frac {B_{ij}}{T}}+{\frac {C_{ij}}{T^{2}}}+D_{ij}\ln {\left({T}\right)}+E_{ij}T^{F_{ij}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a8efec6b7fcaf9c194cdd741b4e3176d782e8d78) | | (1.4) |
For derivative calculations, it becomes convenient to further compartmentalize the general NRTL equation by abstracting the summation terms. While this does introduce an additional substitution for chain rule differentiation, it does make the final equation more readable.
![{\displaystyle \ln(\gamma _{i})={\frac {S_{1_{ij}}}{S_{2_{ik}}}}+\sum _{j=1}^{n}{\frac {x_{j}G_{ij}}{S_{2_{jk}}}}{\left({\tau _{ij}-{\frac {S_{1_{jm}}}{S_{2_{jk}}}}}\right)}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bb6da010088d19416cf192d21b77e06ddcc53862) | | (1.5) |
with
Sum type 1:
![{\displaystyle S_{1_{ij}}=\displaystyle \sum _{j=1}^{n}{x_{j}\tau _{ji}G_{ji}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/37909e15d70718265549f705320bab3df8fbc393) | | (1.6) |
Sum type 2:
![{\displaystyle S_{2_{ij}}=\displaystyle \sum _{j=1}^{n}{x_{j}G_{ji}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2e9afdb559846a29db0bdbfe69905c76e668bd38) | | (1.7) |
Derivatives of system variables
The system variables are intrinsic to a system being observed or predicted. These include several directly measurable variables and many indirectly measurable variables. Only the directly measurable system variables need to be calculated for the NRTL model.
- Directly measurable system variables:
Temperature
These, along with the composition derivatives, are the primary derivatives of interest, the reason for which will become obvious when reading the pressure derivatives section.
Adjustable parameters
Since the adjustable parameters are treated as constants during typical VLE calculations, their derivatives are straight forward degenerate solutions:
Non-randomness
![{\displaystyle {\begin{matrix}\displaystyle \left({\frac {\partial \alpha _{ij_{0}}}{\partial T}}\right)_{P,{\vec {N}}}=0&\displaystyle \left({\frac {\partial \alpha _{ij_{1}}}{\partial T}}\right)_{P,{\vec {N}}}=0\end{matrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c1a29231d08c43444ba85c753754e3dfe79c60ad) | | (2.1-1) |
Interaction
![{\displaystyle {\begin{matrix}\displaystyle \left({\frac {\partial A_{ij}}{\partial T}}\right)_{P,{\vec {N}}}=0&\displaystyle \left({\frac {\partial B_{ij}}{\partial T}}\right)_{P,{\vec {N}}}=0&\displaystyle \left({\frac {\partial C_{ij}}{\partial T}}\right)_{P,{\vec {N}}}=0&\displaystyle \left({\frac {\partial D_{ij}}{\partial T}}\right)_{P,{\vec {N}}}=0&\displaystyle \left({\frac {\partial E_{ij}}{\partial T}}\right)_{P,{\vec {N}}}=0&\displaystyle \left({\frac {\partial F_{ij}}{\partial T}}\right)_{P,{\vec {N}}}=0\end{matrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0f1849b51c55d90d8445d54847e36d92ea8f5a82) | | (2.1.1-1) |
Non-randomness term
The non randomness term has a similarly simple series of temperature derivatives, being linear with respect to temperature.
First order:
![{\displaystyle \left({\frac {\partial \alpha _{ij}}{\partial T}}\right)_{P,{\vec {N}}}=\alpha _{ij_{1}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/92519a8eef4560b69cff79e2c5b1604f1ecb78e5) | | (2.1.2-1) |
Second and higher order:
![{\displaystyle \left({\frac {\partial ^{n}\alpha _{ij}}{\partial T^{n}}}\right)_{P,{\vec {N}}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cfd1d90b767fe5521ef77cc3c9ab7b36912170d6) | | (2.1.2-2) |
Interaction term
The interaction term can be conveniently expressed in terms of the derivative order, which becomes more or less obvious after a couple of iterations:
First order:
![{\displaystyle \left({\frac {\partial \tau _{ij}}{\partial T}}\right)_{P,{\vec {N}}}=-{\frac {B_{ij}}{T^{2}}}-{\frac {2C_{ij}}{T^{3}}}+{\frac {D_{ij}}{T}}+E_{ij}F_{ij}T^{F_{ij}-1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7b6264bbb400c184da190017a95a5736282fe937) | | (2.1.3-1) |
Second order:
![{\displaystyle \left({\frac {\partial ^{2}\tau _{ij}}{\partial T^{2}}}\right)_{P,{\vec {N}}}={\frac {2B_{ij}}{T^{3}}}+{\frac {6C_{ij}}{T^{4}}}-{\frac {D_{ij}}{T^{2}}}+E_{ij}F_{ij}\left(F_{ij}-1\right)T^{F_{ij}-2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/726459b5f43b630cda4dc8e57fc35b98642573e6) | | (2.1.3-2) |
Higher order: {{NumBlk||
|2.1.3.-3}
First order
Interaction energy term
![{\displaystyle \left({\frac {\partial G_{ij}}{\partial T}}\right)_{P,{\vec {N}}}=-\left(\alpha _{ij}\left({\frac {\partial \tau _{ij}}{\partial T}}\right)_{P,{\vec {N}}}+\tau _{ij}\left({\frac {\partial \alpha _{ij}}{\partial T}}\right)_{P,{\vec {N}}}\right){\text{exp}}\left({-\alpha _{ij}\tau _{ij}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ffe79bc8f90137adfcdac8e0ecb3eb777239d69b) | | (2.1.4.1-1) |
![{\displaystyle \left({\frac {\partial G_{ij}}{\partial T}}\right)_{P,{\vec {N}}}=-\left(\alpha _{ij}\left({\frac {\partial \tau _{ij}}{\partial T}}\right)_{P,{\vec {N}}}+\tau _{ij}\left({\frac {\partial \alpha _{ij}}{\partial T}}\right)_{P,{\vec {N}}}\right)G_{ij}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4079099a85961961cea823aef1b4a9d9a995c207) | | (2.1.4.1-2) |
Sum type 1
![{\displaystyle \left({\frac {\partial S_{1_{ij}}}{\partial T}}\right)_{P,{\vec {N}}}=\displaystyle \sum _{j=1}^{n}{x_{j}\left(\tau _{ji}\left({\frac {\partial G_{ji}}{\partial T}}\right)_{P,{\vec {N}}}+G_{ji}\left({\frac {\partial \tau _{ji}}{\partial T}}\right)_{P,{\vec {N}}}\right)}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3b90c02a05cfaba5ea3b0275eb05c5c239f983ef) | | (2.1.4.2-1) |
Sum type 2
![{\displaystyle \left({\frac {\partial S_{1_{ij}}}{\partial T}}\right)_{P,{\vec {N}}}=\displaystyle \sum _{j=1}^{n}{x_{j}\left({\frac {\partial G_{ji}}{\partial T}}\right)_{P,{\vec {N}}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d530ac183f08370e8178f20b90e1106859121d9f) | | (2.1.4.3-1) |
Main equation
![{\displaystyle \left({\frac {\partial \ln {\gamma _{i}}}{\partial T}}\right)_{P,{\vec {N}}}={\frac {S_{2_{ik}}\left({\frac {\partial S_{1_{ij}}}{\partial T}}\right)_{T,{\vec {N}}}-S_{1_{ij}}\left({\frac {\partial S_{2_{ik}}}{\partial T}}\right)_{T,{\vec {N}}}}{{S_{2_{ik}}}^{2}}}+\sum _{j=1}^{n}{x_{j}{\frac {S_{2_{jk}}\left(G_{ij}\left(\left({\frac {\partial \tau _{ij}}{\partial T}}\right)_{P,{\vec {N}}}-{\frac {S_{2_{jk}}\left({\frac {\partial S_{1_{jm}}}{\partial T}}\right)_{P,{\vec {N}}}-S_{1_{jm}}\left({\frac {\partial S_{2_{jk}}}{\partial T}}\right)_{P,{\vec {N}}}}{{S_{2_{jk}}}^{2}}}\right)+\left(\tau _{ij}-{\frac {S_{1_{jm}}}{S_{2_{jk}}}}\right)\left({\frac {\partial G_{ij}}{\partial T}}\right)_{P,{\vec {N}}}\right)+G_{ij}\left(\tau _{ij}-{\frac {S_{1_{jm}}}{S_{2_{jk}}}}\right)\left({\frac {\partial S_{2_{jk}}}{\partial T}}\right)_{P,{\vec {N}}}}{{S_{2_{jk}}}^{2}}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6fd5d84f7ff6ac6ba449fc4379cf56be6eb5f614) | | (2.1.4.4-1) |
Second order
Third order
Pressure
These are the least interesting, though still very important and useful, derivatives. Since the pressure derivatives are evaluated at constant temperature and composition and none of the terms involved are explicit functions of pressure, all of the derivative expressions are essentially derivatives of a constant scalar term, which is zero. All of the partial derivatives are listed below for sake of complete coverage of the topic.
Adjustable parameters
The pressure derivatives for the adjustable parameters are analogous to their respective temperature derivatives.
![{\displaystyle \left({\frac {\partial ^{n}\phi }{\partial P^{n}}}\right)_{T,{\vec {N}}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1a44b17f7b4f9020c07a761fcf342d8a1c65eb35) | | (2.2.1-1) |
where
is an adjustable parameter.
Non-randomness term
Unlike the temperature derivatives, the pressure derivatives of the non-randomness term are all zero.
![{\displaystyle \left({\frac {\partial ^{n}\alpha _{ij}}{\partial P^{n}}}\right)_{T,{\vec {N}}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/95e69c51cc63f0dfd3fc38d83e14ac3b9dadea3f) | | (2.2.2-1) |
Interaction term
![{\displaystyle \left({\frac {\partial ^{n}\tau _{ij}}{\partial P^{n}}}\right)_{T,{\vec {N}}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bad1091325fbfdc3986fac500c3e4f4a7b42568b) | | (2.2.3-1) |
Interaction energy term
![{\displaystyle \left({\frac {\partial ^{n}G_{ij}}{\partial P^{n}}}\right)_{T,{\vec {N}}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/52f251b8611a24cc2831cfa72a1b8ef8420183e2) | | (2.2.4-1) |
Sum type 1
![{\displaystyle \left({\frac {\partial ^{n}S_{1_{ij}}}{\partial P^{n}}}\right)_{T,{\vec {N}}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/798ae362e298887de329148fb2d17d9ba52964a6) | | (2.2.5-1) |
Sum type 2
![{\displaystyle \left({\frac {\partial ^{n}S_{2_{ij}}}{\partial P^{n}}}\right)_{T,{\vec {N}}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6914ab62d94969a4f53e049ff5cc96d0d81d1888) | | (2.2.1-6) |
Main equation
This is the most important derivative resulting from the pressure derivatives due to its use in simplifying calculations.
![{\displaystyle \left({\frac {\partial ^{n}\ln {\gamma _{i}}}{\partial P^{n}}}\right)_{T,{\vec {N}}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b88b252a72cf0d40c93ba694287db386d5e2d0bd) | | (2.2.7-1) |
Composition
Adjustable parameters
![{\displaystyle \left({\frac {\partial ^{n}N\phi }{\partial {N_{i}}^{n}}}\right)_{T,P,N_{j\neq i}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/509c37a3f86d4334dc5bafa63c191a85865152ca) | | (2.3.1-1) |
where
is an adjustable parameter.
Non-randomness term
![{\displaystyle \left({\frac {\partial ^{n}N\alpha _{ij}}{\partial {N_{k}}^{n}}}\right)_{T,P,N_{m\neq k}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7c96483317caa4e83fe467e5dfbfea3b1930fb34) | | (2.3.2-1) |
Interaction term
![{\displaystyle \left({\frac {\partial ^{n}N\tau _{ij}}{\partial {N_{k}}^{n}}}\right)_{T,P,N_{m\neq k}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a632e0060b63c2e1ed86fa559e33bebee98f4659) | | (2.3.3-1) |
Interaction energy term
![{\displaystyle \left({\frac {\partial ^{n}NG_{ij}}{\partial {N_{k}}^{n}}}\right)_{T,P,N_{m\neq k}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f1a4ab018cd8dc67b8440e5f6faa991d4dea3fcb) | | (2.3.4-1) |
Sum type 1
![{\displaystyle \left({\frac {\partial NS_{1_{ij}}}{\partial N_{k}}}\right)_{T,P,N_{m\neq k}}=\tau _{ki}G_{ki}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/91635e94da93572674a7f3a539a7391085c78dd6) | | (2.3.5-1) |
Sum type 2
![{\displaystyle \left({\frac {\partial NS_{1_{ij}}}{\partial N_{k}}}\right)_{T,P,N_{m\neq k}}=G_{ki}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c6142ab310d7113b479f0438b52748209ead7489) | | (2.3.6-1) |
Main equation
First order
Second order
Third order
Derivatives of adjustable parameters
The adjustable parameters are fit to experimental data for known compositions and are coefficients for temperature dependent parameters of the general NRTL equation. These parameters are considered variables during fit optimization calculations only. During the calculations for vapor-liquid equilibrium, these parameters are considered as constants.
- Non-randomness parameters:
, scalar
, linearly proportional
- Interaction parameters:
, scalar
, inversely proportional
, inverse square proportional
, logarithmically proportional
, power proportional
, power term exponent
Non-randomness parameters
Scalar
Linear
Interaction parameters
Scalar
Inversely proportional
Inverse square proportional
Logarithmically proportional
Power proportional
Power term exponent
References