A number of industrially important processes, such as distillation, absorption, and extraction, bring into contact two phases that are not at equilibrium. The rate at which a species is transferred from one phase to the other depends on the departure of the system from equilibrium. The quantitative treatment of these rate processes requires knowledge of the equilibrium states of the system. Apart from these cases, vapor-liquid equilibrium (VLE) relationships in multicomponent systems are needed for the solution of many other classes of engineering problems, such as the computation of evaporation rates in spray combustion applications.

For the multicomponent particle, the vaporization rate of component is computed from Equation 12–162 or Equation 12–163 and requires knowledge of the component concentration or the mass fraction at the droplet surface.

The concentration or mass fraction of species at the particle surface is computed from the mole fraction , which can be obtained from the general expression for two phase equilibrium, equating the fugacity of the liquid and vapor mixture components [525]:

(12–174)

where is the mole fraction, is the fugacity coefficient for the species in the mixture, and is the absolute pressure. The superscripts and are the vapor and the liquid phase variables, respectively. The fugacity coefficients account for the nonideality in the gas and liquid mixture. The fugacity of the liquid phase can be calculated from the pure component’s saturation pressure [610]:

(12–175)

Here, is the fugacity coefficient for pure at the saturation pressure; is the activity coefficient for species in the mixture, and accounts for the nonideality in the liquid phase; is the particle surface temperature. We assume perfect thermal conductivity inside the particle, so the particle temperature is used instead; is the universal gas constant; is the molar volume of the liquid. The exponential term is the Poynting correction factor and accounts for the compressibility effects within the liquid. Except at high pressures, the Poynting factor is usually negligible. Under low pressure conditions where the gas phase may be assumed to be ideal, and . Furthermore, if the liquid is also assumed to be ideal, , then Equation 12–174 reduces to Raoult’s law,

(12–176)

Raoult’s law is the default vapor-liquid equilibrium expression used in the Ansys Fluent multicomponent droplet model. However, there is a UDF hook available for user-defined vapor-liquid equilibrium models.

While Raoult’s law represents the simplest form of the VLE equation, keep in mind that it is of limited use, as the assumptions made for its derivation are usually unrealistic. The most critical assumption is that the liquid phase is an ideal solution. This is not likely to be valid, unless the system is made up of species of similar molecular sizes and chemical nature, such as in the case of benzene and toluene, or n-heptane and n-hexane.

For higher pressures, especially near or above the critical point of the components, real gas effects must be considered. Most models describing the fugacity coefficients use a cubic equation of state with the general form:

(12–177)

where is the molar volume. For in-cylinder applications, the Peng-Robinson equation of state is often used [510], where , , and :

(12–178)

This equation defines the compressibility

(12–179)

The implementation of the Peng-Robinson equation of state in Ansys Fluent uses this expression for both phases, the particle liquid and the vapor phase. The parameters and are determined by the composition using a simple mixing law:

(12–180)

where is the number of components in the mixture. The pure component parameters can be obtained using the relationship with the Peng-Robinson constants:

(12–181)

where is the critical temperature, is the critical pressure and is the acentric factor of the component .

The fugacities of the components depend on the compressibility of the liquid and vapor phase:

(12–182)

where is the residual Helmholtz energy, which is a function of the compressibility:

(12–183)

In summary, the vapor mole fraction , the pressure , and the compressibilities of the vapor () and the liquid () phase at the surface of the particle are determined from the liquid particle mole fraction of the components and the particle temperature . The surface vapor concentrations are calculated using the following equation:

(12–184)