Spreading with evaporation and condensation in one-component fluids

Transient coarsening and the motility of optically heated Janus colloids in a binary liquid mixture

A gold-capped Janus particle suspended in a near-critical binary liquid mixture can self-propel under illumination. We have immobilized such a particle in a narrow channel and carried out a combined experimental and theoretical study of the non-equilibrium dynamics of a binary solvent around it - lasting from the very moment of switching illumination on until the steady state is reached. In the theoretical study we use both a purely diffusive and a hydrodynamic model, which we solve numerically. Our results demonstrate a remarkable complexity of the time evolution of the concentration field around the colloid. This evolution is governed by the combined effects of the temperature gradient and the wettability, and crucially depends on whether the colloid is free to move or is trapped. For the trapped colloid, all approaches indicate that the early time dynamics is purely diffusive and characterized by composition layers travelling with constant speed from the surface of the colloid into the bulk of the solvent. Subsequently, hydrodynamic effects set in. Anomalously large nonequilibrium fluctuations, which result from the temperature gradient and the vicinity of the critical point of the binary liquid mixture, give rise to strong concentration fluctuations in the solvent and to permanently changing coarsening patterns not observed for a mobile particle. The early time dynamics around initially still Janus colloids produces a force which is able to set the Janus colloid into motion. The propulsion due to this transient dynamics is in the direction opposite to that observed after the steady state is attained.

Illumination-induced motion of a Janus nanoparticle in binary solvents

Using a fluid particle dynamics method we numerically investigate the motion of a spherical Janus particle suspended in a binary liquid mixture, which emerges under heating of one-half of a colloid surface. The method treats simultaneously the flow of the solvent and the motion of the particle, hence, the velocity of the particle can be computed directly. Our approach accounts for a phenomenon of critical adsorption, therefore, a particle that is adsorptionwise nonneutral is always completely covered by an adsorption layer (droplet). In order to establish the mechanism of self-propulsion, we study systematically various combinations of adsorption preference on both hemispheres of the Janus colloid as function of the heating power for symmetric and nonsymmetric binary solvents and for various particle sizes in three spatial dimensions. Only for a particle for which the heated hemisphere is neutral whereas the other hemisphere prefers one of the two components of the mixture does the reversal of the direction of motion occur. The particle self-propels much faster in nonsymmetric binary solvents. Self-propulsion originates from a gradient of mechanical stress, in a way similar to the Marangoni effect. This stress is not localized at the edge but distributed within the whole droplet. We compare our findings with the experimental observations and other theoretical results.

Liquid phase stabilization versus bubble formation at a nanoscale curved interface

We investigate the nature of vapor bubble formation near a nanoscale-curved convex liquid-solid interface using two models: an equilibrium Gibbs model for homogenous nucleation, and a nonequilibrium dynamic van der Waals-diffuse-interface model for phase change in an initially cool liquid. Vapor bubble formation is shown to occur for sufficiently large radius of curvature and is suppressed for smaller radii. Solid-fluid interactions are accounted for and it is shown that liquid-vapor interfacial energy, and hence Laplace pressure, has limited influence over bubble formation. The dominant factor is the energetic cost of creating the solid-vapor interface from the existing solid-liquid interface, as demonstrated via both equilibrium and nonequilibrium arguments.

Ballistic heat transport in laser generated nano-bubbles

Nanobubbles generated by laser heated plasmonic nanoparticles are of interest for biomedical and energy harvesting applications. Of utmost importance is the maximal size of these transient bubbles. Here, we report hydrodynamic phase field simulations of the dynamics of laser induced nanobubbles, with the aim to understand which physical processes govern their maximal size. We show that the nanobubble maximal size and lifetime are to a large extent controlled by the ballistic thermal flux which is present inside the bubble. Taking into account this thermal flux, we can reproduce the fluence dependence of the maximal nanobubble radius as reported experimentally. We also discuss the influence of the laser pulse duration on the number of nanobubbles generated and their maximal size. These studies represent a significant step toward the optimization of the nanobubble size, which is of crucial importance for photothermal cancer therapy applications.

Density functional theory of gas-liquid phase separation in dilute binary mixtures

We examine statics and dynamics of phase-separated states of dilute binary mixtures using density functional theory. In our systems, the difference of the solvation chemical potential between liquid and gas [Formula: see text] (the Gibbs energy of transfer) is considerably larger than the thermal energy [Formula: see text] for each solute particle and the attractive interaction among the solute particles is weaker than that among the solvent particles. In these conditions, the saturated vapor pressure increases by [Formula: see text], where [Formula: see text] is the solute density added in liquid. For [Formula: see text], phase separation is induced at low solute densities in liquid and the new phase remains in gaseous states, even when the liquid pressure is outside the coexistence curve of the solvent. This explains the widely observed formation of stable nanobubbles in ambient water with a dissolved gas. We calculate the density and stress profiles across planar and spherical interfaces, where the surface tension decreases with increasing interfacial solute adsorption. We realize stable solute-rich bubbles with radius about 30 nm, which minimize the free energy functional. We then study dynamics around such a bubble after a decompression of the surrounding liquid, where the bubble undergoes a damped oscillation. In addition, we present some exact and approximate expressions for the surface tension and the interfacial stress tensor.

Thermal singularity and contact line motion in pool boiling: Effects of substrate wettability

The dynamic van der Waals theory [Phys. Rev. E 75, 036304 (2007)] is employed to model the growth of a single vapor bubble in a superheated liquid on a flat homogeneous substrate. The bubble spreading dynamics in the pool boiling regime has been numerically investigated for one-component van der Waals fluids close to the critical point, with a focus on the effect of the substrate wettability on bubble growth and contact line motion. The substrate wettability is found to control the apparent contact angle and the rate of bubble growth (the rate of total evaporation), through which the contact line speed is determined. An approximate expression is derived for the contact line speed, showing good agreement with the simulation results. This demonstrates that the contact line speed is primarily governed by (1) the circular shape of interface (for slow bubble growth), (2) the constant apparent contact angle, and (3) the constant bubble growth rate. It follows that the contact line speed has a sensitive dependence on the substrate wettability via the apparent contact angle which also determines the bubble growth rate. Compared to hydrophilic surfaces, hydrophobic surfaces give rise to a thinner shape of bubble and a higher rate of total evaporation, which combine to result in a much faster contact line speed. This can be linked to the earlier formation of a vapor film and hence the onset of boiling crisis.

Nanobubbles around plasmonic nanoparticles: Thermodynamic analysis

We describe the dynamics of vapor nanobubbles in water, on the basis of simulations of a hydrodynamics phase-field model. This situation is relevant to recent experiments, where a water nanobubble is generated around a nanoparticle immersed in water, and heated by an intense laser pulse. We emphasize the importance of nanoscale effects in the dynamics of the nanobubble. We first analyze the evolution of the temperature inside the bubble. We show that the temperature drops by hundredths of kelvins in a few picoseconds, just after nanobubble formation. This is the result of the huge drop of the thermal boundary conductance between the nanoparticle and the fluid accompanying vaporization. Subsequently, the temperature inside the vapor is almost homogeneous and the temperature gradient is concentrated in the liquid, whose thermodynamic state locally follows the saturation line. We discuss also the evolution of the pressure inside the vapor nanobubble. We show that nanobubble generation is accompanied by a pressure wave propagating in the liquid at a velocity close to the liquid speed of sound. The internal pressure inside the vapor just after its formation largely exceeds Laplace pressure and quickly relaxes as a result of the damping generated by the viscous forces. All these considerations shed light on the thermodynamics of the nanobubbles generated experimentally.

Quantitative phase-field modeling for wetting phenomena

A new phase-field model is developed for studying partial wetting. The introduction of a third phase representing a solid wall allows for the derivation of a new surface tension force that accounts for energy changes at the contact line. In contrast to other multi-phase-field formulations, the present model does not need the introduction of surface energies for the fluid-wall interactions. Instead, all wetting properties are included in a unique parameter known as the equilibrium contact angle θeq. The model requires the solution of a single elliptic phase-field equation, which, coupled to conservation laws for mass and linear momentum, admits the existence of steady and unsteady compact solutions (compactons). The representation of the wall by an additional phase field allows for the study of wetting phenomena on flat, rough, or patterned surfaces in a straightforward manner. The model contains only two free parameters, a measure of interface thickness W and β, which is used in the definition of the mixture viscosity μ=μlϕl+μvϕv+βμlϕw. The former controls the convergence towards the sharp interface limit and the latter the energy dissipation at the contact line. Simulations on rough surfaces show that by taking values for β higher than 1, the model can reproduce, on average, the effects of pinning events of the contact line during its dynamic motion. The model is able to capture, in good agreement with experimental observations, many physical phenomena fundamental to wetting science, such as the wetting transition on micro-structured surfaces and droplet dynamics on solid substrates.

Single-bubble dynamics in pool boiling of one-component fluids

We numerically investigate the pool boiling of one-component fluids with a focus on the effects of surface wettability on the single-bubble dynamics. We employed the dynamic van der Waals theory [Phys. Rev. E 75, 036304 (2007)], a diffuse-interface model for liquid-vapor flows involving liquid-vapor transition in nonuniform temperature fields. We first perform simulations for bubbles on homogeneous surfaces. We find that an increase in either the contact angle or the surface superheating can enhance the bubble spreading over the heating surface and increase the bubble departure diameter as well and therefore facilitate the transition into film boiling. We then examine the dynamics of bubbles on patterned surfaces, which incorporate the advantages of both hydrophobic and hydrophilic surfaces. The central hydrophobic region increases the thermodynamic probability of bubble nucleation while the surrounding hydrophilic region hinders the continuous bubble spreading by pinning the contact line at the hydrophobic-hydrophilic intersection. This leads to a small bubble departure diameter and therefore prevents the transition from nucleate boiling into film boiling. With the bubble nucleation probability increased and the bubble departure facilitated, the efficiency of heat transfer on such patterned surfaces is highly enhanced, as observed experimentally [Int. J. Heat Mass Transfer 57, 733 (2013)]. In addition, the stick-slip motion of contact line on patterned surfaces is demonstrated in one-component fluids, with the effect weakened by surface superheating.

Kinetics of nanobubble generation around overheated nanoparticles

We report on the formation and growth of nanobubbles around laser-heated gold nanoparticles in water. Using a hydrodynamic free-energy model, we show that the temporal evolution of the nanobubble radius is asymmetrical: the expansion is found to be adiabatic, while the collapse is best described by an isothermal evolution. We unveil the critical role of the thermal boundary resistance in the kinetics of formation of the nanobubbles: close to the vapor production threshold, nanobubble generation is very long, yielding optimal conditions for laser-energy conversion. Furthermore, the long appearance times allow nanoparticle melting before the onset of vaporization.

Molecular dynamics simulations for the motion of evaporative droplets driven by thermal gradients along nanochannels

For a one-component fluid on a solid substrate, a thermal singularity may occur at the contact line where the liquid-vapor interface intersects the solid surface. Physically, the liquid-vapor interface is almost isothermal at the liquid-vapor coexistence temperature in one-component fluids while the solid surface is almost isothermal for solids of high thermal conductivity. Therefore, a temperature discontinuity is formed if the two isothermal interfaces are of different temperatures and intersect at the contact line. This leads to the so-called thermal singularity. The localized hydrodynamics involving evaporation/condensation near the contact line leads to a contact angle depending on the underlying substrate temperature. This dependence has been shown to lead to the motion of liquid droplets on solid substrates with thermal gradients (Xu and Qian 2012 Phys. Rev. E 85 061603). In the present work, we carry out molecular dynamics (MD) simulations as numerical experiments to further confirm the predictions made from our previous continuum hydrodynamic modeling and simulations, which are actually semi-quantitatively accurate down to the small length scales in the problem. Using MD simulations, we investigate the motion of evaporative droplets in one-component Lennard-Jones fluids confined in nanochannels with thermal gradients. The droplet is found to migrate in the direction of decreasing temperature of solid walls, with a migration velocity linearly proportional to the temperature gradient. This agrees with the prediction of our continuum model. We then measure the effect of droplet size on the droplet motion. It is found that the droplet mobility is inversely proportional to a dimensionless coefficient associated with the total rate of dissipation due to droplet movement. Our results show that this coefficient is of order unity and increases with the droplet size for the small droplets (~10 nm) simulated in the present work. These findings are in semi-quantitative agreement with the predictions of our continuum model. Finally, we measure the effect of liquid-vapor coexistence temperature on the droplet motion. Through a theoretical analysis on the size of the thermal singularity, it can be shown that the droplet mobility decreases with decreasing coexistence temperature. This is observed in our MD simulations.

Hydrodynamics of Leidenfrost droplets in one-component fluids

Using the dynamic van der Waals theory [Phys. Rev. E 75, 036304 (2007)], we numerically investigate the hydrodynamics of Leidenfrost droplets under gravity in two dimensions. Some recent theoretical predictions and experimental observations are confirmed in our simulations. A Leidenfrost droplet larger than a critical size is shown to be unstable and break up into smaller droplets due to the Rayleigh-Taylor instability of the bottom surface of the droplet. Our simulations demonstrate that an evaporating Leidenfrost droplet changes continuously from a puddle to a circular droplet, with the droplet shape controlled by its size in comparison with a few characteristic length scales. The geometry of the vapor layer under the droplet is found to mainly depend on the droplet size and is nearly independent of the substrate temperature, as reported in a recent experimental study [Phys. Rev. Lett. 109, 074301 (2012)]. Finally, our simulations demonstrate that a Leidenfrost droplet smaller than a characteristic size takes off from the hot substrate because the levitating force due to evaporation can no longer be balanced by the weight of the droplet, as observed in a recent experimental study [Phys. Rev. Lett. 109, 034501 (2012)].

Thermal singularity and droplet motion in one-component fluids on solid substrates with thermal gradients

Using a continuum model capable of describing the one-component liquid-gas hydrodynamics down to the contact line scale, we carry out numerical simulation and physical analysis for the droplet motion driven by thermal singularity. For liquid droplets in one-component fluids on heated or cooled substrates, the liquid-gas interface is nearly isothermal. Consequently, a thermal singularity occurs at the contact line and the Marangoni effect due to temperature gradient is suppressed. Through evaporation or condensation in the vicinity of the contact line, the thermal singularity makes the contact angle increase with the increasing substrate temperature. This effect on the contact angle can be used to move the droplets on substrates with thermal gradients. Our numerical results for this kind of droplet motion are explained by a simple fluid dynamical model at the droplet length scale. Since the mechanism for droplet motion is based on the change of contact angle, a separation of length scales is exhibited through a comparison between the droplet motion induced by a wettability gradient and that by a thermal gradient. It is shown that the flow field at the droplet length scale is independent of the statics or dynamics at the contact line scale.

Droplet motion in one-component fluids on solid substrates with wettability gradients

Droplet motion on solid substrates has been widely studied not only because of its importance in fundamental research but also because of its promising potentials in droplet-based devices developed for various applications in chemistry, biology, and industry. In this paper, we investigate the motion of an evaporating droplet in one-component fluids on a solid substrate with a wettability gradient. As is well known, there are two major difficulties in the continuum description of fluid flows and heat fluxes near the contact line of droplets on solid substrates, namely, the hydrodynamic (stress) singularity and thermal singularity. To model the droplet motion, we use the dynamic van der Waals theory [Phys. Rev. E 75, 036304 (2007)] for the hydrodynamic equations in the bulk region, supplemented with the boundary conditions at the fluid-solid interface. In this continuum hydrodynamic model, various physical processes involved in the droplet motion can be taken into account simultaneously, e.g., phase transitions (evaporation or condensation), capillary flows, fluid velocity slip, and substrate cooling or heating. Due to the use of the phase field method (diffuse interface method), the hydrodynamic and thermal singularities are resolved automatically. Furthermore, in the dynamic van der Waals theory, the evaporation or condensation rate at the liquid-gas interface is an outcome of the calculation rather than a prerequisite as in most of the other models proposed for evaporating droplets. Numerical results show that the droplet migrates in the direction of increasing wettability on the solid substrates. The migration velocity of the droplet is found to be proportional to the wettability gradients as predicted by Brochard [Langmuir 5, 432 (1989)]. The proportionality coefficient is found to be linearly dependent on the ratio of slip length to initial droplet radius. These results indicate that the steady migration of the droplets results from the balance between the (conservative) driving force due to the wettability gradient and the (dissipative) viscous drag force. In addition, we study the motion of droplets on cooled or heated solid substrates with wettability gradients. The fast temperature variations from the solid to the fluid can be accurately described in the present approach. It is observed that accompanying the droplet migration, the contact lines move through phase transition and boundary velocity slip with their relative contributions mostly determined by the slip length. The results presented in this paper may lead to a more complete understanding of the droplet motion driven by wettability gradients with a detailed picture of the fluid flows and phase transitions in the vicinity of the moving contact line.

Thermohydrodynamics of boiling in a van der Waals fluid

We present a modeling approach that enables numerical simulations of a boiling Van der Waals fluid based on the diffuse interface description. A boundary condition is implemented that allows in and out flux of mass at constant external pressure. In addition, a boundary condition for controlled wetting properties of the boiling surface is also proposed. We present isothermal verification cases for each element of our modeling approach. By using these two boundary conditions we are able to numerically access a system that contains the essential physics of the boiling process at microscopic scales. Evolution of bubbles under film boiling and nucleate boiling conditions are observed by varying boiling surface wettability. We observe flow patters around the three-phase contact line where the phase change is greatest. For a hydrophilic boiling surface, a complex flow pattern consistent with vapor recoil theory is observed.

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.

Phase separation in thermal systems: a lattice Boltzmann study and morphological characterization

We investigate thermal and isothermal symmetric liquid-vapor separations via a fast Fourier transform thermal lattice Boltzmann (FFT-TLB) model. Structure factor, domain size, and Minkowski functionals are employed to characterize the density and velocity fields, as well as to understand the configurations and the kinetic processes. Compared with the isothermal phase separation, the freedom in temperature prolongs the spinodal decomposition (SD) stage and induces different rheological and morphological behaviors in the thermal system. After the transient procedure, both the thermal and isothermal separations show power-law scalings in domain growth, while the exponent for thermal system is lower than that for isothermal system. With respect to the density field, the isothermal system presents more likely bicontinuous configurations with narrower interfaces, while the thermal system presents more likely configurations with scattered bubbles. Heat creation, conduction, and lower interfacial stresses are the main reasons for the differences in thermal system. Different from the isothermal case, the release of latent heat causes the changing of local temperature, which results in new local mechanical balance. When the Prandtl number becomes smaller, the system approaches thermodynamical equilibrium much more quickly. The increasing of mean temperature makes the interfacial stress lower in the following way: σ=σ(0)[(T(c)-T)/(T(c)-T(0))](3/2), where T(c) is the critical temperature and σ(0) is the interfacial stress at a reference temperature T(0), which is the main reason for the prolonged SD stage and the lower growth exponent in the thermal case. Besides thermodynamics, we probe how the local viscosities influence the morphology of the phase separating system. We find that, for both the isothermal and thermal cases, the growth exponents and local flow velocities are inversely proportional to the corresponding viscosities. Compared with the isothermal case, the local flow velocity depends not only on viscosity but also on temperature.

Phase-field model of solid-liquid phase transition with density difference and latent heat in velocity and elastic fields

We present a phase-field model of solid-liquid transitions with inhomogeneous temperature in one-component systems, including hydrodynamics and elasticity. Our model can describe plastic deformations at large elastic strains. We use it to investigate the melting of a solid domain, accounting for the latent heat effect, where there appears a velocity field in liquid and an elastic field in solid. We present simulation results in two dimensions for three cases of melting. First, a solid domain is placed on a heated wall, which melts mostly near the solid-liquid-wall contact region. Second, a solid domain is suspended in a warmer liquid under shear flow, which rotates as a whole because of elasticity and melts gradually. Cooling of the surrounding liquid is accelerated by convection. Third, a solid rod is under high compression in liquid, where slips appear from the solid-liquid interface, leading to a plastic deformation. Subsequently, melting starts in the plastically deformed areas, eventually resulting in the fracture of the rod into pieces. In these phase-transition processes, the interface temperature is kept nearly equal to the coexisting temperature T(cs)(p) away from the heated wall, but this local equilibrium is not attained near the the contact region. We also examine a first-order liquid-liquid phase transition under heating from a boundary in one-component liquids.