Dewetting acceleration by evaporation

Abstract

Dewetting of a substrate, i.e. the liquid film retraction under partial wetting conditions, has been studied extensively over the past decades. This theoretical work deals with the dewetting phenomenon in the presence of liquid evaporation into the pure vapour atmosphere driven by substrate superheating with respect to the saturation temperature corresponding to the vapour pressure. The dynamic profile of the vapour–liquid interface is analysed numerically by using the generalised lubrication approximation that can be applied to an interface slope larger than its conventional counterpart. This approximation uses a phenomenological argument to increase the precision of calculations. For small slopes, it reduces to the conventional lubrication approach. Several nanoscale effects are included in the theory to obtain more precise results: the Kelvin effect, the hydrodynamic slip, the Marangoni effect, the vapour recoil and the interfacial thermal resistance. Their relative importance is discussed. The results include the film shape evolution, the apparent contact angle and the contact line speed, all defined by the substrate superheating and its wetting properties. These dependencies agree with existing experimental results. It is found that the dewetting is accelerated by evaporation. The numerical results are compared to those of the multiscale theory that uses two main parameters, the Voinov angle (the apparent angle obtained within the microregion problem) and the Voinov length. The latter is studied for several fluids as a function of superheating, and an approximate expression for it is proposed. A good agreement between the numerics and the multiscale theory shows the validity of the latter. This means that the dewetting speed in the presence of evaporation can be calculated with an analytical expression that requires only the Voinov angle as a numerical input. Such results are important for many applications, for instance, bubble growth in boiling or film dynamics in heat pipes.

Time-averaged approach to the dewetting problem at evaporation

“…Another equation determines the CL speed U cl (which is a part of the problem). The relation between U cl and the evaporation flux at CL [3]…”

“…( 11). The variation in V turns out to be small during the film evaporation [3]. In applications (e.g., in pulsating heat pipe modeling), the focus is primarily on Ca and θ app , which experience slow variations during the dewetting process.…”

“…In these cases, the liquid film is surrounded by the pure vapor of the same fluid. The numerical results of [3] indicated that evaporation accelerates the dewetting; the contact line recedes faster when the solid substrate is heated. This letter is structured as follows.…”

“…Until recently, the dewetting studies [1,2] concerned the adiabatic conditions (in the absence of net heat and mass transfer between liquid, solid and ambient gas). The work of [3] investigated dewetting in the presence of phase change over the liquid film (evaporation/condensation). The co-occurrence of the two processes is frequently encountered in studies of heat exchange systems, such as (a) E-mail: vadim.nikolayev@cea.fr (corresponding author) Fig.…”

“…This letter is structured as follows. First, we will briefly summarize the problem in [3] because these results will help to understand the multiscale analysis discussed next. Based on it, new ("time-averaged") approach will be finally introduced.…”