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 415 records · Page 23

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.

Tryggvason, Gretar↗

Experiments on annular liquid jet instability and on the formation of liquid shells

An annular jet flow of liquid surrounding a flow of gas at its core is extremely unstable. Experiments are described in which such a flow is generated by an annular nozzle operated at fairly specific conditions. It is shown that periodic, axisymmetric oscillations arise spontaneously within the cylindrical sheet emerging from the nozzle and grow with such rapidity along the axial dimension that a sealing-off and encapsulation of the core gas occurs within a few jet diameters. This is closely followed by a pinchoff of the liquid between adjacent bubbles. The liquid shells set free thereby assume spherically symmetric form under capillary forces, and each contains a precisely uniform measure of gas and of liquid on account of the extremely high frequency-stability of the process. Description is given of the fluid dynamic processes by which the shells are formed, and mention is made of exploiting the instability for the production of rigid shells for technological applications.

Kendall, J. M.↗

Comparison of Oxygen Liquefaction Gaseous Oxygen Injection via Ullage and Dip Tube

The Cryogenic Fluid In-situ Liquefaction for Landers (CryoFILL) testing demonstrated constant and transient oxygen liquefaction in a scaled lander tank integrated with an industrial cryocooler. Early in the test series, two methods for injecting gaseous oxygen into the test tank were ran under similar constant liquefaction conditions for comparison. Gaseous oxygen was either injected directly from the lid into the ullage or bubbled up through the liquid using a dip tube that extends to the bottom of the test tank. The direct comparison of these gaseous oxygen injection methods for constant liquefaction indicated ullage injection provided a greater liquefaction rate than injection via a dip tube. Injection through the ullage was subsequently selected as the nominal injection method for the CryoFILL test series. The comparison of the two injection methods will be highlighted with a discussion of what may cause the variation in liquefaction rate.

Cryogenic Fluid Management↗

Effect of Baffle on Gravity-Gradient-Excited Slosh Waves and Spacecraft Moment and Angular-Momentum Fluctuations in Microgravity

The dynamical behavior of fluids affected by the asymmetric gravity gradient acceleration has been investigated. In particular, the effects of surface tension on partially filled rotating fluids applicable to a full-scale Gravity Probe-B Spacecraft dewar tank with and without baffles are studied. Results of slosh wave excitation along the liquid-vapor interface induced by gravity gradient acceleration indicate that the gravity gradient acceleration is equivalent to the combined effect of a twisting force and a torsional moment acting on the spacecraft. The results are clearly seen from one-up one-down and one-down one-up oscillations in the cross-section profiles of two bubbles in the vertical (r, z)-plane of the rotating dewar, and from the eccentric contour of the bubble rotating around the axis of the dewar in a horizontal (r, theta)-plane. As the viscous force, between liquid and solid interface, greatly contributes to the damping of slosh wave excitation, a rotating dewar with baffles provides more areas of liquid-solid interface than that of a rotating dewar without baffles. Results show that the damping effect provided by the baffles reduces the amplitude of slosh wave excitation and lowers the degree of asymmetry in liquid-vapor distribution. Fluctuations of angular momentum and fluid moment caused by the slosh wave excited by gravity gradient acceleration with and without baffle boards are also investigated. It is also shown that the damping effect provided by the baffles greatly reduces the amplitudes of angular momentum and fluid moment fluctuations.

Hung, R. J.↗

Coupled Microbial-Conversion and Computational-Fluid-Dynamics (CFD) Models for Butanediol Production in Micro-Aerated Reactors

Microbial conversion of substrates to macromolecules has been widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries. These bioreactions are also being investigated in the production of low-value commodities such as biofuels [1]. Gas-liquid mass-transfer and transport-reaction coupling are important challenges when designing and scaling up these reactor systems. Experiments in wellmixed small-scale reactors have enabled characterization of microbial reactivity, while their coupling with macroscale transport remains relatively unexplored. In this work, we use a coupled metabolic-CFD model to study the action of a genetically engineered microbe Zymomonas mobilis [2] on sugars to produce 2,3-Butanediol (BDO). BDO is an important hydrocarbon intermediate that can be catalytically upgraded to several fuels and chemicals [3]. An important aspect to this particular microbial conversion is the need for micro-aerated environments as opposed to traditional aerobic fermentation. Slight variations in oxygen concentration can result in competing reaction pathways that disable BDO production. Hence, gas-liquid mass transfer and transport need to be optimized in large-scale reactors to maximize BDO production, for which CFD is a valuable tool. The aerobic-fermentation CFD model previously developed by the authors [5] for simulating bubble-column and airlift reactors at scale was used in this study. The Reynolds-averaged mass, momentum and energy transport equations for interpenetrating gas and liquid phase are solved in this model along with the transport and interphase mass transfer of oxygen. Our previous work used a phenomenological model for microbial oxygen uptake that neglected microbial growth and other reaction pathways. In this work, a detailed metabolic model enabled prediction of product formation and inhibition pathways. In order to manage computational cost, we used a subcycling technique [6] that takes advantage of the clear separation in transport (~ 200 sec) and reaction (~ 2-3 hours) timescales. The CFD model is first solved to steady state, after which the metabolic model is advanced at every cell in the computational domain using the local oxygen concentration. The CFD model is then run to achieve a new steady state that provides a new oxygen distribution for the metabolic model. This process, where reaction and fluid updates are interleaved together, is iterated until reactants are completely exhausted. This work will examine the performance of different reactor designs such as bubble column and airlift reactors at scale (250-500 m3). Oxygen mass-transfer coefficient and distribution are critically analyzed among reactors, and optimization studies pertaining to aeration is presented. Furthermore, it has been observed in experiments that high BDO production may be achieved by manipulating the aerobic environment over the course of reaction, such that oxygen concentration is high during the growth phase, and very low as sugar is depleted. This characteristic will be addressed by our simulations for which a timedependent scheduling strategy for aeration is presented that maximizes BDO production. [1] Humbird, D., Davis, R., and McMillan, J., Aeration costs in stirred-tank and bubble column bioreactors, Biochemical Engineering Journal, 127, 161—166, 2017 [2] Yang, S., Mohagheghi, A., Franden, M. A., Chou, Y.-C., Chen, X., Dowe, N., Himmel, M. E., and Zhang, M., Metabolic engineering of zymomonas mobilis for 2, 3-butanediol production from lignocellulosic biomass sugars. Biotechnology for biofuels, 9(1):189, 2016 [3] Kim, S. J., Sim, H. J., Kim, J. W., Lee, Y. G., Park, Y. C., and Seo, J. H., Enhanced production of 2,3-butanediol from xylose by combinatorial engineering of xylose metabolic pathway and cofactor regeneration in pyruvate decarboxylase-deficient Saccharomyces cerevisiae. Bioresource Technology, 245:1551–1557, 2017 [4] 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

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Diabatic visualization shows effects of micro-fins on evaporation of R410A: Smooth, axial micro-fin, and helical micro-fin tubes

A novel approach to visualize the flow boiling inside a clear micro-fin tube under diabatic conditions is presented. Transparent smooth, axial micro-fin, and helical micro-fin tubes are made by 3D printing. The 3D printed tube is placed inside a glass tube and heated by the transparent secondary fluid flowing between the two tubes for providing the evaporation conditions and transparency. R410A flow boiling and the flow patterns in the three geometries captured with a high speed camera are compared. The experimental results show that micro-fin geometry influences flow behavior. Bubbles are mainly generated in the groove region due to a higher superheat than the fin region. Some of the liquid refrigerant is trapped in the grooves of the upper part of the micro-fin tube when the slug flow is present, and bubbles are generated in this thin liquid layer. In addition, the bubbles in the helical micro-fin tube are easier to merge because of the obstacle of the fin geometry for the bubble flowing path. As the two bubble merge, the conversion of surface free energy causes the bubble to travel with a higher velocity.

3D printed↗

Heat pipe systems using new working fluids

The performance of a heat pipe system is greatly improved by the use of a dilute aqueous solution of about 0.0005 and about 0.005 moles per liter of a long chain alcohol as the working fluid. The surface tension-temperature gradient of the long-chain alcohol solutions turns positive as the temperature exceeds a certain value, for example about 40.degree. C. for n-heptanol solutions. Consequently, the Marangoni effect does not impede, but rather aids in bubble departure from the heating surface. Thus, the bubble size at departure is substantially reduced at higher frequencies and, therefore, increases the boiling limit of heat pipes. This feature is useful in microgravity conditions. In addition to microgravity applications, the heat pipe system may be used for commercial, residential and vehicular air conditioning systems, micro heat pipes for electronic devices, refrigeration and heat exchangers, and chemistry and cryogenics.

Chao, David F.↗

A phase-field method for boiling heat transfer

Here we present a phase field method for heat transfer in two-phase flow with boiling. The vapor/liquid interface evolution is modeled by the Cahn-Hilliard equation. The phase change rate is determined by accounting for the heat conduction balance on either vapor or liquid side of the interfacial area, depending on which side the temperature is assumed to be maintained at the saturation temperature during boiling. The velocity correction scheme proposed by Dong & Shen [27] is extended to solve the Navier-Stokes equations for a non-solenoidal velocity field, and the entropy viscosity method is employed for stabilization. The phase change model is verified by two-dimensional simulations of a vapor bubble growing in super-heated liquid and in film boiling. In both cases, mesh independence of the results is systematically performed. Subsequently, the method is applied to predict the growth of three-dimensional vapor bubble in a rectangular microchannel with boiling flow, achieving good agreement with experimental measurements and available simulation results using the level-set method. The numerical experiments demonstrate that the required mesh resolution for the phase field method is comparable with that of volume of fluid (VoF) and level-set methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The structure and dynamics of bubble-type vortex breakdown

A unique discrete form of the Navier-Stokes equations for unsteady, three-dimensional, incompressible flow has been used to study vortex breakdown numerically. A Burgers-type vortex was introduced along the central axis of the computational domain, and allowed to evolve in space and time. By varying the strength of the vortex and the free stream axial velocity distribution, using a previously developed Rossby number criterion as a guide, the location and size of the vortex breakdown region was controlled. While the boundaries of the vortex breakdown bubble appear to be nominally symmetric, the internal flow field is not. Consequently, the mechanisms for mixing and entrainment required to sustain the bubble region are different from those suggested by earlier axisymmetric models. Results presented in this study, for a Reynolds number of 200, are in good qualitative agreement with higher Reynolds number experimental observations, and a variety of plots have been presented to help illuminate the fluid physics.

Spall, R. E.↗

Open source implementation of glued sphere discrete element method and nonspherical biomass fast pyrolysis simulation

Here, in this research, a glued-sphere Discrete Element Method (DEM) method was implemented in the open-source, Computational Fluid Dynamics (CFD) software MFiX. The implementation was verified using a cylinder-wall collision and then validated by simulating the packing and fluidization of non-spherical particles. The validated code was applied to simulate fast pyrolysis of non-spherical biomass particles in a bubbling fluidized bed. The glued sphere occupancy ratio was proposed to quantify the quality of shape resolution using glued sphere. Shape resolution showed significant influence on the packing height in the simulation of particle packing and an occupancy ratio of 80 % was recommended. Its influence is minor in fully fluidized bed but can be 8 times higher in fixed fluidized bed. Three tested heat transfer models predicted similar yields of elongated biomass fast pyrolysis. The solver developed in this research can be used to simulate other multiphase reacting flows involving non-spherical particles.

42 ENGINEERING↗

Three-Dimensional Structure and Modeling of a Normal Bifurcated Shock Train from Experimental Measurements

An asymmetric, bifurcated normal shock train in a Mach 2 constant area, rectangular duct is investigated in the University of Michigan Direct Connect Isolator facility. High-speed schlieren imaging, wall static pressure measurements, surface oil flow visualization, and particle image velocimetry of a shock train are synthesized into a three-dimensional representation of shock train structure. This visualization is then used to inform the underlying flow physics of the distributed fluid dynamical processes along the structure. The detailed, three-dimensional morphology of the flow profile entering the shock train is shown to have a significant impact on the separated flow morphology of the shock train. This results in skewed supersonic core flow with alternating separation bubbles in the surrounding boundary layers. Finally, the efficacy of the prevailing pseudo-shock models in the literature is analyzed with the support of the available flow measurements. None are found to simultaneously model both the pressure and streamwise Mach number profile of the UMDCI pseudo-shock.

Hypersonics↗

Experiments and CFD simulation of mass transfer and hydrodynamics in a cylindrical bubble column

Flow visualization, Particle Image Velocimetry (PIV) and dissolved oxygen charging experiments are carried out in a small-scale cylindrical air-water bubble column at different superficial gas velocities. Void fraction profile, total gas holdup, liquid phase velocity field and volumetric mass transfer coefficient are obtained. Details of each experiment are given, with a new procedure for obtaining a 2D void fraction profile from flow visualization by image processing and time averaging which could provide useful information for CFD model validation. The experimental results are analyzed and compared with simulations where good agreement is found for both local and global flow patterns. The averaged governing equations and constitutive relations used in the simulation are discussed briefly, with emphasis on the proper comparison between simulations and experiments based on the derivation of the two-fluid model. Finally, by conducting three sets of experiments related to phase distribution, velocity field and mass transfer in the same apparatus, a validated CFD model for multiphase mass transfer is established and can be used in the future design of a multi-phase reactor.

42 ENGINEERING↗

Nonlinear bubble competition of the multimode ablative Rayleigh–Taylor instability and applications to inertial confinement fusion

The self-similar nonlinear evolution of the multimode ablative Rayleigh–Taylor instability (RTI) and the ablation-generated vorticity effect are studied for a range of initial conditions. We show that, unlike classical RTI, the nonlinear multimode bubble-front evolution remains in the bubble competition regime due to ablation-generated vorticity, which accelerates the bubbles, thereby preventing a transition into the bubble-merger regime. We develop an analytical bubble competition model to describe the linear and nonlinear stages of ablative RTI. Here, we show that vorticity inside the multimode bubbles is most significant at small scales with large initial perturbation. Since these small scales persist in the bubble competition regime, the self-similar growth coefficient ab can be enhanced by up to 30% relative to ablative bubble competition without vorticity effects. We use the ablative bubble competition model to explain the hydrodynamic stability boundary observed in OMEGA low-adiabat implosion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental study of bubble cavities attached to a rotating shaft in a reservoir

Bubble cavities formed by air entrainment and attached to a rotating shaft in an oil reservoir were studied. The cavities appear to the unaided eye as toroidal. High speed photography, however, reveals the individuality of the bubble cavities and their near solid body rotational characteristics. The cavities are distorted by the rotation effects but remain attached and tend to merge because of edge effects in the axial direction. The flow field within the reservoir is influenced by the unusual character of the two phase fluid found there; the vorticity is readily visualized. Other examples of vapor entrapment at the inlet of an eccentric rotor are also discussed. A simplified analytical method is provided, and a numerical analysis is being investigated. Vapor (void) entrainment and generation can significantly alter leakage rates and stability of seals, bearings, and dampers. Recognition of these effects in the component design systems will result only after detailed studies of the above phenomena.

Hendricks, R. C.↗

Three-Dimensional Structure and Modeling of a Normal Bifurcated Shock Train from Experimental Measurements

An asymmetric, normal bifurcated shock train in a Mach 2 constant area, rectangular duct is investigated in the University of Michigan Direct Connect Isolator facility. High-speed schlieren imaging, wall static pressure measurements, surface oil flow visualization, and particle image velocimetry of the shock train are synthesized into a three-dimensional representation. This visualization is then used to inform the underlying flow physics of the distributed fluid dynamical processes along the shock train. The detailed, three-dimensional morphology of the flow profile entering the shock train is shown to have a significant impact on the separated flow morphology of the shock train. This results in a skewed supersonic core flow with alternating separation bubbles in the surrounding boundary layers. Finally, the efficacy of the prevailing pseudo-shock models in the literature is analyzed with the support of the available flow measurements. None are found to model both the pressure and streamwise Mach number profile of the UMDCI pseudo-shock. Note: This presentation is accompanied by the video presentation that can be downloaded for viewing and is formatted as an mp4 with a runtime of 10 min 15 secs.

Hypersonics↗

Development of Hairpin Vortices in Turbulent Spots and End-Wall Transition

The end-stage phase of boundary layer transition is characterized by the development of hairpin-like vortices which evolve rapidly into patches of turbulent behavior. In general, the characteristics of the evolution form this hairpin stage to the turbulent stage is poorly understood, which has prompted the present experimental examination of hairpin vortex development and growth processes. Two topics of particular relevance to the workshop focus will be covered: 1) the growth of turbulent spots through the generatio and amalgamation of hairpin-like vortices, and 2) the development of hairpin vortices during transition in an end-wall junction flow. Brief summaries of these studies are described below. Using controlled generation of hairpin vortices by surface injection in a critical laminar boundary layer, detailed flow visualization studies have been done of the phases of growth of single hairpin vortices, from the initial hairgin generation, through the systematic generation of secondary hairpin-like flow structures, culminating in the evolution to a turbulent spot. The key to the growth process is strong vortex-surface interactions, which give rise to strong eruptive events adjacent to the surface, which results in the generation of subsequent hairpin vortex structures due to inviscid-viscuous interactions between the eruptive events and the free steam fluid. The general process of vortex-surface fluid interaction, coupled with subsequent interactions and amalgamation of the generated multiple hairpin-type vortices, is demonstrated as a physical mechanism for the growth and development of turbulent spots. When a boundary layer flow along a surface encounters a bluff body obstruction extending from the surface (such as cylinder or wing), the strong adverse pressure gradients generated by these types of flows result in the concentration of the impinging vorticity into a system of discrete vortices near the end-wall juncture of the obstruction, with the extensions of the vortices engirdling the obstruction to form "necklace" or "horseshoe" vortices. Recent hydrogen bubble and particle image visualization have shown that as Reynolds number is increased for a laminar approach flow, the flow will become critical, and a destabilization of the necklace vortices results in the development of an azimuthal waviness, or "kinks", in the vortices. These vortex kinks are accentuated by Biot-Savart effects, causing portions of a distorted necklace vortex to make a rapid approach to the surface, precipitating processes of localized, three-dimensional surface interactions. These interactions result in the rapid generation, focussing, and ejection of thin tongues of surface fluid, which rapidly roll-over and appear as hairpin vortices in the junction region. Subsequent amalgamation of these hairpin vortices with the necklace vortices produces a complex transitional-type flow. A presentation of key results from both these studies will be done, emphasizing both the ubiquity of such hairpin-type flow structures in manifold transitional-type flows, and the importance of vortex-surface interactions n the development of hairpin vortices.

Smith, Charles R.↗

Low cost silicon solar array project: Feasibility of low-cost, high-volume production of silane and pyrolysis of silane to semiconductor-grade silicon

Silicon epitaxy analysis of silane produced in the Process Development Unit operating in a completely integrated mode consuming only hydrogen and metallurgical silicon resulted in film resistivities of up to 120 ohms cm N type. Preliminary kinetic studies of dichlorosilane disproportionation in the liquid phase have shown that 11.59% SiH4 is formed at equilibrium after 12 minutes contact time at 56 C. The fluid-bed reactor was operated continuously for 48 hours with a mixture of one percent silane in helium as the fluidizing gas. A high silane pyrolysis efficiency was obtained without the generation of excessive fines. Gas flow conditions near the base of the reactor were unfavorable for maintaining a bubbling bed with good heat transfer characteristics. Consequently, a porous agglomerate formed in the lower portion of the reactor. Dense coherent plating was obtained on the silicon seed particles which had remained fluidizied throughout the experiment.

Breneman, W. C.↗