Table of Contents
Journal of Thermodynamics
Volume 2016 (2016), Article ID 3806364, 21 pages
Research Article

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

1Laboratoire d’Energétique et de Mécanique Théorique et Appliquée, CNRS-UMR 7563, Université de Lorraine, 54500 Nancy, France
2Laboratory of Fluid Dynamics and Seismics, Moscow Institute of Physics and Technology (State University), Dolgoprudny 141700, Russia

Received 3 October 2015; Revised 4 January 2016; Accepted 5 January 2016

Academic Editor: Mohammad Al-Nimr

Copyright © 2016 Mikhail Panfilov and Alexandre Koldoba. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


To describe phase equilibrium of mixtures, we develop a nonclassical approach based on using 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 the mathematical theory of this new approach for pure cores and for mixtures, 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.

1. 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 [14]. 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 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 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 phi-gamma 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.

2. 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.

2.1. 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, 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 iswhere is the number of moles of fluid.

Proof of (1a). The most general conditions of equilibrium for simple fluid (without capillarity) areUsing definition (2) and the equivalence, , we obtain which yields the well known explicit relationship for the chemical potential:Substituting (5) in the equilibrium conditions (3) we obtain (1a).

In system ((1a)–(1c)) the chemical potentials and are known at reference pressures and .

Then, for any , (1a) represents the nonlinear system of one equation with respect to which is the pressure of gas-liquid 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.

2.2. 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.

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

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 and molar volumes and (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: , , and .

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 , while the true liquid corresponds to .

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 , while the pseudogas is the prolongation of the true gas into the domain of volumes lower than .

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 , 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: . 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)) becomeswhich is the algebraic equation with respect to the equilibrium pressure .

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 .

2.4. 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) leads to orwhich is the system of two nonlinear equations for two equilibrium molar volumes: and .

For a cubic EOS, as the Van-der-Waals equation, , 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 .

Figure 2: Graphical technique of determining the equilibrium line in the case of using the uniform (cubic) EOS for gas and liquid.

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 . 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.

2.5. 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 . 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:At fixed ,

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

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 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 and .

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 two-phase equilibrium is impossible for an incompressible liquid. Indeed, any monotonic function cannot cross twice the vertical line .

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 and (Figure 3);(c)the curve for pseudogas is situated on the left from that of pseudoliquid between and , according to condition (13b).(ii)The equilibrium pressure is determined graphically in such a way that the surface of the area between two curves between points and (Figure 3) is identical to the area of the surface between points and .

In the classical case considered in the previous section the equilibrium equations (12) are solved with respect to two different volumes and . In contrast to it, the solution of (10) is searched directly with respect to the unique value of pressure .

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.

2.6. 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, is an arbitrary constant pressure value (different from in (1a)–(1c) and (5)), while 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: , 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.

Figure 4: Solvability of (17).

Two solutions and transform to a single value if the exponential function touches the straight line , which means the identical derivatives of and : . 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 crosses the line twice at two points and . When , the first pressure tends to and the second pressure 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 is determined from (17), the equilibrium pressure, , is found as the solution to the following equation resulting from (14):

The value of is the nontrivial solution of this equation, different from .

So, for given , one can get and . Since we know the gas and liquid EOS we can get and and, finally, the envelope of the phase diagram (this being the line joining all the dew points and all the bubble points).

3. Determination of the Coefficients

3.1. 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 (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 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 formwhere the new parameter should be adjusted by using the data on the equilibrium gas volume, .

The virial EOS gives more accurate results for equilibrium gas volume than the ideal gas EOS.

Let us assume that the equilibrium pressure , the equilibrium molar liquid, and gas volumes and are given as experimental data. For pure components, they can be found in the literature or calculated by Peng-Robinson EOS.

The objective is to determine the values of , , , and . Parameter is selected in a free way (in the examples below it was chosen as ). Then the system of (17), (19), and (20) and the second equation in (15) are sufficient to determine , , , and through , , and . This system is

Introducing new dimensionless variables, we can transform system (21) to one nonlinear equation which determines the function , where and are the given parameters,and an explicit relationship for function :

If the gas phase is ideal, then ; then we obtain for and instead of (23)

Proof. (1) Eliminating and from the last two equations (21), we obtain from the first two equations (2) Let us eliminate from the second relationship, which gives (24). Substituting it into the first equation, we deduce (23).

Let us note that the functions and are universal and do not depend on individual properties of chemical components.

The solution of (25) obtained numerically is shown in Figure 5.

Figure 5: Functions and defined by (25).

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: ; then .

Note that the boundary values of these functions which correspond to the formal critical point are , , and .

Using and , we can determine parameters , , , and :

3.2. Example of Fitting

As the example of application of this technique, we present the results obtained for Decane for four temperatures:  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.

The values of parameter are calculated from its definition; parameter is calculated as the solution to (25). The reference pressure was selected as . The obtained values of the coefficient of compressibility are (MPa−1) as follows: and the reference molar liquid volume , calculated from (27), is (m3/mol):

The results are presented in Figure 6.

Figure 6: curves for Decane by dissymmetric model (the black curves) and through Peng-Robinson EOS (the red dashed curves).

In Figure 7, we present the diagrams for several pure components calculated by the dissymmetric model and Peng-Robinson EOS, for  K.

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

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.

4. 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 “” 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) is the chemical potential of component in the mixture if the mixture behaves as an ideal gas.(v) 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 .(viii) is the partial molar volume of component in phase .(ix) is the total molar volume of phase .(x) and are two pseudocritical pressures of pure component .(xi) is the pressure of gas-liquid coexistence for pure component .(xii) is the molar volume of pure component in equilibrium gas state.(xiii) is 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:

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 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 variables: phase concentrations and (due to (30d)), pressure , and temperature . Such a system is unclosed and contains independent variables. We select them as , , and 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, , 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 material balance equations:where is the gas molar fraction.Then the system of equations (30a) and (31) contains unknown variables: 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: (), and the disappearing gas phase has some limit values of phase concentrations.Then the system of equations (30a) contains unknown variables: limit concentrations and pressure , which is called the bubble-pressure, . Its is clear that 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 () and so on, similar to bubble point.

4.2. 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: , where 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 : , 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, ).

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 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.

4.3. Multicomponent Phase Diagrams

The expected behaviour of the phase diagrams of mixtures is easily revealed from the equilibrium equations (30a). Let us consider a two-component mixture (component 1 is light), which satisfies the condition of the ideal mixing both in liquid and in gas. Then , , and system (30a) becomes

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 , then it will continuously vary between and 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.

Figure 8: Example of consistency of the same EOS for gas for individual components: graphical determination of the equilibrium pressures (a) and (b).

Figures 8(a) and 8(b) illustrate the graphical determination of the equilibrium pressures and for two pure components. Pressure is determined in such a way that the yellow surface would be equal to the yellow surface . Pressure 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.

4.4. 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 , or for all . From (34), we obtain the following equilibrium equations for :

Using (35), we obtain the expression of the activity coefficients through the phase concentrations:

4.5. 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 ():where , , , , and 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 .

The equilibrium equations (39) become ()

In dimensionless form these equations take the following form:where

4.6. 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, and , as the functions of pressure .

This system has an explicit analytical solution for phase concentrations and in two cases, the ideal mixing in both phases and 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 ):

For nonideal mixing in liquid, let us assume the hypothesis of moderate nonideality: the interactions between different molecules and identical molecules are nearly equal: , which leads to . Then -values in (46) are calculated aswhere Coefficients and are calculated by using formulae (47). This solution is developed in Appendix B.

4.7. Analytical Solutions for Phase Concentrations in Three-Component Case

In a three-component system, three equations of phase equilibrium contain four phase concentrations , , , and , which means one concentration is independent. We will select as this independent parameter. Then the equilibrium system determines three other concentrations as the functions of the equilibrium pressure, temperature, and . 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 ():

For nonideal mixing in liquid, where is calculated through (40) in which one should use instead of , for all ; , , and 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 considered as independent in the equilibrium equations:where is the gas molar fraction. We obtain two equations with two unknown variables: and .

The envelopes of the two-phase domain are calculated by using the same equations in which for the bubble curve (transition to single-phase liquid) and for the dew curve (transition to single-phase gas), which results from (54) if and , respectively.

5. 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 ,” where is the total concentration of the lightest component. At high , 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]:(i)2-Methylpentane (C6H14) and Heptane (C7H16).(ii)2-Methylpentane and Octane (C8H18).(iii)Heptane and Octane.(iv)Butane (C4H10) and Decane (C10H22).(v)Pentane (C5H12) and Benzene (C6H6).

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.

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

For mixture, the phase diagrams (the envelops of the two-phase zone) were constructed by using the model presented in this paper in terms of versus the molar volume (see Figure 9(a)) and versus the overall concentration of the light component (Figure 9(b)). Diagram was compared with the Peng-Robinson model, as the experimental data on such diagrams practically do not exist. For diagrams , we used experimental data.

Figure 9: Binary phase diagrams for 2-Methylpentane (C6H14)-Octane (C8H18) at  K.

The comparison of calculations by using the dissymmetric model with experimental data for Butane (C4H10)-Decane (C10H22) mixture is presented in Figure 10.

Figure 10: Binary phase diagrams for Butane (C4H10)-Decane (C10H22) mixture at  K.

The case of nonideal mixing is illustrated in Figure 11 for the mixture Pentane-Benzene, for which there exist experimental data.

Figure 11: Binary phase diagrams for Pentane (C5H12)-Benzene (C6H6) mixture at  K.

Two other cases of comparison with experimental data are shown in Figure 17 for Heptane and Octane (a) and for 2-Methylpentane and Heptane (b).

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.

5.2. 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 (), 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.

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

The first three Figures 12, 13, 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.

Figure 12: Ternary phase diagram for Propane-Pentane-Decane mixture, for bar and  K: (a) the dissymmetric model and (b) Peng-Robinson EOS.
Figure 13: Ternary phase diagram for 2-Methylpentane-Heptane-Octane mixture, for bar and  K: (a) the dissymmetric model and (b) Peng-Robinson EOS.
Figure 14: Ternary phase diagram for Hydrogen Sulfide (H2S)-Propane (C3H8)-Pentane (C5H12), for bar and  K: (a) the dissymmetric model and (b) the Peng-Robinson EOS.

Four other Figures 15 and 16 show the comparison of the results obtained by the dissymmetric model with experimental data.

Figure 15: Ternary phase diagrams calculated by the dissymmetric model (lines) and experimental data (crosslets): (a) mixture of Nitrogen (N2)-Argon (Ar)-Methane (CH4) at  K and bar; experimental data from [16]; (b) mixture of Chlorodifluoromethane (R22)-Nitrogen (N2)-Carbon Dioxide (CO2) at  K and bar; experimental data from [17].
Figure 16: Ternary phase diagrams calculated by the dissymmetric model (lines) and experimental data (crosslets): (a) mixture of Dichloromethane (R30)-Tetrachloromethane (R10)-Trichloromethane (R20) at  K and bar; experimental data from [18]; (b) mixture of Chlorodifluoromethane (R142b)-1,1,1-Trichloroethane (R140a)-1,1-Dichloro-1-fluoroethane (R141b) at  K and bar; experimental data from [19].
Figure 17: Binary phase diagrams for Heptane (C7H16)-Octane (C8H18) (a) and 2-Methylpentane (C6H14)-Octane (b) with capillary pressure and without it at  K.

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).

6. Application to Capillary Condensation

6.1. 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.

The general equilibrium equations ((30a)–(30d)) for multicomponent mixtures becomewhere is the gas pressure.

For virial EOS for gas and exponential EOS for liquid, we obtain, similar to (43),where is the capillary number.

6.2. Calculation of the Phase Composition at Nonzero Capillary Pressure

For ternary mixture, the following can be obtained for -values (), similar to (52):

Then the phase concentrations are calculated by using the same formulae (51).

Parameter is zero for ideal mixture. For nonideal solution, it is calculated through (40), in which one should use instead of , for all ; , , and are calculated by using formulae (51) but with -values corresponding to ideal mixing, that is, by using (60), in which .

The examples of calculation of the phase diagrams at nonzero capillary pressure are given in Figure 17.

The values of we have used are proper to natural porous media: for pure components the curvature radius of gas-liquid interface is assumed to be of the order of pore radius (~200), and the surface tension of each component was calculated using the Brock and Bird Model. The capillary pressure of the mixture was estimated as the mean value ().

6.3. Kelvin’s Equation of Capillary Condensation for One-Component Case

For one-component fluid, the equilibrium equations (58) determine a unique equilibrium pressure in gas phase :

At given EOS for pseudogas, , and for pseudoliquid , (61) represents a nonlinear equation with respect to the equilibrium pressure .

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 intersects the hyperbolic curve of the ideal gas only at one point : . Then we obtain the following from (61) by integrating: . This yieldswhere ; the capillary number .

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 on the horizontal axis corresponds to ).

Figure 18: Graphical interpretation of (62). The physical nontrivial solution corresponds to point .

Point has no physical meaning, as the corresponding value ; then . This means that the molar volume of the coexisting liquid is higher than that of gas (), which is physically impossible. So the true solution corresponds to point . We obtain the characteristic pressure which is lower than , 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 , which has a unique solution , or . Hence, we can consider pressure to be the pressure of gas-liquid coexistence at zero capillarity. It should also be noted again here that the case of cannot be described by the model of incompressible liquid and ideal gas.

When the gas pressure is close to , or (i.e., the capillarity is sufficiently low), then (62) gives , orwhich 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 . However, (62) has no solutions if , as in this case the straight line 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.

7. 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 ( and ) and the area between these curves located between points and must be equivalent to the area between points and , as seen in Figure 3(a). This graphical construction determines the equilibrium pressure 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 two-component 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.


A. Proof of (30a)

(1.) The following relationship for the chemical potential of component in an arbitrary phase through the chemical potential of the ideal gas is known in classical thermodynamics [14]:where is an arbitrary reference value of pressure which may be different for various components. Then

(2.) For a pure component we obtain

(3.) Let us subtract these two relationships from each other:

(4.) Using (5), we obtain

(5.) Substituting the last relationship in the general equilibrium equation, , we obtain

(6.) Selecting the reference pressures as the second pseudocritical pressures of individual components, , we obtain (30a), because at the point where the difference between the pressure, temperature, and molar volumes of pseudogas and pseudoliquid disappears.

B. Proof of (47) and (48)

Let us assume that is a small parameter. Then one represents the solution in the following form:where and are the values which correspond to the ideal mixing and and are the small deviations of order .

All the values are of order and disappear in zero approximation. Substituting these expressions in (43) and applying the method of perturbation, one obtains the following for the zero approximation:And, for the terms of order ,where is defined through (40) in which ().

Relationships (B.2) give (47) for , while (B.1) and (B.3) yield , which leads to (48).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


The research was partially supported by the Groupement de Recherche MOMAS: Modélisation Mathématique et Simulations Numériques liées aux études d’entreposage souterrain de déchets radioactifs (CNRS, ANDRA, EDF, CEA, BRGM).


  1. A. I. Brusilovsky, Phase Transition Under Exploitation of Oil and Gas Reservoirs, Graal, Moscow, Russia, 2002 (Russian).
  2. A. Firoozabadi, Thermodynamics of Hydrocarbon Reservoirs, McGraw-Hill, New York, NY, USA, 1999.
  3. M. Hillert, Phase Equilibria, Phase Diagrams and Phase Transformations, Cambridge University Press, Cambridge, UK, 2007.
  4. S. M. Walas, Phase Equilibria in Chemical Engineering, Elsevier, Butterworth, Malaysia, 1985.
  5. M. L. Michelsen, “The isothermal flash problem. Part I. Stability,” Fluid Phase Equilibria, vol. 9, no. 1, pp. 1–19, 1982. View at Publisher · View at Google Scholar · View at Scopus
  6. M. L. Michelsen, “The isothermal flash problem. Part II: phase split calculation,” Fluid Phase Equilibria, vol. 9, no. 1, pp. 21–40, 1982. View at Publisher · View at Google Scholar · View at Scopus
  7. M. L. Michelsen, “Simplified flash calculations for cubic equations of state,” Industrial & Engineering Chemistry Process Design and Development, vol. 25, no. 1, pp. 184–188, 1986. View at Publisher · View at Google Scholar · View at Scopus
  8. B. H. Jensen and A. Fredenslund, “A simplified flash procedure for multicomponent mixtures containing hydrocarbons and one non-hydrocarbon using two-parameter cubic equations of state,” Industrial and Engineering Chemistry Research, vol. 26, no. 10, pp. 2129–2134, 1987. View at Publisher · View at Google Scholar · View at Scopus
  9. Y. Li and R. T. Johns, “Rapid flash calculations for compositional simulation,” SPE Reservoir Evaluation & Engineering, vol. 9, no. 5, pp. 521–529, 2006. View at Publisher · View at Google Scholar · View at Scopus
  10. S. Oladyshkin and M. Panfilov, “Limit thermodynamic model for compositional gas-liquid systems moving in a porous medium,” Transport in Porous Media, vol. 70, no. 2, pp. 147–165, 2007. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  11. S. Oladyshkin, S. Skachkov, I. Panfilova, and M. Panfilov, “Upscaling fractured media and streamline HT-splitting in compositional reservoir simulation,” Oil & Gas Science and Technology, vol. 62, no. 2, pp. 137–146, 2007. View at Publisher · View at Google Scholar
  12. A. P. Plres, P. G. Bedrikovetsky, and A. A. Shapiro, “Analytical modeling for two-phase EOR processes: splitting between hydrodynamics and thermodynamics,” in Proceedings of the SPE Annual Technical Conference and Exhibition, SPE Paper 89919 MS, pp. 573–583, Houston, Tex, USA, September 2004. View at Scopus
  13. A. V. Koldoba and E. V. Koldoba, “Model equation of state and Gibbs energy for numerical modeling of the problems of multicomponent flow with phase transition in porous media,” Geochemistry, vol. 4, pp. 573–576, 2004 (Russian). View at Google Scholar · View at Scopus
  14. A. A. Shapiro and E. H. Stenby, “Kelvin equation for a non-ideal multicomponent mixture,” Fluid Phase Equilibria, vol. 134, no. 1-2, pp. 87–101, 1997. View at Publisher · View at Google Scholar · View at Scopus
  15. H. Orbey and S. I. Sandler, Modeling Vapor-Liquid Equilibria: Cubic Equations of State and Their Mixing Rules, Cambridge University Press, Cambridge, UK, 1998.
  16. R. C. Miller, A. J. Kidnay, and M. J. Hiza, “Liquid-vapor equilibria at 112.00 K for systems containing nitrogen, argon, and methane,” AIChE Journal, vol. 19, no. 1, pp. 145–151, 1973. View at Publisher · View at Google Scholar
  17. X. Wang, Y. Wang, J. Shi, and B. C.-Y. Lu, “Isothermal vapor-liquid equilibria at elevated pressures for the systems containing nitrogen, carbon dioxide, and chlorodifluoromethane,” Journal of Chemical and Engineering Data, vol. 36, no. 4, pp. 436–439, 1991. View at Publisher · View at Google Scholar · View at Scopus
  18. G. Rulewicz, H. Schuberth, and E. Leibnitz, “Untersuchungen über das Dampf-Flüssigkeits-Phasengleichgewicht des Systems Methylenchlorid/Chloroform/Tetrachlorkohlenstoff bei 45°C unter Verwendung einer gaschromatographischen Analysenmethode,” Journal für Praktische Chemie, vol. 37, no. 3-4, pp. 122–136, 1968. View at Publisher · View at Google Scholar
  19. Y. W. Kang and Y. Y. Lee, “Vapor-liquid equilibria for the systems composed of 1-chloro-1,1-difluoroethane, 1,1-dichloro-1-fluoroethane, and 1,1,1-trichloroethane at 50.1°C,” Journal of Chemical and Engineering Data, vol. 41, no. 2, pp. 303–305, 1996. View at Publisher · View at Google Scholar · View at Scopus
  20. H. Renon and J. M. Prausnitz, “Local compositions in thermodynamic excess functions for liquid mixtures,” AIChE Journal, vol. 14, no. 1, pp. 135–144, 1968. View at Publisher · View at Google Scholar
  21. J. Gmehling, U. Onken, and W. Arlt, Chemistry Data Series: Vapor-Liquid Equilibrium Data Collection, Dechema, 1977.