Next Article in Journal
Study on Treatment of Tiny Pollution Water with PAC-HUM System in Kuitun River
Next Article in Special Issue
Permeation of Ternary Mixture Containing H2S, CO2 and CH4 in Aquivion® Perfluorosulfonic Acid (PFSA) Ionomer Membranes
Previous Article in Journal
Magnesium-Doped Sr2(Fe,Mo)O6−δ Double Perovskites with Excellent Redox Stability as Stable Electrode Materials for Symmetrical Solid Oxide Fuel Cells
Previous Article in Special Issue
Purification of Adeno-Associated Virus (AAV) Serotype 2 from Spodoptera frugiperda (Sf9) Lysate by Chromatographic Nonwoven Membranes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Rigid Adsorbent Lattice Fluid Model: Thermodynamic Consistency and Relationship to the Real Adsorbed Solution Theory

School of Engineering, University of Edinburgh, The King’s Buildings, Mayfield Road, Edinburgh EH9 3FB, UK
Membranes 2022, 12(10), 1009; https://doi.org/10.3390/membranes12101009
Submission received: 30 August 2022 / Revised: 13 October 2022 / Accepted: 14 October 2022 / Published: 18 October 2022

Abstract

:
The Rigid Adsorbent Lattice Fluid model has been shown to comply with all the requirements for thermodynamic consistency in the case of an adsorbent that does not undergo structural changes. This is achieved by introducing a correction to the reduced density function that multiplies the combinatorial term. A procedure to calculate the predicted adsorbed mixture activity coefficients has been presented that allows the production of excess Gibbs energy plots at a constant reduced grand potential. The predicted nonideality is structurally consistent with the Non-Ideal Adsorbed Solution Theory of Myers in terms of both its dependence on concentration and reduced grand potential. The ability to generate excess Gibbs energy values allows linking the new Rigid Adsorbent Lattice Fluid model to the traditional Real Adsorbed Solution Theory providing an alternative approach to predicting multicomponent adsorption based solely on pure component data.

1. Introduction

A new thermodynamic framework for pure and mixed gas adsorption, the Rigid Adsorbent Lattice Fluid (RALF) model was presented in [1], adapting to a non-distributing solid the Non-Equilibrium Lattice Fluid (NELF) approach for polymers developed and applied extensively by Professor Sarti and his colleagues at the University of Bologna [2,3,4,5], including a recent extensive review [6]. Since its introduction, which considered the adsorption of gases and vapours in the zeolite silicalite, the RALF model has been used to correlate adsorption in breathing metal–organic frameworks [7]; to demonstrate that it allows predicting different types of isotherms depending on the physical parameters of the adsorbates [8]; and more recently, it was shown to capture the behaviour of ZSM-25 in dynamic experiments [9] as well as water adsorption with stepped behaviour [10]. The emphasis of these publications [1,7,8,9,10] has been to investigate the applicability of RALF to a range of challenging adsorption systems as well as test the prediction capability for multicomponent adsorption, thus focusing primarily on the phase equilibria of the adsorbates.
The purpose of the present contribution was to study in greater detail the case of the frozen solid, where the adsorbent does not undergo structural changes, and focus on:
1. The requirements for thermodynamic consistency, which as will be shown lead to a small modification to the adsorbed phase expression for multicomponent mixtures.
2. Demonstrate how it is possible to determine the adsorbed phase excess Gibbs energy from RALF, thus providing the means to bridge this approach to the classical Real Adsorbed Solution Theory (RAST).
It is hoped that this contribution will be of interest to both adsorption and membrane communities that have developed thermodynamic approaches for phase equilibria between the gas (and liquid) and solid phases with limited attempts at reconciling the apparently different approaches. One would hope as well that this will not lead to a sense of bewilderment in the reader [6] as in both cases the assumption made is that the solid phase (membrane or adsorbent) is considered to be a lattice fluid.

2. RALF Model for a Frozen Adsorbent

A basic assumption in the thermodynamics of adsorption is that the adsorbent does not undergo structural changes [11] and in the RALF framework this corresponds to the frozen solid limit [1]. The volume of the system, V , is assumed to be the volume occupied by the solid, including the micropores, V S , [12] and
V = V S = m S ρ S = j m j ρ
where ρ S is the density of the solid including the micropores, ρ is the density of the adsorbed phase, including the adsorbates. m S is the mass of the solid, while m j are the masses of all components, including the solid. In the frozen solid limit V S is a constant [1].
The RALF model was derived from the Lattice Fluid model of Sanchez and Lacombe [13,14,15] and the following expression for the residual Gibbs energy for fluid was obtained in [1]
G R ( T , P , N ) R T = r N [ ρ ˜ T ˜ + ( 1 ρ ˜ ) ln ( 1 ρ ˜ ) ρ ˜ + 1 ] + N ρ ˜ j x j ln ( ϕ j x j ) + N ( z 1 ln z )
with → ρ ˜ = ρ ρ * → → T ˜ = T T * → → z = P V N R T = r P ˜ ρ ˜ T ˜ → → P ˜ = P P * . where ρ * is the characteristic density; T * is the characteristic temperature, related to an energy; P * is the characteristic pressure, related to an energy density; ϕ j are the volume fractions of each specie; r j are the number of lattice sites occupied by each specie while r is the same for the mixture.
The following relationship links the characteristic parameters
P * v * = R T *
Equation (2) has a + 1 in the first parenthesis [16] included so that the first term cancels in the limit of zero pressure. The combinatorial term is pre-multiplied by ρ ˜ to ensure that at zero pressure the ideal gas state is recovered [1]. While these are minor changes from the NELF formulation, the basis for the RALF model is to define Equation (2) for the adsorbed phase. First the volume occupied by adsorbed molecules at close packing, v i * , is corrected to take into account confinement constraints.
v i A * = ( 1 + ξ i A ) v i *       and       ρ i A * = ρ i * 1 + ξ i A
The second important modification is that the combinatorial contribution needs to reflect the fact that due to the rigid nature of the solid it does not contribute to the combinatorial term, but has the effect of reducing the volume available to the molecules in the lattice. Therefore, in the original RALF framework, the residual Gibbs energy of the adsorbed phase was specified as
G A R ( T , P , N ) R T = r N [ ρ ˜ T ˜ + ( 1 ρ ˜ ) ln ( 1 ρ ˜ ) ρ ˜ + 1 ] + N ρ ˜ i x i ln ϕ i ( 1 ϕ S ) x i + N ( z 1 ln z )
Here, ( 1 ϕ S ) is highlighted to identify the additional deviation from the case of a polymer.
In the model sums over the index i include only the adsorbates, while sums over j include also the solid, which is the last component. In the LF model component, 0 is used for the vacancies.
What is not immediately apparent is the fact that for multicomponent mixtures Equation (5) does not reduce to the Ideal Adsorbed Solution Theory (IAST) [17], which is one of the requirements for thermodynamic consistency [18]. The issue comes from the combinatorial term, which does not cancel as in the limit of zero pressure the reduced density is ρ ˜ = ρ ˜ S = ρ S ρ S * and not zero.
To ensure thermodynamic consistency the residual Gibbs energy of the adsorbed phase has to be corrected so that the combinatorial term is zero at zero pressure and a suitable expression is:
G A R ( T , P , N ) R T = r N [ ρ ˜ T ˜ + ( 1 ρ ˜ ) ln ( 1 ρ ˜ ) ρ ˜ + 1 ] + N [ ρ ˜ ρ ˜ S ( 1 + ln ρ ˜ ρ ˜ S ) ] i x i ln ϕ i ( 1 ϕ S ) x i + + N ( z 1 ln z )
The term multiplying the combinatorial part ensures that both the residual energy and the chemical potential of the adsorbates have the correct limit at zero pressure since
lim P 0 ρ ˜ ρ ˜ S ( 1 + ln ρ ˜ ρ ˜ S ) = 0
This modification does not affect the single adsorbate case as the combinatorial term for a pure adsorbate is zero given that for a single adsorbate ϕ 1 = ( 1 ϕ S ) . For mixtures at high adsorbed concentrations the correction is small, thus all numerical results presented in the original RALF model [1] are not affected by this correction.
The mathematical structure of the correction term leads to an equation of state term given by
z A E o S 1 = r [ ρ ˜ T ˜ ln ( 1 ρ ˜ ) ρ ˜ 1 ] + ( ρ ˜ ρ ˜ S ) i x i ln ϕ i x i ( 1 ϕ S )  
for which the combinatorial term of the adsorbed phase cancels in the limit of zero pressure.
To apply Equation (6), mixing rules for the model parameters need to be specified. These are [1]
r i v * = r i 0 v i A *
which conserves the close-packed molecular volume of each component;
j r j 0 N j = r N
which preserves the number of pair interactions in the close-packed state; and
P * = j k ϕ j ϕ k P j k *
with → P j k * = P k j * = ( 1 κ k j ) P k * P j * → and κ k k = 0 . which is the classical quadratic mixing rule [19].
The reduced residual chemical potentials can be obtained by derivation
μ k R R T = 1 R T ( G A R N k ) T , P , N j k = ln φ k         μ S R m R T = 1 R T ( G A R m S ) T , P , N j S
where φ k is the fugacity coefficient of adsorbate k . Note that an asymmetric convention is used with moles for the adsorbates and mass for the adsorbent.
The expressions for the reduced residual chemical potential of the adsorbates in the case of a frozen solid are:
μ k R R T = ρ ˜ T ˜ r k ( 2 j ϕ j P k j * P * 1 ) + [ ( 1 ρ ˜ ) ln ( 1 ρ ˜ ) ρ ˜ + 1 ] r k 0 + r k r ( z A E o S 1 ) + [ ρ ˜ ρ ˜ S ( 1 + ln ρ ˜ ρ ˜ S ) ] ( ln r k r ( 1 ϕ S ) + 1 r k r ( 1 ϕ S ) ) ln z
Equation (13) is all that is needed for phase equilibrium calculations but to establish the connection between RALF and RAST what is important is the expression for the reduced residual chemical potential of the solid, which for a frozen solid is given by:
μ S R m R T = ρ ˜ ϕ S ρ S v S * [ ( 1 ρ ˜ ) ln ( 1 ρ ˜ ) ρ ˜ + 1 ] ρ ˜ T ˜ [ 2 j ϕ j P S j * P * 1 ] ρ ˜ ϕ S ρ S v * + n ( ϕ S 1 ) ( z A E o S 1 ) n + P R T ρ S
where n = N m S is the total absolute adsorbed amount in moles/kg of adsorbent and P R T ρ S = n z .
From Equation (14) it is possible to define also the reference state for the solid as that of the solid phase without adsorbates at the same pressure of the system [20] as in this limit ϕ S = 1 ; ρ ˜ = ρ ˜ S ; v * = v S * and N = 0 , therefore
μ S 0 R m R T = 1 v S * ρ S * [ ( 1 ρ ˜ S ) ln ( 1 ρ ˜ S ) ρ ˜ S + 1 ρ ˜ S T ˜ S ] + P R T ρ S
These last two expressions lead to the reduced grand potential [20]
ψ = μ S 0 R m R T μ S R m R T = 1 v S * ρ S * [ ( 1 ρ ˜ S ) ln ( 1 ρ ˜ S ) ρ ˜ S + 1 ρ ˜ S T ˜ S ] ρ ˜ ϕ S ρ S v S * [ ( 1 ρ ˜ ) ln ( 1 ρ ˜ ) ρ ˜ + 1 ] + ρ ˜ T ˜ [ 2 j ϕ j P S j * P * 1 ] ρ ˜ ϕ S ρ S v * n ( ϕ S 1 ) ( z A E o S 1 ) + n
which is a fundamental quantity [11] in the thermodynamics of the adsorbed phase as it allows to define the reference state relative to which activity coefficients can be obtained. One should note that Equation (15) is only a function of composition and the total adsorbed amount and does not depend explicitly on the system pressure.
For completeness and to be able to reproduce all the terms needed in the calculations Table 1 includes all the steps and expressions required to calculate the reduced residual potentials of the adsorbates and the reduced grand potential in the RALF model for the frozen solid. In what follows the numerical calculations will be based on the parameters for silicalite obtained following the procedure in [1] from experimental data including Henry law constants, adsorption energies and pure component isotherms for several adsorbates including normal alkanes [21,22,23].

3. Thermodynamic Consistency of an Adsorbed Phase

Talu and Myers [18] cover in some detail the requirements for consistency of adsorption experimental data. The requirements for experimental data apply equally to the results from a model and are the following [18]:
  • Single-gas adsorption isotherms should reduce to Henry’s law at the limit of zero pressure.
  • Multicomponent isotherms should display continuity with single-gas isotherms.
  • At fixed temperature and pressure, thermodynamically consistent x-y diagrams should intersect the predictions from the IAST at least once.
  • In the limit of zero pressure, the IAST should be obtained.
  • Activity coefficients in the adsorbed phase are a function of composition and the reduced grand potential.
The last condition requires a method to determine the activity coefficients of the adsorbed phase from the RALF model and this will be covered in a separate section.
Condition 1 is fulfilled as the RALF model reduces to a Henry law constant at low pressure [1] or
n 1 = N 1 m S = K P ( T ) P
At infinite dilution and in the limit of low pressure, for a single adsorbate, it is possible to obtain [1]
ln K P 1 = ln 1   k g ρ S R T + ρ ˜ S R T r 1 0 v 1 * 2 P 1 S * [ ( 1 ρ ˜ S ) ln ( 1 ρ ˜ S ) ρ ˜ S + 1 ] r 1 0 + r 1 0 v 1 * v S * [ ln ( 1 ρ ˜ S ) ρ ˜ S + 1 ]
Given that the RALF model is formulated in terms of thermodynamically consistent mixing rules, there is continuity when one component in the mixture becomes very dilute, i.e., a ternary mixture reduces correctly to the binary mixtures and a binary mixture reduces correctly to the single component isotherms, thus condition 2 is automatically fulfilled.
Condition 3 stems from the integral consistency test and the Gibbs–Duhem equation which requires that at least at one point the activity coefficients of the two components are the same [18]. As a result, the binary selectivity is equal to that of the IAST at least at one point. Figure 1 shows the difference between the gas phase mole fraction calculated from RALF and IAST for the binary ethane/n-butane in silicalite using the RALF parameters reported in [1] at 300 K and 350 kPa. The two models cross at x = 0.574 . The IAST prediction was obtained using a nested-loop calculation [24] where the reduced grand potential of the ideal mixture, ψ I A S T , requires the solution of
1 i = 1 N c f i 0 ( ψ I A S T ) x i φ i F P = 0
where φ i F is the fugacity coefficient of each component in the fluid phase and f i 0 is the reference fugacity calculated from the pure component isotherm. In the RALF model, this last term requires the intermediate calculation of the adsorbed amount for each pure adsorbate from the knowledge of ψ I A S T . The multicomponent calculations for the RALF model require the solution of the equilibrium relationships given in [1]. For convenience, the fluid phase was described using the Sanchez–Lacombe equation of state modified to be consistent at zero pressure [1].
Condition 4 is fulfilled and this can be demonstrated considering the limit as pressure goes to zero of Equation (13).
μ k R 0 R T = ρ ˜ S R T r k 0 v k * 2 P k S * + [ ( 1 ρ ˜ S ) ln ( 1 ρ ˜ S ) ρ ˜ S + 1 ] r k 0 r k 0 v k * v S * [ ln ( 1 ρ ˜ S ) ρ ˜ S + 1 ] ln z = ln ( 1   k g i n i K P k P )
When this is combined with the equilibrium relationship one obtains the correct limiting behaviour of an ideal mixture in the Henry law region
n k = K P k P y k
As a final internal thermodynamic consistency check, one can also verify that the reduced grand potential obtained from Equation (16) is consistent with the Gibbs adsorption isotherm [11,25]:
ψ = 0 P i n i A d l n f i
It is relatively straightforward to test this relationship for a single component isotherm and this is a useful check especially in verifying the correct implementation of the RALF model.

4. RALF and RAST

In this section, a procedure to obtain the activity coefficients of the adsorbed components from the RALF model will be derived. As the RALF model is a combination of an energy term of a regular solution with the combinatorial term of an athermal solution [1], one would expect that a simple excess Gibbs energy expression should suffice. The picture in an adsorbed phase is complicated by the fact that the excess Gibbs energy will be a function also of the reduced grand potential. Myers has shown that to a very good approximation g E x models from liquid phase correlations [19] can be used for the concentration dependence of adsorbed phase provided that an additional term is included [26,27,28]
g E x A R T = ( 1 e C ψ ) g E x R T = i x i ln γ i A
One can see that the additional term ensures that at low pressure the IAST is recovered in compliance with condition 4.
To calculate the activity coefficients of the RALF model one has to consider the definition of the fugacity of the adsorbed phase from RAST, which requires the reference state of the pure components at the same reduced grand potential of the mixture. Therefore, the two adsorbed phase fugacities can be equated
γ k A φ k A 0 P k 0 = φ k A P
This expression can be rearranged to a more convenient form
ln γ k A = μ k R R T + ln z ln n R T ρ S f k A 0
With the pure component reference fugacity f k A 0 = φ k A 0 P k 0 . Equation (25) depends only on the adsorbed phase concentrations and the total adsorbed amount as is the case for the reduced grand potential.
It is now possible to calculate the excess Gibbs energy by assigning the adsorbed phase mole fractions and the reduced grand potential. The procedure is as follows
6.
Determine the total adsorbed amount from Equation (16).
7.
Repeat this for the pure component case in order to determine the pure component adsorbed amounts at the same reduced grand potential of the mixture.
8.
From the pure component isotherm, calculate the reference pressure and fugacity corresponding to the adsorbed amounts obtained in step 2.
9.
Calculate the activity coefficients of all components from Equation (25).
Steps 1 and 2 are fairly easy to implement as each corresponds to the solution of a single nonlinear equation with a function that increases monotonically with the number of moles adsorbed. The physically valid range of the adsorbed amounts will lie between 0 and the finite saturation capacity values for the mixture and for each pure component, respectively. The saturation capacities are obtained by inspection and from the definition of the close-packed state:
n S a t i x i M w i ρ i A * = 1 ρ S ( 1 ρ S ρ S * )
For a single adsorbent [1]
n 1 S a t = 1 M w 1 ρ 1 A * ρ S ( 1 ρ S ρ S * )
Step 3 is slightly more complicated as the adsorption isotherm is an implicit function of pressure, but the adsorbed amount increases monotonically with pressure and this makes the numerical solution robust.
Spanning the compositions at a fixed reduced grand potential is the most convenient way to see the composition dependence of the excess Gibbs energy of the adsorbed phase. An alternative is to assign increasing values of the total adsorbed amount and determine at a fixed mole fraction the dependence with respect to the reduced grand potential. In this case, Step 1 is replaced by direct computation of the reduced grand potential of the mixture having specified the total amount adsorbed and composition of the adsorbates.
While a detailed investigation of all possible cases that can be generated using the RALF model including the use of the mixture parameters is beyond the scope of this study, here it is possible to investigate the behaviour of the binary mixtures of alkanes in silicalite, which should be represented well by the Non-Ideal Adsorbed Solution Theory, NIAST, of Myers [28]. NIAST is a predictive model based on the simplest expression for activity coefficients, which at a fixed temperature can be written as
g E x A R T = 4 δ 12 x 1 x 2 ( 1 e C 12 ψ )
With
δ 12 = | ψ 1 0.5 n 1 S a t ψ 2 0.5 n 2 S a t |   ;   C 12 = e δ 12 n 12 S a t   ;   1 n 12 S a t = 1 2 ( 1 n 1 S a t + 1 n 2 S a t )
ψ k 0.5 is the reduced grand potential of pure component k at an adsorbed concentration of 0.5 n k S a t . Therefore, the NIAST parameters can be calculated from the pure component isotherms of the RALF model providing an alternative means of predicting the binary mixtures.
As the deviations from ideality are mild in the alkanes silicalite systems considered, calculations were carried out on the methane/n-butane mixture which shows the largest deviations. Figure 2 shows the excess Gibbs energy predicted at three mole fractions (0.25, 0.5 and 0.75) as a function of the reduced grand potential of the mixture at 300 K.
As can be seen from Figure 2 the RALF model shows a more complex behaviour than the NIAST approach, as the excess Gibbs energy includes some asymmetry that is more pronounced as the adsorbed amount increases (at the same reduced grand potential the curves at 0.25 and 0.75 should overlap in the case of symmetry), and there is also a small region at low reduced grand potential (low pressure and adsorbed amount region) where positive deviations are predicted. Overall, though negative deviations are predicted where typical experimental measurements are carried out. To see more clearly the exponential dependence with respect to the reduced grand potential, Figure 3 shows the normalised excess Gibbs energy for a mole fraction of 0.5.
Structurally the RALF model contains the correct dependence on the reduced grand potential but the predictions are of smaller deviations from ideality at lower pressures compared to the NIAST formulation.
Figure 4 shows the comparison at high adsorbed amounts for ψ = 150 mol/kg. It is very interesting to note that both NIAST and RALF models predict the same infinite dilution activity coefficient for n-butane in methane (region close to x = 1 ), while RALF predicts a slightly smaller infinite dilution activity coefficient at the other limit, with a resulting asymmetric excess Gibbs energy that can be matched well with a two-parameter Margules binary activity coefficient expression [29] (dotted line). Including the effect of the reduced grand potential the approximate expression is given by
g E x A R T = ( A 21 x 1 + A 12 x 2 ) x 1 x 2 ( 1 e C ψ )
Matching the excess Gibbs energy at high reduced grand potential A 21 = 0.156 and A 12 = 0.229 .
At lower reduced grand potential values, the excess Gibbs energy predicted from RALF becomes nearly symmetric as shown in Figure 5 for ψ = 10 mol/kg. The dotted line close to the RALF data points is calculated using Equation (29) and the deviation is small.
Structurally RALF is a rather complex model that produces the correct mixture behaviour for the systems considered with similar quantitative results when compared to the NIAST of Myers [28]. The key advantage of RALF is that it is formulated as a multicomponent model and is not limited to binary systems [1].

5. Conclusions

A correction to the reduced density function that multiplies the combinatorial term in RALF was shown to lead to a formulation that passes all thermodynamic consistency checks. As the Sanchez–Lacombe equation of state was derived starting from the close-packed limit it was important to carry out these checks and identify the incorrect zero pressure limit of the original formulation.
A procedure to calculate the predicted adsorbed mixture activity coefficients has been presented that allows for the production of excess Gibbs energy plots at a constant reduced grand potential. This calculation does not require the knowledge of the pressure of the system as the problem can be solved by specifying the mole fractions of the adsorbed phase and the total adsorbed amount.
Results from the NIAST model of Myers indicate that the RALF model for frozen solids contains the correct structure in terms of the dependence of non-ideality with respect to composition and the reduced grand potential. In addition to providing the means to analyse the RALF model, the procedure presented can also be used to link the new framework to the traditional non-ideal adsorbed phase approach if an excess Gibbs energy generated from RALF can be matched to a suitable RAST model.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The development of the RALF model began through an Advanced Study Institute Senior Fellowship visit to the Alma Mater Studiorum Università di Bologna. Discussions and follow-up exchanges on the lattice fluid model with several colleagues in Bologna (De Angelis, Doghieri, Giacinti-Baschetti, Minelli and in particular Sarti) are gratefully acknowledged. This paper is dedicated to Giulio C. Sarti as a contribution to this Special Issue in his honour.

Conflicts of Interest

The author declares no conflict of interest.

Glossary

f fugacity (Pa)
G R residual Gibbs energy (J)
G A R residual Gibbs energy of the adsorbed phase (J)
K dimensionless Henry law constant (–)
K P Henry law constant (mol kg–1 Pa–1)
n i adsorbed amount of component i (mol kg–1)
n i s a t saturation capacity of component i (mol kg–1)
N number of moles (mol)
N j number of moles of species j (mol)
N r size of lattice in mole equivalents (mol)
m j mass of species j (kg)
m S mass of solid (kg)
M w i molecular mass of species i (kg mol–1)
P pressure (Pa)
P 0 reference pressure (Pa)
P ˜ reduced pressure (–)
P * characteristic pressure of the mixture (Pa)
P i * characteristic pressure of component i pure (Pa)
P i A * characteristic pressure of component i in the adsorbed phase (Pa)
P S * characteristic pressure of the solid (Pa)
P i j * pair characteristic pressure (Pa)
r average number of mers in a molecule (–)
r j number of mers in molecule j in the mixture (–)
r j 0 number of mers in molecule j pure (–)
R ideal gas constant (J mol–1 K–1)
T temperature (K)
T ˜ reduced temperature (–)
T * characteristic temperature of the mixture (K)
T i * characteristic temperature of component i pure (K)
T S * characteristic temperature of the solid (K)
v ˜ reduced molar volume (–)
v * average close-packed volume of mers in a mixture (m3 mer-mol–1)
v j * close-packed volume of mers molecule j pure (m3 mer-mol–1)
v j A * close-packed volume of mers molecule j in the adsorbed phase (m3 mer-mol–1)
V volume of the lattice (m3)
V * close-packed volume for the mixture (m3)
V S volume of the solid (m3)
w j mass fraction of species j (–)
x i mole fraction of species i in the adsorbed phase (–)
y i mole fraction of species i in the fluid phase (–)
z compressibility factor (–)
z E o S compressibility factor derived from the Helmholtz energy (–)

Greek Letters

β T isothermal solid compressibility (Pa–1)
Δ H 0 adsorption enthalpy at zero loading (J mol–1)
Δ U 0 adsorption energy at zero loading (J mol–1)
φ i fugacity coefficient of species i (–)
ϕ j volume fraction in the lattice occupied by species j at close-packing (–)
ϕ S volume fraction in the lattice occupied by the solid at close-packing (–)
ϕ j L volume fraction in the lattice occupied by species j (–)
κ i j pair interaction coefficient (–)
μ i R residual chemical potential of species i in the adsorbed phase (J mol–1)
μ i F R residual chemical potential of species i in the fluid phase (J mol–1)
μ i R residual chemical potential of species i at infinite dilution in the adsorbed phase (J mol–1)
μ S R m residual chemical potential of the solid on a mass basis (J kg–1)
μ S 0 R m residual chemical potential of the solid without adsorbates (J kg–1)
ρ ˜ reduced mass density (–)
ρ * average close-packed mass density in a mixture (kg m–3)
ρ j * close-packed mass density of molecule j (kg m–3)
ρ j A * close-packed mass density of molecule j in the adsorbed phase (kg m–3)
ρ S mass density of the solid (kg m–3)
ξ j A volume correction due to confinement constraints (–)
Ψ reduced grand potential (mol kg–1)

References

  1. Brandani, S. The rigid adsorbent lattice fluid model for pure and mixed gas adsorption. AIChE J. 2018, 65, 1304–1314. [Google Scholar] [CrossRef] [Green Version]
  2. Doghieri, F.; Sarti, G.C. Nonequilibrium lattice fluids:  A predictive model for the solubility in glassy polymers. Macromolecules 1996, 29, 7885–7896. [Google Scholar] [CrossRef]
  3. Sarti, G.C.; Doghieri, F. Predictions of the solubility of gases in glassy polymers based on the NELF model. Chem. Eng. Sci. 1998, 53, 3435–3447. [Google Scholar] [CrossRef]
  4. Baschetti, M.G.; Doghieri, A.F.; Sarti, G.C. Solubility in Glassy Polymers: Correlations through the Nonequilibrium Lattice Fluid Model. Ind. Eng. Chem. Res. 2001, 40, 3027–3037. [Google Scholar] [CrossRef]
  5. De Angelis, M.G.; Sarti, G.C.; Doghieri, F. Correlations between Penetrant Properties and Infinite Dilution Gas Solubility in Glassy Polymers: NELF Model Derivation. Ind. Eng. Chem. Res. 2007, 46, 7645–7656. [Google Scholar] [CrossRef]
  6. Minelli, M.; Sarti, G.C. 110th Anniversary: Gas and Vapor Sorption in Glassy Polymeric Membranes—Critical Review of Different Physical and Mathematical Models. Ind. Eng. Chem. Res. 2019, 59, 341–365. [Google Scholar] [CrossRef]
  7. Verbraeken, M.C.; Brandani, S. Predictions of Stepped Isotherms in Breathing Adsorbents by the Rigid Adsorbent Lattice Fluid. J. Phys. Chem. C 2019, 123, 14517–14529. [Google Scholar] [CrossRef]
  8. Verbraeken, M.C.; Brandani, S. A priori predictions of type I and type V isotherms by the rigid adsorbent lattice fluid. Adsorption 2019, 26, 989–1000. [Google Scholar] [CrossRef] [Green Version]
  9. Verbraeken, M.C.; Mennitto, R.; Georgieva, V.M.; Bruce, E.L.; Greenaway, A.G.; Cox, P.A.; Min, J.G.; Hong, S.B.; Wright, P.A.; Brandani, S. Understanding CO2 adsorption in a flexible zeolite through a combination of structural, kinetic and modelling techniques. Sep. Purif. Technol. 2020, 256, 117846. [Google Scholar] [CrossRef]
  10. Brandani, S.; Mangano, E.; Santori, G. Water Adsorption on AQSOA-FAM-Z02 Beads. J. Chem. Eng. Data 2022, 67, 1723–1731. [Google Scholar] [CrossRef]
  11. Myers, A.L.; Monson, P.A. Physical adsorption of gases: The case for absolute adsorption as the basis for thermodynamic analysis. Adsorption 2014, 20, 591–622. [Google Scholar] [CrossRef]
  12. Brandani, S.; Mangano, E.; Sarkisov, L. Net, excess and absolute adsorption and adsorption of helium. Adsorption 2016, 22, 261–276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Sanchez, I.C.; Lacombe, R.H. An elementary molecular theory of classical fluids. Pure fluids. J. Phys. Chem. 1976, 80, 2352–2362. [Google Scholar] [CrossRef]
  14. Lacombe, R.H.; Sanchez, I.C. Statistical thermodynamics of fluid mixtures. J. Phys. Chem. 1976, 80, 2568–2580. [Google Scholar] [CrossRef]
  15. Sanchez, I.C.; Lacombe, R.H. Statistical Thermodynamics of Polymer Solutions. Macromolecules 1978, 11, 1145–1156. [Google Scholar] [CrossRef] [Green Version]
  16. Neau, E. A consistent method for phase equilibrium calculation using the Sanchez–Lacombe lattice–fluid equation-of-state. Fluid Phase Equilibria 2002, 203, 133–140. [Google Scholar] [CrossRef]
  17. Myers, A.L.; Prausnitz, J.M. Thermodynamics of mixed-gas adsorption. AIChE J. 1965, 11, 121–127. [Google Scholar] [CrossRef]
  18. Talu, O.; Myers, A.L. Rigorous thermodynamic treatment of gas adsorption. AIChE J. 1988, 34, 1887–1893. [Google Scholar] [CrossRef]
  19. Prausnitz, J.M.; Lichtenthaler, R.N.; de Azevedo, E.G. Molecular Thermodynamics of Fluid-Phase Equilibria, 3rd ed.; Prentice Hall PTR: Upper Saddle River, NJ, USA, 1999. [Google Scholar]
  20. Brandani, S.; Mangano, E.; Luberti, M. Net, excess and absolute adsorption in mixed gas adsorption. Adsorption 2017, 23, 569–576. [Google Scholar] [CrossRef] [Green Version]
  21. Golden, T.; Sircar, S. Gas Adsorption on Silicalite. J. Colloid Interface Sci. 1994, 162, 182–188. [Google Scholar] [CrossRef]
  22. Hufton, J.R.; Danner, R.P. Chromatographic study of alkanes in silicalite: Equilibrium properties. AIChE J. 1993, 39, 954–961. [Google Scholar] [CrossRef]
  23. Abdul-Rehman, H.B.; Hasanain, M.A.; Loughlin, K.F. Quaternary, ternary, binary, and pure component sorption on zeolites. Light alkanes on Linde S-115 silicalite at moderate to high pressures. Ind. Eng. Chem. Res. 1990, 29, 1525–1535. [Google Scholar] [CrossRef]
  24. Mangano, E.; Friedrich, D.; Brandani, S. Robust algorithms for the solution of the ideal adsorbed solution theory equations. AIChE J. 2014, 61, 981–991. [Google Scholar] [CrossRef] [Green Version]
  25. Ruthven, D.M. Principles of Adsorption and Adsorption Processes; Wiley: New York, NY, USA, 1984. [Google Scholar]
  26. Talu, O.; Li, J.; Myers, A.L. Activity coefficients of adsorbed mixtures. Adsorption 1995, 1, 103–112. [Google Scholar] [CrossRef] [Green Version]
  27. Siperstein, F.R.; Myers, A.L. Mixed-gas adsorption. AIChE J. 2001, 47, 1141–1159. [Google Scholar] [CrossRef]
  28. Myers, A.L. Prediction of Adsorption of Nonideal Mixtures in Nanoporous Materials. Adsorption 2005, 11, 37–42. [Google Scholar] [CrossRef]
  29. Smith, J.M.; Van Ness, H.C.; Abbott, M.M. Introduction to Chemical Engineering Thermodynamics, 5th ed.; McGraw-Hill: New York, NY, USA, 1996. [Google Scholar]
Figure 1. Difference in gas phase mole fraction of ethane calculated from RALF and IAST as a function of adsorbed phase mole fraction. Ethane/n-butane binary mixture at 300 K and 350 kPa.
Figure 1. Difference in gas phase mole fraction of ethane calculated from RALF and IAST as a function of adsorbed phase mole fraction. Ethane/n-butane binary mixture at 300 K and 350 kPa.
Membranes 12 01009 g001
Figure 2. Reduced excess Gibbs energy as a function of the reduced grand potential at mole fractions of methane of 0.25 (dotted line), 0.5 (continuous line) and 0.75 (dashed line). Methane/n-butane binary mixture.
Figure 2. Reduced excess Gibbs energy as a function of the reduced grand potential at mole fractions of methane of 0.25 (dotted line), 0.5 (continuous line) and 0.75 (dashed line). Methane/n-butane binary mixture.
Membranes 12 01009 g002
Figure 3. Normalised reduced excess Gibbs energy as a function of the reduced grand potential at a mole fraction of methane of 0.5. RALF predictions for the binary methane/n-butane (continuous line) are compared to simple exponential decay curves with C = 0.066 (dotted line) and the predicted NIAST C 12 = 0.328 (dashed line).
Figure 3. Normalised reduced excess Gibbs energy as a function of the reduced grand potential at a mole fraction of methane of 0.5. RALF predictions for the binary methane/n-butane (continuous line) are compared to simple exponential decay curves with C = 0.066 (dotted line) and the predicted NIAST C 12 = 0.328 (dashed line).
Membranes 12 01009 g003
Figure 4. Normalised reduced excess Gibbs energy as a function of the mole fraction of methane at a fixed reduced grand potential ψ = 150 mol/kg. RALF prediction (continuous line) and two-parameter Margules equation (dotted line) are compared with NIAST predictions (dashed line) based on RALF pure component isotherms.
Figure 4. Normalised reduced excess Gibbs energy as a function of the mole fraction of methane at a fixed reduced grand potential ψ = 150 mol/kg. RALF prediction (continuous line) and two-parameter Margules equation (dotted line) are compared with NIAST predictions (dashed line) based on RALF pure component isotherms.
Membranes 12 01009 g004
Figure 5. Normalised reduced excess Gibbs energy as a function of the mole fraction of methane at a fixed reduced grand potential ψ = 10 mol/kg. RALF prediction (continuous line) and two-parameter Margules equation (dotted line) are compared with NIAST predictions (dashed line) based on RALF pure component isotherms.
Figure 5. Normalised reduced excess Gibbs energy as a function of the mole fraction of methane at a fixed reduced grand potential ψ = 10 mol/kg. RALF prediction (continuous line) and two-parameter Margules equation (dotted line) are compared with NIAST predictions (dashed line) based on RALF pure component isotherms.
Membranes 12 01009 g005
Table 1. RALF frozen solid multicomponent steps and expressions.
Table 1. RALF frozen solid multicomponent steps and expressions.
Sequence
1 Specify   parameters :   T i * ;   P i * ;   ρ i * ,   M w i ;   T S * ;   P S * ;   ρ S *   and   κ i j
2 Specify   variables :   T ;   P ;   n i ;   m S
3 v S * = R T S * P S *
4 P i A * = P i * 1 + ξ i A   ;   v i A * = R T i * P i A * ;   ρ i A * = ρ i * 1 + ξ i A ;   r i A 0 = M w i v i A * ρ i A *
5 Define   m i = n i M w i   and   m T = j m j
6 ρ * = m T i m i ρ i A * + m S ρ S *   ;   ϕ i = m i ρ i A * ρ * m T   ;   ϕ S = m S ρ S * ρ * m T   ;   ρ = m T m S ρ S
7 P * = j k ϕ j ϕ k P j k *   with   P j k * = P k j * = ( 1 κ k j ) P k * P j *   and   κ k k = 0 .
8 1 v * = j ϕ j v j *   ;   T * = P * v * R ;   T ˜ = T T * ;   P ˜ = P P *   and   ρ ˜ = ρ ρ *
9 r i = r i A 0 v i A * v *   ;   1 r = i ϕ i r i   ;   r N = m S ρ S ρ ˜ v *   and   z = r P ˜ ρ ˜ T ˜
10 x i = r ϕ i r i   and   z A E o S 1 = r [ ρ ˜ T ˜ ln ( 1 ρ ˜ ) ρ ˜ 1 ] + ( ρ ˜ ρ ˜ S ) j x j ln ϕ j x j ( 1 ϕ S )  
μ k R R T = ρ ˜ T ˜ r k ( 2 j ϕ j P k j * P * 1 ) + [ ( 1 ρ ˜ ) ln ( 1 ρ ˜ ) ρ ˜ + 1 ] r k 0 + r k r ( z A E o S 1 ) ln z + [ ρ ˜ ρ ˜ S ( 1 + ln ρ ˜ ρ ˜ S ) ] ( ln r k r ( 1 ϕ S ) + 1 r k r ( 1 ϕ S ) )   Ψ = μ S 0 R m R T μ S R m R T = 1 v S * ρ S * [ ( 1 ρ ˜ S ) ln ( 1 ρ ˜ S ) ρ ˜ S + 1 ρ ˜ S T ˜ S ] ϕ S ρ S ρ ˜ v S * [ ( 1 ρ ˜ ) ln ( 1 ρ ˜ ) ρ ˜ + 1 ] + ρ ˜ T ˜ [ 2 j ϕ j P S j * P * 1 ] ϕ S ρ S ρ ˜ v * n ( ϕ S 1 ) ( z E o s 1 ) + n
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Brandani, S. The Rigid Adsorbent Lattice Fluid Model: Thermodynamic Consistency and Relationship to the Real Adsorbed Solution Theory. Membranes 2022, 12, 1009. https://doi.org/10.3390/membranes12101009

AMA Style

Brandani S. The Rigid Adsorbent Lattice Fluid Model: Thermodynamic Consistency and Relationship to the Real Adsorbed Solution Theory. Membranes. 2022; 12(10):1009. https://doi.org/10.3390/membranes12101009

Chicago/Turabian Style

Brandani, Stefano. 2022. "The Rigid Adsorbent Lattice Fluid Model: Thermodynamic Consistency and Relationship to the Real Adsorbed Solution Theory" Membranes 12, no. 10: 1009. https://doi.org/10.3390/membranes12101009

APA Style

Brandani, S. (2022). The Rigid Adsorbent Lattice Fluid Model: Thermodynamic Consistency and Relationship to the Real Adsorbed Solution Theory. Membranes, 12(10), 1009. https://doi.org/10.3390/membranes12101009

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop