Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reaction diffusion systems”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

SBML Level 3: an extensible format for the exchange and reuse of biological models

Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction-based models and packages that extend the core with features suited to other model types including constraint-based models, reaction-diffusion models, logical network models, and rule-based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single-cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Micro-continuum approach for mineral precipitation

Abstract Rates and extents of mineral precipitation in porous media are difficult to predict, in part because laboratory experiments are problematic. It is similarly challenging to implement numerical methods that model this process due to the need to dynamically evolve the interface of solid material. We developed a multiphase solver that implements a micro-continuum simulation approach based on the Darcy–Brinkman–Stokes equation to study mineral precipitation. We used the volume-of-fluid technique in sharp interface implementation to capture the propagation of the solid mineral surface. Additionally, we utilize an adaptive mesh refinement method to improve the resolution of near interface simulation domain dynamically. The developed solver was validated against both analytical solution and Arbitrary Lagrangian–Eulerian approach to ensure its accuracy on simulating the propagation of the solid interface. The precipitation of barite (BaSO 4 ) was chosen as a model system to test the solver using variety of simulation parameters: different geometrical constraints, flow conditions, reaction rate and ion diffusion. The growth of a single barite crystal was simulated to demonstrate the solver’s capability to capture the crystal face specific directional growth.

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

GP-BayesOpInf

SAND2025-01851O GP-BayesOpInf is a software tool that uses algorithms to combine Gaussian process regression, principal component analysis, and linear Bayesian inference to produce a probabilistic reduced-order model for time-dependent systems. Numerical examples include the compressible Euler equations for an ideal gas, a heat diffusion process with a nonlinear reaction term, and a set of ordinary differential equations describing a compartmental model in epidemiology. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗
NASA NTRS1988

Preparation for microgravity - The role of the Microgravity Material Science Laboratory

Experiments at the NASA Lewis Research Center's Microgravity Material Science Laboratory using physical and mathematical models to delineate the effects of gravity on processes of scientific and commercial interest are discussed. Where possible, transparent model systems are used to visually track convection, settling, crystal growth, phase separation, agglomeration, vapor transport, diffusive flow, and polymer reactions. Materials studied include metals, alloys, salts, glasses, ceramics, and polymers. Specific technologies discussed include the General Purpose furnace used in the study of metals and crystal growth, the isothermal dendrite growth apparatus, the electromagnetic levitator/instrumented drop tube, the high temperature directional solidification furnace, the ceramics and polymer laboratories and the center's computing facilities.

Johnston, J. Christopher↗
NASA NTRS1985

Amorphous metallic films in silicon metallization systems

Diffusion barrier research was focussed on lowering the chemical reactivity of amorphous thin films on silicon. An additional area of concern is the reaction with metal overlays such as aluminum, silver, and gold. Gold was included to allow for technology transfer to gallium arsenide PV cells. Amorphous tungsten nitride films have shown much promise. Stability to annealing temperatures of 700, 800, and 550 C were achieved for overlays of silver, gold, and aluminum, respectively. The lower results for aluminum were not surprising because there is an eutectic that can form at a lower temperature. It seems that titanium and zirconium will remove the nitrogen from a tungsten nitride amorphous film and render it unstable. Other variables of research interest were substrate bias and base pressure during sputtering.

So, F.↗
NASA NTRS1975

International Symposium on Combustion, 15th, Tokyo, Japan, August 25-31, 1974, Proceedings

Recent theoretical and experimental studies concerned with detonation and pressure wave combustion interaction, fire and explosion research and safety, heterogeneous combustion, flame-flow interactions, kinetics of elementary reactions, pollution control in and by combustion systems, and ignition are presented. Some of the topics covered include critical power density for direct initiation of unconfined gaseous detonations, extinction of laminar diffusion flames for liquid fuels, combustion of bulk titanium in oxygen, flame propagation in small spheres of unconfined and slightly confined flammable mixtures, kinetics of the reaction of nitric oxide with hydrogen, production of chemi-ions and formation of CH and CH2 radicals in methane-oxygen and ethylene-oxygen flames, NOx emission characteristics in two-stage combustion, and spherical ignition of oxyhydrogen behind a reflected shock wave. Individual items are announced in this issue.

Source record↗
NASA NTRS2013

Calcium-Magnesium-Aluminosilicate (CMAS) Reactions and Degradation Mechanisms of Advanced Environmental Barrier Coatings

The thermochemical reactions between calcium-magnesium-aluminosilicate- (CMAS-) based road sand and several advanced turbine engine environmental barrier coating (EBC) materials were studied. The phase stability, reaction kinetics and degradation mechanisms of rare earth (RE)-silicates Yb2SiO5, Y2Si2O7, and RE-oxide doped HfO2 and ZrO2 under the CMAS infiltration condition at 1500 C were investigated, and the microstructure and phase characteristics of CMAS-EBC specimens were examined using Scanning Electron Microscopy (SEM) and X-ray Diffraction (XRD). Experimental results showed that the CMAS dissolved RE-silicates to form crystalline, highly non-stoichiometric apatite phases, and in particular attacking the silicate grain boundaries. Cross-section images show that the CMAS reacted with specimens and deeply penetrated into the EBC grain boundaries and formed extensive low-melting eutectic phases, causing grain boundary recession with increasing testing time in the silicate materials. The preliminary results also showed that CMAS reactions also formed low melting grain boundary phases in the higher concentration RE-oxide doped HfO2 systems. The effect of the test temperature on CMAS reactions of the EBC materials will also be discussed. The faster diffusion exhibited by apatite and RE-doped oxide phases and the formation of extensive grain boundary low-melting phases may limit the CMAS resistance of some of the environmental barrier coatings at high temperatures.

Ahlborg, Nadia L.↗
DOE OSTI2020

Boosting the sodium storage behaviors of carbon materials in ether-based electrolyte through the artificial manipulation of microstructure

The porous carbon blacks rationally designed by a facile yet efficient NH 3 thermal etching route have been investigated as anode materials in an ether-based electrolyte for sodium-ion batteries. The as-synthesized CBN35 carbon black with a 35% weight loss after NH3 thermal etching exhibited a large specific charge capacity of 352 mAh g -1 at 50 mA g -1 and a superior rate capability of 101 mAh g -1 at 16000 mA g -1 , due to its highest microporosity, an appropriate surface area, a desirable microstructure, and a promising hybrid intercalation mechanism. Impressively, even cycled at 1600 mA g -1 over 3200 cycles, an outstanding reversible capacity of 103 mAh g -1 with a negligible 0.0162% capacity loss per cycle can still be achieved. Based on the multimodal characterizations including the structural probes of phase evolution for carbon materials, the electrochemical techniques, and the surface-sensitive XAS measurements, the exceptional electrochemical properties should stem from several merits of modified carbon black system. While the particular microporous structure provides relatively more accessible sodium storage sites, a novel hybrid intercalation mechanism in ether-based electrolyte would incorporate the sodium ion insertion into the disordered structure with the solvated sodium ion species co-intercalation into the graphitic phase. In addition to the diffusion-controlled redox reactions, the noticeable surface-induced pseudocapacitive reactions also significantly contribute to the charge storage upon sodiation and guarantee the rapid migrations of sodium ions/solvated compounds. In conclusion, this system further features a controlled emergence of a robust but thin solid electrolyte interphase layer, which could suppress the side reactions of active electrode with reactive electrolyte, maintain the fragile porous structure upon cycling, and facilitate the migrations of sodium ions and solvated sodium ion compounds.

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

Symmetrization of conservation laws with entropy for high-temperature hypersonic computations

Results of Hughes, France, and Mallet are generalized to conservation law systems taking into account high-temperature effects. Symmetric forms of different equation sets are derived in terms of entropy variables. First, the case of a general divariant gas is studied; it can be specialized to the usual Navier-Stokes equations, as well as to situations where the gas is vibrationally excited, and undergoes equilibrium chemical reactions. The case of gas in thermochemical nonequilibrium is considered next. Transport phenomena, and in particular mass diffusion, are examined in the framework of symmetric advective-diffusive systems.

Chalot, F.↗
DOE OSTI2022

A pseudo-two-dimensional (P2D) model for FeS 2 conversion cathode batteries

Conversion cathode materials are gaining interest for secondary batteries due to their high theoretical energy and power density. However, practical application as a secondary battery material is currently limited by practical issues such as poor cyclability. To better understand these materials, we have, for this study, developed a pseudo-two-dimensional model for conversion cathodes. We apply this model to FeS 2 – a material that undergoes intercalation followed by conversion during discharge. The model is derived from the half-cell Doyle–Fuller–Newman model with additional loss terms added to reflect the converted shell resistance as the reaction progresses. We also account for polydisperse active material particles by incorporating a variable active surface area and effective particle radius. Using the model, we show that the leading loss mechanisms for FeS 2 are associated with solid-state diffusion and electrical transport limitations through the converted shell material. The polydisperse simulations are also compared to a monodisperse system, and we show that polydispersity has very little effect on the intercalation behavior yet leads to capacity loss during the conversion reaction. Finally, we provide the code as an open-source Python Battery Mathematical Modeling (PyBaMM) model that can be used to identify performance limitations for other conversion cathode materials.

25 ENERGY STORAGE↗
NASA NTRS2002

Laminar Soot Processes (LSP)

This is the final report of a research program considering the structure and the soot surface reaction properties of laminar nonpremixed (diffusion) flames. The study was limited to ground-based measurements of buoyant laminar jet diffusion flames at pressures of 0.1-1.0 atm. The motivation for the research is that soot formation in flames is a major unresolved problem of combustion science that influences the pollutant emissions, durability and performance of power and propulsion systems, as well as the potential for developing computational combustion. The investigation was divided into two phases considering the structure of laminar soot-containing diffusion flames and the soot surface reaction properties (soot surface growth and oxidation) of these flames, in turn. The first phase of the research addressed flame and soot structure properties of buoyant laminar jet diffusion flames at various pressures. The measurements showed that H, OH and O radical concentrations were generally in superequilibrium concentrations at atmospheric pressure but tended toward subequilibrium concentrations as pressures decreased. The measurements indicated that the original fuel decomposed into more robust compounds at elevated temperatures, such as acetylene (unless the original fuel was acetylene) and H, which are the major reactants for soot surface growth, and that the main effect of the parent fuel on soot surface growth involved its yield of acetylene and H for present test conditions. The second phase of the research addressed soot surface reaction properties, e.g., soot surface growth and surface oxidation. It was found that soot surface growth rates in both laminar premixed and diffusion flames were in good agreement, that these rates were relatively independent of fuel type, and that these rates could be correlated by the Hydrogen-Abstraction/Carbon-Addition (HACA) mechanisms of Colket and Hall (1994), Frenklach et al. (1990,1994), and Kazakov et al. (1995). It was also found that soot surface oxidation rates were relatively independent of fuel type, were not correlated with O2, CO2, H2O and O collision rates but were correlated with the collision rates of OH with a collision efficiency of 0.14, in agreement with the early measurements in premixed flames of Neoh et al. (1980), after allowing for oxidation by O2 via the classical rate expression of Nagle and Strickland-Constable (1962).

Dai, Z.↗
DOE OSTI2023

MFiX: Fractional-Step Method Implementation

A comprehensive, multiphase computational fluid dynamics (CFD) simulation solves several coupled transport equations including continuity, momentum, species, and energy. Chemical reactions further couple these equations through heats of reaction and rates of formation of products and rates of destruction of reactants. A fractional-step method separates changes attributed to chemical reactions from transport phenomena like convection and diffusion. When the governing equations are split into the transport and reacting components, efficient and independent methodologies can be exploited to solve the different systems. Specifically, discretization of field variable transport equations results in large, sparse matrices which are loosely coupled. These systems are solved in succession using iterative techniques that take advantage of the matrix structure. In contrast, chemical reactions tightly couple field variables locally within the domain (e.g., within a single computational cell) resulting in low dimensional but dense, nonlinear systems that are better solved using direct integration techniques.

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

Time-lapse approach to correct deficiencies of 2D soil zymography

Membrane zymography is commonly used in rhizosphere ecology to map enzyme activities in intact soil samples and plant roots. The method consists of incubating a membrane saturated with an enzyme-specific fluorogenic substrate on the soil/root surface followed by measurements of fluorescence intensity of the product in the membrane. The traditional zymography is based on assumptions that fluorescence on membrane images is linearly increasing with time during zymography and an increase in product content is numerically equal to enzyme activity in soil below the membrane. These assumptions are unlikely to hold in experimental settings. In this work we introduce a new zymography technique, time-lapse zymography (TLZ); the approach that eliminates the need for assumptions of the traditional zymography and provides more realistic estimates of enzymatic activities. We assessed the new technique in a series of laboratory and modeling experiments, including quantification of the fluorescent product diffusion (e.g. MUF: 4-methylumbelliferone) from plant roots, enzyme activity measurements on the roots, and HYDRUS-2D & HP2 software calibration with obtained data. The calibrated model was used to analyze the processes governing spatial and temporal dynamics of MUF contents in the membrane. The results indicated that the enzyme diffusion within the membrane-soil system was negligible, and measured zymograms were adequately reproduced solely by accounting for substrate and product diffusions and for catalytic enzyme reaction described by the Michaelis-Menten equation. TLZ enabled identifying and using linear parts on MUF time series and accounting for MUF losses in each zymogram pixel, considerably improving the accuracy as compared to the traditional method. Results demonstrated that enzymatic activity from only a thin soil layer (~0.2 mm) is adequately represented in zymograms.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1988

Preparation for microgravity: The role of the microgravity materials science laboratory

A laboratory dedicated to ground based materials processing in preparation for space flight was established at the NASA Lewis Research Center. Experiments are performed to delineate the effects of gravity on processes of both scientific and commercial interest. Processes are modeled physically and mathematically. Transport model systems are used where possible to visually track convection, settling, crystal growth, phase separation, agglomeration, vapor transport, diffusive flow, and polymers reactions. The laboratory contains apparatus which functionally duplicates apparatus available for flight experiments and other pieces instrumented specifically to allow process characterization. Materials addressed include metals, alloys, salts, glasses, ceramics, and polymers. The Microgravity Materials Science Laboratory is staffed by engineers and technicians from a variety of disciplines and is open to users from industry and academia as well as the government. Examples will be given of the laboratory apparatus typical experiments and results.

Johnston, J. Christopher↗
DOE OSTI2024

Geochemical-hydromechanical couplings during water-serpentinized harzburgite interactions at 20° C and 30 bars

Interaction of groundwater with the serpentinized harzburgites of the upper mantle plays an important role in the geologic carbon cycle and serpentinization processes. Here, in this study, an intact serpentinized harzburgite with more than 60% lizardite from the Semail ophiolite in Oman was reacted with water at near neutral pH and low PCO 2 to observe its dissolution behavior and reaction path within the CaO-MgO-SiO 2 -CO 2 -H 2 O system at 20 °C and 30 bars. In the experiment, mass (i.e. H 2 O, ions) transfer is by diffusion and the selected temperature and pressure conditions mimic water-rock interactions at the shallow portions of most ultramafic rock hosted aquifers. Equilibrium activity-activity diagrams and thermodynamic data were used to compare the experimentally determined reaction path with the theoretical stability boundaries for minerals involved to investigate the evolution of bulk chemical composition, mineralogy, and water composition. The experiments demonstrated that (i) the water-rock interactions increased the pH of the aqueous phase from 5.9 to 7.5 over a period of 18 weeks; (ii) the dissolution of lizardite, orthopyroxene (enstatite), and clinopyroxene (diopside and augite) increased the concentration of Mg, Ca, and Si in the aqueous phase; (iii) the dissolution was incongruent with respect to Mg and Si, favoring Si release at pH > 6 due to (a) breaking of more reactive cation‑oxygen bonds that is consistent with the stoichiometry of magnesium for proton exchange reaction which favors a surface that is enriched in Mg at basic conditions, and (b) preferential dissolution of clinopyroxene Mg 1.30 Ca 0.52 Fe 0.06 Si 2 O 6 ; (iv) the aqueous phase was undersaturated with respect to carbonate (e.g. calcite, magnesite, and hydromagnesite) and hydrous (e.g. lizardite, chrysotile, brucite, and talc) minerals; (v) the bulk chemical composition and mineralogy of the intact serpentinized harzburgite matrix did not change; (vi) the thermodynamic data can successfully predict water-intact serpentinized harzburgite behavior if the water chemical composition can be constrained; and (vii) no detectable macroscopic form of damage (e.g. microcracks forming as a result of differential stress fields due to changes in volume during chemical reaction) was observed.

Carbonation↗
NASA NTRS1994

Soot Formation in Hydrocarbon/Air Laminar Jet Diffusion Flames

Soot processes within hydrocarbon/air diffusion flames are important because they affect the durability and performance of propulsion systems, the hazards of unwanted fires, the pollutant and particulate emissions from combustion processes, and the potential for developing computational combustion. Motivated by these observations, this investigation involved an experimental study of the structure and soot properties of round laminar jet diffusion flames, seeking an improved understanding of soot formation (growth and nucleation) within diffusion flames. The present study extends earlier work in this laboratory concerning laminar smoke points (l) and soot formation in acetylene/air laminar jet diffusion flames (2), emphasizing soot formation in hydrocarbon/air laminar jet diffusion flames for fuels other than acetylene. In the flame system, acetylene is the dominant gas species in the soot formation region and both nucleation and growth were successfully attributed to first-order reactions of acetylene, with nucleation exhibiting an activation energy of 32 kcal/gmol while growth involved negligible activation energy and a collision efficiency of O.53%. In addition, soot growth in the acetylene diffusion flames was comparable to new soot in premixed flame (which also has been attributed to first-order acetylene reactions). In view of this status, a major issue is the nature of soot formation processes in diffusion flame involving hydrocarbon fuels other than acetylene. In particular, information is needed about th dominant gas species in the soot formation region and the impact of gas species other than acetylene on soot nucleation and growth.

Sunderland, P. B.↗
DOE OSTI2023

Julia as a unifying end-to-end workflow language on the Frontier exascale system

We evaluate Julia as a single language and ecosystem paradigm powered by LLVM to develop workflow components for high-performance computing. We run a Gray-Scott, 2-variable diffusion-reaction application using a memory-bound, 7-point stencil kernel on Frontier, the US Department of Energy’s first exascale supercomputer. We evaluate the performance, scaling, and trade-offs of (i) the computational kernel on AMD’s MI250x GPUs, (ii) weak scaling up to 4,096 MPI processes/GPUs or 512 nodes, (iii) parallel I/O writes using the ADIOS2 library bindings, and (iv) Jupyter Notebooks for interactive analysis. Results suggest that although Julia generates a reasonable LLVM-IR, a nearly 50% performance difference exists vs. native AMD HIP stencil codes when running on the GPUs. As expected, we observed near-zero overhead when using MPI and parallel I/O bindings for system-wide installed implementations. Consequently, Julia emerges as a compelling high-performance and high-productivity workflow composition language, as measured on the fastest supercomputer in the world.

Godoy, William↗