Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion−reaction modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

43 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Numerical investigation of microbial quorum sensing under various flow conditions

Microorganisms efficiently coordinate phenotype expressions through a decision-making process known as quorum sensing (QS). We investigated QS amongst distinct, spatially distributed microbial aggregates under various flow conditions using a process-driven numerical model. Model simulations assess the conditions suitable for QS induction and quantify the importance of advective transport of signaling molecules. In addition, advection dilutes signaling molecules so that faster flow conditions require higher microbial densities, faster signal production rates, or higher sensitivities to signaling molecules to induce QS. However, autoinduction of signal production can substantially increase the transport distance of signaling molecules in both upstream and downstream directions. We present empirical approximations to the solutions of the advection–diffusion–reaction equation that describe the concentration profiles of signaling molecules for a wide range of flow and reaction rates. These empirical relationships, which predict the distribution of dissolved solutes along pore channels, allow to quantitatively estimate the effective communication distances amongst multiple microbial aggregates without further numerical simulations.

59 BASIC BIOLOGICAL SCIENCES↗
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 OSTI2020

A mass-transfer particle-tracking method for simulating transport with discontinuous diffusion coefficients

The problem of a spatially discontinuous diffusion coefficient (D(x)) is one that may be encountered in hydrogeologic systems due to natural geological features or as a consequence of numerical discretization of flow properties. To date, mass-transfer particle-tracking (MTPT) methods, a family of Lagrangian methods in which diffusion is jointly simulated by random walk and diffusive mass transfers, have been unable to solve this problem. This manuscript presents a new mass-transfer (MT) algorithm that enables MTPT methods to accurately solve the problem of discontinuous D(x). To achieve this, we derive a semi-analytical solution to the discontinuous D(x) problem by employing a predictor-corrector approach, and we use this semi-analytical solution as the weighting function in a reformulated MT algorithm. Furthermore, this semi-analytical solution is generalized for cases with multiple 1D interfaces as well as for 2D cases, including a 2 × 2 tiling of 4 subdomains that corresponds to a numerically-generated diffusion field. The solutions generated by this new mass-transfer algorithm closely agree with an analytical 1D solution or, in more complicated cases, trusted numerical results, demonstrating the success of our proposed approach.

42 ENGINEERING↗
DOE OSTI2020

Persistent compositions of non-stoichiometric compounds with low bulk diffusivity: A theory and application to Nb3Sn superconductors

Non-stoichiometric compounds may develop a composition gradient when they are formed by reactive diffusional processes. This paper reports an interesting phenomenon that in compounds with low bulk diffusivities, which rely mainly on grain boundary diffusion for their growth, the final bulk compositions may be far from equilibrium, with a very low bulk diffusivity leading to fixed (persistent) compositions – one such example is Nb 3 Sn, a superconductor. We investigated the microchemistry at the reactive interface using atom probe tomography to clarify the diffusion reaction mechanism for this low-bulk-diffusivity case and thus propose a theory for what determines the compound composition, using Nb 3 Sn as an example for concreteness. Using certain approximations, we derive an explicit analytical equation that illustrates what factors determine its composition profile. We compare our model with the known facts of Nb 3 Sn and see good agreement. In particular, this model predicts that internal oxidation may lead to higher Sn contents than conventional, non-oxidized Nb 3 Sn. Our measurements show that this is indeed true, and that the internally-oxidized Nb 3 Sn also has higher upper critical fields, achieving up to 28.2 T at 4.2 K. We discuss the general applicability of this model to non-stoichiometric compounds with low bulk diffusivity, and propose it as a tool to help in the design and processing of such materials for compositional control.

36 MATERIALS SCIENCE↗
DOE OSTI2021

DRACO: An Overview [Slides]

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

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

Direct simulations of chemically reacting turbulent mixing layers

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

Riley, J. J.↗
NASA NTRSDate not supplied

Extinction of the Stagnation Point Diffusion Flame: Effect of Conductive Heat Loss into Solid Interior

The planned microgravity experiment Growth and Extinction Limit (GEL) is designed to study the material flammability and extinction of thermally thick solid with an emphasis to quantify the effect of conductive heat loss into the solid interior by measuring the in-depth temperature profile. To obtain a more detailed understanding of the extinction phenomena, a previously developed one-dimensional axisymmetric stagnation-point diffusion flame model is expanded to account for the excessive heat loss into the solid beyond the steady burning state. This is important because for thick solids, this is important because the burning period can have a long and slow solid transient with a quasi-steady gas -phase. A non-dimensional excessive conductive heat loss parameter Ψ is defined in the paper and extensive numerical computations have been performed to determine the extinction boundary and the characteristics of near-limit flames. The classical U-shaped flammability boundary in the oxygen vs. stretch rate plot now becomes a series of U-shaped boundaries as a function of Ψ. For a given oxygen ambient, extinction boundary is a function of Ψ. The boundary consists of a blowoff branch and a radiative quenching branch. This is similar to the steady burning flammability boundary (Ψ =0) in the oxygen vs stretch rate plot. Extensive computed data including flame and pyrolysis temperatures, burning rate, flame standoff distance, and reaction rate are presented along the extinction boundaries. In particular, we note that the non-dimensional gas reaction rate is the ratio of diffusion reaction time to the reaction diffusion time, so it has the classical Damhkoler number character. By examining the maximum non-dimensional reaction rates, we observe that along the extinction boundaries, they are approximate constant across the entire range of stretch rates (from quenching to blow-off). These critical values are of the order unity but may vary in different oxygen ambient. The implication will be discussed in the presentation together with the near-limit flame structures.

Chengyao Li↗