Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “drops & bubbles”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

18 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Pressure fields produced by single-bubble collapse near a corner

Damage produced by repeated bubble collapse to neighboring rigid objects in hydraulic systems is an important consequence of cavitation. Although bubble collapse near a single wall has received significant attention in the past, few studies exist on the dynamics of bubbles collapsing near a corner, i.e., two flat rigid surfaces intersecting at a right angle. In this work we quantify the pressure fields produced by a single bubble collapsing near two perpendicular rigid walls. Using a high-order accurate shock- and interface-capturing method to solve the three-dimensional compressible Navier-Stokes equations for gas and liquid flows, we simulate the dynamics of a single bubble collapsing at different initial stand-off distances from the two walls. In contrast to a bubble collapsing near a single wall, the collapse of bubbles within a critical stand-off distance is not symmetric about the bisecting plane due to the interaction between the bubble and the second wall. The second wall affects the pressure produced during the collapse in the following ways: (i) For bubbles initially located sufficiently close to both walls, the reentrant jet produced during collapse no longer points in the direction normal to the closest wall but at an angle toward the corner, (ii) the part of the emitted shock with the highest amplitude propagates in line with the jet, and (iii) the bubble migrates in that same direction during its collapse with a dependence on the stand-off distance, consistent with predictions made using Kelvin impulse. The location of maximum pressure along the walls is measured for the different initial stand-off distances. Using acoustic arguments, we find a semiempirical relationship to predict the initial stand-off distances for which the maximum pressure occurs in the corner. Here we find that when the bubble is sufficiently close to equidistant from each boundary, the maximum pressure is observed in the corner due to the water-hammer and implosion shocks reflecting off the boundaries and intersecting in the corner. We also show that when bubbles are initially attached to either wall the wall pressure produced can be significantly increased compared to bubbles detached from either wall.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Ray-tracing image simulations of transparent objects with complex shape and inhomogeneous refractive index

Optical images of transparent three-dimensional objects can be different from a replica of the object’s cross section in the image plane, due to refraction at the surface or in the body of the object. Simulations of the object’s image are thus needed for the visualization and validation of physical models. We report ray tracing image simulations that achieved high physical fidelity, reproducing optical behaviors and image features not rendered in previous studies. We replicated brightfield microscopy images of drops with complex shapes, and images of pressure and shock waves traveling inside them. For high physical fidelity, the simulations must replicate the spatial and angular distribution of illumination rays, and both the experiment and the simulation must be designed for accurate optical modeling. The simulations are highly sensitive to the properties of the drops and can be used to diagnose and refine fluid dynamics models. The simulated images can also be optimized to extract multiple 3D properties from experimental images. Compared to specialized single-shot 3D imaging methods, this approach has the advantage that it preserves the experimental simplicity, the high resolution, and the visual interpretability characteristic to basic optical imaging. The techniques introduced here are directly applicable to optical microscopy, so they can be used in other fields, such as microfluidics and biology, to expand the type and the accuracy of three-dimensional information that can be extracted from basic optical images.

Cavitation↗
DOE OSTI2023

Rocket drops: The self-propulsion of supercooled freezing drops

Isolated supercooled water drops have been observed to move spontaneously while freezing in vacuum. This motion is caused by an increase of the evaporation rate during the early stages of freezing, which transfers momentum to the drops, similar to rocket propulsion. As in other cases of self-propulsion, symmetry-breaking is necessary, and occurs when ice nucleation occurs away from the center of the drop. The self-propulsion velocity was modeled analytically for a simplified case, and numerically for two experiments. Here, the model predicts peak velocities on the order of 1 m/s in vacuum, and a drop kinematics similar to that observed experimentally. In air, the self-propulsion velocity is expected to be much smaller but may be detectable for micron-sized drops.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Ephemeral antibubbles: Spatiotemporal evolution from direct numerical simulations

Antibubbles, which consist of a shell of a low-density fluid inside a high-density fluid, have several promising applications. We show, via extensive direct numerical simulations (DNSs), in both two and three dimensions, that the spatiotemporal evolution of antibubbles can be described naturally by the coupled Cahn-Hilliard-Navier-Stokes (CHNS) equations for a binary fluid. Our DNSs capture elegantly the gravity-induced thinning and breakup of an antibubble via the time evolution of the Cahn-Hilliard scalar-order-parameter field $\phi$, which varies continuously across interfaces, so we do not have to enforce complicated boundary conditions at the moving antibubble interfaces. To ensure that our results are robust, we supplement our CHNS simulations with sharp-interface volume-of-fluid DNSs. We track the thickness of the antibubble and calculate the dependence of the lifetime of an antibubble on several parameters; we show that our DNS results agree with various experimental results; in particular, the velocity with which the arms of the antibubble retract after breakup scales as σ 1/2 , where σ is the surface tension.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Competing mechanisms at vibrated interfaces of density-contrast fluids

Fluid-fluid interfacial instability and subsequent fluid mixing are ubiquitous in nature and engineering. The hydrodynamic instability of fluid interfaces has long centered on the pressure gradient-driven long-wavelength Rayleigh-Taylor instability and the resonance-induced short-wavelength Faraday instability. However, neither instability alone can explain the dynamics when both mechanisms are present. We identify a previously unseen multi-modal instability emerging from their coexistence. When the denser fluid is polydimethylsiloxane, the mixed region at a high density contrast (Atwood number = 0.9) spans a vibration amplitude range approximately twice the gravitational acceleration. Using Floquet stability analysis, we show how vibrations govern transitions between the RT and Faraday instabilities, leading to contention between these instabilities rather than resonant enhancement. Here, the initial transient growth is represented by the exponential modal growth of the most unstable Floquet exponent, along with its accompanying periodic behavior. Direct numerical simulations validate these findings and track interface breakup into the multiscale and nonlinear regimes. Specifically, we show that growing RT modes nonlinearly suppresses Faraday responses even when the initial growth rate of the Faraday instability is 3.63 times that of RT, so a bidirectional competition hinders their sustained coexistence.

Direct numerical simulations↗
DOE OSTI2023

Sizing and Selection of Pressure Relief Valves for High-Pressure Thermal–Hydraulic Systems

This study covers the critical concerns related to the sizing, selection, installation, maintenance, and testing of pressure safety valves (PSVs). The aim is to ensure the safety of pressurized systems, hydrostatic transmission systems, and hydraulic plants, including process plants, thermal power plants, and nuclear reactor systems. PSVs are devices that ensure the safety and reliability of pressurized vessels, lines, and systems during overpressure events. The task of selecting which PSV features are of greatest value for a specific purpose is complex—especially in the design of a high-pressure experimental thermal–hydraulic facility for hydrostatic and transient testing of the reactor system—when the systems are in the design and development phases and require qualification and demonstration to prove that they have reached a given level of technological readiness. The present study highlights the required steps for users to follow the associated rules, guidelines, and recommendations. As a part of this research, case studies are presented to help readers better understand the applicable strategy and standards. A discussion and a review of PSV performance degradation and failure are summarized to provide a better understanding of varied process applications and conditions, including fluid flow dynamics, boundary-layer formation and pressure drops, gas bubble formation and collapse, geometric configurations, inlet/outlet piping, abrupt pressure fluctuations, and acoustic resonance. Moreover, this study discusses the servicing and testing of PSVs in a multiphase pressurized system. Overall, it provides a basic overview of how PSVs ensure the safety of pressurized systems, supported by case studies and industrial practices.

02 PETROLEUM↗
DOE OSTI2024

Impulse-driven release of gas-encapsulated drops

Gas-encapsulated drops, much like antibubbles, are drops enclosed in a bubble within a liquid. They show potential as payload carriers in fluid transport and mixing techniques where sound waves can be leveraged to induce the collapse of the gas core and the subsequent release of the drop. Here, the interaction of millimetre-sized gas-encapsulated drops with impulsive laser-induced shock waves is investigated to gain fundamental insights into the release process. Experimental synchrotron X-ray phase contrast imaging, which allows the drop dynamics to be visualised inside the encapsulating bubble, is complemented by numerical simulations to study the intricate physics at play. Three drop dynamical release regimes are discovered, namely the drop impact , partial deposition and jet impact regimes. The regime type is mainly dependent on the shape of the bubble interface impacting the drop and the associated Weber and Reynolds numbers. The drop dynamics of the drop impact and partial deposition regimes show similarities with the canonical configuration of drops impacting flat liquid surfaces, whereas the jet impact regime resembles binary drop collisions, which allows existing scaling laws to be applied to describe the underlying processes. The release of the drop is investigated numerically. The time evolution of the drop dissemination within the surrounding liquid discloses enhanced mixing for dynamics involving high Weber and Reynolds numbers such as the drop impact and jet impact regimes.

42 ENGINEERING↗
DOE OSTI2024

Kelvin-Helmholtz instability as one of the key features for fast and efficient emulsification by hydrodynamic cavitation

The paper investigates the oil–water emulsification process inside a micro-venturi channel. More specifically, the possible influence of Kelvin-Helmholtz instability on the emulsification process. High-speed visualizations were conducted inside a square venturi constriction with throat dimensions of 450 µm by 450 µm, both under visible light and X-Rays. We show that cavity shedding caused by the instability results in the formation of several cavity vortices. Their rotation causes the deformation of the oil stream into a distinct wave-like shape, combined with fragmentation into larger drops due to cavitation bubble collapse. Later on, the cavity collapse further disperses the larger drops into a finer emulsion. Thus, it turns out that the Kelvin-Helmholtz instability is similarly characteristic for hydrodynamic cavitation emulsification inside a microchannel as is the Rayleigh-Taylor instability for acoustically driven emulsion formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Shock-induced bubble jets: a dual perspective of bubble collapse and interfacial instability theory

Interactions between shock waves and gas bubbles in a liquid can lead to bubble collapse and high-speed liquid jet formation, relevant to biomedical applications such as shock wave lithotripsy and targeted drug delivery. This study reveals a complex interplay between acceleration-induced instabilities that drive jet formation and radial accelerations causing overall bubble collapse under shock wave pressure. Using high-speed synchrotron X-ray phase contrast imaging, the dynamics of micrometre-sized air bubbles interacting with laser-induced underwater shock waves are visualised. These images offer full optical access to phase discontinuities along the X-ray path, including jet formation, its propagation inside the bubble, and penetration through the distal side. Jet formation from laser-induced shock waves is suggested to be an acceleration-driven process. A model predicting jet speed based on the perturbation growth rate of a single-mode Richtmyer–Meshkov instability shows good agreement with experimental data, despite uncertainties in the jet-driving mechanisms. The jet initially follows a linear growth phase, transitioning into a nonlinear regime as it evolves. To capture this transition, a heuristic model bridging the linear and nonlinear growth phases is introduced, also approximating jet shape as a single-mode instability, again matching experimental observations. Upon piercing the distal bubble surface, jets can entrain gas and form a toroidal secondary bubble. Linear scaling laws are identified for the pinch-off time and volume of the ejected bubble relative to the jet’s Weber number, characterising the balance of inertia and surface tension. At low speeds, jets destabilise due to capillary effects, resulting in ligament pinch-off.

Drops and Bubbles: Bubble dynamics↗
DOE OSTI2024

Muon Flux Variance from Severe Atmospheric Conditions

We analyze data from the European Centre for Medium-Range Weather Forecasts (ECMWF) and the QuarkNet Research Group at Fermilab to study the variance in muon flux due to extreme weather conditions. We considered three high-intensity storms for each season: winter, spring, and summer. The most significant observation was during winter, with a ~125% increase in muon flux corresponding to a temperature drop of ~20K. Additionally, we are developing CosmicWatch muon detectors to be used in tests for the MMBC/MAMBA Bubble Chamber. A pair of these detectors will be employed to trigger on muons, providing extremely precise and inexpensive tracking granularity in the form of a muon telescope.

Aamir, Yusuf↗
DOE OSTI2025

Active microrheology and dynamic phases for pattern-forming systems with competing interactions

We consider the driven dynamics of a probe particle moving through an assembly of particles with competing long-range repulsive and short-range attractive interactions, which form crystal, stripe, labyrinth, and bubble states as the ratio of attraction to repulsion is varied. We show that the probe particle exhibits a depinning-like threshold from an elastic regime, where the probe particle is trapped by interactions with the other particles, to a plastic flow regime, where the probe particle can break bonds in the surrounding medium. For a fixed particle density, the depinning threshold and sliding velocity of the probe particle vary nonmonotonically as the attraction term is increased. A velocity minimum appears near the crystal to stripe crossover, and there is a significant increase in the depinning threshold in the bubble regime when the probe particle is strongly confined inside the bubbles. For fixed attractive interaction but increasing particle density, the behavior is also nonmonotonic, and there are jumps and drops in the velocity and depinning threshold corresponding to points at which the system transitions between different structures. There are also several distinct flow states that can be characterized by the amount of plastic deformation induced by the probe particle in the surrounding medium. Each flow state generates a different amount of effective drag on the probe particle, and there can be jumps in the velocity-force curve at transitions between the states. In conclusion, we also find that when oriented stripes are present, the probe particle can move along the stripe in an edge transport state that has a finite Hall angle.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Methane foam performance in oil-wet unconsolidated porous media: A systematic experimental investigation at reservoir conditions

Methane foam performance was systematically evaluated using high-pressure and high-temperature experiments performed on oil-wet proppant packs. The sensitivities of foam performance to various key factors, including surfactant concentration, foam quality, total injection rate, salinity, operating pressure, and initial oil saturation were examined. The results showed that even though the foamability and foam strength in oil-wet proppant packs were adversely affected by wettability and presence of oil, the steady-state pressure drops across the proppant packs were insensitive to the quantity of oil initially present. The presence of oil above the critical oil saturation impeded the onset of foam generation. This suggested that the foamability of surfactants depends on the synergetic effects of surfactant’s ionic nature and foam parameters, such as foam quality and total injection rate, on their ability to reduce the oil saturation to a threshold value below which favorable bubble generation sites could be created. It was observed that, for low injection rate, anionic surfactant was ineffective at low salinity conditions, while amphoteric agent performed well at all salinities. Additionally, an increase in salinity resulted in enhanced pressure drop across the proppant packs for both anionic and amphoteric surfactants. Further, the results also established that too high operating pressures are somewhat detrimental to the foam strength and stability, irrespective of the ionic nature of foaming agents and the wettability of the porous medium. For the anionic surfactant, the increase in the total flow rate was found to mitigate the adverse effects of low salt concentration and operating pressure. Results from quality-sensitivity foam tests identified two regions for the foam behavior, named low-quality and high-quality regimes, separated by the transition foam quality. Finally, we discuss the implications of the results for the design and effective implementation of enhanced oil recovery (EOR) schemes that deploy hydrocarbon-based foams.

02 PETROLEUM↗
DOE OSTI2026

Rapid Electrochemical Carburization of Mild Steel

Next-generation iron production generates iron with 0 wt.% C, in contrast to contemporary blast furnace production of iron with 1 - 5 wt.% C. This makes new iron feedstocks difficult to use as a drop-in replacement in current steelmaking. We report results of an emissions-negative approach of carburizing iron or mild steel using molten salt electrochemistry and CO2. Electrochemical carburization of mild steel via molten carbonate salt electrolyte and CO2 bubbling shows rapid carbide development tunable by applied electric potential. Mild steel rods, serving as the cathode, show pearlite development at the surface after exposure to applied voltage in a LiKCO3 molten carbonate salt bath at 715 degrees C. The thickness of the pearlitic layer could be varied by changing the applied potential, with a strong difference in microstructure observed between an applied potential of -1.8 V and -2.4 V. This rapid carburization has multiple applications for enhancing steel production by upgrading the carbon content of low-carbon feedstocks.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Enhanced rare earth element recovery with cross-linked glutaraldehyde-lanthanide binding peptides in foam-based separations

Lanthanide Binding Tag (LBT) peptides that coordinate selectively with lanthanide ions can be used to replace the energy intensive processes used for the separation of rare earth elements (REEs). These surface-active biomolecules, once selectively complexed with the trivalent REE cations, can adsorb to air/aqueous interfaces of bubbles for foam-based REEs recovery. Glutaraldehyde, an organic compound that is a homobifunctional crosslinker for proteins and peptides, can be used to enhance the adsorption and interfacial stabilization of lanthanide-bound peptides films. The stability of the interfacial cross-linked films was tested by measuring their dilational and shear surface rheological properties. Surface activity of the adsorbed species was analyzed using pendant drop tensiometry, while surface density and molecular arrangement were determined using x-ray reflectivity and x-ray fluorescence near total reflection. Glutaraldehyde cross-linked REE-peptide complexes enhance the adsorption of lanthanides to air-water interfaces, resulting in thicker interfacial structures. Subsequently, these thicker layers enhance the dilational and shear interfacial rheological properties. The interfacial film stabilization and REEs extraction promoted by the cross-linker presented in this work provides an approach to integrate glutaraldehyde as a substitute of common foam stabilizers such as polymers, surfactants, and particles to optimize the recovery of REEs when using biomolecules as extractants.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Jetted and Turbulent Stellar Deaths: New LVK-detectable Gravitational-wave Sources

Upcoming LIGO–Virgo–KAGRA (LVK) observing runs are expected to detect a variety of inspiralling gravitational-wave (GW) events that come from black hole and neutron star binary mergers. Detection of noninspiral GW sources is also anticipated. We report the discovery of a new class of noninspiral GW sources—the end states of massive stars—that can produce the brightest simulated stochastic GW burst signal in the LVK bands known to date, and could be detectable in LVK run A+. Some dying massive stars launch bipolar relativistic jets, which inflate a turbulent energetic bubble—cocoon—inside of the star. We simulate such a system using state-of-the-art 3D general relativistic magnetohydrodynamic simulations and show that these cocoons emit quasi-isotropic GW emission in the LVK band, ~10–100 Hz, over a characteristic jet activity timescale ~10–100 s. Our first-principles simulations show that jets exhibit a wobbling behavior, in which case cocoon-powered GWs might be detected already in LVK run A+, but it is more likely that these GWs will be detected by the third-generation GW detectors with an estimated rate of ~10 events yr –1 . The detection rate drops to ~1% of that value if all jets were to feature a traditional axisymmetric structure instead of a wobble. Accompanied by electromagnetic emission from the energetic core-collapse supernova and the cocoon, we predict that collapsars are powerful multimessenger events.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Batch VUV4 characterization for the SBC-LAr10 scintillating bubble chamber

The Scintillating Bubble Chamber (SBC) collaboration purchased 32 Hamamatsu VUV4 silicon photomultipliers (SiPMs) for use in SBC-LAr10, a bubble chamber containing 10 kg of liquid argon. A dark-count characterization technique, which avoids the use of a single-photon source, was used at two temperatures to measure the VUV4 SiPMs breakdown voltage (V BD ), the SiPM gain (g SiPM ), the rate of change of g SiPM with respect to voltage (m), the dark count rate (DCR), and the probability of a correlated avalanche (P CA ) as well as the temperature coefficients of these parameters. A Peltier-based chilled vacuum chamber was developed at Queen's University to cool down the Quads to 233.15 ± 0.2 K and 255.15 ± 0.2 K with average stability of ±20 mK. An analysis framework was developed to estimate V BD to tens of mV precision and DCR close to Poissonian error. The temperature dependence of V BD was found to be 56 ± 2 mV K -1 , and m on average across all Quads was found to be (459 ± 3(stat.)±23(sys.))× 10 3 e- PE -1 V -1 . The average DCR temperature coefficient was estimated to be 0.099 ± 0.008 K -1 corresponding to a reduction factor of 7 for every 20 K drop in temperature. The average temperature dependence of P CA was estimated to be 4000 ± 1000 ppm K -1 . P CA estimated from the average across all SiPMs is a better estimator than the P CA calculated from individual SiPMs, for all of the other parameters, the opposite is true. All the estimated parameters were measured to the precision required for SBC-LAr10, and the Quads will be used in conditions to optimize the signal-to-noise ratio.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗