Search NASA⌕ Search

SEARCH · Search NASA

Results for “reaction-diffusion system”

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

Accelerating network layouts using graph neural networks

Graph layout algorithms used in network visualization represent the first and the most widely used tool to unveil the inner structure and the behavior of complex networks. Current network visualization software relies on the force-directed layout (FDL) algorithm, whose high computational complexity makes the visualization of large real networks computationally prohibitive and traps large graphs into high energy configurations, resulting in hard-to-interpret “hairball” layouts. Here we use Graph Neural Networks (GNN) to accelerate FDL, showing that deep learning can address both limitations of FDL: it offers a 10 to 100 fold improvement in speed while also yielding layouts which are more informative. We analytically derive the speedup offered by GNN, relating it to the number of outliers in the eigenspectrum of the adjacency matrix, predicting that GNNs are particularly effective for networks with communities and local regularities. Finally, we use GNN to generate a three-dimensional layout of the Internet, and introduce additional measures to assess the layout quality and its interpretability, exploring the algorithm’s ability to separate communities and the link-length distribution. The novel use of deep neural networks can help accelerate other network-based optimization problems as well, with applications from reaction-diffusion systems to epidemics.

97 MATHEMATICS AND COMPUTING↗
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

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 OSTI2023

Probing the Structure and Dynamics of Fluid Mixtures in Porous Materials Through Ultrafast Vibrational Spectro-Microscopy and Many-Body Molecular Dynamics

This research program focused on the characterization of both molecular structure and dynamics of aqueous solutions in various metal-organic frameworks (MOFs) and organic polymers that have recently been proposed for applications in water treatment technologies. This was accomplished by integrating ultrafast vibrational spectroscopy and microscopy, which are sensitive to local environments, molecular orientation, and dynamical couplings, with many-body molecular dynamics (MB-MD) simulations, which enable realistic modeling of aqueous systems and porous materials under different thermodynamic conditions. Our studies allowed for identifying the physical mechanisms and characterizing the underlying molecular interactions between water and host materials that determine the adsorption and transport processes of water in various MOFs and β-cyclodextrin polymers – prototypical examples of porous materials and polymers, that have recently been proposed for technological applications in water purification. By combining bulk- and surface-sensitive, spatially resolved ultrafast vibrational spectroscopy with MB-MD simulations, it was possible to gain broad insight into the structure and mobility of water under heterogeneous confinement at various length scales. In particular, measurements of surface- and bulk-sensitive vibrational spectra at single crystal level, integrated with MB-MD simulations were used to characterize domain-specific, structure–binding affinity relationships. Besides providing specific information about the behavior of aqueous solutions in prototypical MOFs and organic polymers for potential applications in water treatment technologies, our studies also provided fundamental insights into the properties of aqueous solutions in confined heterogeneous environments, which have broad implications in many areas, ranging from heterogeneous catalysis to oil recovery, ion transport processes, and reaction-diffusion processes in crowded systems.

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

Long-time integration of parametric evolution equations with physics-informed DeepONets

Ordinary and partial differential equations (ODEs/PDEs) play a paramount role in analyzing and simulating complex dynamic processes across all corners of science and engineering. In recent years machine learning tools are aspiring to introduce new effective ways of simulating such equations, however existing approaches are not able to reliably return stable and accurate predictions across long temporal horizons. We aim to address this challenge by introducing an effective framework for learning evolution operators that map random initial conditions to associated ODE/PDE solutions within a short time interval. Such operators can be parametrized by deep neural networks that are trained in an entirely self-supervised manner without requiring one to generate any paired input-output observations. Global long-time predictions across a range of initial conditions can be then obtained by iteratively evaluating the trained model using each prediction as the initial condition for the next evaluation step. Here, this introduces a new approach to temporal domain decomposition that is shown to be effective in performing accurate long-time simulations for a wide range of parametric ODE and PDE systems, from wave propagation, to reaction-diffusion dynamics and stiff chemical kinetics, introducing a new way of rapidly emulating non-equilibrium processes in science and engineering.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Multi-resolution partial differential equations preserved learning framework for spatiotemporal dynamics

Traditional data-driven deep learning models often struggle with high training costs, error accumulation, and poor generalizability in complex physical processes. Physics-informed deep learning (PiDL) addresses these challenges by incorporating physical principles into the model. Most PiDL approaches regularize training by embedding governing equations into the loss function, yet this depends heavily on extensive hyperparameter tuning to weigh each loss term. To this end, we propose to leverage physics prior knowledge by “baking” the discretized governing equations into the neural network architecture via the connection between the partial differential equations (PDE) operators and network structures, resulting in a PDE-preserved neural network (PPNN). This method, embedding discretized PDEs through convolutional residual networks in a multi-resolution setting, largely improves the generalizability and long-term prediction accuracy, outperforming conventional black-box models. The effectiveness and merit of the proposed methods have been demonstrated across various spatiotemporal dynamical systems governed by spatiotemporal PDEs, including reaction-diffusion, Burgers’, and Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2023

Reaction-diffusive dynamics of number-conserving dissipative quantum state preparation

The use of dissipation for the controlled creation of nontrivial quantum many-body correlated states is of much fundamental and practical interest. What is the result of imposing number conservation, which, in closed system, gives rise to diffusive spreading? For this work, we investigate this question for a paradigmatic model of a two-band system, with dissipative dynamics aiming to empty one band and to populate the other, which had been introduced before for the dissipative stabilization of topological states. Going beyond the mean-field treatment of the dissipative dynamics, we demonstrate the emergence of a diffusive regime for the particle and hole density modes at intermediate length- and timescales, which, interestingly, can only be excited in nonlinear response to external fields. We also identify processes that limit the diffusive behavior of this mode at the longest length- and timescales. Strikingly, we find that these processes lead to a reaction-diffusion dynamics governed by the Fisher-Kolmogorov-Petrovsky-Piskunov equation, making the designed dark state unstable towards a state with a finite particle and hole density.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2023

A new approach for simulating inhomogeneous chemical kinetics

Abstract In this paper, inhomogeneous chemical kinetics are simulated by describing the concentrations of interacting chemical species by a linear expansion of basis functions in such a manner that the coupled reaction and diffusion processes are propagated through time efficiently by tailor-made numerical methods. The approach is illustrated through modelling $$\alpha$$ α - and $$\gamma$$ γ -radiolysis in thin layers of water and at their solid interfaces from the start of the chemical phase until equilibrium was established. The method’s efficiency is such that hundreds of such systems can be modelled in a few hours using a single core of a typical laptop, allowing the investigation of the effects of the underlying parameter space. Illustrative calculations showing the effects of changing dose-rate and water-layer thickness are presented. Other simulations are presented which show the approach’s capability to solve problems with spherical symmetry (an approximation to an isolated radiolytic spur), where the hollowing out of an initial Gaussian distribution is observed, in line with previous calculations. These illustrative simulations show the generality and the computational efficiency of this approach to solving reaction-diffusion problems. Furthermore, these example simulations illustrate the method’s suitability for simulating solid-fluid interfaces, which have received a lot of experimental attention in contrast to the lack of computational studies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Reversible Reactions, Mesh Size, and Segmental Dynamics Control Penetrant Diffusion in Ethylene Vitrimers

The diffusion of two aromatic dyes with nearly identical sizes was measured in ethylene vitrimers with precise linker lengths and borate ester cross-links using fluorescence recovery after photobleaching (FRAP). One dye possessed a reactive hydroxyl group, while the second was inert. The reaction of the hydroxyl group with the network is slow relative to the hopping times of the dye, resulting in a large slowdown by a factor of 50 for a reactive probe molecule. A kinetic model was fit to the fluorescence intensity data to determine rate constants for the reversible reaction of the dye from the network, which confirms the role of slow reaction kinetics. A second network cross-linker was also investigated with a substituted boronic ester showing ∼10,000 times faster exchange kinetics. In this system, the two dyes show the same diffusion coefficient, as the reaction is no longer the rate-limiting step. The role of dense meshes on small and large dyes is also discussed in the context of the existing theories. Finally, these results highlight the potential of dynamic networks to control penetrant transport through synergistic effects of the mesh size, dynamic bond kinetics, and penetrant–network interactions.

confinement↗