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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Coupled Experimental/Computational Investigation of the Dynamics of Interacting Magnetized Plasmas

The interaction, or interpenetration, of magnetized plasmas of different density and/or pressure occurs in a wide variety of natural and man-made systems. Such systems include extragalactic jets propagating into the intergalactic medium, solar coronal mass ejections into background solar wind, compact toroid (CT) fueling of magnetic fusion plasmas, and jets of capsule shell impurities into DT fusion fuel, which can lead to enhanced impurity mix in inertial fusion implosions. These plasmas may take the form of jets, with open, helical magnetic structure, or plasma “bubbles” with closed magnetic fields (B-fields), such as spheromaks or CT’s. Such structures, both open and closed B-field cases, can transport heat, particles and magnetic flux or magnetic helicity into background plasma regions. For example, the origin of extragalactic magnetic fields may be due, at least in part, to transport by astrophysical jets. The goal of this proposed work was to elucidate the detailed plasma and magnetic field dynamics of high-density plasma jets (open B-field) and bubbles (closed B-field) propagating into lower density background magnetized plasma through controlled laboratory experiments and closely coupled nonlinear MHD modeling. These experiments were conducted in the HelCat (Helicon-Cathode) linear plasma device at the University of New Mexico (UNM). Plasma jets and bubbles were launched via an existing compact coaxial plasma gun, mounted on the HelCat device. This gun produced plasmas tens of cm in scale and lasting tens of microseconds, thereby allowing detailed multipoint, space- and time-resolved measurements to be made routinely. The experiments were directly modeled using the extended magnetohydrodynamic (XMHD) PERSEUS code, developed at Cornell University [23,24]. Both experimental and numerical modeling work are ongoing. The main results to date are reported here. Additional supplemental funding for one year (8/1/2019 – 7/31/2020) supported numerical investigation of photoionization processes important in many low temperature plasmas, including the HelCat device. Initial results of this modeling work is also reported.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Large Eddy Simulation Study of Flow Turbulence, Alumina Transport, and Bath Temperature Evolution in Conventional Aluminum-Smelting Cell Using OpenFOAM

In this study, a Large Eddy Simulation (LES) of the aluminum-smelting process is performed using OpenFOAM. To understand the coupled behavior of heat transfer, mass transfer, and flow of the smelting process, a multi-physics computational fluid dynamics (CFD) model based on the Eulerian–Eulerian multi-fluid approach is adopted. The model accounts for CO 2 bubble and magnetohydrodynamics (MHD)-driven flow, along with alumina dissolution, transport, and bath temperature evolution. The simulation predictions show small-scale turbulent vortical structures in the anode–cathode space caused by combined effect of MHD and CO 2 bubble-bath interactions and relatively large-scale asymmetric vortices in the inter-anode space caused by the CO 2 bubble-bath interactions. The vortex formation at the edges of the anodes evidently aids in transporting alumina from the central channel to the bottom of the anodes and prevents accumulation of gas bubbles in the periphery of the anode bottom. Symmetric bath cold spots are observed in the vicinity of the feeder. Cold spots are also observed in the anode–cathode distance space below the anode bottom due to the transport of undissolved solid to this region by the flow. The findings from the work are useful in developing and designing alumina-feeding strategy leading to reduced anode effects and smooth operation of the cell. The work also highlights the important flow structures in conventional aluminum-smelting cell.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Comparison of Interface Growth Models Applied to Rayleigh–Taylor and Richtmyer–Meshkov Instabilities

Sophisticated numerical models that contain fluid interfaces rely upon interface evolution models to approximate the transition to turbulence near interfaces, in the presence of Rayleigh–Taylor (RTI) or Richtmyer–Meshkov (RMI) instability. Semi-analytical models have been developed in recent decades to predict the interface growth from an initial state into the nonlinear regime. Two of these models are considered in this study. They are the Goncharov and the z-models. Both of these interface models have strengths and weaknesses, which are examined here. Both of them have been implemented in the xRAGE compressible flow solver, which models fluid interfaces. The flow solver provides the fluids acceleration as a body force to the interface model. The purpose of such interface model is to evolve the early times of interface position as a subgrid model within a compressible flow simulation in order to then initialize a turbulence model. In this work, the interface models are assessed and compared for their evolution of RTI and RMI. The z-model performed better than the Goncharov model for all cases that were explored.

42 ENGINEERING↗
NASA NTRS1985

Formation of metallic and metallic-glass hollow spheres and their solidification characteristics

Various metals and metallic glass systems have bene processed into hollow spheres with sizes ranging from 3 mm to 440 microns in diameter. The technique for the formation of the large hollow spheres, in general, is based on the fluid-dynamic instability of a hollow annular jet. A refined technique has also been developed for microshell formation, in which discrete bubbles are injected into the stream of the molten material and individually 'flushed' out at a frequency related to the Rayleigh jet instability. The surfaces of those spheres of all sizes exhibit a range of contrasting solidification behaviors and characteristics. Metal shells of varying materials, sizes, aspect ratios, sphericity and concentricity have many useful and novel applications.

Lee, M. C.↗
NASA NTRS1990

A mixed element for laminated plates and shells

Formulation and numerical evaluation of a simple shear-flexible four-noded quadrilateral laminated composite plate/shell finite-element is presented. The element developed is based on a generalized mixed variational principle with independently assumed displacement and laminate internal strain fields. For the latter, a layer-number-independent polynomial interpolation is utilized, together with a judicious selection of the in-plane (spanwise) distributions for the strain components. This was facilitated by the use of a set of bubble functions as additional kinematic degrees of freedom, which was shown to be crucial for eliminating the locking phenomenon. For dynamic applications, a simplified lumped mass matrix is employed. Finally, an extensive number of critical test cases are given to access the element's preformance and to illustrate its effectiveness in static as well as vibration problems for anisotropic laminated plates and shells.

Wilt, T. E.↗
DOE OSTI2021

Heat transfer from glass melt to cold cap: Computational fluid dynamics study of cavities beneath cold cap

Efficient glass production depends on the continuous supply of heat from the glass melt to the floating layer of batch, or cold cap. Computational fluid dynamics (CFD) are employed to investigate the formation and behavior of gas cavities that form beneath the batch by gases released from the collapsing primary foam bubbles, ascending secondary bubbles, and in the case of forced bubbling, from the rising bubbling gas. The gas phase fraction, temperature, and velocity distributions below the cold cap are used to calculate local and average heat transfer rates as a function of the bubbling rate. It is shown that the thickness of the cavities is nearly independent of the cold cap shape and the amount of foam evolved during batch conversion. Furthermore, it is ~7 mm and up to ~15 mm for the cases without and with forced bubbling used to promote circulation within the melt, respectively. Using computed velocity and temperature profiles, the melting rate of the simulated high-level nuclear waste glass batch was estimated to increase with the bubbling rate to the power of ~0.3 to 0.9, depending on the flow pattern. The simulation results are in good agreement with experimental data from laboratory- and pilot-scale melter tests.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2002

A Comparative Study of Three Methodologies for Modeling Dynamic Stall

During the past two decades, there has been an increased reliance on the use of computational fluid dynamics methods for modeling rotors in high speed forward flight. Computational methods are being developed for modeling the shock induced loads on the advancing side, first-principles based modeling of the trailing wake evolution, and for retreating blade stall. The retreating blade dynamic stall problem has received particular attention, because the large variations in lift and pitching moments encountered in dynamic stall can lead to blade vibrations and pitch link fatigue. Restricting to aerodynamics, the numerical prediction of dynamic stall is still a complex and challenging CFD problem, that, even in two dimensions at low speed, gathers the major difficulties of aerodynamics, such as the grid resolution requirements for the viscous phenomena at leading-edge bubbles or in mixing-layers, the bias of the numerical viscosity, and the major difficulties of the physical modeling, such as the turbulence models, the transition models, whose both determinant influences, already present in static maximal-lift or stall computations, are emphasized by the dynamic aspect of the phenomena.

Sankar, L.↗
NASA NTRS1998

Dryout and Rewetting in the Pool Boiling Experiment Flown on STS-72 (PBE-2 B) and STS-77 (PBE-2 A)

Experiments were conducted in the microgravity of space in which a pool of liquid (R-113), initially at a precisely defined pressure and temperature, is subjected to a step imposed heat flux from a semi-transparent thin-film heater forming part of one wall of the container such that boiling is initiated and maintained for a defined period of time at a constant pressure level. A total of nine tests were conducted at three levels of heat flux and three levels of subcooling in each of the two space experiments in a GAS canister on the STS-77, -72, respectively. Three (3) modes of propagation of boiling across the heater surface and subsequent vapor bubble growths were observed, in addition to the two (2) modes observed in the previous microgravity pool boiling space flights on STS-47, -57, and -60. Of particular interest were the extremely dynamic or "explosive" growths, which were determined to be the consequence of the large increase in the liquid-vapor interface area associated with the appearance of a corrugated or rough interface. Predictions of circumstances for its onset have been carried out. Assumptions were necessary regarding the character of disturbances necessary for the instabilities to grow. Also, a new vapor bubble phenomena was observed in which small vapor bubbles migrated toward a larger bubble, eventually coalescing with this larger bubble. The heat transfer was enhanced approximately 30% as a result of these migrating bubbles, which is believed to be a vapor bubble manifestation of Marangoni convection and/or molecular momentum effects, sometimes referred to as vapor recoil. The circumstances of heat flux and liquid subcooling necessary to produce heater surface dryout for an initially stagnant liquid subjected to an imposed heat flux have been more closely identified.

Merte, Herman, Jr.↗
NASA NTRS1978

Dynamic stall experiments on the NACA 0012 airfoil

The flow over a NACA 0012 airfoil undergoing large oscillations in pitch was experimentally studied at a Reynolds number of and over a range of frequencies and amplitudes. Hot-wire probes and surface-pressure transducers were used to clarify the role of the laminar separation bubble, to delineate the growth and shedding of the stall vortex, and to quantify the resultant aerodynamic loads. In addition to the pressure distributions and normal force and pitching moment data that have often been obtained in previous investigations, estimates of the unsteady drag force during dynamic stall have been derived from the surface pressure measurements. Special characteristics of the pressure response, which are symptomatic of the occurrence and relative severity of moment stall, have also been examined.

Mcalister, K. W.↗
NASA NTRS1999

Dynamic Nucleation of Supercooled Melts and Measurement of the Surface Tension and Viscosity

We investigate the phenomenon of acoustic pressure-induced nucleation by using a novel approach involving the large amplitude resonant radial oscillations and collapse of a single bubble intentionally injected into a supercooled liquid. Using a combination of previously developed and proven techniques, the bubble is suspended in a fluid host by an ultrasonic field which supplies both the levitation capability as well as the forcing of the radial oscillations. We observe the effects of an increase in pressure (due to bubble collapse) in a region no larger than 100 microns within the supercooled melt to rigorously probe the hypothesis of pressure-induced nucleation of the solid phase. The use of single bubbles operating in narrow temporal and spatial scales will allow the direct and unambiguous correlation between the origin and location of the generation of the disturbance and the location and timing of the nucleation event. In a companion research effort, we are developing novel techniques for the non-contact measurements of the surface tension and viscosity of highly viscous supercooled liquids. Currently used non-invasive methods of surface tension measurement for the case of undercooled liquids generally rely of the quantitative determination of the resonance frequencies of drop shape oscillations, of the dynamics of surface capillary waves, or of the velocity of streaming flows. These methods become quickly ineffective when the liquid viscosity rises to a significant value. An alternate and accurate method which would be applicable to liquids of significant viscosity is therefore needed. We plan to develop such a capability by measuring the equilibrium shape of levitated undercooled melt droplets as they undergo solid-body rotation. The experimental measurement of the characteristic point of transition (bifurcation point) between axisymmetric and two-lobed shapes will be used to calculate the surface tension of the liquid. Such an approach has already been validated through the experimental verification of numerical modeling results. The experimental approach involves levitation, melting, and solidification of undercooled droplets using a hybrid ultrasonic-electrostatic technique in both a gaseous as well as a vacuum environment. A shape relaxation method will be investigated in order to derive a reliable method to measure the viscosity of undercooled melts. The analysis of the monotonic relaxation to equilibrium shape of a drastically deformed and super-critically damped free drop has been used to derive interfacial tension of immiscible liquid combinations where one of the component has high viscosity. A standard approach uses the initial elongation of a droplet through shear flows, but an equivalent method could involve the initial deformation of a drop levitated in a gas by ultrasonic radiation pressure, electric stresses, or even solid body rotation. The dynamic behavior of the free drop relaxing back to equilibrium shape will be modeled, and its characteristic time dependence should provide a quantitative means to evaluate the liquid viscosity.

Trinh, E. H.↗
NASA NTRS1985

Analysis of an oscillatory oil squeeze film containing a central gas bubble

A squeeze-film damper, consisting of two circular plates, having only normal oscillatory relative motion is considered. The liquid lubricant between the plates is assumed to contain a single central gas bubble. The effect of the bubble on the damper performance is analyzed. Comparison is made with the performance of a pure liquid damper. Substantial deviations in peak dynamic pressures are predicted which explain discrepancies between experimental and theoretical results reported in the literature.

Haber, S.↗
NASA NTRS1991

Accretion dynamics and polarized X-ray emission of magnetized neutron stars

The basic ideas of accretion onto magnetized neutron stars are outlined. These are applied to a simple model of the structure of the plasma mound sitting at the magnetic poles of such a star, in which upward diffusion of photons is balanced by their downward advection. This steady flow model of the plasma's dynamical state is used to compute the emission of polarized X-raysfrom the optically thick, birefringent medium. The linear polarization of the continuum radiation emerging from the quasi-static mound is found to be as much as 40 percent at some rotation phases, but is insensitive to the geometry of the accretion flow. The role of the accretion shock, whose detailed polarimetric and spectral characteristics have yet to be calculated, is emphasized as the final determinant of the properties of the emerging X-rays. Some results describing the fully time dependent dynamics of the flow are also presented. In particular, steady flow onto a neutron star is shown to exhibit formation of 'photon bubbles', regions of greatly reduced plasma density filled with radiation which form and rise on millisecond time scale. The possible role of these complex structures in the flow for the formation of the emergent spectrum is briefly outlined.

Arons, Jonathan↗
NASA NTRS1993

Liquid Oxygen (LO2) propellant conditioning concept testing

Marshall Space Flight Center (MSFC) and industry contractors have undertaken activities to develop a simplified liquid oxygen (LO2) propellant conditioning concept for future expendable launch vehicles. The objective of these activities is to reduce operations costs and timelines and to improve safety of these vehicles. The approach followed has been to identify novel concepts through system level studies and demonstrate the feasibility of these concepts through small-scale and full-scale testing. Testing will also provide data for design guidelines and validation of analytical models. Four different concepts are being investigated: no-bleed, low-bleed, use of a recirculation line, and helium (He) bubbling. This investigation is being done under a Joint Institutional Research and Development (JIRAD) program currently in effect between MSFC and General Dynamics Space Systems (GDSS). A full-scale test article, which is a facsimile of a propellant feed duct with an attached section to simulate heat input from a LO2 turbopump, will be tested at the Cold Flow Facility at MSFC's West Test Area. Liquid nitrogen (LN2), which has similar properties to LO2, will be used in place of LO2 for safety and budget reasons. Work to date includes design and fabrication of the test article, design of the test facility and initial fabrication, development of a test matrix and test procedures, initial predictions of test output, and heat leak calibration and heat exchanger tests on the test article. The tests for all propellant conditioning concepts will be conducted in the summer of 1993, with the final report completed by October, 1993.

Perry, Gretchen L. E.↗
DOE OSTI2025

In-situ ion irradiation of a spent UO2 fuel: evolution of fission products and nanograins with radiation dose

This study investigates the behavior of fission gas bubbles, five metal precipitates (5MPs) (Mo, Ru, Rh, Tc, Pd) and nanograins in spent uranium dioxide (UO2) fuel under various ion irradiation doses, temperatures, and flux conditions. Utilizing in-situ ion irradiation and advanced transmission electron microscopy, we analyzed the evolution of fission gas bubbles, 5MPs and nanograins in UO2 samples from the Belgium Reactor 3 (BR-3). Our findings reveal significant shrinkage of fission gas bubbles and 5MPs with increasing irradiation dose, accompanied by a decrease in pair density. The study also reveals nanograin coarsening under irradiation. We demonstrate that ion irradiation induces a homogeneous re-solution process where individual atoms are ejected from bubbles and precipitates, leading to their dissolution and subsequent re-precipitation in the matrix. This study provides critical insights into the dynamic behavior of fission products and nanograins under irradiation, facilitating the development of predictive models and contributing to the optimization of nuclear fuel performance and safety.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Helium bubble size effects on the surface morphological response of plasma-facing tungsten

Abstract We report a simulation study on the effects of helium (He) bubble size on the morphological evolution and pattern formation on the surface of tungsten used as a plasma-facing component (PFC) in nuclear fusion devices. We have carried out a systematic investigation based on self-consistent dynamical simulations of surface morphological evolution according to an atomistically-informed, 3D continuum-scale model that captures well the relevant length and time scales of surface nanostructure formation in PFC tungsten. The model accounts for PFC surface diffusion, driven by the biaxial compressive stress originating from the over-pressurized He bubbles in the near-surface region of PFC tungsten as a result of He plasma exposure, combined with the formation of self-interstitial atoms in tungsten that diffuse toward the PFC surface and the flux of surface adatoms generated as a result of surface vacancy-adatom pair formation upon He implantation; this transport of surface adatoms contributes to the anisotropic growth of surface nanostructural features due to the different rates of adatom diffusion along and across step edges of islands on the tungsten surface. Our detailed analysis reveals that varying the average He bubble size plays an important role in the PFC surface growth kinetics as well as the resulting surface topography. Specifically, we find that an increase in the He bubble size leads to a deceleration in the growth rate of the tungsten nanotendrils that emanate from the PFC surface. We also find that the separation distance between the resulting surface features increases with increasing He bubble size, as well as over time. This coarsening effect is a thermally activated process resulting in an accurate description of the temperature dependence of the average surface feature separation by an Arrhenius relation.

36 MATERIALS SCIENCE↗
NASA NTRS2004

Effects of Gravity on Bubble Formation in an Annular Jet

The effects of gravity on the bubble formation in an annular jet were studied. The experiments were conducted in the 2.2-second drop tower at the NASA Glenn Research Center. Terrestrial gravity experiments were conducted at the Fluid Dynamics Research Laboratory at the University of Oklahoma. Stainless steel tubing with inner diameters of 1/8" (gas inner annulus) and 5/16" (liquid outer annulus) served as the injector. A rectangular test section, 6" x 6" x 14" tall, made out of half-inch thick Lexan was used. Images of the annular jet were acquired using a high-speed camera. The effects of gravity and varying liquid and gas flow rates on bubble size, wavelength, and breakup length were documented. In general, the bubble diameter was found to be larger in terrestrial gravity than in microgravity for varying Weber numbers (0.05 - 0.16 and 5 - 11) and liquid flow rates (1.5 ft/s - 3.0 ft/s). The wavelength was found to be larger in terrestrial gravity than in microgravity, but remained constant for varying Weber numbers. For low Weber numbers (0.05 - 0.16), the breakup length in microgravity was significantly higher than in terrestrial gravity. Comparison with linear stability analysis showed estimated bubble sizes within 9% of experimental bubble sizes. Bubble size compared to other terrestrial gravity experiments with same flow conditions showed distinct differences in bubble size, which displayed the importance of injector geometry on bubble formation.

Koepp, R. A.↗
DOE OSTI2020

Comparison of Surface Tension Generation Methods in Smoothed Particle Hydrodynamics for Dynamic Systems

Developing robust numerical models of dynamic surface tension dominated multiphase systems is an ongoing challenge, especially in scenarios with large density and viscosity ratios. This is critical to the design and understanding of various physical and engineering systems, such as fluidized beds, fuel injectors, and drug delivery schemes. Much of the computational work in surface tension dominated multiphase flows has employed the continuum surface force method (CSF) of Brackbill et al. [1], which recasts surface tension from a surface force to a volumetric force that can be imposed in the vicinity of an interface. The CSF method produces accurate results across a variety of systems, however it relies on the identication of surface normals, which can be unreliable under certain conditions. Alternative methods of simulating surface tension have been proposed. Here the advantages and disadvantages of the CSF method in comparison to a pairwise forces (PF) method proposed by Tartakovsky and Panchenko [2] are explored. The CSF and PF methods are used in a smoothed particle hydrodynamics (SPH) framework to model dynamic systems. Results are compared to existing test cases from the literature and to analytic solutions derived from fundamental normal mode behavior of bubbles and droplets. A physical system for which the PF method is more stable and physically appropriate than the CSF method is identied.

smoothed particle hydrodynamics, multiphase flow↗
DOE OSTI2023

Compressible Rayleigh–Taylor instability subject to isochoric initial background state

The effect of compressibility on the single-mode Rayleigh–Taylor instability is examined using two (2D) and three-dimensional (3D) direct numerical simulations. To isolate compressibility from background stratification effects, this work employs a constant density profile on each side of the interface. The numerical simulations are performed at various Reynolds numbers using the gas kinetic method for static Mach numbers up to M = 0.4. The most important finding is that compressibility acting in isolation enhances the instability and perturbations grows faster with increasing Mach number, unlike previous results with background isothermal state, which show suppression of the instability at higher static Mach numbers. In addition, compressibility is also shown to increase the bubble-spike asymmetry. While the instability grows faster for the 3D case, the findings are qualitatively similar in 2D and 3D. Furthermore, the dynamical reasons underlying the effect of compressibility are elucidated by examining the evolution of vorticity and turbulent kinetic energy transport equations.

42 ENGINEERING↗