Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion−reaction modeling”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A mineral precipitation model based on the volume of fluid method

A novel volume of fluid method is presented for mineral precipitation coupled with fluid flow and reactive transport. The approach describes the fluid-solid interface as a smooth transitional region, which is designed to provide the same precipitation rate and viscous drag force as a sharp interface. Specifically, the governing equation of mineral precipitation is discretized by an upwind scheme, and a rigorous effective viscosity model is derived around the interface. The model is validated against analytical solutions for mineral precipitation in channel and ring-shaped structures. It also compares well with interface tracking simulations of advection-diffusion-reaction problems. Here, the methodology is finally employed to model mineral precipitation in fracture networks, which is challenging due to the low porosity and complex geometry. Compared to other approaches, the proposed model has a concise algorithm and contains no free parameters. In the modeling, only the pore space requires meshing, which improves the computational effciency especially for low-porosity media.

58 GEOSCIENCES↗
DOE OSTI2022

Reaction–drift–diffusion models from master equations: application to material defects

We present a general method to produce well-conditioned continuum reaction–drift–diffusion equations directly from master equations on a discrete, periodic state space. We assume the underlying data to be kinetic Monte Carlo models (i.e. continuous-time Markov chains) produced from atomic sampling of point defects in locally periodic environments, such as perfect lattices, ordered surface structures or dislocation cores, possibly under the influence of a slowly varying external field. Our approach also applies to any discrete, periodic Markov chain. Here, the analysis identifies a previously omitted non-equilibrium drift term, present even in the absence of external forces, which can compete in magnitude with the reaction rates, thus being essential to correctly capture the kinetics. To remove fast modes which hinder time integration, we use a generalized Bloch relation to efficiently calculate the eigenspectrum of the master equation. A well conditioned continuum equation then emerges by searching for spectral gaps in the long wavelength limit, using an established kinetic clustering algorithm to define a proper reduced, Markovian state space.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Bayesian learning with Gaussian processes for low-dimensional representations of time-dependent nonlinear systems

This work presents a data-driven method for learning low-dimensional time-dependent physics-based surrogate models whose predictions are endowed with uncertainty estimates. We use the operator inference approach to model reduction that poses the problem of learning low-dimensional model terms as a regression of state space data and corresponding time derivatives by minimizing the residual of reduced system equations. Standard operator inference models perform well with accurate training data that are dense in time, but producing stable and accurate models when the state data are noisy and/or sparse in time remains a challenge. Another challenge is the lack of uncertainty estimation for the predictions from the operator inference models. Our approach addresses these challenges by incorporating Gaussian process surrogates into the operator inference framework to (1) probabilistically describe uncertainties in the state predictions and (2) procure analytical time derivative estimates with quantified uncertainties. The formulation leads to a generalized least-squares regression and, ultimately, reduced-order models that are described probabilistically with a closed-form expression for the posterior distribution of the operators. The resulting probabilistic surrogate model propagates uncertainties from the observed state data to reduced-order predictions. Furthermore, we demonstrate the method is effective for constructing low-dimensional models of two nonlinear partial differential equations representing a compressible flow and a nonlinear diffusion–reaction process, as well as for estimating the parameters of a low-dimensional system of nonlinear ordinary differential equations representing compartmental models in epidemiology.

Data-driven model reduction↗
DOE OSTI2022

Diffusion-limited hydrolysis in polymeric materials

Polymer degradation under aggressive environmental stressors often develops heterogeneities due to diffusion-limited reaction phenomena. This is well established for diffusion-limited oxidation (DLO), which is known to occur for most polymeric materials at elevated temperatures but has been less summarized for the conditions of diffusion-limited hydrolysis (DLH). Here, an overview of hydrolysis for several materials and a computational model, analogous to the underlying equations for DLO, able to define this diffusion-reaction system is presented. A systematic study of the influence of various parameters, such as water diffusivity, reaction rate and order, and a more in-depth focus on residual isocyanate hydrolysis in a polymeric methylene diphenyl diisocyanate (pMDI) based polyurethane (PU) foam is given. For this system, we present experimental data for model ‘input’ parameters and discuss predictions for different conditions. We conceptually compare the behavior of diffusion-limited oxidation to that of diffusion-limited hydrolysis (DLH). With the mathematical framework and key material properties presented herein, any DLH phenomena following Fickian diffusion behavior can be understood, modeled, and predicted.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Assimilating partial observation to enhance feedback control of stochastic dynamical systems

Here, in this paper, we present a novel methodology to tackle feedback optimal control problems in scenarios where the exact state of the controlled process is unknown. It integrates data assimilation techniques and optimal control solvers to manage partial observation of the state process, a common occurrence in practical scenarios. Traditional stochastic optimal control methods assume full state observation, which is often not feasible in real-world fluid dynamics control problems. Our approach underscores the significance of utilizing observational data to inform control policy design. Specifically, we introduce a kernel learning backward stochastic differential equation (SDE) filter to enhance data assimilation efficiency and propose a sample-wise stochastic optimization method within the stochastic maximum principle framework. We demonstrate the efficacy and accuracy of our method in the control of advection-diffusion-reaction flow problem and the Dubins airplane maneuvering problem with model uncertainty.

data driven↗
DOE OSTI2020

Defect cluster and nonequilibrium gas bubble associated growth in irradiated UMo fuels – a cluster dynamics and phase field model

Irradiation examination shows that gas bubble swelling kinetics is much faster after irradiation-induced recrystallization than that prior recrystallization in UMo fuels. It implies that gas bubbles in coarse grains and small recrystallized grains have different growth behavior. In this work, a phase-field model of gas bubble evolution integrating microstructure dependent cluster dynamics has been developed, for the first time, to study the gas bubble swelling behavior in the recrystallization zone of UMo fuels. Generation, diffusion, reaction, sink, emission and clustering of vacancies and interstitials are described by the cluster dynamics model while a phase-field model is used to describe the evolution of non-equilibrium gas bubbles including nucleation and growth. With the coupled model, the effect of defect generation rate, clustering rate, interstitial emission and sink rates on grain boundaries on the gas bubble evolution are systematically simulated. A set of model parameters (defect generation rate, clustering rate, interstitial emission and sink rates) is determined by comparing measured and simulated gas bubble swelling kinetics. The results demonstrate that interstitial clustering is one of the important physical mechanisms which results in a fast gas bubble swelling kinetics in the recrystallization zone. The developed model can also be extended to study the associated growth of defect and second phase precipitates often observed in irradiated materials.

Hu, Shenyang↗
DOE OSTI2023

CompLaB v1.0: a scalable pore-scale model for flow, biogeochemistry, microbial metabolism, and biofilm dynamics

Abstract. Microbial activity and chemical reactions in porous media depend on the local conditions at the pore scale and can involve complex feedback with fluid flow and mass transport. We present a modeling framework that quantitatively accounts for the interactions between the bio(geo)chemical and physical processes and that can integrate genome-scale microbial metabolic information into a dynamically changing, spatially explicit representation of environmental conditions. The model couples a lattice Boltzmann implementation of Navier–Stokes (flow) and advection–diffusion-reaction (mass conservation) equations. Reaction formulations can include both kinetic rate expressions and flux balance analysis, thereby integrating reactive transport modeling and systems biology. We also show that the use of surrogate models such as neural network representations of in silico cell models can speed up computations significantly, facilitating applications to complex environmental systems. Parallelization enables simulations that resolve heterogeneity at multiple scales, and a cellular automaton module provides additional capabilities to simulate biofilm dynamics. The code thus constitutes a platform suitable for a range of environmental, engineering and – potentially – medical applications, in particular ones that involve the simulation of microbial dynamics.

58 GEOSCIENCES↗
NASA NTRS2016

Materials Outgassing Rate Decay in Vacuum at Isothermal Conditions

As a laboratory for scientific research, the International Space Station has been in Low Earth Orbit for nearly 20 years and is expected to be on-orbit for another 10 years. The ISS has been maintaining a relatively pristine contamination environment for science payloads. Materials outgassing induced contamination is currently the dominant source for sensitive surfaces on ISS and modeling the outgassing rate decay over a 20 to 30 year period is challenging. Materials outgassing is described herein as a diffusion-reaction process using ASTM E 1559 rate data. The observation of -1/2 (diffusion) or non-integers (reaction limited) as rate decay exponents for common ISS materials indicate classical reaction kinetics is unsatisfactory in modeling materials outgassing. Non-randomness of reactant concentrations at the interface is the source of this deviation from classical reaction kinetics. A diffusion limited decay was adopted as the result of the correlation of the contaminant layer thicknesses on returned ISS hardware, the existence of high outgassing silicone exhibiting near diffusion limited decay, and the confirmation of non-depleted material after ten years in the Low Earth Orbit.Keywords: Materials Outgassing, ASTM E 1559, Reaction Kinetics, Diffusion, Space Environments Effects, Contamination

Huang, Alvin Y.↗
NASA NTRS2014

RITRACKS: A Software for Simulation of Stochastic Radiation Track Structure, Micro and Nanodosimetry, Radiation Chemistry and DNA Damage for Heavy Ions

The code RITRACKS (Relativistic Ion Tracks) has been developed over the last few years at the NASA Johnson Space Center to simulate the effects of ionizing radiations at the microscopic scale, to understand the effects of space radiation at the biological level. The fundamental part of this code is the stochastic simulation of radiation track structure of heavy ions, an important component of space radiations. The code can calculate many relevant quantities such as the radial dose, voxel dose, and may also be used to calculate the dose in spherical and cylindrical targets of various sizes. Recently, we have incorporated DNA structure and damage simulations at the molecular scale in RITRACKS. The direct effect of radiations is simulated by introducing a slight modification of the existing particle transport algorithms, using the Binary‐Encounter‐Bethe model of ionization cross sections for each molecular orbitals of DNA. The simulation of radiation chemistry is done by a step‐by‐step diffusion‐reaction program based on the Green's functions of the diffusion equation]. This approach is also used to simulate the indirect effect of ionizing radiation on DNA. The software can be installed independently on PC and tablets using the Windows operating system and does not require any coding from the user. It includes a Graphic User Interface (GUI) and a 3D OpenGL visualization interface. The calculations are executed simultaneously (in parallel) on multiple CPUs. The main features of the software will be presented.

Plante, I↗
DOE OSTI2021

Analyzing the 2019 Chilean social outbreak: Modelling Latin American economies

In this work, we propose a quantitative model for the 2019 Chilean protests. We utilize public data for the consumer price index, the gross domestic product, and the employee and per capita income distributions as inputs for a nonlinear diffusion-reaction equation, the solutions to which provide an in-depth analysis of the population dynamics. Specifically, the per capita income distribution stands out as a solution to the extended Fisher-Kolmogorov equation. According to our results, the concavity of employee income distribution is a decisive input parameter and, in contrast to the distributions typically observed for Chile and other countries in Latin America, should ideally be non-negative. Based on the results of our model, we advocate for the implementation of social policies designed to stimulate social mobility by broadening the distribution of higher salaries.

Curilef, Sergio (ORCID:0000000283183206)↗
DOE OSTI2025

Nonintrusive projection-based reduced order modeling using stable learned differential operators

Nonintrusive projection-based reduced order models (ROMs) are essential for dynamics prediction in multi-query applications where underlying governing equations are known but the access to the source of the underlying full order model (FOM) is unavailable; that is, FOM is a glass-box. This article proposes a learn-then-project approach for nonintrusive model reduction. In the first step of this approach, high-dimensional stable sparse learned differential operators (S-LDOs) are determined using the generated data. In the second step, the ordinary differential equations, comprising these S-LDOs, are used with suitable dimensionality reduction and low-dimensional subspace projection methods to provide equations for the evolution of reduced states. This approach allows easy integration into the existing intrusive ROM framework to enable nonintrusive model reduction while allowing the use of Petrov–Galerkin projections. The applicability of the proposed approach is demonstrated for Galerkin and LSPG projection-based ROMs through four numerical experiments: 1-D scalar advection, 1-D Burgers, 2-D scalar advection and 1-D scalar advection–diffusion–reaction equations. In conclusion, the results indicate that the proposed nonintrusive ROM strategy provides accurate and stable dynamics prediction.

42 ENGINEERING↗
NASA NTRS1987

Optimal cure cycle design of a resin-fiber composite laminate

A unified computed aided design method was studied for the cure cycle design that incorporates an optimal design technique with the analytical model of a composite cure process. The preliminary results of using this proposed method for optimal cure cycle design are reported and discussed. The cure process of interest is the compression molding of a polyester which is described by a diffusion reaction system. The finite element method is employed to convert the initial boundary value problem into a set of first order differential equations which are solved simultaneously by the DE program. The equations for thermal design sensitivities are derived by using the direct differentiation method and are solved by the DE program. A recursive quadratic programming algorithm with an active set strategy called a linearization method is used to optimally design the cure cycle, subjected to the given design performance requirements. The difficulty of casting the cure cycle design process into a proper mathematical form is recognized. Various optimal design problems are formulated to address theses aspects. The optimal solutions of these formulations are compared and discussed.

Hou, Jean W.↗
DOE OSTI2021

Encapsulation of highly viscous CO 2 capture solvents for enhanced capture kinetics: Modeling investigation of mass transfer mechanisms

The encapsulation of highly viscous liquid-like Nanoparticle Organic Hybrid Materials (NOHMs) inside a gas permeable polymer to form SIPs (Solvent Impregnated Polymers) significantly enhanced the CO 2 capture kinetics of NOHMs, leading to a remarkable 50-fold increase in CO 2 flux compared to the neat NOHMs. To understand the mechanism for enhanced CO 2 mass transfer within these hybrid materials, kinetic modeling of CO 2 uptake into SIPs containing polyethylenimine functionalized NOHMs, denoted NPEI-SIP, was conducted. CO 2 mass transfer into NPEI-SIP films was found to conform, both qualitatively and quantitatively, with a diffusion-controlled moving front model. The diffusion-controlled model was also used to simulate CO 2 uptake within a fixed bed containing polydisperse NPEI-SIP particles, and this model accurately predicted experimentally measured breakthrough curves at 25 °C and 50 °C. The 50-fold increase in gas flux was shown to be a consequence of the very large CO 2 permeability within the polymer–solvent composite. Furthermore, the increase in gas flux is also dependent on the diffusion–reaction regime in which the chemical solvent operates, and the largest improvement will occur when immobilizing solvents, such as NOHMs, which operate in the instantaneous-reaction regime. The CO 2 capacity (~3 mol CO 2 /kg) and saturation time (~5 min) of 430 μm SIP particles were comparable to popular CO 2 chemisorption materials such as amine grafted silicates, in spite of the slow kinetics of NOHM-I-PEI in liquid form (CO 2 saturation time ~24 h for a 1 mm thin film).

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Generalized Grain-Scale Model for the Non-Plasma and Plasma-Assisted Hydrogen Direct Reduction of Iron Ore

Direct Reduction of Iron ore using hydrogen (H-DRI) is a promising pathway towards efficient steelmaking and accurate predictive models are a necessity for scale-up and optimization of this technology. However, accurate models of this process remain limited because existing models oversimplify grain-scale phenomena, such as nonlinearity inside grain, self-sufficient porosity, surface reactions, and the role of plasma species. These phenomena are important for flash steelmaking and plasma-assisted H-DRI processes. To address this need, we present a phenomenological model for simulating H-DRI at the scale of a single micron-sized grain of the iron ore. We call this the Transient Reactive Grain Model (TRGM). TRGM incorporates key physical process: gas species transport, a chemical kinetics of material conversion, nanopore structural evolution and, adsorption-desorption surface kinetics at the reactive nanopore surface. The important contribution of this work is that the model provides a dependence on different reductant species, specifically hydrogen atoms versus molecules, so that role of hydrogen plasma reduction can be clarified compared to the use of pure hydrogen gas reduction. TRGM predictions agree well with experimental data for both molecular H2 reduction of Fe2O3 and plasma hydrogen reduction of Fe3O4. Results reveal species concentration gradients with a diffuse reaction zone, and enhanced hydrogen diffusion at the grain outer surface due to evolving porosity. These findings challenge common assumptions in existing models, including sharp reaction fronts, quasi-steady diffusion and kinetics, and the neglect of surface chemistry. As a generalized grain-scale model for H-DRI processes, TRGM has practical applications in flash steelmaking and in-flight reduction using both molecular and plasma hydrogen.

08 HYDROGEN↗
DOE OSTI2023

On the influence of over-parameterization in manifold based surrogates and deep neural operators

Constructing accurate and generalizable approximators (surrogate models) for complex physico-chemical processes exhibiting highly non-smooth dynamics is challenging. The main question is what type of surrogate models we should construct and should these models be under-parameterized or over-parameterized. In this work, we propose new developments and perform comparisons for two promising approaches: manifold-based polynomial chaos expansion (m-PCE) and the deep neural operator (DeepONet), and we examine the effect of over-parameterization on generalization. While m-PCE enables the construction of a mapping by first identifying low-dimensional embeddings of the input functions, parameters, and quantities of interest (QoIs), a neural operator learns the nonlinear mapping via the use of deep neural networks. Here, we demonstrate the performance of these methods in terms of generalization accuracy by solving the 2D time-dependent Brusselator reaction-diffusion system with uncertainty sources, modeling an autocatalytic chemical reaction between two species. We first propose an extension of the m-PCE by constructing a mapping between latent spaces formed by two separate embeddings of the input functions and the output QoIs. To further enhance the accuracy of the DeepONet, we introduce weight self-adaptivity in the loss function. We demonstrate that the performance of m-PCE and DeepONet is comparable for cases of relatively smooth input-output mappings. However, when highly non-smooth dynamics is considered, DeepONet shows higher approximation accuracy. We also find that for m-PCE, modest over-parameterization leads to better generalization, both within and outside of distribution, whereas aggressive over-parameterization leads to over-fitting. In contrast, an even highly over-parameterized DeepONet leads to better generalization for both smooth and non-smooth dynamics. Furthermore, we compare the performance of the above models with another recently proposed operator learning model, the Fourier Neural Operator, and show that its over-parameterization also leads to better generalization. Taken together, our studies show that m-PCE can provide very good accuracy at very low training cost, whereas a highly over-parameterized DeepONet can provide better accuracy and robustness to noise but at higher training cost. In both methods, the inference cost is negligible.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

DRACO: An Overview [Slides]

DRACO (Diffusion ReACtiOn) is a diffusion and chemistry code designed to: 1) Operate on 3D with an unstructured grid defining an arbitrary geometry of interacting parts. 2) Generate its own meshes and use meshes created by other software. 3)Model the transport of any number of diffusing quantities: Concentrations, pressures, temperature, etc. 4) Allow diffusion coefficients to depend in an arbitrary way on concentration, temperature, position, time, etc. 5) Model general chemistry between concentrations with arbitrary reaction rates. 6) Allow arbitrary initial conditions, boundary conditions, and sources/sinks. 7) Allow all of the above to be specified by the user.

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

Direct simulations of chemically reacting turbulent mixing layers

The report presents the results of direct numerical simulations of chemically reacting turbulent mixing layers. The work consists of two parts: (1) the development and testing of a spectral numerical computer code that treats the diffusion reaction equations; and (2) the simulation of a series of cases of chemical reactions occurring on mixing layers. The reaction considered is a binary, irreversible reaction with no heat release. The reacting species are nonpremixed. The results of the numerical tests indicate that the high accuracy of the spectral methods observed for rigid body rotation are also obtained when diffusion, reaction, and more complex flows are considered. In the simulations, the effects of vortex rollup and smaller scale turbulence on the overall reaction rates are investigated. The simulation results are found to be in approximate agreement with similarity theory. Comparisons of simulation results with certain modeling hypotheses indicate limitations in these hypotheses. The nondimensional product thickness computed from the simulations is compared with laboratory values and is found to be in reasonable agreement, especially since there are no adjustable constants in the method.

Riley, J. J.↗