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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Towards Probablistic Assessment of Hypobaric Decompression Sickness Treatment

INTRODUCTION: Pressure, oxygen (O2), and time are the pillars to effective treatment of decompression sickness (DCS). The NASA DCS Treatment Model links a decrease in computed bubble volume to the resolution of a symptom. The decrease in volume is realized in two stages: a) during the Boyle's Law compression and b) during subsequent dissolution of the gas phase by the O2 window. METHODS: The cumulative distribution of 154 symptoms that resolved during repressurization was described with a log-logistic density function of pressure difference (deltaP as psid) associated with symptom resolution and two other explanatory variables. The 154 symptoms originated from 119 cases of DCS during 969 exposures in 47 different altitude tests. RESULTS: The probability of symptom resolution [P(symptom resolution)] = 1 / (1+exp(- (ln(deltaP) - 1.682 + 1.089×AMB - 0.00395×SYMPTOM TIME) / 0.633)), where AMB is 1 when the subject ambulated as part of the altitude exposure or else 0 and SYMPTOM TIME is the elapsed time in min from start of the altitude exposure to recognition of a DCS symptom. The P(symptom resolution) was estimated from computed deltaP from the Tissue Bubble Dynamics Model based on the "effective" Boyle's Law change: P2 - P1 (deltaP, psid) = P1×V1/V2 - P1, where V1 is the computed volume of a spherical bubble in a unit volume of tissue at low pressure P1 and V2 is computed volume after a change to a higher pressure P2. V2 continues to decrease through time at P2, at a faster rate if 100% ground level O2 was breathed. The computed deltaP is the effective treatment pressure at any point in time as if the entire deltaP was just from Boyle's Law compression. DISCUSSION: Given the low probability of DCS during extravehicular activity and the prompt treatment of a symptom with options through the model it is likely that the symptom and gas phase will resolve with minimum resources and minimal impact on astronaut health, safety, and productivity.

Conkin, J.↗
NASA NTRS1980

Dissolution of multicomponent bubbles

The behavior of an isolated, stationary, multicomponent gas bubble in a glassmelt containing several dissolved gases is considered. The relevant mass-transport equations are formulated and calculations are performed for the case of two diffusing gases using a quasi-stationary model and a numerical solution of the exact mass-transfer equations. The results obtained from these two approaches are compared. The factors which govern the dissolution or growth of a bubble are thermodynamic and kinetic in origin. The tendency of a bubble to grow or shrink at long times is controlled by departure from overall equilibrium, whereas the short-time bubble dynamics may be dominated by kinetic effects. As a result of the existence of these dual influences, maxima and/or minima occur in the functional dependence of the bubble radius on time.

Weinberg, M. C.↗
DOE OSTI2023

Gas-Liquid Flow Modeling for Renewable Fuels Production

Aerobic/anaerobic and gas fermentation pathways have emerged as promising new technologies for the generation of renewable fuels/chemicals from biomass derived sugars, and mixtures of greenhouse/energy rich gas streams (CO2/CH4/H2/CO) via microbial action. Example pathways include sugars-to-ethanol conversion, biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the scale-up, design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of at-scale (~500 m3) bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for gas uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM [1], to perform these simulations, which has been validated against small-scale reactors in our previous work [2]. There is currently a knowledge-gap regarding bubble-size distributions when using gas mixtures with vastly different properties, which can have a significant impact overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a population balance model (PBM) for bubble-size-distribution modeling that is validated against small-scale experiments in our solver with an uncertainty quantification study for bubble coalescence and break-up model parameters. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient are analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors. This work will also examine the performance of different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor mass-transfer coefficient, gas hold-up, and dissolved gas distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented. [1] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998. [2] Rahimi, M., Sitaraman, H., Humbird, D. and Stickel, J., Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution, Chemical Engineering Research and Design, 139: 283-295.

BIOMASS FUELS↗
NASA NTRS1988

Elasticity effects on cavitation in a squeeze film damper undergoing noncentered circular whirl

Elasticity of the liner and its effects on cavitation were numerically determined for a squeeze film damper subjected to dynamic loading. The loading was manifested as a prescribed motion of the rotor undergoing noncentered circular whirl. The boundary conditions were implemented using Elrod's algorithm which conserves lineal mass flux through the moving cavitation bubble as well as the oil film region of the damper. Computational movies were used to analyze the rapidly changing pressures and vapor bubble dynamics throughout the dynamic cycle for various flexibilities in the damper liner. The effects of liner elasticity on cavitation were only noticeable for the intermediate and high values of viscosity used in this study.

Brewe, David E.↗
NASA NTRS2002

Length Scale and Gravity Effects on Microgravity Boiling Heat Transfer

Boiling is a complex phenomenon where hydrodynamics, heat transfer, mass transfer, and interfacial phenomena are tightly interwoven. An understanding of boiling and critical heat flux in microgravity environments is of importance to space based hardware and processes such as heat exchange, cryogenic fuel storage and transportation, electronic cooling, and material processing due to the large amounts of heat that can be removed with relatively little increase in temperature. Although research in this area has been performed in the past four decades, the mechanisms by which heat is removed from surfaces in microgravity are still unclear. In earth gravity, buoyancy is an important parameter that affects boiling heat transfer through the rate at which bubbles are removed from the surface. A simple model describing the bubble departure size based on a quasistatic force balance between buoyancy and surface tension is given by the Fritz [I] relation: Bo(exp 1/2) = 0.0208 theta where Bo is the ratio between buoyancy and surface tension forces. For small, rapidly growing bubbles, inertia associated with the induced liquid motion can also cause bubble departure. In microgravity, the magnitude of effects related to natural convection and buoyancy are small and physical mechanisms normally masked by natural convection in earth gravity such as Marangoni convection can substantially influence the boiling and vapor bubble dynamics. CHF (critical heat transfer) is also substantially affected by microgravity. In 1 g environments, Bo has been used as a correlating parameter for CHF. Zuber's CHF model for an infinite horizontal surface assumes that vapor columns formed by the merger of bubbles become unstable due to a Helmholtz instability blocking the supply of liquid to the surface. The jets are spaced lambda(sub D) apart, where lambda(sub D) = 2pi square root of 3[(sigma)/(g(rho(sub l) - rho(sub v)](exp 1/2) = 2pi square root of 3 L Bo(exp -1/2) = square root of 3 lambda(sub c) and is the wavelength that amplifies most rapidly. The critical wavelength, lambda(sub c), is the wavelength below which a vapor layer underneath a liquid layer is stable. For heaters with Bo smaller than about 3 (heaters smaller than lambda(sub D)), the above model is not applicable, and surface tension effects dominate. Bubble coalescence is thought to be the mechanism for CHF under these conditions. Small Bo can result by decreasing the size of a heater in earth gravity, or by operating a large heater in a lower gravity environment. In the microgravity of space, even large heaters can have low Bo, and models based on Helmholtz instability should not be applicable. The macrolayer model of Haramura and Katto is dimensionally equivalent to Zuber's model and has the same dependence on gravity, so it should not be applicable as well. The goal of this work is to determine how boiling heat transfer mechanisms in a low-g environment are altered from those at higher gravity levels. Boiling data using a microheater array was obtained under gravity environments ranging from 1.8 g to 0.02 g with heater sizes ranging from 2.7 mm to 1 mm. The boiling behavior for 2.7 mm at 0.02 g looked quite similar to boiling on the 1 mm heater at 1 g-the formation of a large primary bubble surrounded by smaller satellite bubbles was observed under both conditions. The similarity suggests that for heaters smaller than some fraction of I(sub c), coalescence and surface tension dominate boiling heat transfer. It also suggests that microgravity boiling can be studied by studying boiling on very small heaters.

Kim, Jungho↗
NASA NTRS2003

Mathematical model of diffusion-limited evolution of multiple gas bubbles in tissue

Models of gas bubble dynamics employed in probabilistic analyses of decompression sickness incidence in man must be theoretically consistent and simple, if they are to yield useful results without requiring excessive computations. They are generally formulated in terms of ordinary differential equations that describe diffusion-limited gas exchange between a gas bubble and the extravascular tissue surrounding it. In our previous model (Ann. Biomed. Eng. 30: 232-246, 2002), we showed that with appropriate representation of sink pressures to account for gas loss or gain due to heterogeneous blood perfusion in the unstirred diffusion region around the bubble, diffusion-limited bubble growth in a tissue of finite volume can be simulated without postulating a boundary layer across which gas flux is discontinuous. However, interactions between two or more bubbles caused by competition for available gas cannot be considered in this model, because the diffusion region has a fixed volume with zero gas flux at its outer boundary. The present work extends the previous model to accommodate interactions among multiple bubbles by allowing the diffusion region volume of each bubble to vary during bubble evolution. For given decompression and tissue volume, bubble growth is sustained only if the bubble number density is below a certain maximum.

Evaluation Studies↗
DOE OSTI2024

Molecular-Resolution Electron Imaging of Defects and Dynamics at the Ice-Water Interface

Water crystallization into hexagonal ice (type I h ) is one of the most critical processes relevant to the Earth’s environment and human activities. However, despite recent breakthroughs in imaging non-equilibrium condensed ice structures, the ice-water interface has never been imaged at a molecular resolution. This is primarily due to the low stability of the hydrogen bonds in ice under high-resolution microscopy conditions and a lack of methods to prepare compatible samples. Here, this presentation describes the first molecular-resolution imaging of ice crystallized from liquid water and the ice-water interface using high-resolution transmission electron microscopy (HRTEM). By encapsulating deionized (DI) water between two amorphous carbon (a-C) TEM grids and subsequently freezing it with liquid N 2 on a cryo sample stage, we generated two types of ice: non-equilibrium, condensed ice from the atmosphere and encapsulated ice from the DI water (Fig. 1A). Condensed ice usually shows irregular, spherulitic shapes (Fig. 1B). Selected area electron diffraction (SAED) shows that they are a mixture of cubic and hexagonal crystals (Fig. 1C). On the contrary, encapsulated ice forms thin films that contain large-area single-crystalline regions of hexagonal ice oriented along the [0001] zone axis (Fig. 1D). Differential electron energy-loss spectroscopy (EELS) confirmed the high purity of the encapsulated ice free from organic contaminations that are common in other encapsulation methods for HRTEM such as graphene liquid cells. These single-crystalline areas are robust under the electron beam up to ~100 e/Å 2 s. Aberration-corrected HRTEM imaging in these areas achieved a line resolution of ~1.3 Å (Fig. 1E and F). This platform allows us to study near-equilibrium ice structures and dynamics at an unprecedented spatial resolution (Fig. 2). For example, we discovered subdomain-rich regions near the defective crystal edges despite the structure appearing single-crystalline according to diffraction criteria. These subdomains connect via low-angle grain boundaries with flat energy landscapes as a function of tilt angles (according to simulations), showing the high tolerance of ice to defect structures. When we tuned the sample temperature and electron flux rate, we observed radiolysis-controlled bubble generation and dissolution in ice single crystals near a steady state of bubble dynamics. Furthermore, rich beam-induced melting and recrystallization dynamics were observed at the ice-water interface with lattice resolution. These data represent the first observation of the ice-water phase transformation at the sub-nanometer level. In summary, the methods developed in this work enabled molecular-resolution observations of ice and the ice-water interface and shed light on the microstructures and phase transformation pathways. Finer control on the temperature, electron irradiation profile, and imaging detector could eventually lead to real-time observation of ice nucleation in water and address long-standing questions in the nucleation pathways.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Tracking Local pH Dynamics during Water Electrolysis via In-Line Continuous Flow Raman Spectroscopy

The performance of electrochemical devices, which play a critical role in decarbonization efforts, is often governed by proton-coupled electron transfer reactions at the electrode–electrolyte interface. These reactions are highly sensitive to the complex and dynamic microenvironment present at the electrode surface. However, characterizing this environment─particularly monitoring interfacial pH and its evolution under reaction conditions─remains challenging, necessitating the development of advanced analytical tools. Here, in this study, we introduce in-line continuous flow Raman spectroscopy (CFRS) as a spectroelectrochemical platform for quantifying interfacial pH swings generated during water-splitting. By monitoring phosphate ion speciation and controlling the hydrodynamics with a flow cell, we measure pH swings as a function of current density, flow rate, and distance from the electrode. Comparison with theoretical models reveals the impact of bulk pH, boundary layer thickness, and bubble dynamics at high current densities. Collectively, these findings establish CFRS as a platform for quantitatively investigating pH dynamics, offering critical insights for advancing electrochemical energy conversion technologies.

Marquez, Raul A. [Univ. of Texas, Austin, TX (Unit↗
NASA NTRS2007

Potential Fifty Percent Reduction in Saturation Diving Decompression Time Using a Combination of Intermittent Recompression and Exercise

Conventional saturation decompression protocols use linear decompression rates that become progressively slower at shallower depths, consistent with free gas phase control vs. dissolved gas elimination kinetics. If decompression is limited by control of free gas phase, linear decompression is an inefficient strategy. The NASA prebreathe reduction program demonstrated that exercise during O2 prebreathe resulted in a 50% reduction (2 h vs. 4 h) in the saturation decompression time from 14.7 to 4.3 psi and a significant reduction in decompression sickness (DCS: 0 vs. 23.7%). Combining exercise with intermittent recompression, which controls gas phase growth and eliminates supersaturation before exercising, may enable more efficient saturation decompression schedules. A tissue bubble dynamics model (TBDM) was used in conjunction with a NASA exercise prebreathe model (NEPM) that relates tissue inert gas exchange rate constants to exercise (ml O2/kg-min), to develop a schedule for decompression from helium saturation at 400 fsw. The models provide significant prediction (p < 0.001) and goodness of fit with 430 cases of DCS in 6437 laboratory dives for TBDM (p = 0.77) and with 22 cases of DCS in 159 altitude exposures for NEPM (p = 0.70). The models have also been used operationally in over 25,000 dives (TBDM) and 40 spacewalks (NEPM). The standard U.S. Navy (USN) linear saturation decompression schedule from saturation at 400 fsw required 114.5 h with a maximum Bubble Growth Index (BGI(sub max)) of 17.5. Decompression using intermittent recompression combined with two 10 min exercise periods (75% VO2 (sub peak)) per day required 54.25 h (BGI(sub max): 14.7). Combined intermittent recompression and exercise resulted in a theoretical 53% (2.5 day) reduction in decompression time and theoretically lower DCS risk compared to the standard USN decompression schedule. These results warrant future decompression trials to evaluate the efficacy of this approach.

Gernhardt, Michael I.↗
DOE OSTI2024

Modelling the ejection of primary aerosols during the fast pyrolysis of biomass anisotropic particles

A model for the fast pyrolysis of anisotropic biomass particles is presented which considers bubbling dynamics within the liquid intermediate phase (metaplast) and aerosol ejection from this phase. The model employs the population balance equation and the method of moments to estimate the production rate and resultant size distribution of aerosol ejections, incorporating a detailed CRECK reaction mechanism, and considers the effect of anisotropic biomass microstructure on the intraparticle transport of mass and energy. Here, this study investigates the impact of particle size, heating rate (heat transfer coefficient), and lignocellulosic composition on aerosol ejection. The model predicts that, at high heating rates (convective heat transfer coefficient of 359 W/m 2 .K), aerosols can contribute over 20% to the heavy fraction yield in bio-oil for small particles (1 mm diameter, 4 mm length). The model can predict aerosol size distribution and surface area, indicating an average size of 20 μm for bubbles and 5 μm for aerosols during increased bubble production and aerosol ejection rates. These findings are consistent with prior experimental results and provide essential information for future modeling of extra-particle reactions of the aerosols as they progress through the reactor.

09 BIOMASS FUELS↗
DOE OSTI2021

Favorable morphology and electronic conductivity of functional sublayers for highly efficient water splitting electrodes

Low electronic conductivities and improper morphologies of anode electrodes greatly limit the reaction area, catalyst utilization and efficiency in proton exchange membrane water electrolyzers. In this study, conductive sublayers with different conductivities and morphologies were introduced into anode electrodes in membrane-based water electrolyzers. In-situ and ex-situ investigation results showed that conductive sublayers (Au mesh and carbon nanotube (CNT) film) augmented the sheet conductivity of anode electrodes by up to 4000 times (from 2000 to 0.5 ohm square -1 ), and the ohmic resistance of water electrolyzers was reduced to 1/3 when inserting conductive sublayers. In addition, CNT film provided a higher electrochemical active area than Au mesh, because of favorable morphologies (large porosity and surface area) of CNT fibers on CNT films. Therefore, the current density of water splitting was increased by 3 times (from 4.55 mA cm -2 to 14.83 mA cm -2 ) at 2.5 V compared to a conventional anode electrode. Visualizations on bubble dynamics showed improved performances with conductive sublayers; this was mainly due to greatly increased number of reaction sites, highly spread reaction area (from 50 to 1000 um), and reduced activation overpotential. Therefore, a balance between high electronic conductivity and nanoporous morphology is essential to the anode electrode for larger reaction sites and areas in highly efficient water electrolyzers.

25 ENERGY STORAGE↗
NASA NTRS1974

Non-Newtonian effects on flow-generated cavitation and on cavitation in a pressure field

Observations are presented which show that the stresses in a flow field of very dilute polymer are not well enough described by the Navier-Stokes equations to accurately predict cavitation. The contitutive equation for the particular polymer and concentration used is needed. The second-order fluid form in which accelerations are relatively important appears capable of explaining observed cavitation suppression by changing the pressure field due to flow. Bubble dynamics in stationary dilute polymer solutions are also examined and found to be little different from those in water.

Ellis, A. T.↗
NASA NTRS1976

Development of a high capacity bubble domain memory element and related epitaxial garnet materials for application in spacecraft data recorders. Item 2: The optimization of material-device parameters for application in bubble domain memory elements for spacecraft data recorders

Bubble domain materials and devices are discussed. One of the materials development goals was a materials system suitable for operation of 16 micrometer period bubble domain devices at 150 kHz over the temperature range -10 C to +60 C. Several material compositions and hard bubble suppression techniques were characterized and the most promising candidates were evaluated in device structures. The technique of pulsed laser stroboscopic microscopy was used to characterize bubble dynamic properties and device performance at 150 kHz. Techniques for large area LPE film growth were developed as a separate task. Device studies included detector optimization, passive replicator design and test and on-chip bridge evaluation. As a technology demonstration an 8 chip memory cell was designed, tested and delivered. The memory elements used in the cell were 10 kilobit serial registers.

Besser, P. J.↗
NASA NTRS1980

Conceptual design of two-phase fluid mechanics and heat transfer facility for spacelab

Five specific experiments were analyzed to provide definition of experiments designed to evaluate two phase fluid behavior in low gravity. The conceptual design represents a fluid mechanics and heat transfer facility for a double rack in Spacelab. The five experiments are two phase flow patterns and pressure drop, flow boiling, liquid reorientation, and interface bubble dynamics. Hardware was sized, instrumentation and data recording requirements defined, and the five experiments were installed as an integrated experimental package. Applicable available hardware was selected in the experiment design and total experiment program costs were defined.

North, B. F.↗
NASA NTRS1983

Space Processing Applications Rocket (SPAR) project SPAR 7

The postflight reports of each of the Principal Investigators of three selected science payloads are presented in addition to the engineering report as documented by the Marshall Space Flight Center (MSFC). Pertinent portions of ground-based research leading to the ultimate selection of the flight sample composition are described including design, fabrication and testing. Containerless processing technology, containerless processing bubble dynamics, and comparative alloy solidification are the experiments discussed.

Poorman, R. M.↗
NASA NTRS1988

Mass transport phenomena between bubbles and dissolved gases in liquids under reduced gravity conditions

The experimental and analytical work that was done to establish justification and feasibility for a shuttle middeck experiment involving mass transfer between a gas bubble and a liquid is described. The experiment involves the observation and measurement of the dissolution of an isolated immobile gas bubble of specified size and composition in a thermostatted solvent liquid of known concentration in the reduced gravity environment of earth orbit. Methods to generate and deploy the bubble were successful both in normal gravity using mutually buoyant fluids and under reduced gravity conditions in the NASA Lear Jet. Initialization of the experiment with a bubble of a prescribed size and composition in a liquid of known concentration was accomplished using the concept of unstable equilibrium. Subsequent bubble dissolution or growth is obtained by a step increase or decrease in the liquid pressure. A numerical model was developed which simulates the bubble dynamics and can be used to determine molecular parameters by comparison with the experimental data. The primary objective of the experiment is the elimination of convective effects that occur in normal gravity.

Dewitt, Kenneth J.↗