Search NASA⌕ Search

SEARCH · Search NASA

Results for “complex fluid”

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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Understanding the breakdown of classic two-phase theory and spray atomization at engine-relevant conditions

We present a generalized framework for multi-component liquid injections to understand and predict the breakdown of classic two-phase theory and spray atomization at engine-relevant conditions. The analysis focuses on the thermodynamic structure and the immiscibility state of representative gas-liquid interfaces. The most modern form of Helmholtz energy mixture state equation is utilized which exhibits a unique and physically-consistent behavior over the entire two-phase regime of fluid densities. It is combined with generalized models for non-linear Gradient Theory and for liquid injections to quantify multi-component two-phase interface structures in global thermal equilibrium. Then, the Helmholtz free energy is minimized which determines the interfacial species distribution as a consequence. This minimal free energy state is demonstrated to validate the underlying assumptions of classic two-phase theory and spray atomization. However, under certain engine-relevant conditions for which corroborating experimental data is presented, this requirement for interfacial thermal equilibrium becomes unsustainable. A rigorously derived probability density function quantifies the ability of the interface to develop internal spatial temperature gradients in the presence of significant temperature differences between injected liquid and ambient gas. Then, the interface can no longer be viewed as an isolated system at minimal free energy. Instead, the interfacial dynamics become intimately connected to those of the separated homogeneous phases. Hence, the interface transitions toward a state in local equilibrium whereupon it becomes a dense-fluid mixing layer. A new conceptual view of a transitional liquid injection process emerges from a transition time scale analysis. Close to the nozzle exit, the two-phase interface still remains largely intact and more classic two-phase processes prevail as a consequence. Further downstream, however, the transition to dense- fluid mixing generally occurs before the liquid length is reached. The significance of the presented modeling expressions is established by a direct comparison to a reduced model, which utilizes widely-applied approximations but fundamentally fails to capture the physical complexity discussed in this paper.

42 ENGINEERING↗
DOE OSTI2016

A Study of Physics-Informed Deep Learning for System Fluid Dynamics Closures

The sub-grid-scale (SGS) physics models, or so-called closure relations (CR), are essential in thermal-fluid modeling and simulation codes, ranging from large-eddy simulation to system-level thermo-fluid dynamic problems. While details vary, all thermo-fluid codes solve coupled mass-momentum-energy conservation partial differential equations (PDE). Traditionally, the analytical form of closure models has to be known, and then we calibrate the model using data from a set of applicable experiments. Characteristically, it takes long effort to gain insights and mechanistic understanding needed for developing closure models. Such a model is a compact form of data. The long time required for new model development limits the pace of model application while dealing with new geometries, new coolants or system conditions, particularly in newly designed systems. Alternative to mechanistic and semi-analytical models, it is possible to discover the pattern of closure relations by applying the modern machine learning methods to the body of relevant experimental data. Deep learning (DL) is a relatively new branch in the machine learning family, and a major breakthrough has been achieved since Hinton introduces a fast algorithm to efficiently train a feed-forward neural network (NN) with multilayer perceptrons. Generally, any NN with more than two layers is considered as deep learning, and it contains a large number of hyper-parameters to achieve pattern recognition or regression for complex datasets. Hornik shows that multilayer NN is compatible with universal approximation theorem. There is no constraint for a multilayer NN to capture the properties of any measurable information. This capability makes DL attractive for closure models in a system thermal-hydraulic code. As system geometry and/or operating conditions change, new physical experiments (increasingly sophisticated) are conducted to support evaluation and revision of closure models. Such a so-called data-driven modeling relaxes the structural form of system codes to directly assimilate data instead of hard-wiring the model to pre-calibrated correlations. The objective of this research is to enable cost-effective development of a DL-based closure that reflects the data. The previous work on applying machine learning techniques to fluid dynamics focused on statistical learning with basic NN from direct numerical simulations (DNS) dataset. In the present work, we show that DL-based closure models can cause the fluid calculation to diverge if multilayer NN is used. To address the issue, we investigate conditions by which the DL-based closure models work compatibly, stably, and effectively with PDE-constrained forward prediction problems. Optimally meeting the compatibility conditions constitutes a physics-informed deep learning (PIDL) strategy. We refer to this process as PDE-constrained machine learning or PDE-ML. To develop the technical basis for PDE-ML using DL-based fluid dynamics closure, numerical experiments are performed for manufactured problems of increasing complexity. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Induced seismicity in geologic carbon storage

Geologic carbon storage, as well as other geo-energy applications, such as geothermal energy, seasonal natural gas storage and subsurface energy storage imply fluid injection and/or extraction that causes changes in rock stress field and may induce (micro)seismicity. If felt, seismicity has a negative effect on public perception and may jeopardize wellbore stability and damage infrastructure. Thus, induced earthquakes should be minimized to successfully deploy geo-energies. However, numerous processes may trigger induced seismicity, which contribute to making it complex and translates into a limited forecast ability of current predictive models. We review the triggering mechanisms of induced seismicity. Specifically, we analyze (1) the impact of pore pressure evolution and the effect that properties of the injected fluid have on fracture and/or fault stability; (2) non-isothermal effects caused by the fact that the injected fluid usually reaches the injection formation at a lower temperature than that of the rock, inducing rock contraction, thermal stress reduction and stress redistribution around the cooled region; (3) local stress changes induced when low-permeability faults cross the injection formation, which may reduce their stability and eventually cause fault reactivation; (4) stress transfer caused by seismic or aseismic slip; and (5) geochemical effects, which may be especially relevant in carbonate-containing formations. We also review characterization techniques developed by the authors to reduce the uncertainty in rock properties and subsurface heterogeneity both for the screening of injection sites and for the operation of projects. Based on the review, we propose a methodology based on proper site characterization, monitoring and pressure management to minimize induced seismicity.

58 GEOSCIENCES↗
DOE OSTI2015

HYDRATE v1.5 OPTION OF TOUGH+ v1.5

HYDRATE v1.5 is a numerical code that for the simulation of the behavior of hydrate-bearing geologic systems, and represents the third update of the code since its first release [Moridis et al., 2008]. It is an option of TOUGH+ v1.5 [Moridis and Pruess, 2014], a successor to the TOUGH2 [Pruess et al., 1999, 2012] family of codes for multi-component, multiphase fluid and heat flow developed at the Lawrence Berkeley National Laboratory. HYDRATE v1.5 needs the TOUGH+ v1.5 core code in order to compile and execute. It is written in standard FORTRAN 95/2003, and can be run on any computational platform (workstation, PC, Macintosh) for which such compilers are available. By solving the coupled equations of mass and heat balance, the fully operational TOUGH+HYDRATE code can model the non-isothermal gas release, phase behavior and flow of fluids and heat under conditions typical of common natural CH4-hydrate deposits (i.e., in the permafrost and in deep ocean sediments) in complex geological media at any scale (from laboratory to reservoir) at which Darcy's law is valid. TOUGH+HYDRATE v1.5 includes both an equilibrium and a kinetic model of hydrate formation and dissociation. The model accounts for heat and up to four mass components, i.e., water, CH4, hydrate, and water-soluble inhibitors such as salts or alcohols. These are partitioned among four possible phases (gas phase, liquid phase, ice phase and hydrate phase). Hydrate dissociation or formation, phase changes and the corresponding thermal effects are fully described, as are the effects of inhibitors. The model can describe all possible hydrate dissociation mechanisms, i.e., depressurization, thermal stimulation, salting-out effects and inhibitor-induced effects.

Moridis, GeorgeJ↗
DOE OSTI2015

RealGasBrine v1.0 option of TOUGH+ v1.5

RealGasBrine v1.0 is a numerical code that for the simulation of the behavior of gas-bearing porous and/fractured geologic media. It is an option of TOUGH+ v1.5 [Moridis, 2014], a successor to the TOUGH2 [Pruess et al., 1999; 2012] family of codes for multi-component, multiphase fluid and heat flow developed at the Lawrence Berkeley National Laboratory. RealGasBrine v1.0 needs the TOUGH+ v1.5 core code in order to compile and execute. It is written in standard FORTRAN 95/2003, and can be run on any computational platform (workstation, PC, Macintosh) for which such compilers are available. RealGasBrine v1.0 describes the non-isothermal two- (for pure water) or three-phase (for brine) flow of an aqueous phase and a real gas mixture in a gas-bearing medium, with a particular focus in ultra-tight (such as tight-sand and shale gas) systems. Up to 12 individual real gases can be tracked, and salt can precipitate as solid halite. The capabilities of the code include coupled flow and thermal effects, real gas behavior, Darcy and non-Darcy flow, several isotherm options of gas sorption onto the grains of the porous media, complex fracture descriptions, gas solubility into water, and geomechanical effects on flow properties. RealGasBrine v1.0 allows the study of flow and transport of fluids and heat over a wide range of time frames and spatial scales not only in gas reservoirs, but also in any problem involving the flow of gases in geologic media, including the geologic storage of greenhouse gas mixtures, the behavior of geothermal reservoirs with multi-component condensable (H2O and CO2) and non-condensable gas mixtures, the transport of water and released H2 in nuclear waste storage applications, etc.

Moridis, GeorgeJ↗
DOE OSTI2024

Vidyut3d: A Non-Equilibrium Plasma Modeling Tool [SWR-24-101]

Vidyut3d is a massively-parallel plasma-fluid solver for low-temperature plasmas (LTPs) that supports both local field (LFA) and local mean energy (LMEA) approximations, as well as complex gas and surface-phase chemistry. The solver supports 2D and 3D domains, and uses AMReX's adaptive mesh refinement capabilities to increase the grid resolution around complex structures (e.g. streamer heads and sheaths) while maintaining a tractable problem size. Vidyut specializes in simulating various types of gas-phase discharges, as well as plasma/surface interactions and surface chemistry (e.g. for plasma-mediated catalysis applications). The solver also supports hybrid CPU/GPU parallelization strategies, and has demonstrated excellent scaling on various HPC architectures for problem sizes consisting of O(100 M) control volumes.

Sitaraman, Hariswaran↗
DOE OSTI2014

User's Manual for the Hydrate v1.5 Option of TOUGH+ v1.5: A Code for the Simulation of System Behavior in Hydrate-Bearing Geologic Media

HYDRATE v1.5 is a numerical code that for the simulation of the behavior of hydrate-bearing geologic systems, and represents the third update of the code since its first release [Moridis et al., 2008]. It is an option of TOUGH+ v1.5 [Moridis, 2014], a successor to the TOUGH2 [Pruess et al., 1999] family of codes for multi-component, multiphase fluid and heat flow developed at the Lawrence Berkeley National Laboratory. HYDRATE v1.5 needs the TOUGH+ v1.5 core code in order to compile and execute. It is written in standard FORTRAN 95/2003, and can be run on any computational platform (workstation, PC, Macintosh) for which such compilers are available. By solving the coupled equations of mass and heat balance, the fully operational TOUGH+HYDRATE code can model the non-isothermal gas release, phase behavior and flow of fluids and heat under conditions typical of common natural CH4-hydrate deposits (i.e., in the permafrost and in deep ocean sediments) in complex geological media at any scale (from laboratory to reservoir) at which Darcy’s law is valid. TOUGH+HYDRATE v1.5 includes both an equilibrium and a kinetic model of hydrate formation and dissociation. The model accounts for heat and up to four mass components, i.e., water, CH4, hydrate, and water-soluble inhibitors such as salts or alcohols. These are partitioned among four possible phases (gas phase, liquid phase, ice phase and hydrate phase). Hydrate dissociation or formation, phase changes and the corresponding thermal effects are fully described, as are the effects of inhibitors. The model can describe all possible hydrate dissociation mechanisms, i.e., depressurization, thermal stimulation, salting-out effects and inhibitor-induced effects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Modeling of Chemical Slurry Rheology in DWPF Sludge Batch (SB) 10 Simulants

The Defense Waste Processing Facility (DWPF) treats high-activity radionuclides from sludge through a process called vitrification. This process converts radioactive liquid waste currently stored in tank farms into a solid glass form that is suitable for long-term storage and disposal. Due to the complexities involved in vitrifying this waste within each operation of the Chemical Processing Cell (CPC), waste rheology is studied to characterize the fluid-mechanical properties as it passes through the CPC and into the Melter. To better understand the waste and validate flow behavior, slurry rheology of simulants that represents the waste was studied at various acid stoichiometry percentages and solids concentrations to determine the simulant’s yield stress and viscosity. This research work has been supported by the DOE-FIU Science & Technology Workforce Development Initiative, an innovative program developed by the U.S. Department of Energy’s Office of Environmental Management (DOE-EM) and Florida International University’s Applied Research Center (FIU-ARC). During the spring of 2022, a DOE Fellow intern, Brendon Cintas, spent 10 weeks doing a summer internship at Savannah River National Laboratory (SRS) under the supervision and guidance of Dan Lambert, Chemical Flowsheet Development. The intern’s project was initiated on June 6, 2022, and continued through August 11, 2022 with the objective of assisting scientists at SRNL’s Rheology and Grout Laboratory at Aiken Country Technology Lab (ACTL) better understand the sludge composition on the rheology of a simulant slurry using a HAAKE RheoStress 6000 rheometer and extrapolate the results to the real-waste data.

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

Gas Diffusion in Fluids Containing Bubbles

Mathematical model describes movement of gases in fluid containing many bubbles. Model makes it possible to predict growth and shrink age of bubbles as function of time. New model overcomes complexities involved in analysis of varying conditions by making two simplifying assumptions. It treats bubbles as point sources, and it employs approximate expression for gas concentration gradient at liquid/bubble interface. In particular, it is expected to help in developing processes for production of high-quality optical glasses in space.

Zak, M.↗
NASA NTRS1993

Convection equation modeling: A non-iterative direct matrix solution algorithm for use with SINDA

The determination of the boundary conditions for a component-level analysis, applying discrete finite element and finite difference modeling techniques often requires an analysis of complex coupled phenomenon that cannot be described algebraically. For example, an analysis of the temperature field of a coldplate surface with an integral fluid loop requires a solution to the parabolic heat equation and also requires the boundary conditions that describe the local fluid temperature. However, the local fluid temperature is described by a convection equation that can only be solved with the knowledge of the locally-coupled coldplate temperatures. Generally speaking, it is not computationally efficient, and sometimes, not even possible to perform a direct, coupled phenomenon analysis of the component-level and boundary condition models within a single analysis code. An alternative is to perform a disjoint analysis, but transmit the necessary information between models during the simulation to provide an indirect coupling. For this approach to be effective, the component-level model retains full detail while the boundary condition model is simplified to provide a fast, first-order prediction of the phenomenon in question. Specifically for the present study, the coldplate structure is analyzed with a discrete, numerical model (SINDA) while the fluid loop convection equation is analyzed with a discrete, analytical model (direct matrix solution). This indirect coupling allows a satisfactory prediction of the boundary condition, while not subjugating the overall computational efficiency of the component-level analysis. In the present study a discussion of the complete analysis of the derivation and direct matrix solution algorithm of the convection equation is presented. Discretization is analyzed and discussed to extend of solution accuracy, stability and computation speed. Case studies considering a pulsed and harmonic inlet disturbance to the fluid loop are analyzed to assist in the discussion of numerical dissipation and accuracy. In addition, the issues of code melding or integration with standard class solvers such as SINDA are discussed to advise the user of the potential problems to be encountered.

Schrage, Dean S.↗
NASA NTRS1996

A New Methodology for Turbulence Modelers Using DNS Database Analysis

Many industrial applications in such fields as aeronautical, mechanical, thermal, and environmental engineering involve complex turbulent flows containing global separations and subsequent reattachment zones. Accurate prediction of this phenomena is very important because separations influence the whole fluid flow and may have an even bigger impact on surface heat transfer. In particular, reattaching flows are known to be responsible for large local variations of the local wall heat transfer coefficient as well as modifying the overall heat transfer. For incompressible, non-buoyant situations, the fluid mechanics have to be accurately predicted in order to have a good resolution of the temperature field.

Parneix, S.↗
NASA NTRS1999

Three-Dimensional Imaging

This chapter is concerned with three-dimensional imaging of fluid flows. Although relatively young, this field of research has already yielded an enormous range of techniques. These vary widely in cost and complexity, with the cheapest light sheet systems being within the budgets of most laboratories, and the most expensive Magnetic Resonance Imaging systems available to a select few. Taking the view that the most likely systems to be developed are those using light sheets, the authors will relate their knowledge and experience of such systems. Other systems will be described briefly and references provided. Flows are inherently three-dimensional in structure; even those generated around nominally 2-D surface geometry. It is becoming increasingly apparent to scientists and engineers that the three-dimensionalities, both large and small scale, are important in terms of overall flow structure and species, momentum, and energy transport. Furthermore, we are accustomed to seeing the world in three dimensions, so it is natural that we should wish to view, measure and interpret flows in three-dimensions. Unfortunately, 3-D images do not lend themselves to convenient presentation on the printed page, and this task is one of the challenges facing us.

Kelso, R. M.↗
NASA NTRS2002

Computational Systems for Multidisciplinary Applications

In this paper, we briefly describe our efforts to develop complex simulation systems. We focus first on four key infrastructure items: enterprise computational services, simulation synthesis, geometry modeling and mesh generation, and a fluid flow solver for arbitrary meshes. We conclude by presenting three diverse applications developed using these technologies.

Soni, Bharat↗
NASA NTRS1978

Pipe Line Control

The array of tanks, pipes and valves in the photo below is a petroleum tank farm in Georgia, part of a petrochemical pipe line system that moves refined petroleum products from Texas and Louisiana to the mid-Eastern seaboard. The same pipes handle a number of different products, such as gasoline, kerosene, jet fuel or fuel oil. The fluids are temporarily stored in tanks, pumped into the pipes in turn and routed to other way stations along the pipe line. The complex job of controlling, measuring and monitoring fuel flow is accomplished automatically by a computerized control and communications system which incorporates multiple space technologies.

Source record↗
NASA NTRSDate not supplied

Parametric Analysis of the Charge-Hold-Vent Method for Cryogenic Propellant Tank Chill Down

In the absence of external heat exchangers, the on-orbit transfer of cryogenic propellants requires the receiver tank to first be quenched to a sufficiently low energy state to allow for a continuous no-vent fill to avoid unnecessary venting of liquid. One proposed method for tank chilldown that minimizes the potential for venting liquid is the charge hold vent (CHV) method. CHV follows a cyclic process that gradually removes thermal energy from the receiver tank by injecting liquid with the vent valve closed and allowing the fluid and wall to reach near-thermal equilibrium before venting the superheated vapor. However, the CHV method must be optimized to minimize complexity, mass, and time. This paper presents a modular CHV analytical model used to quantify the number of cycles and propellant mass consumed based on first principles. The model is used to examine the effect of eight parameters: receiver tank material, volume, mass, maximum expected operating pressure, and initial pressure, liquid injection pressure and temperature, and the target temperature. The model is validated against the only two available CHV datasets. Based on results, the tank mass-to-volume ratio is the most important factor in determining the number of CHV cycles and thus degree of difficulty in tank chilldown. The model can easily be used for early-stage design, sizing, and analysis of cryogenic propellant transfer systems.

Tank Chilldown↗
NASA NTRSDate not supplied

Parametric Analysis of the Charge-Hold-Vent Method for Cryogenic Propellant Tank Chill Down

Abstract. In the absence of external heat exchangers, the on-orbit transfer of cryogenic propellants requires the receiver tank to first be quenched to a sufficiently low energy state to allow for a continuous no-vent fill to avoid unnecessary venting of liquid. One proposed method for tank chilldown that minimizes the potential for venting liquid is the charge hold vent (CHV) method. CHV follows a cyclic process that gradually removes thermal energy from the receiver tank by injecting liquid with the vent valve closed and allowing the fluid and wall to reach near-thermal equilibrium before venting the superheated vapor. However, the CHV method must be optimized to minimize complexity, mass, and time. This paper presents a modular CHV analytical model used to quantify the number of cycles and propellant mass consumed based on first principles. The model is used to examine the effect of eight parameters: receiver tank material, volume, mass, maximum expected operating pressure, and initial pressure, liquid injection pressure and temperature, and the target temperature. The model is validated against the only two available CHV datasets. Based on results, the tank mass-to-volume ratio is the most important factor in determining the number of CHV cycles and thus degree of difficulty in tank chilldown. The model can easily be used for early-stage design, sizing, and analysis of cryogenic propellant transfer systems. Keywords: tank chilldown, charge-hold-vent, no-vent fill

Tank Chilldown↗
NASA NTRS2024

A Coordinated Mineralogical and Oxygen Isotopic Study of Refractory Ca-Al-Rich Inclusions

Refractory Ca-Al-rich inclusions (CAIs) in primitive chondritic meteorites represent the oldest solids that formed in our Solar System. Coordinated analyses using multiple ion and electron-beam techniques are vital to fully understand complex nebular and asteroidal processes that CAIs may have experienced in the context of the formation and early evolution of our Solar System. We report the results of a coordinated mineralogical and oxygen isotopic study of unique grossite bearing CAIs from reduced CV (Vigarano type) carbonaceous chondrites using transmission electron microscopy and secondary ion mass spectrometry. Grossite (CaAl 4 O 7 ) occurs in cores of nodules in fine grained CAIs that likely represent aggregates of nebular condensates. In many occurrences, grossite has been partially replaced by hercynite [(Fe,Mg,Zn)Al 2 O 4 ], which displays complex microstructures and compositions, and magnetite nanoparticles. This grossite alteration likely occurred via crystallographically controlled metasomatic reactions with Zn, Fe-rich aqueous fluids on the parent asteroid. The constituent phases of grossite bearing CAIs show heterogeneous oxygen isotopic compositions . That is, grossite and perovskite are 16 O poor and melilite is variably 16 O-depleted in contrast to solar like, 16 O-rich hibonite and spinel. The observed oxygen isotopic distribution is interpreted as a result of mineralogically controlled oxygen isotopic exchange with an 16 O-poor fluid on the parent asteroid. Collectively, we conclude that the presence of limited fluids played an important role in preferential alteration of grossite to hercynite and magnetite and various degrees of 16 O depletion in grossite, perovskite, and melilite during thermal metamorphism on the CV chondrite parent asteroid.

Jangmi Han↗
DOE OSTI2021

A Study on Computational Cost Reduction of Simulations of Phase-Change Material (PCM) Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.

25 ENERGY STORAGE↗