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.
Computations of Boiling in Microgravity
The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work is limited to very simple models. In this project, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. The proposed work is based on previously funded NASA work that allowed us to develop a two-dimensional numerical method for boiling flows and to demonstrate the ability of the method to simulate film boiling. While numerical simulations of multi-fluid flows have been advanced in a major way during the last five years, or so, similar capability for flows with phase change are still in their infancy. Although the feasibility of the proposed approach has been demonstrated, it has yet to be extended and applied to fully three-dimensional simulations. Here, a fully three-dimensional, parallel, grid adaptive code will be developed. The numerical method will be used to study nucleate boiling in microgravity, with particular emphasis on two aspects of the problem: 1) Examination of the growth of bubbles at a wall nucleation site and the instabilities of rapidly growing bubbles. Particular emphasis will be put on accurately capturing the thin wall layer left behind as a bubble expands along a wall, on computing instabilities on bubble surfaces as bubbles grow, and on quantifying the effects of both these phenomena on heat transfer; and 2) Examination of the effect of shear flow on bubble growth and heat transfer.
Computations of Boiling in Microgravity
The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.
Theoretical and experimental comparison of vapor cavitation in dynamically loaded journal bearings
Vapor cavitation for a submerged journal bearing under dynamically loaded conditions was investigated. The observation of vapor cavitation in the laboratory was done by high-speed photography. It was found that vapor cavitation occurs when the tensile stress applied to the oil exceeded the tensile strength of the oil or the binding of the oil to the surface. The theoretical solution to the Reynolds equation is determined numerically using a moving boundary algorithm. This algorithm conserves mass throughout the computational domain including the region of cavitation and its boundaries. An alternating direction implicit (MDI) method is used to effect the time march. A rotor undergoing circular whirl was studied. Predicted cavitation behavior was analyzed by three-dimensional computer graphic movies. The formation, growth, and collapse of the bubble in response to the dynamic conditions is shown. For the same conditions of dynamic loading, the cavitation bubble was studied in the laboratory using high-speed photography.
The nuclear superbubble of NGC 3079
We have used the Hawaii Imaging Fabry-Perot Interferometer (HIFI) at the CFH 3.6 m telescope to map H-alpha + (N II) lambda-lambda 6548, 6583 emission-line profiles across the entire edge-on, nearby SBc galaxy NGC 3079, with resolution 70 km/s and subarcsecond sampling. Blue and red long-slit spectra were obtained with the Double Spectrograph on the Palomar 5 meter telescope to provide additional emission-line diagnostics. A spectacular, line emitting bubble of diameter 13 sec (approximately 1.1 kpc) is observed immediately east of the nucleus. Its unusual gaseous excitation (e.g., (N II) lambda(6583)/H-alpha greater than 1) suggests that shocks are important. Extremely violent gas motions that range over 2000 km/s are detected across the bubble and diametrically opposite on the west side of the nucleus. Nonrotational motions are also found in the inner galaxy disk. The superbubble of NGC 3079 is the most powerful example known of a wind-blown bubble, and an excellent laboratory to study wind dynamics. The dimensions and energies of the bubble imply that is likely to be in the blowout phase and partially ruptured. The predicted rate of kinetic energy output from the central starburst appears sufficient to power most of this outflow. It is possible that a central active galactic nucleus also contributes to the outflow.
Directional locking and hysteresis in stripe- and bubble-forming systems on one-dimensional periodic substrates with a rotating drive
We examine the dynamics of a two-dimensional stripe, bubble, and crystal forming system interacting with a periodic one-dimensional substrate under an applied drive that is rotated with respect to the substrate periodicity direction 𝑥. We find that the stripes remain strongly directionally locked to the 𝑥 direction for an extended range of drives before undergoing motion parallel to the drive. In some cases, the stripes break apart at the unlocking transition, but they can dynamically reform into stripes aligned perpendicular to the 𝑥 direction, producing hysteresis in the directional locking and unlocking transitions. In contrast, moving anisotropic crystal and bubble phases exhibit weaker directional locking and reduced or no hysteresis. The hysteresis occurs in regimes where the particle rearrangements occur and is most pronounced near the stripe phase. In conclusion, we also show that for varied substrate strength, substrate spacing, and particle density, a number of novel dynamical patterns can form that include a combination of stripe, bubble, and crystal morphologies.
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.
In-situ irradiation of uranium carbide
Uranium carbide (UC) is a leading candidate fuel for Generation IV reactors due to its high uranium density and thermal conductivity. However, its irradiation performance—particularly gas bubble swelling and defect dynamics—remains poorly characterized. Using in-situ transmission electron microscopy (TEM), we irradiated UC with 300 keV Xe + and 1 MeV Kr 2+ ions at temperatures up to 900 °C to quantify swelling behavior and dislocation loop evolution. The swelling remained below 0.6 % across all temperatures, suggesting the dimensional stability of UC under irradiation at these temperatures. Dislocation loops grew faster in UC than in UO 2 or UN, correlating with its lower homologous temperature. Notably, nanograin structures emerged in thin regions of the lamellar, mirroring phenomena previously observed in UO 2 and ZrC. These results address critical knowledge gaps in the radiation tolerance of UC and provide insight into its suitability for advanced reactor systems.
A improved equation of state for Xe gas bubbles in gamma U-Mo fuels
A monolithic fuel design based on a U-Mo alloy has been selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). An issue with U-Mo monolithic fuel is the large amount of swelling that takes place during operation. The accurate prediction of fuel evolution under irradiation requires implementation of correct thermodynamic properties into mesoscale and continuum level fuel performance modeling codes. However, the thermodynamic properties of the fission gas bubbles (such as the relationship among bubble size, equilibrium Xe concentration, and bubble pressure) are not well known. This work studies Xe bubbles in ?U-Mo from a diameter of 3 nm up to 8.5 nm and from 400 K up to 700 K. The energetic relationship of Xe bubbles with regard to voids and Xe substitutional atoms is described. The transition is also determined for when a bubble becomes over-pressurized. Finally, an equation of state is fit to the pressure as a function of molar volume and temperature.
Long-lived oscillations of metastable states in neutral atom systems
Metastable states arise in a range of quantum systems and can be observed in various dynamical scenarios, including decay, bubble nucleation, and long-lived oscillations. The phenomenology of metastable states has been examined in quantum many-body systems, notably in one-dimensional (1D) ferromagnetic Ising spin systems and superfluids. Here, in this paper, we study long-lived oscillations of metastable and ground states in 1D antiferromagnetic neutral atom chains with long-range Rydberg interactions. We use a staggered local detuning field to achieve confinement. Using theoretical and numerical models, we identify novel spectral signatures of quasiparticle oscillations distinct to antiferromagnetic neutral atom systems and interpret them using a classical energy model of short-range meson repulsion. Finally, we evaluate the experimental accessibility of our proposed setup on current neutral-Atom platforms and discuss experimental feasibility and constraints.
Probability of Decompression Sickness and Venous Gas Emboli from 49 NASA Hypobaric Chamber Tests with Reference to Exploration Atmosphere
Introduction: Decompression sickness (DCS) is a complex biophysical event; it combines human perception of pain, for instance, and the presence of a gas phase in the tissues. Living tissues are complex and dynamic. Micronuclei and later bubbles may or may not form given what appears to be the same conditions. Even when bubbles grow, symptoms may or may not develop under what appears to be the same conditions. Therefore, at this time it is appropriate to consider DCS as a probabilistic rather than a deterministic event. Methods: Probabilistic models about hypobaric DCS and venous gas emboli (VGE) require a large amount of quality research data, a definition of decompression dose using physical and physiologic variables, and a flexible analytical approach that can quantify the association between each outcome and all covariates of interest (assuming independence between DCS and VGE) and then ultimately be extended to acknowledge dependencies between DCS and VGE. Our DCS and VGE data are from 1,031 hypobaric decompressions from 1983 to 2016. A total of 577 humans participated in 49 hypobaric chamber tests to evaluate denitrogenation procedures used by astronauts in the Space Shuttle and International Space Station programs. We defined decompression dose as the ratio of computed nitrogen tension in a theoretical 360-minute half-time compartment to ambient pressure, which accounts for denitrogenation and exposure pressure as well as explanatory variables such as age, sex, body mass index, and the presence or absence of ambulation as part of exercise at the exposure pressure. A parametric survival model, using a log-logistic distribution, was used to quantify the time to development of DCS, VGE, and Grade IV VGE. Results: Our survival estimates are applicable to simple hypobaric decompressions, such as depressurizations in 5 to 30 minutes to exposure pressures between 4 and 10 pounds per square inch absolute (psia) and after minutes to hours of denitrogenation, either under resting or exercise conditions to accelerated denitrogenation. The regressions are applicable to exposures between 2 to 6 hours and under conditions of ambulation or no ambulation as part of exercise at the test pressure. We estimate that an exposure to 4.3 psia with simulated extravehicular activity (EVA) that includes ambulation after equilibration to the exploration atmosphere at 8.2 psia with a 34% oxygen atmosphere will result in 3.1% DCS (1.8% to 5.2 95% confidence interval), 23.2% VGE (16.7 to 31.2%), and 8.5% Grade IV VGE (4.7 to 14.7%) in equal samples of men and women exposed for 6 hours. Discussion: Probabilistic models for DCS, VGE, and Grade IV VGE can be used to inform those that plan future EVAs. Their applications are useful to quantify the risk of DCS and VGE in astronauts that perform EVAs in low-pressure space suits while in space or while exploring the surfaces of the moon or Mars.
Shallow Bubble Collapse and AESOP: New Ejecta Mechanisms and Specialized Diagnostics
Presentation to be given on May 18 at the MSTS Works in Progress Seminar, virtually. The presentation will also be archived on the Science & Technology website.
Dynamics of Aqueous Foam Drops
We develop a model for the nonlinear oscillations of spherical drops composed of aqueous foam. Beginning with a simple mixture law, and utilizing a mass-conserving bubble-in-cell scheme, we obtain a Rayleigh-Plesset-like equation for the dynamics of bubbles in a foam mixture. The dispersion relation for sound waves in a bubbly liquid is then coupled with a normal modes expansion to derive expressions for the frequencies of eigenmodal oscillations. These eigenmodal (breathing plus higher-order shape modes) frequencies are elicited as a function of the void fraction of the foam. A Mathieu-like equation is obtained for the dynamics of the higher-order shape modes and their parametric coupling to the breathing mode. The proposed model is used to explain recently obtained experimental data.
The effect of composition on helium bubble nucleation in complex concentrated alloys
Molecular dynamics simulations of helium implantation have been performed in the MoNbTaTi complex concentrated alloy (CCA) to determine the impact of composition on helium clustering and bubble nucleation. It was found that compositions with lower interstitial atomic volume resulted in smaller cluster sizes, higher helium-to-vacancy ratios, and larger migration barriers. As the atomic volume decreases, it is harder to accommodate additional helium in the lattice and for helium to diffuse between interstitial sites, resulting in smaller cluster sizes. In particular, increasing molybdenum directly correlated with a decreased interstitial atomic volume and increased migration barrier. Niobium exhibited the opposite trend, where more niobium in the material resulted in larger interstitial atomic volumes and cluster sizes but lower migration barriers. The atomic volume was determined to be an indicator of the susceptibility of the material to helium bubble growth and could potentially be used as a metric to screen compositions for resistance to helium damage. In this way, properties such as available interstitial volumes could be used to screen materials more rapidly given the vast compositional space of CCAs.
Predicting critical heat flux using localized sensing at invisible vapor-liquid interfaces
Predicting critical heat flux (CHF) in two-phase electronics cooling systems remains a significant challenge due to the sudden onset of boiling crisis and the difficulty in directly visualizing vapor-liquid interfaces. Existing sensing methods rely on lagging temperature measurements, optically accessible systems, or spatially averaged signals that cannot pinpoint CHF initiation at localized high-heat-flux regions. Here, we report a planar capacitive sensing approach that enables real-time, localized detection of vapor-liquid interface dynamics for CHF prediction in boiling heat transfer. The capacitive sensor exploits the dielectric constant difference between liquid and vapor phases to capture bubble nucleation, growth, and departure dynamics with a temporal resolution down to 2 ms. The capacitive sensing reveals distinct signals across boiling regimes: from high-frequency fluctuations during strong nucleate boiling to low-frequency fluctuations with increased amplitudes when approaching CHF. The multi-sensor array experiments demonstrate real-time localized sensing, where each sensor responds exclusively to boiling in its immediate vicinity without crosstalk from neighboring regions. This non-intrusive sensing approach provides predictive rather than lagging sensing signals of CHF occurrence, offering predictive diagnosis of two-phase liquid cooling for the thermal management of high-power-density electronics.
Water masers in W49N - The youngest stellar jet?
Observations by Gwinn et al. (1992) of the proper motions of water masers in W49N show that they have an elongated distribution expanding from a common center. Features with high space velocity only occur far from the center, while low-velocity features occur at all distances. It is proposed that water masers in star-forming regions occur in expanding shells swept up by high-velocity winds from young, massive stars during the early phases of the expansion. In W49N, confinement of the bubble by a density distribution with an axial cavity can explain both the velocity field and the shape of the maser distribution. A fully dynamical calculation of the expanding bubble is presented which provides a satisfactory fit for the observations and suggests that this system is only about 250 yr old. Thus these observations may show the very first stages of the formation of a jet from a young stellar object.
Dynamics of surface tension in microgravity environment
Time-dependent evolutions of the profile of free surface (bubble shapes) for a cylindrical container partially filled with a Newtonian fluid of constant density, rotating about its axis of symmetry, in low and microgravity environments, have been studied. Numerical computations of the dynamics of bubble shapes have been carried out (1) linear time-dependent functions of spin-up and spin-down in low and microgravity environments, (2) linear time-dependent functions of increasing and decreasing gravity environment at high and low rotating cylinder speeds, (3) time-dependent step functions of spin-up and spin-down in a low-gravity environment, and (4) sinusoidal function oscillation of the gravity environment at high and low rotating cylinder speeds.
Energy-conserving physics for nonhydrostatic dynamics in mass coordinate models
Abstract. Motivated by reducing errors in the energy budget related to enthalpy fluxes within the Energy Exascale Earth System Model (E3SM), we study several physics–dynamics coupling approaches. Using idealized physics, a moist rising bubble test case, and the E3SM's nonhydrostatic dynamical core, we consider unapproximated and approximated thermodynamics applied at constant pressure or constant volume. With the standard dynamics and physics time-split implementation, we describe how the constant-pressure and constant-volume approaches use different mechanisms to transform physics tendencies into dynamical motion and show that only the constant-volume approach is consistent with the underlying equations. Using time step convergence studies, we show that the two approaches both converge but to slightly different solutions. We reproduce the large inconsistencies between the energy flux internal to the model and the energy flux of precipitation when using approximate thermodynamics, which can only be removed by considering variable latent heats, both when computing the latent heating from phase change and when applying this heating to update the temperature. Finally, we show that in the nonhydrostatic case, for physics applied at constant pressure, the general relation that enthalpy is locally conserved no longer holds. In this case, the conserved quantity is enthalpy plus an additional term proportional to the difference between hydrostatic pressure and full pressure.