Enhancement of boiling heat transfer on biphilic (mixed-wettability) surfaces faces a sudden reversal at low pressures, which is brought about by excessive contact-line spreading across the wetting heterogeneities. We employ the diffuse-interface approach to numerically study bubble expansion on a heating surface that consists of opposing wettabilities. The results show a dramatic shift in the dynamics of a traversing contact line across the wettability divide under different gravities, which correspond to variable bubble growth rates. Specifically, it is found that the contact-line propagation tends to follow closely the rapidly expanding bubble at low gravity, with only a brief interruption at the border between the hydrophobic and hydrophilic sections of the surface. Only when the bubble growth becomes sufficiently weakened at high gravity does the contact line get slowed down drastically to the point of being nearly immobilized at the edge of the hydrophilic surface. The following bubble expansion, which faces strong limitations in the direction parallel to the surface, features a consistent apparent contact angle at around 66.4 degrees, regardless of the wettability combination. A simple theoretical model based on the force-balance analysis is proposed to describe the physical mechanism behind such a dramatic transition in the contact-line behavior.

In this paper, we numerically study pool boiling of a binary (water and nitrogen) mixture on a surface endowed with a combination of hydrophobicity and hydrophilicity (i.e., the so called biphilic surface). Here we adopt a numerical approach based on the phase field theory, where the vapor-liquid interface is assumed to be of a finite thickness (hence diffusive in nature) and requires no explicit tracking schemes. The theoretical modeling of two-phase heat and mass transfer in water diluted with nitrogen demonstrates the signiant impact of impurities on bubble dynamics. The simulations show that locally high concentrations of nitrogen gas within the vapor bubble is essential to weakening the condensation effect, which results in sustained bubble growth and ultimately (partial) departure from the surface under the artificially enlarged gravity. Simply increasing the solubility of nitrogen in water, however, turns out to be counterproductive because possible re-dissolution of the aggregated nitrogen by the bulk water could deprive the bubble of vital gas contents, leading instead to continuous bubble shrinkage and collapse. Additionally, it is found that with the significant accumulation of nitrogen, the bubble interface is increasingly dominated by a strong interfacial thermocapillary flow due to the Marangoni effect.

For phase-change cooling schemes for electronics, quick activation of nucleate boiling helps safeguard the electronics components from thermal shocks associated with undesired surface superheating at boiling incipience, which is of great importance to the long-term system stability and reliability. Previous experimental studies show that bubble nucleation can occur surprisingly early on mixed-wettability surfaces. In this paper, we report unambiguous evidence that such unusual bubble generation at extremely low temperatures-even below the boiling point-is induced by a significant presence of incondensable gas retained by the hydrophobic surface, which exhibits exceptional stability even surviving extensive boiling deaeration. By means of high-speed imaging, it is revealed that the consequently gassy boiling leads to unique bubble behaviour that stands in sharp contrast with that of pure vapour bubbles. Such findings agree qualitatively well with numerical simulations based on a diffuse-interface method. Moreover, the simulations further demonstrate strong thermocapillary flows accompanying growing bubbles with considerable gas contents, which is associated with heat transfer enhancement on the biphilic surface in the low-superheat region.

Multicomponent, multiphase, compressible flows are very important in real life, as well as in scientific research, while their modeling is in an early stage. In this paper, we propose a diffuse interface model for compressible binary mixtures, based on the balance of mass, momentum, energy, and the second law of thermodynamics. We show both analytically and numerically that this model is able to describe the phase equilibrium for a real binary mixture (CO2 + ethanol is considered in this paper) very well by adjusting the parameter which measures the attraction force between molecules of the two components in the model. We also show that the calculated surface tension of the CO2 + ethanol mixture at different concentrations match measurements in the literature when the mixing capillary coefficient is taken to be the geometric mean of the capillary coefficient of each component. Three different cases of two droplets in a shear flow, with the same or different concentration, are simulated, showing that the higher concentration of CO2 the smaller the surface tension and the easier the drop deforms.

Axisymmetric rapid expansion of supercritical carbon dioxide is investigated in this article. The extended generalized Bender equation of state is used to give a good description of the fluids over a wide range of pressure and temperature conditions. The locations of Mach disks are analyzed and compared with an experimental correlation for the case where there is no plate positioned in front of the nozzle exit. It is found that the disagreement between our numerical results and the experimental formula is very small when the pressure ratio is small, and increases as the pressure ratio increases. It is also found that with different equations of state, the predicted positions of Mach disks do not differ a lot, but the temperature profiles in the chamber differ a lot. The case where there is a plate positioned in front of the nozzle exit is also studied in this article. A universal similarity solution is obtained.

We numerically study the thermohydrodynamics of boiling for a CO2 + ethanol mixture on lyophilic and lyophobic surfaces in both closed and open systems, based on a diffuse interface model for a two-component system. The corresponding wetting boundary conditions for an isothermal system are proposed and verified in this paper. New phenomena due to the addition of another component, mainly the preferential evaporation of the more volatile component, are observed. In the open system and the closed system, the physical process shows very different characteristics. In the open system, except for the movement of the contact line, the qualitative features are rather similar for lyophobic and lyophilic surfaces. In the closed system, the vortices that are observed on a lyophobic surface are not seen on a lyophilic surface. More sophisticated wetting boundary conditions for nonisothermal, two-component systems might need to be further developed, taking into account the variations of density, temperature, and surface tension near the wall, while numerical results show that the boundary conditions proposed here also work well even in boiling, where the temperature is nonuniform.

In this paper, we numerically study particle formation in the rapid expansion of supercritical solution (RESS) process in a two dimensional, axisymmetric geometry, for a benzoic acid + CO2 system. The fluid is described by the classical Navier-Stokes equation, with the thermodynamic pressure being replaced by a generalized pressure tensor. Homogenous particle nucleation, transport, condensation and coagulation are described by a general dynamic equation, which is solved using the method of moments. The results show that the maximal nucleation rate and number density occurs near the nozzle exit, and particle precipitation inside the nozzle might not be ignored. Particles grow mainly across the shocks. Fluid in the shear layer of the jet shows a relatively low temperature, high nucleation rate, and carries particles with small sizes. On the plate, particles within the jet have smaller average size and higher geometric mean, while particles outside the jet shows a larger average size and a lower geometric mean. Increasing the preexpansion temperature will increase both the average particle size and standard deviation. The preexpansion pressure does not show a monotonic dependency with the average particle size. Increasing the distance between the plate and the nozzle exit might decrease the particle size. For all the cases in this paper, the average particle size on the plate is on the order of tens of nanometers.