Stationary longitudinal vortical rolls emerge in katabatic and anabatic Prandtl slope flows at shallow slopes as a result of an instability when the imposed surface buoyancy flux relative to the background stratification is sufficiently large. Here, we identify the self-pairing of these longitudinal rolls as a unique flow structure. The topology of the counter-rotating vortex pair bears a striking resemblance to speaker-wires and their interaction with each other is a precursor to further destabilization and breakdown of the flow field into smaller structures. On its own, a speaker-wire vortex retains its unique topology without any vortex reconnection or breakup. For a fixed slope angle $\alpha =3^{\circ }$ and at a constant Prandtl number, we analyse the saturated state of speaker-wire vortices and perform a bi-global linear stability analysis based on their stationary state. We establish the existence of both fundamental and subharmonic secondary instabilities depending on the circulation and transverse wavelength of the base state of speaker-wire vortices. The dominance of subharmonic modes relative to the fundamental mode helps to explain the relative stability of a single vortex pair compared to the vortex dynamics in the presence of two or an even number of pairs. These instability modes are essential for the bending and merging of multiple speaker-wire vortices, which break up and lead to more dynamically unstable states, eventually paving the way for transition towards turbulence. This process is demonstrated via three-dimensional flow simulations with which we are able to track the nonlinear temporal evolution of these instabilities.
more »
« less
Impact of inlet gas turbulence on the formation, development and breakup of interfacial waves in a two-phase mixing layer
Understanding the development and breakup of interfacial waves in a two-phase mixing layer between the gas and liquid streams is paramount to atomization. Due to the velocity difference between the two streams, the shear on the interface triggers a longitudinal instability, which develops to interfacial waves that propagate downstream. As the interfacial waves grow spatially, transverse modulations arise, turning the interfacial waves from quasi-two-dimensional to fully three-dimensional. The inlet gas turbulence intensity has a strong impact on the interfacial instability. Therefore, parametric direct numerical simulations are performed in the present study to systematically investigate the effect of the inlet gas turbulence on the formation, development and breakup of the interfacial waves. The open-source multiphase flow solver, PARIS, is used for the simulations and the mass–momentum consistent volume-of-fluid method is used to capture the sharp gas–liquid interfaces. Two computational domain widths are considered and the wide domain will allow a detailed study of the transverse development of the interfacial waves. The dominant frequency and spatial growth rate of the longitudinal instability are found to increase with the inlet gas turbulence intensity. The dominant transverse wavenumber, determined by the Rayleigh–Taylor instability, scales with the longitudinal frequency, so it also increases with the inlet gas turbulence intensity. The holes formed in the liquid sheet are important to the disintegration of the interfacial waves. The hole formation is influenced by the inlet gas turbulence. As a result, the sheet breakup dynamics and the statistics of the droplets formed also change accordingly.
more »
« less
- Award ID(s):
- 1942324
- NSF-PAR ID:
- 10329037
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 921
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Aerobreakup of liquid drops are important to many droplet applications, such as fuel injection. When a liquid drop is subjected to a gas stream of high velocity, the drop can deform and break into small droplets. The drop aerobreakup is controlled by multiple dimensionless parameters. The Weber number (We) has been commonly used to characterize the different breakup regimes. While the effects of Weber and Ohnesorge numbers on the aerobreakup of a drop in unbounded domain have been extensively studied, the effect of the Reynolds number (Re) based on gas properties are less understood and will be investigated by 2D axis-symmetric and 3D detailed numerical simulations in the present paper. Attention will be focused on the moderate We regime, where the drop mostly breaks in the bag mode. In many previous studies for millimeter drops, Re is too large to be relevant. However, for applications where drops are small and the relative velocity is high, Re can be quite small when the drop breaks. Parametric simulations of Re and We are performed to systematically investigate the effect of Re on the drop aerobreakup dynamics. The simulations are performed using the Basilisk solver, where the mass-momentum consistent VOF method is used to capture the interfacial dynamics on an adaptive mesh. The reduced Re is found to induce significant changes in the drop acceleration, deformation, bag morphology, and the bag breakup dynamics, which in turn lead to significant variation in the size and spatial distributions of the children droplets formed.more » « less
-
null (Ed.)Abstract ‘Horizontal convection’ (HC) is the generic name for the flow resulting from a buoyancy variation imposed along a horizontal boundary of a fluid. We study the effects of rotation on three-dimensional HC numerically in two stages: first, when baroclinic instability is suppressed and, second, when it ensues and baroclinic eddies are formed. We concentrate on changes to the thickness of the near-surface boundary layer, the stratification at depth, the overturning circulation and the flow energetics during each of these stages. Our results show that, for moderate flux Rayleigh numbers ( $O(1{0}^{11} )$ ), rapid rotation greatly alters the steady-state solution of HC. When the flow is constrained to be uniform in the transverse direction, rapidly rotating solutions do not support a boundary layer, exhibit weaker overturning circulation and greater stratification at all depths. In this case, diffusion is the dominant mechanism for lateral buoyancy flux and the consequent buildup of available potential energy leads to baroclinically unstable solutions. When these rapidly rotating flows are perturbed, baroclinic instability develops and baroclinic eddies dominate both the lateral and vertical buoyancy fluxes. The resulting statistically steady solution supports a boundary layer, larger values of deep stratification and multiple overturning cells compared with non-rotating HC. A transformed Eulerian-mean approach shows that the residual circulation is dominated by the quasi-geostrophic eddy streamfunction and that the eddy buoyancy flux has a non-negligible interior diabatic component. The kinetic and available potential energies are greater than in the non-rotating case and the mixing efficiency drops from ${\sim }0. 7$ to ${\sim }0. 17$ . The eddies play an important role in the formation of the thermal boundary layer and, together with the negatively buoyant plume, help establish deep stratification. These baroclinically active solutions have characteristics of geostrophic turbulence.more » « less
-
Experiments were carried out to observe the flow inside counterflow atomizers over a range of operating conditions and fluid properties. Liquids used were water and propylene glycol, while the gas was either air or helium. Liquid flow rates ranged from 10 ml/min to 40 ml/min, with gas liquid ratio (GLR) ranging from 0.1 to 0.6. The primary experiments used the 7-BM line of the Advanced Photon Source in Argonne National Laboratories with a 2.6 mm atomizer produced from (Poly)Ethyl-Ether-Ketone (PEEK). The X-Ray beam was operated in phase contrast mode, leading to interference patterns near the gas-liquid interface and enabling a qualitative understanding of the flow structure. Complementary optical work applied laser shadowgraphy to a 1 mm orifice atomizer constructed with quartz capillary tubing. A diffuse pulsed Nd:YAG laser backlight captured instantaneous gas-liquid interface positions in the internal flow. With both techniques, two distinct flow behaviors are observed corresponding to low and high GLR values. At low GLR, the inertia of the injected gas is insufficient to penetrate the liquid downflow. The gas stream entering the mixing chamber in the upstream direction is immediately deflected by the denser liquid and enters the discharge tube around a central liquid jet, which is sheared and accelerated by the surrounding gas, leading to breakup. A distinct frequency of jet breakup is observed inside the discharge tube, with the liquid jet oscillating and fragmenting against the walls. The situation at high GLR is quite different, however, as the incoming gas stream asymmetrically penetrates upstream into the mixing chamber, taking the form of a high-speed jet confined along one wall, and displaying a flapping instability as it encounters the liquid flowing downstream. This flapping causes violent mixing, resulting in a highly disturbed interface, along with the generation of liquid ligaments and gas bubbles. This two-phase mixture enters the discharge tube with no liquid jet formation evident for this case. The transition between these two regimes is explored by changing the liquid viscosity and gas molar mass, and weak sensitivity to fluid properties is observed. Further, quantitative image analysis techniques applied to the low and high GLR cases allow extraction of the frequencies of the liquid jet in the discharge tube at low GLR, as well as the flapping mode at high GLR.more » « less
-
more » « less
Abstract Collisional growth of cloud droplets is an essential yet uncertain process for drizzle and precipitation formation. To improve the quantitative understanding of this key component of cloud‐aerosol‐turbulence interactions, observational studies of collision‐coalescence in a controlled laboratory environment are needed. In an existing convection‐cloud chamber (the Pi Chamber), collisional growth is limited by low liquid water content and short droplet residence times. In this work, we use numerical simulations to explore various configurations of a convection‐cloud chamber that may intensify collision‐coalescence. We employ a large‐eddy simulation (LES) model with a size‐resolved (bin) cloud microphysics scheme to explore how cloud properties and the intensity of collision‐coalescence are affected by the chamber size and aspect ratio, surface roughness, side‐wall wetness, side‐wall temperature arrangement, and aerosol injection rate. Simulations without condensation and evaporation within the domain are first performed to explore the turbulence dynamics and wall fluxes. The LES wall fluxes are used to modify the Scalar Flux‐budget Model, which is then applied to demonstrate the need for non‐uniform side‐wall temperature (two side walls as warm as the bottom and the two others as cold as the top) to maintain high supersaturation in a tall chamber. The results of LES with full cloud microphysics reveal that collision‐coalescence is greatly enhanced by employing a taller chamber with saturated side walls, non‐uniform side‐wall temperature, and rough surfaces. For the conditions explored, although lowering the aerosol injection rate broadens the droplet size distribution, favoring collision‐coalescence, the reduced droplet number concentration decreases the frequency of collisions.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2021.481
