No Arabic abstract
Remaining within the pure hydrodynamic approach, we formulate a self-consistent model for simulating the dynamic behavior of matter passing through metastable states in the two-phase liquid-vapor region of the phase diagram. The model is based on the local criterion of explosive boiling, derived by applying the theory of homogeneous bubble nucleation in superheated liquids. Practical application of the proposed model is illustrated with hydrodynamic simulations of a volumetrically uniformly heated planar layer of fused silica SiO2. Implications for experimentally measurable quantities are briefly discussed. A newly developed equation of state, based on the well known QEOS model and capable of handling homogeneous mixtures of elements, was used in the numerical simulations.
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 study an instability of thin liquid-vapor layers bounded by rigid parallel walls from both below and above. In this system, the interfacial instability is induced by lateral vapor pressure fluctuation, which is in turn attributed to the effect of phase change: evaporation occurs at a hotter portion of the interface and condensation at a colder one. The high vapor pressure pushes the interface downward and the low one pulls it upward. A set of equations describing the temporal evolution of the interface of the liquid-vapor layers is derived. This model neglects the effect of mass loss or gain at the interface and guarantees the mass conservation of the liquid layer. The result of linear stability analysis of the model shows that the presence of the pressure dependence of the local saturation temperature mitigates the growth of long-wave disturbances. The thinner vapor layer enhances the vapor pressure effect. We find the stability criterion, which suggests that only slight temperature gradients are sufficient to overcome the gravitational effect for a water/vapor system. The same holds for the Rayleigh-Taylor unstable case, with a possibility that the vapor pressure effect may be weakened if the accommodation coefficient is below a certain critical value.
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.
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.
We study the time-dependent neutralization of a slow highly charged ion that penetrates a hexagonal hollow-centred graphene nanoflake. To compute the ultrafast charge transfer dynamics, we apply an effective Hubbard nanocluster model and use the method of nonequilibrium Green functions (NEGF) in conjunction with an embedding self-energy scheme which allows one to follow the temporal changes of the number of electrons in the nanoflake. We perform extensive simulations of the charge transfer dynamics for a broad range of ion charge states and impact velocities. The results are used to put forward a simple semi-analytical model of the neutralization dynamics that is in very good agreement with transmission experiments, in which highly charged xenon ions pass through sheets of single-layer graphene.