A new numerical procedure to determine the VLE curve
Introduction
The knowledge of phase transitions of pure or mixed fluids is essential in many practical applications. Engineers and applied physicists and chemists commonly use an equation of state (EOS) to evaluate the vapour-liquid equilibrium (VLE) properties. However, this calculation is no easy task numerically, and the shortcut techniques that have been developed are not applicable to every EOS (Wisniak et al., 1998).
As is well known, to determine the VLE curve of a pure fluid, it is necessary to solve the equationswhere the pressure, P, and the chemical potential, μ, are analytical functions of temperature, T, and density, ρ.
One way to obtain the VLE properties from a given EOS is through the equality of the areas for each temperature. This is equivalent to the equality of the chemical potentials corresponding to the gas and liquid phases at a given temperature. In practice however, the accurate determination of the position of the straight section of the isotherm for which the pressures and areas must be equal is no trivial problem, because the isotherm has a steep slope for densities (molar volumes) corresponding to the liquid phase. A small deviation in the equality of the areas involves small deviations in the densities ρL of the liquid phase and ρV of the gas phase. Because of the steepness of the slope of the isotherm in the liquid region, the deviation in the determination of PL is extraordinarily sensitive to the accuracy of the determination of the liquid density, ρL, and hence, so is the accuracy in the determination of the mechanical equilibrium between the two phases.
This problem may be avoided by solving simultaneously the system of equations (1), to which end one would need analytical expressions not only for the EOS but also for the chemical potential. Whereas there have been numerous EOSs in the literature, analytical expressions for the chemical potential are more difficult to find (Mulero et al., 1999).
Various EOSs have been constructed on the basis of intensive computer simulation data for Lennard–Jones (LJ) fluids (Johnson et al., 1993, Kolafa and Nezbeda, 1994, Cuadros et al., 1996). But the complexity of the corresponding and expressions makes it very difficult to solve system (1) analytically. Instead, it is necessary to use numerical techniques. Thus, to determine each pair of points , one must solve system (1) for a given temperature, with the inherent associated problem of choosing the initial value, especially in the neighbourhood of the critical point.
We present here an easier and faster method of solving system (1). This method yields practically the entire VLE curve, without having to solve the system point by point. It is based on numerically solving a system of nonlinear differential equations using the program mathematica (Wolfram, 1991). This program offers to the user an interactive tool for numerical calculations and a versatile programming language for fast and accurate solutions to the most diverse problems. For these reasons, it is becoming ever more commonly used by scientists, engineers, and higher education students.
The main goal of the present paper is to describe a new procedure to calculate the coexistence curve between the liquid and gas phases of a LJ system. The method could be extrapolated to any system for which one knows the EOS and the expression
Section snippets
Theoretical basis of the method
For any T below the critical temperature TC, the functions and satisfy Eq. (1). Since P and μ are known analytically, one then differentiates Eq. (1) with respect to T:and hencewhere and denote the derivatives of ρV and ρL with respect to T.
From Eq. (2), one obtains the system of two first-order nonlinear differential equations
Algorithm
- 1.
Initiate T0.
- 2.
Define the functions and
- 3.
Search for a numerical solution, using Newton's method, to the simultaneous equations and
- 4.
Assign the solution's values to ρV0 and ρL0.
- 5.
Find a numerical solution to the ordinary differential equations and , for the functions ρV and ρL with the independent variable T in the range 0.7 to 1.35.
As we noted above, to solve the system of two nonlinear differential equations
Applications to calculating the VLE curve
For our calculations, we chose three EOS models: Johnson et al., 1993, Kolafa and Nezbeda, 1994, and Cuadros et al., 1996. All the model EOSs will be expressed in reduced LJ units.
Conclusions
- 1.
We have described a new approach to calculating the VLE curve, system (1), which is mathematically equivalent to the classical requirement of solving the system that results from differentiating with respect to T.
- 2.
If the EOS is sufficiently regular, with the imposition of determined initial conditions , the solution of the system (2) exists and is unique.
- 3.
The complexity of the EOS expressions makes it practically impossible to obtain analytical formulae for ρV and ρL We therefore
Acknowledgements
The authors express their gratitude to the Consejerı́a de Educación y Juventud of the Junta de Extremadura and the Fondo Social Europeo for the financial aid granted through the Project IPR98B004. Also, M.I. Parra would like to express her gratitude to the Socrates–Erasmus program for its support through a grant, as well as for facilitating her academic visit to the Institute of Mathematics of the Politecnika Zielonogórska.
References (10)
- et al.
Fluid Phase Equilibria
(1994) - et al.
Chem. Eng. Sci.
(1998) - et al.
J. Chem. Phys.
(1969) - et al.
Phase Transitions
(1996) Ordinary Differential Equations
(1964)
Cited by (7)
Simple approach to approximate predictions of the vapor-liquid equilibrium curve near the critical point and its application to Lennard-Jones fluids
2012, Physics Letters, Section A: General, Atomic and Solid State PhysicsCitation Excerpt :The result of this procedure was to remove a small amount of data ‘scatter’ of graphs illustrating the percentage deviation of the values being compared (see Figs. 2, 4). Argon and methane are the subject of numerous scientific studies [16–20], providing coverage of both theoretical and experimental studies. It has been payed much attention to the thermodynamic properties, due to widespread use of this element both in industry and science.
Many-body and quantum effects in the surface tension and surface energy of liquid neon and argon using the Fowler's approximation
2012, Chemical PhysicsCitation Excerpt :Although the RDF expression of Cuadros et al. [47] is a complex analytical model, but it is (such as the RDF equation of Morsali et al. [40]) in good agreement with the real RDF obtained in computer simulations [9,47]. Therefore, the difference between the calculated values in this work and those of Mulero et al. [9] can be due to the difference between the experimental densities [46] and the density equations of Okrasinski et al. [48,49]. In the other words, the Okrasinski et al. expressions produce the density values somewhat bigger than those of Ref. [46].
Modeling and simulation of supercritical extraction columns using computer algebra
2005, Journal of Supercritical FluidsThe slowly varying density approximation for the surface tension and surface energy of Lennard-Jones fluids revisited
2003, Journal of Physics Condensed MatterFowler's approximation for the surface tension and surface energy of Lennard-Jones fluids revisited
2003, Journal of Physics Condensed Matter