Long-wave-instability-induced pattern formation in an evaporating sessile or pendent liquid layer

We investigate the nonlinear dynamics and stability of an evaporating liquid layer subject to vapor recoil, capillarity, thermocapillarity, ambient cooling, viscosity, and negative or positive gravity combined with buoyancy effects in the lubrication approximation. Using linear theory, we identify the mechanisms of finite-time rupture, independent of thermocapillarity and direction of gravity, and predict the effective growth rate of an interfacial perturbation which reveals competition among the mechanisms. A stability diagram is predicted for the onset of long-wave (LW) evaporative convection. In the two-dimensional simulation, we observe well-defined capillary ridges on both sides of the valley under positive gravity and main and secondary droplets under negative gravity, while a ridge can be trapped in a large-scale drained region in both cases. Neglecting the other non-Boussinesq effects, buoyancy does not have a significant influence on interfacial evolution and rupture time but makes contributions to the evaporation-driven convection and heat transfer. The average Nusselt number is found to increase with a stronger buoyancy effect. The flow field and interface profile jointly manifest the LW Marangoni-Rayleigh-Benard convection under positive gravity and the LW Marangoni convection under negative gravity. In the three-dimensional simulation of moderate evaporation with a random perturbation, the rupture patterns are characterized by irregular ridge networks with distinct height scales for positive and negative gravity. A variety of interfacial and internal dynamics are displayed, depending on evaporation conditions, gravity, Marangoni effect, and ambient cooling. Reasonable agreement is found between the present results and the reported experiments and simulations. The concept of dissipative compacton also sheds light on the properties of interfacial fractalization.