No Arabic abstract
The injection of liquid fuel at supercritical pressures is a relevant topic in combustion, but usually overlooked. In the past, the wrong assumption whereby the liquid experiments a fast transition to a supercritical state was made, thus neglecting any role of two-phase interface dynamics in the early stages of the atomization process. However, recent studies have shown that local thermodynamic phase equilibrium and mixing between the involved species allow the coexistence of both phases in this pressure range. In this work, a Volume-of-Fluid method adapted to variable-density real fluids is used to solve the low-Mach-number governing equations coupled with a thermodynamic model based on the Soave-Redlich-Kwong equation of state. The mixing process, interface thermodynamics and early deformation of a cool liquid jet composed of n-decane surrounded by a hotter gas composed of oxygen at 150 bar are analyzed. Although heat conducts from the hotter gas into the liquid, net condensation can provide the proper local energy balance at high pressures. Then, vaporization and condensation may happen simultaneously at different interface locations. As pressure increases, liquid and gas mixtures become more alike in the vicinity of the interface. Thus, a combination of low surface tension force and gas-like liquid viscosities causes an early growth of surface instabilities. Early results indicate some similarity with high-Weber-number incompressible flows. The role of vortex dynamics on the interface deformation is analyzed by using the dynamical vortex identification method.
Numerical heat and mass transfer analysis of a configuration where a cool liquid hydrocarbon is suddenly introduced to a hotter gas at supercritical pressure shows that a well-defined phase equilibrium can be established before substantial growth of typical hydrodynamic instabilities. The equilibrium values at the interface quickly reach near-steady values. Sufficiently thick diffusion layers form quickly around the liquid-gas interface (e.g., 3-10 microns for the liquid phase and 10-30 microns for the gas phase in 10-100 microseconds), where density variations become increasingly important with pressure as mixing of species is enhanced. While the hydrocarbon vaporizes and the gas condenses for all analyzed pressures, the net mass flux across the interface reverses as pressure is increased, showing that a clear vaporization-driven problem at low pressures may present condensation at higher pressures. This is achieved while heat still conducts from gas to liquid. Analysis of fundamental thermodynamic laws on a fixed-mass element containing the diffusion layers proves the thermodynamic viability of the obtained results.
We study the Richtmyer--Meshkov (RM) instability of a relativistic perfect fluid by means of high order numerical simulations with adaptive mesh refinement (AMR). The numerical scheme adopts a finite volume Weighted Essentially Non-Oscillatory (WENO) reconstruction to increase accuracy in space, a local space-time discontinuous Galerkin predictor method to obtain high order of accuracy in time and a high order one-step time update scheme together with a cell-by-cell space-time AMR strategy with time-accurate local time stepping. In this way, third order accurate (both in space and in time) numerical simulations of the RM instability are performed, spanning a wide parameter space. We present results both for the case in which a light fluid penetrates into a higher density one (Atwood number $A>0$), and for the case in which a heavy fluid penetrates into a lower density one (Atwood number $A<0$). We find that, for large Lorentz factors gamma_s of the incident shock wave, the relativistic RM instability is substantially weakened and ultimately suppressed. More specifically, the growth rate of the RM instability in the linear phase has a local maximum which occurs at a critical value of gamma_s ~ [1.2,2]. Moreover, we have also revealed a genuine relativistic effect, absent in Newtonian hydrodynamics, which arises in three dimensional configurations with a non-zero velocity component tangent to the incident shock front. In this case, the RM instability is strongly affected, typically resulting in less efficient mixing of the fluid.
We report on a new class of electromagnetically-driven fluid interface instability. Using the optical radiation pressure of a cw laser to bend a very soft near-critical liquid-liquid interface, we show that it becomes unstable for sufficiently large beam power P, leading to the formation of a stationary beam-centered liquid micro-jet. We explore the behavior of the instability onset by tuning the interface softness with temperature and varying the size of the exciting beam. The instability mechanism is experimentally demonstrated. It simply relies on total reflection of light at the deformed interface whose condition provides the universal scaling relation for the onset Ps of the instability.
In this letter, we provide fundamental insights into the dynamic transcritical transition process using molecular dynamics simulations. A transcritical region, which covers three different fluid states, was discovered as a substitute for the traditional interface. The physical properties, such as temperature and density, exhibited a highly non-linear distribution in the transcritical region. Meanwhile, the surface tension was found to exist throughout the transcritical region, and the magnitude was directly proportional to $ - rho cdot { abla ^2}rho $
Flows forced by a precessional motion can exhibit instabilities of crucial importance, whether they concern the fuel of a flying object or the liquid core of a telluric planet. So far, stability analyses of these flows have focused on the special case of a resonant forcing. Here, we address the instability of the flow inside a precessing cylinder in the general case. We first show that the base flow forced by the cylinder precession is a superposition of a vertical or horizontal shear flow and an infinite sum of forced modes. We then perform a linear stability analysis of this base flow by considering its triadic resonance with two free Kelvin modes. Finally, we derive the amplitude equations of the free Kelvin modes and obtain an expression of the instability threshold and growth rate.