skip to main content


Title: Direct Numerical Simulation of Bubble Formation Through a Submerged “Flute” With Experimental Validation
Abstract Direct numerical simulation (DNS) is often used to uncover and highlight physical phenomena that are not properly resolved using other computational fluid dynamics methods due to shortcuts taken in the latter to cheapen computational cost. In this work, we use DNS along with interface tracking to take an in-depth look at bubble formation, departure, and ascent through water. To form the bubbles, air is injected through a novel orifice geometry not unlike that of a flute submerged underwater, which introduces phenomena that are not typically brought to light in conventional orifice studies. For example, our single-phase simulations show a significant leaning effect, wherein pressure accumulating at the trailing nozzle edges leads to asymmetric discharge through the nozzle hole and an upward bias in the flow in the rest of the pipe. In our two-phase simulations, this effect is masked by the surface tension of the bubble sitting on the nozzle, but it can still be seen following departure events. After bubble departure, we observe the bubbles converge toward an ellipsoidal shape, which has been validated by experiments. As the bubbles rise, we note that local variations in the vertical velocity cause the bubble edges to flap slightly, oscillating between relatively low and high velocities at the edges.  more » « less
Award ID(s):
2107901
NSF-PAR ID:
10315390
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Fluids Engineering
Volume:
144
Issue:
2
ISSN:
0098-2202
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Extensive research has been conducted to resolve small-scale microlayer and bubble nucleation and departure processes in flow boiling, building on controlled pool boiling studies. Large-scale two-phase flow structures, such as Taylor bubbles, are known to locally modify transport due to their wakes and varying surrounding liquid film thickness. However, the effect of interaction of such large-scale flow processes with bubble nucleation is not yet well characterized. Wakes may drive premature nucleating bubble departure, or conversely, suppress boiling due to boundary layer quenching, significantly affecting overall heat transfer. To explore such phenomena, a two-phase flow boiling visualization facility is developed to collect simultaneous high-speed visualization and infrared (IR) thermal imaging temperature distribution data. The test cell channel is 420 mm long with a 10 mm × 10 mm internal square-cross section. A transparent conductive indium tin oxide (ITO) coated sapphire window serves as a heater and IR interface for measuring the internal wall temperature. The facility is charged with a low boiling point fluid (HFE7000) to reduce uncertainties from heat loss to the laboratory environment. Vertical saturated flow boiling wake-nucleation interaction experiments are performed for varying liquid volume flow rates (0.5 − 1.5 L min-1, laminar-to-turbulent Re) and heat fluxes (0 − 100 kW m-2). Discrete vapor slugs are injected to explore interactions with nucleate boiling processes. By measuring film heater power, surface temperature distributions, and pressures, local instantaneous heat transfer coefficients (HTC) can be obtained. Results will be applied to assess simulations at matched conditions for void fraction, and size statistics of flow structures. 
    more » « less
  2. Abstract

    The ignition of plasmas in liquids has applications from medical instrumentation to manipulation of liquid chemistry. Formation of plasmas directly in a liquid often requires prohibitively large voltages to initiate breakdown. Producing plasma streamers in bubbles submerged in a liquid with higher permittivity can significantly lower the voltage needed to initiate a discharge by reducing the electric field required to produce breakdown. The proximity of the bubble to the electrodes and the shape of the bubbles play critical roles in the manner in which the plasma is produced in, and propagates through, the bubble. In this paper, we discuss results from a three-dimensional direct numerical simulation (DNS) used to investigate the shapes of bubbles formed by injection of air into water. Comparisons are made to results from a companion experiment. A two-dimensional plasma hydrodynamics model was then used to capture the plasma streamer propagation in the bubble using a static bubble geometry generated by the DNS The simulations showed two different modes for streamer formation depending on the bubble shape. In an elliptical bubble, a short electron avalanche triggered a surface ionization wave (SIWs) resulting in plasma propagating along the surface of the bubble. In a circular bubble, an electron avalanche first traveled through the middle of the bubble before two SIWs began to propagate from the point closest to the grounded electrode where a volumetric streamer intersected the surface. In an elliptical bubble approaching a powered electrode in a pin-to-pin configuration, we experimentally observed streamer behavior that qualitatively corresponds with computational results. Optical emission captured over the lifetime of the streamer curve along the path of deformed bubbles, suggesting propagation of the streamer along the liquid/gas boundary interface. Plasma generation supported by the local field enhancement of the deformed bubble surface boundaries is a mechanism that is likely responsible for initiating streamer formation.

     
    more » « less
  3. Needle-free jet injections are actuated by a pressure impulse that can be delivered by different mechanisms to generate high-speed jets (Vj ∼ 102 m/s). During filling and transportation of disposable cartridges and ampoules, bubbles can form, which can be problematic especially for viscous fluids. Here, we report on the effect of location and size of entrapped air pockets in cartridges used in spring-powered jet injections. As air bubbles pass through the orifice, they undergo depressurization, which results in intermittent atomization and spray formation, temporarily increasing the jet dispersion. Atomization and dispersion of the jet can lead to product loss during an injection. We find that the effect of bubble location on the jet exit speed, delivery efficiency, and the projected area of the blebs formed after the injection was statistically significant (p < 0.05). The findings of this study have implications for the development of pre-filled cartridges for jet injection applications. 
    more » « less
  4. The effect of nozzle surface features on the overall atomization behavior of a liquid jet is analyzed in the present computational work by adopting three representative geometries, namely a single X-ray tomography scan of the Engine Combustion Network’s Spray A nozzle (Unprocessed), a spline reconstruction of multiple scans (Educated), and a purely external flow configuration. The latter configuration is often used in fundamental jet atomization studies. Numerically, the two-phase flow is solved based on algebraic volume-of-fluid methodology utilizing the OpenFoam solver, interFoam. Quantitative characterization of the surface features concerning the first two geometries reveals that while both of them have similar levels of cylindrical asymmetries, the nozzle configuration pertaining to the Unprocessed geometry has much larger surface features along the streamwise direction than the Educated geometry. This produces for the Unprocessed configuration a much larger degree of non-axial velocity components in the flow exiting the orifice and also a more pronounced disturbance of the liquid surface in the first few diameters downstream of the nozzle orifice. Furthermore, this heightened level of surface destabilization generates a much shorter intact liquid core length, that is, it produces faster primary atomization. The surprising aspect of this finding is that the differences between the Unprocessed and Educated geometries are of [Formula: see text](1) μm, and they are able to produce [Formula: see text](1) mm effects in the intact liquid core length. In spite of more pronounced atomization for the Unprocessed geometry, the magnitude of the turbulent liquid kinetic energy is roughly the same as the Educated geometry. This highlights the important role of mean field quantities, in particular, non-axial velocity components, in precipitating primary atomization. At the other end of the spectrum, the external-only configuration has the mildest level of surface disturbances in the near field resulting in the longest intact liquid core length. 
    more » « less
  5. The dynamics of a three-phase gas–liquid–liquid multiphase system is examined by direct numerical simulations. The system consists of a continuous liquid phase, buoyant gas bubbles, and smaller heavy drops that fall relative to the continuous liquid. The computational domain is fully periodic, and a force equal to the weight of the mixture is added to keep it in place. The governing parameters are selected so that the terminal Reynolds numbers of the bubbles and the drops are moderate; while the effect of bubble deformability is examined by changing its surface tension, the surface tension for the drops is sufficiently high so they do not deform. One bubble in a “unit cell” and eight freely interacting bubbles are examined. The dependency of the slip velocities, the velocity fluctuations, and the distribution of the dispersed phases on the volume fraction of each phase are examined. It is found that while the distribution of drops around a single bubble in a “unit cell” is uneven and depends on its deformability, the distribution of drops around freely interacting bubbles is relatively uniform for the parameters examined in this study. 
    more » « less