Kelvin’s Dissymmetric Models and Consistency Conditions of Multicomponent Gas-Liquid Equilibrium and Capillary Condensation

To describe phase equilibriumofmixtures, we develop a nonclassical approach based onusing different equations of state for gas and liquid.We show that not all the types of EOS are admissible but only those which verify some specific conditions of consistency.We developed themathematical theory of this new approach for pure cores and formixtures, in presence and absence of capillary forces, which leads to explicit analytical relationships for phase concentrations of chemical components. Several examples of comparison with experimental data for binary and ternary mixtures illustrate the feasibility of the suggested approach.


Introduction
The objective of the present paper is to develop simplified mathematical models of phase equilibria which would significantly reduce computation time, while retaining an acceptable level of precision.This problem is highly pertinent for multicomponent two-phase or three-phase fluids.The results of such simulations provide information about the chemical composition of each phase at any pressure and temperature and the fractional volume of each phase.
Modern thermodynamic models are essentially based on cubic equations of state (EOS) such as those developed by Van-der-Waals, Peng-Robinson, or Redlich-Kwong [1][2][3][4].They ensure high quality in descriptions of gas-liquid coexistence, which determines the success of the modern thermodynamic theory.At the same time, when such models are used to simulate fluid flow in geological oil, gas, or water reservoirs, more than the half of the computation time is spent on the thermodynamic block.The iterative numerical procedure suggested in [5,6] in a form called "the flash calculation" must be done at each time step and at each space point which is highly time-consuming.
Consequently a lot of papers have been published with the objective of reducing the time of flash calculation by reducing the number of primary variables.Michelsen in [7] reduced the number of parameters to three and accepted zero values for all binary interaction parameters (BIPs).By only taking into account the BIPs between a single component and the system of all other components, the number of parameters was reduced to five in [8].More recently, in [9] the authors suggested the split of the BIPs into two parameters using a simple quadratic expression.
Unfortunately the simplest approaches based on empirical Henry's or Raoulth's laws only permit the calculation of the composition of one phase and are not adapted for use with multiphase multicomponent systems.
Another way to reduce the computational time for the general coupled hydrodynamic-thermodynamic problem consists of splitting the thermodynamics and hydrodynamics.An example of such splitting is published in [10,11], where the authors developed the so-called HT-splitting, an asymptotic procedure providing nonclassic thermodynamic equations completely independent of hydrodynamic equations.This 2 Journal of Thermodynamics approach has lead to a new kind of thermodynamic model made up of traditional algebraic relationships and new differential thermodynamic equations.
Completely different technique of splitting the thermodynamics and hydrodynamics was developed in [12] for mono-dimensional two-phase -component flow.By a specific replacement of variables, the  ×  conservation laws are transformed into a reduced ( − 1) × ( − 1) closed auxiliary system containing just thermodynamic variables and one lifting hydrodynamic equation containing the phase saturation.
A cubic EOS has to be excessively complicated in order to be capable of capturing not only the single-phase behaviour of gas and liquid but also the transition two-phase zone.However, just within this zone the cubic EOS is not used: its complicated nonmonotonic behaviour is simply replaced by an isobaric straight line (Maxwell's line) which corresponds to the true gas-liquid equilibrium.Then our logic was the following: it is not necessary to require that an EOS would describe the transition zone; then it is possible to describe liquid and gas by two different EOS, which may represent very simple equations.This becomes the significant source of simplifying equilibrium equations.
The system of equilibrium equations and two different EOS for gas and liquid will be called the dissymmetric equilibrium model.
The approach based on dissymmetric EOS is not new and is frequently applied in chemistry for various purposes.A lot of other papers exist when the dissymmetric models are applied to calculate some thermodynamic parameters of mixtures.Their application to calculate phase equilibria is less known.This was done first probably by Kelvin, who used the model of the ideal gas and an incompressible liquid to describe the equilibrium controlled by capillary forces for a monocomponent fluid.The dissymmetric EOS for two-component fluids was applied in [13] to calculate the two-phase envelop.In [14] the generalization of the Kelvin equation was obtained for multicompositional mixtures by applying the real-gas EOS and an incompressible liquid and by introducing the new concept of a "mixed volume."The resulting equation is valid only near the dew point.
Another kind of dissymmetric models, known as the phigamma approach, has been used for decades in chemical engineering thermodynamics [15].In this approach vapor is described by an EOS, while liquid is described by the activity coefficient.The model used for the activity coefficient represents some semiempirical approximations but never an EOS.
The main difference of the present paper from all other papers in which a kind of dissymmetric approaches was applied consists of the following: (i) Within the present paper we use two different EOS for gas and liquid to calculate phase equilibria.(ii) We have shown that not any couple of EOS is consistent with the two-phase equilibrium.(iii) We have determined the general theoretical constraints to the choice of two EOS consistent with the gas-liquid equilibrium.
In all the papers which used the dissymmetric approach the selection of the individual EOS is arbitrary and not constrained.However, once we tried to develop dissymmetric EOS models and to extend them to multicomponent mixtures, we found that several EOS lead to nonexistence of two-phase states.In particular, it is possible to show that the incompressible liquid is unable to capture two-phase states whatever the EOS for gas is.
The capacity of an EOS to capture two-phase states is called consistency conditions.Thus, the selection of the individual dissymmetric EOS for gas and liquid is not arbitrary but must satisfy the consistency conditions which represent mathematically the conditions of solution existence for the equilibrium equations.
The development of the consistency conditions and their use in order to determine the adjustable parameters of the dissymmetric EOS is the second key element of the paper.
The developed approach gives a form of equilibrium equations which can be solved analytically with respect to phase concentrations even for multicomponent nonideal systems, which significantly reduces the computational time.

Principle Idea: One-Component Fluid
The best way of presenting the principles of our approach is to analyze a simple fluid consisting of one chemical component.

Main Equilibrium
Relationship.The gas-liquid equilibrium may be described using the system of one equilibrium equation, (1a), and two equations of state, (1b) and (1c), which have the following form: where  is the pressure,  is the temperature,  is the molar volume,  is the chemical potential,   and   are two arbitrary reference values of pressure,  eq is the equilibrium pressure at given , and indexes  and  correspond to gas and liquid, respectively.
The definition of the chemical potential, the fluid volume, and the molar volume through the Gibbs energy  is where  is the number of moles of fluid.
Then, for any , (1a) represents the nonlinear system of one equation with respect to  eq which is the pressure of gasliquid coexistence.
The classic approach consists of using the same EOS for all fluids, whereas our approach consists of using different EOS for gas and for liquid.

Pseudoliquid and Pseudogas.
Within the framework of the suggested approach, it is necessary to introduce the following definitions.
The pseudoliquid and the pseudogas are the hypothetical fluids which are authorized to exist at any  and , with each of them being defined only by a formal EOS, (1b) or (1c).
This means that both pseudophases can coexist at any  and , but this kind of coexistence is not proved by any equilibrium conditions.Figure 1 presents two curves, ADE and ABM, which are the graphs of the EOS for pseudogas and pseudoliquid, respectively.
The true coexisting gas and liquid are the fluids defined by the system of two EOS, (1b) and (1c), and the equilibrium equation, (1a), in such a way that this system of three equations can only have one solution for the coexistence pressure  eq and molar volumes  ,eq and  ,eq (at a fixed temperature ) which satisfy the strong inequality: In Figure 1, this corresponds to the horizontal line of Maxwell , which is characterized by three parameters:  eq ,  ,eq , and  ,eq .
The true gas and liquid are the fluids which are determined only by the EOS and exist outside the two-phase coexistence domain; that is, the true gas corresponds to   ≥  ,eq , while the true liquid corresponds to   ≤  ,eq .In Figure 1, this corresponds to the segment DE for the true gas and BA for the true liquid.
Both the true phases and pseudophases are described by the same EOS.This means that a pseudoliquid may be considered as the continuous prolongation of the true liquid into the domain of molar volumes greater than  ,eq , while the pseudogas is the prolongation of the true gas into the domain of volumes lower than  ,eq .

2.
3. Pseudocritical Point.Two curves () described by two independent EOS can cross one another at a point (), which corresponds to pressure  * .At this point the pressure, temperature, and volume of gas and liquid are trivially identical, since the difference between gas and liquid disappears.Consequently, the chemical potentials in  are identical.This can be proved more strictly in three different ways.
Thermodynamic Method.At point  we have two substances, for which all the parameters of state are identical: pressure, volume, temperature, and composition.Then each function of state is also identical for them.The chemical potential is the internal energy reported to one mole, and the internal energy is the function of state.Then the chemical potentials of two substances are identical.
Molecular Physics Method.At point  we have two substances, which have identical molecules (the same composition), identical distances between the molecules (identical pressure and volume), and identical fluctuation energy (identical temperature).Then the internal energy of these two substances is identical.Then their chemical potentials are also identical.
Mathematical Method.Let us use the property of homogeneity of the Gibbs energy.Indeed, the Gibbs free energy depends on fluid pressure , temperature , volume , and the number of moles :  = (, , , ).Being an additive function of the number of particles, the energy satisfies the following homogeneity condition: (, , , ) = (, , , ), for any .Selecting  = 1/, we obtain The Euler theorem for homogeneous functions says that It follows from ( 7) and ( 8) that  = (, , ).Thus, for fixed values of , , , the chemical potential of a fluid, , is a constant value, independently of its phase state.Then we have the following at point  * :   =   .
Respectively, point  * will be called pseudocritical.Pressure  * () is the pseudocritical pressure which is defined as the solution of the following equation: In general case, two monotonically decreasing curves can cross each other at several points:  I * ,  II * , . ... All of them will be called pseudocritical.
If we select for the reference pressures,   =   =  * , then the chemical potentials   ( * , ) and   ( * , ) must be identical.Consequently, the equilibrium equation ((1a)-(1c)) becomes which is the algebraic equation with respect to the equilibrium pressure  eq .The geometrical meaning of the left-hand side is the surface between two curves   () and   () (at a fixed ), calculated between the pseudocritical pressure  * and the pressure of gas-liquid coexistence  eq .

Application to the Case of a Uniform EOS for Both
Phases.In the traditional case, the EOS is identical both for liquid and for gas.Then any point corresponds to (9) and is pseudocritical.It is easy to show that in this case equilibrium equation (10) is independent of the selection of  * .
Let  * be an arbitrary reference pressure located far from the two-phase zone.Let us transform the integration ∫ () in (10) to ∫ ().The integration by parts of (10) which is the system of two nonlinear equations for two equilibrium molar volumes:  ,eq and  ,eq .For a cubic EOS, as the Van-der-Waals equation,  = /( − ) − / 2 , with two empirical coefficients,  and , a curve () below the critical point has the form shown in Figure 2. Then the existence of the nontrivial solution of the first equation in (12) may be interpreted graphically as shown in Figure 2. The left-hand side of this equation is the rectangular area .The right-hand side obviously corresponds to the same area if the grey areas  and CNB eliminate one another, that is, are identical, which determines the unique location of the Maxwell line  and the equilibrium pressure  eq .This graphical technique is known as Maxwell's equal area rule.
The Maxwell rule determines the mathematical necessary condition on the behaviour of the curve () which ensures the existence of two-phase states: it must be nonmonotonic with  minimums and  maximums, where  ≥ 1.Only in this case does system (12) have a nontrivial solution.This property of the curve () represents the condition of consistency between the EOS and the equilibrium equation.Obviously, this property is verified for a cubic EOS with respect to volume.

Consistency Conditions.
For a dissymmetric model, the technique of application of the equilibrium equation (10) and the consistency condition are different and are sufficiently based on using the pseudocritical pressure  * .
The existence of two-phase state imposes some constraints on the form of two individual EOS which should ensure the solvability of (10) with respect to  eq .We will call these constraints the solvability conditions of the equilibrium equation or consistency conditions between the equilibrium equation and EOS.
The necessary conditions of solvability may be formulated in a geometrical way, taking into account the fact that the integral in (10) represents the surface between two curves   () and   () (at a fixed ), as shown in Figure 3(a).
In order to be zero, this surface has to have at least one positive part and one negative part and the areas of these two parts must be identical.Thus, we obtain the necessary and sufficient criteria of solvability or the consistency conditions: The first condition means that the curves   () and   () cross each other twice (points  and   in Figure 3).This is necessary to obtain two segments with different signs of the area between these two curves.
The second condition (which is valid only for nontrue gas and liquids) determines the correct placing of these two curves with respect to one another in order to ensure inequality of (6) for true gas and liquid.
Then the equilibrium equation ( 10) takes the more exact form: If the conditions ((13a) and (13b)) are satisfied, then the equilibrium pressure  eq is determined from (14) in such a way that the negative and positive areas of the surface between the curves   () and   () are identical between points  II * and  eq .
In general case, two curves   () and   () can cross each other at a number of points larger than two.But to ensure solvability, it is sufficient to have only two points.
From this graphical technique it follows that the twophase equilibrium is impossible for an incompressible liquid.Indeed, any monotonic function   () cannot cross twice the vertical line   = const.
From these principles, we deduced the following graphical technique which may be easily used to detect the consistency (or inconsistency) of the selected system of EOS and to determine the equilibrium pressure: (i) Two EOS whose plots  versus  satisfy the following conditions need to be selected: (a) both functions () for gas and liquid should be monotonic at each point ; (b) they cross each other at two different pressures  I * and  II * (Figure 3); (c) the curve () for pseudogas is situated on the left from that of pseudoliquid between  I * and  II * , according to condition (13b).(ii) The equilibrium pressure  eq is determined graphically in such a way that the surface of the area between two curves between points  I * and  eq (Figure 3) is identical to the area of the surface between points  I * and  II * .In the classical case considered in the previous section the equilibrium equations ( 12) are solved with respect to two different volumes  ,eq and  ,eq .In contrast to it, the solution of ( 10) is searched directly with respect to the unique value of pressure  eq .
The graphical areas which eliminate one another are located along axis  (Figure 3(a)) and not axis , which represents reversed Maxwell's rule; that is, the analogue of the Maxwell rule of equal areas turned by 90 ∘ counterclockwise.

Consistent Exponential EOS for Liquid and Ideal EOS
for Gas.Another model which satisfies the consistency conditions is that of the ideal EOS for gas and the exponential EOS for liquid, which represents a theoretical equation based on the assumption that the compressibility coefficient is independent of pressure: where  is the universal gas constant value,  0 is an arbitrary constant pressure value (different from   in (1a)-(1c) and ( 5)), while  0 () and () are two empirical parameters of the model which determine the behaviour of the phase diagram.Parameter  has the meaning of the isothermal compressibility coefficient.Indeed, the exponential equation results immediately from the definition of the compressibility coefficient:  ≡ −(1/)(/), if we assume that  = (), which is valid far from the critical point.Let us prove that the selected model ( 15) is self-consistent, if the parameters verify the following conditions: Indeed, (13a) takes the following dimensionless form through the dimensionless pressure  ≡ / ⊕ : The plot of the left-hand and right-hand sides is shown in Figure 4.
Two solutions  I * and  II * transform to a single value if the exponential function  1 =   touches the straight line  2 = , which means the identical derivatives of  1 and  2 :   = 1.This equation and (17) determine two critical values of  and : which corresponds to the upper estimation in (16).For larger values of , two curves do not cross each other, which means that ( 17) has no solutions.
For smaller values of , the curve  1 =   crosses the line  2 =  twice at two points  I * and  II * .When  → 0, the first pressure  I * tends to 1 and the second pressure  II * tends to infinity.This yields the lower estimation in (16), which corresponds to the incompressible liquid.Thus, for these values of , condition (13a) is satisfied.It is easy to show that (13b) is also satisfied.
Once the pseudocritical pressure  II * is determined from (17), the equilibrium pressure,  eq , is found as the solution to the following equation resulting from (14): The value of  eq is the nontrivial solution of this equation, different from  II * .So, for given , one can get  II * and  eq .Since we know the gas and liquid EOS we can get  ,eq and  ,eq and, finally, the envelope of the phase diagram (this being the line joining all the dew points and all the bubble points).

General Technique of Adjusting Dissymmetric Models.
As shown in Section 2.6, for ideal gas and exponential EOS for liquid, if we know the values of  and  0 (i.e., the parameter ), then it is possible to determine the equilibrium pressure and volumes.In the similar way it is possible to formulate the symmetrical problem: to determine parameters () and  0 () by using the experimental information on the equilibrium pressure and volumes of pure components.The obtained adjusted dissymmetric models will be used further to describe multicomponent mixtures without any additional fitting (Section 4).
We can extend this approach by using the virial EOS for gas in the form where the new parameter  0 should be adjusted by using the data on the equilibrium gas volume,  ,eq .The virial EOS gives more accurate results for equilibrium gas volume than the ideal gas EOS.
Let us assume that the equilibrium pressure  eq , the equilibrium molar liquid, and gas volumes  ,eq and  ,eq are given as experimental data.For pure components, they can be found in the literature or calculated by Peng-Robinson EOS.
Introducing new dimensionless variables, we can transform system (21) to one nonlinear equation which determines the function P(, ), where  and  are the given parameters, and an explicit relationship for function Q(, ): If the gas phase is ideal, then  = ; then we obtain for P() and Q() instead of ( 23) (25) Proof.
(1) Eliminating  0 and  0 from the last two equations (21), we obtain from the first two equations (2) Let us eliminate Q from the second relationship, which gives (24).Substituting it into the first equation, we deduce (23).
Let us note that the functions P(, ) and Q(, ) are universal and do not depend on individual properties of chemical components.
The solution of (25) obtained numerically is shown in Figure 5.
To select the initial approximation for iterative solution of (25), it is possible to use the coarse approximation obtained by fitting the curve shown in Figure 5: P ≈  0.98 ; then Q ≈ (0.02/(1 −  0.98 )) ln .
Note that the boundary values of these functions which correspond to the formal critical point are Q cr = 1, P cr = 1, and  cr = 1.
In Figure 7, we present the diagrams for several pure components calculated by the dissymmetric model and Peng-Robinson EOS, for  = 300 K.
Remarks.(1) The developed technique of determination of the parameters gives unsatisfactory results for temperatures close to the critical point.For the moment, the dissymmetric models cannot pretend to be a good tool to describe the critical point.Indeed, the true critical point is the inflexion point and the point of zero derivative for the true curve ().Such a complex behaviour is impossible for the point of intersection of two monotonic curves which correspond to pseudogas and pseudoliquid.
(2) In the present paper, we use the ideal EOS for gas just to simplify the principle ideas of the suggested approach.Obviously, the ideal EOS can lead to some errors in calculating the equilibrium gas volume, especially at high equilibrium pressures.The virial EOS (20) gives more adequate results.

Equilibrium of Multicomponent Phases
Let us suppose that the fluid consists of  chemical components which are present in gas and liquid.The gas-liquid equilibrium would now determine the composition of each phase and not just the equilibrium pressure and volumes.
The superscript  will be the number of chemical components; the superscript "pure" refers to the pure component considered at the same pressure and temperature as the mixture; the subscript "eq" refers to the gas-liquid equilibrium.Other notations are as follows: (i)    is the mole fraction of component  in phase , which will be also called "the phase concentration." (ii)   is the mole fraction of component  in the entire fluid, which will be also called "the total concentration." (iii)    is the chemical potential of component  in phase .
(iv)  ,ideal is the chemical potential of component  in the mixture if the mixture behaves as an ideal gas.
(v)  , ideal pure is the chemical potential of the pure component  if it behaves as an ideal gas.
(vi)    is the number of moles of component  in phase .
(vii)   is the volume of phase .
=   (, , Equation (30a) results from the equivalence between chemical potentials; (30b) and (30c) are two EOS for pseudogas and pseudoliquid.The proof of (30a) is given in Appendix A.
The molar volumes    are defined if EOS (30b) and (30c) are given; the individual parameters of pure components as their pseudocritical pressures  II, * are the known values.The consistency of system ((30a)-(30d)) is ensured, if the consistency conditions ((13a) and (13b)) are verified for all pure components.We assume that all individual phase diagrams for pure components are well adjusted to experimental data by applying the technique described in Section 3.
System of  equations (30a) contains 2 variables: 2 − 2 phase concentrations    and    (due to (30d)), pressure , and temperature .Such a system is unclosed and contains  independent variables.We select them as , , and −2 phase concentrations called independent.(This means, in particular, that the pressure in the equilibrium is not uniquely defined and can vary within the two-phase region.) Depending on the closure relationships, we obtain three main problems related to multicomponent equilibria: (i) Flash calculation: at given , , and total composition of the fluid, {  } −1 =1 , it is necessary to calculate the composition of both phases and the molar fraction of liquid phase.
In this case, the equilibrium relationships (30a) are completed with  − 1 material balance equations: where   is the gas molar fraction.
Then the system of 2 − 1 equations (30a) and (31) contains 2 − 1 unknown variables: 2 − 2 phase concentrations and saturation   .Therefore it is closed.(ii) Calculation of the bubble-point curve: at given , we assume that only liquid phase exists, while the gas saturation tends to zero.This means that the total composition is equal to the liquid composition:    =   ( = 1, . . .,  − 1), and the disappearing gas phase has some limit values of phase concentrations.Then the system of  equations (30a) contains  unknown variables:  − 1 limit concentrations    and pressure , which is called the bubble-pressure,  bbl .Its is clear that  bbl is the function of the total composition and temperature.(iii) Calculation of the dew-point curve: At given , we assume that only gas phase exists.Then    =   ( = 1, . . .,  − 1) and so on, similar to bubble point.

Equivalent form of the Equilibrium Equations through the Activity Coefficient.
It is useful to represent the partial volumes of components in mixture through the individual volumes of pure components.This can be done by using the concept of ideal mixing, which assumes that the mixture volume is the sum of the volumes of the pure components:  = ∑  =1  ,pure , where  ,pure is the volume occupied by component  in the pure state at the same , .
From this definition, it is possible to obtain the following property of an ideal phase :    =  ,pure  , for any component  (the definitions of variables are given at the beginning of Section 4).
For nonideal mixture, we can introduce the excess in volume,   , which takes into account the nonideality: (for an ideal mixture,   ≡ 0).According to the definition of the Gibbs energy, volume   is related with the excess in the Gibbs energy with respect to that of an ideal mixture,   :   =   /.Then this gives the following result for molar volume from (32): where we have introduced the new variable called the activity coefficient of component ,   .The value   /  is the excess in the chemical potential, which is usually presented in the logarithmic form through   .
Applying (33) to both phases, we obtain the following from (30a): where    and    are the activity coefficients of component  in gas and liquid.We take into account the fact that both activity coefficients tend to 1 when pressure tends to zero.
Two versions of the equilibrium equations, (30a) and (34), are equivalent, but they assume different use of EOS for liquid and gas: for (30a), we should use EOS to determine the partial volumes, while for (34) the EOS should give us the activity coefficient.The activity coefficient may be also determined through several specific models as NRTL expression (Nonrandom Two-Liquid) [20].The NRTL equation is one of the most popular and efficient methods of describing nonideal mixing.This model assumes that, even in multicomponent mixtures, the interactions between molecules can be described as two to two interactions.As a result, the model for multicomponent system contains only binary fitting parameters.The NRTL equation can efficiently represent behaviour of various mixtures, even those of a highly nonideal nature.For a mixture of  components, the NRTL equation is (in terms of the excess in Gibbs energy) Other parameters are and   are dimensionless adjustable interaction parameters,   is the characteristic potential of  −  interaction, and parameter  is related to the nonrandomness of molecule motion in the mixture which varies from about 0.20 to 0.47.

Multicomponent Phase
Compared with the similar equation for pure components, (14), we see that the presence of new term  disturbs the equilibrium.As explained for (14), the term  represents graphically the sum of one negative area shown in Figure 3(a) and one positive area.Term  is also positive for the light component and augments the contribution of positive areas.To equilibrate their role, the negative area must become higher than for the pure light component.This means that the equilibrium pressure of the mixture will be lower than that of the pure light component.Similarly, it is possible to show that it will be higher than the equilibrium pressure for the pure heavy component.As the equilibrium pressure of a mixture depends on the total concentrations  1 , then it will continuously vary between  2 eq and  1 eq when the total composition of mixture changes.
The use of the same ideal gas model for all gaseous individual components significantly reduces the number of parameters which need to be fitted.Let us verify that such a simplification does not lead to some nonphysical situations.
This case is illustrated in Figure 8 for a two-component fluid.The two red curves correspond to two EOS for two pure components in a pseudoliquid state, while the pseudogas state is described by the same EOS of ideal gas for both pure components (the blue curve).Component 1 is light.
Figures 8(a) and 8(b) illustrate the graphical determination of the equilibrium pressures  1 eq and  2 eq for two pure components.Pressure  1 eq is determined in such a way that the yellow surface  would be equal to the yellow surface .Pressure  2 eq results from the equality between the yellow surfaces  and .The results of these constructions are which satisfy the physical meaning: indeed, the equilibrium pressure must be lower for the heavier component, while the equilibrium molar volumes must by lower for the light component.Thus, the use of the ideal EOS for all pure components does not lead to any physical contradiction.

Equilibrium Equations for Ideal Mixing in Gas and Nonideal Mixing in Liquid.
Let the mixing in the liquid phase be nonideal, while the gas is ideal.Then    =  ,pure  , or    = 1 for all .From (34), we obtain the following equilibrium equations for  = 1, . . ., : Using (35), we obtain the expression of the activity coefficients through the phase concentrations: Figure 8: Example of consistency of the same EOS for gas for individual components: graphical determination of the equilibrium pressures  1 eq (a) and  2 eq (b).

Equilibrium for Virial EOS for Gas and Exponential
EOS for Liquid.Let the pure components be a liquid with exponential EOS and a gas described by virial EOS ( = 1, . . ., ): where  ,pure 0 ,   0 ,   ,  II, * , and   0 are the parameters of the model which determine the behaviour of the envelope of the phase diagram.They are selected to comply with experimental data on pure components and satisfy consistency conditions imposed by the envelope (Section 3).For the ideal gas  ,pure 0 ≡ 0.

Analytical Solutions for Phase Concentrations in Two-
Component Case.For a binary mixture, the equilibrium equations (43) take the following form: where the superscript 2 means the number of components but not "square." System (45) determines two concentrations,  1  and  2  , as the functions of pressure .
This system has an explicit analytical solution for phase concentrations  1  (, ) and  1  (, ) in two cases, (1) the ideal mixing in both phases and (2) the moderately nonideal mixing in the liquid phase: where -values are defined as   ≡    /   .For ideal mixing in both phases, -values are calculated directly from (45) (for  = 1, 2): For nonideal mixing in liquid, let us assume the hypothesis of moderate nonideality: the interactions between different molecules and identical molecules are nearly equal:  12 ≃  21 ≃  11 ≃  22 , which leads to   ≃ 0. Then -values in (46) are calculated as where Coefficients  01 and  02 are calculated by using formulae (47).This solution is developed in Appendix B.

Analytical Solutions for Phase Concentrations in Three-
Component Case.In a three-component system, three equations of phase equilibrium contain four phase concentrations  1  ,  2  ,  1  , and  2  , which means one concentration is independent.We will select  1  as this independent parameter.Then the equilibrium system determines three other concentrations as the functions of the equilibrium pressure, temperature, and  1  .The dimensionless equilibrium equations are As said previously, this system has an explicit analytical solution regarding phase concentrations in two cases, (i) ideal mixing everywhere and (ii) moderately nonideal mixing in liquid: For ideal mixing in both phases, the following was obtained for -values ( = 1, 2, 3): For nonideal mixing in liquid, where   0 is calculated through (40) in which one should use   0 instead of    , for all ;   0 ,   0 , and  0 are calculated as for ideal mixing by using formulae (51) and (52).
As mentioned in Section 4.1, the equilibrium equations can be completed with  material balance equations (31), for given total composition of the fluid   .This makes it possible to determine parameter  1  considered as independent in the equilibrium equations:  where   is the gas molar fraction.We obtain two equations with two unknown variables:  1  and   .The envelopes of the two-phase domain are calculated by using the same equations in which  1  =  1 for the bubble curve (transition to single-phase liquid) and  1  =  1 / 1 for the dew curve (transition to single-phase gas), which results from (54) if   = 0 and   = 1, respectively.

Comparison with Experimental
Data and Cubic EOS 5.1.Two-Component Case.The gas-liquid equilibrium for a two-component mixture at constant temperature is usually represented by diagrams " versus  1 ," where  1 is the total concentration of the lightest component.At high  1 , the system is gaseous; at high , it is liquid, while a two-phase domain exists between these two states.The dissymmetric model of pure components was presented by exponential EOS for liquid and ideal gas.
We used the experimental data for several binary mixtures published in [21]: Parameters of the dissymmetric models of pure components obtained by fitting their - curves with the diagrams calculated by using Peng-Robinson EOS are presented in Table 1.
For mixture, the phase diagrams (the envelops of the twophase zone) were constructed by using the model presented in this paper in terms of  versus the molar volume  (see Figure 9 The comparison of calculations by using the dissymmetric model with experimental data for Butane (C 4 H 10 )-Decane (C 10 H 22 ) mixture is presented in Figure 10.
The case of nonideal mixing is illustrated in Figure 11 for the mixture Pentane-Benzene, for which there exist experimental data.As seen, the nonideal dissymmetric model gives satisfactory results with respect to real data.
We underline that, for mixtures, only the parameters of pure components should be fit.

Three-Component
Case.The gas-liquid equilibrium of a three-component mixture at constant temperature and pressure is often represented in an equilateral triangle (ternary diagram).Its height is equal to unity.A mixture composition is represented by a point within the triangle, so that the molar fractions of three components are given by three heights lowered from that point on three sides.
For each total composition ( 1 ,  2 ,  3 ), phase concentrations can be calculated from (51)-(54).We drew the bubble curve and the dew curve by using the following relationships resulting from (51)-( 54): (i) For bubble point curve, (ii) For dew point curve, We have calculated the ternary diagrams for several mixtures by using the dissymmetric model and compared them to experimental data existing in the literature or to other methods of calculation.Parameters of the dissymmetric model obtained by fitting the phase equilibrium of pure components are presented in Table 2.
The first three Figures 12,13 Remark.To use the dissymmetric model, it is necessary to fit only the parameters of pure components, by using the procedure described in Section 3. The consecutive calculation of phase composition is performed by analytical formulae.Therefore, once the dissymmetric model has been adjusted for pure components, it calculates the multicomponent equilibrium instantaneously.
The time of adjustment for a pure component by using algorithm presented in Section 3 is lower than 2 s (for various processors), if the equilibrium pressure and volumes are known (they can be found in various source data).

Equilibrium Equations in Presence of Capillary Pressure.
In the situation when the interface between gas and liquid is curvilinear, the equilibrium gas and liquid pressures are different due to the capillary pressure calculated according to the Laplace equation: where   is the capillary pressure,  is the surface tension, and  is the mean curvature radius of the gas-liquid interface.
where  is the gas pressure.
For virial EOS for gas and exponential EOS for liquid, we obtain, similar to (43), where Ca ≡   / 1 0 is the capillary number.
Due to the presence of the new term (with respect to (10)), in this case consistency can be reached even if we have only one point of intersection between the EOS for pseudoliquid and pseudogas.This for instance is the case of incompressible liquid and ideal gas, for which the vertical line   = const intersects the hyperbolic curve of the ideal gas only at one

Journal of Thermodynamics
where  ≡  eq / I * ; the capillary number Ca ≡ 2/( I * ).This equation has two nontrivial solutions, as shown in Figure 18, where the left-and right-hand sides of (62) are presented.They correspond to points  and  (point  = 1 on the horizontal axis corresponds to  =  I * ).Point  has no physical meaning, as the corresponding value   > 1; then  eq >  I * .This means that the molar volume of the coexisting liquid is higher than that of gas ( ,eq = / eq ), which is physically impossible.So the true solution corresponds to point .We obtain the characteristic pressure   which is lower than  I * , at which the true gas and liquid can coexist.
According to Kelvin, for a plan gas-liquid interface, the capillary pressure is zero; therefore (62) yields ln  =  − 1, which has a unique solution  = 1, or  eq =  I * .Hence, we can consider pressure  I * to be the pressure of gas-liquid coexistence at zero capillarity.It should also be noted again here that the case of   = 0 cannot be described by the model of incompressible liquid and ideal gas.
When the gas pressure is close to  I * , or  → 1 (i.e., the capillarity is sufficiently low), then (62) gives ln  ≈ −Ca, or which is the Kelvin equation for the capillary condensation which relates the pressure in gas under a curved gas-liquid interface with the pressure at a plane interface (strictly speaking, the exponential function in (63) should be expanded into the Taylor series).The case of positive Ca corresponds to a gas bubble in continuous liquid or to a porous medium wetted more by liquid than by gas.The opposite case corresponds to Ca < 0. However, (62) has no solutions if Ca < 0, as in this case the straight line  − Ca − 1 is placed under the logarithmic curve in Figure 18.This means that the liquid drops in continuum gas cannot exist in the equilibrium and are therefore unstable.

Conclusion
As shown, the development of the idea of using different EOS for gas and liquid required us to develop a mathematical tool and to introduce some new concepts such as a pseudogas, pseudoliquid, and pseudocritical point and the consistency conditions.The consistency conditions, which determine the solvability of the equilibrium equations based on two different EOS for gas and liquid, are the main elements of this theory which determine its novelty.We have obtained exact theoretical equations which determine the pseudocritical pressure (6) and the equilibrium pressure (10).For the case of the ideal gas and the exponential EOS for liquid, they have the form (17) and (19), respectively.
In Section 5, we have presented the mathematical technique of using the relationships (17) and (19) to determine the adjustable parameters of two EOS for pure components.
In particular cases, we have shown that the model of the ideal gas for all pure components is consistent and can be used to calculate gas-liquid equilibrium, while the model of incompressible liquid is nonconsistent.
We have proven that the consistency condition means graphically that two curves () for gas and liquid cross each other at two points ( I * and  II * ) and the area between these curves located between points  I * and  II * must be equivalent to the area between points  eq and  I * , as seen in Figure 3(a).This graphical construction determines the equilibrium pressure  eq and represents the analogue of the Maxwell equal area rule but turned counterclockwise.
The extension of this approach to multicomponent case leads to a more general system of equilibrium equations, which can be formulated through partial volumes (30a) or through activity coefficient (34).It was obtained from the fundamental relationships for chemical potentials through the molar volumes.This model does not require any additional adjustment, except the adjustments already done for each individual pure component.We have illustrated the developed approach using the example of ideal mixing in gas and nonideal mixing in liquid.For individual pure components, we used the virial EOS for gas or ideal gas and the exponentially compressible liquid, which satisfy the conditions of self-consistency.We have shown that a dissymmetric model is capable of ensuring explicit analytical relationships for phase concentrations.
In the last section we extended this approach to the equilibrium in presence of capillary forces, which leads to the extension of Kelvin's formula of capillary condensation to multicomponent case.
The objective of the present paper was to illustrate the feasibility of the suggested approach in principle.The comparison to some experimental data and simulations with Peng-Robinson EOS has shown that the developed approach is fundamentally able to predict the equilibrium of twocomponent and three-component mixtures.We used the data for various kinds of components: hydrocarbons (C1-C10), aromatic hydrocarbons (benzene), Oxygen, Nitrogen, Carbon Dioxide, Argon, Hydrogen Sulfide, fluorines, and chlorines.For all these data the results calculated by analytical formulae for phase concentrations are in satisfactory agreement with experimental data, even by using the simplest types of EOS.Obviously, this approach can be improved and advanced through the use of more complex EOS.
The practical advantage of using the developed theory, rather than that based on a uniform EOS for both phases, consists of the fact that this approach highly simplifies the flush calculations by leading to explicit analytical relationships for phase concentrations even for multicomponent nonideal mixtures instead of traditional nonlinear time-consuming algorithms.One efficient application of this approach would be in the field of oil engineering as currently the simulation of the multiphase compositional fluid flow is hindered by the time it takes to carry out the necessary thermodynamic calculations.

Figure 1 :
Figure 1: Illustration of the concept of the true and pseudophases in the diagram -.

Figure 3 :
Figure 3: (a) Graphical interpretation of the consistency conditions: the gray area corresponds to the integral in (10).(b) The resulting phase diagram.

Figure 7 :
Figure 7: - curves for three pure components by dissymmetric model (the black curves) and through Peng-Robinson EOS (the red dashed curves).

(
viii)    ≡   /   is the partial molar volume of component  in phase .(ix)   ≡   /  is the total molar volume of phase .(x)  I, * and  II, * are two pseudocritical pressures of pure component .(xi)   eq is the pressure of gas-liquid coexistence for pure component .(xii)  ,pure ,eq is the molar volume of pure component  in equilibrium gas state.(xiii)  ,pure ,eqis the molar volume of pure component  in equilibrium liquid state.(xiv)  is the activity coefficient.4.1.Equilibrium Equations.We will use the equilibrium equations for a multicomponent system in the following form:
(a)) and  versus the overall concentration of the light component  1 (Figure 9(b)).Diagram - was compared with the Peng-Robinson model, as the experimental data on such diagrams practically do not exist.For diagrams - 1 , we used experimental data.
, and 14 illustrate the comparison of the results obtained by the dissymmetric model with simulations based on Peng-Robinson EOS.The dew point curve is blue, and the bubble point curves are red.The tie lines are dashed.Four other Figures 15 and 16 show the comparison of the results obtained by the dissymmetric model with experimental data.
Figure 6: - curves for Decane by dissymmetric model (the black curves) and through Peng-Robinson EOS (the red dashed curves).
3.2.Example of Fitting.As the example of application of this technique, we present the results obtained for Decane for four temperatures:  = {308.9,360.4,411.9, 463.3}K.The critical pressure and temperature for it are 2.11 MPa and 617.8 K.The exact data on the equilibrium pressure and liquid volume have been calculated by using the Peng-Robinson equation of state.

Table 1 :
Parameters of dissymmetric models for pure components used to calculate binary mixtures.

Table 2 :
Parameters of fitting pure components used to calculate ternary mixtures.