Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Properties”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A review of high order strong stability preserving two-derivative explicit, implicit, and IMEX methods

High order strong stability preserving time discretizations ensure the nonlinear non-inner-product strong stability properties of spatial discretizations suited for the stable simulation of hyperbolic PDEs in a wide variety of application areas including fluid dynamics, magnetohydrodynamics, semiconductor devices, electromagnetics, and astrophysics. Over the past decade multiderivative time-stepping have been increasingly used for the time-evolution hyperbolic PDEs, so that the strong stability properties of these methods have become important. In this work we review sufficient conditions for a two-derivative multistage method to preserve the strong stability properties of spatial discretizations in a forward Euler and different conditions on the second derivative. In particular we present the strong stability preserving theory for explicit and implicit two-derivative Runge–Kutta schemes, including a special condition on the second derivative under which these implicit methods may be unconditionally strong stability preserving. This special condition is natural for the stiff component of wide range of plasma physics problems, and can be useful in the context of strong stability preserving implicit-explicit multi-derivative Runge–Kutta schemes, where the time-step restriction is then independent of the stiff term. Lastly, we present the strong stability preserving theory for implicit-explicit multi-derivative general linear methods, and some novel second and third order methods where the time-step restriction is independent of the stiff term.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Molten Chloride Thermophysical Properties, Chemical Optimization, and Purification

Next-generation concentrating solar power (CSP) technology requires high-temperature heat-transfer fluids and thermal energy storage that can operate in the temperature range of 500°–720°C. This high temperature demands a new chemistry for the heat-transfer fluids that must have high-temperature thermal and chemical stabilities, high energy density, and low corrosion rate on metallic materials used in the heat exchangers, piping, and thermal energy storage tanks. A ternary MgCl 2 -KCl-NaCl chloride-salt system has been proposed to achieve such goals, given its high-temperature thermal stability. The chloride-salt-based systems have been used by the metallurgy industry as the electrolyte to electrolytically extract metallic magnesium. They have also attracted attentions from the nuclear industry for developing the molten salt fast reactors that can operate at close to atmospheric pressure and at higher temperature than traditional light-water reactors. However, the corrosion properties of the ternary salt are not well known. The corrosion mechanism is believed to be correlated to the presence of MgOHCl impurity at the operation temperature. MgOHCl is fundamentally caused by the hygroscopic nature of the MgCl 2 component. Therefore, the first key objective of this project is to understand the nature of the corrosiveness of the ternary chloride-salt system and its prevention by designing new lab-scale purification processes. Because of its importance, a working group called the Chloride Salt Collective across multiple U.S. Department of Energy national laboratories is established to gain experimental agreement on the corrosion behaviors. In addition, water generation during salt dehydration could lead to highly corrosive HCl/water mixture via a hydrolysis process between water and anhydrous MgCl 2 or MgCl 2 hydrates. MgOHCl will eventually be converted to MgO particulates regardless of the chemical path chosen (i.e., by thermal decomposition, chemical purification with an active metal, or electrochemical purification utilizing the reduction power of electron flows). Large MgO particles with high hardness can wear out bearing and bushing and even clog the delicate channels in the primary heat exchanger. Because of this, a second objective is to understand the chemistry behind salt purification and transfer the scientific understanding to a set of engineering know-how for large-scale salt melter design. Because there is limited information on the thermophysical properties and the sensitivity of these properties around the targeted salt composition, a third objective is to measure thermophysical properties relevant to system and component design and support Topic 1 Liquid Pathway development. For example, composition shift due to salt reaction with water can be a potential issue because it can lead to changes of thermophysical properties, such as (1) melting point (or freeze point), which can jeopardize safe operation of the supercritical CO 2 power cycle proposed for next-generation CSP, and (2) heat capacity and density, which can decrease the overall power-generation capability of the plants. The value of the thermal conductivity of the selected salt composition is also a critical criterion for the solar receiver design.

14 SOLAR ENERGY↗
DOE OSTI2020

Seismicity and Stress Associated With a Fluid-Driven Fracture: Estimating the Evolving Geometry

A coupled approach, combining the theory of rate- and state-dependent friction and methods from poroelasticity, forms the basis for a quantitative relationship between displacements and fluid leak-off from a growing fracture and changes in the rate of seismic events in the region surrounding the fracture. Poroelastic Green's functions link fracture aperture changes and fluid flow from the fracture to changes in the stress field and pore pressure in the adjacent formation. The theory of rate- and state-dependent friction provides a connection between Coulomb stress changes and variations in the rate of seismic events. Numerical modeling indicates that the Coulomb stress changes can vary significantly between formations with differing properties. The relationship between the seismicity rate changes and the changes in the formation stresses and fluid pressure is nonlinear, but a transformation produces a quantity that is linearly related to the aperture changes and fluid leak-off from the fracture. The methodology provides a means for mapping changes in seismicity into fracture aperture changes and to image an evolving fracture. An application to observed microseismicity associated with a hydrofracture reveals asymmetric fracture propagation within two main zones, with extended propagation in the upper zone. Furthermore, the time-varying volume of the fracture agrees with the injected volume, given by the integration of rate changes at the injection well, providing validation of the estimated aperture changes.

58 GEOSCIENCES↗
DOE OSTI2020

Confinement-mediated phase behavior of hydrocarbon fluids: Insights from Monte Carlo simulations

The phase behavior of hydrocarbon fluids confined in porous media has been discussed to deviate significantly from that in the bulk environment due to the existence of sub-10nm pores. Though experiments and simulations have measured the bubble/dew points and sorption isotherms of hydrocarbons confined in both natural and synthetic nanopores, the confinement effects in terms of the strength of fluid-pore interactions tuned by surface wettability and chemistry have received comparably less discussion. More significantly, the underlying physics of confinement-induced phenomena remain obfuscated. In this work, we studied the phase behavior and capillary condensation of n-hexane to understand the effects of confinement at the molecular level. To systematically investigate the pore effects, we constructed two types of wall confinements; one is a structureless virtual wall described by the Steele potential and the other one is an all-atom amorphous silica structure with surface modified by hydroxyl groups. Our numerical results demonstrated the importance of fluid-pore interaction, pore size, and pore morphology effects in mediating the pressure-volume-temperature (PVT) properties of hydrocarbons. The most remarkable finding of this work was that the saturation pressure predicted from the van der Waals-type adsorption isothermal loop could be elevated or suppressed relative to the bulk phase, as illustrated in the graphical abstract. As the surface energy (i.e., fluid-pore interaction) decreased, the isothermal vapor pressure increased, indicating a greater preference for the fluid to exist in the vapor state. Sufficient reduction of the fluid-pore interactions could even elevate the vapor pressure above that of the bulk fluid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Physics-based Induced Earthquake Forecasting: Process Understanding, and Hazards Mitigation

Disposal of saltwater co-produced with oil and gas is linked to elevated seismicity in the Central and Midwest US. There is a concern that these events may lead to widespread damage and an overall increase in seismicity. Thus an improved understanding of the spatially and temporally variable deformation and stress field associated with fluid injection operation is critically necessary for evaluating time-varying seismic hazards. Despite the improvements in seismic monitoring capacity and the resulting decrease in the magnitude detection threshold, estimates of induced earthquake probability remain elusive due to insufficient models incapable of accounting for the complex physics governing the process of induced seismicity. The proposed research effort will comprehensively analyze, integrate, and interpret geodetic, injection, and seismic data in the vicinity of the injection sites in Oklahoma to resolve the 4-dimensional distribution of pore pressure and stress in the shallow crust. This project, in particular, is focused on exploring the statistical relation between injection operation and increased earthquake hazard. The amplitude of and the extent to which pore pressure changes are determined by some factors, in particular, the hydrogeological properties of the rocks, such as diffusivity. Thus the available deformation data will be used to constrain hydrogeological properties of the medium, to accurately resolve the evolution of crustal stresses due to fluid injection. Having the time-varying models of stress changes, a statistical framework will be implemented to estimate the time-dependent probability of large earthquakes on the nearby fault systems. These data and models help to improve seismic hazard estimates and aid in constructing operational-induced earthquake forecast models. This information can also be integrated into the updated U.S. National Seismic Hazard Map, which local communities and authorities use in their earthquake risk estimates and mitigation efforts.

58 GEOSCIENCES↗
DOE OSTI2021

Internally cooled membrane-based absorber for dehumidification and water heating: Validated model and simulation study

This report presents a comprehensive model for advanced membrane-based absorber components accounting for three separate fluid streams, one of them providing internal cooling, with complex flow patterns. The model was implemented with various liquid desiccants and their properties (including LiCl, CaCl 2 , and [emim][OMS] ionic liquid [IL]). The model has been validated using data from two laboratory prototype absorbers, one with CaCl 2 as the working fluid and one with an IL as the working fluid. Entropy analysis was further carried out to verify the model and understand effect of internal cooling on the entropy variation of three fluid streams. The developed model and codes are expected to enable detailed configuration optimization and provide in-depth understanding of three-fluid heat and mass exchanger (HMX) performance. The paper also describes parametric studies of the HMX components in utilizing the latent heat removed in space cooling to heat water and numerically explores the unique benefit of its application in a semi-open absorption heat pump water heater. The results show that a single three-fluid HMX has the potential to achieve simultaneous dehumidification and water heating efficiently at cost effectiveness, particularly in hot and humid climates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Relating Hydro–Mechanical and Elastodynamic Properties of Dynamically Stressed Tensile–Fractured Rock in Relation to Applied Normal Stress, Fracture Aperture, and Contact Area

We exploit nonlinear elastodynamic properties of fractured rock to probe the micro-scale mechanics of fractures and understand the relation between fluid transport and fracture aperture under dynamic stressing. Experiments were conducted on rough, tensile-fractured Westerly granite subject to triaxial stresses. We measure fracture permeability for steady-state fluid flow with deionized water. Pore pressure oscillations are applied at amplitudes ranging from 0.2 to 1 MPa at 1 Hz frequency. During dynamic stressing we transmit ultrasonic signals through the fracture using an array of piezoelectric transducers (PZTs) to monitor evolution of interface properties. We examine the influence of fracture aperture and contact area by conducting measurements at effective normal stresses of 10–20 MPa. Additionally, the evolution of contact area with stress is characterized using pressure sensitive film. These experiments are conducted separately with the same fracture and map contact area at stresses from 9 to 21 MPa. The measurements are a proxy for “true” contact area for the fracture surface and we relate them to elastic properties using the calculated PZT sensor footprints via numerical modeling of Fresnel zones. We compare the elastodynamic response of the fracture using the stress-induced changes in ultrasonic wave velocities for transmitter-receiver pairs to image spatial variations in contact properties. We show that nonlinear elasticity and permeability enhancement decrease with increasing normal stress. Additionally, post-oscillation wave velocity and permeability exhibit quick recoveries toward pre-oscillation values. Estimates of fracture contact area (global and local) demonstrate that the elastodynamic and permeability responses are dominated by fracture topology.

58 GEOSCIENCES↗
DOE OSTI2023

Permeability of Single–Layer–Free–Standing Meshes at Varying Capillary Pressure via a Novel Method

The permeability of mesh wicks is important for various applications, including two–phase heat transfer. However, the understanding of the permeability of single–layer, free–standing mesh wicks, with liquid–gas interfaces on both sides, is limited. A novel and simpler method is presented to determine the permeability of a free–standing wick and apply it to a representative mesh. This method involves modifying the capillary pressure via elevation and simultaneously measuring the permeability to determine the permeability–capillary pressure relationship. When applied to a copper mesh with plain weave having undergone surface cleaning, the permeability is found to decrease as capillary pressure for deionized water increases. A dimensional analysis is presented to generalize this data for other mesh sizes with similar weaves and fluids. The behavior of mesh in application is modeled, based on the integration of Darcy's law with an analytic function fit to measured data, and parametric studies are conducted to investigate the superficial velocity of liquids through the mesh under varying driving pressures, transport lengths, and liquid viscosity, based on the obtained capillary pressure–permeability relationship. This study provides valuable insights into the transport properties of mesh wicks, with potential applications in fields such as electronics cooling, electrochemical devices, and fluid purification technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗
DOE OSTI2022

Explaining the lack of power degradation of energy confinement in wide pedestal quiescent H-modes via transport modeling

Wide pedestal quiescent H (WPQH)-mode is an attractive scenario for future burning plasmas as they operate without ELMs. WPQH is characterized by formation of a wider and higher pedestal (than quiescent H-mode), and broadband fluctuations in the pedestal. Unlike conventional H-modes, where the energy confinement time reduces with increasing heating power, the WPQH plasmas reported in this paper do not show power degradation of the energy confinement. As the injected neutral beam power was increased, reduced core (ρ ≤ 0.45) transport calculated by transp, as well as increased core temperatures, pressure gradient and diamagnetic E × B shear rate were observed. The reduction in the heat transport and rapid decrease in the ion temperature gradient scale length suggest the formation of an ion internal transport barrier (ITB) that was accompanied by increased stored energy in the core. Quasilinear turbulent transport modeling using the trapped gyro Landau fluid (tglf) code was used to predict the ITB and its turbulence stability properties. By using profiles and equilibria produced by matching the transp transport fluxes with the tglf transport model within the tgyro transport solver, the energy confinement time captures the experimentally observed insensitivity to the increased PNBI. Linear stability analysis reveals that drift-wave instabilities in the core are stabilized by E × B shear, T i /T e ratio and Shafranov shift; the latter was found to have the strongest effect on the turbulence suppression at the highest heating level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Machine learning of consistent thermodynamic models using automatic differentiation

In this study, we propose a data-driven method to describe consistent equations of state (EOS) for arbitrary systems. Complex EOS are traditionally obtained by fitting suitable analytical expressions to thermophysical data. A key aspect of EOS is that the relationships between state variables are given by derivatives of the system free energy. In this work, we model the free energy with an artificial neural network and utilize automatic differentiation to directly learn the derivatives of the free energy. We demonstrate this approach on two different systems, the analytic van der Waals EOS and published data for the Lennard-Jones fluid, and we show that it is advantageous over direct learning of thermodynamic properties (i.e., not as derivatives of the free energy but as independent properties), in terms of both accuracy and the exact preservation of the Maxwell relations. Furthermore, the method implicitly provides the free energy of a system without explicit integration.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

3D fault architecture controls the dynamism of earthquake swarms

The vibrant evolutionary patterns made by earthquake swarms are incompatible with standard, effectively two-dimensional (2D) models for general fault architecture. Here we leverage advances in earthquake monitoring with a deep-learning algorithm to image a fault zone hosting a 4-year-long swarm in southern California. We infer that fluids are naturally injected into the fault zone from below and diffuse through strike-parallel channels while triggering earthquakes. A permeability barrier initially limits up-dip swarm migration but ultimately is circumvented. This enables fluid migration within a shallower section of the fault with fundamentally different mechanical properties. Our observations provide high-resolution constraints on the processes by which swarms initiate, grow, and arrest. These findings illustrate how swarm evolution is strongly controlled by 3D variations in fault architecture.

58 GEOSCIENCES↗
DOE OSTI2021

Thermal-Hydrological-Mechanical Modelling of Stockton University Reservoir Cooling System, Fine Scale Stress Test Modelling

Mesh, properties, initial conditions, injection/withdrawal rates for modelling thermal, hydrological, and mechanical effects of fluid injection to and withdrawal from ground for Stockton University reservoir cooling system (aquifer storage cooling system), Galloway, New Jersey, for unscheduled two hour injection at 133 % designed capacity, on fine scale grid, with some results. Second simulation of J.T. Smith, E. Sonnenthal, P. Dobson, P. Nico, and M. Worthington, 2021. Thermal-hydrological-mechanical modeling of Stockton University reservoir cooling system, Proceedings of the 46th Workshop on Geothermal Reservoir Engineering, Stanford University, SGP-TR-218, from which Figures 6-9, pertain.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Thermal-Hydrological-Mechanical Modelling of Stockton University Reservoir Cooling System, Large Scale Grid

Mesh, properties, initial conditions, injection/withdrawal rates for modeling thermal, hydrological, and mechanical effects of fluid injection to and withdrawal from ground for Stockton University reservoir cooling system (aquifer storage cooling system), Galloway, New Jersey, on large scale grid, with some results. First simulation of J.T. Smith, E. Sonnenthal, P. Dobson, P. Nico, and M. Worthington, 2021. Thermal-hydrological-mechanical modeling of Stockton University reservoir cooling system, Proceedings of the 46th Workshop on Geothermal Reservoir Engineering, Stanford University, SGP-TR-218, from which Figures 1-5 pertain.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Gravitational vortex disruption in mesofluidic separators

Typically, bump array devices separate submicrometer particles where gravitational effects remain small. However, as particle size increases from a few micrometers to hundreds of micrometers, gravitational effects become important and settling may induce vortex formation. Here, this article shows under laminar flow that these vortices suffice to disrupt flow fields, frustrating particle separation. The limiting diameter, above which vortices do form, is determined, and the dependence of the limiting diameter on viscosity, particle properties, and separator orientation is quantified.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

An advanced membrane-based HMX model and its application in utilizing sensible and latent heat removed in space cooling

This paper presents a comprehensive simulation model for advanced membrane-based heat and mass exchanger (HMX) components which account for three separate fluid streams with complex flow patterns. The model was implemented with various liquid desiccants and their properties (including LiCl, CaCl and ionic liquid (IL)). Preliminary model validation is presented using data from a laboratory prototype absorber. The paper also presents implementation of the model into a sorption system modeling platform and analyzes potential benefit analysis of the HMX component in utilizing the sensible and latent heat removed in space cooling to heat water. The developed model and codes are expected to enable detailed optimization studies and identify novel configurations for the three-fluid HMX maximum performance and cost effectiveness.

Gao, Zhiming↗
DOE OSTI2020

Effects of salinity and pH on the spectral induced polarization signals of graphite particles

The electrical property of micrometre-sized graphite particles was investigated under different particle concentration, particle size, fluid conductivity and pH conditions. Due to its large internal electronic conductivity and ability to polarize under external potential field, significant enhancement of its spectral induced polarization (SIP) responses is observed when graphite is included in sand mixtures. While a small amount of graphite inclusion significantly increases the SIP response of its mixtures with sand, further concentration increase does not necessarily lead to a proportional increase of the SIP response. This is shown to be related to the formation of graphite aggregates at higher concentrations. Changes of fluid salinity have a significant effect on graphite's SIP behaviour. This includes a positive impact on normalized chargeability and imaginary conductivity, but a negative impact on chargeability and relaxation time constant. The effect of pH on the SIP response of graphite is small but shows consistent trend, where pH increase leads to a decrease of both the chargeability and relaxation time constant. The underlying cause of this effect is not clear.

58 GEOSCIENCES↗