Search NASA⌕ Search

SEARCH · Search NASA

Results for “reaction-diffusion equation”

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

A deterministic particle method for one-dimensional reaction-diffusion equations

We derive a deterministic particle method for the solution of nonlinear reaction-diffusion equations in one spatial dimension. This deterministic method is an analog of a Monte Carlo method for the solution of these problems that has been previously investigated by the author. The deterministic method leads to the consideration of a system of ordinary differential equations for the positions of suitably defined particles. We then consider the time explicit and implicit methods for this system of ordinary differential equations and we study a Picard and Newton iteration for the solution of the implicit system. Next we solve numerically this system and study the discretization error both analytically and numerically. Numerical computation shows that this deterministic method is automatically adaptive to large gradients in the solution.

Mascagni, Michael↗

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

Periodic solutions to systems of reaction-diffusion equations

Necessary and sufficient conditions are derived for the existence of temporally periodic 'dissipative structure' solutions in weak diffusion with the reaction rate terms dominant in a generic system of reaction-diffusion differential equations. The enumerator index i of the equations denotes the density or concentration of the ith participating molecular or biological species, and D sub i is the diffusivity constant for the ith species while Q sub i (c), an algebraic function of the n-tuple c, expresses the local rate of production of the ith species due to chemical reactions or biological interactions.

Rosen, G.↗

Numerical Determination of Critical Conditions for Thermal Ignition

The determination of ignition or thermal explosion in an oxidizing porous body of material, as described by a dimensionless reaction-diffusion equation of the form .tu = .2u + .e-1/u over the bounded region O, is critically reexamined from a modern perspective using numerical methodologies. First, the classic stationary model is revisited to establish the proper reference frame for the steady-state solution space, and it is demonstrated how the resulting nonlinear two-point boundary value problem can be reexpressed as an initial value problem for a system of first-order differential equations, which may be readily solved using standard algorithms. Then, the numerical procedure is implemented and thoroughly validated against previous computational results based on sophisticated path-following techniques. Next, the transient nonstationary model is attacked, and the full nonlinear form of the reaction-diffusion equation, including a generalized convective boundary condition, is discretized and expressed as a system of linear algebraic equations. The numerical methodology is implemented as a computer algorithm, and validation computations are carried out as a prelude to a broad-ranging evaluation of the assembly problem and identification of the watershed critical initial temperature conditions for thermal ignition. This numerical methodology is then used as the basis for studying the relationship between the shape of the critical initial temperature distribution and the corresponding spatial moments of its energy content integral and an attempt to forge a fundamental conjecture governing this relation. Finally, the effects of dynamic boundary conditions on the classic storage problem are investigated and the groundwork is laid for the development of an approximate solution methodology based on adaptation of the standard stationary model.

Luo, W.↗

Finite flame thickness effects on Kolmogorov-Petrovsky-Piskunov turbulent burning velocities

KPP (Kolmogorov-Petrovsky-Piskunov) solutions of the reaction-diffusion equation have application in various physical phenomena occurring in biology, ecology, and reacting flows. In particular, these solutions are commonly used in turbulent combustion to scale turbulent burning velocities. Subject to certain conditions on reaction rate profile through the flame brush and turbulent diffusivity, this theory relates the turbulent burning velocity to the derivative of the reaction rate ($\tilde{ω}$) at the leading edge of the flame brush ($d\tilde{ω}$/$d\tilde{c}$ | $\tilde{c} =$0 ). Such waves are often referred to as “pulled fronts.” However, turbulent flames never actually satisfy the KPP conditions for a pulled front, as the turbulent flame brush, parametrized here by the thickness δ t , consists of an ensemble of laminar flamelets of thickness δ, where ε=δ/δ t $\ll$ 1 is very small, but nonzero, and $d\tilde{ω}$ /$d\tilde{c}$ tends to zero at the brush leading edge for high activation energy, combustion-type kinetics. Here, this paper analyzes these effects on KPP wave solutions, parametrized by ε=δ/δ t and Zeldovich number Ze focusing on whether turbulent flames retain their pulled front character and what the correction to the KPP wave speed is. Variational solutions of the reaction-diffusion equation show that the solution can be expanded in powers of 1/|ln⁡ε|. Both numerical and asymptotic results are presented, showing that the wave still exhibits pulled front solutions but with significant corrections to the KPP result. The leading order correction is of the form |ln⁡ε| –2 and independent of Ze. Higher order corrections are function of both ε and Ze. However, the dominant factor influencing the wave speed correction is due to the finite ε, with Ze exhibiting a weaker effect.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Asymptotic analysis of dissipative waves with applications to their numerical simulation

Various problems involving the interplay of asymptotics and numerics in the analysis of wave propagation in dissipative systems are studied. A general approach to the asymptotic analysis of linear, dissipative waves is developed. It was applied to the derivation of asymptotic boundary conditions for numerical solutions on unbounded domains. Applications include the Navier-Stokes equations. Multidimensional traveling wave solutions to reaction-diffusion equations are also considered. A preliminary numerical investigation of a thermo-diffusive model of flame propagation in a channel with heat loss at the walls is presented.

Hagstrom, Thomas↗

Multi-dimensional initial-boundary value problems with strong nonlinearities

Initial boundary value problems for the general scalar singularly perturbed parabolic equation in a cylindrical domain are studied. The one-dimensional reaction-diffusion equation is considered in the absence of convection. Then, the general one-dimensional and N-dimensional reaction-convection-diffusion equations are analyzed.

Howes, F. A.↗

Fast estimation of reaction rates in spherical and non-spherical porous catalysts

Here we present a methodology for modeling multi-step reaction rates in porous catalyst particles for use in CFD-DEM and two fluid models. Single-step effectiveness factors based on a Thiele modulus, while useful, cannot accurately capture the cascading reaction systems common in high temperature vapor-phase chemical reactors like fluidized catalytic cracking units and catalytic biomass fast pyrolysis systems. Instead, multi-step effectiveness vectors derived from steady-state solutions to the governing reaction-diffusion equations are needed. Solutions for various catalyst shapes are presented, including spheres, cylinders, and prisms. Computational challenges inherent in repeated evaluation of reaction rates with diffusion limitations are discussed, and an efficient implementation based on pre-computed lookup tables is proposed and demonstrated on a simulation of a fluidized bed reactor. Open-source code is provided for the compilation of reaction rate tables for use in ODE, DEM, and two-fluid models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

The thermal stability of coronal loops by nonlinear diffusion asymptotics

A nonlinear reaction-diffusion equation and some additional constraints are derived which describe the time-dependent behavior of the temperature structure of the plasma in coronal loops. The equation is analyzed using nonlinear diffusion asymptotics, in particular singular perturbation techniques, and the results are interpreted in the context of the physical problem of the thermal stability and temporal behavior of the plasma. The results are consistent with the possibility of cyclic thermal behavior of the plasma, as suggested by Kuin and Martens (1982).

Pakkert, J. W.↗

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↗

Discovering a reaction–diffusion model for Alzheimer’s disease by combining PINNs with symbolic regression

Misfolded tau proteins play a critical role in the progression and pathology of Alzheimer's disease. Recent studies suggest that the spatio-temporal pattern of misfolded tau follows a reaction-diffusion type equation. However, the precise mathematical model and parameters that characterize the progression of misfolded protein across the brain remain incompletely understood. Here, we use deep learning and artificial intelligence to discover a mathematical model for the progression of Alzheimer's disease using longitudinal tau positron emission tomography from the Alzheimer's Disease Neuroimaging Initiative database. Specifically, we integrate physics informed neural networks (PINNs) and symbolic regression to discover a reaction-diffusion type partial differential equation for tau protein misfolding and spreading. First, we demonstrate the potential of our model and parameter discovery on synthetic data. Then, we apply our method to discover the best model and parameters to explain tau imaging data from 46 individuals who are likely to develop Alzheimer's disease and 30 healthy controls. Our symbolic regression discovers different misfolding models f(c) for two groups, with a faster misfolding for the Alzheimer's group, f(c) = 0.23c 3 – 1.34c 2 + 1.11c, than for the healthy control group, f(c) = –c 3 + 0.62c 2 + 0.39c. Our results suggest that PINNs, supplemented by symbolic regression, can discover a reaction-diffusion type model to explain misfolded tau protein concentrations in Alzheimer's disease. Furthermore, we expect our study to be the starting point for a more holistic analysis to provide image-based technologies for early diagnosis, and ideally early treatment of neurodegeneration in Alzheimer's disease and possibly other misfolding-protein based neurodegenerative disorders.

60 APPLIED LIFE SCIENCES↗

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

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↗

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↗

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↗

Numerical studies of the thermal design sensitivity calculation for a reaction-diffusion system with discontinuous derivatives

The aim of this study is to find a reliable numerical algorithm to calculate thermal design sensitivities of a transient problem with discontinuous derivatives. The thermal system of interest is a transient heat conduction problem related to the curing process of a composite laminate. A logical function which can smoothly approximate the discontinuity is introduced to modify the system equation. Two commonly used methods, the adjoint variable method and the direct differentiation method, are then applied to find the design derivatives of the modified system. The comparisons of numerical results obtained by these two methods demonstrate that the direct differentiation method is a better choice to be used in calculating thermal design sensitivity.

Hou, Jean W.↗

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↗