Ivana Seric

Direct numerical simulation of variable surface tension flows using a Volume-of-Fluid method

Abstract We develop a general methodology for the inclusion of a variable surface tension coefficient into a Volume-of-Fluid based Navier–Stokes solver. This new numerical model provides a robust and accurate method for computing the surface gradients directly by finding the tangent directions on the interface using height functions. The implementation is applicable to both temperature and concentration dependent surface tension coefficient, along with the setups involving a large jump in the temperature between the fluid and its surrounding, as well as the situations where the concentration should be strictly confined to the fluid domain, such as the mixing of fluids with different surface tension coefficients. We demonstrate the applicability of our method to the thermocapillary migration of bubbles and the coalescence of drops characterized by a different surface tension coefficient.

Exploiting the Marangoni Effect To Initiate Instabilities and Direct the Assembly of Liquid Metal Filaments

Utilization of the Marangoni effect in a liquid metal is investigated, focusing on initiating instabilities to direct material assembly via the Rayleigh-Plateau instability. Thin (2 nm) copper (Cu) films are lithographically patterned onto thick (12 nm) nickel (Ni) strips to induce a surface energy gradient at the maximum wavelength of the filament instability predicted by Rayleigh-Plateau instability analysis. The pattern is irradiated with an 18 ns pulsed laser such that the pattern melts and the resultant Ni-Cu surface tension gradient induces Marangoni flows due to the difference in surface energies. The experimental results, supported by extensive direct numerical simulations, demonstrate that the Marangoni flow exceeds the capillary flow induced by the initial geometry, guiding instabilities such that final nanoparticle location is directed toward the regions of higher surface energy (Ni regions). Our work shows a route for manipulation, by means of the Marangoni effect, to direct the evolution of the surface instabilities and the resulting pattern formation.

Instability of Gravity Driven Flow of Liquid Crystal Films

This paper discusses modeling of spreading nematic liquid crystal films. We concentrate on gravity driven spreading and consider various instabilities which occur during the spreading. We find that nematic character of the spreading film leads to stronger instabilities of the film fronts, and that it also leads to surface instabilities. We also present results of physical experiments involving spreading nematic films and find good agreement with the theoretical and computational predictions.

Influence of thermal effects on stability of nanoscale films and filaments on thermally conductive substrates

We consider films and filaments of nanoscale thickness on thermally conductive substrates exposed to external heating. Particular focus is on metal films exposed to laser irradiation. Due to short length scales involved, the absorption of heat in the metal is directly coupled to the film evolution, since the absorption length and the film thickness are comparable. Such a setup requires self-consistent consideration of fluid mechanical and thermal effects. We approach the problem via Volume-of-Fluid based simulations that include destabilizing liquid metal-solid substrate interaction potentials. These simulations couple fluid dynamics directly with the spatio-temporal evolution of the temperature field both in the fluid and in the substrate. We focus on the influence of the temperature variation of material parameters, in particular of surface tension and viscosity. Regarding variation of surface tension with temperature, the main finding is that while Marangoni effect may not play a significant role in the considered setting, the temporal variation of surface tension (modifying normal stress balance) is significant and could lead to complex evolution including oscillatory evolution of the liquid metal-air interface. Temperature variation of film viscosity is also found to be relevant. Therefore, the variations of surface tensions and viscosity could both influence the emerging wavelengths in experiments. In contrast, the filament geometry is found to be much less sensitive to a variation of material parameters with temperature.