Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive transport”

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

Alquimia v1.0: a generic interface to biogeochemical codes – a tool for interoperable development, prototyping and benchmarking for multiphysics simulators

Alquimia v1.0 is a generic interface to geochemical solvers that facilitates development of multiphysics simulators by enabling code coupling, prototyping and benchmarking. The interface enforces the function arguments and their types for setting up, solving, serving up output data and carrying out other common auxiliary tasks while providing a set of structures for data transfer between the multiphysics code driving the simulation and the geochemical solver. Alquimia relies on a single-cell approach that permits operator splitting coupling and parallel computation. We describe the implementation in Alquimia of two widely used open-source codes that perform geochemical calculations: PFLOTRAN and CrunchFlow. We then exemplify its use for the implementation and simulation of reactive transport in porous media by two open-source flow and transport simulators: Amanzi and ParFlow. We also demonstrate its use for the simulation of coupled processes in novel multiphysics applications including the effect of multiphase flow on reaction rates at the pore scale with OpenFOAM, the role of complex biogeochemical processes in land surface models such as the E3SM Land Model (ELM) and the impact of surface–subsurface hydrological interactions on hydrogeochemical export from watersheds with the Advanced Terrestrial Simulator (ATS). These applications make it apparent that the availability of a well-defined yet flexible interface has the potential to improve the software development workflow, freeing up resources to focus on advances in process models and mechanistic understanding of coupled problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2012

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗
NASA NTRS2016

Biological Potential in Serpentinizing Systems

Generation of the microbial substrate hydrogen during serpentinization, the aqueous alteration of ultramafic rocks, has focused interest on the potential of serpentinizing systems to support biological communities or even the origin of life. However the process also generates considerable alkalinity, a challenge to life, and both pH and hydrogen concentrations vary widely across natural systems as a result of different host rock and fluid composition and differing physical and hydrogeologic conditions. Biological potential is expected to vary in concert. We examined the impact of such variability on the bioenergetics of an example metabolism, methanogenesis, using a cell-scale reactive transport model to compare rates of metabolic energy generation as a function of physicochemical environment. Potential rates vary over more than 5 orders of magnitude, including bioenergetically non-viable conditions, across the range of naturally occurring conditions. In parallel, we assayed rates of hydrogen metabolism in wells associated with the actively serpentinizing Coast Range Ophiolite, which includes conditions more alkaline and considerably less reducing than is typical of serpentinizing systems. Hydrogen metabolism is observed at pH approaching 12 but, consistent with the model predictions, biological methanogenesis is not observed.

the potential of serpentinizing systems to support↗
NASA NTRS2017

Investigating Weathering of Basaltic Materials in Gale Crater, Mars: A Combined Laboratory, Modeling and Terrestrial Field Approach

Recent observations from Gale Crater, Mars document past aqueous alteration both in the formation of the Stimson sandstone unit, as well as in the formation of altered fractures within that unit. Geochemical and mineralogical data from Curiosity also suggest Fe-rich amorphous weathering products are present in most samples measured to date. Here we interpret conditions of possible past weathering in Gale Crater using a combination of field, laboratory, and modeling work. In order to better understand secondary Fe-rich phases on Mars, we are examining formation of weathering products in high Fe and Mg and low Al serpentine soils in the Klamath Mountains, CA. We have isolated potential weathering products from these soils, and are analyzing them using synchrotron μXRF and μXRD as well as FullPat for a direct comparison to analyses from Gale Crater. In order to interpret the implications of the persistence of potential secondary Fe-containing phases on Mars, we are also measuring the dissolution rates of the secondary weathering products allophane, Fe-rich allophane, and hisingerite. Ongoing dissolution experiments of these materials suggest that they dissolve significantly more rapidly than more crystalline secondary minerals with similar chemical compositions. Finally, to quantify the specific conditions of past aqueous alteration in Gale Crater we are performing reactive transport modeling of a range of possible past environmental conditions. Specifically, we are testing the conditions under which a Stimson unit-like material forms from a parent material similar to Rocknest or Bagnold eolian deposits, and the conditions under which observed altered fracture zones form from a Stimson unit-like parent material. Our modeling results indicate that the formation of the Stimson unit is consistent with leaching of an eolian deposit with a solution of pH = 6-8, and that formation of the altered fracture zones is consistent with leaching with a very acidic (pH = 2-3) high sulfate solution containing Ca. These results suggest circumneutral pH conditions during authigenesis or early diagenesis in the Stimson formation sediments followed by diagenetic alteration by very acidic solutions along fracture zones.

Hausrath, Elisabeth↗
NASA NTRS2017

Energy in a Planetary Context

The potential present day habitability of solar system bodies beyond Earth is limited to subsurface environments, where the availability of energy in biologically useful form is a paramount consideration. Energy availability is commonly quantified in terms of molar Gibbs energy changes for metabolisms of interest, but this can provide an incomplete and even misleading picture. A second aspect of life's requirement for energy - the rate of delivery, or power - strongly influences habitability, biomass abundance, growth rates, and, ultimately, rates of evolution. We are developing an approach to quantify metabolic power, using a cell-scale reactive transport model in which physical and chemical environmental parameters are varied. Simultaneously, we evaluate cell-specific energy flux requirements and their dependence on environmental "extremes". Comparison of metabolic power supply and demand provides a constraint on how biomass abundance varies across a range of environmental parameters, and thereby a prediction of the relative habitability of different environments. We are evaluating the predictive capability of this approach through comparison to observed distributions of microbial abundance in a range of subsurface (predominantly serpentinizing) systems.

MARINE GEOMICROBIOLOGY↗
DOE OSTI2022

Geologic Disposal Safety Assessment (GDSA) Biosphere Model Development

The Spent Fuel and Waste Science and Technology Campaign of the U.S. Department of Energy Office of Nuclear Energy, Office of Spent Fuel and Waste Disposition is conducting research and development on geologic disposal of spent nuclear fuel and high-level nuclear waste. This work includes the Geologic Disposal Safety Assessment (GDSA) program which is charged with development of generic deep geologic repository concepts and system performance assessment models. One part of the GDSA framework is the development of a biosphere model capable of assessing doses to potential receptors exposed to radionuclides released from geologic disposal sites. As part of the GDSA framework, a biosphere model compatible with the PFLOTRAN massively parallel subsurface flow and reactive transport code is under development. The PFLOTRAN model provides the radionuclide source term for the biosphere model. The GDSA Biosphere model then assesses the potential movement of radionuclides through the surface biosphere and the subsequent exposure to a human receptor living in the biosphere. The biosphere model includes pathways originating from the groundwater as well as pathways originating from surface water bodies that have a water exchange with a contaminated groundwater body. The pathways for human exposure include consumption of drinking water, irrigated crops, meat animals, aquatic vegetation, and animals, etc.; external exposure from irrigated ground surfaces, surface water bodies, recreational activities, etc.; and inadvertent exposures such as ingestion of contaminated soils or shower water, etc. The GDSA Biosphere model was designed to be flexible and generic in order to accommodate a variety of different sites and climate states. This presentation will present the on the purpose, design, and development progress of the GDSA Biosphere Model.

GDSA, biosphere, repository↗
DOE OSTI2020

Pump-and-Treat Impacts on Aquifer Conditions and the Development of System-Wide Performance Strategies - 20377

The pump-and-treat (P and T) system on the Hanford Site's Central Plateau currently treats multiple contaminants within the 200-ZP-1 operable unit, including nitrate through its biological treatment processes. Although it was initially designed to support the 200-ZP-1 Remedial Action Objectives, the P and T system has evolved in use to address contaminants from other regions of the aquifer as well. The P and T influents are biogeochemically distinct, treating groundwaters with elevated concentrations of cyanide, uranium, technetium-99, and nitrate. Given the varying geochemical signatures from current and potential future sources, operational conditions, and amendments considered for treatment processes, there is a need to identify the impacts that the P and T effluent has on the aquifer, and how the aquifer affects P and T performance. To this end, laboratory efforts have focused on characterizing P and T effluent and its biogeochemical effects on the sediments through a series repacked column experiments. The P and T effluent is rich in particulate matter. A systematic inspection of this particulate matter in the P and T effluent showed significant levels of manganese/iron oxides, calcite crystals, relatively large organic matter molecules, and biomass. Therefore, the column experiments were designed and conducted for four different treatment systems: artificial groundwater, filtered P and T effluent, non-filtered P and T effluent, and chlorinated P and T effluent. Different levels of microbiologically- and geochemically-induced changes in the columns occurred and modified the flow systems. For example, physical filtration occurred in the sediments where the particulate material had accumulated in the top portion of the column. Reactive transport modeling of the laboratory experiments is being carried out so that the reaction network can be extrapolated to the field scale. These simulations will provide predictive evaluations of system performance based on P and T operational changes and aquifer mass loading and capacity estimates. This information can be used to address combined subsurface and operational challenges, such as the need for frequent well redevelopment to maintain injection capacity. This information also provides the technical basis for decisions related to the P and T system operations that support remedy optimization efforts and short- and long-term remedy decisions related to multiple operable units on the Hanford Site Central Plateau. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Nano- to microscale three-dimensional morphology relevant to transport properties in reactive porous composite paint films

The quantitative evaluation of the three-dimensional (3D) morphology of porous composite materials is important for understanding mass transport phenomena, which further impact their functionalities and durability. Reactive porous paint materials are composites in nature and widely used in arts and technological applications. In artistic oil paintings, ambient moisture and water and organic solvents used in conservation treatments are known to trigger multiple physical and chemical degradation processes; however, there is no complete physical model that can quantitatively describe their transport in the paint films. In the present study, model oil paints with lead white (2PbCO 3 ·Pb(OH) 2 ) and zinc white (ZnO) pigments, which are frequently found in artistic oil paintings and are associated with the widespread heavy metal soap deterioration, were studied using synchrotron X-ray nano-tomography and unilateral nuclear magnetic resonance. This study aims to establish a relationship among the paints’ compositions, the 3D morphological properties and degradation. This connection is crucial for establishing reliable models that can predict transport properties of solvents used in conservation treatments and of species involved in deterioration reactions, such as soap formation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

TReactMech v4.217

TReactMech couples geomechanical processes (poroelasticity, failure, and inelastic strain) with multiphase nonisothermal flow (derived from TOUGH2) and reactive geochemical transport. At its core is the reactive-transport code TOUGHREACT v4.13. TReactMech is an efficient hybrid parallel simulator, solving the geomechanics using finite elements and MPI/PetSc, the multiphase flow using integrated finite difference and MPI/PETSc, and the reactive chemistry using OpenMP. The advantages of TReactMech are in its multiphase flow capabilities (e.g., supercritical CO2, supercritical water, air) and parallel geomechanics including full 3-D stress tensor, shear and tensile failure, coupled to porosity and permeability changes. It is backwardly compatible with TOUGH2 and TOUGHREACT v4.13, allowing for easier transitions between the codes. TReactMech can be used to simulate many natural and engineered subsurface systems, including geothermal reservoirs, borehole heat exchangers, geologic carbon sequestration, geologic storage of nuclear waste, groundwater resources, weathering, sediment diagenesis, seafloor hydrothermal circulation, hydrofracturing in unconventional reservoirs, and injection/production-induced surface deformation.

Sonnenthal, Eric↗
DOE OSTI2026

Concerted Electron-Ion Transport by Polyacrylonitrile Elucidated with Reactive Deep Learning Potentials

Charge transport in polymers, such as polyacrylonitrile (PAN), is crucial for electronics and energy storage. For instance, PAN can transport cations e.g., Li + , by facilitating dynamic cation-nitrile coordination in batteries. However, little is known regarding the underlying role of complex reactive polymer configurations. Herein, we develop a deep-learning potential, trained on ab initio energies and forces of nonequilibrium reactive PAN configurations, to unravel the kinetics of PAN cyclization initiated by a nucleophile (OH – dissociated from LiOH) attacking the terminal nitrile carbon. We find, based on the reaction free-energetics, rates, and charge analysis, that the nucleophile attack producing the first ring is the rate-limiting step, which subsequently triggers Li + -coupled electron transfer along the PAN backbone, causing ∼10 4 times faster sequential ring-formation of the remaining nitriles. PAN’s extended configurations, where dipolar and H-bonding interactions are minimal, enable such rapid kinetics. By validating our computational findings with IR and NMR experiments, we establish a pathway for designing reactive polymers with enhanced charge transport for energy applications.

Chahal-Crockett, Rajni [Oak Ridge National Laborat↗
NASA NTRS1972

Transport phenomena of reactive fluid flow in heterogeneous combustion processes.

A previously developed computer program was used to model two transient hybrid combustion processes involving tubes of solid Plexiglas. In the first study, representing combustion of a hybrid rocket, the oxidizing gas was oxygen, and calculations were continued sufficiently long to obtain steady-state values. Systematic variations were made in reaction rate constant, mass flow rate, and pressure, alternatively using constant and temperature dependent regression rate models for the fuel surface. Consistent results were obtained, as is evidenced by the values for the mass function of the reaction product and the flame temperature, for which plots are supplied. In the second study, fire initiation in a duct was studied, with an air mixture as the oxidizing gas. It was demonstrated that a satisfactory flame spread mechanism could be reproduced on the computer. In both of the above applications, the general, transient, two-dimensional conservation equations were represented, together with chemical reactions, solid-fuel interface conditions, and heat conduction in the solid fuel.

Hung, W. S. Y.↗
NASA NTRS1992

Reactive collision terms for fluid transport theory

A modified Bhatnagar-Gross-Krook (BGK) collision method is used for the derivation of reactive collision terms to allow for complete higher-moment approximations in fluid transport theory. The general reactive collision terms for the 8-, 10- and 13-moment approximation for binary reactions are derived. Then, reactive collision terms for specific chemistry are derived with the assumptions: (1) all reactants and products have 13-moment distribution functions; (2) simple, classical assumptions provide sufficient models for reaction dynamics; and (3) the reaction rate is energy-independent. The derived specific formulas reflect two extremes in reactive collision dynamics: direct and indirect reaction mechanisms. Finally, considerations for correct use of the reactive collision terms are discussed in the context of space plasma environments.

Eccles, J. V.↗
DOE OSTI2026

Computer simulation of reaction and transport of core-scale serpentinization of Fe-bearing olivine and evolution of geological hydrogen

We present reactive flow and transport simulations of core-scale serpentinization of forsteritic, Fe-bearing olivine under hydrothermal conditions. The model captures fluid flow, evolving porosity–permeability, and redox-controlled H 2 production in a cylindrical core over 55 days at 245 °C and 37 bar. Key processes include olivine dissolution, precipitation of serpentine and magnetite, and oxidation of Fe2+ to Fe 3+ . Here, results show strong coupling among flow velocity, alteration front propagation, and spatially heterogeneous H 2 generation. Elevated H2 concentrations align with inlet-localized magnetite precipitation, consistent with experiments. Secondary mineral formation reduces porosity and permeability, altering transport pathways. Reactive flow at the fluid–solid interface governs H 2 generation rates and distribution. Dissolved SiO2 promotes serpentine and talc formation by suppressing brucite, while bicarbonate extends reaction duration and moderates surface complexation effects.

Geologic hydrogen↗
DOE OSTI2021

Simulation of hardened cement degradation and estimation of uncertainty in predicted failure times with peridynamics

Modeling the degradation of cement-based infrastructure due to aqueous environmental conditions continues to be a challenge. In order to develop a capability to predict concrete infrastructure failure due to chemical degradation, here we created a chemomechanical model of the effects of long-term water exposure on cement paste. The model couples the mechanical static equilibrium balance with reactive–diffusive transport and incorporates fracture and failure via peridynamics (a meshless simulation method). The model includes fundamental aspects of degradation of ordinary Portland cement (OPC) paste, including the observed softening, reduced toughness, and shrinkage of the cement paste, and increased reactivity and transport with water induced degradation. This version of the model focuses on the first stage of cement paste decalcification, the dissolution of portlandite. Given unknowns in the cement paste degradation process and the cost of uncertainty quantification (UQ), we adopt a minimally complex model in two dimensions (2D) in order to perform sensitivity analysis and UQ. We calibrate the model to existing experimental data using simulations of common tests such as flexure, compression and diffusion. Then we calculate the global sensitivity and uncertainty of predicted failure times based on variation of eleven unique and fundamental material properties. We observed particularly strong sensitivities to the diffusion coefficient, the reaction rate, and the shrinkage with degradation. Also, the predicted time of first fracture is highly correlated with the time to total failure in compression, which implies fracture can indicate impending degradation induced failure; however, the distributions of the two events overlap so the lead time may be minimal. Extension of the model to include the multiple reactions that describe complete degradation, viscous relaxation, post-peak load mechanisms, and to three dimensions to explore the interactions of complex fracture patterns evoked by more realistic geometry is straightforward and ongoing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Automated Symbolic Upscaling: 2. Model Generation for Extended Applicability Regimes

Abstract In this second part of the two paper series, we detail an algorithmic procedure for systematically implementing the generalized closure form strategy presented in Part 1. This strategy extends the applicability of homogenized models with respect to classical homogenization theory, as demonstrated in Part 1 where upscaled models are rigorously derived in moderately reactive physical regimes. After encoding the algorithm into Symbolica, an automated upscaling framework, we upscale two reactive mass transport problems and numerically validate the resulting nonlinear homogenized models by showing the absolute error estimates predicted by homogenization theory are satisfied. In both problems, nontrivial closure forms and closure problems are automatically formulated using the encoded strategy with no human interaction, nor prior knowledge regarding the closure required for the systems. We hope these demonstrations spark further interest in automated analytical frameworks for multiscale modeling, as such capabilities are invaluable for generating rigorous multiscale models of complex phenomena in porous media.

Pietrzyk, Kyle↗
DOE OSTI2024

Molecular Dynamics Simulation of Complex Reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) Software Package

Simulating chemically reactive phenomena such as proton transport on nanosecond to microsecond and beyond time scales is a challenging task. Ab initio methods are unable to currently access these time scales routinely, and traditional molecular dynamics methods feature fixed bonding arrangements that cannot account for changes in the system’s bonding topology. The Multiscale Reactive Molecular Dynamics (MS-RMD) method, as implemented in the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) software package for the LAMMPS molecular dynamics code, offers a method to routinely sample longer time scale reactive simulation data with statistical precision. RAPTOR may also be interfaced with enhanced sampling methods to drive simulations toward the analysis of reactive rare events, and a number of collective variables (CVs) have been developed to facilitate this. Key advances to this methodology, including GPU acceleration efforts and novel CVs to model water wire formation are reviewed, along with recent applications of the method which demonstrate its versatility and robustness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗