Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid 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.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. I. Dynamics

We study the effect of bubble breakup and coalescence on fine-scale dynamics of bubbly turbulent flows using direct numerical simulations. We perform two different simulations of dilute bubbly flows of void fraction 0.5%: one with bubbles breaking up and coalescing and the other without these physical processes. The volume of the fluid method is used for simulating bubbles undergoing breakup and coalescence while the bubbles are treated as rigid spheres in the immersed boundary method simulation. The energy spectrum in both types of simulation, consistent with previous studies, exhibits a -3 slope. We follow a single infinitesimal fluid element as it evolves to understand velocity gradient dynamics using conditional mean trajectories. We note finite-time divergence when the fluid element evolves under the action of inertial and pressure dynamics. The inertial, pressure, and viscous velocity gradient dynamics, when considered individually, produce the same results for bubble-induced turbulence (BIT) as with the classical homogeneous isotropic turbulence (HIT). Yet when the overall velocity gradient dynamics is considered, BIT results in non-cyclic trajectories moving toward stable node and unstable saddle while classical HIT shows cyclic behavior in their trajectories that move toward the origin. Interestingly, both the volume of fluid and immersed boundary simulations produce similar results. Therefore, there are two main takeaways from this research. First, new velocity gradient models are needed for BIT as their velocity gradients behavior is entirely different from the HIT. Second, we can neglect the bubble topology, breakup, and coalescence while studying or modeling the fine-scale dynamics of BIT.

42 ENGINEERING↗
DOE OSTI2026

Two-phase flow numerical analysis of electrode geometry for alkaline water electrolyzers

Hydrogen is a promising component of a future energy-secure and efficient economy, but its competitiveness depends on reducing production costs. One strategy is to operate alkaline water electrolyzers at higher current densities to increase output. However, this intensifies performance losses due to gas bubble accumulation, which blocks transport pathways and deactivates electrochemically active surfaces. Enhancing bubble evacuation through electrode design is therefore essential. Previous studies have explored various approaches — such as modifying surface morphology, applying sonication or pressure modulation, and introducing surfactants — but these efforts have addressed a limited range of conditions due to the complexity of two-phase flow and electrode geometries. Experiments have also largely been focused on either cell level improvements, which lack the information necessary to isolate each contributing factor, or on modified geometries that are not relevant to practical cell operation. From a modeling perspective, conventional Eulerian multiphase models do not track the complex gas–liquid interfacial dynamics and often neglect surface tension and contact angle effects, reducing their predictive accuracy. To provide insights on the effects of different electrode geometries on the performance of alklaine water electrolyzers this work employs an immersed boundary volume-of-fluid method to simulate bubble behavior in 3D porous electrodes. Multiple base electrode geometries, typically used in practice, with varying porosity are evaluated under a constant surface gas generation rate. Simulation data is analyzed to quantify electrode gas coverage, bubble size dynamics and other relevant metrics. Results show that porosity strongly influences bubble accumulation on electrode surfaces, with higher porosity reducing gas coverage, and its not strictly dependent on the electrode geometry. However, the electrode’s base geometry significantly affects gas accumulation at the separator gap, independent of porosity. A foam electrode geometry resulted in the lowest gas coverage of all electrodes with a median volumetric gas coverage of 11%, but at the cost of a 70% reduction in active area compared with the largest surface area electrode, while gyroid electrodes showed the best trade-off between gas coverage, particularly at the separator surface, and electrochemically active area. In conclusion, the results highlight the need for holistic electrode design strategies.

Alkaline↗
DOE OSTI2022

Vorticity-induced flow-focusing leads to bubble entrainment in an inkjet printhead: Synchrotron x-ray and volume-of-fluid visualizations

The oscillatory flows present in an inkjet printhead can lead to strong deformations of the air-liquid interface at the nozzle exit. Such deformations may lead to an inward directed air jet with bubble pinch-off and the subsequent entrainment of an air bubble, which is highly detrimental to the stability of inkjet printing. Understanding the mechanisms of bubble entrainment is therefore crucial to improving print stability. Here, in the present work, we use ultrafast x-ray phase-contrast imaging and direct numerical simulations based on the volume-of-fluid method to study the mechanisms underlying the bubble entrainment in a piezoacoustic printhead. We first demonstrate good agreement between experiments and numerics. We then show the different classes of bubble pinch-off obtained in experiments, and that those were also captured numerically. The numerical results are then used to show that the baroclinic torque, which is generated at the gas-liquid interface due to the misalignment of density and pressure gradients, results in a flow-focusing effect that drives the formation of the air jet from which a bubble can pinch off.

42 ENGINEERING↗
DOE OSTI2020

Shapes of rotating normal fluid He-3 versus superfluid He-4 droplets in molecular beams

Previous single-pulse extreme ultraviolet and X-ray coherent diffraction studies revealed that superfluid 4He droplets obtained in free jet expansion acquire sizable angular momentum, resulting in significant centrifugal distortion. Similar experiments with normal fluid 3He droplets may help elucidating the origin of the of the large degree of rotational excitation and highlight similarities and differences of dynamics in normal and superfluid droplets. Here, we present the first comparison of the shapes of isolated 3He and 4He droplets following expansion of the corresponding fluids in vacuum at temperatures as low as ~ 2 K. Large 3He and 4He droplets with average radii of ~160 nm and ~350 nm, respectively, were produced. We find that the majority of the 3He droplets in the beam correspond to rotating oblate spheroids with reduced average angular momentum ($\Lambda$) and reduced angular velocities ($\Omega$) similar to that of 4He droplets. Given the different physical nature of 3He and 4He, this similarity in $\Lambda$ and $\Omega$ may be surprising and suggest that similar mechanisms induce rotation regardless of the isotope. We hypothesized that the observed distribution of droplet sizes and angular momenta stem from processes in the dense region close to the nozzle. In this region, the significant velocity spread and collisions between the droplets induce excessive rotation followed by droplet fission. The process may repeat itself several times before the droplets enter the collision-fee high vacuum region further downstream.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

DEMONSTRATION OF A DATA-DRIVEN PHYSICS-BASED APPROACH FOR COMPUTATIONALLY EFFICIENT CFD PREDICTION OF TWO-PHASE BUBBLY FLOW

To realize efficient computational fluid dynamics (CFD) prediction of two-phase flow, a multi-scale physics-guided data-driven approach, Feature Similarity Measurement (FSM) technique was developed for error estimation in two-phase flow simulation using coarse-mesh CFD, to achieve a comparable accuracy as fine-mesh simulations with fast-running feature. By defining physics-guided parameters and variable gradients as physical features, FSM has the capability to capture the underlying local patterns in the coarse-mesh CFD simulation. Massive low-fidelity data and respective high-fidelity data are used to explore the underlying information relevant to the main simulation errors and the effects of phenomenological scaling. By learning from previous simulation data, a surrogate model using deep feedforward neural network (DFNN) can be developed and trained to estimate the simulation error of coarse-mesh CFD. In a demonstration case of two-phase bubbly flow, the DFNN model well captured and corrected the unphysical “peaks” in the velocity and void fraction profiles near the wall in the coarse-mesh configuration, even for extrapolative predictions. The research documented supports the feasibility of the physics-guided deep learning methods for coarse mesh CFD simulations which has a potential for the efficient industrial design.

42 ENGINEERING↗
DOE OSTI2021

Uncertainty quantification for Multiphase-CFD simulations of bubbly flows: a machine learning-based Bayesian approach supported by high-resolution experiments

In this paper, we developed a machine learning-based Bayesian approach to inversely quantify and reduce the uncertainties of multiphase computational fluid dynamics (MCFD) simulations for bubbly flows. The proposed approach is supported by high-resolution two-phase flow measurements, including those by double-sensor conductivity probes, high-speed imaging, and particle image velocimetry. Local distributions of key physical quantities of interest (QoIs), including the void fraction and phasic velocities, are obtained to support the Bayesian inference. In the process, the epistemic uncertainties of the closure relations are inversely quantified while the aleatory uncertainties from stochastic fluctuations of the system are evaluated based on experimental uncertainty analysis. The combined uncertainties are then propagated through the MCFD solver to obtain uncertainties of the QoIs, based on which probability-boxes are constructed for validation. The proposed approach relies on three machine learning methods: feedforward neural networks and principal component analysis for surrogate modeling, and Gaussian processes for model form uncertainty modeling. The whole process is implemented within the framework of an open-source deep learning library PyTorch with graphics processing unit (GPU) acceleration, thus ensuring the efficiency of the computation. The results demonstrate that with the support of high-resolution data, the uncertainties of MCFD simulations can be significantly reduced. The proposed approach has the potential for other applications that involve numerical models with empirical parameters.

42 ENGINEERING↗
DOE OSTI2020

Probing interfacial momentum closures in two-phase bubbly flow with machine learning-aided methods

Computational fluid dynamics (CFD) approach has already reached a high level of maturity for single-phase flows, however the development of closure models for two-phase flow requires additional attention. Multiphase CFD (M-CFD) methods resolve the conservation equations for mass, momentum and energy while differing in the approaches and strategies adopted in the physical closure models. The most widely adopted framework for M-CFD is the Eulerian-Eulerian two-fluid approach which assumes that all phases are co-existing inside each computational cell. For each fluid, the full set of conservation equations is solved; therefore, each fluid has a different velocity field. For adiabatic two-phase flow, the mechanisms of the interfacial momentum transfer are modeled by the interfacial forces representing different physical mechanisms. One of the crucial issues in the development and application of two-fluid model is the understanding of the interfacial momentum closures which determines the bubble distribution and migration behaviors. Dedicated experiments are performed to support the physical understanding and drive the closures’ development. However, limitations exist due to the uncertainties in the experimental measurement and the simplified analytical assumptions which have difficulties on representing the complex non-linear flow fields. In this paper, a data-driven approach, Feature Similarity Measurement (FSM), is developed and proposed to resolve the challenges of modeling the interfacial forces closures. Case study is performed with two-phase flow scenarios where the high-fidelity experimental data is available. Within the Eulerian-Eulerian two-fluid framework, only momentum equations for gas and liquid phases are solved and reduced-order interfacial momentum closures are aided with FSM. Predictions of void fraction and velocity fields are analyzed and demonstrate the potential of machine learning-driven interfacial forces closures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Motion of Air Bubbles in a Cement Slurry

The dynamics of air (gas) bubbles in a column of cement slurry is examined numerically. The air injected at the bottom of a laboratory-scale column through a porous distributor plate spatially distributes and migrates as a swarm of bubbles throughout the slurry toward the freeboard. The two-phase system of the cement slurry and the air bubbles is modeled using the conservation equations of mass and linear momentum in the framework of the volume-of-fluid (VOF) approach. The cement slurry is modeled using the Herschel–Bulkley and Bingham fluid models. Results show that the mean Sauter diameter and the mean rise velocity of the bubbles decrease with the gas flow rate. Meanwhile, it is found that the rising of the bubbles is controlled by breakup events, along with relatively weak path instabilities of the bubbles resulting in relatively straight trajectories, independent of the gas flow rate. The extent of the yielded region appears larger for the Herschel–Bulkley model compared to the Bingham fluid model (by approximately 10%).

36 MATERIALS SCIENCE↗
DOE OSTI2020

Performance of Closed Loop Venturi Aspirated Aeration System: Experimental Study and Numerical Analysis with Discrete Bubble Model

In wastewater treatment plants, aeration plays a significant role as it encourages aerobic respiration of microbes, which are necessary to break down carbonaceous matter in the waste stream. This process can account for the majority of energy use in wastewater treatment plants. The aeration process is also necessary in odor control in lagoons and in the aquaculture industry. Generally, the aeration process is accomplished with compressors or blowers which can be of low efficiency due to ideal gas laws. This study introduces the idea of increasing aeration efficiency by looping water from a reservoir through a piping network which includes a venturi aspirator at its inlet. For this purpose, an experimental study has been conducted in a laboratory setup with a pump which pulls water from a tank, passes it through a Venturi aspirator and a helical piping network intended to increase bubble residence time, before depositing it back into the bulk fluid tank. This same process is modeled computationally, using a discrete bubble method (DBM), with good agreement with experiments. The overall purpose here is to determine the optimal configurations for standard aeration efficiency (SAE) and the standard oxygen transfer rate (SOTR). A parametric study has been implemented using the DBM based on different hydraulic and flow parameters. The model is also used to predict the SAE of a hypothetical aeration system. Results indicate that it is possible to achieve SAE values in the range of surface aerators or submersed jet aerators using the proposed aeration system with less complex components, thereby decreasing overall costs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Benchmark of interfacial area concentration approaches for the two-fluid model in gas-dispersed condensing flow

Accurate predictions of interfacial area concentration (IAC) are required for closure of the condensation phase change term. The interfacial area transport equation (IATE) can dynamically predict interfacial area concentration and is increasingly being used in the benchmarking of phase change flows over static based correlations which require a flow regime map. For higher flow regimes, the two-group two-fluid IATE model can separately model the transport of spherical bubbles and cap/slug bubbles. Recently, a modified formulation of the one-group and two-group two-fluid model with interfacial area transport equation was proposed for gas-dispersed condensing flows. Three datasets have been used for validation consisting of condensing two-phase flow data in the unheated section of a vertical annulus with conditions spanning a range of pressures, inlet subcoolings and mass fluxes for benchmarking the various IAC approaches. For bubbly flows, the coupled void transport - IATE simulation results in good agreement in the predictions of the void fraction and the interfacial area concentration under nearly all sixty conditions. Considering the proposed two-group two-fluid IATE model formulation, the group-2 Nusselt number correlation is investigated to account for the appropriate heat transfer length scale for group-2 bubbles. Validation of the proposed two-group two-fluid model against the new dataset shows good agreement with experimental data with significant improvement over the one-group IATE model for conditions with substantial group-2 void fraction throughout the simulation domain.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A Numerical Analysis of The Pool Boiling Process in Flooded Evaporators Using Different Refrigerants

The pool boiling process can be observed in several energy conversion processes including commercial and industrial refrigeration, industrial air-cooling operations, and power generation. The process becomes more involved when pool boiling in a tube bundle is considered. In the current study, a numerical model is developed to predict the key performance parameters of a flooded evaporator while considering a range of working fluids. A kettle reboiler configuration was considered, and a performance model was developed to account for boiling on individual tubes, merging of vapor bubbles, and movement under gravity. A volume of fluid (VOF) model was used to deal with the different phases in the simulation. Various fluids were considered in this study including HFE-7000, HFE-7300, and water. The trajectories of the bubbles were tracked, and the resulting information has been summarized in physical measurable quantities.

Zhang, Mingkan↗
DOE OSTI2023

Heat transfer enhancement due to cold cap motion from bubbling in a waste glass melter

In this study, a computational fluid dynamics (CFD) model was developed to model the motion of a solid cold cap in a waste glass melter. Forced convection bubblers at the base of the melter release air into the molten glass, which forms large bubbles that travel upwards to the cold cap and augment heat transfer from the glass to the cold cap. The CFD model employs the Navier- Stokes equations to solve for the fluctuating flowfield using a rigid body motion dynamic fluid body interaction module. This allows for movement of the floating body in response to the bubbling forces calculated at each time step. Further, the heat flux delivered to the cold cap by the convective bubbling is studied as a function of the bubbling rate. Results for the moving cold cap are compared with the computed heat flux trends for a stationary cold cap. The heat flux delivered to the cold cap from the glass is 25% higher for the case with the moving cold cap. The heat flux was found to be proportional to v0.6 as opposed to v0.9 (where v is the normalized bubbling rate) for the stationary cold cap.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Machine learning accelerated turbulence modeling of transient flashing jets

We report modeling the sudden depressurization of superheated liquids through nozzles is a challenge because the pressure drop causes rapid flash boiling of the liquid. The resulting jet usually demonstrates a wide range of structures, including ligaments and droplets, due to both mechanical and thermodynamic effects. As the simulation comprises increasingly numerous phenomena, the computational cost begins to increase. One way to moderate the additional cost is to use machine learning surrogacy for specific elements of the calculation. This study presents a machine learning-assisted computational fluid dynamics approach for simulating the atomization of flashing liquids accounting for distinct stages, from primary atomization to secondary breakup to small droplets using the Σ - Y model coupled with the homogeneous relaxation model. Notably, the models for thermodynamic non-equilibrium (HRM) and Σ - Y are coupled, for the first time, with a deep neural network that simulates the turbulence quantities, which are then used in the prediction of superheated liquid jet atomization. The data-driven component of this method is used for turbulence modeling, avoiding the solution of the two-equation turbulence model typically used for Reynolds-averaged Navier-Stokes simulations for these problems. Both the accuracy and speed of the hybrid approach are evaluated, demonstrating adequate accuracy and at least 25% faster computational fluid dynamics simulations than the traditional approach. This acceleration suggests that perhaps additional components of the calculation could be replaced for even further benefit. Published under an exclusive license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Development of machine learning framework for interface force closures based on bubble tracking data

Interfacial force closures in the two-fluid model play a critical role for the predictive capabilities of void fraction distribution. However, the practices of interfacial force modeling have long been challenged by the inherent physical complexity of the two-phase flows. The rapidly expanding computational capabilities in the recent years have made high-fidelity data from the interface-captured direct numerical simulation become more available, and hence potential for data-driven interfacial force modeling has prevailed. In this work, we established a data-driven modeling framework integrated to the HZDR multiphase Eulerian-Eulerian framework for computational fluid dynamics simulations. The data-driven framework is verified in a benchmark problem, where a feedforward neural network managed to capture the non-linear mapping between bubble Reynolds number and drag coefficient and reproduce the void distribution resulting from the baseline model in the test case. The second focus is on utilizing the bubble tracking data set to form a closure for the bubble drag in the turbulent bubbly flow, in which the drag coefficient is set to be correlated with the bubble Reynolds number and the Eötvös number. Pseudo-steady state filtering in the Frenet Frame was carried out to obtain the drag coefficient from the turbulent bubbly flow data. The performance of the data-driven drag model is also examined through a case study, where improvement of model’s prediction near-wall is regarded necessary. In conclusion, discussion and further plans of investigation are provided.

42 ENGINEERING↗
DOE OSTI2020

A Scintillating Xenon Bubble Chamber for Dark Matter Detection. Final Report

This report describes the progress in the search for particle dark matter achieved under DOE award DE-SC0012161, as well as the invention of a new dark matter and neutrino detection technique. Work supported by this award follows in three distinct thrusts: (1) the successful completion of the Generation-1 Direct Detection experiment PICO-60, which set the world-leading limit on the spin-dependent coupling of dark matter to protons and achieved a precise understanding of the response of bubble chambers to nuclear-recoil signals and electron-recoil backgrounds; (2) the construction of the Generation-2 Direct Detection experiment LZ, which will soon be the world's most sensitive dark matter detector; and (3) the demonstration of the first scintillating bubble chamber. The development of a scintillating liquid noble bubble chamber was the primary goal of the proposed work, with the aim of combining the excellent background rejection of a PICO-style bubble chamber with the event-by-event energy resolution of a liquid-noble detector such as LZ. This work produced the first ever observation of coincident scintillation and bubble nucleation by a nuclear recoil in a superheated fluid and also revealed an unexpected benefit to the use of noble liquids in a bubble chamber, namely the ability to increase the degree of superheat by an order of magnitude beyond that achievable in PICO-style detectors while maintaining PICO's world-leading background rejection. This discovery has led to new projects in the US and Canada, developing noble liquid bubble chambers both as detectors for low-mass dark matter and as detectors of coherent elastic neutrino-nucleus scattering (CEvNS) by neutrinos produced at nuclear reactors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

System and method of mapping topology

A system for determining a fingerprint of a structure is provided. The system includes a plurality of granules inserted in a structure having a plurality of fissures, fractures, and cracks (collectively apertures), each granule comprising a membrane, and at least one bubble of compressed gas formed in the membrane, the membrane selectively dissolving in presence of a predetermined fluid and thereby selectively bursting the at least one bubble, thereby generating a concussing vibration, at least at least i) three detection devices for two-dimensional mapping or ii) four detection devices for three-dimensional mapping placed proximate to the structure according to a predetermined placement schedule, and a computing device comprising a processor configured to receive data from the at least three or four detection devices and to determine location of the at least one bubble of each of the plurality of the granules at the time of bursting by triangulating the concussive vibration in order to determine location of the at least one bubble.

Pyrak-Nolte, Laura↗
DOE OSTI2021

The benefits of diligence: how precise are predicted gravitational wave spectra in models with phase transitions?

Models of particle physics that feature phase transitions typically provide predictions for stochastic gravitational wave signals at future detectors and such predictions are used to delineate portions of the model parameter space that can be constrained. The question is: how precise are such predictions? Uncertainties enter in the calculation of the macroscopic thermal parameters and the dynamics of the phase transition itself. We calculate such uncertainties with increasing levels of sophistication in treating the phase transition dynamics. Currently, the highest level of diligence corresponds to careful treatments of the source lifetime; mean bubble separation; going beyond the bag model approximation in solving the hydrodynamics equations and explicitly calculating the fraction of energy in the fluid from these equations rather than using a fit; and including fits for the energy lost to vorticity modes and reheating effects. The lowest level of diligence incorporates none of these effects. We compute the percolation and nucleation temperatures, the mean bubble separation, the fluid velocity, and ultimately the gravitational wave spectrum corresponding to the level of highest diligence for three explicit examples: SMEFT, a dark sector Higgs model, and the real singlet-extended Standard Model (xSM). In each model, we contrast different levels of diligence in the calculation and find that the difference in the final predicted signal can be several orders of magnitude. Our results indicate that calculating the gravitational wave spectrum for particle physics models and deducing precise constraints on the parameter space of such models continues to remain very much a work in progress and warrants care.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Evaluation of the Barracuda Software Package for Simulating Bubble Motion in Vibrating Liquid-Filled Containers

The commercial software package Barracuda, developed by CPFD Software for simulating particle-laden fluid flows, is evaluated as a means to simulate the motion of bubbles in vibrating liquid-filled containers. Demonstration simulations of bubbles rising due to buoyancy forces in a cylinder filled with silicone oil and angled at 0, 30, 45, and 60 degrees from the vertical were performed by CPFD Software. The results of these simulations are discussed, and the capabilities of Barracuda for simulating bubble motion are assessed. It was determined that at present Barracuda does not meet the needs of the desired application. Further developments that would enable its use for this application are highlighted.

42 ENGINEERING↗