Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bubble dynamics”

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.

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Insights into the rapid two-phase transport dynamics in different structured porous transport layers of water electrolyzers through high-speed visualization

In proton exchange membrane electrolyzer cells (PEMECs), maintaining efficient two-phase transport is one of the most important functions of porous transport layers (PTLs). To enhance the two-phase transport in PTLs, thin/titanium liquid/gas diffusion layers (TT-LGDLs) are introduced in PEMECs, and their difference from the conventional Ti felt PTLs are analyzed in-situ through high-speed and microscale visualization and electrochemical characterizations. The visualization results show that unfavorable large slugs can be greatly reduced in the PEMEC with a TT-LGDL compared to the PEMEC with a Ti felt PTL. More importantly, the recovery capability of water starvation with different PTLs is studied. After water starvation, the PEMEC with the TT-LGDL can recover the water starvation much more rapidly than the Ti felt cell, benefiting from its short and straight-through flow paths. Furthermore, the TT-LGDL tends to generate oxygen bubbles that are almost six times smaller and 236 times more frequently than the Ti felt PTL, indicating significantly boosted removal efficiency of produced oxygen and PEMEC performance. Finally, this study offers new insights into the dynamic processes of two-phase transport and the recovery capability of water starvation for different PTLs, which will provide valuable guidance for further optimization of PTLs and performance enhancement of PEMECs.

08 HYDROGEN↗
DOE OSTI2026

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2023

Dynamics of oil–water interface at the beginning of the ultrasonic emulsification process

A lot of effort has been dedicated in recent years towards understanding the basics of cavitation induced emulsification, mainly in the form of single cavitation bubbles. Regarding bulk acoustic emulsification, a lot less research has been done. In our here presented work we utilize advanced high-speed observation techniques in visible light and X-Rays to build upon that knowledge and advance the understanding of bulk emulsion preparation. During research we discovered that emulsion formation has an acute impact on the behavior of the interface and more importantly on its position relative to the horn, hence their interdependence must be carefully studied. We did this by observing bulk emulsification with 2 cameras simultaneously and corroborating these measurements with observation under X-Rays. Since the ultrasonic horns location also influences interface behavior, we shifted its initial position to different locations nearer to and further away from the oil–water interface in both phases. We found that a few millimeters distance between the horn and interface is not enough for fine emulsion formation, but that they must be completely adjacent to each other, with the horn being located inside the oil–water interface. We also observed some previously undiscovered phenomena, such as the splitting of the interface to preserve continuous emulsion formation, climbing of the interface up the horn and circular interface protrusions towards the horn forming vertical emulsion streams. Interestingly, no visible W/O emulsion was ever formed during our experiments, only O/W regardless of initial horn position.

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

Experimental and numerical investigation of sands and Geldart A biomass co-fluidization

This study investigated the fluidization of sands and small Geldart A biomass mixtures. The mixture fluidized like Geldart A type particles with a uniform bed expansion regime before bubbling. The video recorded color distance between pure sands and sands–biomass mixtures was used to estimate the sands–biomass mixing. The coarse-grained computational fluid dynamics–discrete element method with a hybrid drag model which couples the Syamlal–O'Brien drag and a filtered drag can capture the mixing while the simulation with Gidaspow drag predicted a segregated bed. The simulations were further validated with experimental measured pressure drops. The time averaged pressure drop equals the weight of the bed material, however, its fluctuation is about three times of the bed material fluctuation.

42 ENGINEERING↗
NASA NTRS1996

Multiphase Flow: The Gravity of the Situation

A brief survey is presented of flow patterns in two-phase, gas-liquid flows at normal and microgravity, the differences between them being explored. It seems that the flow patterns in zero gravity are in general much simpler than those in normal gravity with only three main regimes (namely bubbly, slug and annular flows) being observed. Each of these three regimes is then reviewed, with particular reference to identification of areas of study where investigation of flows at microgravity might not only be interesting in themselves, but also throw light on mechanisms at normal earth gravity. In bubbly flow, the main area of interest seems to be that of bubble coalescence. In slug flow, the extension of simple displacement experiments to the zero gravity case would appear to be a useful option, supplemented by computational fluid dynamics (CFD) studies. For annular flow, the most interesting area appears to be the study of the mechanisms of disturbance waves; it should be possible to extend the region of investigation of the onset and behavior of these waves to much low gas velocities where measurements are clearly much easier.

Hewitt, Geoffrey F.↗
NASA NTRS2014

OASIS Observation and Analysis of Smectic Islands in Space

The Observation and Analysis of Smectic Islands in Space (OASIS) project comprises a series of experiments that will probe the interfacial and hydrodynamic behavior of freely suspended liquid crystal films in space. These are the thinnest known stable condensed phase structures, making them ideal for studies of fluctuation and interface phenomena. The experiments seek to verify theories of coarsening dynamics, hydrodynamic flow, relaxation of hydrodynamic perturbations, and hydrodynamic interactions of a near two-dimensional structure. The effects of introducing islands or droplets on a very thin bubble will be studied, both as controllable inclusions that modify the flow and as markers of flow.

liquid crystals↗
NASA NTRS2008

Non-steady Reconnection in Global Simulations of Magnetosphere Dynamics

To analyze the non-steady magnetic reconnection during quasi-steady solar wind driving we employed high resolution global MHD model BATSRUS with non-MHD corrections in diffusion regions around the reconnection sites. To clarify the role of small-scale non-MHD effects on the global magnetospheric dynamic we performed simulations with different models of dissipation. We found that magnetopause surface is not in steady state even during extended periods of steady solar wind conditions. The so-called tilted reconnection lines become unstable due to formation of pressure bubbles, strong core field flux tubes, vortices, and traveling magnetic field cavities. Non-steady dayside reconnection results in formation of flux tubes with bended axis magnetically connecting magnetic field cavities generated at flanks and strong core segments formed near the subsolar region. We found that the rate of magnetic flux loading to the tail lobes is not very sensitive to the dissipation mechanism and details of the dayside reconnection. On the other hand the magnetotail reconnection rate, the speed of the reconnection site retreat and the global magnetotail dynamics strongly depend on the model of dissipation. THEMIS and Cluster observations are consistent with signatures predicted by simulations.

Kuznetsova, M. M.↗
NASA NTRS1988

The effects of compressibility on dynamic stall

In this paper, typical computational predictions and experimental measurements of compressible flow past an airfoil at dynamic stall conditions are studied and compared to develop an insight into the effect of compressibility on dynamic stall. The dependency of the critical Mach number on airfoil leading edge curvature, camber, and angle of attack is investigated. Evidence is presented to show that a local region of supersonic flow occurs on an oscillating airfoil, even for a freestream Mach number as low as 0.2, if the boundary layer remains attached and the angle of attack is sufficiently high; that a shock terminates this local supersonic bubble; and that the vorticity that this shock generates grows rapidly and becomes very unstable as the angle of attack increases beyond the value at which the maximum local flow speed first exceeds the speed of sound. It is suggested that these shock-induced effects compete with the dynamic viscous effects occurring in the boundary layer in determining the onset of separation, which can lead to premature dynamic stall and can significantly reduce the maximum dynamic lift that can otherwise be obtained.

Fung, K.-Y.↗
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 OSTI2021

Optical Imaging on Z LDRD: Design and Development of Self-Emission and Debris Imagers

We present an overview of the design and development of optical self-emission and debris imaging diagnostics for the Z Machine at Sandia National Laboratories. These diagnostics were designed and implemented to address several gaps in our understanding of visibly emitting phenomenon on Z and the post-shot debris environment. Optical emission arises from plasmas that form on the transmission line that delivers energy to Z loads and on the Z targets themselves; however, the dynamics of these plasmas are difficult to assess without imaging data. Addressing this, we developed a new optical imager called SEGOI (Self-Emission Gated Optical Imager) that leverages the eight gated optical imagers and two streak cameras of the Z Line VISAR system. SEGOI is a low cost, side-on imager with a 1 cm field of view and 30-50 µm spatial resolution, sensitive to green light (540-600 nm). This report outlines the design considerations and development of this diagnostic and presents an overview of the first diagnostic data acquired from four experimental campaigns. SEGOI was fielded on power flow experiments to image plasmas forming on and between transmission lines, on an inertial confinement fusion experiment called the Dynamic Screw Pinch to image low density plasmas forming on return current posts, on an experiment designed to measure the magneto Rayleigh-Taylor instability to image the instability bubble trajectory and self-emission structures, and finally on a Magnetized Liner Inertial Fusion (MagLIF) experiment to image the emission from the target. The second diagnostic developed, called DINGOZ (Debris ImagiNG on Z), was designed to improve our understanding of the post-shot debris environment. DINGOZ is an airtight enclosure that houses electronics and batteries to operate a high-speed (10-400 kfps) camera in the Z Machine center section. We report on the design considerations of this new diagnostic and present the first high-speed imaging data of the post-shot debris environment on Z.

42 ENGINEERING↗
NASA NTRS2003

Model of Fluidized Bed Containing Reacting Solids and Gases

A mathematical model has been developed for describing the thermofluid dynamics of a dense, chemically reacting mixture of solid particles and gases. As used here, "dense" signifies having a large volume fraction of particles, as for example in a bubbling fluidized bed. The model is intended especially for application to fluidized beds that contain mixtures of carrier gases, biomass undergoing pyrolysis, and sand. So far, the design of fluidized beds and other gas/solid industrial processing equipment has been based on empirical correlations derived from laboratory- and pilot-scale units. The present mathematical model is a product of continuing efforts to develop a computational capability for optimizing the designs of fluidized beds and related equipment on the basis of first principles. Such a capability could eliminate the need for expensive, time-consuming predesign testing.

Bellan, Josette↗
NASA NTRS2021

A Model for the Coupled Eruption of a Pseudostreamer and Helmet Streamer

A highly important aspect of solar activity is the coupling between eruptions and the surrounding coronal magnetic-field topology, which determines the trajectory and morphology of the event and can even lead to sympathetic eruptions from multiple sources. In this paper, we report on a numerica lsimulation of a new type of coupled eruption, in which a coronal jet initiated by a large pseudostreamer filament eruption triggers a streamer-blowout coronal mass ejection (CME) from the neighboring helmet streamer. Our configuration has a large opposite-polarity region positioned between the polarcoronal hole and a small equatorial coronal hole, forming a pseudostreamer flanked by the coronalholes and the helmet streamer. Further out, the pseudostreamer stalk takes the shape of an extendedarc in the heliosphere. We energize the system by applying photospheric shear along a section of the polarity inversion line within the pseudostreamer. The resulting sheared-arcade filament channel develops a flux rope that eventually erupts as a classic coronal-hole-type jet. However, the enhanced breakout reconnection above the channel as the jet is launched progresses into the neighboring helmet streamer, partially launching the jet along closed helmet streamer field lines and blowing out thestreamer top to produce a classic bubble-like CME. This CME is strongly deflected from the jet’s initial trajectory and contains a mixture of open and closed magnetic field lines. We present the detailed dynamics of this new type of coupled eruption, its underlying mechanisms and the implications of thiswork for the interpretation of in-situ and remote-sensing observations.

P F Wyper↗
DOE OSTI2022

Interface Capturing Simulations on Pool Boiling Performance With Multiple Nucleation Sites

High-fidelity pool boiling simulations can provide a numerical database for improving mechanistic boiling models by allowing for specific evaluation of interactions among bubbles. Previously published pool boiling simulations investigated two nucleation sites in which bubble growth at one site suppressed nucleation at the other site. Based on previous study results, more complicated interface-capturing simulations on pool boiling were conducted using PHASTA code with locally refined unstructured mesh.First, different boundary conditions (BCs) were assessed to support robustness and reproducibility of the boiling model. Then, a scale study was conducted at a larger domain with nine nucleation sites, four of which were activated. Involving more nucleation sites increased the complexity of bubble interactions from surrounding sites. Finally, bubble departure behavior influenced by wall heat flux was investigated. When heat flux was increased, the order of bubble departure changed, but diagonal bubbles always departed after one another. The departure time interval between the first and second bubble reduced as heat flux increased. The corresponding frequency was almost linearly proportional to the heat flux. In addition, bubble departure behavior was found to be greatly influenced by the nucleation site pattern. Multiple nucleation sites resulted in superimposed inhibitive effects from surrounding sites to each bubble, which extensively delayed the departure.The work presented here provides new insight on the fundamental understanding of boiling phenomena, contributes to the development of a 3D multiphase computational fluid dynamics (M-CFD) model, and provides a more comprehensive database for data-driven pool boiling studies.

Fan, Joy↗
DOE OSTI2020

Chapter 21. Atomic-Level Dislocation Dynamics in Irradiated Metals

Primary damage and microstructure evolution in structural nuclear materials operating under conditions of a high flux of energetic atomic particles and high temperature and stress lead to the formation of a high concentration, non-homogeneous distribution of defect clusters in the form of dislocation loops, voids, gas-filled bubbles and radiation-induced precipitates of nanometer scale. They cause changes in many material properties. Being obstacles to dislocation glide, they strongly affect the mechanical properties, with an increase in yield and flow stresses and a reduction in ductility. Atomic-scale computer simulations can provide details of how these effects are influenced by the obstacle structure, applied stress, strain rate and temperature. Processes such as obstacle cutting, transformation, absorption and drag are observed. Some recent results for body-centered and face-centered cubic metals are described in this review and, where appropriate, comparisons are drawn with predictions based on the elasticity theory of crystal defects. Perspectives on how to use this information at higher scales and in particular in mesoscale, dislocation dynamics simulations are also discussed.

Osetskiy, Yury N.↗
DOE OSTI2024

The Onset of Nucleate Boiling Prediction in a Mini Channel with the Eulerian Multiphase Flow

The prediction of boiling incipience is a critical issue for a reactor. Vapor lock causes operating instability and undesirable temperature rise, especially for small channels. This study uses the Eulerian multiphase flow boiling model to predict the onset of nucleate boiling in a millimeter scale rectangular channel. Two wall boiling models, namely the RPI (Rensselaer Polytechnic Institute) and the non-equilibrium sub cooled model are applied. Computational fluid dynamics is employed for this study to investigate the impact of operating conditions like liquid flow rate, operating pressure, heat flux, velocity profile, turbulence model, and inlet liquid temperature. Here, the effect of fundamental boiling parameters on the surface wall temperature and vapor volume fraction is also studied. The parameters under consideration are the bubble departure diameter, bubble departure frequency, nucleate site density, quenching time period, and interface heat transfer coefficients. A series of CFD calculations is conducted by varying the considered variables systematically in a wide range of flow conditions covering laminar, transition, and turbulent flows. The sub-component heat fluxes like the liquid convective, the vapor convective, the evaporative, and the quenching terms are monitored to infer the boiling dynamics in the transition region from the single-phase flow to the nucleate boiling zone. It is found that the CFD approach to detecting the boiling incipience point agrees reasonably with available experimental data. However, its limitations, like inaccuracy for lower flow rates and early transition to boiling flow, are also noticed due to the nature of the considered RPI model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

risingMicrobubbleLattice

This code is a simulation case to be run with OpenFOAM, an open-source computational fluid dynamics software. A gmsh mesh file is also included. Specifically, this simulation demonstrates the transport of a single microbubble rising through an ordered lattice due to an applied flow. The bubble deforms as it squeezes through the pores of the lattice. Is

Guo, Jack [Lawrence Livermore National Laboratory ↗
DOE OSTI2020

FICH: Feature instability cavitation history [Slides]

We seek to improve our ability to predict and understand the Richtmyer-Meshkov (RM) unstable dynamics of explosively shocked isolated features interfaced with vacuum or gases, and an HE products cavity. Explosively shocked features form as 2D or 3D RM phenomena. Liquid cavitation cause cavitation bubbles when shockwaves collide, and cavitation waves from wave interactions with isolated features form that source mass into the RM unstable sheets or spikes, and the cavitation waves may slice through bulk material into the HE products cavity.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

SNICAR-ADv4: a physically based radiative transfer model to represent the spectral albedo of glacier ice

Abstract. Accurate modeling of cryospheric surface albedo is essential for our understanding of climate change as snow and ice surfaces regulate the global radiative budget and sea-level through their albedo and mass balance. Although significant progress has been made using physical principles to represent the dynamic albedo of snow, models of glacier ice albedo tend to be heavily parameterized and not explicitly connected with physical properties that govern albedo, such as the number and size of air bubbles, specific surface area (SSA), presence of abiotic and biotic light absorbing constituents (LACs), and characteristics of any overlying snow. Here, we introduce SNICAR-ADv4, an extension of the multi-layer two-stream delta-Eddington radiative transfer model with the adding–doubling solver that has been previously applied to represent snow and sea-ice spectral albedo. SNICAR-ADv4 treats spectrally resolved Fresnel reflectance and transmittance between overlying snow and higher-density glacier ice, scattering by air bubbles of varying sizes, and numerous types of LACs. SNICAR-ADv4 simulates a wide range of clean snow and ice broadband albedo (BBA), ranging from 0.88 for (30 µm) fine-grain snow to 0.03 for bare and bubble-free ice under direct light. Our results indicate that representing ice with a density of 650 kg m−3 as snow with no refractive Fresnel layer, as done previously, generally overestimates the BBA by an average of 0.058. However, because most naturally occurring ice surfaces are roughened “white ice”, we recommend modeling a thin snow layer over bare ice simulations. We find optimal agreement with measurements by representing cryospheric media with densities less than 650 kg m−3 as snow and larger-density media as bubbly ice with a Fresnel layer. SNICAR-ADv4 also simulates the non-linear albedo impacts from LACs with changing ice SSA, with peak impact per unit mass of LACs near SSAs of 0.1–0.01 m2 kg−1. For bare, bubble-free ice, LACs actually increase the albedo. SNICAR-ADv4 represents smooth transitions between snow, firn, and ice surfaces and accurately reproduces measured spectral albedos of a variety of glacier surfaces. This work paves the way for adapting SNICAR-ADv4 to be used in land ice model components of Earth system models.

54 ENVIRONMENTAL SCIENCES↗