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
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
- PAR ID:
- 10315390
- 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
-
-
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
-
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
-
Heat transfer due to the convective boiling mechanism in the microchannel plays an important role in heat transfer during boiling. Therefore, it is relevant to find ways to manipulate the vapor bubbles such that convection heat transfer is enhanced. This numerical study investigates the effects of different geometrical parameters on bubble movement through a micro tapered gap. The objective is to identify an optimal configuration such that the bubble moves at the fastest possible speed when it travels through the micro gap. To conduct this research a model is created using ANSYS-Fluent which uses the Volume of Fluid (VOF) interface tracking method. The multiphase VOF model tracks the air-water interface. A bubble is generated inside the microchannel in which fluid is flowing. The overall domain of the model consists of the surface at the bottom, having an orifice through which the air bubble is generated. Three different cases of an angled tapered surface are created 5°, 10°, and 15°. The airflow rate is kept constant throughout each simulation. Simulation results show the impact of the tapered angle on the bubble’s flow movement and flow direction. Liquid and air velocity contours can be used to analyze the flow. The impact of the taper angles on the movement and flow direction of the air bubble is discussed. It is observed that the performed simulations help to better understand the experimental observation of bubble motion; the simulations give clear evidence of the fluid dynamic behavior along the tapered microchannel.more » « less
-
We study bubble break-up in homogeneous and isotropic turbulence by direct numerical simulations of the two-phase incompressible Navier–Stokes equations. We create the turbulence by forcing in physical space and introduce the bubble once a statistically stationary state is reached. We perform a large ensemble of simulations to investigate the effect of the Weber number (the ratio of turbulent and surface tension forces) on bubble break-up dynamics and statistics, including the child bubble size distribution, and discuss the numerical requirements to obtain results independent of grid size. We characterize the critical Weber number below which no break-up occurs and the associated Hinze scale $$d_h$$ . At Weber number close to stable conditions (initial bubble sizes $$d_0\approx d_h$$ ), we observe binary and tertiary break-ups, leading to bubbles mostly between $$0.5d_h$$ and $$d_h$$ , a signature of a production process local in scale. For large Weber numbers ( $$d_0> 3d_h$$ ), we observe the creation of a wide range of bubble radii, with numerous child bubbles between $$0.1d_h$$ and $$0.3d_h$$ , an order of magnitude smaller than the parent bubble. The separation of scales between the parent and child bubble is a signature of a production process non-local in scale. The formation mechanism of these sub-Hinze scale bubbles relates to rapid large deformation and successive break-ups: the first break-up in a sequence leaves highly deformed bubbles which will break again, without recovering a spherical shape and creating an array of much smaller bubbles. We discuss the application of this scenario to the production of sub-Hinze bubbles under breaking waves.more » « less