GENERALIZATION OF THERMODYNAMIC LANGMUIR ISOTHERMS FOR MIXED-GAS ADSORPTION EQUILIBRIA

Information

  • Patent Application
  • 20250035572
  • Publication Number
    20250035572
  • Date Filed
    July 22, 2022
    2 years ago
  • Date Published
    January 30, 2025
    25 days ago
Abstract
A system and method for estimating an adsorption equilibria for one or more gases from pure component adsorption isotherms includes providing one or more processors, a memory communicably coupled to the one or more processors and an output device communicably coupled to the one or more processors, calculating an adsorption of each gas on a constant monolayer adsorption surface is calculated using the one or more processors and the generalized Langmuir isotherm equations (24)-(26) or equation (27), providing the adsorption of each gas to the output device, and developing a chemical process or product is developed using the adsorption of each gas.
Description
TECHNICAL FIELD OF THE INVENTION

The present invention relates in general to the field of thermodynamic modeling, and more particularly, to a generalization of thermodynamic Langmuir isotherms for mixed-gas adsorption equilibria.


BACKGROUND OF THE INVENTION

Without limiting the scope of the invention, its background is described in connection with Langmuir isotherm modeling.


The chemical process industry is among the most essential manufacturing sectors around the globe. The United States alone has a chemical product output of greater than $979 billion in 2019. [1] Over the years, the energy demands to operate these chemical processes grow rapidly due to an increase in chemical products demand. The separation of chemical products consumes nearly 45-55% of the industrial energy and 10-15% of the overall energy produced in the United States. [2, 3] To date, distillation, membrane-based separation, and adsorptive separation are among the widely used separation technologies. [4] However, after decades of optimization and improvements, distillation processes remain the dominant separation technologies and account for more than 50% of overall separation energy consumed due to low energy efficiency. [4, 5] Membrane-based separation is an emerging and energy-efficient technology but its industrial applications continue to be limited by low selectivity, low capacity, low stability, and elevated fouling tendency. [4, 6] In contrast, adsorptive separation is relatively established, energy-efficient, and heavily practiced in the air separation industry and others. [7]


In spite of the high industrial potential of adsorptive separation, lack of experimental measurement of mixed-gas adsorption equilibrium data hampers its wide application in industry. [8, 9] Therefore, it is imperative that robust and accurate thermodynamic models be developed and validated to reliably extrapolate and predict mixed-gas adsorption equilibria from pure component adsorption isotherms. [10] Among a wide range of models and theories proposed to estimate mixed-gas adsorption equilibria are Vacancy Solution Theory (VST) [11], Loading Ratio Correlation (LRC) [12], the extended Langmuir (eL) [13], Ideal Adsorbed Solution Theory (IAST) [14], etc. Although VST predicts azeotropic adsorption behavior, it is thermodynamically inconsistent and the model parameters lack physical significance. [7,15] The LRC model can reasonably predict mixed-gas adsorption equilibria for ideal systems. [16] However, its empiricism and failure to predict adsorption azeotrope diminish LRC's applicability to model and design adsorbers. [7, 17] Similarly, the widely practiced eL model is thermodynamically inconsistent and fails to address the pressure dependence as it presents a constant selectivity at all pressures for a given mixed-gas adsorption system. [7, 18, 19] The benchmark to predict mixed-gas adsorption equilibria from pure component adsorption isotherms and the only thermodynamically consistent model is IAST. However, IAST does not account for surface heterogeneity and it fails to predict nonideal mixed-gas adsorption equilibria. [10. 20] In short, the state-of-the-art adsorption isotherm models are either thermodynamically inconsistent or incapable of describing the nonideal behavior of mixed-gas adsorption equilibria.


Recently, for the first time, Kaur et al. [21] proposed an adsorption Nonrandom Two-Liquid (aNRTL) activity coefficient model to account for the adsorbent surface heterogeneity and its underlying adsorbate-adsorbent interactions in the adsorbate phase. Subsequently, substituting adsorbate and vacant site concentrations in the classical Langmuir (cL) adsorption rate formulation with activities using aNRTL model, Chang et al. [22] presented a thermodynamic Langmuir (tL) isotherm to capture the adsorbent surface heterogeneity for pure component adsorption and isosteric heat of adsorption. [23] Moreover, Tun and Chen [10] implemented IAST in combination with the aNRTL model and showed the feasibility to reliably correlate and predict binary and multicomponent gas adsorption equilibria from pure component isotherms at constant spreading pressure. However, the adsorbate phase activity coefficients in the IAST thermodynamic framework are functions of spreading pressure. [24] As spreading pressure is not a measurable property, it is believed that no theoretical investigations have adequately addressed the spreading pressure dependence [25] for the adsorbate phase activity coefficients.


What is needed is a new Langmuir model that allows an accurate account of both surface loading dependence and adsorbate phase composition dependence for mixed-gas adsorption equilibria.


SUMMARY OF THE INVENTION

A simple and yet thermodynamically consistent multicomponent adsorption isotherm expression to estimate mixed-gas adsorption equilibria from pure component adsorption isotherms is described herein. The generalized Langmuir (gL) adsorption isotherm treats the adsorbent as an integral part of the adsorption system, assumes competitive adsorption of multiple adsorbates on adsorbent surface, tracks the available adsorbent vacant site surface area, and computes adsorbate and vacant site activities with an area-based aNRTL activity coefficient model. Devoid of the needs to consider the system spreading pressure constraint and the spreading pressure dependence of the adsorbate phase activity coefficients as imposed by IAST, the gL model allows an accurate account of both surface loading dependence and adsorbate phase composition dependence for mixed-gas adsorption equilibria.


In one embodiment, a computerized method for estimating an adsorption equilibria for one or more gases from pure component adsorption isotherms include providing one or more processors, a memory communicably coupled to the one or more processors and an output device communicably coupled to the one or more processors, calculating, using the one or more processors, an adsorption of each gas on a constant monolayer adsorption surface and the generalized Langmuir isotherm equations:







θ
i

=




n
i



A
i



A
o


=



K
i
o



y
i


P




γ
i



γ
ϕ



q
i



+




j
=
1

n





γ
i



q
j




γ
j



q
i





K
j
o



y
j


P












θ
ϕ

=




n
ϕ



A
ϕ



A
o


=

1

1
+




j
=
1

n





γ
ϕ



q
j



γ
j




K
j
o



y
j


P












q
i

=


A
i


A
ϕ






where: θi is an adsorbate phase area fraction covered with the gas i,

    • niAi is an occupied area for the gas i,
    • Kio is an intrinsic adsorption equilibrium constant of the gas i,
    • yi is a gas phase mole fraction of gas i,
    • P is a gas vapor pressure,
    • γi is an activity coefficient of the gas i,
    • γϕ is an activity coefficient of vacant sites,
    • qi is a ratio of an effective area of the gas i (Ai) and an effective area of a phantom molecule ϕ (Aϕ),
    • n is a number of the one or more gases,
    • θϕ is an adsorbate phase vacant site area fraction, and
    • nϕAϕ is a vacant area for the phantom molecule ϕ.


The adsorption of each gas is provided to the output device, and a chemical process or product is developed using the adsorption of each gas. The method can be implemented by an apparatus, system, computer, or non-transitory computer readable medium encoded with a computer program for execution by a processor that performs the steps of the method.


In one aspect, the generalized Langmuir isotherm equations reduce to








n
i


n
i
0


=



K
i
o



y
i


P


1
+




i
=
1

n



K
i
o



y
i


P








when (1) the adsorbate and vacant site effective areas are the same A1=A2= . . . =Ai=Aϕ, or equivalently, the saturation loadings of adsorbates and phantom molecule are same n10=n20= . . . =n10=nϕ0, and (2) the adsorbate phase activity coefficients are unity γiϕ=1. In another aspect, the one or more gases comprise a mixed gas having two or more components. In another aspect, the constant monolayer adsorption surface comprises activated carbon, LiLSX or Zeolite H-mordenite. In another aspect, the gas i comprises CH4, C2H4, C2H6, C3H6, N2, O2, CO2, H2S, or C3H8. In another aspect, the one or more gasses comprise a mixed gas selected from CH4—C2H4, CH4—C2H6, C2H4—C2H6, C2H4—C3H6, and C2H6—C3H6. In another aspect, the one or more gasses comprise a mixed gas selected from N2 and O2. In another aspect, the one or more gasses comprise a mixed gas selected from H2S—CO2, C3H8—H2S, and C3H8—CO2.


In another embodiment, an apparatus, system or computer includes at least one input/output interface, a data storage, and one or more processors communicably coupled to the at least one input/output interface and the data storage. The one or more processors calculate an adsorption of each gas i on a constant monolayer adsorption surface Ao using generalized Langmuir isotherm equations:







θ
i

=




n
i



A
i



A
o


=



K
i
o



y
i


P




γ
i



γ
ϕ



q
i



+




j
=
1

n





γ
i



q
j




γ
j



q
i





K
j
o



y
j


P












θ
ϕ

=




n
ϕ



A
ϕ



A
o


=

1

1
+




j
=
1

n





γ
ϕ



q
j



γ
j




K
j
o



y
j


P












q
i

=


A
i


A
ϕ






where: θi is an adsorbate phase area fraction covered with the gas i,

    • niAi is an occupied area for the gas i,
    • Kio is an intrinsic adsorption equilibrium constant of the gas i,
    • yi is a gas phase mole fraction of gas i,
    • P is a gas vapor pressure,
    • γi is an activity coefficient of the gas i,
    • γϕ is an activity coefficient of vacant sites,
    • qi is a ratio of an effective area of the gas i (Ai) and an effective area of a phantom molecule ϕ (Aϕ),
    • n is a number of the one or more gases, θϕ is an adsorbate phase vacant site area fraction, and
    • nϕAϕ is a vacant area for the phantom molecule ϕ.


      The adsorption of each gas i is provided to the output device, and a chemical process or a product is developed using the adsorption of each gas i.


In one aspect, the generalized Langmuir isotherm equations reduce to








n
i


n
i
0


=



K
i
o



y
i


P


1
+




i
=
1

n



K
i
o



y
i


P








when (1) the adsorbate and vacant site effective areas are the same A1=A2= . . . =Ai=Aϕ, or equivalently, the saturation loadings of adsorbates and phantom molecule are same n10=n20= . . . =ni0=nϕ0, and (2) the adsorbate phase activity coefficients are unity γiϕ=1. In another aspect, the one or more gases comprise a mixed gas having two or more components. In another aspect, the constant monolayer adsorption surface comprises activated carbon, LiLSX or Zeolite H-mordenite. In another aspect, the gas i comprises CH4, C2H4, C2H6, C3H6, N2, O2, CO2, H2S, or C3H8. In another aspect, the one or more gasses comprise a mixed gas selected from CH4—C2H4, CH4—C2H6, C2H4—C2H6, C2H4—C3H6, and C2H6—C2H6. In another aspect, the one or more gasses comprise a mixed gas selected from N2 and O2. In another aspect, the one or more gasses comprise a mixed gas selected from H2S—CO2, C3H8—H2S, and C3H8—CO2.


In another embodiment, a method of adsorbing one or more gases includes providing a vessel containing a constant monolayer adsorption surface, introducing one or more gases into the vessel, wherein the adsorption of each gas on the constant monolayer adsorption surface is determined by the generalized Langmuir isotherm equations:







θ
i

=




n
i



A
i



A
o


=



K
i
o



y
i


P




γ
i



γ
ϕ



q
i



+




j
=
1

n





γ
i



q
j




γ
j



q
i





K
j
o



y
j


P












θ
ϕ

=




n
ϕ



A
ϕ



A
o


=

1

1
+




j
=
1

n





γ
ϕ



q
j



γ
j




K
j
o



y
j


P












q
i

=


A
i


A
ϕ






where: θi is an adsorbate phase area fraction covered with the gas i,

    • niAi is an occupied area for the gas i,
    • Kio is an intrinsic adsorption equilibrium constant of the gas i,
    • yi is a gas phase mole fraction of gas i,
    • P is a gas vapor pressure,
    • γi is an activity coefficient of the gas i,
    • γϕ is an activity coefficient of vacant sites,
    • qi is a ratio of an effective area of the gas i (Ai) and an effective area of a phantom molecule ϕ (Aϕ),
    • n is a number of the one or more gases,
    • θϕ is an adsorbate phase vacant site area fraction, and
    • nϕAϕ is a vacant area for the phantom molecule ϕ.


The adsorption of each gas i is provided to the output device, and a chemical process or a product is developed using the adsorption of each gas i. A product can be produced in accordance with this method.


In one aspect, the generalized Langmuir isotherm equations reduce to








n
i


n
i
0


=



K
i
o



y
i


P


1
+




i
=
1

n



K
i
o



y
i


P








when (1) the adsorbate and vacant site effective areas are the same A1=A2= . . . =Ai=Aϕ, or equivalently, the saturation loadings of adsorbates and phantom molecule are same n10=n20= . . . =ni0=nϕ0, and (2) the adsorbate phase activity coefficients are unity γiϕ=1. In another aspect, the one or more gases comprise a mixed gas having two or more components. In another aspect, the constant monolayer adsorption surface comprises activated carbon, LiLSX or Zeolite H-mordenite. In another aspect, the gas i comprises CH4, C2H4, C2H6, C3H6, N2, O2, CO2, H2S, or C3H8. In another aspect, the one or more gasses comprise a mixed gas selected from CH4—C2H4, CH4—C2H6, C2H4—C2H6, C2H4—C3H6, and C2H6—C3H6. In another aspect, the one or more gasses comprise a mixed gas selected from N2 and O2. In another aspect, the one or more gasses comprise a mixed gas selected from H2S—CO2, C3H8—H2S, and C3H8—CO2.





BRIEF DESCRIPTION OF THE DRAWINGS

For a more complete understanding of the features and advantages of the present invention, reference is now made to the detailed description of the invention along with the accompanying figures and in which:



FIGS. 1A to 1C depict the data sets and model results for five pure component adsorption isotherms: (FIG. 1A) CH4, C2H4, C2H6, and C3H6 on activated carbon at 323 K [28], (FIG. 1B) N2 and O2 on LiLSX at 303.15 K [29], (FIG. 1C) CO2, H2S, and C3H8 on zeolite H-mordenite at 303.15 K [30];



FIGS. 2A to 2E depict the experimental data and model results for the five binary adsorption equilibrium compositions on activated carbon at 323 K and 0.1 bar [28]: (FIG. 2A) CH4—C2H4, (FIG. 2B) CH4—C2H6, (FIG. 2C) C2H4—C2H6, (FIG. 2D) C2H4—C3H6, and (FIG. 2E) C2H6—C3H6;



FIG. 3 depicts the experimental data and model results for the mixed-gas adsorption of N2—O2 binary mixture on LiLSX at 1.013 bar and 6.08 bar at 303.15 K [29];



FIG. 4A depicts the experimental data and the model results for mixed-gas adsorption of H2S CO2 binary mixture on zeolite H-mordenite at 303.15 K and 0.156 bar [30];



FIGS. 4B and 4C depicts the adsorption azeotropic behaviors of (FIG. 4B) C3H8—H2S and (FIG. 4C) C3H8—CO2 binary gas mixtures at 303.15 K on zeolite H-mordenite at 0.081 bar and 0.41 bar respectively [30];



FIG. 5A depicts the overall surface loading jumps of N2—O2 on LiLSX when the system pressure jumps from 1.013 bar to 6.08 bar at 303.15 K [29];



FIGS. 5B and 5C depict the adsorbent surface area tracking for (FIG. 5B) C3H8—H2S on zeolite H-mordenite at 0.081 bar and 303.15 K [30], and (FIG. 5C) C3H8—CO2 on zeolite H-mordenite at 0.41 bar and 303.15 K [30];



FIGS. 6A and 6B depict the mixed-gas adsorption equilibria phase diagram for (FIG. 6A) N2—O2 on LiLSX at 1.013 bar and 6.08 bar at 303.15 K [29], and (FIG. 6B) C3H8—H2S at 0.081 bar and C3H8—CO2 at 0.41 bar on zeolite H-mordenite at 303.15 K [30];



FIG. 7A depicts the activity coefficient of binary mixtures N2—O2 on LiLSX at 1.013 bar and 6.08 bar at 303.15 K [29];



FIGS. 7B and 7C depict the activity coefficient of binary mixtures (FIG. 7B) C3H8—H2S on zeolite H-mordenite at 0.081 bar and 303.15 K [30], and (FIG. 7C) C3H8—CO2 on zeolite H-mordenite at 0.41 bar and 303.15 K [30];



FIGS. 8A to 8C depict the ternary mixed-gas adsorption equilibria: (FIG. 8A) CH, (∘)-C2H, (Δ)-C2H, (x) on activated carbon at 323 K and 0.1 bar [28], (FIG. 8B) C2H4(∘)-C2H, (Δ)-C3H6 (x) on activated carbon at 323 K and 0.1 bar [28], (FIG. 8C) CO2 (∘)-H2S (Δ)-C3H8 (x) on zeolite H-mordenite at 303.15 K and 0.007-0.1342 bar [30], where the blue, black, and red colors represent eL, IAST, and gL predictions, respectively;



FIG. 9 is a block diagram of an apparatus or system suitable for performing the methods described herein;



FIG. 10 is a flow chart of a method for estimating an adsorption equilibria for one or more gases from pure component adsorption isotherms; and



FIG. 11 is a flow chart of a method for adsorbing one or more gases.





DETAILED DESCRIPTION OF THE INVENTION

While the making and using of various embodiments of the present invention are discussed in detail below, it should be appreciated that the present invention provides many applicable inventive concepts that can be embodied in a wide variety of specific contexts. The specific embodiments discussed herein are merely illustrative of specific ways to make and use the invention and do not delimit the scope of the invention.


To facilitate the understanding of this invention, a number of terms are defined below. Terms defined herein have meanings as commonly understood by a person of ordinary skill in the areas relevant to the present invention. Terms such as “a”, “an” and “the” are not intended to refer to only a singular entity but include the general class of which a specific example may be used for illustration. The terminology herein is used to describe specific embodiments of the invention, but their usage does not limit the invention, except as outlined in the claims.


A thermodynamically consistent model to predict mixed-gas adsorption equilibria from pure gas adsorption isotherms is described herein. A generalization of thermodynamic Langmuir isotherm for pure component adsorption, the model assumes competitive adsorption of multiple adsorbates on adsorbent surface and it applies an area-based adsorption Nonrandom Two-Liquid activity coefficient model in the activity coefficient calculations for the adsorbate phase. The resulting generalized Langmuir (gL) isotherm properly captures both surface loading dependence and adsorbate phase composition dependence for mixed-gas adsorption equilibria. The model is validated with accurate representations of gas adsorption equilibrium data for wide varieties of unary, binary, and ternary gas systems. The model results are further compared with those calculated from extended Langmuir isotherm and Ideal Adsorbed Solution Theory.


Model Formulation
Generalized Langmuir Isotherm for Pure Component Adsorption

Starting from the fundamental adsorption and desorption reactions of pure adsorbate gas A on an adsorbent surface containing vacant sites S:











A

(
g
)


+
S


AS




(
1
)







The adsorption reaction with rate constant ka results in occupied sites denoted with AS. In contrast, the desorption reaction having rate constant kd results in pure gas A and vacant sites S. Once the adsorption equilibrium has been achieved, the rates of adsorption and desorption become equal as shown in Eq. (2).












k
a



P
[
S
]


=

k
d


[
AS
]




(
2
)







Here [S] and [AS] denote the classical Langmuir (cL) site concentrations of vacant sites and occupied sites, respectively. The occupied sites can be expressed in the amount adsorbed for adsorbate component 1, n1. The vacant sites can be represented as (n10−n1), where n10 is the saturation amount adsorbed for component 1. The apparent adsorption equilibrium constant K1 can be expressed as a function of the adsorbate phase site concentrations as:










K
1

=



k
o


k
d


=



[
AS
]


P
[
S
]


=



n
1


P

(


n
1
0

-

n
1


)


=


x
1


P

(

1
-

x
1


)









(
3
)







where x1 is the ratio of n1 and n10. Simplification of Eq. (3) for x1 yields the classical Langmuir isotherm shown in Eq. (4):











n
1


n
1
0


=



K
1


P


1
+


K
1


P








(
4
)








The classical Langmuir isotherm disregards the adsorbent surface heterogeneity upon assuming the rates of adsorption and desorption are proportional to the site concentrations. Chang et al. [22] proposed the thermodynamic Langmuir isotherm to address the surface heterogeneity by substituting the site concentrations with the site activities expressed as the product of site concentration and site activity coefficient as shown in Eq. (5), leading to the thermodynamic Langmuir isotherm expressed as Eq. (6):










K
1
o

=



x
1



γ
1




γ
ϕ



x
ϕ


P







(
5
)















n
1


n
1
0


=



K
1
o


P




γ
1


γ
ϕ


+


K
1
o


P







(
6
)







where K1o is the intrinsic adsorption equilibrium constant of adsorbate component 1, x1 is the ratio of n1 and n10, and xϕ, representing the vacant site fraction on the adsorbent surface, is calculated as (1−x1). y1 is the activity coefficient of component 1 on the occupied sites and γϕ is the activity coefficient of “phantom” molecule ϕ on the vacant sites. The reference state for adsorbate component 1 is the adsorbate phase fully occupied with component 1 while the reference state for the phantom molecule ϕ is the absorbate phase with vacant sites only.


In addition to addressing the surface heterogeneity, this work further proposes that there is a constant total adsorbent surface area, Ao, which is covered with adsorbate component 1 with the “effective” molecular area, A1. The adsorbate phase area fraction covered with component 1, θ1, is a function of Ao, A1, and n1, expressed as:










n
T

=


n
1

+

n
ϕ






(
7
)














x

1

=


n
1


n
T






(
8
)













x
ϕ

=


n
ϕ


n
T






(
9
)













θ
1

=



n
1



A
1



A
o







(
10
)














θ
ϕ

=



n
ϕ



A
ϕ



A
o







(


11

)













A
o

=




n
1



A
1


+


n
ϕ



A
ϕ



=



n
1
0



A
1


=


n
ϕ
0



A
ϕ









(
12
)








where nϕ and nϕ0 are the remaining amount and the maximum amount of vacant sites, respectively, Aϕ is the “effective” area of phantom molecule ϕ, with nitrogen chosen as the model molecule for ϕ. x1 is the ratio of n1 and nT referred as the adsorbate phase mole fraction of component 1, and xϕ is the ratio of nϕ and nT and referred as the adsorbate phase mole fraction of vacant sites. While Ao remains constant, the occupied site area fraction, θ1, and the vacant site area fraction, θϕ, change with respect to the loading of component 1. Upon considering both the site activities and the site area fractions, the generalized Langmuir isotherm is expressed in Eqs. (13)-(14) for adsorbate component 1 and vacant site area fractions:










θ
1

=




n
1



A
1



A
o


=



K
1
o


P




γ
1



γ
ϕ



q
1



+


K
1
o


P









(
13
)














θ
ϕ

=




n
ϕ



A
ϕ



A
o


=

1

1
+




γ
ϕ



q
1



γ
1




K
1
o


P









(
14
)














q
1

=


A
1


A
ϕ







(
15
)








Here, q1 is the ratio of effective area of component 1 and the effective area of phantom molecule ϕ and it should depend on the adsorbent surface characteristics, adsorption temperature, and the local minimal energy adsorbate molecular configuration on the adsorbent surface. [26] γ1 and γϕ are the activity coefficients of component 1 and phantom molecule, respectively, which are functions of x1, xϕ, A1, and Aϕ. Eq. (13) reduces to the thermodynamic Langmuir (tL) isotherm, Eq. (6), when the effective area of adsorbate is same as the effective area of phantom molecule. Eq. (13) further reduces to the classical Langmuir isotherm, Eq. (4), when both the effective area of adsorbate and the effective area of phantom molecule are same and the activity coefficients are unity.


Generalized Langmuir Isotherm for Mixed-Gas Adsorption

The generalized Langmuir isotherm for pure component adsorption can be readily extended for mixed-gas adsorption. For a competitive adsorption of multiple adsorbates on a constant monolayer adsorption surface area, A°, the corresponding generalized Langmuir isotherm equations for mixed-gas adsorption equilibria can be expressed as:










K
i
o

=



x
i



γ
i




γ
ϕ



x
ϕ



Py
i








(
16
)














n
T

=





i
=
1

n


n
i


+

n
ϕ







(
17
)














x
i

=


n
i


n
T







(
18
)















x

ϕ

=


n
ϕ


n
T







(
19
)


















i
=
1

n


x
i


+

x
ϕ


=
1





(
20
)














θ
i

=



n
i



A
i



A
o






(
21
)













θ
ϕ

=



n
ϕ



A
ϕ



A
o







(
22
)














A
o

=






i
=
1

n



n
i



A
i



+


n
ϕ



A
ϕ



=



n
i
0



A
i


=


n
ϕ
0



A
ϕ









(
23
)








where γi is the gas phase mole fraction of adsorbate component i, n is the number of adsorbates, nT is the sum of the total amount adsorbed of all adsorbates and the remaining amount of vacant sites. θi is the adsorbate phase area fraction covered with component i, and θϕ is the adsorbate phase vacant site area fraction. Eq. (16) defines the intrinsic adsorption equilibrium constant of component i while Eqs. (17)-(20) express the adsorbate phase mole fractions of component i and vacant sites and their relationships. Eqs. (21) and (22) track the adsorbate phase area fractions occupied by adsorbate component i and vacant sites, respectively. Eq. (23) shows, while the total adsorbent surface area, Ao, remains constant, the occupied area for component i, niAi and the are for vacant sites, nϕAϕ, vary with loading and adsorbate phase compositions.


Rearrangement of Eqs. (16)-(20) in combination with Eqs. (21)-(22) results in the generalized Langmuir isotherm equations for mixed-gas adsorption equilibria as presented in Eqs. (24) and (25).










θ
i

=




n
1



A
1



A
o


=



K
i
o



y
i


P




γ
i



γ
ϕ



q
i



+





j
=
1


n





γ
i



q
j




γ
j



q
i





K
j
o



y
j


P










(
24
)














θ
ϕ

=




n
ϕ



A
ϕ



A
o


=

1

1
+





j
=
1


n





γ
ϕ



q
j



γ
j




K
j
o



y
j


P










(
25
)














q
i

=


A
i


A
ϕ







(
26
)








Here, qi is the ratio of effective area of component i and effective area of phantom molecule, γi and γϕ represent the activity coefficients of component i and vacant sites, respectively, which are functions of xi's, xϕ, Ai's, and Aϕ. In summary, solving n number of Eq. (24) for n number of adsorbates coupled with Eq. (25) for nϕ along with an appropriate activity coefficient model results in the generalized Langmuir isotherm for mixed-gas adsorption equilibria. The generalized Langmuir isotherm reduces to the extended Langmuir isotherm, shown in Eq. (27), if the following two conditions are satisfied: 1) the adsorbate and vacant site effective areas are the same, i.e., A1=A2= . . . =Ai=Aϕ, or equivalently, the saturation loadings of adsorbates and phantom molecule are same, i.e., n10=n20= . . . =ni=n00, and 2) the adsorbate phase activity coefficients are unity. Devoid of the simplifying assumptions of the extended Langmuir isotherm, the generalized Langmuir isotherm provides a thermodynamically consistent generalization of the extended Langmuir expression.











n
i


n
i
0


=



K
i
o



y
i


P


1
+




i
=
1

n



K
i
o



y
i


P








(
27
)







Area-Based Multicomponent aNRTL Model


The original adsorption NRTL model derived from the two-fluid theory does not take into account the effective areas of adsorbates and phantom molecule on the adsorbate phase. [21] To make aNRTL model consistent with the generalized Langmuir model formulation, the aNRTL model formulation has been revisited to incorporate the adsorbate area fractions. The derivation of area-based aNRTL model is provided below.


Consider an adsorbent surface containing energetically non-interacting multiple adsorption sites where each site “0” is surrounded with the molecules of n number of adsorbates. To address the local composition of multicomponent adsorbate phase, Nonrandom two-liquid theory [36] has been implied while incorporating the size of adsorbates. So, the sum of local area fractions of n adsorbates around the adsorbent site “0” shall be unity, written as:













i
=
1

n


θ

i

0



=
1




(
28
)







Here, θi0 is the local area fractions of ith adsorbate near the adsorbent site “0” that can be expressed in the form of area fraction and mole fraction as follow:










θ

i

0


=




θ
i



exp



(

-


α


g

i

0



RT


)






k
=
1

n



θ

k





exp



(

-


α


g

k

0



RT


)




=




θ
i



exp



(


-
α





g

i

0


-

g
00


RT


)






k
=
1

n



θ

k





exp



(


-
α





g

k

0


-

g
00


RT


)




=



θ
i



G

i

0







k
=
1

n



θ

k






G

k

0












(
29
)














θ
i

=



x
i



q
i






i
=
1

n



x
i



q
i









(
30
)








Here, gi indicates the interaction energy between the ith adsorbate and the adsorption site. α denotes the non-randomness factor and qi is the effective area of the ith adsorbate.


To derive the expression for molar excess Gibbs energy for area based aNRTL model, the multicomponent adsorbate system is treated with the two-liquid theory. [37] Starting from the sum of changes of residual Gibbs energy of the adsorbate phase:










g
mix

=




j
=
1

n


=


x
j



q
j



g

(
r
)










(
31
)














g
E

=




j
=
1

n



x
j




q
j

(





i
=
1

n



θ

i

0




g

i

0




-

g

j

0



)








(
32
)








Here g is the residual Gibbs energy. Now, using Eq. (29) and (32) results in:










g
E

=






j
=
1


n



x
j



q
j









i
=
1


n



g

i

0






θ
i



exp



(

-


α


g

i

0



RT


)







k
=
1


n



θ

k





exp



(

-


α


g

k

0



RT


)






-





j
=
1


n



x
j



q
j



g

j

0










(
33
)














g
E

=






i
=
1


n




x
i



q
i



exp




(

-


α


g

i

0



RT


)

·

[






j
=
1


n



x
j



q
j



g

i

0




-





j
=
1


n



x
j



q
j



g

j

0





]








k
-
1

n



x
k



q
k



exp



(

-


α


g

k

0



RT


)









(
34
)








Eq. (34) can be simplified to Eq. (35):











g
E

RT

=




i
=
1

n



x
i



q

i





G

i

0








j
=
1

n



x
j




q
j

[



g

i

0


-

g

j

0



RT

]







k
=
1

n



x
k



q
k



G

k

0












(
35
)








Simplification of Eq. (35) results in the final expression of Gibbs energy for multicomponent adsorbate systems, written as:











g
E

RT

=




i
=
1

n



x
i



q

i









j
=
1

n



x
j



q

j





τ

ij









k
=
1

n



x
k



q
k



G

k

i












(
36
)








Afterwards, the activity coefficients are related to the partial molar excess Gibbs energy expressed as follow:










ln



γ
i


=







n
i





[



n
T



g
E


RT

]


T
,
π
,

n

j

i










(
37
)








Rearranging Eq. (36) and then substituting in Eq. (37) results in:












n
T



g
E


RT

=



p



n
p



q
p






j



n
j



q
j



τ

pj








k



n
k



q
k



G
kp











(
38
)


















n
T



g
E


RT

=





n
i



q
i






j



n
j



q
j



τ
ij






k



n
k



q
k



G
ki









1
st



term


+






p

i




n
p



q
p






j



n
j



q
j



τ
pj






k



n
k



q
k



G
kp










2
nd



term







(
39
)







Here, p=j=k=1, 2, 3, . . . , i, . . . , m. Now, taking the derivative of Eq. (39) with respect to ng at constant T, π, and nj≠i;















n
i





[



n
T



g
E


RT

]



1

st





term



=



-

n
i





q
i
2

·




j



n
j



q
j



τ
ij





(



k



n
k



q
k



G
ki



)

2




+


q
i

·




j



n
j



q
j



τ
ij






k



n
k



q
k



G
ki










(
40
)







Upon simplification, Eq. (40) can be written as:















n
i





[



n
T



g
E


RT

]



1

st





term



=



-

n
i





q
i
2

·




j



n
j



q
j



τ
ij





(



k



n
k



q
k



G
ki



)

2




+


q
i

·




j



n
j



q
j



τ
ij



G
ij






k



n
k



q
k



G
kj










(
41
)







Now, taking the partial derivative of 2nd term as expressed in Eq. (41):















n
i





[



n
T



g
E


RT

]



2

nd





term



=




p

i




n
p



q
p




(





k



n
k



q
k




G
kp

·

q
i




τ
pi



-



j



n
j



q
j




τ
pj

·

q
i




G
ip






(



k



n
k



q
k



G
kp



)

2


)







(
42
)







Upon simplification of Eq. (42), the resulting expressions are:















n
i





[



n
T



g
E


RT

]



2
nd



term



=


q
i







p

i



[




j

i




n
p



q
p




G

p

0


·

n
j




q
j




(



τ
pi



G

j

0



-


τ
pj



G

i

0




)



]




(



k



n
k



q
k



G

k

0




)

2







(
43
)


















n
i





[



n
T



g
E


RT

]



2
nd



term



=



q
i





j



n
j



q
j






p



n
p



q
p




G
pj

(


τ
pi

-


τ
pj



G
ij



)





(



k



n
k



q
k



G

k
j




)

2





+


q
i





j



n
j



q
j





n
i



q
i



τ
ij



G
ij
2




(



k



n
k



q
k



G
kj



)

2










(
44
)







Now, combining Eq. (41) and (44) results in:










ln



γ
i


=


q
i





j



n
j




q
j


[




τ
ij



G
ij





k



n
k



q
k



G
kj




+




p



n
p



q
p




G
pj

(


τ
pi

-


τ
pj



G
ij



)





(



k



n
k



q
k



G
kj



)

2



]








(
45
)







Simplification of Eq. (45) provides the final expression of activity coefficient of component i in the multicomponent gas adsorption equilibria system.










ln



γ
i


=







n
i





[



n
T



g
E


RT

]


T
,
π
,

n

j

i





=


q
i







j
=
1

n



x
j
2



q
j
2




τ
ij


[


G
ij

-
1

]





(




k
=
1

n



x
k



q
k



G
kj



)

2








(
46
)







Here, qi=Ai/Aϕ is the ratio of effective adsorbate area of component i and that of the model molecule (nitrogen) for phantom molecule, and xi is the adsorbate phase mole fraction of component i.










τ
ij

=


-

τ
ji


=



g

i

0


-

g

j

0



RT






(
47
)













G
ij

=


1

G
ji


=

exp

(


-
α



τ
ij


)






(
48
)







It is worth noting that the gL isotherm treats the adsorbent as a part of the adsorption system. In other words, the adsorbate phase of a single component gas adsorption system is treated as a binary system of adsorbate component 1 and phantom molecule ϕ for adsorbent vacant sites with the composition of (x1, xϕ). Likewise, the adsorbate phase of a binary gas adsorption system is treated as a ternary system of two adsorbates and phantom molecule with the composition of (x1, x2, xϕ). For a mixed-gas adsorption system with n components, adsorbent vacant sites or phantom molecule ϕ included in the mixture, the excess Gibbs energy expression for the adsorbate phase is:











g
E

RT

=




i
=
1

n



x
i



q
i







j
=
1

n



x
j



q
j



τ
ij







k
=
1

n



x
k



q
k



G
ki










(
49
)







The corresponding activity coefficient expression is given in Eq. (50):










ln



γ
i


=


q
i








j
=
1


n



x
j
2



q
j
2




τ
ij

[


G
ij

-
1

]





(




k
=
1

n



x
k



q
k



G
kj



)

2







(
50
)













τ
ij

=


-

τ
ji


=



g

i

0


-

g

j

0



RT






(
51
)













G
ij

=

exp



(

-

ατ
ij


)






(
52
)







where α is the nonrandomness factor fixed at 0.3 and τij's are the adjustable binary interaction parameters for the i-j pair, with τij=−τji. In short, the area-based adsorption NRTL activity coefficients are functions of xi's, xϕ, Ai's, and Aϕ, and τij's.


Results and Discussion
Pure Component Adsorption Isotherm

Accurate representation of experimental pure component adsorption isotherm data is the first step towards successful modeling of mixed-gas adsorption equilibria. The classical Langmuir (cL), thermodynamic Langmuir (tL), and generalized Langmuir (gL) models are used herein to represent pure component adsorption isotherms. Requiring no activity coefficient calculations, cL makes use of two model parameters (Ki and ni0). In contrast, both tL and gL isotherms require three model parameters (τ, Kio, and ni0). For internal consistency, the original aNRTL model [21] is used in the activity coefficient calculation for tL and the area-based aNRTL model presented above is used with gL. Note that, with literature reported values for Ao and Ai the ni0 parameter in gL can be calculated per Eq. (23) and leave gL with only two adjustable model parameters, i.e., τand Kio. The effective areas in gL, Ai's, are notably dependent on adsorbent, adsorption temperature, and the configuration of adsorbate molecule on the adsorbent surface. [26] The recommended effective areas of adsorbates and adsorbent surface areas are available in the literature. [26] To regress the adsorption isotherm parameters, a Maximum Likelihood Principle-based objective function is minimized, [27] given as follow:









Obj
=



i



(



n
i
calc

-

n
i
expt



σ
expt


)

2






(
53
)







where σ indicates the standard deviation of the experimental data set to 0.05 mol/kg while nicalc and niexpt denote the calculated and experimentally measured adsorption amounts. To evaluate the performance of isotherm fitting, Root Mean Square Error (RMSE) is used, and it is expressed as:









RMSE
=





i



(


n
i
calc

-

n
i
expt


)

2


N






(
54
)







where N represents the total number of data points.


Nine pure component adsorption isotherm data sets are selected to support subsequent mixed-gas adsorption equilibria studies and they have been used to identify the model parameters of cL, tL, and gL pure component isotherms. These data include adsorption isotherms of CH4, C2H4, C2H6, and C3H6 on activated carbon at 323 K, [28] N2 and O2 on LiLSX at 303.15 K [29], and CO2, H2S, and C3H8 on zeolite H-mordenite at 303.15 K [30]. The graphical representation of these data sets and the model results are shown in the FIGS. 1A to 1C. The corresponding regression results are summarized in Table 1, Table 2, and Table 3 for cL, tL, and gL, respectively. Note that the identification of tL isotherm parameters for systems involving only data in the Henry's law region, i.e., at low relative pressures, is more challenging due to lack of information on saturation capacities. For example, the much higher uncertainty in estimating ni0 parameter of tL for O2 on LiLSX, as shown in Table 2, can be attributed to the fact that the experimental data are limited to the Henry's law region.









TABLE 1







Regressed parameters for classical Langmuir isotherm















T
K
ni0
RMSE



Adsorbent
Adsorbate
(K)
(bar−1)
(mol/kg)
(mol/kg)
Ref.
















Activated
CH4
323
0.402 ± 0.265
1.44 ± 0.70
0.003
[28]


Carbon
C2H4
323
2.265 ± 0.196
2.14 ± 0.09
0.041



C2H6
323
2.511 ± 0.271
2.13 ± 0.10
0.046



C3H6
323
9.020 ± 0.366
2.96 ± 0.04
0.143


LiLSX
N2
303.15
0.435 ± 0.029
2.54 ± 0.07
0.013
[29]



O2
303.15
0.046 ± 0.029
3.70 ± 0.07
0.011


Zeolite H-
CO2
303.15
7.459 ± 0.054
2.41 ± 0.01
0.129
[30]


mordenite
H2S
303.15
102.193 ± 3.584 
2.14 ± 0.02
0.120



C3H8
303.15
33.413 ± 0.754 
1.12 ± 0.01
0.107
















TABLE 2







Regressed parameters for thermodynamic Langmuir isotherm
















T
Kio
nio

RMSE



Adsorbent
Adsorbate
(K)
(bar−1)
(mol/kg)
τ
(mol/kg)
Ref.

















Activated
CH4
323
0.402 ±
1.44 ±
−0.02 ±
0.003
[28]


Carbon


0.292
0.78
8.06



C2H4
323
0.148 ±
6.73 ±
−2.18 ±
0.014





0.398
6.63
0.79



C2H6
323
0.193 ±
6.13 ±
−2.21 ±
0.021





0.480
5.55
0.71



C3H6
323
0.942 ±
6.67 ±
−2.52 ±
0.028





0.732
1.95
0.39


LiLSX
N2
303.15
0.344 ±
2.83 ±
−0.91 ±
0.006
[29]





0.047
0.17
0.16



O2
303.15
0.0033 ±
15.01 ±
−1.77 ±
0.011





0.045
85.14
4.56


Zeolite H-
CO2
303.15
1.632 ±
4.14 ±
−2.42 ±
0.016
[30]


Mordenite


0.544
0.44
0.11



H2S
303.15
75.082 ±
2.78 ±
−3.09 ±
0.022





10.181
0.11
0.19



C3H8
303.15
41.038 ±
1.27 ±
−3.02 ±
0.050





4.537
0.05
0.11
















TABLE 3







Regressed parameters for generalized Langmuir isotherm

















T
Kio
nio

RMSE
Ai



Adsorbent
Adsorbate
(K)
(bar−1)
(mol/kg)
τ
(mol/kg)
(nm2)
Ref.


















Activated
CH4
323
0.017 ±
7.27 **
−1.99 ±
0.003
0.160 33
[28]


Carbon


0.015

0.56



C2H4
323
0.192 ±
5.17 **
−2.14 ±
0.014
0.225 33





0.094

0.17



C2H6
323
0.196 ±
5.17 **
−1.85 ±
0.021
0.225 33





0.104

0.17



C3H6
323
2.121 ±
3.88 **
−1.59 ±
0.046
0.300 32





1.793

0.34


LiLSX
N2
303.15
0.344 ±
2.83 ±
−0.91 ±
0.006
0.162 26
[29]





0.047
0.17
0.16



O2
303.15
0.062 ±
3.25 **
−0.02 ±
0.012
0.141 26





0.024

1.59


Zeolite H-
CO2
303.15
1.619 ±
4.14 ±
−2.42 ±
0.016
0.163 26
[30]


Mordenite


0.544
0.44
0.11



H2S
303.15
26.553 ±
3.21 **
−3.10 ±
0.028
0.210 35





27.009

0.51



C2H
303.15
0.877 ±
2.25 **
−3.74 +
0.113
0.300 32





2.477

1.97





** nio is calculated from Ao






The results show both tL and gL perform much better than cL. In addition, while tL and gL show very similar RMSE's in fitting the data, some of the values of tL isotherm parameter ni0 seem unreasonable. For example, the regressed tL ni0 value for CH4 adsorbed on activated carbon is exceedingly small and only a fraction (˜20%) of that of C2H4. In contrast, given a constant Ao for an adsorbent, the values of gL isotherm parameter ni0 are inversely proportional to the values of Ai. Note that the adsorbate saturation loadings, i.e., ni0's, on activated carbon shown in Table 3 are calculated from the reported adsorbent surface area of






700





m
2

g

[

28
,
31

]





and the corresponding adsorbate effective areas in the literature [32, 33]. For LiLSX, Ai is known for N2 [26], ni0 for N2 is fitted to the isotherm data, and Ao is then calculated to be






276





m
2

g

.





For zeolite H-mordenite, Ai is known for CO2 [26], ni0 for CO2 is fitted to the isotherm data, and Ao is then calculated to be







406




m
2

g


,




consistent with the literature value of









429




m
2

g

.





[
34
]







Binary Mixed-Gas Adsorption Equilibria

Ten binary mixed-gas adsorption systems have been investigated with eL, IAST, and gL. The regressed gL parameters and the corresponding RMSE's for all the binary mixed-gas adsorption systems with the three models are reported in Table 4. For internal model consistency, the cL model is used to calculate pure component adsorption isotherms in the eL calculations for mixed-gas adsorption equilibria. The tL model is used to calculate pure component adsorption isotherms and the corresponding spreading pressures in the IAST calculations. In contrast, the gL isotherm is used in both pure component adsorption isotherms and mixed-gas adsorption equilibria.









TABLE 4







Model results for binary mixed-gas adsorption equilibria














T
P
IAST
eL
gL
















Adsorbent
Adsorbates
(K)
(bar)
RMSE
RMSE
τ12
RMSE
Ref.


















Activated
CH4(1)—C2H4(2)
323
0.1
0.059
0.114
0
0.062
[28]


Carbon
CH4(1)—C2H6(2)
323
0.1
0.161
0.233
0
0.163



C2H4(1)—C2H6(2)
323
0.1
0.041
0.061
−0.54 ±
0.044








0.00



C2H4(1)—C3H6(2)
323
0.1
0.080
0.196
−1.34 ±
0.084








0.00



C2H6(1)—C3H6(2)
323
0.1
0.040
0.128
−2.36 ±
0.036








0.01


LiLSX
N2(1)—O2(2)
303.15
1.013
0.013
0.004
−1.83 ±
0.005
[29]








0.00



N2(1)—O2(2)
303.15
6.080
0.015
0.028
−1.83 ±
0.001








0.00


Zeolite H-
H2S (1)—CO2(2)
303.15
0.156
0.049
0.044
−2.95 ±
0.013
[30]


mordenite





0.00



C3H8 (1)—H2S(2)
303.15
0.081
0.145
0.181
−3.93 ±
0.038








0.00



C3H8 (1)—CO2(2)
303.15
0.410
0.162
0.212
−4.87 ±
0.085








0.00









Five binary mixed-gas adsorption of hydrocarbons on activated carbon at 323 K and 0.1 bar [28] are investigated. gL estimations accurately match the experimental data and outperform eL predictions for the five binary systems of CH4—C2H4, CH4—C2H6, C2H4—C2H6, C2H4—C3H6, and C2H6—C3H6. RMSEs of IAST and gL are same for all five binary systems on activated carbon.



FIGS. 2A to 2E show the experimental data [28] and the model results for the five binary systems on activated carbon at 323 K and 0.1 bar [28]: (FIG. 2A) CH4—C2H4, (FIG. 2B) CH4—C2H6, (FIG. 2C) C2H4—C2H6, (FIG. 2D) C2H4—C3H6, and (FIG. 2E) C2H6—C3H6. The eL predictions deviate much from the experimental adsorbate phase composition while the IAST and gL results match or are close to the data. Note that the adsorbate phase mole fractions, x, calculated from gL with Eqs. (18)-(19) are referred as “true” adsorbate phase mole fractions and they can be transformed to “apparent” adsorbate phase mole fractions, xi, upon renormalization with xϕ excluded. “Apparent” adsorbate phase mole fractions are used to generate all the figures herein.



FIG. 3 presents the experimental data [29] and the model results for the mixed-gas adsorption of N2—O2 binary mixture on LiLSX at pressure of 1.013 bar and 6.08 bar at 303.15 K.29 While eL satisfactorily predicts the mixed-gas adsorption equilibria at 1.013 bar, eL does not address the system pressure effect and its predictions for 6.08 bar remain unchanged from those for 1.013 bar and significantly depart from the experimental data. Although the IAST predictions show correct pressure dependence, the IAST predictions for both 1.013 bar and 6.08 bar deviate from the experimental data. Interestingly, gL successfully correlates the adsorption equilibria data at 6.08 bar with τ12=−1.83, and then accurately predicts the adsorption equilibria at 1.013 bar.



FIG. 4A shows the experimental data [30] and the model results for mixed-gas adsorption of H2S—CO2 binary mixture on zeolite H-mordenite at 303.15 K and 0.156 bar. Neither IAST nor eL predict the adsorption behavior while gL accurately correlates the mixed-gas adsorption equilibria with τ12=−2.95. FIG. 4B and FIG. 4C show the adsorption azeotropic behaviors of C3H8—H2S and C3H8—CO2 binary gas mixtures at 303.15 K on zeolite H-mordenite at 0.081 bar and 0.41 bar, respectively. For these azeotropes-forming systems, both IAST and eL fail to predict the azeotropic behaviors while gL reliably captures the adsorption azeotropes with τ12=−3.93 for the C3H8—H2S binary and τ12=−4.87 for the C3H8—CO2 binary.


Estimate Surface Loading from Adsorption Equilibria


The unique feature of gL to consider adsorbent as a part of adsorption system makes the tracking of adsorbent surface loading possible. For the adsorption of N2—O2 binary on LiLSX at 1.013 bar and 6.08 bar and at 303.15 K [29], FIG. 5A shows the overall surface loading jumps when the system pressure changes from 1.013 bar (black lines denoted as 502) to 6.08 bar (red lines denoted as 504). It further shows that at low pressure, i.e., 1.013 bar, N2 adsorption and O2 adsorption are relatively independent of each other. As the gas phase N2 mole fraction increases, the N2 adsorbate phase area fraction increases linearly, albeit with a slight positive deviation, while the O2 adsorbate phase area fraction declines linearly. At high pressure, i.e., 6.08 bar, the adsorption competition between O2 and N2 strengthens. The slight positive deviation from the linear relationship for the N2 adsorbate phase area fraction vs the gas phase N2 mole fraction is magnified while the O2 adsorbate phase area fraction shows negative deviation from the linear relationship.


For the two binary mixtures exhibiting azeotropic behavior, i.e., C3H8—H2S binary on zeolite H-mordenite at 0.081 bar and 303.15 K [30] and C3H8—CO2 binary on zeolite H-mordenite at 0.41 bar and 303.15 K [30], FIG. 5B and FIG. 5C show both H2S and CO2 adsorptions drop as their gas phase mole fractions drop while, responsible for the azeotropic behavior, the C3H8 adsorbate phase area fractions show logit S-shape behavior as the gas phase C3H8 mole fractions increase. Interestingly, both these two azeotrope-forming systems exhibit a maximum total absorbent surface coverage at the azeotropes.


The mixed-gas adsorption equilibria phase diagram analogous to P-xy phase diagram in vapor-liquid equilibria can be generated from gL. FIG. 6A shows the θ-x′y phase diagrams of N2—O2 binary at 1.013 bar (lower pair of black lines denoted as 602) and 6.08 bar (upper pair of red lines denoted as 604) on LiLSX at 303.15 K. Here θ is the overall surface loading, i.e., (θ12). The θ-x′y phase diagram of N2—O2 binary at 1.013 bar shows nearly “Raoult's law” type phase diagram and the nonideality is much enhanced for the system at 6.08 bar. FIG. 6B shows the θ-x′y phase diagrams of C3H8—H2S binary at 0.081 bar (upper pair of black lines denoted as 606) and C3H8—CO2 binary at 0.41 bar (lower pair of red lines denoted as 608) on zeolite H-mordenite at 303.15 K. At azeotropes, the systems exhibit maximum overall surface loading and the apparent adsorbate phase component mole fractions equal to the gas phase component mole fractions.


Estimate Activity Coefficients from Adsorption Equilibria


Besides tracking the adsorbent surface loading, gL captures the adsorbent surface heterogeneity by tracking the adsorbate phase activity coefficients of adsorbates and vacant sites. As previously mentioned, full adsorbent surface coverage of pure adsorbate molecules is defined as the reference state for adsorbates. For the adsorption of N2—O2 binary on LiLSX at 1.013 (black lines denoted as 702) and 6.08 bar (red lines denoted as 704), and 303.15 K [29], FIG. 7A shows the activity coefficients of adsorbates and vacant sites as function of apparent adsorbate phase N2 mole fraction. With τ≈0 for O2, activity coefficients for O2 and vacant sites remain close to unity until the adsorbent surface is significantly covered with N2. On the other hand, activity coefficient for N2 increases and approaches unity as the N2 adsorbate phase area fraction increases with increasing apparent adsorbate phase mole fraction and system pressure.



FIG. 7B and FIG. 7C show the activity coefficients of adsorbates and vacant sites as function of apparent adsorbate phase C3H8 mole fraction for the C3H8—H2S binary on zeolite H-mordenite at 0.081 bar and 303.15 K [30] and the C3H8—CO2 binary at 0.41 bar and 303.15 K [30]. With t=−3.74 for C3H8, the activity coefficient of C3H8 at infinite dilution is very small for both binary systems, indicative of strong attractive interaction between the adsorbate and the adsorbent. The activity coefficient of C3H8 then increases sharply and reaches a plateau as its apparent adsorbate phase mole fraction increases. As expected, the corresponding activity coefficients of CO2 and H2S in their respective binary adsorption systems with C3H8 are close to unity (˜0.8) when their apparent adsorbate phase mole fractions are unity. The activity coefficients then drop to around 0.1 as their respective apparent adsorbate phase mole fraction approaches zero. The activity coefficients of vacant sites remain relatively unchanged since the overall surface loadings, θ's, stay relatively constant. The vacant sites activity coefficients hover around 0.4 to 0.5 for the C3H8—H2S binary and drop from 0.7 to 0.3 for the C3H8—CO2 binary as the C3H8 apparent adsorbate phase mole fractions vary from zero to unity.


Prediction of Mixed-Gas Adsorption Equilibria

To examine the predictive capability of gL, the mixed-gas adsorption equilibria is predicted for CH4—C2H4—C2H6 ternary and C2H4—C2H6—C3H6 ternary on activated carbon at 323 K and 0.1 bar [28] and CO2— H2S—C3H8 ternary on zeolite H-mordenite at 303.15 K and 0.007-0.1342 bar [30]. The predictions with eL, IAST, and gL are based on the pure component adsorption isotherm parameters reported in Table 1, Table 2, and Table 3 and the regressed τij's reported in Table 4. The parity plots of above mentioned mixed-gas adsorption equilibria predictions for apparent adsorbate phase compositions are presented in FIGS. 8A-8C. FIG. 8A depicts CH4 (∘)-C2H4(Δ)-C2H, (x) on activated carbon at 323 K and 0.1 bar [28]. FIG. 8B depicts C2H4(∘)-C2H, (Δ)-C3H6(x) on activated carbon at 323 K and 0.1 bar [28]. FIG. 8C depicts CO2 (∘)-H2S (Δ)-C3H8 (x) on zeolite H-mordenite at 303.15 K and 0.007-0.1342 bar [30]. The blue, black, and red colors represent eL, IAST, and gL predictions, respectively.


The average relative deviation (ARD %) results are reported in Table 5.










ARD


%

=



1

0

0

N






i
=
1

N




"\[LeftBracketingBar]"




x

i
,
calc



-

x

i
,
expt





x

i
,
expt






"\[RightBracketingBar]"








(
55
)














TABLE 5







Average relative deviation (ARD %) of ternary mixed-gas


adsorption equilibria systems











Adsorbents
Adsorbates
IAST
eL
gL





Activated Carbon at 323K and
CH4(1)—C2H4(2)—C2H6(3)
27.2%
47.1%
27.5%


0.1 bar 28
C2H4(1)—C2H6(2)—C3H6(3)
17.6%
37.7%
18.6%


Zeolite H-mordenite at 303.15K
CO2(1)—H2S(2)—C3H8(3)
48.8%
27.0%
19.3%


and 0.007-0.1342 bar 30









For mixed-gas adsorption equilibria for CH4—C2H4—C2H6 ternary and C2H4—C2H6-C3H6 ternary on activated carbon at 323 K and 0.1 bar [28], the eL predictions have the highest ARD's of 47% and 38%, respectively. IAST and gL provide very similar predictions for CH4—C2H4-C2H6 ternary and C2H4—C2H6—C3H6 ternary on activated carbon with ARD around 27% and 18%, respectively. In case of mixed-gas adsorption for CO2— H2S—C3H8 ternary on zeolite H-mordenite at 303.15 K and 0.007-0.1342 bar [30], the IAST predictions yield the highest ARD of nearly 49%, followed by the eL predictions with ARD of 27%. The gL predictions yield the best predictions with ARD of 19%.


Comparison of IAST-aNRTL and generalized Langmuir capabilities


In addition to the comparison with eL and IAST, the gL results have also been checked against our prior investigation for mixed-gas adsorption equilibria based a modified IAST and the original adsorption NRTL model [10]. Specifically, Table 6 reports the RMSEs for IAST-aNRTL and gL for five binary systems: N2—O2 binary at 1.013 bar and at 6.08 bar on LiLSX at 303.15 K [29], H2S—CO2 binary at 0.156 bar, C3H8—H2S binary at 0.081 bar, and C3H8—CO2 binary at 0.41 bar on zeolite H-mordenite at 303.15 K. [30] Overall, the gL results for the mixed-gas adsorption equilibria are slightly better than the IAST-aNRTL results. More importantly, the gL results are thermodynamically consistent while the IAST-aNRTL results are valid at constant spreading pressures even though the adsorption equilibria are at constant system pressures. Moreover, gL does not require the computationally expansive calculations on spreading pressures as required for IAST-aNRTL.









TABLE 6







Regressed IAST-aNRTL and generalized Langmuir parameters


for binary mixed-gas adsorption equilibria













T
P
IAST-aNRTL 10
gL
















Adsorbent
Adsorbates
(K)
(bar)
τ12
RMSE
τ12
RMSE
Ref.


















LiLSX
N2(1)—O2(2)
303.15
1.013
−0.84 ±
0.006
−1.83 ±
0.005
[29]






0.04

0.00



N2(1)—O2(2)
303.15
6.080
−0.84 ±
0.005
−1.83 ±
0.001






0.04

0.00


Zeolite H-
H2S (1)—CO2(2)
303.15
0.156
−1.60 ±
0.020
−2.95 ±
0.013
[30]


Mordenite



0.31

0.00



C3H8(1)—H2S(2)
303.15
0.081
−3.22 ±
0.042
−3.93 ±
0.038






0.02

0.00



C3H8 (1)—CO2(2)
303.15
0.410
−3.77 ±
0.063
−4.87 ±
0.085






0.02

0.00









CONCLUSION

A simple, robust, and thermodynamically consistent adsorption model for pure component adsorption isotherms and mixed-gas adsorption equilibria has been described herein. Treating the adsorbent surface as an integral part of adsorption systems, the proposed generalized Langmuir isotherm accurately correlates and predicts mixed-gas adsorption equilibria without the need to compute “spreading pressure” as required for Ideal Adsorbed Solution Theory. Tested with mixed-gas adsorption equilibria data for ten binary and three ternary systems, the generalized Langmuir isotherm outperforms both extended Langmuir and Ideal Adsorbed Solution Theory in predicting mixed-gas adsorption equilibria. As the generalized Langmuir isotherm tracks the surface loading in terms of occupied and vacant sites surface area fractions and the surface heterogeneity in terms of the adsorbate phase activity coefficients, it is shown for the first-time complete phase diagrams for mixed-gas adsorption equilibria analogous to the TPxy phase diagrams for vapor-liquid equilibria. Additionally, the generalized Langmuir isotherm predictions for mixed-gas adsorption equilibria can be used to elucidate the spreading pressure dependence of adsorbate phase activity coefficients in the context of Ideal Adsorbed Solution Theory.


Some other embodiments of the present invention will now be described with respect to FIGS. 9-11. FIG. 9 is a block diagram of an apparatus, system or computer 900, such as a workstation, laptop, desktop, tablet computer, mainframe, or other single or distributed computing platform suitable for performing the methods described herein. Note that the components can be integrated into a single device or communicably coupled to one another via a network. The apparatus, system or computer 900 includes one or more processors 902, a memory or data storage 904, and one or more communication interfaces or input/output interfaces 906, which can be communicably coupled to one or more output device(s) 908 (e.g., printer, internal or external data storage device, display or monitor, remote database, remote computer, etc.) via a network or communications link 910 (e.g., wired, wireless, optical, etc.). The one or more output device(s) can be integrated into the computer 900 as indicated by the dashed line 912


The apparatus, system or computer 900 can be used to estimate an adsorption equilibria for one or more gases from pure component adsorption isotherms. The one or more processors calculate an adsorption of each gas i on a constant monolayer adsorption surface Ao using generalized Langmuir isotherm equations:







θ
i

=




n
i



A
i



A
o


=



K
i
o



y
i


P




γ
i



γ
ϕ



q
i



+




j
=
1

n





γ
i



q
j




γ
j



q
i





K
j
o



y
j


P












θ
ϕ

=




n
ϕ



A
ϕ



A
o


=

1

1
+




j
=
1

n





γ
ϕ



q
j



γ
j




K
j
o



y
j


P












q
i

=


A
i


A
ϕ






where: θi is an adsorbate phase area fraction covered with the gas i,

    • niAi is an occupied area for the gas i,
    • Kio is an intrinsic adsorption equilibrium constant of the gas i,
    • yi is a gas phase mole fraction of gas i,
    • P is a gas vapor pressure,
    • γi is an activity coefficient of the gas i,
    • γϕ is an activity coefficient of vacant sites,
    • qi is a ratio of an effective area of the gas i (Ai) and an effective area of a phantom molecule ϕ (Aϕ),
    • n is a number of the one or more gases,
    • θϕ is an adsorbate phase vacant site area fraction, and
    • nϕAϕ is a vacant area for the phantom molecule ϕ.


The adsorption of each gas i is provided to the output device 908, and a chemical process or a product is developed using the adsorption of each gas i.


In one aspect, the generalized Langmuir isotherm equations reduce to








n
i



n
i
0




=



K
i
o



y
i


P


1
+




i
=
1

n



K
i
o



y
i


P








when (1) the adsorbate and vacant site effective areas are the same A1=A2= . . . =Ai=Aϕ, or equivalently, the saturation loadings of adsorbates and phantom molecule are same n10=n20= . . . =ni0=nϕ0, and (2) the adsorbate phase activity coefficients are unity γiϕ=1. In another aspect, the one or more gases comprise a mixed gas having two or more components. In another aspect, the constant monolayer adsorption surface comprises activated carbon, LiLSX or Zeolite H-mordenite. In another aspect, the gas i comprises CH4, C2H4, C2H6, C3H6, N2, O2, CO2, H2S, or C3H8. In another aspect, the one or more gasses comprise a mixed gas selected from CH4—C2H4, CH4—C2H6, C2H4—C2H6, C2H4—C2H6, and C2H6—C3H6. In another aspect, the one or more gasses comprise a mixed gas selected from N2 and O2. In another aspect, the one or more gasses comprise a mixed gas selected from H2S—CO2, C3H8—H2S, and C3H8—CO2.



FIG. 10 is a flow chart depicting a computerized method 1000 for estimating an adsorption equilibria for one or more gases from pure component adsorption isotherms. One or more processors, a memory communicably coupled to the one or more processors and an output device communicably coupled to the one or more processors are provided in block 1002. An adsorption of each gas on a constant monolayer adsorption surface is calculated using the one or more processors and the generalized Langmuir isotherm equations (24)-(26) or equation (27) in block 1004, namely:







θ
i

=




n
i



A
i



A
o


=



K
i
o



y
i


P




γ
i



γ
ϕ



q
i



+




j
=
1

n





γ
i



q
j




γ
j



q
i





K
j
o



y
j


P












θ
ϕ

=




n
ϕ



A
ϕ



A
o


=

1

1
+




j
=
1

n





γ
ϕ



q
j



γ
j




K
j
o



y
j


P












q
i

=


A
i


A
ϕ






where: θi is an adsorbate phase area fraction covered with the gas i,

    • niAi is an occupied area for the gas i,
    • Kio is an intrinsic adsorption equilibrium constant of the gas i,
    • yi is a gas phase mole fraction of gas i,
    • P is a gas vapor pressure,
    • γi is an activity coefficient of the gas i,
    • γϕ is an activity coefficient of vacant sites,
    • qi is a ratio of an effective area of the gas i (Ai) and an effective area of a phantom molecule ϕ (Aϕ),
    • n is a number of the one or more gases,
    • θϕ is an adsorbate phase vacant site area fraction, and
    • nϕAϕ is a vacant area for the phantom molecule ϕ.


The adsorption of each gas is provided to the output device in block 1006, and a chemical process or product is developed using the adsorption of each gas in block 1008. The method 1000 can be implemented by the apparatus 900 or by a non-transitory computer readable medium encoded with a computer program for execution by a processor that performs the steps of the method 1000.


In one aspect, the generalized Langmuir isotherm equations reduce to








n
i


n
i
0


=



K
i
o



y
i


P


1
+




i
=
1

n



K
i
o



y
i


P








when (1) the adsorbate and vacant site effective areas are the same A1=A2= . . . =Ai=Aϕ, or equivalently, the saturation loadings of adsorbates and phantom molecule are same n10=n20= . . . =ni0?=nϕ0, and (2) the adsorbate phase activity coefficients are unity γiϕ=1. In another aspect, the one or more gases comprise a mixed gas having two or more components. In another aspect, the constant monolayer adsorption surface comprises activated carbon, LiLSX or Zeolite H-mordenite. In another aspect, the gas i comprises CH4, C2H4, C2H4, C2H6, N2, O2, CO2, H2S, or C3H8. In another aspect, the one or more gasses comprise a mixed gas selected from CH4—C2H4, CH4—C2H6, C2H4—C2H6, C2H4—C3H6, and C2H6—C3H6. In another aspect, the one or more gasses comprise a mixed gas selected from N2 and O2. In another aspect, the one or more gasses comprise a mixed gas selected from H2S—CO2, C3H8—H2S, and C3H8—CO2.



FIG. 11 is a flow chart depicting a method 1100 of adsorbing one or more gases. A vessel containing a constant monolayer adsorption surface is provided in block 1102. One or more gases are introduced into the vessel in block 1104, wherein the adsorption of each gas on the constant monolayer adsorption surface is determined by the generalized Langmuir isotherm equations (24)-(26) or equation (27), namely:







θ
i

=




n
i



A
i



A
o


=



K
i
o



y
i


P




γ
i



γ
ϕ



q
i



+




j
=
1

n





γ
i



q
j




γ
j



q
i





K
j
o



y
j


P












θ
ϕ

=




n
ϕ



A
ϕ



A
o


=

1

1
+




j
=
1

n





γ
ϕ



q
j



γ
j




K
j
o



y
j


P












q
i

=


A
i


A
ϕ






where: θi is an adsorbate phase area fraction covered with the gas i,

    • niAi is an occupied area for the gas i,
    • Kio is an intrinsic adsorption equilibrium constant of the gas i,
    • yi is a gas phase mole fraction of gas i,
    • P is a gas vapor pressure,
    • γi is an activity coefficient of the gas i,
    • γϕ is an activity coefficient of vacant sites,
    • qi is a ratio of an effective area of the gas i (Ai) and an effective area of a phantom molecule ϕ (Aϕ),
    • n is a number of the one or more gases,
    • θϕ is an adsorbate phase vacant site area fraction, and
    • nϕAϕ is a vacant area for the phantom molecule ϕ.


      A product can be produced in accordance with the method 1100.


In one aspect, the generalized Langmuir isotherm equations reduce to








n
i



n
i
0




=



K
i
o



y
i


P


1
+




i
=
1

n



K
i
o



y
i


P








when (1) the adsorbate and vacant site effective areas are the same A1=A2= . . . =Ai=Aϕ, or equivalently, the saturation loadings of adsorbates and phantom molecule are same n10=n20= . . . =ni0=nϕ0, and (2) the adsorbate phase activity coefficients are unity γiϕ=1. In another aspect, the one or more gases comprise a mixed gas having two or more components. In another aspect, the constant monolayer adsorption surface comprises activated carbon, LiLSX or Zeolite H-mordenite. In another aspect, the gas i comprises CH4, C2H4, C2H6, C3H6, N2, O2, CO2, H2S, or C3H8. In another aspect, the one or more gasses comprise a mixed gas selected from CH4—C2H4, CH4—C2H6, C2H4—C2H6, C2H4—C3H6, and C2H6—C3H6. In another aspect, the one or more gasses comprise a mixed gas selected from N2 and O2. In another aspect, the one or more gasses comprise a mixed gas selected from H2S—CO2, C3H6—H2S, and C3H8—CO2.


It will be understood that particular embodiments described herein are shown by way of illustration and not as limitations of the invention. The principal features of this invention can be employed in various embodiments without departing from the scope of the invention. Those skilled in the art will recognize or be able to ascertain using no more than routine experimentation, numerous equivalents to the specific procedures described herein. Such equivalents are considered to be within the scope of this invention and are covered by the claims.


All publications and patent applications mentioned in the specification are indicative of the level of skill of those skilled in the art to which this invention pertains. All publications and patent applications are herein incorporated by reference to the same extent as if each individual publication or patent application was specifically and individually indicated to be incorporated by reference.


The use of the word “a” or “an” when used in conjunction with the term “comprising” in the claims and/or the specification may mean “one,” but it is also consistent with the meaning of “one or more,” “at least one,” and “one or more than one.” The use of the term “or” in the claims is used to mean “and/or” unless explicitly indicated to refer to alternatives only or the alternatives are mutually exclusive, although the disclosure supports a definition that refers to only alternatives and “and/or.” Throughout this application, the term “about” is used to indicate that a value includes the inherent variation of error for the device, the method being employed to determine the value, or the variation that exists among the study subjects.


As used in this specification and claim(s), the words “comprising” (and any form of comprising, such as “comprise” and “comprises”), “having” (and any form of having, such as “have” and “has”), “including” (and any form of including, such as “includes” and “include”) or “containing” (and any form of containing, such as “contains” and “contain”) are inclusive or open-ended and do not exclude additional, unrecited features, elements, components, groups, integers, and/or steps, but do not exclude the presence of other unstated features, elements, components, groups, integers and/or steps. In embodiments of any of the compositions and methods provided herein, “comprising” may be replaced with “consisting essentially of” or “consisting of”. As used herein, the term “consisting” is used to indicate the presence of the recited integer (e.g., a feature, an element, a characteristic, a property, a method/process step or a limitation) or group of integers (e.g., feature(s), element(s), characteristic(s), property(ies), method/process steps or limitation(s)) only. As used herein, the phrase “consisting essentially of” requires the specified features, elements, components, groups, integers, and/or steps, but do not exclude the presence of other unstated features, elements, components, groups, integers and/or steps as well as those that do not materially affect the basic and novel characteristic(s) and/or function of the claimed invention.


The term “or combinations thereof” as used herein refers to all permutations and combinations of the listed items preceding the term. For example, “A, B, C, or combinations thereof” is intended to include at least one of: A, B, C, AB, AC, BC, or ABC, and if order is important in a particular context, also BA, CA, CB, CBA, BCA, ACB, BAC, or CAB. Continuing with this example, expressly included are combinations that contain repeats of one or more item or term, such as BB, AAA, AB, BBC, AAABCCCC, CBBAAA, CABABB, and so forth. The skilled artisan will understand that typically there is no limit on the number of items or terms in any combination, unless otherwise apparent from the context.


As used herein, words of approximation such as, without limitation, “about”, “substantial” or “substantially” refers to a condition that when so modified is understood to not necessarily be absolute or perfect but would be considered close enough to those of ordinary skill in the art to warrant designating the condition as being present. The extent to which the description may vary will depend on how great a change can be instituted and still have one of ordinary skill in the art recognize the modified feature as still having the required characteristics and capabilities of the unmodified feature. In general, but subject to the preceding discussion, a numerical value herein that is modified by a word of approximation such as “about” may vary from the stated value by at least ±1, 2, 3, 4, 5, 6, 7, 10, 12 or 15%.


All of the compositions and/or methods disclosed and claimed herein can be made and executed without undue experimentation in light of the present disclosure. While the compositions and methods of this invention have been described in terms of preferred embodiments, it will be apparent to those of skill in the art that variations may be applied to the compositions and/or methods and in the steps or in the sequence of steps of the method described herein without departing from the concept, spirit and scope of the invention. All such similar substitutes and modifications apparent to those skilled in the art are deemed to be within the spirit, scope and concept of the invention as defined by the appended claims.


To aid the Patent Office, and any readers of any patent issued on this application in interpreting the claims appended hereto, applicants wish to note that they do not intend any of the appended claims to invoke paragraph 6 of 35 U.S.C. § 112, U.S.C. § 112 paragraph (f), or equivalent, as it exists on the date of filing hereof unless the words “means for” or “step for” are explicitly used in the particular claim.


For each of the claims, each dependent claim can depend both from the independent claim and from each of the prior dependent claims for each and every claim so long as the prior claim provides a proper antecedent basis for a claim term or element.


REFERENCES



  • 1. Garside, M. U.S. Chemical industry-Statistics & Facts; Statista: 2020.

  • 2. Sholl, D. S.; Lively, R. P., Seven chemical separations to change the world. Nature 2016, 532 (7600), 435-7.

  • 3. Murnen, H. Separation Processes with Hope for Purification Energy Efficiency. www.chemical-materials.elsevier.com/chemical-manufacturing-excellence/separation-processes-for-energy-efficiency/.

  • 4. Laboratory, O. R. N. Materials for Separation Technologies: Energy and Emission Reduction Opportunities; 2005.

  • 5. Lively, R. P.; Realff, M. J., On thermodynamic separation efficiency: Adsorption processes. AIChE Journal 2016, 62 (10), 3699-3705.

  • 6. Shao, L., Grand Challenges in Emerging Separation Technologies. Frontiers in Environmental Chemistry 2020, 1.

  • 7. Yang, R. T., Gas separation by adsorption processes. Butterworth-Heinemann: 2013.

  • 8. Walton, K. S.; Sholl, D. S., Predicting multicomponent adsorption: 50 years of the ideal adsorbed solution theory. AIChE Journal 2015, 61 (9), 2757-2762.

  • 9. Walton, K. S., 110th Anniversary: Commentary: Perspectives on Adsorption of Complex Mixtures. Industrial & Engineering Chemistry Research 2019, 58 (37), 17100-17105.

  • 10. Tun, H.; Chen, C. C., Prediction of mixed-gas adsorption equilibria from pure component adsorption isotherms. AIChE Journal 2020, 66 (7).

  • 11. Suwanayuen, S.; Danner, R. P., Vacancy solution theory of adsorption from gas mixtures. AIChE Journal 1980, 26 (1), 76-83.

  • 12. Wang, Z.; Li, Y.; Guo, P.; Meng, W., Analyzing the Adaption of Different Adsorption Models for Describing the Shale Gas Adsorption Law. Chemical Engineering & Technology 2016, 39 (10), 1921-1932.

  • 13. Markham, E. C.; Benton, A. F., The Adsorption of Gas Mixtures by Silica. Journal of the American Chemical Society 1931, 53 (2), 497-507.

  • 14. Myers, A. L.; Prausnitz, J. M., Thermodynamics of mixed-gas adsorption. AIChE Journal 1965, 11 (1), 121-127.

  • 15. Sorial, G. A.; Granville, W. H.; Daly, W. O., Adsorption equilibria for oxygen and nitrogen gas mixtures on 5 A molecular sieves. Chemical Engineering Science 1983, 38 (9), 1517-1523.

  • 16. Do, D. D., Adsorption Analysis: Equilibria and Kinetics. Imperial college press London: 1998; Vol. 2.

  • 17. Wu, Q.; Zhou, L.; Wu, J.; Zhou, Y., Adsorption Equilibrium of the Mixture CH4+N2+H2 on Activated Carbon. Journal of Chemical & Engineering Data 2005, 50 (2), 635-642.

  • 18. Tien, C., Introduction to adsorption: Basics, analysis, and applications. Elsevier: 2018.

  • 19. Hermann, M.; Niessen, W.; Karge, H. G., Sorption Kinetics of n-Paraffins in MFI-type Zeolites Investigated by micro-FTIR. In Fundamentals of Adsorption: Proceedings of the Fifth International Conference on Fundamentals of Adsorption, LeVan, M. D., Ed. Springer US: Boston, MA, 1996; pp 377-384.

  • 20. Harlick, P. J. E.; Tezel, F. H., Adsorption of carbon dioxide, methane and nitrogen: pure and binary mixture adsorption for ZSM-5 with SiO2/Al2O3 ratio of 280. Separation and Purification Technology 2003, 33 (2), 199-210.

  • 21. Kaur, H.; Tun, H.; Sees, M.; Chen, C.-C., Local composition activity coefficient model for mixed-gas adsorption equilibria. Adsorption 2019, 25 (5), 951-964.

  • 22. Chang, C.-K.; Tun, H.; Chen, C.-C., An activity-based formulation for Langmuir adsorption isotherm. Adsorption 2019, 26 (3), 375-386.

  • 23. Tun, H.; Chen, C.-C., Isosteric heat of adsorption from thermodynamic Langmuir isotherm. Adsorption 2021.

  • 24. Talu, O.; Li, J.; Myers, A. L., Activity coefficients of adsorbed mixtures. Adsorption 1995, 1 (2), 103-112.

  • 25. Talu, O.; Myers, A. L., Rigorous thermodynamic treatment of gas adsorption. AIChE Journal 1988, 34 (11), 1887-1893.

  • 26. Sing, K. S. W., Assessment of Surface Area by Gas Adsorption. In Adsorption by Powders and Porous Solids, Rouquerol, F.; Rouquerol, J.; Sing, K. S. W.; Llewellyn, P.; Maurin, G., Eds. Academic Press: Oxford, 2014; pp 237-268.

  • 27. Britt, H. I.; Luecke, R. H., The Estimation of Parameters in Nonlinear, Implicit Models. Technometrics 1973, 15 (2), 233-247.

  • 28. Costa, E.; Calleja, G.; Marron, C.; Jimenez, A.; Pau, J., Equilibrium adsorption of methane, ethane, ethylene, and propylene and their mixtures on activated carbon. Journal of Chemical & Engineering Data 1989, 34 (2), 156-160.

  • 29. Wu, C.-W.; Kothare, M. V.; Sircar, S., Equilibrium Adsorption Isotherms of Pure N2 and O2 and Their Binary Mixtures on LiLSX Zeolite: Experimental Data and Thermodynamic Analysis. Industrial & Engineering Chemistry Research 2014, 53 (17), 7195-7201.

  • 30. Talu, O.; Zwiebel, I., Multicomponent adsorption equilibria of nonideal mixtures. AIChE Journal 1986, 32 (8), 1263-1276.

  • 31. Costa, E.; Sotelo, J. L.; Calleja, G.; Marrón, C., Adsorption of binary and ternary hydrocarbon gas mixtures on activated carbon: Experimental determination and theoretical prediction of the ternary equilibrium data. AIChE Journal 1981, 27 (1), 5-12.

  • 32. Bezus, A.; Dreving, V.; Kiselev, A., Isotherms and heats of adsorption of propane and propylene on graphitized carbon black-energy of adsorption forces. Kolloidnyi Zhurnal 1961, 23 (4), 389-&.

  • 33. Livingston, H., The cross-sectional areas of molecules adsorbed on solid surfaces. Journal of colloid science 1949, 4 (5), 447-458.

  • 34. Cintas, E. N. Technoeconomic assessment of thermochemical biorefinery based on DME as a platform chemical. University of Seville, 2016.

  • 35. Harris, B. L.; Emmett, P., Adsorption Studies. Physical Adsorption of Nitrogen, Toluene, Benzene, Ethyl Iodide, Hydrogen Sulfide, Water Vapor, Carbon Disulfide, and Pentane on Various Porous and Non-Porous Solids. The Journal of Physical Chemistry 1949, 53 (6), 811-825.

  • 36. Renon, H.; Prausnitz, J. M., Local compositions in thermodynamic excess functions for liquid mixtures. AIChE Journal 1968, 14 (1), 135-144.

  • 37. Scott, R. L., Corresponding States Treatment of Nonelectrolyte Solutions. The Journal of Chemical Physics 1956, 25 (2), 193-205.


Claims
  • 1. A computerized method for estimating an adsorption equilibria for one or more gases from pure component adsorption isotherms comprising: providing one or more processors, a memory communicably coupled to the one or more processors and an output device communicably coupled to the one or more processors;calculating, using the one or more processors, an adsorption of each gas i on a constant monolayer adsorption surface Ao using generalized Langmuir isotherm equations:
  • 2. The method of claim 1, wherein the generalized Langmuir isotherm equations reduce to
  • 3. The method of claim 1, wherein the one or more gases comprise a mixed gas having two or more components.
  • 4. The method of claim 1, wherein the constant monolayer adsorption surface comprises activated carbon, LiLSX or Zeolite H-mordenite.
  • 5. The method of claim 1, wherein the gas i comprises CH4, C2H4, C2H6, C3H6, N2, O2, CO2, H2S, or C3H8.
  • 6. The method of claim 1, wherein the one or more gasses comprise a mixed gas selected from CH4—C2H4, CH4—C2H6, C2H4—C2H6, C2H4—C3H6, and C2H6—C3H6.
  • 7. The method of claim 1, wherein the one or more gasses comprise a mixed gas selected from N2 and O2.
  • 8. The method of claim 1, wherein the one or more gasses comprise a mixed gas selected from H2S—CO2, C3H8—H2S, and C3H8—CO2.
  • 9. The chemical process or the product developed in accordance with claim 1.
  • 10. A system for estimating an adsorption equilibria for one or more gases from pure component adsorption isotherms comprising: at least one input/output interface;a data storage;one or more processors communicably coupled to the at least one input/output interface and the data storage, wherein the one or more processors calculate an adsorption of each gas i on a constant monolayer adsorption surface Ao using generalized Langmuir isotherm equations:
  • 11. The system of claim 10, wherein the generalized Langmuir isotherm equations reduce to
  • 12. The system of claim 10, wherein the one or more gases comprise a mixed gas having two or more components.
  • 13. The system of claim 10, wherein the constant monolayer adsorption surface comprises activated carbon, LiLSX or Zeolite H-mordenite.
  • 14. The system of claim 10, wherein the gas i comprises CH4, C2H4, C2H6, C3H6, N2, O2, CO2, H2S, or C3H8.
  • 15. The system of claim 10, wherein the one or more gasses comprise a mixed gas selected from CH4—C2H4, CH4—C2H6, C2H4—C2H6, C2H4—C3H6, and C2H6—C3H6.
  • 16. The system of claim 10, wherein the one or more gasses comprise a mixed gas selected from N2 and O2.
  • 17. The system of claim 10, wherein the one or more gasses comprise a mixed gas selected from H2S—C2, C3H8—H2S, and C3H8—CO2.
  • 18. A computer program embodied on a non-transitory computer readable storage medium that is executed using one or more processors for estimating an adsorption equilibria for one or more gases from pure component adsorption isotherms comprising: a code segment for calculate an adsorption of each gas i on a constant monolayer adsorption surface Ao using generalized Langmuir isotherm equations:
  • 19. A method of adsorbing one or more gases comprising: providing a vessel containing a constant monolayer adsorption surface Ao;introducing the one or more gasses into the vessel; andwherein the adsorption of each gas i on the constant monolayer adsorption surface Ao is determined by generalized Langmuir isotherm equations:
  • 20. The method of claim 19, wherein the generalized Langmuir isotherm equations reduce to
  • 21. The method of claim 19, wherein the one or more gases comprise a mixed gas having two or more components.
  • 22. The method of claim 19, wherein the constant monolayer adsorption surface comprises activated carbon, LiLSX or Zeolite H-mordenite.
  • 23. The method of claim 19, wherein the gas i comprises CH4, C2H4, C2H6, C3H6, N2, O2, CO2, H2S, or C3H8.
  • 24. The method of claim 19, wherein the one or more gasses comprise a mixed gas selected from CH4—C2H4, CH4—C2H6, C2H4—C2H6, C2H4—C3H6, and C2H6—C3H6.
  • 25. The method of claim 19, wherein the one or more gasses comprise a mixed gas selected from N2 and O2.
  • 26. The method of claim 19, wherein the one or more gasses comprise a mixed gas selected from H2S—C2, C3H8—H2S, and C3H8—CO2.
  • 27. A product produced in accordance with claim 19.
CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims priority to U.S. Provisional Application Ser. No. 63/226,303, filed Jul. 28, 2021, the entire contents of which are incorporated herein by reference. This application is related to PCT Application Serial No. PCT/US20/45586, filed Aug. 10, 2020, which claims priority to U.S. Provisional Application Ser. No. 62/860,319, filed Jun. 12, 2019. The entire contents of the foregoing applications are incorporated herein by reference.

STATEMENT OF FEDERALLY FUNDED RESEARCH

This invention was made with government support under DE-EE0007888 awarded by the U.S. Department of Energy. The government has certain rights in the invention.

PCT Information
Filing Document Filing Date Country Kind
PCT/US2022/037972 7/22/2022 WO
Provisional Applications (1)
Number Date Country
63226303 Jul 2021 US