Marangoni Convection of Self-Rewetting Fluid Layers with a Deformable Interface in a Square Enclosure and Driven by Imposed Nonuniform Heat Energy Fluxes

• Published in: Energies (Basel) Vol. 18; no. 13; p. 3563
• Main Authors: Elbousefi, Bashir, Schupbach, William, Premnath, Kannan N.
• Format: Journal Article
• Language: English
• Published: Basel MDPI AG 01.07.2025
• Subjects: Analysis; Behavior; Boundary conditions; Distribution (Probability theory); Electric properties; Fluid dynamics; Heat transfer; Interfaces; Investigations; lattice Boltzmann method; Marangoni effects; Methods; Microelectromechanical systems; Microfluidics; Numerical analysis; self-rewetting fluids; Simulation methods; surface-tension-driven flow; Surfactants; thermocapillary convection; Nunavut
• ISSN: 1996-1073; 1996-1073
• DOI: 10.3390/en18133563

Fluids that exhibit self-rewetting properties, such as aqueous long-chain alcohol solutions, display a unique quadratic relationship between surface tension and temperature and are marked by a positive gradient. This characteristic leads to distinctive patterns of thermocapillary convection and associated interfacial dynamics, setting self-rewetting fluids apart from normal fluids (NFs). The potential to improve heat transfer using self-rewetting fluids (SRFs) is garnering interest for use in various technologies, including low-gravity conditions and microfluidic systems. Our research aims to shed light on the contrasting behaviors of SRFs in comparison to NFs regarding interfacial transport phenomena. This study focuses on the thermocapillary convection in SRF layers with a deformable interface enclosed inside a closed container modeled as a square cavity, which is subject to nonuniform heating, represented using a Gaussian profile for the heat flux variation on one of its sides, in the absence of gravity. To achieve this, we have enhanced a central-moment-based lattice Boltzmann method (LBM) utilizing three distribution functions for tracking interfaces, computing two-fluid motions with temperature-dependent surface tension and energy transport, respectively. Through numerical simulations, the impacts of several characteristic parameters, including the viscosity and thermal conductivity ratios, as well as the surface tension–temperature sensitivity parameters, on the distribution and magnitude of the thermocapillary-driven motion are examined. In contrast to that in NFs, the counter-rotating pair of vortices generated in the SRF layers, due to the surface tension gradient at the interface, is found to be directed toward the SRF layers’ hotter zones. Significant interfacial deformations are observed, especially when there are contrasts in the viscosities of the SRF layers. The thermocapillary convection is found to be enhanced if the bottom SRF layer has a higher thermal conductivity or viscosity than that of the top layer or when distributed, rather than localized, heating is applied. Furthermore, the higher the magnitude of the effect of the dimensionless quadratic surface tension sensitivity coefficient on the temperature, or of the effect of the imposed heat flux, the greater the peak interfacial velocity current generated due to the Marangoni stresses. In addition, an examination of the Nusselt number profiles reveals significant redistribution of the heat transfer rates in the SRF layers due to concomitant nonlinear thermocapillary effects.