No Arabic abstract
Motived by recent ground-based and microgravity experiments investigating the interfacial dynamics of a volatile liquid (FC-72, $Pr=12.34$) contained in a heated cylindrical cell, we numerically study the thermocapillary-driven flow in such an evaporating liquid layer. Particular attention is given to the prediction of the transition of the axisymmetric flow to fully three-dimensional patterns when the applied temperature increases. The numerical simulations rely on an improved one-sided model of evaporation by including heat and mass transfer through the gas phase via the heat transfer Biot number and the evaporative Biot number. We present the axisymmetric flow characteristics, show the variation of the transition points with these Biot numbers, and more importantly elucidate the twofold role of the latent heat of evaporation in the stability; evaporation not only destabilizes the flow but also stabilizes it, depending upon the place where the evaporation-induced thermal gradients come into play. We also show that buoyancy in the liquid layer has a stabilizing effect, though its effect is insignificant. At high Marangoni numbers, the numerical simulations revealed smaller-scale thermal patterns formed on the surface of a thinner evaporating layer, in qualitative agreement with experimental observations. The present work helps to gain a better understanding of the role of a phase change in the thermocapillary instability of an evaporating liquid layer.
In this work we consider a new class of oscillatory instabilities that pertain to thermocapillary destabilization of a liquid film heated by a solid substrate. We assume the substrate thickness and substrate-film thermal conductivity ratio are large so that the effect of substrate thermal diffusion is retained at leading order in the long-wave approximation. As a result, system dynamics are described by a nonlinear partial differential equation for the film thickness that is nonlocally coupled to the full substrate heat equation. Perturbing about a steady quiescent state, we find that its stability is described by a non-self adjoint eigenvalue problem. We show that, under appropriate model parameters, the linearized eigenvalue problem admits complex eigenvalues that physically correspond to oscillatory (in time) instabilities of the thin film height. As the principal results of our work, we provide a complete picture of the susceptibility to oscillatory instabilities for different model parameters. Using this description, we conclude that oscillatory instabilities are more relevant experimentally for films heated by insulating substrates. Furthermore, we show that oscillatory instability where the fastest-growing (most unstable) wavenumber is complex, arises only for systems with sufficiently large substrate thicknesses.
We study the conductive and convective states of phase-change of pure water in a rectangular container where two opposite walls are kept respectively at temperatures below and above the freezing point and all the other boundaries are thermally insulating. The global ice content at the equilibrium and the corresponding shape of the ice-water interface are examined, extending the available experimental measurements and numerical simulations. We first address the effect of the initial condition, either fully liquid or fully frozen, on the system evolution. Secondly, we explore the influence of the aspect ratio of the cell, both in the configurations where the background thermal-gradient is antiparallel to the gravity, namely the Rayleigh-Benard (RB) setting, and when they are perpendicular, i.e., vertical convection (VC). We find that for a set of well-identified conditions the system in the RB configuration displays multiple equilibrium states, either conductive rather than convective, or convective but with different ice front patterns. The shape of the ice front appears to be always determined by the large scale circulation in the system. In RB, the precise shape depends on the degree of lateral confinement. In the VC case the ice front morphology is more robust, due to the presence of two vertically stacked counter-rotating convective rolls for all the studied cell aspect-ratios.
We consider the unsteady regimes of an acoustically-driven jet that forces a recirculating flow through successive reflections on the walls of a square cavity. The specific question being addressed is to know whether the system can sustain states of low-dimensional chaos when the acoustic intensity driving the jet is increased, and, if so, to characterise the pathway and underlying physical mechanisms. We adopt two complementary approaches, both based on data extracted from numerical simulations: (i) We first characterise successive bifurcations through the analysis of leading frequencies. Two successive phases in the evolution of the system are singled out in this way, both leading to potentially chaotic states. The two phases are separated by a drastic simplification of the dynamics that immediately follows the emergence of intermittency. The second phase also features a second intermediate state where the dynamics is simplified due to frequency-locking. (ii) Nonlinear time series analysis enables us to reconstruct the attractor of the underlying dynamical system, and to calculate its correlation dimension and leading Lyapunov exponent. Both these quantities bring confirmation that the state preceding the dynamic simplification that initiates the second phase is chaotic. Poincare maps further reveal that this chaotic state in fact results from a dynamic instability of the system between two non-chaotic states respectively observed at slightly lower and slightly higher acoustic forcing.
In a cylindrical container filled with an eutectic GaInSn alloy, an electro-vortex flow (EVF) is generated by the interaction of a non-uniform current with its own magnetic field. In this paper, we investigate the EVF phenomenon numerically and experimentally. Ultrasound Doppler Velocimetry (UDV) is applied to measure the velocity field in a cylindrical vessel. Second, we enhance an old numerical solver by taking into account the effect of Joule heating, and employ it for the numerical simulation of the EVF experiment. Special focus is laid on the role of the magnetic field, which is the combination of the current induced magnetic field and the external geomagnetic field. For getting a higher computational efficiency, the so-called parent-child mesh technique is applied in OpenFOAM when computing the electric potential, the current density and the temperature in the coupled solid-liquid conductor system. The results of the experiment are in good agreement with those of the simulation. This study may help to identify the factors that are essential for the EVF phenomenon, and for quantifying its role in liquid metal batteries.
A thin liquid film falling on a uniformly heated horizontal plate spreads into fingering ripples that can display a complex dynamics ranging from continuous waves, nonlinear spatially localized periodic wave patterns (i.e. rivulet structures) to modulated nonlinear wavetrain structures. Some of these structures have been observed experimentally, however conditions under which they form are still not well understood. In this work we examine profiles of nonlinear wave patterns formed by a thin liquid film falling on a uniformly heated horizonal plate. In this purpose, the Benney model is considered assuming a uniform temperature distribution along the film propagation on the horizontal surface. It is shown that for strong surface tension but relatively small Biot number, spatially localized periodic-wave structures can be analytically obtained by solving the governing equation under appropriate conditions. In the regime of weak nonlinearity, a multiple-scale expansion combined with the reductive perturbation method leads to a complex Ginzburg-Landau equation, the solutions of which are modulated periodic pulse trains which amplitude, width and period are expressed in terms of characteristic parameters of the model.