Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reaction-diffusion models”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Spiral waves are stable in discrete element models of two-dimensional homogeneous excitable media

The spontaneous breakup of a single spiral wave of excitation into a turbulent wave pattern has been observed in both discrete element models and continuous reaction-diffusion models of spatially homogeneous 2D excitable media. These results have attracted considerable interest, since spiral breakup is thought to be an important mechanism of transition from the heart rhythm disturbance ventricular tachycardia to the fatal arrhythmia ventricular fibrillation. It is not known whether this process can occur in the absence of disease-induced spatial heterogeneity of the electrical properties of the ventricular tissue. Candidate mechanisms for spiral breakup in uniform 2D media have emerged, but the physical validity of the mechanisms and their applicability to myocardium require further scrutiny. In this letter, we examine the computer simulation results obtained in two discrete element models and show that the instability of each spiral is an artifact resulting from an unphysical dependence of wave speed on wave front curvature in the medium. We conclude that spiral breakup does not occur in these two models at the specified parameter values and that great care must be exercised in the representation of a continuous excitable medium via discrete elements.

NASA Discipline Cardiopulmonary↗
NASA NTRS1996

Effects of Material Degradation on the Structural Integrity of Composite Materials: Experimental Investigation and Modeling of High Temperature Degradation Mechanisms

It has previously been demonstrated that simple coupled reaction-diffusion models can approximate the aging behavior of PMR-15 resin subjected to different oxidative environments. Based on empirically observed phenomena, a model coupling chemical reactions, both thermal and oxidative, with diffusion of oxygen into the material bulk should allow simulation of the aging process. Through preliminary modeling techniques such as this it has become apparent that accurate analytical models cannot be created until the phenomena which cause the aging of these materials are quantified. An experimental program is currently underway to quantify all of the reaction/diffusion related mechanisms involved. The following contains a summary of the experimental data which has been collected through thermogravimetric analyses of neat PMR-15 resin, along with analytical predictions from models based on the empirical data. Thermogravimetric analyses were carried out in a number of different environments - nitrogen, air and oxygen. The nitrogen provides data for the purely thermal degradation mechanisms while those in air provide data for the coupled oxidative-thermal process. The intent here is to effectively subtract the nitrogen atmosphere data (assumed to represent only thermal reactions) from the air and oxygen atmosphere data to back-figure the purely oxidative reactions. Once purely oxidative (concentration dependent) reactions have been quantified it should then be possible to quantify the diffusion of oxygen into the material bulk.

Cunningham, Ronan A.↗
DOE OSTI2026

Reaction Diffusion Modelling of 3D Pillar Electrodes in Single-Catalyst CO 2 Reduction Cascades

Effective electrochemical CO 2 reduction to liquid fuels requires that the local catalytic environment facilitates the desired reactivity, yet a microscopic understanding of this environment is difficult to achieve from experiment alone. In this work, a 3D reaction-diffusion model was developed to explore the effects of electrode surface area and local geometry on the performance of a heterogeneous catalyst that performs a two-step CO 2 reduction cascade reaction to CO and then CH 3 OH under aqueous conditions. Kinetic parameters for the model were inspired by experimental results using a cobalt phthalocyanine (CoPc) catalyst. Three-dimensional architectures composed of arrays of square pillars with varying dimensions and either smooth or periodically modulated surfaces were tested, revealing the extent to which geometry modulates the performance of the cascade reactions. Although structural variations modulate local concentration gradients, we find that electrochemically active surface area predominantly governs the overall cascade reaction. Moreover, the results suggest that supersaturation of CO, with concentrations up to ten-fold higher than the equilibrium solubility limit, might be critical for more efficient conversion to CH 3 OH. For any given geometry, the spatially averaged ratio of [CO] to [CO 2 ] is dictated by the electrochemically active surface area and determines the yield of CH 3 OH. For a fixed surface area, geometries that spatially confine the electrolyte yield moderate local [CO] to [CO 2 ] ratios within small volumes. In contrast, less confining geometries result in a broader distribution of local ratios spread over larger volumes, with both configurations yielding the same spatially averaged [CO] to [CO 2 ] ratio. These insights provide valuable design principles—highlighting the critical importance of surface area and possibly CO supersaturation—for engineering advanced electrode architectures that leverage intermediate trapping and CO supersaturation to enhance overall performance in tandem CO 2 reduction systems.

COMSOL↗
DOE OSTI2025

Modeling diffusion and depletion in high-aspect-ratio atomic layer deposition processes: Process parameters and manufacturing impacts

Atomic layer deposition (ALD) is a powerful technique for modifying the surface chemistry and properties of substrates with complex and nonplanar topologies. However, achieving uniform and conformal deposition on ultrahigh-aspect-ratio substrates remains challenging, typically requiring large quantities of precursors and long exposure times. Furthermore, process optimization is often performed empirically and involves substantial trial and error. In this work, we perform a combined experimental and computational study of ALD Al 2 O 3 infiltration into silica aerogel monoliths (aspect ratio >10 5 ). A reaction-diffusion model is used to explore the effects of key processing parameters, namely, exposure time per dose, precursor source temperature, number of aerogels in the reactor, and reactor volume. The model is based on quasi-static mode ALD, where the dosed precursor is held in the chamber for a fixed period of time before purging. We analyze the trade-offs between process throughput and precursor utilization for each of these parameters. Furthermore, we investigate the co-optimization and interactions between multiple process parameters, demonstrating the potential for further improvements. Furthermore, this physics-based model can be used to identify a set of process parameters for high-aspect-ratio ALD that meet specific manufacturing objective functions, including throughput, cost, and sustainability.

Aerogel↗
DOE OSTI2025

Atomic layer deposition of nanofilms on porous polymer substrates: Strategies for success

Atomic layer deposition (ALD) is a versatile technique for engineering the surfaces of porous polymers, imbuing the flexible, high-surface-area substrates with inorganic and hybrid material properties. Previously reported enhancements include fouling resistance, electrical conductance, thermal stability, photocatalytic activity, hydrophilicity, and oleophilicity. However, there are many poorly understood phenomena that introduce challenges in applying ALD to porous polymers. In this paper, we address five common challenges and ways to overcome them: (1) entrapped precursor, (2) embrittlement, (3) film fracture, (4) deformation, and (5) pore collapse. These challenges are often interrelated and can exacerbate one another. To investigate these phenomena, we applied various ALD chemistries to porous polymers including polyethersulfone, polysulfone, polyvinylidene fluoride, and polycarbonate track-etched membranes. Reaction-diffusion modeling revealed why certain precursors and processing conditions result in embrittling subsurface material growth, entrapment of unreacted precursors, and nongrowth. We quantify the limits of ALD processing temperatures that are dictated by thermal expansion mismatch and can lead to fractured ALD films. The results herein allow us to make recommendations to avoid, mitigate, or overcome the difficulties encountered when performing ALD and plasma-enhanced ALD on porous polymers. We intend this article to serve as a “lessons learned” guide informed by previous experience to provide a better understanding of the difficulties and limitations of ALD on porous polymers and knowledge-based guidelines for successful depositions. This knowledge can accelerate future research and help experimentalists navigate and troubleshoot as they expose porous polymers to reactive precursor vapors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Sequential Infiltration Synthesis of Multinary Metal Oxides: The Role of Precursor Diffusion and Reactivity

Understanding volatile metal–organic precursor transport in polymer matrices during sequential infiltration synthesis (SIS) is critical to the design of hybrid materials of tailored composition and uniform structure. Here, we investigate the diffusion and reaction behavior of diethylzinc (DEZ) in poly(methyl methacrylate) (PMMA) matrices with variable indium and zinc incorporation. Unlike conventional SIS where precursors adduct or react directly with the polymer backbone, DEZ shows minimal interaction with PMMA. However, DEZ may interact strongly with In–OH or Zn–OH moieties that are previously incorporated into a PMMA thick film. Using spectroscopic ellipsometry, X-ray photoelectron spectroscopy depth profiling, and EDX line scans, we quantify Zn incorporation as a function of exposure time, polymer thickness, and previous infiltration steps. The Zn distribution evolves from surface-localized to deep within the ∼1 μm thick film with increasing DEZ exposures up to 30 min. A reaction-diffusion model incorporating a hindering factor to capture diffusivity decay due to product accumulation reproduces the Zn depth profile and yields physically meaningful parameters to guide future growth. Furthermore, this study provides a framework for interpreting metal–organic precursor infiltration behavior and highlights the utility of model-guided SIS process design.

diethylzinc↗
DOE OSTI2026

The Cascade Effectiveness of 3-Terminal Tandem Photocathode Architectures as Applied to CO 2 Reduction

Cascade catalysis for photoelectrochemical CO 2 reduction (CO 2 R) decouples the overall reaction into sequential steps occurring on separately optimized catalysts (for example, Ag and Cu) between which an intermediate species such as CO is transferred. A 3-terminal tandem (3TT) photovoltaic architecture advantageously holds two different catalytic regions at different potentials under a single illumination source, but its overall efficiency is low. Using a stochastic reaction-diffusion model, we have examined 3TT photocathode design principles, focusing on the coupling of surface chemistry and transport of CO both inside the boundary layer and outward, into the bulk electrolyte. We find that ensuring that the lateral diffusion distance within the boundary layer is short compared to the boundary layer thickness and controlling bulk flow are key, with interdigitated designs showing an overall conversion efficiency improvement by 2 orders of magnitude compared to the side-by-side case without flow. These findings can be adapted for other cascaded architectures.

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

Curvature effects in the reaction-diffusion system governing frontal polymerization

Frontal polymerization is a fast and energy-efficient curing process of thermoset polymers involving the self-sustained propagation of a sharp front that transforms a liquid monomer or gel into a polymer. The temperature, speed, and geometry of the polymerization front are closely intertwined quantities that can affect the curing process and the mechanical properties of the polymer. This work investigates the effect of the front curvature on the front propagation speed based on a reaction-diffusion model. The proposed analytical relations between front speed and curvature are verified by finite element simulations and validated through experiments. The models show that a concave front has a higher front temperature and speed than its convex counterpart. When the curvature of a convex front reaches a critical value, the reaction can no longer be sustained, and the front is quenched. The relation between the speed of a steadily propagating front and its geometry is derived, showing good agreement with experimental observations. This study also provides analytical estimates of the thickness of the polymerization front and insight into the relationship between the front curvature and the kinetics of the exothermic reaction.

42 ENGINEERING↗
DOE OSTI2025

Developing the Science Basis for Understanding Polymer Encapsulant Degradation Mechanisms: DuraMAT 2.0 Final Project Report

Polymeric encapsulants are essential materials in photovoltaic modules, protecting sensitive electronics from the environment while providing mechanical integrity to the multilayered assembly. However, these polymeric materials are susceptible to degradation processes driven by the ingress of environmental species, ultraviolet radiation, thermal stresses, and mechanical loading. In this study, we employ a combined atomistic simulation and accelerated aging experimental approach to study the molecular-scale mechanisms of encapsulant degradation. Classical molecular dynamics simulations quantify the diffusion of environmental and degradation species through the polymer matrix, producing composition-specific diffusion coefficients. Reactive simulations characterize activation energy barriers and reaction rate constants for key chemical pathways. In parallel, thermal-desorption analyses coupled with mass spectrometry monitor the emergence and concentration profiles of degradation products under controlled stressor conditions. By integrating simulation and experiment, we establish quantitative correlations between polymer composition, species diffusivity, and chemical reactivity. We anticipate that these relations and quantitative values could serve as high-fidelity inputs to reaction-diffusion models, enabling physics-informed lifetime predictions and guiding the design of more durable encapsulant materials for solar energy applications.

36 MATERIALS SCIENCE↗
NASA NTRS1983

Numerical solution of the problem of flame propagation by the use of the random element method

A numerical, grid-free algorithm is presented for one-dimensional reaction-diffusion model of laminar flame propagation in premixed gases. It is based on the random element method we developed for the analysis of diffusional processes. The effect of combustion is taken into account by applying the principle of fractional steps to separate the process of diffusion, modeled by the random walk of computational elements, from the exothermic effects of chemical reaction, monitoring their strength. The validity of the algorithm is demonstrated by application to flame propagation problems for which exact solutions exist. The flame speed evaluated by its use oscillates around the exact value at a relatively small amplitude, while the temperature and species concentration profiles are self-correcting in their convergence to the exact solution. A satisfactory resolution is obtained by the use of quite a small number of computational elements which automatically adjust their distribution of fit sharp gradients.

Ghoniem, A. F.↗
DOE OSTI2026

Efficient data-driven regression for reduced-order modeling of spatial pattern formation

We present an efficient data-driven regression approach for constructing reduced-order models (ROMs) of reaction-diffusion systems exhibiting pattern formation. The ROMs are learned non-intrusively from available training data of physically accurate numerical simulations. The method can be applied to general nonlinear systems through the use of polynomial model form, while not requiring knowledge of the underlying physical model, governing equations, or numerical solvers. The process of learning ROMs is posed as a low-cost least-squares problem in a reduced-order subspace identified via Proper Orthogonal Decomposition (POD). Numerical experiments on classical pattern-forming systems–including the Schnakenberg and Mimura–Tsujikawa models–demonstrate that higher-order surrogate models significantly improve prediction accuracy while maintaining low computational cost. The proposed method provides a flexible, non-intrusive model reduction framework, well suited for the analysis of complex spatio-temporal pattern formation phenomena.

Data-driven modeling↗
NASA NTRS1990

Computation of the unsteady facilitated transport of oxygen in hemoglobin

The transport of a reacting permeant diffusing through a thin membrane is extended to more realistic dissociation models. A new nonlinear analysis of the reaction-diffusion equations, using implicit finite-difference methods and direct block solvers, is used to study the limits of linearized and equilibrium theories. Computed curves of molecular oxygen permeating through hemoglobin solution are used to illustrate higher-order reaction models, the effect of concentration boundary layers at the membrane interfaces, and the transient buildup of oxygen flux.

Davis, Sanford↗
NASA NTRS1994

Microwave heating and joining of ceramic cylinders: A mathematical model

A thin cylindrical ceramic sample is placed in a single mode microwave applicator in such a way that the electric field strength is allowed to vary along its axis. The sample can either be a single rod or two rods butted together. We present a simple mathematical model which describes the microwave heating process. It is built on the assumption that the Biot number of the material is small, and that the electric field is known and uniform throughout the cylinder's cross-section. The model takes the form of a nonlinear parabolic equation of reaction-diffusion type, with a spatially varying reaction term that corresponds to the spatial variation of the electromagnetic field strength in the waveguide. The equation is analyzed and a solution is found which develops a hot spot near the center of the cylindrical sample and which then propagates outwards until it stabilizes. The propagation and stabilization phenomenon concentrates the microwave energy in a localized region about the center where elevated temperatures may be desirable.

Booty, Michael R.↗
DOE OSTI2025

Steric effects of central dogma processes on the compaction and segregation of bacterial nucleoids

The bacterial cytoplasm is characterized by a distinctive membrane-less organelle, the nucleoid, which harbors the chromosomal DNA. Here, we investigate the steric effects of dynamic processes associated with transcription and translation on the structure of this organelle using coarse-grained molecular dynamics simulations that incorporate out-of-equilibrium reactions. Our model captures the scale of the entire cell and incorporates a reaction-diffusion system for ribosomes and polyribosomes, coupling their nonequilibrium kinetics to DNA excluded volume interactions. Our findings demonstrate that out-of-equilibrium reactions increase the size of the nucleoid and the number of ribosomes in this subcellular region. In addition, we show that nucleoid size is proportional to transcriptional activity. Our model reproduces the time-dependent change in nucleoid size observed in rifampicin treatment experiments, where the pool of polyribosomes is depleted. Furthermore, we find that these active processes are essential for complete sister chromosome separation and correct nucleoid positioning within the cell. Overall, our study reveals the effects of the central dogma processes on the internal organization and localization of bacterial nucleoids.

Chang, Mu-Hung [Univ. of Tennessee, Knoxville, TN ↗
DOE OSTI2025

Gray-Scott reaction-diffusion

This dataset consists of 1000 records in a single HDF5 file generated from the simulation code available at: https://github.com/lezahlie/greyscott_simulation. More about the Gray-Scott model is explained here https://visualpde.com/nonlinear-physics/gray-scott.html.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗
NASA NTRS1987

Numerical methods for one-dimensional reaction-diffusion equations arising in combustion theory

A review of numerical methods for one-dimensional reaction-diffusion equations arising in combustion theory is presented. The methods reviewed include explicit, implicit, quasi-linearization, time linearization, operator-splitting, random walk and finite-element techniques and methods of lines. Adaptive and nonadaptive procedures are also reviewed. These techniques are applied first to solve two model problems which have exact traveling wave solutions with which the numerical results can be compared. This comparison is performed in terms of both the wave profile and computed wave speed. It is shown that the computed wave speed is not a good indicator of the accuracy of a particular method. A fourth-order time-linearized, Hermitian compact operator technique is found to be the most accurate method for a variety of time and space sizes.

Ramos, J. I.↗
DOE OSTI2025

Morphogenic Growth 3D Printing

Inspired by nature's morphogenesis, a new 3D printing process –growth printing (GP)– takes advantage of a self‐propagating curing front to produce 3D polymeric parts following a growth‐like development plan. The propagation of the curing front is driven by the exothermic polymerization of dicyclopentadiene (DCPD), which transforms the liquid resin into a stiff polymer as it propagates at 1 mm s −1 . GP is triggered when a heated initiator contacts the uncured liquid resin in an open container. The initiator nucleates the frontal polymerization reaction and the isotropic radial propagation of the growth front. Simultaneously, the initiator is moved up across the free surface of the resin, pulling the cured object out of the uncured resin. The motion trajectory of the initiator with respect to the free resin surface controls the growth morphology of the 3D part. An inverse design algorithm is developed to produce 3D parts by modeling the reaction‐diffusion‐driven solidification process. This process has substantial energy savings and high printing speeds.

3D printing↗