Generalized corresponding states model for bulk and interfacial properties in pure fluids and fluid mixtures

Predrying transition on a hydrophobic surface: statics and dynamics

For one-component fluids, we examine the predrying phase transition between a thin and thick low-density layer in liquid on a wall repelling the fluid. This is the case of a hydrophobic wall for water. A predrying line starts from the coexistence curve and ends at a surface critical point in the phase diagram. We calculate this line numerically using the van der Waals model and analytically using the free-energy expansion up to the quartic order. We also examine the predrying dynamics of a layer created on a hydrophobic spot on a heterogeneous wall. It is from a thin to thick layer during decompression and from a thick to thin layer during compression. Upon the transition, a liquid region above the film is cooled for decompression and heated for compression due to latent heat convection, and a small pressure pulse is emitted from the film into the liquid.

Solvation effects in phase transitions in soft matter

Phase transitions in polar binary mixtures can be drastically altered by even a small amount of salt. This is because the preferential solvation strongly depends on the ambient composition. Together with a summary of our research on this problem, we present some detailed results on the role of antagonistic salt composed of hydrophilic and hydrophobic ions. These ions tend to segregate at liquid-liquid interfaces and selectively couple to water-rich and oil-rich composition fluctuations, leading to mesophase formation. In our two-dimensional simulation, the coarsening of the domain structures can be stopped or slowed down, depending on the interaction parameter (or the temperature) and the salt density. We realize stripe patterns at the critical composition and droplet patterns at off-critical compositions. In the latter case, charged droplets emerge with considerable size dispersity in a percolated region. We also give the structure factors among the ions, accounting for the Coulomb interaction and the solvation interaction mediated by the composition fluctuations.

Precipitation in aqueous mixtures with addition of a strongly hydrophilic or hydrophobic solute

We examine phase separation in aqueous mixtures due to preferential solvation with a low-density solute (hydrophilic ions or hydrophobic particles). For hydrophilic ions, preferential solvation can stabilize water domains enriched with ions. This precipitation occurs above a critical solute density n(p) in wide ranges of the temperature and the average composition, where the mixture solvent would be in a one-phase state without solute. The volume fraction of precipitated domains tends to zero as the average solute density n is decreased to np or as the interaction parameter χ is decreased to a critical value χ(p). If we start with one-phase states with n>n(p) or χ>χ(p), precipitation proceeds via homogeneous nucleation or via heterogeneous nucleation, for example, around suspended colloids. In the latter case, colloid particles are wrapped by thick wetting layers. We also predict a first-order prewetting transition for n or χ slightly below np or χ(p) for neutral colloids.

A new crossover sine model based on trigonometric model and its application to the crossover lattice equation of state

In this study, a new crossover sine model (CSM) n was developed from a trigonometric model [M. E. Fisher, S. Zinn, and P. J. Upton, Phys. Rev. B 59, 14533 (1999)]. The trigonometric model is a parametric formulation model that is used to represent the thermodynamic variables near a critical point. Although there are other crossover models based on this trigonometric model, such as the CSM and the analytical sine model, which is an analytic formulation of the CSM, the new sine model (NSM) employs a different approach from these two models in terms of the connections between the parametric variables of the trigonometric model and thermodynamic variables. In order to test the performance of the NSM, the crossover lattice equation of state [M. S. Shin, Y. Lee, and H. Kim, J. Chem. Thermodyn. 40, 174 (2008)] was applied using the NSM for correlations of various pure fluids and fluid mixtures. The results showed that over a wide range of states, the crossover lattice fluid (xLF)/NSM yields the saturated properties of pure fluids and the phase behavior of binary mixtures more accurately than the original lattice equation of state. Moreover, a comparison with the crossover lattice equation of state using the CSM (xLF/CSM) showed that the new model presents good correlation results that are comparable to the xLF/CSM.

A new model approach for the near-critical point region: 1. Construction of the generalized van der Waals equation of state

To date, it has been considered that all classical equations of state (EOS) have failed to describe the properties of fluids near the critical region, where the density fluctuations have a significant influence on fluid properties. In this paper, we suggest a newly constructed equation for fluid states, the generalized van der Waals (GvdW) EOS with the highly simplified Dieterici's form P = [RT/(V - b)] - a(b/V)c by a new model potential construction describing intermolecular interactions. On the basis of the model potential construction, it is shown that the a, b, and c parameters have physical interpretations as an internal pressure, a void volume, and a dimensionless value that represents an inharmonic intermolecular cell potential, respectively. As an illustration of our model approach, we initially apply it to near the critical point (cp) region, where all classical EOS descriptions have been incorporated with experimental thermodynamic data, and we obtain a table of three parameters for 12 pure normal fluids, which precisely describes thermodynamic critical values. On the basis of the basic relations between pressure and volume at the critical point, we express the corresponding EOS in terms of the c parameter, and by this means, we also obtain a theoretical vapor-liquid equilibrium (VLE) line, which closely coincides with the experimental data for several pure normal fluids near the critical region. As a result, we show that thermodynamic properties near the critical region can be described analytically by only three parameters. In addition, to validate our EOS for the temperature-differential derivatives, we show that the calculated isochoric heat capacity (Cv) of saturated argon closely coincides with the experimental data. Moreover, the possibility of a precise description with respect to the entire fluid region is also argued, in terms of the physical cases from the triple point to the ideal gas region.

Master crossover functions for one-component fluids

By introducing three well-defined dimensionless numbers, we establish the link between the scale dilatation method able to estimate master (i.e., unique) singular behaviors of the one-component fluid subclass and the universal crossover functions recently estimated [Garrabos and Bervillier, Phys. Rev. E 74, 021113 (2006)] from the bounded results of the massive renormalization scheme applied to the Phi(d)(4)(n) model of scalar order parameter (n=1) and three dimensions (d=3), representative of the Ising-like universality class. The master (i.e., rescaled) crossover functions are then able to fit the singular behaviors of any one-component fluid without adjustable parameter, using only one critical energy scale factor, one critical length scale factor, and two dimensionless asymptotic scale factors, which characterize the fluid critical interaction cell at its liquid-gas critical point. An additional adjustable parameter accounts for quantum effects in light fluids at the critical temperature. The effective extension of the thermal field range along the critical isochore where the master crossover functions seems to be valid corresponds to a correlation length greater than three times the effective range of the microscopic short-range molecular interaction.

Master crossover behavior of parachor correlations for one-component fluids

The master asymptotic behavior of the usual parachor correlations, expressing surface tension sigma as a power law of the density difference rho(L)-rho(V) between coexisting liquid and vapor, is analyzed for a series of pure compounds close to their liquid-vapor critical point, using only four critical parameters (beta(c))-1 , alpha(c) , Z(c) , and Y(c) , for each fluid. This is accomplished by the scale dilatation method of the fluid variables where, in addition to the energy unit (beta(c))-1 and the length unit alpha(c) , the dimensionless numbers Z(c) and Y(c) are the characteristic scale factors of the ordering field along the critical isotherm and of the temperature field along the critical isochore, respectively. The scale dilatation method is then formally analogous to the basic system-dependent formulation of the renormalization theory. Accounting for the hyperscaling law delta-1/delta+1=eta-2/2d , we show that the Ising-like asymptotic value pi(a) of the parachor exponent is unequivocally linked to the critical exponents eta or delta by pi(a)/d-1=2/d-(2-eta)=delta+1/d (here d=3 is the space dimension). Such mixed hyperscaling laws combine either the exponent eta or the exponent delta , which characterizes bulk critical properties of d dimension along the critical isotherm or exactly at the critical point, with the parachor exponent pi(a) which characterizes interfacial properties of d-1 dimension in the nonhomogeneous domain. Then we show that the asymptotic (symmetric) power law [abstract; see text] is the two-dimensional critical equation of state of the liquid-gas interface between the two-phase system at constant total (critical) density rho(c) . This power law complements the asymptotic (antisymmetric) form [abstract; see text] of the three-dimensional critical equation of state for a fluid of density rho not equal to rho_(c) and pressure p not equal to p_(c) , maintained at constant (critical) temperature T=T_(c)} [mu_(rho)(mu_(rho,c)) is the specific (critical) chemical potential; p_(c) is the critical pressure; and T_(c) is the critical temperature]. We demonstrate the existence of the related universal amplitude combination [abstract; see text] = universal constant, constructed with the amplitudes D_(rho)(sigma) and D_(rho)(c) , separating then the respective contributions of each scale factor Y_(c) and Z_(c) , characteristic of each thermodynamic path, i.e., the critical isochore and the critical isotherm (or the critical point), respectively. The main consequences of these theoretical estimations are discussed in light of engineering applications and process simulations where parachor correlations constitute one of the most practical methods for estimating surface tension from density and capillary rise measurements.

Dynamic van der Waals theory

We present a dynamic van der Waals theory starting with entropy and energy functional with gradient contributions. The resultant hydrodynamic equations contain the stress arising from the density gradient. It provides a general scheme of two-phase hydrodynamics involving the gas-liquid transition in nonuniform temperature. Some complex hydrodynamic processes with evaporation and condensation are examined numerically. They are (i) adiabatically induced spinodal decomposition, (ii) piston effect with a bubble in liquid, (iii) temperature and velocity profiles around a droplet in heat flow, (iv) efficient latent heat transport at small liquid densities (the mechanism of heat pipes), (v) boiling in gravity with continuous bubble formation and rising, and (vi) spreading and evaporation of liquid on a heated boundary wall.

Study on Vapor−Liquid Equilibria and Surface Tensions for Nonpolar Fluids by Renormalization Group Theory and Density Gradient Theory

An equation of state (EOS) applicable for both the uniform and nonuniform fluids is established by using the density-gradient theory (DGT). In the bulk phases, the EOS reduces to statistical associating fluid theory (SAFT). By combining the EOS with the renormalization group theory (RGT), the vapor-liquid-phase equilibria and surface tensions for 10 nonpolar chainlike fluids are investigated from low temperature up to the critical point. The obtained results agree well with the experimental data.

Mean crossover functions for uniaxial three-dimensional Ising-like systems

We give simple expressions for the mean of the max and min bounds of the critical-to-classical crossover functions, previously calculated [Bagnuls and Bervillier, Phys. Rev. E 65, 066132 (2002)] within the massive renormalization scheme of the Phi(d)4(n) model in three dimensions (d = 3) and scalar order parameter (n = 1) of the Ising-like universality class. The mean functions are determined relying on the properties of the theoretical functions in the two limiting three-dimensional (3D) Ising-like and mean-field-like descriptions close to the Wilson-Fisher fixed point and to the Gaussian fixed point, respectively. Such descriptions correspond to the preasymptotic domains near each fixed point where a Wegner expansion restricted to two terms (leading and first confluent terms) is valid. The Ising-like preasymptotic domain description includes the correlations between parameters due to the error-bar determination of the exponents and amplitude combinations very close to the Wilson-Fisher fixed point. Adding the equivalent description of the mean field preasymptotic domain close to the Gaussian fixed point leads to define each mean crossover function with three calculated parameters. Fixing a unique value of one parameter whatever the selected mean crossover function, we use this parameter as a relative sensor to estimate the dominant nature, either (Ising-like) critical, or (mean-field-like) classical, of the crossover behavior. Finally, we obtain an explicit criterion to measure the extension of the Ising-like preasymptotic domain which can then permit to coherently account for measurements performed in systems where the asymptotical approach to the critical point remains finite, using a well-controlled number of system-dependent parameters (like in the subclass of one-component fluids).

Droplet motion with phase change in a temperature gradient

We examine the droplet motion in one-component fluids in a small temperature gradient by solving linearized hydrodynamic equations supplemented with appropriate surface boundary conditions. We show that the velocity field and the temperature around the droplet are strongly influenced by first-order phase transition taking place at the interface. Latent heat released or absorbed at the interface drastically changes the hydrodynamic flow around the droplet. As a result, the temperature becomes almost homogeneous inside the droplet and the Marangoni effect arising from the surface tension gradient is much suppressed. The droplet velocity is also much decelerated.

Ion-induced nucleation in polar one-component fluids

We present a Ginzburg-Landau theory of ion-induced nucleation in a gas phase of polar one-component fluids, where a liquid droplet grows with an ion at its center. By calculating the density profile around an ion, we show that the solvation free energy is larger in gas than in liquid at the same temperature on the coexistence curve. This difference much reduces the nucleation barrier in a metastable gas.

Predicting Mixture Phase Equilibria and Critical Behavior Using the SAFT-VRX Approach

The SAFT-VRX equation of state combines the SAFT-VR equation with a crossover function that smoothly transforms the classical equation into a nonanalytical form close to the critical point. By a combinination of the accuracy of the SAFT-VR approach away from the critical region with the asymptotic scaling behavior seen at the critical point of real fluids, the SAFT-VRX equation can accurately describe the global fluid phase diagram. In previous work, we demonstrated that the SAFT-VRX equation very accurately describes the pvT and phase behavior of both nonassociating and associating pure fluids, with a minimum of fitting to experimental data. Here, we present a generalized SAFT-VRX equation of state for binary mixtures that is found to accurately predict the vapor-liquid equilibrium and pvT behavior of the systems studied. In particular, we examine binary mixtures of n-alkanes and carbon dioxide + n-alkanes. The SAFT-VRX equation accurately describes not only the gas-liquid critical locus for these systems but also the vapor-liquid equilibrium phase diagrams and thermal properties in single-phase regions.

Simplified crossover droplet model for adsorption of pure fluids in slit pores

We present a generalized crossover (GC) model for the excess adsorption of pure fluids at a flat solid-liquid interface, which reproduces scaling behavior of the excess adsorption in the critical region and is reduced to the classical, van der Waals-type analytical model far away from the bulk critical point. In developing this model, we used the density-functional theory (DFT) approach for the order parameter profile calculations with a generalized corresponding states model for the local free-energy density. The GC DFT model well represents the available experimental adsorption data for Kr/graphite, C2H4/graphite, C3H8/graphite, CO2/silica, and SF6/graphite systems in the entire density range 0 < rho < or = 3rhoc and temperatures up to 1.7Tc. In the critical region 0.5 rhoc < r < or = 1.5rhoc and T < or = 1.15Tc, the GC DFT model is consistent with the predictions of the asymptotic renormalization-group crossover model for the critical adsorption in a semi-infinite system developed earlier. For the excess adsorption on the critical isochore, both theories predict a scaling-law behavior Gamma proportional tau(-nu+beta), but fail to reproduce a "critical depletion" of the excess adsorption along the critical isochore of the SF6/graphite system near Tc. We show that an anomalous decrease of adsorption observed in this system at tau = T/Tc - 1 < 10(-2) can be explained by finite-size effect and develop a simplified crossover droplet (SCD) model for the excess adsorption in a slit pore. With the effective size of the pore of L = 50 nm, the SCD model reproduces all available experimental data for SF6/graphite, including the critical isochore data where tau-->0, within experimental accuracy. At L >> xib (where xib is a bulk correlation length) the SCD model is transformed into the GC DFT model for semi-infinite systems. Application of the SCD model to the excess adsorption of carbon dioxide on the silica gel is also discussed.