Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Critical Fluids”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Interfacial tension and vapor-liquid equilibria in the critical region of mixtures

In the critical region, the concept of two-scale-factor universality can be used to accurately predict the surface tension between near-critical vapor and liquid phases from the singularity in the thermodynamic properties of the bulk fluid. In the present work, this idea is generalized to binary mixtures and is illustrated using the data of Hsu et al. (1985) for CO2 + n-butane. The pressure-temperature-composition-density data for coexisting, near-critical phases of the mixtures are fitted with a thermodynamic potential comprised of a sum of a singular term and nonsingular terms. The nonuniversal amplitudes characterizing the singular term for the mixtures are obtained from the amplitudes for the pure components by interpolation in a space of thermodynamic 'field' variables. The interfacial tensions predicted for the mixtures from the singular term are within 10 percent of the data on three isotherms in the pressure range (Pc - P)/Pc of less than 0.5. This difference is comparable to the combined experimental and model errors.

Moldover, Michael R.↗
DOE OSTI2025

A quantitative risk assessment framework for fault reactivation in underground hydrogen storage: Coupled simulation and deep learning approach

Underground hydrogen storage (UHS) is emerging as a critical solution for large-scale energy storage. However, like all subsurface fluid injection activities, UHS poses the risk of injection-induced fault reactivation. Accurate risk assessment is essential to ensuring the safety and efficiency of UHS operations. This study presents the development of deep-learning surrogate models for fault reactivation prediction in UHS, trained on a comprehensive database of fully coupled fluid flow-geomechanics simulations. Our findings reveal that analytical models often yield unreliable estimates, with errors up to 54% in the allowable injection pressure, potentially leading to a 40% reduction in UHS operational capacity. The developed surrogate models were incorporated into a quantitative risk assessment (QRA) framework, enabling probabilistic evaluation of fault reactivation risk while accounting for uncertainties in the input variables. Site-specific features, such as horizontal stress gradients, fault’s dip and strike angles, and operational parameters like bottom-hole injection pressure and well-fault distance, were identified as the primary drivers of fault reactivation across various stress regimes. Whereas other hydraulic, geological, and poroelastic reservoir properties were found to have a secondary impact. Notably, we observed that the risk of fault reactivation for a critically oriented fault with a static friction coefficient greater than 0.55 remains below 10% in a normal faulting stress regime. However, the risk significantly increases as the stress regime transitions from normal to strike-slip and ultimately to reverse faulting conditions. These findings underscore the importance of rigorous site characterization and comprehensive QRA evaluations to optimize UHS performance and minimize geomechanical risks.

25 ENERGY STORAGE↗
NASA NTRSDate not supplied

Modeling of Cryogenic Heated-Tube Flow Boiling Experiments of Hydrogen and Helium With GFSSP

Accurate modeling of cryogenic boiling heat transfer is vital for the development of extended-duration space missions. Such missions may require the transfer of cryogenic propellants from in-space storage depots or the cooling of nuclear reactors. Purdue University in collaboration with NASA has assembled a database of cryogenic flow boiling data points from heated-tube experiments dating back to 1959, which has been used to develop new flow boiling correlations specifically for cryogens. Computational models of several of these experiments have been constructed in the Generalized Fluid System Simulation Program (GFSSP), a network flow code developed at NASA’s Marshall Space Flight Center. The new Purdue-developed correlations cover the full boiling curve: onset of nucleate boiling, nucleate boiling, critical heat flux, and film boiling. These correlations have been coded into GFSSP user subroutines. The fluids modeled are hydrogen and helium. Predictions of wall temperature, heat transfer coefficient, and pressure drop are presented and compared to the test data.

multi-phase flow↗
NASA NTRSDate not supplied

Modeling of Cryogenic Heated-Tube Flow Boiling Experiments of Nitrogen and Methane with the Generalized Fluid System Simulation Program

Accurate modeling of cryogenic boiling heat transfer is vital for the development of extended-duration space missions. Such missions may require the transfer of cryogenic propellants from in-space storage depots or the cooling of nuclear reactors. Purdue University in collaboration with NASA has assembled a database of cryogenic flow boiling data points from heated-tube experiments dating back to 1959, which has been used to develop new flow boiling correlations specifically for cryogens. Computational models of several of these experiments have been constructed in the Generalized Fluid System Simulation Program (GFSSP), a network flow code developed at NASA’s Marshall Space Flight Center. The new Purdue-developed correlations cover the full boiling curve: onset of nucleate boiling, nucleate boiling, critical heat flux, and film boiling. These correlations have been coded into GFSSP user subroutines. The fluids modeled are methane and nitrogen. Predictions of wall temperature and heat transfer coefficient are presented and compared to the test data.

Multi-Phase Flow↗
DOE OSTI2021

Computational Fluid Dynamics Simulation of Compressible Non-Newtonian Biomass in a Compression-Screw Feeder

Compression-screw feeders play a critical role in biorefineries to transport lignocellulosic biomass feedstocks from the feeding hoppers to biomass-conversion reactors in order to pretreat and convert biomass to hydrocarbon liquid biofuels and other power and energy resources. One of the main challenges in the operation of screw feeder is plugging and jamming of compressed biomass with high concentration of insoluble solids. The focus of this paper is to numerically investigate the screw feeder at these challenging operating conditions and help with the optimization of the screw feeder design to avoid operation failure. In this work a customized CFD model based on open-source OpenFOAM package [1] was developed to simulate the concentrated biomass as a highly viscous non-Newtonian fluid in the screw feeder. The biomass is modeled as a single-phase compressible Bingham fluid with a plastic viscosity as well as a density-dependent yield stress. The compressibility formulation (pressure-dependent density) and the density-dependent yield stress formulation in the governing equations follow the suggestions from a recent study by Duncan et al. [2]. A pilot-scale hopper/screw feeding system at NREL [3] is used to compare the experimental observations with our simulation results. The auger is 280 mm long and tapered with outer diameter changing from 80 mm to 35 mm. The auger rotates from 10 to 60 rpm in a conical throat which contains anti-rotational bars. The simulations predicted the required torque for the screw feeder and the pressure increase at the exit for biomass feedstocks with various fluid viscosity properties and auger rotating speeds. The analysis of the stress forces helped to identify the critical conditions were the screw feeder excessive wear or jamming could occur.

biofuels↗
NASA NTRSDate not supplied

Modeling of Cryogenic Heated-Tube Flow Boiling Experiments of Hydrogen and Helium With the Generalized Fluid System Simulation Program

Accurate modeling of cryogenic boiling heat transfer is vital for the development of extended-duration space missions. Such missions may require the transfer of cryogenic propellants from in-space storage depots or the cooling of nuclear reactors. Purdue University in collaboration with NASA has assembled a database of cryogenic flow boiling data points from heated-tube experiments dating back to 1959, which has been used to develop new flow boiling correlations specifically for cryogens. Computational models of several of these experiments have been constructed in the Generalized Fluid System Simulation Program (GFSSP), a network flow code developed at NASA’s Marshall Space Flight Center. The new Purdue-developed universal correlations cover the full boiling curve: onset of nucleate boiling, nucleate boiling, critical heat flux, and film boiling. These correlations have been coded into GFSSP user subroutines. The fluids modeled in this study are liquid hydrogen and liquid helium. Predictions of local wall temperature and pressure drop are presented and compared to the test data.

multi-phase flow↗
DOE OSTI2025

Thermal shock resistance of lightweight cements developed for geothermal conditions

Cements are a critical component in well construction, as they act to prevent well fluid and gas escape, prevent corrosion of the casing, and strengthen the wellbore to prevent deformation. Under the high temperature/pressure conditions common in geothermal systems, the injection of cold water for energy production is expected to induce cyclic damage to the borehole cement through the rapid temperature fluctuations. These “thermal shocks” are expected to cause casing shrinkage, annulus formation, and cement tensile stresses. To understand the effect of cold water injection on the wellbore environment, a set of rock-cement-steel samples were created to simulate the structure of a geothermal well. Lightweight thermally-insulating cement blends were tested under thermal shock conditions in this study. In each test, the samples were pressurized to an effective pressure of ~3.5 MPa and placed at high temperatures. Thermal shocks were performed by injecting cold water (~10-15 °C) through the samples at a constant rate while keeping the samples at high temperatures until the sample temperature stopped decreasing and deformation ceased. Eight thermal shock tests were conducted with each sample – two at 100 °C and six at 200 °C. Post-tests analysis was then conducted by cutting open each sample to examine the damage in each component of the simulated wellbore. Experimental results suggest that all samples experienced similar degrees of axial and lateral contraction during cold water injection, but for the most part this contraction is recoverable when injection halts. Post-test analysis revealed that fly ash cenosphere pre-treatment had the best effect on improving thermal shock resistance in the cement blends. Thermomechanical modeling of likely stress paths experienced by the cements during heating/cooling cycles shows that elasto-plastic cement constitutive behavior results in most plastic strain occurring during the initial heating steps, with mostly elastic strain occurring during the thermal shock cycles. In conclusion, this agrees with experimental evidence, suggesting that cement damage from shocking occurs via other mechanisms such as chemical alteration, corrosion, and fatigue.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Predicted Heat Flux Performance of Actively Cooled Tungsten-Armored Graphitic Foam Monoblocks

Tungsten (W)–armored graphitic foam monoblocks were developed for applications requiring high-Z plasma-facing material in long-pulse fusion experiments and ultimately deuterium-tritium fusion reactors. The monoblocks are an integrated material system combining the advantages of a chemical vapor deposited (CVD) W coating with a high-conductivity graphitic foam. The W is a high-melting-point, high-Z material with low tritium retention. The graphitic foam coupled to a swirl tube serves as a high-thermal-conductivity heat sink that cannot melt, although it can sublime at much higher temperatures than copper melts. Together, they comprise a robust plasma-facing component (PFC) weighing roughly 5% of an all-W component or 17% of a traditional W-coated copper heat sink.A single-channel mock-up consisting of four graphitic foam monoblocks equipped with a water-cooled swirl tube was fabricated for eventual testing in the 60-kW, EB-60, rastered electron beam at the Applied Research Laboratory of The Pennsylvania State University. Two monoblocks have a thin 50-μm-thick coating of pure W chemically vapor deposited over NbC and pure Nb interlayers. Two others have a 2-mm-thick pure W coating CVD on graphitic monoblocks using the same interlayers. The mock-up will be cooled with available 10 m/s, 0.7 MPa water with a 22°C inlet temperature and subjected to varying uniform heat loads up to 20 MW/m 2 . It is equipped with type-K thermocouples at various depths, and calibrated infrared thermography and spot pyrometry will be used to characterize the heated surface. Real-time water calorimetry will be used to ascertain the absorbed steady-state power and infer the heat flux during testing. Since testing cannot be done under prototypic divertor flow conditions, it is necessary to predict the thermal response of this novel PFC system and investigate the power sharing between radiation and convection at divertor heat flux levels and its inherent ability to avoid critical heat flux. Results are reported for predictions obtained from computational fluid dynamics models up to 30 MW/m 2 of steady-state uniform heat flux. Leading-edge heat loads of 30 MW/m 2 on a 2-mm-wide side strip were also investigated to ascertain if coating delamination is likely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Novel Multiphysics Multiscale Multiporosity Shale Gas Transport Model for Geomechanics/Flow Coupling in Steady and Transient States

Summary A novel multiphysics multiscale multiporosity shale gas transport (M3ST) model was developed to investigate shale gas transport in both transient and steady states. The microscale model component contains a kerogen domain and an inorganic matrix domain, and each domain has its own geomechanical and gas transport properties. Permeabilities of various shale cores were measured in the laboratory using a pulse decay permeameter (PDP) with different pore pressure and confining stress combinations. The PDP-measured apparent permeability as a function of pore pressure under two effective stresses was fitted using the microscale M3ST model component based on nonlinear least squares fitting (NLSF), and the fitted model parameters were able to provide accurate model predictions for another effective stress. The parameters and petrophysical properties determined in the steady state were then used in the transient-state, continuum-scale M3ST model component, which performed history matching of the evolutions of the upstream and downstream gas pressures. In addition, a double-exponential empirical model was developed as a powerful alternative to the M3ST model to fit laboratory-measured apparent permeability under various effective stresses and pore pressures. The developed M3ST model and the research findings in this study provided critical insights into the role of the multiphysics mechanisms, including geomechanics, fluid dynamics and transport, and the Klinkenberg effect on shale gas transport across different spatial scales in both steady and transient states.

Engineering↗
DOE OSTI2021

Computational Fluid Dynamics Modeling to Simulate a Combined Reforming Process for Syngas and Hydrogen Production

An Oxygen Transport Membrane (OTM) combined reforming technology for producing syngas and hydrogen integrates the advantages of multiple processes—steam methane reforming (SMR), autothermal reforming (ATR), an air separation unit (ASU)—into a single integrated technology. The OTM consists of a primary reforming tube, in which desulfurized natural gas is partially reformed by steam at high pressure in the presence of a metal catalyst. This process is followed in series by a ceramic OTM with a secondary reformer, in which residual methane reforms and O 2 - ions react with a portion of the CO and H 2 fuel to provide the heat to support both primary and secondary reforming. Although the OTM combined reformer technology for syngas and H 2 production has been substantially developed in the last decade, several challenges that affect the overall production efficiency and reliability are yet to be fully understood, addressed, and resolved. Therefore, developing Computational Fluid Dynamics (CFD) models that incorporate fluid dynamics, mass transport, kinetics, heat transport, and structural mechanics is critical to understanding and minimizing the probability of tube failures during the startup and operation. In this report, an exhaustive literature review was performed to survey the current state of technology for producing syngas and H 2 using either conventional or renewable energy sources. The feedstocks reviewed include natural gas and coal for the conventional technologies, whereas biomass, solar, wind, and nuclear energy for the renewable technologies. The existing industry-grade COMSOL multiphysics models of OTM were upgraded for the latest software release. In addition, they were improved to help achieve grid and solver independence and were successfully ported on the ORNL high-performance computing clusters to speed up their run times. A 42% reduction in the simulation run time was achieved. A new higher-fidelity CFD model of an OTM tube was developed in the StarCCM+ simulation platform. This new model was designed to simulate various physics using first principles, e.g., turbulent flow, heat transfer, and chemical reactions while avoiding unnecessary simplifications. The resulting predictions were qualitatively assessed and provided useful insights into the multiphysics complexity of an OTM tube.

08 HYDROGEN↗
DOE OSTI2022

Capillary Pressure to Relative Permeability: Task 7.3.4 Relative Permeability Analysis Final Report

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability relationship can be determined. Relative permeability is a numerical description of the interplay between two or more fluids and the porous media they flow through. It is a critical parameter for various tools used to characterized subsurface multiphase flow systems, such as numerical simulation for oil and gas development, carbon sequestration, and groundwater contamination remediation. Therefore, it is critical to get a good statistic distribution of relative permeability in the porous media under study. Empirical relationships for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility required to match laboratory derive relative permeability relationships. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability relationships. Existing laboratory-measured relative permeability data along with mercury intrusion capillary (MICP) data coupled with interfacial tension and contact angle measurements were used to determine the efficacy of this approach to relative permeability curve creation. The relative permeability relationships determined with this method were fit to existing laboratory data to elucidate common fitting parameters that were then used to create relative permeability relationships from MICP data that does not have an associated laboratory-measured relative permeability relationship.

58 GEOSCIENCES↗
NASA NTRS1984

Statistical mechanics of light elements at high pressure. VI - Liquid-state calculations with Thomas-Fermi-Dirac theory

A model free energy is developed for hydrogen-helium mixtures based on solid-state Thomas-Fermi-Dirac calculations at pressures relevant to the interiors of giant planets. Using a model potential similar to that for a two-component plasma, effective charges for the nuclei (which are in general smaller than the actual charges because of screening effects) are parameterized, being constrained by calculations at a number of densities, compositions, and lattice structures. These model potentials are then used to compute the equilibrium properties of H-He fluids using a charged hard-sphere model. The results find critical temperatures of about 0 K, 500 K, and 1500 K, for pressures of 10, 100, and 1000 Mbar, respectively. These phase separation temperatures are considerably lower (approximately 6,000-10,000 K) than those found from calculations using free electron perturbation theory, and suggest that H-He solutions should be stable against phase separation in the metallic zones of Jupiter and Saturn.

Macfarlane, J. J.↗
NASA NTRS1987

Three-step labyrinth seal for high-performance turbomachines

A three-step labyrinth seal with 12, 11, and 10 labyrinth teeth per step, respectively, was tested under static (nonrotating) conditions. The configuration represented the seal for a high-performance turbopump (e.g., the space shuttle main engine fuel pump). The test data included critical mass flux and pressure profiles over a wide range of fluid conditions at concentric, partially eccentric, and fully eccentric seal positions. The seal mass fluxes (leakage rates) were lower over the entire range of fluid conditions tested than those for data collected for similar straight and three-step cylindrical seals, and this conformed somewhat to expectations. However, the pressure profiles for the eccentric positions indicated little, if any, direct stiffness for this configuration in contrast to significant direct stiffness reported for the straight and three-step cylindrical seals over the range of test conditions. Seal dynamics depend on geometric configuration, inlet and exit parameters, fluid phase, and rotation. The method of corresponding states was applied to the mass flux data, which were found to have a pressure dependency for helium.

Hendricks, Robert C.↗
NASA NTRS1988

Convective effects on directional solidification of a simulated metal alloy

The first significant results of a ground-based experimental program which supports a low gravity space processing Spacelab experiment are reported. The phenomena which precipitate pluming and thus freckling in a metal alloy analog (ammonium chloride and water) are studied in detail and the sequential events leading to massive channeling and convection are optically documented. The pluming is shown to be other than a random burst of unstable fluid from a preferred channel but rather a natural occurrence resulting from a fundamental (Rayleigh-Benard) fluid dynamic instability at the density inversion interface. This extrapolates to critical size parameters appropriate to processing of actual metal alloys.

Mccay, T. D.↗
NASA NTRS1995

Durability of Silica-Protected Ceramics in Combustion Atmospheres

This overview of oxidation and corrosion effects on silica-protected ceramics and composites is divided into two parts: (1) External effects and (2) Internal effects. In regard to the first part, the literature indicates silica the best of the common protective oxides in pure oxygen. However in actual applications, there are four major degradation routes which must be considered: (1) Formation of SiO(g); (2) water vapor enhancement of oxidation rate of Si(OH)4 volatiles; (3) low level metal cation enhancement of oxidation rate; and (4) molten salt and slag corrosion. Several protective coating concepts have been discussed, which may minimize these effects. Internal effects deal with a microcrack providing a path for oxygen ingress to the oxygen sensitive fiber and fiber coating. This is a critical area of study for theses materials. Possible solutions involve a fluid glass to fill the cracks and/or oxygen getting along the walls of the crack.

Jacobson, Nathan S.↗
NASA NTRS2003

Magnetic Control of Solutal Buoyancy Driven Convection

Volumetric forces resulting from local density variations and gravitational acceleration cause buoyancy induced convective motion in melts and solutions. Solutal buoyancy is a result of concentration differences in an otherwise isothermal fluid. If the fluid also exhibits variations in magnetic susceptibility with concentration then convection control by external magnetic fields can be hypothesized. Magnetic control of thermal buoyancy induced convection in ferrofluids (dispersions of ferromagnetic particles in a carrier fluid) and paramagnetic fluids have been demonstrated. Here we show the nature of magnetic control of solutal buoyancy driven convection of a paramagnetic fluid, an aqueous solution of Manganese Chloride hydrate. We predict the critical magnetic field required for balancing gravitational solutal buoyancy driven convection and validate it through a simple experiment. We demonstrate that gravity driven flow can be completely reversed by a magnetic field but the exact cancellation of the flow is not possible. This is because the phenomenon is unstable. The technique can be applied to crystal growth processes in order to reduce convection and to heat exchanger devices for enhancing convection. The method can also be applied to impose a desired g-level in reduced gravity applications.

Ramachandran, N.↗
NASA NTRS1994

Oscillatory Thermocapillary Flow Experiment (OTFE)

An experiment was performed on oscillatory thermocapillary flow in the Glovebox aboard the USML-1 Spacelab which was launched in July, 1992. Cylindrical containers of 1 and 3 em in diameter were used. Silicone oils of 2 and 5 cSt viscosity were the test fluids. The fluid was heated by a cylindrical heater placed along the centerline of the container. The diameter of the heater was 10% of the container diameter. The fluid motion was studied by flow visualization. Although oscillations were observed briefly, bubbles generated in the fluid during the experiment disturbed the flow substantially so that the critical temperature differences could not be determined.

Kamotani, Y.↗
NASA NTRS2000

Extracellular calcium elicits a chemokinetic response from monocytes in vitro and in vivo

Recruitment of macrophages to sites of cell death is critical for induction of an immunologic response. Calcium concentrations in extracellular fluids vary markedly, and are particularly high at sites of injury or infection. We hypothesized that extracellular calcium participates in modulating the immune response, perhaps acting via the seven-transmembrane calcium-sensing receptor (CaR) on mature monocytes/macrophages. We observed a dose-dependent increase in monocyte chemotaxis in response to extracellular calcium or the selective allosteric CaR activator NPS R-467. In contrast, monocytes derived from mice deficient in CaR lacked the normal chemotactic response to a calcium gradient. Notably, CaR activation of monocytes bearing the receptor synergistically augmented the transmigration response of monocytes to the chemokine MCP-1 in association with increased cell-surface expression of its cognate receptor, CCR2. Conversely, stimulation of monocytes with MCP-1 or SDF-1alpha reciprocally increased CaR expression, suggesting a dual-enhancing interaction of Ca(2+) with chemokines in recruiting inflammatory cells. Subcutaneous administration in mice of Ca(2+), MCP-1, or (more potently) the combination of Ca(2+) and MCP-1, elicited an inflammatory infiltrate consisting of monocytes/macrophages. Thus extracellular calcium functions as an ionic chemokinetic agent capable of modulating the innate immune response in vivo and in vitro by direct and indirect actions on monocytic cells. Calcium deposition may be both consequence and cause of chronic inflammatory changes at sites of injury, infection, and atherosclerosis.

Non-NASA Center↗