A note on the standard diffusion curve of TAP analysis

Toby Isaac

DOI: https://doi.org/10.48550/arXiv.2302.03772

2023-02-08

Abstract:The standard diffusion curve used in models of TAP reactors, as it is usually defined, is numerically unstable for small values. We use a functional equation satisfied by the curve to define a numerically stable way of computing it for all values.

Mathematical Software

What problem does this paper attempt to address?

An Alternate Mathematical Model for the Tap Technique for Numerical Integration

DZ Wang

DOI: https://doi.org/10.1002/jccs.200300081

2003-01-01

Journal of the Chinese Chemical Society

Abstract:The mathematical modeling of the reactor in the TAP technique describes it as composed of zones of different diffusivity, and it makes the model linear by imposing boundary conditions of the continuity of concentration and flux at zone boundaries so as to be able to use a Fourier transformation technique to solve the equations. An alternative description of the reactor that uses a spatially non-constant diffusivity to obviate the need to impose auxiliary boundary conditions in the regions where diffusivities change is proposed. The reactor equation is the mass balance equation without the simplification to the case of constant diffusivity, that is, without the diffusivity taken out of the differential operator. Computations based on it are compared with those of the linear model. The computed response curves of both models are similar, which indicates that the alternate model is equally good as a description of the reaction system. The alternate model is useful for modeling reactor systems with a variable cross-sectional area or with temperature non-uniformity.
Reaction Diffusion TAP

M. Andrecut

2023-09-07

Abstract:The recently introduced Theory of the Adjacent Possible (TAP) is a model of combinatorial innovation aiming to explain the "hockey-stick" upward trend of human technological evolution, where an explosion in the number of produced items with increasing complexity suddenly occurs. In addition, the TAP model was also used to explain the rapidly emerging biological complexity. Inspired by TAP here we propose a reaction-diffusion system aiming to extend the model in both space and time. We show that the new model exhibits similar characteristics to the TAP model, like the sudden increase in the production of items, after a longer period of slow growth. The new model also exhibits wave propagation of "innovation", resulting in self-sustained complex interference patterns.

Populations and Evolution,Data Analysis, Statistics and Probability
TUNED FINITE-DIFFERENCE DIFFUSION OPERATORS

Jason Maron,Mordecai-Mark Mac Low

DOI: https://doi.org/10.1088/0067-0049/182/1/468

2009-05-01

The Astrophysical Journal Supplement Series

Abstract:Finite-difference simulations of fluid dynamics and magnetohydrodynamics generally require an explicit diffusion operator, either to maintain stability by attenuating grid-scale structure, or to implement physical diffusivities such as viscosity or resistivity. If the goal is stability only, the diffusion must act at the grid scale, but should affect structure at larger scales as little as possible. For physical diffusivities the diffusion scale depends on the problem, and diffusion may act at larger scales as well. Diffusivity can undesirably limit the computational time step in both cases. We construct tuned finite-difference diffusion operators that minimally limit the time step while acting as desired near the diffusion scale. Such operators reach peak values at the diffusion scale rather than at the grid scale, but behave as standard operators at larger scales. These operators will be useful for simulations with high magnetic diffusivity or kinematic viscosity such as in the simulation of astrophysical dynamos with magnetic Prandtl number far from unity, or for numerical stabilization using hyperdiffusivity.

astronomy & astrophysics
Boundary conditions in TAP reactor modeling

Lei Gao,Xueliang Zhao,Dezheng Wang

DOI: https://doi.org/10.1016/S0167-2991(06)81687-3

2006-01-01

Abstract:The two boundary conditions (BC) in TAP reactor modeling: the use of the delta function to approximate the pulse input and the assumption of zero concentration (infinitely large pumping speed) at the reactor outlet, are discussed. When the delta function approximation is used, curve fitting should be limited to the latter 3/4 part of the response curve for curves of width < 3 times pulse width, while for curves with widths > 4 times pulse width, it is a fair approximation. The infinitely large pumping speed BC is not a good assumption even for a 1,500 ls(-1) pumping speed and broad outlet pulses with widths > 1000 ms.
Thermally Induced Diffusion of Chemicals under Steady-State Heat Transfer in Saturated Porous Media

Huaxiang Yan,Majid Sedighi,Haijian Xie

DOI: https://doi.org/10.1016/j.ijheatmasstransfer.2020.119664

IF: 5.2

2020-01-01

International Journal of Heat and Mass Transfer

Abstract:Thermally induced diffusion of chemicals is observed in many natural phenomena and engineering processes. We present analytical solutions to the non-isothermal diffusion of chemicals under combined effects of molecular and thermal diffusion processes and under variable diffusion coefficient with temperature. The generalized integral transform technique (GITT) is adopted to successfully derive the analytical solutions to the one-dimensional boundary problem under Neumann and Dirichlet outflow conditions. The accuracy of the analytical model developed is examined against three benchmarks including two experimental datasets. The analytical models are applied to develop an understanding of how the chemical transfer under non-isothermal conditions is governed by a combination of the Soret effect and temperature dependency of diffusion coefficient. We demonstrate that the mass transfer components due to thermal diffusion (or Soret effect) and the temperature dependency of diffusion coefficients are equally important to develop an accurate assessment of non-isothermal diffusion problem. Under the conditions of the case study considered, we show that a threshold for the Soret coefficient exists at which the temperature dependency of chemical diffusion coefficient can induce equally or larger influences on the chemical compared to the influence of thermal diffusion. The analytical solutions proposed also provide modelling benchmarks to examine the accuracy of alternative multiphysics numerical models for studying the non-isothermal diffusion problems in porous media. (C) 2020 Published by Elsevier Ltd.
The measurement of gas diffusivity in porous materials by temporal analysis of products (TAP)

Wenjun Li,Lanying Xie,Lei Gao,Xueliang Zhao,Rongrong Hu,Yi Cheng,Dezheng Wang

DOI: https://doi.org/10.1016/j.cattod.2006.07.020

IF: 5.3

2007-01-01

Catalysis Today

Abstract:The pulsing of argon in a temporal analysis of products (TAP) reactor and reactor modeling of the response curves were used to measure the effective intraparticle diffusivities in porous materials. The diffusivity that can be measured is limited: (1) at the low end by intraparticle diffusion being too slow such that just a small fraction of the pulse gets into the pores to give an indistinguishable tail, which only measures that the diffusivity is smaller than an upper limit and (2) at the high end by intraparticle diffusion being too fast such that it gives a constant concentration in the pores, which only measures that the diffusivity is larger than a lower limit. The limits and range are slightly different for different particle and bed dimensions. A 9mm long packed bed has a sensitive range of about 300-fold where there are discernible changes in the normalized pulse shape due to diffusivity changes. If small particles of about 50μm are used, the range is from 1×10−10 to 3×10−8m2/s, and if large particles of about 500μm are used, the range is from 2×10−9 to 5×10−7m2/s.
Numerical Simulation of Nanoparticle Synthesis in Diffusion Flame Reactor

Minfzhou Yu,Jianzhong Lin,Tatleung Chan

DOI: https://doi.org/10.1016/j.powtec.2007.03.037

IF: 5.64

2008-01-01

Powder Technology

Abstract:A computational model combining the fluid dynamics with the particle kinetics was employed to study TiO2 nanoparticle synthesis in a diffusion flame reactor. A one-step chemical kinetics approach was used to model titanium tetraisopropoxide (TTIP) decomposition that leads to homogeneous nucleation and particle formation. The momentum, heat, and mass transfer, Brownian coagulation and diffusion, surface growth, coalescence and thermophoresis have been taken into account. Based on the particle size distributions, an efficient quadrature method of moments was allowed to approximate the general dynamics equation of particle, and the eddy dissipation concept (EDC) combustion model was used to estimate the flame temperature field. Excellent agreements between the model predictions and experimental data, with respect to the flame temperature distribution and particle kinetics, are reached. By taking the particle size and surface area as independent variables, the full distributions of volume equivalent diameters, evolution of the agglomerate number, the geometric standard deviation based on volume and agglomerate fractal nature, mean primary particle size and the number of primary particles per agglomerate are revealed. The variation of oxygen flow rate has an important influence on the temperature distribution and hence on the particle kinetics accordingly.
Numerical study of anisotropic diffusion in Turing patterns based on Finsler geometry modeling

Gildas Diguet,Madoka Nakayama,Sohei Tasaki,Fumitake Kato,Hiroshi Koibuchi,Tetsuya Uchimoto

DOI: https://doi.org/10.1103/physreve.109.014213

IF: 2.707

2024-01-12

Physical Review E

Abstract:We numerically study the anisotropic Turing patterns (TPs) of an activator-inhibitor system described by the reaction-diffusion (RD) equation of Turing, focusing on anisotropic diffusion using the Finsler geometry (FG) modeling technique. In FG modeling, the diffusion coefficients are dynamically generated to be direction dependent owing to an internal degree of freedom (IDOF) and its interaction with the activator and inhibitor. Because of this dynamical diffusion coefficient, FG modeling of the RD equation sharply contrasts with the standard numerical technique in which direction-dependent coefficients are manually assumed. To find the solution of the RD equations in FG modeling, we use a hybrid numerical technique combining the Metropolis Monte Carlo method for IDOF updates and discrete RD equations for steady-state configurations of the activator-inhibitor variables. We find that the newly introduced IDOF and its interaction are a possible origin of spontaneously emergent anisotropic patterns of living organisms, such as zebra and fishes. Moreover, the IDOF makes TPs controllable by external conditions if the IDOF is identified with the direction of cell diffusion accompanied by thermal fluctuations. https://doi.org/10.1103/PhysRevE.109.014213 ©2024 American Physical Society

physics, fluids & plasmas, mathematical
Uniform Convergent Tailored Finite Point Method for Advection–Diffusion Equation with Discontinuous, Anisotropic and Vanishing Diffusivity

Min Tang,Yihong Wang

DOI: https://doi.org/10.1007/s10915-016-0254-1

2016-01-01

Journal of Scientific Computing

Abstract:We propose two tailored finite point methods for the advection–diffusion equation with anisotropic tensor diffusivity. The diffusion coefficient can be very small in one direction in some part of the domain and be discontinuous at the interfaces. When flows advect from the vanishing-diffusivity region towards the non-vanishing diffusivity region, standard numerical schemes tend to cause spurious oscillations or negative values. Our proposed schemes have uniform convergence in the vanishing diffusivity limit, even when the solution exhibits interface and boundary layers. When the diffusivity is along the coordinates, the positivity and maximum principle can be proved. We use the value as well as their derivatives at the grid points to construct the scheme for nonaligned case, which makes it can achieve good accuracy and convergence for the derivatives as well, even when there exhibit boundary or interface layers. Numerical experiments are presented to show the performance of the proposed scheme.
Optimal finite-differences discretization for the diffusion equation from the perspective of large-deviation theory

Naftali R. Smith

DOI: https://doi.org/10.1088/1742-5468/ad363f

2024-04-08

Abstract:When applying the finite-differences method to numerically solve the one-dimensional diffusion equation, one must choose discretization steps $\Delta x$, $\Delta t$ in space and time, respectively. By applying large-deviation theory on the discretized dynamics, we analyze the numerical errors due to the discretization, and find that the (relative) errors are especially large in regions of space where the concentration of particles is very small. We find that the choice $\Delta t = {\Delta x}^2 / (6D)$, where $D$ is the diffusion coefficient, gives optimal accuracy compared to any other choice (including, in particular, the limit $\Delta t \to 0$), thus reproducing the known result that may be obtained using truncation error analysis. In addition, we give quantitative estimates for the dynamical lengthscale that describes the size of the spatial region in which the numerical solution is accurate, and study its dependence on the discretization parameters. We then turn to study the advection-diffusion equation, and obtain explicit expressions for the optimal $\Delta t$ and other parameters of the finite-differences scheme, in terms of $\Delta x$, $D$ and the advection velocity. We apply these results to study large deviations of the area swept by a diffusing particle in one dimension, trapped by an external potential $\sim |x|$. We extend our analysis to higher dimensions by combining our results from the one dimensional case with the locally one-dimension method.

Statistical Mechanics,Computational Physics
Numerical Recovery of the Diffusion Coefficient in Diffusion Equations from Terminal Measurement

Bangti Jin,Xiliang Lu,Qimeng Quan,Zhi Zhou

2024-05-17

Abstract:In this work, we investigate a numerical procedure for recovering a space-dependent diffusion coefficient in a (sub)diffusion model from the given terminal data, and provide a rigorous numerical analysis of the procedure. By exploiting decay behavior of the observation in time, we establish a novel H{ö}lder type stability estimate for a large terminal time $T$. This is achieved by novel decay estimates of the (fractional) time derivative of the solution. To numerically recover the diffusion coefficient, we employ the standard output least-squares formulation with an $H^1(\Omega)$-seminorm penalty, and discretize the regularized problem by the Galerkin finite element method with continuous piecewise linear finite elements in space and backward Euler convolution quadrature in time. Further, we provide an error analysis of discrete approximations, and prove a convergence rate that matches the stability estimate. The derived $L^2(\Omega)$ error bound depends explicitly on the noise level, regularization parameter and discretization parameter(s), which gives a useful guideline of the \textsl{a priori} choice of discretization parameters with respect to the noise level in practical implementation. The error analysis is achieved using the conditional stability argument and discrete maximum-norm resolvent estimates. Several numerical experiments are also given to illustrate and complement the theoretical analysis.

Numerical Analysis
Identification of Concentration and Residence Time Distributions with One-Dimension Diffusion Model at Close Boundary Conditions

Li Li,Yang Guangyu,Wang Jinfu

DOI: https://doi.org/10.3969/j.issn.2095-2783.2006.01.008

2006-01-01

Abstract:The paper identified the concentration and residence time distributions by analytically solving Laplace transformation of a one-dimension diffusion model.The model Was built-up at the close boundary conditions with the pulse input of the tracer into reactor.The response curves of concentration at the outlet of reactor were presented at the different values of Peclet number.Furthermore,the changes of the concentration with time and location in the reactor were described for a fixed Peclet number.The comparisons of the response curves of the concentration were made for the open and close boundary conditions.Also,the characteristics of residence time distribution were discussed.
A High Accuracy Compact Difference Scheme for Convection Diffusion Reaction Equation on Non-uniform Grid

田芳,田振夫

DOI: https://doi.org/10.3969/j.issn.1005-3085.2009.02.004

2009-01-01

Abstract:A high accuracy compact finite difference scheme on the non-uniform grid is proposed to solve the convection-diffusion-reaction equation. Based on the Taylor series expansion, we first constructed the approximate expressions for the 1st-order and 2nd-order derivative; then, by rewriting the model equation in the form of the convection diffusion equation and utilizing the model equation, we get the finite difference scheme and illustrate by four numerical examples. The numerical results show that the presented scheme has many advantages such as yielding more accurate numerical solutions, having high resolution for the large gradient changes of the unknown quantity, being suitable for both convection-dominant flows and diffusion-dominant flows and so on.
A Diffusion Scale Approximation for Stochastic Fission Chains with State Feedback

Chen Dubi,Anil K. Prinja

DOI: https://doi.org/10.1080/00295639.2024.2324532

IF: 1.2

2024-04-06

Nuclear Science and Engineering

Abstract:Modeling and simulation of stochastic fission chains, often referred to as zero power reactor noise or stochastic transport, is a central topic in nuclear science and engineering, with important practical applications in reactor control and measurements. The common mathematical setting for studying reactor noise is through branching processes, typically modeled using forward and backward master equations. Because of the high complexity of the problem, reactor noise is typically studied under the point model approximation and in a linear setting, which neglects any reactivity feedbacks. On the other hand, since reactivity feedbacks are dominant in power reactors, stochastic models and simulations in a nonlinear setting are of great interest. For this reason, there is a constant interest in new mathematical frameworks that will enable modeling and simulation in a complete fashion. In the present study, we look at a branching process with a negative feedback, realized by a linear increase in the per capita death rate. In particular, the outline of the study is to compare two models for the above process. The first is the exact model, based on a forward master equation, and the second is the diffusion scale approximation (or the functional central limit approximation), realized in an Ito-type stochastic differential equation (SDE). The comparison between the two models is done in two levels. First, we show that the dynamics of the first and second moments of both models are governed by the same equations. Then, by realization of each model, we contract a Monte Carlo simulation for sampling the population size distribution. In particular, we show that once the population size is sufficiently large, the SDE approximation allows us to perform accurate simulations in a critical setting, with a dramatic reduction of the CPU time.

nuclear science & technology
Application of the Discontinuity Factor Theory to Accelerate Routine Rod Worth Calculations for Modular HTRs

Zhou Xuhua,Li Fu,Wang Dengying

DOI: https://doi.org/10.1016/j.nucengdes.2008.10.018

IF: 1.7

2008-01-01

Nuclear Engineering and Design

Abstract:A modular high temperature gas cooled reactor (HTR) is equipped with strong absorbers and many void areas in the control rod regions which are located in the graphitic side reflector. When using diffusion calculations, a special skill is required to homogenize these ex-core absorber regions. These difficulties arise from the fact that these absorber regions do not contain fission sources while the boundaries of the absorber regions are exposed to very strong neutron currents. This paper shows that a mere conventional homogenization technique is not enough to obtain acceptable results. It is shown that a discontinuity factor-corrected diffusion-solution can be found which yields very good, satisfactory results for fast, whole core diffusion calculations. The verification is done by comparing time-consuming transport calculations [SN-codes or multi-group MCNP-codes] with a fast running diffusion code, corrected by the newly proposed homogenization technique.
An iterative scheme to include turbulent diffusion in advective-dominated transport of delayed neutron precursors

Mathis Caprais,André Bergeron

2024-11-06

Abstract:In this study, the Method of Characteristics (MOC) for Delayed Neutron Precursors (DNPs) is used to solve the precursors balance equation with turbulent diffusion. The diffusivity of DNPs, significantly higher than molecular diffusivity, emerges in turbulent flows from the time-averaging of the DNPs mass balance equation. To integrate this effect within the MOC framework, the advection-reaction component of the DNPs balance equation is solved using the MOC, while the diffusive source is computed from the concentration of the previous iteration. The method is validated on a 2D recirculating pipe reactor with high Reynolds number flow, comparing the MOC with diffusion to a standard finite volume (FV) discretization of the fission products balance equation. Additionally, the impact of the diffusivity term on DNP distributions and reactor reactivity is quantified as a function of the turbulent Schmidt number.

Computational Physics,Fluid Dynamics
A convergent discretisation method for transition path theory for diffusion processes

Nada Cvetković,Tim Conrad,Han Cheng Lie

DOI: https://doi.org/10.1137/20M1329354

2020-04-02

Abstract:Transition path theory (TPT) for diffusion processes is a framework for analysing the transitions of multiscale ergodic diffusion processes between disjoint metastable subsets of state space. Most methods for applying TPT involve the construction of a Markov state model on a discretisation of state space that approximates the underlying diffusion process. However, the assumption of Markovianity is difficult to verify in practice, and there are to date no known error bounds or convergence results for these methods. We propose a Monte Carlo method for approximating the forward committor, probability current, and streamlines from TPT for diffusion processes. Our method uses only sample trajectory data and partitions of state space based on Voronoi tessellations. It does not require the construction of a Markovian approximating process. We rigorously prove error bounds for the approximate TPT objects and use these bounds to show convergence to their exact counterparts in the limit of arbitrarily fine discretisation. We illustrate some features of our method by application to a process that solves the Smoluchowski equation on a triple-well potential.

Numerical Analysis
Beyond the effectiveness factor: Multi-step reactions with intraparticle diffusion limitations

Aaron M. Lattanzi,M. Brennan Pecha,Vivek S. Bharadwaj,Peter N. Ciesielski

DOI: https://doi.org/10.1016/j.cej.2019.122507

IF: 15.1

2020-01-01

Chemical Engineering Journal

Abstract:Though common in many industrial chemical processes, multi-step reactions occurring within porous catalyst particles provide significant challenges to previously established modeling frameworks. In this work, a system of reaction-diffusion equations is solved analytically via eigenanalysis. The intraparticle concentration profiles are integrated over the particle volume to obtain effectiveness factors for each reacting species, referred to herein as a multi-step effectiveness vector (MEV). The MEV solution is tested for the case of uniform flow past a stationary catalyst particle with a state-of-the-art, multi-step reaction scheme for fluid catalytic cracking (FCC). Outputs from particle-resolved direct numerical simulation (DNS) are compared to outputs from ordinary differentials with effective reaction closures coming from the MEV solution, the classic single-step effectiveness factor (SEF) (R. Aris, "The mathematical theory of diffusion and reaction in permeable catalysts: The theory of the steady state," Oxford University Press, 1975), and utilizing the intrinsic rates as-is. When compared to DNS, the MEV solution shows excellent agreement for the formation of all species and quantitatively captures the intraparticle concentration profiles. By contrast, the intrinsic treatment shows an over-prediction for product formation and the SEF treatment shows mixed levels of accuracy that depend upon the reaction coupling. Specifically, the SEF treatment accurately characterizes product formation if the kinetic scheme is weakly coupled but shows significant departure if stronger coupling is present. The computational cost associated with the MEV solution is associated with eigendecomposition of the Thiele matrix, an <math>N×N</math> matrix where N is the number of reacting components. Overall, the method is shown to be relatively easy to implement and may significantly enhance the fidelity of chemical reactor simulations that employ multi-step reaction schemes.

engineering, chemical, environmental
Reversible Diffusion by Thermal Fluctuations

A. Donev,T. G. Fai,and E. Vanden-Eijnden

DOI: https://doi.org/10.48550/arXiv.1306.3158

2014-04-02

Abstract:A model for diffusion in liquids that couples the dynamics of tracer particles to a fluctuating Stokes equation for the fluid is investigated in the limit of large Schmidt number. In this limit, the concentration of tracers is shown to satisfy a closed-form stochastic advection-diffusion equation that is used to investigate the collective diffusion of hydrodynamically-correlated tracers through a combination of Eulerian and Lagrangian numerical methods. This analysis indicates that transport in liquids is quite distinct from the traditional Fickian picture of diffusion. While the ensemble-averaged concentration follows Fick's law with a diffusion coefficient that obeys the Stokes-Einstein relation, each instance of the diffusive mixing process exhibits long-ranged giant fluctuations around its average behavior. We construct a class of mesoscopic models for diffusion in liquids at different observation scales in which the renormalized diffusion coefficient depends on this scale. This indicates that the Fickian diffusion coefficient in liquids is not a material constant, but rather, changes with the scale at which experimental measurements are performed.

Statistical Mechanics,Mesoscale and Nanoscale Physics

A note on the standard diffusion curve of TAP analysis

An Alternate Mathematical Model for the Tap Technique for Numerical Integration

Reaction Diffusion TAP

TUNED FINITE-DIFFERENCE DIFFUSION OPERATORS

Boundary conditions in TAP reactor modeling

Thermally Induced Diffusion of Chemicals under Steady-State Heat Transfer in Saturated Porous Media

The measurement of gas diffusivity in porous materials by temporal analysis of products (TAP)

Numerical Simulation of Nanoparticle Synthesis in Diffusion Flame Reactor

Numerical study of anisotropic diffusion in Turing patterns based on Finsler geometry modeling

Uniform Convergent Tailored Finite Point Method for Advection–Diffusion Equation with Discontinuous, Anisotropic and Vanishing Diffusivity

Optimal finite-differences discretization for the diffusion equation from the perspective of large-deviation theory

Numerical Recovery of the Diffusion Coefficient in Diffusion Equations from Terminal Measurement

Identification of Concentration and Residence Time Distributions with One-Dimension Diffusion Model at Close Boundary Conditions

A High Accuracy Compact Difference Scheme for Convection Diffusion Reaction Equation on Non-uniform Grid

A Diffusion Scale Approximation for Stochastic Fission Chains with State Feedback

Application of the Discontinuity Factor Theory to Accelerate Routine Rod Worth Calculations for Modular HTRs

An iterative scheme to include turbulent diffusion in advective-dominated transport of delayed neutron precursors

A convergent discretisation method for transition path theory for diffusion processes

Beyond the effectiveness factor: Multi-step reactions with intraparticle diffusion limitations

Reversible Diffusion by Thermal Fluctuations