Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reaction-diffusion”

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
DOE OSTI2024

Reaction-diffusion modeling provides insights into biophysical carbon-concentrating mechanisms in land plants

Carbon-concentrating mechanisms (CCMs) have evolved numerous times in photosynthetic organisms. They elevate the concentration of CO 2 around the carbon-fixing enzyme rubisco, thereby increasing CO 2 assimilatory flux and reducing photorespiration. Biophysical CCMs, like the pyrenoid-based CCM (PCCM) of Chlamydomonas reinhardtii or carboxysome systems of cyanobacteria, are common in aquatic photosynthetic microbes, but in land plants appear only among the hornworts. To predict the likely efficiency of biophysical CCMs in C 3 plants, we used spatially resolved reaction-diffusion models to predict rubisco saturation and light use efficiency. We found that the energy efficiency of adding individual CCM components to a C 3 land plant is highly dependent on the permeability of lipid membranes to CO 2 , with values in the range reported in the literature that are higher than those used in previous modeling studies resulting in low light use efficiency. Adding a complete PCCM into the leaf cells of a C 3 land plant was predicted to boost net CO 2 fixation, but at higher energetic costs than those incurred by photorespiratory losses without a CCM. Two notable exceptions were when substomatal CO 2 levels are as low as those found in land plants that already use biochemical CCMs and when gas exchange is limited, such as with hornworts, making the use of a biophysical CCM necessary to achieve net positive CO 2 fixation under atmospheric CO 2 levels. This provides an explanation for the uniqueness of hornworts' CCM among land plants and the evolution of pyrenoids multiple times.

Kaste, Joshua A. M.↗
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

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

FROMP-from-Monomer

This repository provides the 1D reaction-diffusion simulation engine for the multi-scale framework. It takes monomer-level kinetic and thermodynamic descriptors as inputs, approximated from quantum chemical calculations, to simulate and predict the macroscopic propagation behavior of Frontal Ring-Opening Metathesis Polymerization (FROMP). By solving a coupled 1D reaction-diffusion PDE system, this tool predicts front propagation speed, front temperature, and degree of conversion. The solver explicitly captures the three fundamental ROMP steps, initiator activation/inhibition, chain initiation, and propagation, while accounting for high-temperature cycloreversion as a competing pathway. The numerical implementation uses FEniCS to handle the coupled heat diffusion and reaction kinetics.

Chua, Lauren [Massachusetts Inst. of Technology (M↗
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↗
DOE OSTI2024

Mechanochemical topological defects in an active nematic

We propose a reaction-diffusion system that converts topological information of an active nematic into chemical signals. We show that a curvature-activated reaction dipole is sufficient for creating a system that dynamically senses topology by producing a concentration field possessing local extrema coinciding with ±$\frac{1}{2}$ defects. The enabling term is analogous to polarization charge density seen in dielectric materials. We demonstrate the ability of this system to identify defects in both passive and active nematics. Our results illustrate that a relatively simple feedback scheme, expressed as a system of partial differential equations, is capable of producing chemical signals in response to inherently nonlocal structures in anisotropic media. Here, we posit that such coarse-grained systems can help generate testable hypotheses for regulated processes in biological systems, such as morphogenesis, and motivate the creation of bio-inspired materials that utilize dynamic coupling between nematic structure and biochemistry.

42 ENGINEERING↗
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

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 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↗
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

Tracking Dendritic Growth in Hydrogen-Based Hematite Reduction via Computer Vision

The reduction of hematite to metallic iron using hydrogen (H2) as a reducing agent presents a promising pathway for decarbonizing steel production. In this study, we employ a combination of in situ confocal scanning laser microscopy (CSLM) and advanced computer vision techniques to quantitatively analyze dendritic growth of ferrite during H2-based reduction of iron oxide at high temperatures. A workflow integrating Watershed Image Segmentation (WIS) and Lucas-Kanade Optical Flow (LKOF) is developed to extract both global and local kinetic information from time-resolved micrograph sequences. H2 reduction experiments conducted at 1400 degrees C and 1500 degrees C demonstrate a clear correlation between temperature and reduction rate, as evidenced by accuracy of fitted Johnson-Mehl-Avrami-Kolmogorov (JMAK) parameters. Optical flow analysis further elucidates the anisotropic and branched nature of dendritic growth, providing spatially resolved velocity fields that correlate well with global transformation kinetics. The proposed methodology demonstrates strong agreement with experimental measurements and literature values, offering a robust framework for automated image-based analysis to study kinetics through microstructural evolution in the reduction of iron ore, and likely other reaction-diffusion phenomena.

08 HYDROGEN↗
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 OSTI2025

Initial Stage of Nanoscale Imaging in Positive Tone Extreme UV Photoresists: The Influence of the Polymer Sequence

Photolithographic patterning using extreme ultraviolet (EUV, 92.5 eV) light is a radiolytic process that initially forms electrons, radical cations, anions, and neutral radicals in the polymeric photoresist matrix. These species may participate in the chemical reactions that define the ultimate resolution of the printed image, and their concentrations and nanometer-scale stochastic variations in their formation influence printed image quality. Proposals have been made that polymer chain uniformity may be advantageous in reducing stochastics due to spatial inhomogeneities, and this aspect of radiolysis is examined in this work. We have simulated the initial subpicosecond stages of the imaging process for a series of photoresist films that are identical in composition but vary in their polymer chain structures. We use detailed, physically accurate stochastic reaction-diffusion calculations to evaluate the influence of defined sequence and random copolymer structures on radiolytic spur formation, i.e., a cluster of species formed by electron-polymer interactions that defines the initial spatial characteristic of the imaging process. Predictions of electron thermalization in the present work are shown to be consistent with the literature, indicating that our overall computational approach for ultrafast nanoscale processes is sound. The computational results show that the polymer sequence has no significant effect on the spur composition. This suggests that any potential imaging improvements to be gained by sequence control must originate from postimaging lithographic process steps.

Absorption↗
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 OSTI2024

Examining Graphite Degradation in Molten Salt Environments: A Chemical, Physical, and Material Analysis

Molten-salt reactors (MSRs) are Generation IV nuclear reactors that use liquid salt as a coolant and/or fuel. In several MSR designs, graphite serves as a moderator and/or reflector. However, due to limited experimental data and operational experience, our understanding of graphite behavior in molten salt environments remains incomplete. This report aims to identify the degradation mechanisms of nuclear graphite in MSRs, detail the mechanisms of each factor, and provide an initial assessment of their impact on the structural integrity of graphite components. This assessment is based on an extensive literature review and insights from subject matter experts. Furthermore, given the limited data, a modeling strategy using existing Grizzly software is proposed for a more thorough analysis where appropriate. Additionally, it presents mitigation strategies where applicable. The report covers physical degradation mechanisms such as infiltration, erosion, and abrasion, as well as chemical degradation mechanisms including fluorination, intercalation, corrosion, and oxidation. Molten salt can infiltrate the porous structure of graphite, leading to several detrimental effects. Entrapment of fissile products within the graphite pores can cause radiation damage and could pose challenges in the handling and disposal of contaminated components. The differential thermal expansion between the infiltrated salt and graphite, along with internal stress from pressurized molten salt and volumetric heating, can compromise the structural integrity of graphite. To mitigate these effects, employing ultra-fine graphite grades and applying sealants and coatings are effective strategies. A computational model based on coupled solid mechanics and heat transfer phenomena could be used to predict the internal stresses using Grizzly software. In pebble-bed MSRs, graphite fuel pebbles can cause abrasion against reactor components due to friction and wear. The severity of wear is influenced by various factors such as temperature, environment, and the presence of lubricants. Tribological studies reveal that higher temperatures and molten salt environments, such as FLiBe, significantly reduce wear rates compared to dry conditions. Additionally, the chemical composition of the salt can further optimize graphite's tribological performance. Long-term wear effects can be modeled by incorporating surface defects into the geometry and predict stresses under thermal and radiation effects using Grizzly software. Chemical degradation of graphite in a molten salt environment can occur through fluorination and intercalation. Fluorination can occur via replacement of hydrogen or oxygen atoms, or at the active sites, but does not cause structural degradation. Intercalation, on the other hand, can lead to exfoliation, where layers of graphite separate and peel away, damaging the graphite. Protective coatings can enhance graphite's resistance to intercalation. Graphite generally exhibits good chemical stability in molten salt environments, though it can corrode under specific conditions, particularly in the presence of impurities or oxidants. Studies have shown that protective coatings, such as plasma-sprayed partially stabilized zirconia (PSZ), can effectively prevent such degradation. Corrosion behavior varies significantly with different graphite grades and coating applications, underscoring the need for detailed studies on uncoated and coated graphite to understand and mitigate corrosion mechanisms in MSRs. Research indicates that the presence of oxidants and impurities can accelerate graphite degradation in molten salts, making it essential to explore acceptable impurity limits. Oxidation is another critical degradation mechanism, leading to weight loss and structural damage due to the formation of CO and CO 2 from the reaction of carbon atoms with oxygen. This process creates new porosity and compromises graphite's integrity. While extensive research on graphite oxidation has been conducted for gas-cooled reactors, studies specific to MSRs are limited. Findings from the coal industry suggest that molten alkali metal salts can significantly accelerate graphite oxidation, a hypothesis worth exploring for fluoride salts in MSRs. Understanding oxidation behavior in MSRs is vital for developing protective measures. The analysis of post-irradiated graphite from the MSRE experiment demonstrated exceptional chemical compatibility with molten fluoride salt, suggesting that the extent of chemical attack on graphite largely depends on the salt's infiltration capability. Therefore, the use of ultra-fine grade graphite could help mitigate chemical degradation effects. Existing oxidation modeling capabilities in Grizzly, which use reaction-diffusion equations to model graphite-air interactions, could be adapted to simulate the chemical degradation effects of graphite in molten salt environments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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↗
DOE OSTI2025

Instability of steady-state mixed-state symmetry-protected topological order to strong-to-weak spontaneous symmetry breaking

Recent experimental progress in controlling open quantum systems enables the pursuit of mixed-state nonequilibrium quantum phases. We investigate whether open quantum systems hosting mixed-state symmetry-protected topological states as steady states retain this property under symmetric perturbations. Focusing on the decohered cluster state – a mixed-state symmetry-protected topological state protected by a combined strong and weak symmetry – we construct a parent Lindbladian that hosts it as a steady state. This Lindbladian can be mapped onto exactly solvable reaction-diffusion dynamics, even in the presence of certain perturbations, allowing us to solve the parent Lindbladian in detail and reveal previously-unknown steady states. Using both analytical and numerical methods, we find that typical symmetric perturbations cause strong-to-weak spontaneous symmetry breaking at arbitrarily small perturbations, destabilize the steady-state mixed-state symmetry-protected topological order. However, when perturbations introduce only weak symmetry defects, the steady-state mixed-state symmetry-protected topological order remains stable. Additionally, we construct a quantum channel which replicates the essential physics of the Lindbladian and can be efficiently simulated using only Clifford gates, Pauli measurements, and feedback.

Shah, Jeet [University of Maryland, College Park, ↗
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↗