Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid interfaces”

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 811 records · Page 45

Numerical Simulation of No-Vent Chill and Fill of a Large Liquid Hydrogen Tank

To enable the design of future in-space cryogenic propellant transfer vehicles, high accuracy models of the cryogenic propellant transfer process are desired. For in-space refueling applications, accurate simulations are desired of the propellant tank chilldown and no-vent fill (NVF) process in order to estimate the amount of propellant required, to determine the maximum sustainable flow rate to fill the tank, and to establish a timeline for the filling process. NVF is a complex process that involves numerous two-phase heat and mass transfer phenomena, such as boiling at the wall and evaporation and condensation at the liquid/vapor interface. In recent years, progress has been made in developing a numerical modeling technique using the Generalized Fluid System Simulation Program (GFSSP) that has been validated against several small-scale liquid nitrogen NVF tests. This paper presents GFSSP numerical model validation against three liquid hydrogen NVF tests on a large-scale propellant tank. Comparing the GFSSP model to the data for pressure, fill level, and wall temperature, the mean absolute percentage error is below 20% across the range of simulated conditions, indicating very good predictive agreement for this large-scale hydrogen NVF test series.

Alok Majumdar↗

Numerical Simulation of No-Vent Chill and Fill of a Large Liquid Hydrogen Tank

To enable the design of future in-space cryogenic propellant transfer vehicles, high accuracy models of the cryogenic propellant transfer process are desired. For in-space refueling applications, accurate simulations are desired of the propellant tank chilldown and no-vent fill (NVF) process in order to estimate the amount of propellant required, to determine the maximum sustainable flow rate to fill the tank, and to establish a timeline for the filling process. NVF is a complex process that involves numerous two-phase heat and mass transfer phenomena, such as boiling at the wall and evaporation and condensation at the liquid/vapor interface. In recent years, progress has been made in developing a numerical modeling technique using the Generalized Fluid System Simulation Program (GFSSP) that has been validated against several small-scale liquid nitrogen NVF tests. This paper presents GFSSP numerical model validation against three liquid hydrogen NVF tests on a large-scale propellant tank. Comparing the GFSSP model to the data for pressure, fill level, and wall temperature, the mean absolute percentage error is below 20% across the range of simulated conditions, indicating very good predictive agreement for this large-scale hydrogen NVF test series.

Alok Majumdar↗

Flow Field Measurement of Mixing Driven by Buoyancy

Mixing driven by buoyancy-induced flows inside a cavity consists of stretching and folding of an interface. Measurement of the flow field using particle imaging velocimetry shows that during stretching the flow field has a single elliptic point, thus dominated by a single vortex. However, global bifurcation that results in folding introduces a hyperbolic point whereby the flow field degenerates to multiple vortex interactions. The short-lived coherent structure observed during mixing which results in the Rayleigh- Taylor morphology is attributed to vortex interactions. The mixing characteristics of non-homogeneous fluids driven by buoyancy are important towards understanding transport phenomenon in a microgravity environment. Mixing consists of stretching and folding of an interface due to a flow field whose intensity depends on the body force. For miscible liquids, the characteristic of the flow field determines whether mass transport is governed by diffusion or bulk stirring which induces mixing. For technologically important processes, transport of mass is governed by the coupling of the body force to scalar gradients such as concentration and or temperature' 2 3 . In order to lend insight into these classes of problems we consider a model experimental system to study mixing driven by buoyancy-induced flows. The characteristics of mixing is addressed from detail measurements of the flow field using particle imaging velocimetry (PIV), and its corresponding interface dynamics using image processing techniques.

Batur, C.↗

Numerical modeling of spray combustion with an advanced VOF method

This paper summarizes the technical development and validation of a multiphase computational fluid dynamics (CFD) numerical method using the volume-of-fluid (VOF) model and a Lagrangian tracking model which can be employed to analyze general multiphase flow problems with free surface mechanism. The gas-liquid interface mass, momentum and energy conservation relationships are modeled by continuum surface mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed flow regimes. The objectives of the present study are to develop and verify the fractional volume-of-fluid cell partitioning approach into a predictor-corrector algorithm and to demonstrate the effectiveness of the present approach by simulating benchmark problems including laminar impinging jets, shear coaxial jet atomization and shear coaxial spray combustion flows.

Chen, Yen-Sen↗

Effects of Traveling Magnetic Field on Dynamics of Solidification

TMF is based on imposing a controlled phase-shift in a train of electromagnets, forming a stack. Thus, the induced magnetic field can be considered to be travelling along the axis of the stack. The coupling of this traveling wave with an electrically conducting fluid results in a basic flow in a form of a single axisymmetric roll. The magnitude and direction of this flow can be remotely controlled. Furthermore, it is possible to localize the effect of this force field though activating only a number of the magnets. This force field generated in the fluid can, in principle, be used to control and modify convection in the molten material. For example, it can be used to enhance convective mixing in the melt, and thereby modify the interface shape, and macrosegregation. Alternatively, it can be used to counteract thermal and/or solutal buoyancy forces. High frequency TMF can be used in containerless processing techniques, such as float zoning, to affect the very edge of the fluid so that Marangoni flow can be counter balanced. The proposed program consists of basic fundamentals and applications. Our goal in conducting the following experiments and analyses is to establish the validity of TMF as a new tool for solidification processes. Due to its low power consumption and simplicity of design, this tool may find wide spread use in a variety of space experiments. The proposed ground based experiments are intended to establish the advantages and limitations of employing this technique. In the fundamentals component of the proposed program, we will use theoretical tools and experiments with mercury to establish the fundamental aspects of TMF-induced convection through a detailed comparison of theoretical predictions and experimental measurements of flow field. In this work, we will conduct a detailed parametric study involving the effects of magnetic field strength, frequency, wave vector, and the fluid geometry. The applications component of this work will be focused on investigating the effect of TMF on the following solidification and pre-directional solidification processes: (1) Bridgman growth of Ga:Ge with the goal of counteracting the buoyancy-driven convection; (2) Mixing of Pb-Ga and Pb-Sn alloys with the aim of initiating and maintaining a uniform melt prior to solidification processing; and (3) Float Zone growth with the aim of identifying, through simulations and model experiments, conditions needed to counteract Marangoni flow in a microgravity environment. The proposed research has strong relevance to microgravity research and the objectives of the NRA. TMF can provide a unique and accurate mechanism for generation and control of desirable flow patterns for microgravity research. These attributes have significant relevance to 1) Alloy mixing prior to solidification in a microgravity environment. TMF can provide this mixing with a low level of power consumption; (2) TMF can offset the deleterious effects of Marangoni convection in microgravity containerless processing. Thus, TMF can be instrumental in further understanding this phenomena; (3) Generation of controlled flows will allow the investigation of the effect of these flows on growth morphology and growth kinetics; and (4) On Earth, TMF has the potential to significantly counter-balance thermosolutal convection, thereby creating conditions similar to those obtained in microgravity. Once demonstrated, this new tool for use in solidification has the strong potential to find applications in a host of microgravity material research projects.

Mazuruk, Konstantin↗

Determination of Interfacial Rheological Properties through Microgravity Oscillations of Bubbles and Drops

This report summarizes our derivations of analytical expressions for the frequencies and damping constants for small-amplitude axisymmetric shape oscillations of a liquid drop suspended in an immiscible fluid host in microgravity. In particular, this work addresses large Reynolds number shape oscillations and focuses on the surface rheological effects that arise from the presence of insoluble surfactants at the interface. Parameters characterizing viscous effects from the bulk phases, surface viscous effects, Marangoni effects from the surface advection and diffusion of surfactants, and the Gibbs elasticity are all considered and analyzed to determine the relative importance of each contribution. Supplementing the analytical treatment for small-amplitude oscillations, a numerical boundary integral equation formulation is developed for the study of large-amplittide axisymmetric oscillations of a drop in vacuum. The boundary integral formulation is an extension of classical potential flow theory and approximately accounts for viscous effects in the bulk fluid as well as the surface viscous and Marangoni effects resulting from an insoluble surfactant contaminating the interface. Theoretical and numerical results are presented for four distinct cases. These, range from the case when the effects of the surfactants are 'negligible' to 'large' when compared to the viscous effects in the bulk phases. The feasibility of the non-contact measurement of the surface parameters, using experimental observations for the oscillation frequencies and damping constants of drops and bubbles, is discussed.

Nadim, Ali↗

The Reservoir Temperature Estimator (RTEst): A multicomponent geothermometry tool

The Reservoir Temperature Estimator (RTEst) is a multicomponent geothermometry tool for estimating reservoir geochemical parameters including reservoir temperature, CO 2 fugacity, mass of water lost or gained, and a reaction factor. It estimates these parameters and their associated uncertainties by minimizing an objective function that is the weighted sum of squares of the saturation indices of a user-selected set of minerals believed to be equilibrated with the reservoir fluid. RTEst accomplishes these estimates by combining the geochemical modeling capabilities of The Geochemist’s Workbench® with the optimization/parameter estimation resources of PEST®. An included interface aids the user in selecting plausible mineral phases to comprise the objective function and calculates their weighting factors. The working principles of RTEst are described and its efficacy is illustrated by presenting results of its application to various geothermal fields with known conditions. These examples show RTEst can account for the alteration of ascending reservoir fluid by mineral (calcite) re-equilibration with changes in temperature, reconstruct waters with CO 2 loss, correct for the deficit of water and other volatiles (CO 2,gas , H 2 S gas ) from boiling, and determine the amount of mixing of thermal and non-thermal waters. RTEst can use data with basis species below detection limit, missing, or unreliable either by assuming equilibrium with a controlling mineral (fixed-analyte method) or by treating the analyte concentration as an optimization parameter. The inverse of variance weighting method included in RTEst provides more representative results than either the normalization or unit weighting methods. Finally, the ability of RTEst to calculate reservoir temperatures, gas fugacity, and mixing fractions demonstrates its usefulness as a tool for evaluating geothermal systems.

15 GEOTHERMAL ENERGY↗

Liquid Hydrogen Tank Chill and No-Vent Fill Prediction using Computational Fluid Dynamics

Cryogenic tank chill and fill is an important cryogenic fluid management (CFM) technology that supports and enables many of NASA’s long-duration space missions. For no-vent fill, the receiver tank pressure remains below the supply tank pressure during the entire duration of the fill so that the tank does not require venting. This is especially advantageous for tank fill operations in low gravity where the position of the liquid is not always known and venting the tank may cause loss of propellant by venting liquid. In lieu of expensive tests conducted on-orbit, accurate computational models capable of predicting receiver tank pressure during cryogenic propellant tank fill may be used to reduce system and propellant mass as well as mission risk. However, these numerical models must be validated or anchored to test data. This study presents a computational fluid dynamics (CFD) model with conjugate heat transfer that is used to predict a liquid hydrogen tank chill and no-vent fill ground test conducted at Lewis Research Center (now Glenn Research Center) in 1991. The specific test case chosen for model validation implemented an upward-facing jet near the bottom of a cylindrical 34 liter tank. CFD predictions are compared to experimental measurements of tank pressure, fill level, fluid temperatures, and wall temperatures. The CFD results show reasonable agreement to the test data but overpredict the pressure collapse near the end of the fill despite the liquid jet penetrating the liquid-vapor interface. Several sensitivity studies are considered due to notable uncertainties in the experiment.

Computational Fluid Dynamics↗

Learning Coarse-Grained Potentials for Binary Fluids

Even though atomistic and CG models have been used to simulate liquid nanodroplets in vapor, very few rigorous studies of the liquid-liquid interface structure are available, and most of them are limited to planar interfaces. In this work, we evaluate several existing force fields (FF)s, including two atomistic and three CG FFs, with respect to modeling the interface structure and thermodynamic properties of the water-hexane interface. Both atomistic FFs are able to quantitatively reproduce the interfacial tension and the coexisting densities of the experimentally-observed planar interface. We use the atomistic FFs to model water droplets in hexane and use these simulations to test the CG FFs. We found that the tested CG FFs cannot reproduce the interfacial tensions of planar and/or curved interfaces. Finally, we propose a new approach for learning coarse-grained (CG) potentials within the CG SDK (Shinoda-DeVane- Klein) FF framework from atomistic simulation data. We demonstrate that the new potential significantly improves the prediction of both the interfacial tension and structure of water-hexane planar and curved interfaces.

Gao, Peiyuan↗

Morphological Stability and Fluid Dynamics of Vapor Crystal Growth

A fundamental understanding of the conditions under which crystals can retain morphological stability, i.e., shape stability of the advancing interface, during growth from vapors was studied. Morphological stability (MS) is a necessary condition for the growth of homogeneous single crystals required for numerous device applications. For crystallization from melts, the MS concepts are well developed and are essentially based on heat and mass transfer conditions about the advancing interface. For crystallization from vapors, the MS requirements are more complex and not well understood. The added complexity arises from the fact that anisotropies in interfacial kinetics are typically stronger in crystallization from vapors than from melts. These pronounced anisotropies root in the distinctly lower atomic roughness of most vapor-solid interfaces.

Rosenberger, F.↗

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.↗

Growth experiment of narrow band-gap semiconductor PbSnTe single crystals in space (M-1)

An experiment on crystal growth of Pb(1-x)Sn(x)Te in microgravity is planned. This material is an alloy of the compound semiconductors PbTe and SnTe. It is a promising material for infrared diode lasers and detectors in the wavelength region between 6 and 30 micron. Since the electrical properties of Pb(1-x)Sn(x)Te depend greatly on the Pb/Sn ratio and crystalline defects as well as impurity concentration, homogeneous, defect-free, high-quality crystals are anticipated. Although many growth methods, such as the pulling method, the Bridgman method, the vapor growth method, etc., have been applied to the growth of Pb(1-x)Sn(x)Te, large, homogeneous, low-defect-density crystals have not yet been grown on Earth. The unsuccessful results were caused by buoyancy-driven convection in the fluids induced by the specific gravity difference between heated and cooled fluids on Earth. A crystal is grown by cooling the melt from one end of the ampoule. In crystal growth from the melt, about 30 percent of the SnTe in the melt is rejected at the solid-liquid interface during solidification. On Earth, the rejected SnTe is completely mixed with the remaining melt by convection in the melt. Therefore, SnTe concentration in the melt, and accordingly in the crystal, increases as the crystal grows. In the microgravity environment, buoyancy-driven convection is suppressed because the specific gravity difference is negligible. In that case, the rejected SnTe remains at the solid-liquid interface and its concentration increases only at the interface. If the growth rate is higher than the PbTe-SnTe interdiffusion rate, the amount of SnTe which diffuses from the interface into the melt increases as SnTe piles up at the interface, and finally it balances the amount of rejected SnTe during solidification, resulting in steady-state SnTe transportation at the interface. By using this principle, compositionally homogeneous crystals can be grown. Furthermore, low-defect-density crystals will be grown in microgravity, because convection causes crystalline defects by mising hot and cold fluids and generating temperature fluctuations in them.

Yamada, Tomoaki↗

CryoTran user's manual, version 1.0

The development of cryogenic fluid management systems for space operation is a major portion of the efforts of the Cryogenic Fluids Technology Office (CFTO) at the NASA Lewis Research Center. Analytical models are a necessary part of experimental programs which are used to verify the results of experiments and are also used as a predictor for parametric studies. The CryoTran computer program is a bridge to obtain analytical results. The object of CryoTran is to coordinate these separate analyses into an integrated framework with a user-friendly interface and a common cryogenic property database. CryoTran is an integrated software system designed to help solve a diverse set of problems involving cryogenic fluid storage and transfer in both ground and low-g environments.

Cowgill, Glenn R.↗

Growth and Morphology of Supercritical Fluids, a Fluid Physics Experiment Conducted on Mir, Complete

The Growth and Morphology of Supercritical Fluids (GMSF) is an international experiment facilitated by the NASA Glenn Research Center and under the guidance of U.S. principal investor Professor Hegseth of the University of New Orleans and three French coinvestigators: Daniel Beysens, Yves Garrabos, and Carole Chabot. The GMSF experiments were concluded in early 1999 on the Russian space station Mir. The experiments spanned the three science themes of near-critical phase separation rates, interface dynamics in near-critical boiling, and measurement of the spectrum of density fluctuation length scales very close to the critical point. The fluids used were pure CO2 or SF6. Three of the five thermostats used could adjust the sample volume with the scheduled crew time. Such a volume adjustment enabled variable sample densities around the critical density as well as pressure steps (as distinct from the usual temperature steps) applied to the sample. The French-built ALICE II facility was used for these experiments. It allows tightly thermostated (left photograph) samples (right photograph) to be controlled and viewed/measured. Its diagnostics include interferometry, shadowgraph, high-speed pressure measurements, and microscopy. Data were logged on DAT tapes, and PCMCIA cards and were returned to Earth only after the mission was over. The ground-breaking near critical boiling experiment has yielded the most results with a paper published in Physical Review Letters (ref. 1). The boiling work also received press in Science Magazine (ref. 2). This work showed that, in very compressible near-critical two-phase pure fluids, a vapor bubble was induced to temporarily overheat during a rapid heating of the sample wall. The temperature rise in the vapor was 23-percent higher than the rise in the driving container wall. The effect is due to adiabatic compression of the vapor bubble by the rapid expansion of fluid near the boundary during heatup. Thermal diffusivity is low near the critical point, so getting heat out of the compressed bubble is observably slow. This gives the appearance of a backward heat flow, or heat flow from a cold surface to a warm fluid.

Wilkinson, R. Allen↗

Fluid System Angular Momentum and Moment Fluctuations Driven by Gravity Gradient or Gravity Jitter in Microgravity

The dynamical behaviour of fluids affected by the asymmetric gravity gradient acceleration and gravity jitter acceleration, in particular the effect of surface tension on partially-filled rotating fluids with special example applicable to a full-scale Gravity Probe-B Spacecraft dewar tank have been investigated. Results of slosh wave excitation along liquid-vapor interface induced by gravity gradient acceleration indicated that one-up one-down and one-down one-up oscillations of two bubbles in the cross-section of doughnut profiles in the vertical r-z plane of rotating dewar, and an eccentric contour of bubble rotating around the axis of dewar in horizontal r-Theta plane. Results of slosh wave excitation along liquid-vapor interface induced by gravity jitter acceleration indicated equivalent to time-dependent oscillatory forces which push bubble in the combined bubble motion of down-and-up and leftward-and-rightward as the bubble is rotating with respect to rotating dewar axis. Fluctuations of angular momentum, fluid moment and bubble mass center caused by slosh wave excitations driven by gravity gradient acceleration or gravity jitter acceleration are also investigated.

Hung, R. J.↗

An arbitrarily high-order three-dimensional Cartesian-grid method for reconstructing interfaces from volume fraction fields

Here Tthis work describes a newly developed, arbitrarily high-order Cartesian-grid method for reconstructing material interfaces from a volume fraction field. The method begins by identifying all of the grid cells in the volume fraction field that are intersected by the interface and need to be approximated by the reconstruction scheme. Finite-differences are used to calculate the gradient of the volume fraction field and provide an estimate of the surface normal in all of the interfacial grid cells. Groups of connected grid cells are then identified which all have the same dominant component of the normal vector. This grouping by orientation determines the proper dependent variable to use in the surface reconstruction (e.g. for a 2D curve, this step determines if the surface will be approximated by a function of x or y). A cumulative integral over the surface is constructed and fit using b-splines for two-dimensional problems or tensor-product b-splines for three-dimensional problems. This construction allows for the interface to be recovered through application of the second fundamental theorem of calculus. Fitting the cumulative integral with $\mathscr{N}$ th-order b-splines (or tensor-product b-splines) yields an ($\mathscr{N}$-1) th-order convergence rate of the interface shape. Differentiation of the b-spline interface function(s) allows for the high-order approximation of the normal vector and curvature to be obtained directly anywhere along b-spline. Together, the proposed reconstruction technique can achieve arbitrarily high mesh convergence rates. Validation tests are presented with mesh convergence rates ranging from fourth- to tenth-order.

97 MATHEMATICS AND COMPUTING↗

An Interfacial Engineering Approach toward Operation of a Porous Solid Electrolyte CO 2 Electrolyzer

Waste CO 2 can be repurposed as a carbon feedstock for synthesizing valuable chemicals via CO 2 electrolysis. Porous solid electrolyte (PSE) CO 2 electrolysis has been demonstrated as an economically viable method to produce high purity products. This work applies an interfacial engineering approach to determine key factors to improve performance in PSE CO 2 electrolyzers. We standardize the assembly by binding the ionic resin into an ionomer wafer and utilize Computational Fluid Dynamics (CFD) to design gaskets for uniform fluid flow. Here, we employ the distribution of relaxation times (DRT) method to determine that anionic-conducting interfaces are the primary contributor to energy losses. To address this, we demonstrate that enhancing the contact between the cathode and the anion exchange membrane (AEM) and the AEM-ionic resin interface allows for low overpotential in deionized water operation.

09 BIOMASS FUELS↗

Raman scattering applied to gas mixing in the successive channel flow of two gases

The Raman scattering technique has been applied to investigate the effects of fluid undercutting in a simulation of the fluid-piston concept for wind-tunnel design. Measurements of the structure of the driver-driven gas interface and of the effect of undercutting on test duration were performed. Tests were conducted over a Froude number range of 0.0021 to 0.526 at driver-to-driven gas density ratios of 1.103 and 1.573. The effects of inlet design, test-channel inclination angle, and pretest pressurization were also examined. The results are compared with shadowgraph and pressure measurements to illustrate the characteristics of the Raman technique and to explain the discrepancy between the shadowgraph and pressure results.

Hillard, M. E., Jr.↗