Exact solutions of certain nonlinear chemotaxis diffusion reaction equations
MISHRA AJAY; KAUSHAL R S; PRASAD AWADHESH
2016-05-01
Using the auxiliary equation method, we obtain exact solutions of certain nonlinear chemotaxis diffusion reaction equations in the presence of a stimulant. In particular, we account for the nonlinearities arising not only from the density-dependent source terms contributed by the particles and the stimulant but also from the coupling term of the stimulant. In addition to this, the diffusion of the stimulant and the effect of long-range interactions are also accounted for in theconstructed coupled differential equations. The results obtained here could be useful in the studies of several biological systems and processes, e.g., in bacterial infection, chemotherapy, etc.
MacDonald, G.; Mackenzie, J. A.; Nolan, M.; Insall, R. H.
2016-03-01
In this paper, we devise a moving mesh finite element method for the approximate solution of coupled bulk-surface reaction-diffusion equations on an evolving two dimensional domain. Fundamental to the success of the method is the robust generation of bulk and surface meshes. For this purpose, we use a novel moving mesh partial differential equation (MMPDE) approach. The developed method is applied to model problems with known analytical solutions; these experiments indicate second-order spatial and temporal accuracy. Coupled bulk-surface problems occur frequently in many areas; in particular, in the modelling of eukaryotic cell migration and chemotaxis. We apply the method to a model of the two-way interaction of a migrating cell in a chemotactic field, where the bulk region corresponds to the extracellular region and the surface to the cell membrane.
Chavanis, Pierre-Henri
2008-01-01
We consider a generalized class of Keller-Segel models describing the chemotaxis of biological populations (bacteria, amoebae, endothelial cells, social insects,...). We show the analogy with nonlinear mean field Fokker-Planck equations and generalized thermodynamics. As an illustration, we introduce a new model of chemotaxis incorporating both effects of anomalous diffusion and exclusion principle (volume filling). We also discuss the analogy between biological populations described by the Keller-Segel model and self-gravitating Brownian particles described by the Smoluchowski-Poisson system.
Effective Medium Equations for Chemotaxis in Porous Media
Valdes-Parada, F.; Porter, M.; Wood, B. D.; Narayanaswamy, K.; Ford, R.
2008-12-01
Biodegradation is an important mechanism for contaminant reduction in groundwater environments; in fact, in-situ bioremediation and bioaugmentation methods represent alternatives to traditional methods such as pump-and-treat. Chemotaxis has been shown to enhance bacterial transport toward or away from concentration gradients of chemical species in laboratory experiments and may signifficantly increase contaminant flux undergoing degradation at the interfaces of low- and high-permeability regions. In this work, the method of volume averaging is used to upscale the microscale description of chemotactic microbial transport in order to obtain the corresponding macroscale equations for bacteria and the chemoattractant. As a first apprach, cellular growth/death and consumption of the attractant by chemical reaction are assumed negligible with respect to convective and diffusive transport, in both levels of scale. For bacteria, two effective coefficients are introduced, namely a total motility tensor and an effective chemotactic sensitivity tensor. Both coefficients are computed by solving the associated closure problems in a capillary tube. Analysis of breakthrough curves resulting from numerical experiments is also presented.
Global Solutions to the Coupled Chemotaxis-Fluid Equations
Duan, Renjun
2010-08-10
In this paper, we are concerned with a model arising from biology, which is a coupled system of the chemotaxis equations and the viscous incompressible fluid equations through transport and external forcing. The global existence of solutions to the Cauchy problem is investigated under certain conditions. Precisely, for the Chemotaxis-Navier-Stokes system over three space dimensions, we obtain global existence and rates of convergence on classical solutions near constant states. When the fluid motion is described by the simpler Stokes equations, we prove global existence of weak solutions in two space dimensions for cell density with finite mass, first-order spatial moment and entropy provided that the external forcing is weak or the substrate concentration is small. © Taylor & Francis Group, LLC.
Travelling Waves in Hyperbolic Chemotaxis Equations
Xue, Chuan
2010-10-16
Mathematical models of bacterial populations are often written as systems of partial differential equations for the densities of bacteria and concentrations of extracellular (signal) chemicals. This approach has been employed since the seminal work of Keller and Segel in the 1970s (Keller and Segel, J. Theor. Biol. 30:235-248, 1971). The system has been shown to permit travelling wave solutions which correspond to travelling band formation in bacterial colonies, yet only under specific criteria, such as a singularity in the chemotactic sensitivity function as the signal approaches zero. Such a singularity generates infinite macroscopic velocities which are biologically unrealistic. In this paper, we formulate a model that takes into consideration relevant details of the intracellular processes while avoiding the singularity in the chemotactic sensitivity. We prove the global existence of solutions and then show the existence of travelling wave solutions both numerically and analytically. © 2010 Society for Mathematical Biology.
Boundedness in a chemotaxis-haptotaxis model with nonlinear diffusion
Li, Yan; Lankeit, Johannes
2016-05-01
This article deals with an initial-boundary value problem for the coupled chemotaxis-haptotaxis system with nonlinear diffusion under homogeneous Neumann boundary conditions in a bounded smooth domain Ω \\subset {{{R}}n} , n = 2, 3, 4, where χ,ξ and μ are given nonnegative parameters. The diffusivity D(u) is assumed to satisfy D(u)≥slant δ {{u}m-1} for all u > 0 with some δ >0 . It is proved that for sufficiently regular initial data global bounded solutions exist whenever m>2-\\frac{2}{n} . For the case of non-degenerate diffusion (i.e. D(0) > 0) the solutions are classical; for the case of possibly degenerate diffusion (D(0)≥slant 0 ), the existence of bounded weak solutions is shown.
Speed ot travelling waves in reaction-diffusion equations
Benguria, R D; Méndez, V
2002-01-01
Reaction diffusion equations arise in several problems of population dynamics, flame propagation and others. In one dimensional cases the systems may evolve into travelling fronts. Here we concentrate on a reaction diffusion equation which arises as a simple model for chemotaxis and present results for the speed of the travelling fronts. (Author)
Speed ot travelling waves in reaction-diffusion equations
Reaction diffusion equations arise in several problems of population dynamics, flame propagation and others. In one dimensional cases the systems may evolve into travelling fronts. Here we concentrate on a reaction diffusion equation which arises as a simple model for chemotaxis and present results for the speed of the travelling fronts. (Author)
Speed ot travelling waves in reaction-diffusion equations
Benguria, R.D.; Depassier, M.C. [Facultad de Fisica, Pontificia Universidad Catolica de Chile, Avda. Vicuna Mackenna 4860, Santiago (Chile); Mendez, V. [Facultat de Ciencies de la Salut, Universidad Internacional de Catalunya, Gomera s/n 08190 Sant Cugat del Valles, Barcelona (Spain)
2002-07-01
Reaction diffusion equations arise in several problems of population dynamics, flame propagation and others. In one dimensional cases the systems may evolve into travelling fronts. Here we concentrate on a reaction diffusion equation which arises as a simple model for chemotaxis and present results for the speed of the travelling fronts. (Author)
Individual-based models for bacterial chemotaxis in the diffusion asymptotics
Rousset, Mathias
2011-01-01
We discuss velocity-jump models for chemotaxis of bacteria with an internal state that allows the velocity jump rate to depend on the memory of the chemoattractant concentration along their path of motion. Using probabilistic techniques, we provide a pathwise result that shows that the considered process converges to an advection-diffusion process in the (long-time) diffusion limit. We also (re-)prove using the same approach that the same limiting equation arises for a related, simpler process with direct sensing of the chemoattractant gradient. Additionally, we propose a time discretization technique that retains these diffusion limits exactly, i.e., without error that depends on the time discretization. In the companion paper \\cite{variance}, these results are used to construct a coupling technique that allows numerical simulation of the process with internal state with asymptotic variance reduction, in the sense that the variance vanishes in the diffusion limit.
Numerical study of plume patterns in the chemotaxis-diffusion-convection coupling system
Deleuze, Yannick; Thiriet, Marc; Sheu, Tony W H
2015-01-01
A chemotaxis-diffusion-convection coupling system for describing a form of buoyant convection in which the fluid develops convection cells and plume patterns will be investigated numerically in this study. Based on the two-dimensional convective chemotaxis-fluid model proposed in the literature, we developed an upwind finite element method to investigate the pattern formation and the hydrodynamical stability of the system. The numerical simulations illustrate different predicted physical regimes in the system. In the convective regime, the predicted plumes resemble B\\'enard instabilities. Our numerical results show how structured layers of bacteria are formed before bacterium rich plumes fall in the fluid. The plumes have a well defined spectrum of wavelengths and have an exponential growth rate, yet their position can only be predicted in very simple examples. In the chemotactic and diffusive regimes, the effects of chemotaxis are investigated. Our results indicate that the chemotaxis can stabilize the overa...
Multinomial diffusion equation
Balter, Ariel; Tartakovsky, Alexandre M.
2011-06-01
We describe a new, microscopic model for diffusion that captures diffusion induced fluctuations at scales where the concept of concentration gives way to discrete particles. We show that in the limit as the number of particles N→∞, our model is equivalent to the classical stochastic diffusion equation (SDE). We test our new model and the SDE against Langevin dynamics in numerical simulations, and show that our model successfully reproduces the correct ensemble statistics, while the classical model fails.
Ibrahim, Moustafa; Saad, Mazen
2015-01-01
This paper is devoted to the mathematical analysis of a degenerate nonlinear parabolic equation. This kind of equations stems either from the modeling of a compressible two phase flow in porous media or from the modeling of a chemotaxis-fluid process. In the degenerate equation, the strong nonlinearities are technically difficult to be controlled by the degenerate dissipative term because the equation itself presents degenerate terms of order 0 and of order 1. In the case of the sub-quadratic...
Planar Gradient Diffusion System to Investigate Chemotaxis in a 3D Collagen Matrix.
Stout, David A; Toyjanova, Jennet; Franck, Christian
2015-01-01
The importance of cell migration can be seen through the development of human life. When cells migrate, they generate forces and transfer these forces to their surrounding area, leading to cell movement and migration. In order to understand the mechanisms that can alter and/or affect cell migration, one can study these forces. In theory, understanding the fundamental mechanisms and forces underlying cell migration holds the promise of effective approaches for treating diseases and promoting cellular transplantation. Unfortunately, modern chemotaxis chambers that have been developed are usually restricted to two dimensions (2D) and have complex diffusion gradients that make the experiment difficult to interpret. To this end, we have developed, and describe in this paper, a direct-viewing chamber for chemotaxis studies, which allows one to overcome modern chemotaxis chamber obstacles able to measure cell forces and specific concentration within the chamber in a 3D environment to study cell 3D migration. More compelling, this approach allows one to successfully model diffusion through 3D collagen matrices and calculate the coefficient of diffusion of a chemoattractant through multiple different concentrations of collagen, while keeping the system simple and user friendly for traction force microscopy (TFM) and digital volume correlation (DVC) analysis. PMID:26131645
Reaction Diffusion and Chemotaxis for Decentralized Gathering on FPGAs
Bernard Girau
2009-01-01
and rapid simulations of the complex dynamics of this reaction-diffusion model. Then we describe the FPGA implementation of the environment together with the agents, to study the major challenges that must be solved when designing a fast embedded implementation of the decentralized gathering model. We analyze the results according to the different goals of these hardware implementations.
Marangoni-driven chemotaxis, chemotactic collapse, and the Keller-Segel equation
Shelley, Michael; Masoud, Hassan
2013-11-01
Almost by definition, chemotaxis involves the biased motion of motile particles along gradients of a chemical concentration field. Perhaps the most famous model for collective chemotaxis in mathematical biology is the Keller-Segel model, conceived to describe collective aggregation of slime mold colonies in response to an intrinsically produced, and diffusing, chemo-attractant. Heavily studied, particularly in 2D where the system is ``super-critical'', it has been proved that the KS model can develop finite-time singularities - so-called chemotactic collapse - of delta-function type. Here, we study the collective dynamics of immotile particles bound to a 2D interface above a 3D fluid. These particles are chemically active and produce a diffusing field that creates surface-tension gradients along the surface. The resultant Marangoni stresses create flows that carry the particles, possibly concentrating them. Remarkably, we show that this system involving 3D diffusion and fluid dynamics, exactly yields the 2D Keller-Segel model for the surface-flow of active particles. We discuss the consequences of collapse on the 3D fluid dynamics, and generalizations of the fluid-dynamical model.
Diffusion equations and turbulent transport
One scrutinized transport equations differing essentially in form from the classical diffusion one. Description of diffusion under strong nonequilibrium and turbulence involved application of equations that took account of transport nonlocality and memory effects. One analyzed ways to derive the mentioned equations starting from quasi-linear approximation and up to equations with fractional derivatives. One points out the generality of the applied theoretical concepts in spite of the essential difference of the exact physical problems. One demonstrated the way of application of the theoretical and probabilistic ideas
Diffusion equations and turbulent transport
Diffusion equations are considered that differ substantially in structure from classical ones. A description of diffusion under strongly nonequilibrium conditions in a highly turbulent plasma requires the use of equations that take into account memory effects and the nonlocal nature of transport. Different methods are developed for constructing such equations, ranging from those in the quasilinear approximation to those with fractional derivatives. It is emphasized that the theoretical concepts underlying the equations proposed are common for a very wide variety of specific physical problems. The ways of applying theoretical probabilistic ideas are demonstrated
The generalized Airy diffusion equation
Frank M. Cholewinski
2003-08-01
Full Text Available Solutions of a generalized Airy diffusion equation and an associated nonlinear partial differential equation are obtained. Trigonometric type functions are derived for a third order generalized radial Euler type operator. An associated complex variable theory and generalized Cauchy-Euler equations are obtained. Further, it is shown that the Airy expansions can be mapped onto the Bessel Calculus of Bochner, Cholewinski and Haimo.
A coupled chemotaxis-fluid model: Global existence
Liu, Jian-Guo
2011-09-01
We consider a model arising from biology, consisting of chemotaxis equations coupled to viscous incompressible fluid equations through transport and external forcing. Global existence of solutions to the Cauchy problem is investigated under certain conditions. Precisely, for the chemotaxis-Navier- Stokes system in two space dimensions, we obtain global existence for large data. In three space dimensions, we prove global existence of weak solutions for the chemotaxis-Stokes system with nonlinear diffusion for the cell density.© 2011 Elsevier Masson SAS. All rights reserved.
Parabolic equations in biology growth, reaction, movement and diffusion
Perthame, Benoît
2015-01-01
This book presents several fundamental questions in mathematical biology such as Turing instability, pattern formation, reaction-diffusion systems, invasion waves and Fokker-Planck equations. These are classical modeling tools for mathematical biology with applications to ecology and population dynamics, the neurosciences, enzymatic reactions, chemotaxis, invasion waves etc. The book presents these aspects from a mathematical perspective, with the aim of identifying those qualitative properties of the models that are relevant for biological applications. To do so, it uncovers the mechanisms at work behind Turing instability, pattern formation and invasion waves. This involves several mathematical tools, such as stability and instability analysis, blow-up in finite time, asymptotic methods and relative entropy properties. Given the content presented, the book is well suited as a textbook for master-level coursework.
Perturbative linearization of reaction-diffusion equations
We develop perturbative expansions to obtain solutions for the initial-value problems of two important reaction-diffusion systems, namely the Fisher equation and the time-dependent Ginzburg-Landau equation. The starting point of our expansion is the corresponding singular-perturbation solution. This approach transforms the solution of nonlinear reaction-diffusion equations into the solution of a hierarchy of linear equations. Our numerical results demonstrate that this hierarchy rapidly converges to the exact solution
A System of Non-linear Partial Differential Equations Modeling Chemotaxis with Sensitivity Functions
Post, Katharina
1999-01-01
Wir betrachten ein System nichtlinearer parabolischer partieller Differentialgleichungen zur Modellierung des biologischen Phänomens Chemotaxis, das unter anderem in Aggregationsprozessen in Lebenszyklen bestimmter Einzeller eine wichtige Rolle spielt. Unser Chemotaxismodell benutzt Sensitivitäts funktionen, die die vorkommenden biologischen Prozesse genauer spezifizieren. Trotz der durch die Sensitivitätsfunktionen eingebrachten, zusätzlichen Nichtlinearitäten in den Gleichungen erhalten w...
Fractional diffusion equations coupled by reaction terms
Lenzi, E. K.; Menechini Neto, R.; Tateishi, A. A.; Lenzi, M. K.; Ribeiro, H. V.
2016-09-01
We investigate the behavior for a set of fractional reaction-diffusion equations that extend the usual ones by the presence of spatial fractional derivatives of distributed order in the diffusive term. These equations are coupled via the reaction terms which may represent reversible or irreversible processes. For these equations, we find exact solutions and show that the spreading of the distributions is asymptotically governed by the same the long-tailed distribution. Furthermore, we observe that the coupling introduced by reaction terms creates an interplay between different diffusive regimes leading us to a rich class of behaviors related to anomalous diffusion.
The Nonlinear Convection—Reaction—Diffusion Equation
ShiminTANG; MaochangCUI; 等
1996-01-01
A nonlinear convection-reaction-diffusion equation is used as a model equation of the El Nino events.In this model,the effects of convection,turbulent diffusion,linear feed-back and nolinear radiation on the anomaly of Sea Surface Temperature(SST) are considered.In the case of constant convection,this equation has exact kink-like travelling wave solutions,which can be used to explain the history of an El Nino event.
Fractional Diffusion Limit for Collisional Kinetic Equations
Mellet, Antoine
2010-08-20
This paper is devoted to diffusion limits of linear Boltzmann equations. When the equilibrium distribution function is a Maxwellian distribution, it is well known that for an appropriate time scale, the small mean free path limit gives rise to a diffusion equation. In this paper, we consider situations in which the equilibrium distribution function is a heavy-tailed distribution with infinite variance. We then show that for an appropriate time scale, the small mean free path limit gives rise to a fractional diffusion equation. © 2010 Springer-Verlag.
Diffusive limits for linear transport equations
The authors show that the Hibert and Chapman-Enskog asymptotic treatments that reduce the nonlinear Boltzmann equation to the Euler and Navier-Stokes fluid equations have analogs in linear transport theory. In this linear setting, these fluid limits are described by diffusion equations, involving familiar and less familiar diffusion coefficients. Because of the linearity extant, one can carry out explicitly the initial and boundary layer analyses required to obtain asymptotically consistent initial and boundary conditions for the diffusion equations. In particular, the effects of boundary curvature and boundary condition variation along the surface can be included in the boundary layer analysis. A brief review of heuristic (nonasymptotic) diffusion description derivations is also included in our discussion
Global Attractors for a Nonclassical Diffusion Equation
Chun You SUN; Su Yun WANG; Cheng Kui ZHONG
2007-01-01
We prove the existence of global attractors in H10 (Ω) for a nonclassical diffusion equation.Two types of nonlinearity f are considered: one is the critical exponent, and the other is the polynomial growth of arbitrary order.
Markowich, Peter
2010-06-01
We study the system ct + u · ∇c = ∇c -nf(c) nt + u · ∇n = ∇n m - ∇ · (n×(c) ∇c) ut + u·∇u + ∇P - η∇u + n∇φ/ = 0 ∇·u = 0. arising in the modelling of the motion of swimming bacteria under the effect of diffusion, oxygen-taxis and transport through an incompressible fluid. The novelty with respect to previous papers in the literature lies in the presence of nonlinear porous-medium-like diffusion in the equation for the density n of the bacteria, motivated by a finite size effect. We prove that, under the constraint m ε (3/2, 2] for the adiabatic exponent, such system features global in time solutions in two space dimensions for large data. Moreover, in the case m = 2 we prove that solutions converge to constant states in the large-time limit. The proofs rely on standard energy methods and on a basic entropy estimate which cannot be achieved in the case m = 1. The case m = 2 is very special as we can provide a Lyapounov functional. We generalize our results to the three-dimensional case and obtain a smaller range of exponents m ε (m*, 2] with m* > 3/2, due to the use of classical Sobolev inequalities.
Linearization of Systems of Nonlinear Diffusion Equations
KANG Jing; QU Chang-Zheng
2007-01-01
We investigate the linearization of systems of n-component nonlinear diffusion equations; such systems have physical applications in soil science, mathematical biology and invariant curve flows. Equivalence transformations of their auxiliary systems are used to identify the systems that can be linearized. We also provide several examples of systems with two-component equations, and show how to linearize them by nonlocal mappings.
Diffusion phenomenon for linear dissipative wave equations
Said-Houari, Belkacem
2012-01-01
In this paper we prove the diffusion phenomenon for the linear wave equation. To derive the diffusion phenomenon, a new method is used. In fact, for initial data in some weighted spaces, we prove that for {equation presented} decays with the rate {equation presented} [0,1] faster than that of either u or v, where u is the solution of the linear wave equation with initial data {equation presented} [0,1], and v is the solution of the related heat equation with initial data v 0 = u 0 + u 1. This result improves the result in H. Yang and A. Milani [Bull. Sci. Math. 124 (2000), 415-433] in the sense that, under the above restriction on the initial data, the decay rate given in that paper can be improved by t -γ/2. © European Mathematical Society.
Reaction diffusion equations with boundary degeneracy
Huashui Zhan
2016-03-01
Full Text Available In this article, we consider the reaction diffusion equation $$ \\frac{\\partial u}{\\partial t} = \\Delta A(u,\\quad (x,t\\in \\Omega \\times (0,T, $$ with the homogeneous boundary condition. Inspired by the Fichera-Oleinik theory, if the equation is not only strongly degenerate in the interior of $\\Omega$, but also degenerate on the boundary, we show that the solution of the equation is free from any limitation of the boundary condition.
A Diffusion Equation on Fractals in Random Media
DeLIU; HouqiangLI; 等
1997-01-01
The proper transformation method is used in this paper to extend the diffusion equation on Euclidean space to the standard diffusion equation on fractals,By this standard diffusion equation,it is proved that the fractal Brownian particle moving belongs to the anomalous diffusion.At the same time,the generalized diffusion equation and its asymptotic solution is discussed.
Relativistic diffusion equation from stochastic quantization
Kazinski, P O
2007-01-01
The new scheme of stochastic quantization is proposed. This quantization procedure is equivalent to the deformation of an algebra of observables in the manner of deformation quantization with an imaginary deformation parameter (the Planck constant). We apply this method to the models of nonrelativistic and relativistic particles interacting with an electromagnetic field. In the first case we establish the equivalence of such a quantization to the Fokker-Planck equation with a special force. The application of the proposed quantization procedure to the model of a relativistic particle results in a relativistic generalization of the Fokker-Planck equation in the coordinate space, which in the absence of the electromagnetic field reduces to the relativistic diffusion (heat) equation. The stationary probability distribution functions for a stochastically quantized particle diffusing under a barrier and a particle in the potential of a harmonic oscillator are derived.
Logarithmic diffusion and porous media equations: a unified description
Pedron, I. T.; Mendes, R. S.; Buratta, T. J.; L. C. Malacarne; Lenzi, E. K.
2005-01-01
In this work we present the logarithmic diffusion equation as a limit case when the index that characterizes a nonlinear Fokker-Planck equation, in its diffusive term, goes to zero. A linear drift and a source term are considered in this equation. Its solution has a lorentzian form, consequently this equation characterizes a super diffusion like a L\\'evy kind. In addition is obtained an equation that unifies the porous media and the logarithmic diffusion equations, including a generalized dif...
Wang, Chi-Jen [Iowa State Univ., Ames, IA (United States)
2013-01-01
In this thesis, we analyze both the spatiotemporal behavior of: (A) non-linear “reaction” models utilizing (discrete) reaction-diffusion equations; and (B) spatial transport problems on surfaces and in nanopores utilizing the relevant (continuum) diffusion or Fokker-Planck equations. Thus, there are some common themes in these studies, as they all involve partial differential equations or their discrete analogues which incorporate a description of diffusion-type processes. However, there are also some qualitative differences, as shall be discussed below.
Existence of global solutions for a chemotaxis-fluid system with nonlinear diffusion
Chung, Yun-Sung; Kang, Kyungkeun
2016-04-01
We consider a coupled system consisting of the Navier-Stokes equations and a porous medium type of Keller-Segel system that model the motion of swimming bacteria living in fluid and consuming oxygen. We establish the global-in-time existence of weak solutions for the Cauchy problem of the system in dimension three. In addition, if the Stokes system, instead Navier-Stokes system, is considered for the fluid equation, we prove that bounded weak solutions exist globally in time.
A Convergent Reaction-Diffusion Master Equation
Isaacson, Samuel A
2012-01-01
The reaction-diffusion master equation (RDME) is a lattice stochastic reaction-diffusion model that has been used to study spatially distributed cellular processes. The RDME has been shown to have the drawback of losing bimolecular reactions in the continuum limit that the lattice spacing approaches zero (in two or more dimensions). In this work we derive a new convergent RDME (CRDME) that eliminates this problem. The CRDME is obtained by finite volume discretization of a spatially-continuous stochastic reaction-diffusion model. We demonstrate the empirical numerical convergence of reaction time statistics associated with the CRDME. Although the reaction time statistics of the RDME diverge as the lattice spacing approaches zero, we show they approach those of the CRDME for sufficiently large lattice spacings or slow bimolecular reaction rates. As such, the RDME may be interpreted as an approximation to the CRDME in several asymptotic limits.
Entropy methods for diffusive partial differential equations
Jüngel, Ansgar
2016-01-01
This book presents a range of entropy methods for diffusive PDEs devised by many researchers in the course of the past few decades, which allow us to understand the qualitative behavior of solutions to diffusive equations (and Markov diffusion processes). Applications include the large-time asymptotics of solutions, the derivation of convex Sobolev inequalities, the existence and uniqueness of weak solutions, and the analysis of discrete and geometric structures of the PDEs. The purpose of the book is to provide readers an introduction to selected entropy methods that can be found in the research literature. In order to highlight the core concepts, the results are not stated in the widest generality and most of the arguments are only formal (in the sense that the functional setting is not specified or sufficient regularity is supposed). The text is also suitable for advanced master and PhD students and could serve as a textbook for special courses and seminars.
Self-similar dynamics of bacterial chemotaxis
Ngamsaad, Waipot
2012-01-01
We investigate the pattern formation of colony generated by chemotactic bacteria through a continuum model. In a simplified case, the dynamics of system is governed by a density-dependent convection-reaction-diffusion equation, $u_t = (u^{m})_{xx} - 2\\kappa(u^m)_{x}+ u - u^{m}$. This equation admits the analytical solutions that show the self-similarity of the bacterial colony's morphogenesis. In addition, we found that the colony evolves long time as the sharp traveling wave. The roles of chemotaxis on the regulation of pattern formation in these results are also discussed.
Simplification and Application of Boyd Membrane-diffusion Equation
甄捷
2004-01-01
The havirng-been-used-for-50-year Boyd membrane diffusion Equation - In( 1 - F) = R t can be deduced into F = kt through using Maclanrin expansion equation and the Lagerange remainders. The latter is a simple membrane diffusion equation, which is available to judge if the exchanging course of the resin obeys the rules of membrane-diffusion mechanism more conveniently.
Turbulence and diffusion scaling versus equations
Bakunin, Oleg G
2008-01-01
This book is an introduction to the multidisciplinary field of anomalous diffusion in complex systems, with emphasis on the scaling approach as opposed to techniques based on the quantitative analysis of underlying transport equations. Typical examples of such systems are turbulent plasmas, convective rolls, zonal flow systems and stochastic magnetic fields. From the more methodological point of view, the approach relies on the general use of correlations estimates, quasilinear equations and continuous time random walk techniques. Yet, the mathematical descriptions are not meant to become a fixed set of recipes but rather develop and strengthen the reader's physical intuition and understanding on the underlying mechanisms involved. Most of the material stems from class-tested lectures, where graduate students where assumed to have a working knowledge of classical physics, fluid dynamics and plasma physics but otherwise no prior knowledge of the subject matter is assumed from the side of the reader.
Differential constraints and exact solutions of nonlinear diffusion equations
The differential constraints are applied to obtain explicit solutions of nonlinear diffusion equations. Certain linear determining equations with parameters are used to find such differential constraints. They generalize the determining equations used in the search for classical Lie symmetries
Differential constraints and exact solutions of nonlinear diffusion equations
Kaptsov, Oleg V.; Verevkin, Igor V.
2002-01-01
The differential constraints are applied to obtain explicit solutions of nonlinear diffusion equations. Certain linear determining equations with parameters are used to find such differential constraints. They generalize the determining equations used in the search for classical Lie symmetries.
Nikitin, A. G.
2005-01-01
Group classification of systems of two coupled nonlinear reaction-diffusion equation with a general diffusion matrix started in papers math-ph/0411027,, math-ph/0411028 is completed in present paper where all non-equivalent equations with triangular diffusion matrix are classified. In addition, symmetries of diffusion systems with nilpotent diffusion matrix and additional first order derivative terms are described.
The solution of integral equation of plasma diffusion problems
Two topics which concern with the solution of one dimensional diffusion equation are discussed in this paper. The solution of the one dimensional diffusion equation with one initial condition is discussed in the first topic, and the solution of the same equation with two initial conditions is discussed in the second topic. Both solutions are presented in integral forms. (author)
Neutron transport equation - indications on homogenization and neutron diffusion
In PWR nuclear reactor, the practical study of the neutrons in the core uses diffusion equation to describe the problem. On the other hand, the most correct method to describe these neutrons is to use the Boltzmann equation, or neutron transport equation. In this paper, we give some theoretical indications to obtain a diffusion equation from the general transport equation, with some simplifying hypothesis. The work is organised as follows: (a) the most general formulations of the transport equation are presented: integro-differential equation and integral equation; (b) the theoretical approximation of this Boltzmann equation by a diffusion equation is introduced, by the way of asymptotic developments; (c) practical homogenization methods of transport equation is then presented. In particular, the relationships with some general and useful methods in neutronic are shown, and some homogenization methods in energy and space are indicated. A lot of other points of view or complements are detailed in the text or the remarks
Finite formulation of one-group neutron diffusion equation
The objective of this paper is to introduce a new numerical method for neutron diffusion equation. This method produces the final algebraic form of the neutron diffusion equation by classifying the neutronic variables and using two kinds of cell complexes (primal and dual cells) without using the conventional differential form of the neutron diffusion equation. The convergence of the method is exactly same as finite element method with linear interpolation. (author)
On an explicit finite difference method for fractional diffusion equations
S. B. Yuste; Acedo, L.
2003-01-01
A numerical method to solve the fractional diffusion equation, which could also be easily extended to many other fractional dynamics equations, is considered. These fractional equations have been proposed in order to describe anomalous transport characterized by non-Markovian kinetics and the breakdown of Fick's law. In this paper we combine the forward time centered space (FTCS) method, well known for the numerical integration of ordinary diffusion equations, with the Grunwald-Letnikov defin...
Voter Model Perturbations and Reaction Diffusion Equations
Cox, J Theodore; Perkins, Edwin
2011-01-01
We consider particle systems that are perturbations of the voter model and show that when space and time are rescaled the system converges to a solution of a reaction diffusion equation in dimensions $d \\ge 3$. Combining this result with properties of the PDE, some methods arising from a low density super-Brownian limit theorem, and a block construction, we give general, and often asymptotically sharp, conditions for the existence of non-trivial stationary distributions, and for extinction of one type. As applications, we describe the phase diagrams of three systems when the parameters are close to the voter model: (i) a stochastic spatial Lotka-Volterra model of Neuhauser and Pacala, (ii) a model of the evolution of cooperation of Ohtsuki, Hauert, Lieberman, and Nowak, and (iii) a continuous time version of the non-linear voter model of Molofsky, Durrett, Dushoff, Griffeath, and Levin. The first application confirms a conjecture of Cox and Perkins and the second confirms a conjecture of Ohtsuki et al in the ...
Diffusion equations and the time evolution of foreign exchange rates
We investigate which type of diffusion equation is most appropriate to describe the time evolution of foreign exchange rates. We modify the geometric diffusion model assuming a non-exponential time evolution and the stochastic term is the sum of a Wiener noise and a jump process. We find the resulting diffusion equation to obey the Kramers–Moyal equation. Analytical solutions are obtained using the characteristic function formalism and compared with empirical data. The analysis focus on the first four central moments considering the returns of foreign exchange rate. It is shown that the proposed model offers a good improvement over the classical geometric diffusion model.
The Non-Classical Boltzmann Equation, and Diffusion-Based Approximations to the Boltzmann Equation
Frank, Martin; Larsen, Edward W; Vasques, Richard
2014-01-01
We show that several diffusion-based approximations (classical diffusion or SP1, SP2, SP3) to the linear Boltzmann equation can (for an infinite, homogeneous medium) be represented exactly by a non-classical transport equation. As a consequence, we indicate a method to solve diffusion-based approximations to the Boltzmann equation via Monte Carlo, with only statistical errors - no truncation errors.
Partial differential equations and non-diffusive structures
In this paper we give a short introduction to open problems and recent studies of classes of partial differential equations, which—in contrast to reaction–diffusion systems—describe phenomena with local interactions. Partial differential equations coupled with ordinary differential equations, models of transport type and hyperbolic systems are discussed with respect to their pattern forming behaviour. (open problem)
Solution of diffusive wave equation using fem for flood forecasting
Flood forecasting can be predicted numerically by the solution of 1D (One-Dimensional) unsteady diffusive wave equation. This research paper presents the development of a fem (Finite Element Model) for flood routing using diffusive wave with and without lateral inflow/outflow. FEM is based on two-step semi implicit Taylor Galerkin technique. The accuracy of the model has been verified by comparing computed results with available solution of the diffusive wave problems available in open literature. Numerical results demonstrate that the technique is an efficient and accurate tool to simulate diffusive wave equation for outflow hydrograph for temporal variation of lateral inflow and outflow. (author)
Numerical solution of the diffusion equation. Part 1
The finite difference method as an approximation to the analytic solution of the diffusion equation is described briefly. For higher dimension this finite difference method will result in linear equation system which can be further solved using various techniques. Some techniques of solving this linear equation are also discussed. The main problem encountered is not that of replacing the governing partial differential equation by the finite difference equations, but that of finding a solution to the system of linear equation formed by the approximation. (author)
Two lines in improved neutron diffusion theory extending the efficiency of finite-difference diffusion codes to the field of optically small systems, are here reviewed. The firs involves the nodal solution for tensorial diffusion equation in slab geometry and tensorial formulation in parallelepiped and cylindrical gemometry; the dependence of critical eigenvalue from small slab thicknesses is also analitically investigated and finally a regularized tensorial diffusion equation is derived for slab. The other line refer to diffusion models formally unchanged with respect to the classical one, but where new size-dependent RTGB definitions for diffusion parameters are adopted, requiring that they allow to reproduce, in diffusion approach, the terms of neutron transport global balance; the trascendental equation for the buckling, arising in slab, sphere and parallelepiped geometry from the above requirement, are reported and the sizedependence of the new diffusion coefficient and extrapolated end point is investigated
LAGRANGE STABILITY IN MEAN SQUARE OF STOCHASTIC REACTION DIFFUSION EQUATIONS
无
2006-01-01
This work is devoted to the discussion of stochastic reaction diffusion equations and some new theorems on Lagrange stability in mean square of the solution are established via Lyapunov method which is nothing to be done in the past.
Biomixing by chemotaxis and enhancement of biological reactions
Kiselev, Alexander
2011-01-01
Many processes in biology involve both reactions and chemotaxis. However, to the best of our knowledge, the question of interaction between chemotaxis and reactions has not yet been addressed either analytically or numerically. We consider a model with a single density function involving diffusion, advection, chemotaxis, and absorbing reaction. The model is motivated, in particular, by studies of coral broadcast spawning, where experimental observations of the efficiency of fertilization rates significantly exceed the data obtained from numerical models that do not take chemotaxis (attraction of sperm gametes by a chemical secreted by egg gametes) into account. We prove that in the framework of our model, chemotaxis plays a crucial role. There is a rigid limit to how much the fertilization efficiency can be enhanced if there is no chemotaxis but only advection and diffusion. On the other hand, when chemotaxis is present, the fertilization rate can be arbitrarily close to being complete provided that the chemo...
Solutions of fractional diffusion equations by variation of parameters method
Mohyud-Din Syed Tauseef
2015-01-01
Full Text Available This article is devoted to establish a novel analytical solution scheme for the fractional diffusion equations. Caputo’s formulation followed by the variation of parameters method has been employed to obtain the analytical solutions. Following the derived analytical scheme, solution of the fractional diffusion equation for several initial functions has been obtained. Graphs are plotted to see the physical behavior of obtained solutions.
Measure solutions for the Smoluchowski coagulation-diffusion equation
Norris, James
2014-01-01
A notion of measure solution is formulated for a coagulation-diffusion equation, which is the natural counterpart of Smoluchowski's coagulation equation in a spatially inhomogeneous setting. Some general properties of such solutions are established. Sufficient conditions are identified on the diffusivity, coagulation rates and initial data for existence, uniqueness and mass conservation of solutions. These conditions impose no form of monotonicity on the coagulation kernel, which may depend o...
Hyperbolic reaction-diffusion equations for a forest fire model
Méndez López, Vicenç; Llebot, Josep Enric,
1997-01-01
Forest fire models have been widely studied from the context of self-organized criticality and from the ecological properties of the forest and combustion. On the other hand, reaction-diffusion equations have interesting applications in biology and physics. We propose here a model for fire propagation in a forest by using hyperbolic reaction-diffusion equations. The dynamical and thermodynamical aspects of the model are analyzed in detail.
The diffusion equation and the steady state. Chapter 2
We shall now study the equations that govern the neutron field in a reactor. These equations are based on the concept of local neutron balance, which takes into account the reaction rates in an element of volume and the net leakage rates out of the volume. The reaction rates are written in terms of the local cross sections, assumed known from a preprocessed database (e.g., ENDF/B-VI). The starting equation is the Maxwell-Boltzmann transport equation, in its integro-differential form. The various approximations required to go from the transport equation to the neutron diffusion equation will be presented first, because all finite-reactor calculations are based on the diffusion approximation. We shall then discuss the multi-group formalism of the diffusion equations and study the mathematical properties of this equation in steady state. This preliminary step will allow us to derive in a more accurate way, in the next chapter, the reactor point-kinetics equations. In the diffusion approximation, neutrons diffuse from regions of high concentration to regions of low concentration, just as heat diffuses from regions of high temperature to those of low temperature, or, rather, as gas molecules diffuse to reduce spatial variations in concentration. While it is sufficiently accurate to treat the transport of gas molecules as a diffusion process, this approach is too limiting for neutron transport. In contrast to a gas, where collisions are very frequent, the cross sections for the interaction of neutrons with nuclei are relatively small, as we saw in chapter 1 (of the order of barns, i.e., 10-24cm2) . This implies that neutrons traverse appreciable distances (of the order of a centimetre) between collisions. This relatively long neutron mean free path, together with the heterogeneity of the physical medium, requires that a more complete treatment be carried out, taking account of variations in the angular distribution of neutron speed in the vicinity of highly absorbing
Notes on Stefan-Maxwell Equation versus Grahan's Diffusion Law
无
2000-01-01
Certain prerequisite information on the component fluxes is necessary for solution of the Stefan-Maxwell equation in multicomponent diffusion systems and the Graham's law of diffusion and effusion is often resorted for this purpose. This article addresses solution of the Stefan-Maxwell equation in binary gas systems and explores the necessary conditions for definite solution of concentration profiles and pertinent component fluxes. It is found that there are multiple solutions for component fluxes in contradiction to what specified by the Graham's law of diffusion. The theorem of minimum entropy production in the non-equilibrium thermodynamics is believed instructive in determining the stable steady state solution out of infinite multiple solutions possible under the specified conditions. It is suggested that only when the boundary condition of component concentration is symmetrical in an isothermal binary system, the counter-diffusion becomes equimolar. The Graham's law of diffusion seems not generally valid for the case of isothermal ordinary diffusion.
Some Remarks on Some Strongly Coupled Reaction-Diffusion Equations
Diagana, Toka
2003-01-01
The primary goal of this paper is to characterize solutions to coupled reaction-diffusion systems. Indeed, we use operators theory to show that under suitable assumptions, then the solutions to the reaction-diffusion equations exist. As applications, we consider a mathematical model arising in Biology and in Chemistry.
A Fluctuating Lattice Boltzmann Method for the Diffusion Equation
Wagner, Alexander J
2016-01-01
We derive a fluctuating lattice Boltzmann method for the diffusion equation. The derivation removes several shortcomings of previous derivations for fluctuating lattice Boltzmann methods for hydrodynamic systems. The comparative simplicity of this diffusive system highlights the basic features of this first exact derivation of a fluctuating lattice Boltzmann method.
Chechkin, A. V.; Gorenflo, R.; I. M. Sokolov
2002-01-01
We propose diffusion-like equations with time and space fractional derivatives of the distributed order for the kinetic description of anomalous diffusion and relaxation phenomena, whose diffusion exponent varies with time and which, correspondingly, can not be viewed as self-affine random processes possessing a unique Hurst exponent. We prove the positivity of the solutions of the proposed equations and establish the relation to the Continuous Time Random Walk theory. We show that the distri...
Ranjit Kumar
2012-09-01
Travelling and solitary wave solutions of certain coupled nonlinear diffusion-reaction equations have been constructed using the auxiliary equation method. These equations arise in a variety of contexts not only in biological, chemical and physical sciences but also in ecological and social sciences.
Higher Order and Fractional Diffusive Equations
D. Assante
2015-07-01
Full Text Available We discuss the solution of various generalized forms of the Heat Equation, by means of different tools ranging from the use of Hermite-Kampé de Fériet polynomials of higher and fractional order to operational techniques. We show that these methods are useful to obtain either numerical or analytical solutions.
From Newton's Equation to Fractional Diffusion and Wave Equations
Vázquez Luis
2011-01-01
Full Text Available Fractional calculus represents a natural instrument to model nonlocal (or long-range dependence phenomena either in space or time. The processes that involve different space and time scales appear in a wide range of contexts, from physics and chemistry to biology and engineering. In many of these problems, the dynamics of the system can be formulated in terms of fractional differential equations which include the nonlocal effects either in space or time. We give a brief, nonexhaustive, panoramic view of the mathematical tools associated with fractional calculus as well as a description of some fields where either it is applied or could be potentially applied.
The numerical simulation of convection delayed dominated diffusion equation
Mohan Kumar P. Murali
2016-01-01
Full Text Available In this paper, we propose a fitted numerical method for solving convection delayed dominated diffusion equation. A fitting factor is introduced and the model equation is discretized by cubic spline method. The error analysis is analyzed for the consider problem. The numerical examples are solved using the present method and compared the result with the exact solution.
EXPONENTIAL ATTRACTOR FOR A CLASS OF NONCLASSICAL DIFFUSION EQUATION
尚亚东; 郭柏灵
2003-01-01
In this paper,we consider the asymptotic behavior of solutions for a class of nonclassical diffusion equation.We show the squeezing property and the existence of exponential attractor for this equation.We also make the estimates on its fractal dimension and exponential attraction.
Semianalytic Solution of Space-Time Fractional Diffusion Equation
A. Elsaid
2016-01-01
Full Text Available We study the space-time fractional diffusion equation with spatial Riesz-Feller fractional derivative and Caputo fractional time derivative. The continuation of the solution of this fractional equation to the solution of the corresponding integer order equation is proved. The series solution of this problem is obtained via the optimal homotopy analysis method (OHAM. Numerical simulations are presented to validate the method and to show the effect of changing the fractional derivative parameters on the solution behavior.
Exact solutions for logistic reaction-diffusion equations in biology
Broadbridge, P.; Bradshaw-Hajek, B. H.
2016-08-01
Reaction-diffusion equations with a nonlinear source have been widely used to model various systems, with particular application to biology. Here, we provide a solution technique for these types of equations in N-dimensions. The nonclassical symmetry method leads to a single relationship between the nonlinear diffusion coefficient and the nonlinear reaction term; the subsequent solutions for the Kirchhoff variable are exponential in time (either growth or decay) and satisfy the linear Helmholtz equation in space. Example solutions are given in two dimensions for particular parameter sets for both quadratic and cubic reaction terms.
Linear fractional diffusion-wave equation for scientists and engineers
Povstenko, Yuriy
2015-01-01
This book systematically presents solutions to the linear time-fractional diffusion-wave equation. It introduces the integral transform technique and discusses the properties of the Mittag-Leffler, Wright, and Mainardi functions that appear in the solutions. The time-nonlocal dependence between the flux and the gradient of the transported quantity with the “long-tail” power kernel results in the time-fractional diffusion-wave equation with the Caputo fractional derivative. Time-nonlocal generalizations of classical Fourier’s, Fick’s and Darcy’s laws are considered and different kinds of boundary conditions for this equation are discussed (Dirichlet, Neumann, Robin, perfect contact). The book provides solutions to the fractional diffusion-wave equation with one, two and three space variables in Cartesian, cylindrical and spherical coordinates. The respective sections of the book can be used for university courses on fractional calculus, heat and mass transfer, transport processes in porous media and ...
Lattice Boltzmann method for the fractional advection-diffusion equation
Zhou, J. G.; Haygarth, P. M.; Withers, P. J. A.; Macleod, C. J. A.; Falloon, P. D.; Beven, K. J.; Ockenden, M. C.; Forber, K. J.; Hollaway, M. J.; Evans, R.; Collins, A. L.; Hiscock, K. M.; Wearing, C.; Kahana, R.; Villamizar Velez, M. L.
2016-04-01
Mass transport, such as movement of phosphorus in soils and solutes in rivers, is a natural phenomenon and its study plays an important role in science and engineering. It is found that there are numerous practical diffusion phenomena that do not obey the classical advection-diffusion equation (ADE). Such diffusion is called abnormal or superdiffusion, and it is well described using a fractional advection-diffusion equation (FADE). The FADE finds a wide range of applications in various areas with great potential for studying complex mass transport in real hydrological systems. However, solution to the FADE is difficult, and the existing numerical methods are complicated and inefficient. In this study, a fresh lattice Boltzmann method is developed for solving the fractional advection-diffusion equation (LabFADE). The FADE is transformed into an equation similar to an advection-diffusion equation and solved using the lattice Boltzmann method. The LabFADE has all the advantages of the conventional lattice Boltzmann method and avoids a complex solution procedure, unlike other existing numerical methods. The method has been validated through simulations of several benchmark tests: a point-source diffusion, a boundary-value problem of steady diffusion, and an initial-boundary-value problem of unsteady diffusion with the coexistence of source and sink terms. In addition, by including the effects of the skewness β , the fractional order α , and the single relaxation time τ , the accuracy and convergence of the method have been assessed. The numerical predictions are compared with the analytical solutions, and they indicate that the method is second-order accurate. The method presented will allow the FADE to be more widely applied to complex mass transport problems in science and engineering.
Nikitin, A. G.
2004-01-01
Group classification of systems of two coupled nonlinear reaction-diffusion equation with a diagonal diffusion matrix is carried out. Symmetries of diffusion systems with singular diffusion matrix and additional first order derivative terms are described.
Stochastic differential equations and diffusion processes
Ikeda, N
1989-01-01
Being a systematic treatment of the modern theory of stochastic integrals and stochastic differential equations, the theory is developed within the martingale framework, which was developed by J.L. Doob and which plays an indispensable role in the modern theory of stochastic analysis.A considerable number of corrections and improvements have been made for the second edition of this classic work. In particular, major and substantial changes are in Chapter III and Chapter V where the sections treating excursions of Brownian Motion and the Malliavin Calculus have been expanded and refined. Sectio
Nonlinear delay reaction-diffusion equations: exact traveling wave solutions
We consider nonlinear delay reaction-diffusion equations. We present a number of exact traveling wave solutions, which can be represented in terms of elementary functions. We consider equations with quadratic, power, exponential, and logarithmic nonlinearities, as well as more complex equations with the kinetic function depending on one or two arbitrary functions of a single argument. All of the solutions obtained involve several free parameters and can be used to test approximate analytical and numerical methods for solving similar and more complex nonlinear differential-difference equations
Ho, C. -L.; Lee, C.-C.
2015-01-01
We consider solvability of the generalized reaction-diffusion equation with both space- and time-dependent diffusion and reaction terms by means of the similarity method. By introducing the similarity variable, the reaction-diffusion equation is reduced to an ordinary differential equation. Matching the resulting ordinary differential equation with known exactly solvable equations, one can obtain corresponding exactly solvable reaction-diffusion systems. Several representative examples of exa...
Unconditionally stable diffusion-acceleration of the transport equation
The standard iterative procedure for solving fixed-source discrete-ordinates problems converges very slowly for problems in optically large regions with scattering ratios c near unity. The diffusion-synthetic acceleration method has been proposed to make use of the fact that for this class of problems the diffusion equation is often an accurate approximation to the transport equation. However, stability difficulties have historically hampered the implementation of this method for general transport differencing schemes. In this article we discuss a recently developed procedure for obtaining unconditionally stable diffusion-synthetic acceleration methods for various transport differencing schemes. We motivate the analysis by first discussing the exact transport equation; then we illustrate the procedure by deriving a new stable acceleration method for the linear discontinuous transport differencing scheme. We also provide some numerical results
Nonlocalized modulation of periodic reaction diffusion waves: The Whitham equation
Johnson, Mathew A.; Noble, Pascal; Rodrigues, L. Miguel; Zumbrun, Kevin
2013-02-01
In a companion paper, we established nonlinear stability with detailed diffusive rates of decay of spectrally stable periodic traveling-wave solutions of reaction diffusion systems under small perturbations consisting of a nonlocalized modulation plus a localized ( L 1) perturbation. Here, we determine time-asymptotic behavior under such perturbations, showing that solutions consist of a leading order of a modulation whose parameter evolution is governed by an associated Whitham averaged equation.
Nonlocalized modulation of periodic reaction diffusion waves: The Whitham equation
Johnson, Mathew; Rodrigues, L Miguel; Zumbrun, Kevin
2011-01-01
In a companion paper, we established nonlinear stability with detailed diffusive rates of decay of spectrally stable periodic traveling-wave solutions of reaction diffusion systems under small perturbations consisting of a nonlocalized modulation plus a localized perturbation. Here, we determine time-asymptotic behavior under such perturbations, showing that solutions consist to leading order of a modulation whose parameter evolution is governed by an associated Whitham averaged equation.
Reaction-diffusion equations with spatially distributed hysteresis
Gurevich, Pavel; Shamin, Roman; Tikhomirov, Sergey
2012-01-01
The paper deals with reaction-diffusion equations involving a hysteretic discontinuity in the source term, which is defined at each spatial point. In particular, such problems describe chemical reactions and biological processes in which diffusive and nondiffusive substances interact according to hysteresis law. We find sufficient conditions that guarantee the existence and uniqueness of solutions as well as their continuous dependence on initial data.
Traveling waves for a boundary reaction-diffusion equation
Caffarelli, L; Sire, Y
2011-01-01
We prove the existence of a traveling wave solution for a boundary reaction diffusion equation when the reaction term is the combustion nonlinearity with ignition temperature. A key role in the proof is plaid by an explicit formula for traveling wave solutions of a free boundary problem obtained as singular limit for the reaction-diffusion equation (the so-called high energy activation energy limit). This explicit formula, which is interesting in itself, also allows us to get an estimate on the decay at infinity of the traveling wave (which turns out to be faster than the usual exponential decay).
Phasefield theory for fractional diffusion-reaction equations and applications
Imbert, Cyril
2009-01-01
This paper is concerned with diffusion-reaction equations where the classical diffusion term, such as the Laplacian operator, is replaced with a singular integral term, such as the fractional Laplacian operator. As far as the reaction term is concerned, we consider bistable non-linearities. After properly rescaling (in time and space) these integro-differential evolution equations, we show that the limits of their solutions as the scaling parameter goes to zero exhibit interfaces moving by anisotropic mean curvature. The singularity and the unbounded support of the potential at stake are both the novelty and the challenging difficulty of this work.
Continuity, the Bloch-Torrey equation, and Diffusion MRI
Hall, Matt G
2016-01-01
The Bloch equation describes the evolution of classical particles tagged with a magnetisation vector in a strong magnetic field and is fundamental to many NMR and MRI contrast methods. The equation can be generalised to include the effects of spin motion by including a spin flux, which typically contains a Fickian diffusive term and/or a coherent velocity term. This form is known as the Bloch-Torrey equation, and is fundamental to MR modalities which are sensitive to spin dynamics such as diffusion MRI. Such modalities have received a great deal of interest in the research literature over the last few years, resulting in a huge range of models and methods. In this work we make make use of a more general Bloch-Torrey equation with a generalised flux term. We show that many commonly employed approaches in Diffusion MRI may be viewed as different choices for the flux terms in this equation. This viewpoint, although obvious theoretically, is not usually emphasised in the diffusion MR literature and points to inte...
The Nonclassical Diffusion Approximation to the Nonclassical Linear Boltzmann Equation
Vasques, Richard
2015-01-01
We show that, by correctly selecting the probability distribution function $p(s)$ for a particle's distance-to-collision, the nonclassical diffusion equation can be represented exactly by the nonclassical linear Boltzmann equation for an infinite homogeneous medium. This choice of $p(s)$ preserves the $true$ mean-squared free path of the system, which sheds new light on the results obtained in previous work.
Nikitin, A. G.
2004-01-01
Group classification of the generalized complex Ginzburg-Landau equations is presented. An approach to group classification of systems of reaction-diffusion equations with general diffusion matrix is developed.
COUPLED CHEMOTAXIS FLUID MODEL
LORZ, ALEXANDER
2010-06-01
We consider a model system for the collective behavior of oxygen-driven swimming bacteria in an aquatic fluid. In certain parameter regimes, such suspensions of bacteria feature large-scale convection patterns as a result of the hydrodynamic interaction between bacteria. The presented model consist of a parabolicparabolic chemotaxis system for the oxygen concentration and the bacteria density coupled to an incompressible Stokes equation for the fluid driven by a gravitational force of the heavier bacteria. We show local existence of weak solutions in a bounded domain in d, d = 2, 3 with no-flux boundary condition and in 2 in the case of inhomogeneous Dirichlet conditions for the oxygen. © 2010 World Scientific Publishing Company.
High order backward discretization of the neutron diffusion equation
Fast codes capable of dealing with three-dimensional geometries, are needed to be able to simulate spatially complicated transients in a nuclear reactor. We propose a new discretization technique for the time integration of the neutron diffusion equation, based on the backward difference formulas for systems of stiff ordinary differential equations. This method needs to solve a system of linear equations for each integration step, and for this purpose, we have developed an iterative block algorithm combined with a variational acceleration technique. We tested the algorithm with two benchmark problems, and compared the results with those provided by other codes, concluding that the performance and overall agreement are very good. (author)
Kinetic methods for the diffusion equations in nuclear PWR
The neutronic diffusion equations are used to calculate the flux in the core of a nuclear reactor. We can solve these equations with direct methods, but these methods are too expensive for slow transient phenomenons. So, according to a method developed by Ott and Meneley, we split the flux into a product of two functions: the first one is only time dependent and is supposed to change rapidly. The second one is time and space dependent but is supposed to vary slowly with time. Numerical results are given and compared to those obtained by a direct resolution of the equations and by experiments
THE INTERIOR LAYER SOLUTION TO NONLOCAL REACTION DIFFUSION EQUATIONS
无
2011-01-01
An initial boundary value problem of semilinear nonlocal reaction diffusion equations is considered.Under some suitable conditions,using the asymptotic theory,the existence and asymptotic behavior of the interior layer solution to the initial boundary value problem are studied.
Target Pattern Waves in Specific Reaction-Diffusion Equation
ZHOU Tian-Shou; TANG Yun
2002-01-01
This paper uses two-timing to solve a class of reaction-diffusion equations of Oregonator and constructsformally stable target pattern solutions to leading term in small parameter, which can be used to account for someobserved experimental facts. In addition, it also gives the explicit expression of the period corresponding to the targetwaves.
Gradient estimates for a nonlinear diffusion equation on complete manifolds
Wu, Jia-Yong
2010-01-01
Let $(M,g)$ be a complete non-compact Riemannian manifold with the $m$-dimensional Bakry-\\'{E}mery Ricci curvature bounded below by a non-positive constant. In this paper, we give a localized Hamilton-type gradient estimate for the positive smooth bounded solutions to the following nonlinear diffusion equation \\[ u_t=\\Delta u-\
Diffusive Limits of the Master Equation in Inhomogeneous Media
Sattin, F; Salasnich, L
2015-01-01
In inhomogeneous environments several expressions for the flux of a diffusing quantity may apply--from Fick-Fourier's to Fokker-Planck's--depending upon the system studied. The integro-differential Master Equation (ME) provides a fairly generic framework for describing the dynamics of arbitrary systems driven by stochastic rules. Diffusive dynamics does arise as long-wavelength limit of the ME. However, while it is straightforward to obtain a diffusion equation with Fokker-Planck flux, its Fick-Fourier counterpart has never been worked out from the ME. In this work we show under which hypothesis the Fick's flux can actually be recovered from the ME. Analytical considerations are supported by explicit computer models.
Numerical solution of a reaction-diffusion equation
The purpose of the present work to continue the observations and the numerical experiences on a reaction-diffusion model, that is a simplified form of the neutronic flux equation. The model is parabolic, nonlinear, with Dirichlet boundary conditions. The purpose is to approximate non trivial solutions, asymptotically stables for t → ∞, that is solutions that tend to the elliptic problem, in the Lyapunov sense. It belongs to the so-called reaction-diffusion equations of semi linear kind, that is, linear equations in the heat operator and they have a nonlinear reaction function, in this case f (u, a, b) = u (a - b u), being u concentration, a and b parameters. The study of the incidence of these parameters take an interest to the neutronic flux physics. So that we search non trivial, positive and bounded solutions. The used algorithm is based on the concept of monotone and ordered sequences, and on the existence theorem of Amann and Sattinger. (author)
Kinetic equations for diffusion in the presence of entropic barriers.
Reguera, D; Rubí, J M
2001-12-01
We use the mesoscopic nonequilibrium thermodynamics theory to derive the general kinetic equation of a system in the presence of potential barriers. The result is applied to a description of the evolution of systems whose dynamics is influenced by entropic barriers. We analyze in detail the case of diffusion in a domain of irregular geometry in which the presence of the boundaries induces an entropy barrier when approaching the exact dynamics by a coarsening of the description. The corresponding kinetic equation, named the Fick-Jacobs equation, is obtained, and its validity is generalized through the formulation of a scaling law for the diffusion coefficient which depends on the shape of the boundaries. The method we propose can be useful to analyze the dynamics of systems at the nanoscale where the presence of entropy barriers is a common feature. PMID:11736170
Langevin equation with fluctuating diffusivity: A two-state model.
Miyaguchi, Tomoshige; Akimoto, Takuma; Yamamoto, Eiji
2016-07-01
Recently, anomalous subdiffusion, aging, and scatter of the diffusion coefficient have been reported in many single-particle-tracking experiments, though the origins of these behaviors are still elusive. Here, as a model to describe such phenomena, we investigate a Langevin equation with diffusivity fluctuating between a fast and a slow state. Namely, the diffusivity follows a dichotomous stochastic process. We assume that the sojourn time distributions of these two states are given by power laws. It is shown that, for a nonequilibrium ensemble, the ensemble-averaged mean-square displacement (MSD) shows transient subdiffusion. In contrast, the time-averaged MSD shows normal diffusion, but an effective diffusion coefficient transiently shows aging behavior. The propagator is non-Gaussian for short time and converges to a Gaussian distribution in a long-time limit; this convergence to Gaussian is extremely slow for some parameter values. For equilibrium ensembles, both ensemble-averaged and time-averaged MSDs show only normal diffusion and thus we cannot detect any traces of the fluctuating diffusivity with these MSDs. Therefore, as an alternative approach to characterizing the fluctuating diffusivity, the relative standard deviation (RSD) of the time-averaged MSD is utilized and it is shown that the RSD exhibits slow relaxation as a signature of the long-time correlation in the fluctuating diffusivity. Furthermore, it is shown that the RSD is related to a non-Gaussian parameter of the propagator. To obtain these theoretical results, we develop a two-state renewal theory as an analytical tool. PMID:27575079
Langevin equation with fluctuating diffusivity: A two-state model
Miyaguchi, Tomoshige; Akimoto, Takuma; Yamamoto, Eiji
2016-07-01
Recently, anomalous subdiffusion, aging, and scatter of the diffusion coefficient have been reported in many single-particle-tracking experiments, though the origins of these behaviors are still elusive. Here, as a model to describe such phenomena, we investigate a Langevin equation with diffusivity fluctuating between a fast and a slow state. Namely, the diffusivity follows a dichotomous stochastic process. We assume that the sojourn time distributions of these two states are given by power laws. It is shown that, for a nonequilibrium ensemble, the ensemble-averaged mean-square displacement (MSD) shows transient subdiffusion. In contrast, the time-averaged MSD shows normal diffusion, but an effective diffusion coefficient transiently shows aging behavior. The propagator is non-Gaussian for short time and converges to a Gaussian distribution in a long-time limit; this convergence to Gaussian is extremely slow for some parameter values. For equilibrium ensembles, both ensemble-averaged and time-averaged MSDs show only normal diffusion and thus we cannot detect any traces of the fluctuating diffusivity with these MSDs. Therefore, as an alternative approach to characterizing the fluctuating diffusivity, the relative standard deviation (RSD) of the time-averaged MSD is utilized and it is shown that the RSD exhibits slow relaxation as a signature of the long-time correlation in the fluctuating diffusivity. Furthermore, it is shown that the RSD is related to a non-Gaussian parameter of the propagator. To obtain these theoretical results, we develop a two-state renewal theory as an analytical tool.
On analytical solutions for the nonlinear diffusion equation
Ulrich Olivier Dangui-Mbani
2014-09-01
Full Text Available The nonlinear diffusion equation arises in many important areas of nonlinear problems of heat and mass transfer, biological systems and processes involving fluid flow and most of the known exact solutions turn out to be approximate solutions in the form of a series which is the exact solution in the closed form. The approximate results obtained by using Homotopy perturbation transform method (HPTM and have been compared with the exact solutions by using software “mathematica” to show the stability of the solutions of nonlinear equation. The comparisons indicate that there is a very good agreement between the HPTM solutions and exact solutions in terms of accuracy
Simple jumping process with memory: Transport equation and diffusion
Kamińska, A.; Srokowski, T.
2004-06-01
We present a stochastic jumping process, defined in terms of jump-size probability density and jumping rate, which is a generalization of the well-known kangaroo process. The definition takes into account two process values: after and before the jump. Therefore, the process is able to preserve memory about its previous values. It possesses a simple stationary limit. Its master equation is interpreted as the kinetic equation with variable collision rate. The process can be easily applied to model systems which relax to distributions other than Maxwellian. The case of a constant jumping rate corresponds to the diffusion process, either normal or ballistic.
Reaction diffusion equation with spatio-temporal delay
Zhao, Zhihong; Rong, Erhua
2014-07-01
We investigate reaction-diffusion equation with spatio-temporal delays, the global existence, uniqueness and asymptotic behavior of solutions for which in relation to constant steady-state solution, included in the region of attraction of a stable steady solution. It is shown that if the delay reaction function satisfies some conditions and the system possesses a pair of upper and lower solutions then there exists a unique global solution. In terms of the maximal and minimal constant solutions of the corresponding steady-state problem, we get the asymptotic stability of reaction-diffusion equation with spatio-temporal delay. Applying this theory to Lotka-Volterra model with spatio-temporal delay, we get the global solution asymptotically tend to the steady-state problem's steady-state solution.
Algorithm refinement for stochastic partial differential equations I. linear diffusion
Alexander, F J; Tartakovsky, D M
2002-01-01
A hybrid particle/continuum algorithm is formulated for Fickian diffusion in the fluctuating hydrodynamic limit. The particles are taken as independent random walkers; the fluctuating diffusion equation is solved by finite differences with deterministic and white-noise fluxes. At the interface between the particle and continuum computations the coupling is by flux matching, giving exact mass conservation. This methodology is an extension of Adaptive Mesh and Algorithm Refinement to stochastic partial differential equations. Results from a variety of numerical experiments are presented for both steady and time-dependent scenarios. In all cases the mean and variance of density are captured correctly by the stochastic hybrid algorithm. For a nonstochastic version (i.e., using only deterministic continuum fluxes) the mean density is correct, but the variance is reduced except in particle regions away from the interface. Extensions of the methodology to fluid mechanics applications are discussed.
Convergence of the iteration to solve the diffusion equation
Linear boundary value problems of diffusion theory in a symmetric region are analyzed with the help of group theory. If the symmetries of the region commute with the operator of the equation to be solved, the symmetries generate a classification of the solution space. This technique is applied to study the numerical behavior of numerical methods. It is shown that approximate polynomials on the boundary and inside the region must be compatible to get stable results. Compatibility conditions are given. To avoid convergence problems, linear functions on the six faces of a regular hexagon would demand as high as sixth order approximation inside the hexagon. The predicted convergence problems have been observed in the neutron physics code VARIANT, which is a routinely used diffusion equation solver at ANL. (Author)
Support Operators Method for the Diffusion Equation in Multiple Materials
Winters, Andrew R. [Los Alamos National Laboratory; Shashkov, Mikhail J. [Los Alamos National Laboratory
2012-08-14
A second-order finite difference scheme for the solution of the diffusion equation on non-uniform meshes is implemented. The method allows the heat conductivity to be discontinuous. The algorithm is formulated on a one dimensional mesh and is derived using the support operators method. A key component of the derivation is that the discrete analog of the flux operator is constructed to be the negative adjoint of the discrete divergence, in an inner product that is a discrete analog of the continuum inner product. The resultant discrete operators in the fully discretized diffusion equation are symmetric and positive definite. The algorithm is generalized to operate on meshes with cells which have mixed material properties. A mechanism to recover intermediate temperature values in mixed cells using a limited linear reconstruction is introduced. The implementation of the algorithm is verified and the linear reconstruction mechanism is compared to previous results for obtaining new material temperatures.
Wavelets method for the time fractional diffusion-wave equation
Heydari, M.H., E-mail: heydari@stu.yazd.ac.ir [Faculty of Mathematics, Yazd University, Yazd (Iran, Islamic Republic of); The Laboratory of Quantum Information Processing, Yazd University, Yazd (Iran, Islamic Republic of); Hooshmandasl, M.R., E-mail: hooshmandasl@yazd.ac.ir [Faculty of Mathematics, Yazd University, Yazd (Iran, Islamic Republic of); The Laboratory of Quantum Information Processing, Yazd University, Yazd (Iran, Islamic Republic of); Maalek Ghaini, F.M., E-mail: maalek@yazd.ac.ir [Faculty of Mathematics, Yazd University, Yazd (Iran, Islamic Republic of); The Laboratory of Quantum Information Processing, Yazd University, Yazd (Iran, Islamic Republic of); Cattani, C., E-mail: ccattani@unisa.it [Department of Mathematics, University of Salerno, Fisciano (Italy)
2015-01-23
In this paper, an efficient and accurate computational method based on the Legendre wavelets (LWs) is proposed for solving the time fractional diffusion-wave equation (FDWE). To this end, a new fractional operational matrix (FOM) of integration for the LWs is derived. The LWs and their FOM of integration are used to transform the problem under consideration into a linear system of algebraic equations, which can be simply solved to achieve the solution of the problem. The proposed method is very convenient for solving such problems, since the initial and boundary conditions are taken into account automatically. - Highlights: • A new operational matrix of fractional integration for the LWs is derived. • A new method based on the LWs is proposed for the time FDWE. • The paper contains some useful properties of the LWs. • The proposed method can be applied for fractional sub-diffusion systems. • The proposed method can be extended for fourth-order FDWE.
UPWIND SPLITTING SCHEME FOR CONVECTION-DIFFUSION EQUATIONS
无
2000-01-01
A new class of numerical schemes is proposed to solve convection-diffusion equations by combining the upwind technique and the method of operator splitting. For every time step, the multi-dimensional approximation is performed in several independent directions alternatively, while the upwind technique is applied to treat the convection term in every individual direction. This scheme possesses maximum principle. Stability and convergence are analysed by energy method.
HEXAN - a hexagonal nodal code for solving the diffusion equation
This report describes the theory of and provides a user's manual for the HEXAN program, which is a nodal program for the solution of the few-group diffusion equation in hexagonal geometry. Based upon symmetry considerations, the theory provides an analytical solution in a homogeneous node. WWER and HTGR test problem solutions are presented. The equivalence of the finite-difference scheme and the response matrix method is proven. The properties of a symmetric node's response matrix are investigated. (author)
CONVERGENCE OF MARKOV CHAIN APPROXIMATIONS TO STOCHASTIC REACTION DIFFUSION EQUATIONS
Kouritzin, Michael A.; Hongwei Long
2001-01-01
In the context of simulating the transport of a chemical or bacterial contaminant through a moving sheet of water, we extend a well-established method of approximating reaction-diffusion equations with Markov chains by allowing convection, certain Poisson measure driving sources and a larger class of reaction functions. Our alterations also feature dramatically slower Markov chain state change rates often yielding a ten to one-hundred-fold simulation speed increase over the previous version o...
The three-dimensional generalized (self-adjoint) atmospheric diffusion equation can be solved analytically. Where the source strength is located at any where within the region of interest, and with arbitrary power- law function of wind speed and eddy diffusivities. The Greens method approach is utilized, where the Green's function is splitting into two factors, As a result of the hermeticity the solution of the diffusion equation can be expressed exact, as a multiple of source and the two factors of the Green's function, vertical dispersion factor, and cross wind dispersion factor. The two factors of the Green's function are derived analytically in which, the Green's of the diffusion equation is splitting into two pair of two-dimensional differential equations, by using: separation of variables, Bessel equation, and transform the dependent and independent variables, and the cosine transformation. All amenable cases are adopted such as: Neumann (total reflection), Dirichlet (total adsorption), and mixed boundary conditions. Gaussian plume model is extracted in all the previously cases
Reaction rates for a generalized reaction-diffusion master equation
Hellander, Stefan; Petzold, Linda
2016-01-01
It has been established that there is an inherent limit to the accuracy of the reaction-diffusion master equation. Specifically, there exists a fundamental lower bound on the mesh size, below which the accuracy deteriorates as the mesh is refined further. In this paper we extend the standard reaction-diffusion master equation to allow molecules occupying neighboring voxels to react, in contrast to the traditional approach, in which molecules react only when occupying the same voxel. We derive reaction rates, in two dimensions as well as three dimensions, to obtain an optimal match to the more fine-grained Smoluchowski model and show in two numerical examples that the extended algorithm is accurate for a wide range of mesh sizes, allowing us to simulate systems that are intractable with the standard reaction-diffusion master equation. In addition, we show that for mesh sizes above the fundamental lower limit of the standard algorithm, the generalized algorithm reduces to the standard algorithm. We derive a lower limit for the generalized algorithm which, in both two dimensions and three dimensions, is of the order of the reaction radius of a reacting pair of molecules.
Analysis of a mixed space-time diffusion equation
Momoniat, Ebrahim
2015-06-01
An energy method is used to analyze the stability of solutions of a mixed space-time diffusion equation that has application in the unidirectional flow of a second-grade fluid and the distribution of a compound Poisson process. Solutions to the model equation satisfying Dirichlet boundary conditions are proven to dissipate total energy and are hence stable. The stability of asymptotic solutions satisfying Neumann boundary conditions coincides with the condition for the positivity of numerical solutions of the model equation from a Crank-Nicolson scheme. The Crank-Nicolson scheme is proven to yield stable numerical solutions for both Dirichlet and Neumann boundary conditions for positive values of the critical parameter. Numerical solutions are compared to analytical solutions that are valid on a finite domain.
A mixed finite element method for nonlinear diffusion equations
Burger, Martin
2010-01-01
We propose a mixed finite element method for a class of nonlinear diffusion equations, which is based on their interpretation as gradient flows in optimal transportation metrics. We introduce an appropriate linearization of the optimal transport problem, which leads to a mixed symmetric formulation. This formulation preserves the maximum principle in case of the semi-discrete scheme as well as the fully discrete scheme for a certain class of problems. In addition solutions of the mixed formulation maintain exponential convergence in the relative entropy towards the steady state in case of a nonlinear Fokker-Planck equation with uniformly convex potential. We demonstrate the behavior of the proposed scheme with 2D simulations of the porous medium equations and blow-up questions in the Patlak-Keller-Segel model. © American Institute of Mathematical Sciences.
A Nonlinear Singularly Perturbed Problem for Reaction Diffusion Equations with Boundary Perturbation
Jia-qi Mo; Wan-tao Lin
2005-01-01
A nonlinear singularly perturbed problems for reaction diffusion equation with boundary perturbation is considered. Under suitable conditions, the asymptotic behavior of solution for the initial boundary value problems of reaction diffusion equations is studied using the theory of differential inequalities.
Group Foliation Method and Functional Separation of Variables to Nonlinear Diffusion Equations
QU Chang-Zheng; ZHANG Shun-Li
2005-01-01
@@ Generalized functional separation of variables to nonlinear diffusion equations is studied in terms of the extended group foliation method. A complete classification for the nonlinear diffusion equation with source term which admits functional separable solutions is presented.
Interior corner problem for the neutron diffusion equation
The qualitative and quantitative behavior of the solution of the neutron diffusion equation at interior corners of two different materials is investigated. This is done by assuming regionwise constant coefficients and sources, reducing the problem to dimensionless form and solving the Neumann problem over a circular disk. The general analytic solution is given, as well as special solutions. A series of parameter studies have been made which include the extreme cases that can occur in neutron diffusion. The results including eigenvalues, eigenfunctions, Fourier coefficients and flux profiles are tabulated and presented in graphical form. It is shown how to select appropriate function spaces to obtain good numerical methods for coarse network calculations. (33 figures, 25 tables) (U.S.)
THE NONLINEAR NONLOCAL SINGULARLY PERTURBED PROBLEMS FOR REACTION DIFFUSION EQUATIONS
莫嘉琪; 朱江
2003-01-01
A class of nonlinear nonlocal for singularly perturbed Robin initial boundaryvalue problems for reaction diffusion equations is considered. Under suitable conditions,firstly, the outer solution of the original problem is obtained, secondly, using the stretchedvariable, the composing expansion method and the expanding theory of power series theinitial layer is constructed, finally, using the theory of differential inequalities theasymptotic behavior of solution for the initial boundary value problems are studied andeducing some relational inequalities the existence and uniqueness of solution for the originalproblem and the uniformly valid asymptotic estimation is discussed.
Diffusive Wave Approximation to the Shallow Water Equations: Computational Approach
Collier, Nathan
2011-05-14
We discuss the use of time adaptivity applied to the one dimensional diffusive wave approximation to the shallow water equations. A simple and computationally economical error estimator is discussed which enables time-step size adaptivity. This robust adaptive time discretization corrects the initial time step size to achieve a user specified bound on the discretization error and allows time step size variations of several orders of magnitude. In particular, in the one dimensional results presented in this work feature a change of four orders of magnitudes for the time step over the entire simulation.
Dynamic hysteresis modeling including skin effect using diffusion equation model
Hamada, Souad; Louai, Fatima Zohra; Nait-Said, Nasreddine; Benabou, Abdelkader
2016-07-01
An improved dynamic hysteresis model is proposed for the prediction of hysteresis loop of electrical steel up to mean frequencies, taking into account the skin effect. In previous works, the analytical solution of the diffusion equation for low frequency (DELF) was coupled with the inverse static Jiles-Atherton (JA) model in order to represent the hysteresis behavior for a lamination. In the present paper, this approach is improved to ensure the reproducibility of measured hysteresis loops at mean frequency. The results of simulation are compared with the experimental ones. The selected results for frequencies 50 Hz, 100 Hz, 200 Hz and 400 Hz are presented and discussed.
IDENTIFYING AN UNKNOWN SOURCE IN SPACE-FRACTIONAL DIFFUSION EQUATION
杨帆; 傅初黎; 李晓晓
2014-01-01
In this paper, we identify a space-dependent source for a fractional diffusion equation. This problem is ill-posed, i.e., the solution (if it exists) does not depend continu-ously on the data. The generalized Tikhonov regularization method is proposed to solve this problem. An a priori error estimate between the exact solution and its regularized approxi-mation is obtained. Moreover, an a posteriori parameter choice rule is proposed and a stable error estimate is also obtained. Numerical examples are presented to illustrate the validity and effectiveness of this method.
Multigrid solution of diffusion equations on distributed memory multiprocessor systems
The subject is the solution of partial differential equations for simulation of the reactor core on high-performance computers. The parallelization and implementation of nodal multigrid diffusion algorithms on array and ring configurations of the DIRMU multiprocessor system is outlined. The particular iteration scheme employed in the nodal expansion method appears similarly efficient in serial and parallel environments. The combination of modern multi-level techniques with innovative hardware (vector-multiprocessor systems) provides powerful tools needed for real time simulation of physical systems. The parallel efficiencies range from 70 to 90%. The same performance is estimated for large problems on large multiprocessor systems being designed at present. (orig.)
Mimetic discretization of two-dimensional magnetic diffusion equations
Lipnikov, Konstantin; Reynolds, James; Nelson, Eric
2013-08-01
In case of non-constant resistivity, cylindrical coordinates, and highly distorted polygonal meshes, a consistent discretization of the magnetic diffusion equations requires new discretization tools based on a discrete vector and tensor calculus. We developed a new discretization method using the mimetic finite difference framework. It is second-order accurate on arbitrary polygonal meshes and a consistent calculation of the Joule heating is intrinsic within it. The second-order convergence rates in L2 and L1 norms were verified with numerical experiments.
New Parallel Three-level Iterative Method for Diffusion Equation
Tinghuai Ma
2010-01-01
Full Text Available To solve the diffusion equation on parallel computers, we first derived an o(τ2+h6 order implicit finite difference method based on a class of alternating group explicit iterative method. Based on this method, we devised a new alternating group explicit iterative method. Moreover, the absolute stability and convergence of the New Alternating Group Explicit Iterative (N-AGEI method was proved. Finally, the numerical experiments were conducted to verify our method. Both the theoretical analysis and simulation results showed that our proposed difference format had satisfied stability error estimate and convergence.
Poincar\\'e's inequality and diffusive evolution equations
Bjorland, Clayton; Schonbek, Maria E.
2007-01-01
This paper addresses the question of change of decay rate from exponential to algebraic for diffusive evolution equations. We show how the behaviour of the spectrum of the Dirichlet Laplacian in the two cases yields the passage from exponential decay in bounded domains to algebraic decay or no decay at all in the case of unbounded domains. It is well known that such rates of decay exist: the purpose of this paper is to explain what makes the change in decay happen. We also discuss what kind o...
Chaotic dynamics and diffusion in a piecewise linear equation
Shahrear, Pabel, E-mail: pabelshahrear@yahoo.com [Department of Mathematics, Shah Jalal University of Science and Technology, Sylhet–3114 (Bangladesh); Glass, Leon, E-mail: glass@cnd.mcgill.ca [Department of Physiology, 3655 Promenade Sir William Osler, McGill University, Montreal, Quebec H3G 1Y6 (Canada); Edwards, Rod, E-mail: edwards@uvic.ca [Department of Mathematics and Statistics, University of Victoria, P.O. Box 1700 STN CSC, Victoria, British Columbia V8W 2Y2 (Canada)
2015-03-15
Genetic interactions are often modeled by logical networks in which time is discrete and all gene activity states update simultaneously. However, there is no synchronizing clock in organisms. An alternative model assumes that the logical network is preserved and plays a key role in driving the dynamics in piecewise nonlinear differential equations. We examine dynamics in a particular 4-dimensional equation of this class. In the equation, two of the variables form a negative feedback loop that drives a second negative feedback loop. By modifying the original equations by eliminating exponential decay, we generate a modified system that is amenable to detailed analysis. In the modified system, we can determine in detail the Poincaré (return) map on a cross section to the flow. By analyzing the eigenvalues of the map for the different trajectories, we are able to show that except for a set of measure 0, the flow must necessarily have an eigenvalue greater than 1 and hence there is sensitive dependence on initial conditions. Further, there is an irregular oscillation whose amplitude is described by a diffusive process that is well-modeled by the Irwin-Hall distribution. There is a large class of other piecewise-linear networks that might be analyzed using similar methods. The analysis gives insight into possible origins of chaotic dynamics in periodically forced dynamical systems.
Chaotic dynamics and diffusion in a piecewise linear equation
Shahrear, Pabel; Glass, Leon; Edwards, Rod
2015-03-01
Genetic interactions are often modeled by logical networks in which time is discrete and all gene activity states update simultaneously. However, there is no synchronizing clock in organisms. An alternative model assumes that the logical network is preserved and plays a key role in driving the dynamics in piecewise nonlinear differential equations. We examine dynamics in a particular 4-dimensional equation of this class. In the equation, two of the variables form a negative feedback loop that drives a second negative feedback loop. By modifying the original equations by eliminating exponential decay, we generate a modified system that is amenable to detailed analysis. In the modified system, we can determine in detail the Poincaré (return) map on a cross section to the flow. By analyzing the eigenvalues of the map for the different trajectories, we are able to show that except for a set of measure 0, the flow must necessarily have an eigenvalue greater than 1 and hence there is sensitive dependence on initial conditions. Further, there is an irregular oscillation whose amplitude is described by a diffusive process that is well-modeled by the Irwin-Hall distribution. There is a large class of other piecewise-linear networks that might be analyzed using similar methods. The analysis gives insight into possible origins of chaotic dynamics in periodically forced dynamical systems.
Nonlinear diffusion equations as asymptotic limits of Cahn-Hilliard systems
Colli, Pierluigi; Fukao, Takeshi
2016-05-01
An asymptotic limit of a class of Cahn-Hilliard systems is investigated to obtain a general nonlinear diffusion equation. The target diffusion equation may reproduce a number of well-known model equations: Stefan problem, porous media equation, Hele-Shaw profile, nonlinear diffusion of singular logarithmic type, nonlinear diffusion of Penrose-Fife type, fast diffusion equation and so on. Namely, by setting the suitable potential of the Cahn-Hilliard systems, all these problems can be obtained as limits of the Cahn-Hilliard related problems. Convergence results and error estimates are proved.
PERTURBATION FINITE VOLUME METHOD FOR CONVECTIVE-DIFFUSION INTEGRAL EQUATION
GAO Zhi; YANG Guowei
2004-01-01
A perturbation finite volume (PFV) method for the convective-diffusion integral equation is developed in this paper. The PFV scheme is an upwind and mixed scheme using any higher-order interpolation and second-order integration approximations, with the least nodes similar to the standard three-point schemes, that is, the number of the nodes needed is equal to unity plus the face-number of the control volume. For instance, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D linear and nonlinear problems, 2-D and 3-D flow model equations. Comparing with other standard three-point schemes, the PFV scheme has much smaller numerical diffusion than the first-order upwind scheme (UDS). Its numerical accuracies are also higher than the second-order central scheme (CDS), the power-law scheme (PLS) and QUICK scheme.
On the supercritically diffusive magneto-geostrophic equations
Friedlander, Susan; Vicol, Vlad
2011-01-01
We address the well-posedness theory for the magento-geostrophic equation, namely an active scalar equation in which the divergence-free drift velocity is one derivative more singular than the active scalar. In the presence of supercritical fractional diffusion given by (-\\Delta)^\\gamma, where 01/2 the equations are locally well-posed, while for \\gamma<1/2 they are ill-posed, in the sense that there is no Lipschitz solution map. The main reason for the striking loss of regularity when \\gamma goes below 1/2 is that the constitutive law used to obtain the velocity from the active scalar is given by an unbounded Fourier multiplier which is both even and anisotropic. Lastly, we note that the anisotropy of the constitutive law for the velocity may be explored in order to obtain an improvement in the regularity of the solutions when the initial data and the force have thin Fourier support, i.e. they are supported on a plane in frequency space. In particular, for such well-prepared data one may prove the local ex...
Travelling Waves in Hybrid Chemotaxis Models
Franz, Benjamin
2013-12-18
Hybrid models of chemotaxis combine agent-based models of cells with partial differential equation models of extracellular chemical signals. In this paper, travelling wave properties of hybrid models of bacterial chemotaxis are investigated. Bacteria are modelled using an agent-based (individual-based) approach with internal dynamics describing signal transduction. In addition to the chemotactic behaviour of the bacteria, the individual-based model also includes cell proliferation and death. Cells consume the extracellular nutrient field (chemoattractant), which is modelled using a partial differential equation. Mesoscopic and macroscopic equations representing the behaviour of the hybrid model are derived and the existence of travelling wave solutions for these models is established. It is shown that cell proliferation is necessary for the existence of non-transient (stationary) travelling waves in hybrid models. Additionally, a numerical comparison between the wave speeds of the continuum models and the hybrid models shows good agreement in the case of weak chemotaxis and qualitative agreement for the strong chemotaxis case. In the case of slow cell adaptation, we detect oscillating behaviour of the wave, which cannot be explained by mean-field approximations. © 2013 Society for Mathematical Biology.
An inherently parallel method for solving discretized diffusion equations
A Monte Carlo approach to solving linear systems of equations is being investigated in the context of the solution of discretized diffusion equations. While the technique was originally devised decades ago, changes in computer architectures (namely, massively parallel machines) have driven the authors to revisit this technique. There are a number of potential advantages to this approach: (1) Analog Monte Carlo techniques are inherently parallel; this is not necessarily true to today's more advanced linear equation solvers (multigrid, conjugate gradient, etc.); (2) Some forms of this technique are adaptive in that they allow the user to specify locations in the problem where resolution is of particular importance and to concentrate the work at those locations; and (3) These techniques permit the solution of very large systems of equations in that matrix elements need not be stored. The user could trade calculational speed for storage if elements of the matrix are calculated on the fly. The goal of this study is to compare the parallel performance of Monte Carlo linear solvers to that of a more traditional parallelized linear solver. The authors observe the linear speedup that they expect from the Monte Carlo algorithm, given that there is no domain decomposition to cause significant communication overhead. Overall, PETSc outperforms the Monte Carlo solver for the test problem. The PETSc parallel performance improves with larger numbers of unknowns for a given number of processors. Parallel performance of the Monte Carlo technique is independent of the size of the matrix and the number of processes. They are investigating modifications to the scheme to accommodate matrix problems with positive off-diagonal elements. They are also currently coding an on-the-fly version of the algorithm to investigate the solution of very large linear systems
Diffusive Boltzmann equation, its fluid dynamics, Couette flow and Knudsen layers
Abramov, Rafail V
2016-01-01
In the current work we propose a diffusive modification of the Boltzmann equation. This naturally leads to the corresponding diffusive fluid dynamics equations, which we numerically investigate in a simple Couette flow setting. This diffusive modification is based on the assumption of the "imperfect" model collision term, which is unable to track all collisions in the corresponding real gas particle system. The effect of missed collisions is then modeled by an appropriately scaled long-term homogenization process of the particle dynamics. The corresponding diffusive fluid dynamics equations are produced in a standard way by closing the hierarchy of the moment equations using either the Euler or the Grad closure. In the numerical experiments with the Couette flow, we discover that the diffusive Euler equations behave similarly to the conventional Navier-Stokes equations, while the diffusive Grad equations additionally exhibit Knudsen-like velocity boundary layers. We compare the simulations with the correspond...
A granular computing method for nonlinear convection-diffusion equation
Tian Ya Lan
2016-01-01
Full Text Available This paper introduces a method of solving nonlinear convection-diffusion equation (NCDE, based on the combination of granular computing (GrC and characteristics finite element method (CFEM. The key idea of the proposed method (denoted as GrC-CFEM is to reconstruct the solution from coarse-grained layer to fine-grained layer. It first gets the nonlinear solution on the coarse-grained layer, and then the function (Taylor expansion is applied to linearize the NCDE on the fine-grained layer. Switch to the fine-grained layer, the linear solution is directly derived from the nonlinear solution. The full nonlinear problem is solved only on the coarse-grained layer. Numerical experiments show that the GrC-CFEM can accelerate the convergence and improve the computational efficiency without sacrificing the accuracy.
Guiding brine shrimp through mazes by solving reaction diffusion equations
Singal, Krishma; Fenton, Flavio
Excitable systems driven by reaction diffusion equations have been shown to not only find solutions to mazes but to also to find the shortest path between the beginning and the end of the maze. In this talk we describe how we can use the Fitzhugh-Nagumo model, a generic model for excitable media, to solve a maze by varying the basin of attraction of its two fixed points. We demonstrate how two dimensional mazes are solved numerically using a Java Applet and then accelerated to run in real time by using graphic processors (GPUs). An application of this work is shown by guiding phototactic brine shrimp through a maze solved by the algorithm. Once the path is obtained, an Arduino directs the shrimp through the maze using lights from LEDs placed at the floor of the Maze. This method running in real time could be eventually used for guiding robots and cars through traffic.
Local multiplicative Schwarz algorithms for convection-diffusion equations
Cai, Xiao-Chuan; Sarkis, Marcus
1995-01-01
We develop a new class of overlapping Schwarz type algorithms for solving scalar convection-diffusion equations discretized by finite element or finite difference methods. The preconditioners consist of two components, namely, the usual two-level additive Schwarz preconditioner and the sum of some quadratic terms constructed by using products of ordered neighboring subdomain preconditioners. The ordering of the subdomain preconditioners is determined by considering the direction of the flow. We prove that the algorithms are optimal in the sense that the convergence rates are independent of the mesh size, as well as the number of subdomains. We show by numerical examples that the new algorithms are less sensitive to the direction of the flow than either the classical multiplicative Schwarz algorithms, and converge faster than the additive Schwarz algorithms. Thus, the new algorithms are more suitable for fluid flow applications than the classical additive or multiplicative Schwarz algorithms.
From baking a cake to solving the diffusion equation
Olszewski, Edward A.
2006-06-01
We explain how modifying a cake recipe by changing either the dimensions of the cake or the amount of cake batter alters the baking time. We restrict our consideration to the génoise and obtain a semiempirical relation for the baking time as a function of oven temperature, initial temperature of the cake batter, and dimensions of the unbaked cake. The relation, which is based on the diffusion equation, has three parameters whose values are estimated from data obtained by baking cakes in cylindrical pans of various diameters. The relation takes into account the evaporation of moisture at the top surface of the cake, which is the dominant factor affecting the baking time of a cake.
Resolution of the time dependent Pn equations by a Godunov type scheme having the diffusion limit
We consider the Pn model to approximate the transport equation in one dimension of space. In a diffusive regime, the solution of this system is solution of a diffusion equation. We are looking for a numerical scheme having the diffusion limit property: in a diffusive regime, it gives the solution of the limiting diffusion equation on a mesh at the diffusion scale. The numerical scheme proposed is an extension of the Godunov type scheme proposed by L. Gosse to solve the P1 model without absorption term. Moreover, it has the well-balanced property: it preserves the steady solutions of the system. (authors)
Approximate Lie group analysis and solutions of 2D nonlinear diffusion-convection equations
Approximate Lie symmetries of the (2+1)-dimensional nonlinear diffusion equation with a small convection are completely classified. It is known that the invariance principle furnishes a systematic method of solving initial-value problems. The solutions of instantaneous source type of the 2D diffusion-convection equation are obtained for the case of power-law diffusivity, using a symmetry reduction
On the uniqueness of semi-wavefronts for non-local delayed reaction-diffusion equations
Aguerrea, Maitere
2013-01-01
We establish the uniqueness of semi-wavefront solution for a non-local delayed reaction-diffusion equation. This result is obtained by using a generalization of the Diekman-Kaper theory for a nonlinear convolution equation. Several applications to the systems of non-local reaction-diffusion equations with distributed time delay are also considered.
Diffusion-equation representations of landform evolution in the simplest circumstances: Appendix C
Hanks, Thomas C.
2009-01-01
The diffusion equation is one of the three great partial differential equations of classical physics. It describes the flow or diffusion of heat in the presence of temperature gradients, fluid flow in porous media in the presence of pressure gradients, and the diffusion of molecules in the presence of chemical gradients. [The other two equations are the wave equation, which describes the propagation of electromagnetic waves (including light), acoustic (sound) waves, and elastic (seismic) waves radiated from earthquakes; and LaPlace’s equation, which describes the behavior of electric, gravitational, and fluid potentials, all part of potential field theory. The diffusion equation reduces to LaPlace’s equation at steady state, when the field of interest does not depend on t. Poisson’s equation is LaPlace’s equation with a source term.
Luchko, Yuri; Povstenko, Yuriy
2012-01-01
In this paper, the one-dimensional time-fractional diffusion-wave equation with the fractional derivative of order $1 \\le \\alpha \\le 2$ is revisited. This equation interpolates between the diffusion and the wave equations that behave quite differently regarding their response to a localized disturbance: whereas the diffusion equation describes a process, where a disturbance spreads infinitely fast, the propagation speed of the disturbance is a constant for the wave equation. For the time fractional diffusion-wave equation, the propagation speed of a disturbance is infinite, but its fundamental solution possesses a maximum that disperses with a finite speed. In this paper, the fundamental solution of the Cauchy problem for the time-fractional diffusion-wave equation, its maximum location, maximum value, and other important characteristics are investigated in detail. To illustrate analytical formulas, results of numerical calculations and plots are presented. Numerical algorithms and programs used to produce pl...
Travelling waves in hybrid chemotaxis models
Franz, Benjamin; Painter, Kevin J; Erban, Radek
2013-01-01
Hybrid models of chemotaxis combine agent-based models of cells with partial differential equation models of extracellular chemical signals. In this paper, travelling wave properties of hybrid models of bacterial chemotaxis are investigated. Bacteria are modelled using an agent-based (individual-based) approach with internal dynamics describing signal transduction. In addition to the chemotactic behaviour of the bacteria, the individual-based model also includes cell proliferation and death. Cells consume the extracellular nutrient field (chemoattractant) which is modelled using a partial differential equation. Mesoscopic and macroscopic equations representing the behaviour of the hybrid model are derived and the existence of travelling wave solutions for these models is established. It is shown that cell proliferation is necessary for the existence of non-transient (stationary) travelling waves in hybrid models. Additionally, a numerical comparison between the wave speeds of the continuum models and the hybr...
QU Chang-Zheng; ZHANG Shun-Li
2004-01-01
The generalized conditionalsymmetry and sign-invariant approaches are developed to study the nonlinear diffusion equations with x-dependent convection and source terms. We obtain conditions under which the equations admit the second-order generalized conditional symmetries and the first-order sign-invariants on the solutions. Several types of different generalized conditional symmetries and first-order sign-invariants for the equations with diffusion of power law are obtained. Exact solutions to the resulting equations are constructed.
Group Analysis of Variable Coefficient Diffusion-Convection Equations. IV. Potential Symmetries
Ivanova, N M; Sophocleous, C
2007-01-01
This paper completes investigation of symmetry properties of nonlinear variable coefficient diffusion-convection equations of the form $f(x)u_t=(g(x)A(u)u_x)_x+h(x)B(u)u_x$. Potential symmetries of equations from the considered class are found and the connection of them with Lie symmetries of diffusion-type equations is shown. Exact solutions of the Fujita--Storm equation $u_t=(u^{-2}u_x)_x$ are constructed.
The stochastic variational method is applied to particle systems and continuum media. As a brief review of this method, we first discuss the application to particle Lagrangians and derive a diffusion-type equation and the Schrödinger equation with the minimum gauge coupling. We further extend the application of the stochastic variational method to Lagrangians of continuum media and show that the Navier–Stokes, Gross–Pitaevskii and generalized diffusion equations are derived. The correction term for the Navier–Stokes equation is also obtained in this method. We discuss the meaning of this correction by comparing it with the diffusion equation. (paper)
Valentyn Tychynin
2015-01-01
Additional nonlocal symmetries of diffusion-convection equations and the Burgers equation are obtained. It is shown that these equations are connected via a generalized hodograph transformation and appropriate nonlocal symmetries arise from additional Lie symmetries of intermediate equations. Two entirely different techniques are used to search nonlocal symmetry of a given equation: the first is based on usage of the characteristic equations generated by additional operators, another techniqu...
The constructive technique and its application in solving a nonlinear reaction diffusion equation
A mathematical technique based on the consideration of a nonlinear partial differential equation together with an additional condition in the form of an ordinary differential equation is employed to study a nonlinear reaction diffusion equation which describes a real process in physics and in chemistry. Several exact solutions for the equation are acquired under certain circumstances. (general)
Parallel computing for homogeneous diffusion and transport equations in neutronics
Parallel computing meets the ever-increasing requirements for neutronic computer code speed and accuracy. In this work, two different approaches have been considered. We first parallelized the sequential algorithm used by the neutronics code CRONOS developed at the French Atomic Energy Commission. The algorithm computes the dominant eigenvalue associated with PN simplified transport equations by a mixed finite element method. Several parallel algorithms have been developed on distributed memory machines. The performances of the parallel algorithms have been studied experimentally by implementation on a T3D Cray and theoretically by complexity models. A comparison of various parallel algorithms has confirmed the chosen implementations. We next applied a domain sub-division technique to the two-group diffusion Eigen problem. In the modal synthesis-based method, the global spectrum is determined from the partial spectra associated with sub-domains. Then the Eigen problem is expanded on a family composed, on the one hand, from eigenfunctions associated with the sub-domains and, on the other hand, from functions corresponding to the contribution from the interface between the sub-domains. For a 2-D homogeneous core, this modal method has been validated and its accuracy has been measured. (author)
A cladding oxidation model based on diffusion equations
During severe accident in PWRs, the cladding oxidation with steam in the core is very important to the accident process. When oxidation time is long, or oxidation occurs in steam starvation conditions, the parabolic rate correlations based on experiments are restricted, which impacts the prediction of cladding failure, hydrogen production, and temperature. According to Fick's laws, a cladding oxidation model in a wide temperature range based on diffusion equations is developed. The developed oxidation model has a wider applicability than those parabolic rate correlations, and can simulate long-term experiments well. The restricted assumptions of short term oxidation time and enough steam environment in the core implemented by those parabolic rate correlations are removed in the model, therefore this model perfectly fit for long-term and steam starvation conditions which are more realistic during a severe accident. This model also can obtain detailed oxygen distribution in the cladding, which is helpful to simulate the cladding failure in detail and develop advanced cladding failure criteria. (authors)
Regularized lattice Boltzmann model for a class of convection-diffusion equations.
Wang, Lei; Shi, Baochang; Chai, Zhenhua
2015-10-01
In this paper, a regularized lattice Boltzmann model for a class of nonlinear convection-diffusion equations with variable coefficients is proposed. The main idea of the present model is to introduce a set of precollision distribution functions that are defined only in terms of macroscopic moments. The Chapman-Enskog analysis shows that the nonlinear convection-diffusion equations can be recovered correctly. Numerical tests, including Fokker-Planck equations, Buckley-Leverett equation with discontinuous initial function, nonlinear convection-diffusion equation with anisotropic diffusion, are carried out to validate the present model, and the results show that the present model is more accurate than some available lattice Boltzmann models. It is also demonstrated that the present model is more stable than the traditional single-relaxation-time model for the nonlinear convection-diffusion equations. PMID:26565368
Diffusion in periodic potential Langevin versus Fokker-Planck equation approach
Phonon activated diffusion of an interstitial impurity in one dimensional cosine potential is discussed with the use both of Langevin (continuous diffusion model) and Fokker-Planck (jump-diffusion model) equations. Jump rate, jump length and diffusion coefficient as a function of temperature at various barrier heights are calculated. There is some difference between results provided by these two models. Therefore the question arises to what extend these two models of diffusion are equivalent. (author)
Macroscopic dynamics of biological cells interacting via chemotaxis and direct contact
Lushnikov, Pavel M; Alber, Mark
2008-01-01
A connection is established between discrete stochastic model describing microscopic motion of fluctuating cells, and macroscopic equations describing dynamics of cellular density. Cells move towards chemical gradient (process called chemotaxis) with their shapes randomly fluctuating. Nonlinear diffusion equation is derived from microscopic dynamics in dimensions one and two using excluded volume approach. Nonlinear diffusion coefficient depends on cellular volume fraction and it is demonstrated to prevent collapse of cellular density. A very good agreement is shown between Monte Carlo simulations of the microscopic Cellular Potts Model and numerical solutions of the macroscopic equations for relatively large cellular volume fractions. Combination of microscopic and macroscopic models were used to simulate growth of structures similar to early vascular networks.
Deuring, Paul; Eymard,, Robert; Mildner, Marcus
2014-01-01
We consider a time-dependent and a steady linear convection-diffusion equation. These equations are approximately solved by a combined finite element -- finite volume method: the diffusion term is discretized by Crouzeix-Raviart piecewise linear finite elements on a triangular grid, and the convection term by upwind barycentric finite volumes. In the unsteady case, the implicit Euler method is used as time discretization. This scheme is shown to be unconditionally L2-stable, uniformly with re...
Innovation diffusion equations on correlated scale-free networks
Bertotti, M. L.; Brunner, J; Modanese, G.
2016-01-01
We introduce a heterogeneous network structure into the Bass diffusion model, in order to study the diffusion times of innovation or information in networks with a scale-free structure, typical of regions where diffusion is sensitive to geographic and logistic influences (like for instance Alpine regions). We consider both the diffusion peak times of the total population and of the link classes. In the familiar trickle-down processes the adoption curve of the hubs is found to anticipate the t...
New Solution of a Fractional Diffusion-Advection Equation Using Ultradistributions
Rocca, M C; Plastino, A R; De Paoli, A L
2014-01-01
In this paper we exactly solve the fractional diffusion-advection equation. For this purpose we use the Theoy of Ultradistributions of J. Sebastiao e Silva, to give a general solution for this equation. From this solution, we obtain several approximations as limiting cases of various situations of physical and astrophysical interest. One of them involves cosmic rays' diffusion.
A Fractional Diffusion Equation for an n-Dimensional Correlated Levy Walk
Taylor-King, J P; Fedotov, S; Van Gorder, R A
2016-01-01
Levy walks define a fundamental concept in random walk theory which allows one to model diffusive spreading that is faster than Brownian motion. They have many applications across different disciplines. However, so far the derivation of a diffusion equation for an n-dimensional correlated Levy walk remained elusive. Starting from a fractional Klein-Kramers equation here we use a moment method combined with a Cattaneo approximation to derive a fractional diffusion equation for superdiffusive short range auto-correlated Levy walks in the large time limit, and solve it. Our derivation discloses different dynamical mechanisms leading to correlated Levy walk diffusion in terms of quantities that can be measured experimentally.
Modeling bacterial chemotaxis inside a cell
Ouannes, Nesrine; Djedi, Noureddine; Luga, Hervé; Duthen, Yves
2014-01-01
This paper describes a bacterial system that reproduces a population of bacteria that behave by simulating the internal reactions of each bacterial cell. The chemotaxis network of a cell is modulated by a hybrid approach that uses an algebraic model for the receptor clusters activity and an ordinary differential equation for the adaptation dynamics. The experiments are defined in order to simulate bacterial growth in an environment where nutrients are regularly added to it. The results show a...
A generalized Fokker–Planck equation for anomalous diffusion in velocity space
A more general quasi-Fokker–Planck equation is derived to describe particle kinetics in situations when the usual Fokker–Planck equation is not applicable. The equation is valid for an arbitrary value of the transferred in a collision act momentum and for the arbitrary mass ratio of the interacting particles. The only assumption is smallness of the typical velocity of the particles, undergoing diffusion. The developed approach is another tool to study anomalous diffusion, avoiding conventional in such problems fractional differentiation. In this Letter anomalous diffusion in velocity space is considered for hard-sphere, Coulomb and dusty plasma collision models. -- Highlights: ► We expand the Master equation into a series. ► We derive a generalized equation for description of anomalous diffusion. ► We consider special cases of anomalous diffusion for hard-sphere interactions, Coulomb systems and dusty plasmas.
Coupled Oscillators with Chemotaxis
Sawai, S; Sawai, Satoshi; Aizawa, Yoji
1998-01-01
A simple coupled oscillator system with chemotaxis is introduced to study morphogenesis of cellular slime molds. The model successfuly explains the migration of pseudoplasmodium which has been experimentally predicted to be lead by cells with higher intrinsic frequencies. Results obtained predict that its velocity attains its maximum value in the interface region between total locking and partial locking and also suggest possible roles played by partial synchrony during multicellular development.
Air Pollution Steady-State Advection-Diffusion Equation: The General Three-Dimensional Solution
Bardo Bodmann; Tiziano Tirabassi; Marco Túllio Vilhena; Daniela Buske
2012-01-01
Atmospheric air pollution turbulent fluxes can be assumed to be proportional to the mean concentration gradient. This assumption, along with the equation of continuity, leads to the advection-diffusion equation. Many models simulating air pollution dispersion are based upon the solution (numerical or analytical) of the advection-diffusion equation as- suming turbulence parameterization for realistic physical scenarios. We present the general steady three-dimensional solution of the advection-...
Koide, T.; Kodama, T.
2011-01-01
The stochastic variational method is applied to particle systems and continuum mediums. As the brief review of this method, we first discuss the application to particle Lagrangians and derive a diffusion-type equation and the Schr\\"{o}dinger equation with the minimum gauge coupling. We further extend the application of the stochastic variational method to Lagrangians of continuum mediums and show that the Navier-Stokes, Gross-Pitaevskii and generalized diffusion equations are derived. The cor...
Frank, T. D.
2008-02-01
We discuss two central claims made in the study by Bassler et al. [K.E. Bassler, G.H. Gunaratne, J.L. McCauley, Physica A 369 (2006) 343]. Bassler et al. claimed that Green functions and Langevin equations cannot be defined for nonlinear diffusion equations. In addition, they claimed that nonlinear diffusion equations are linear partial differential equations disguised as nonlinear ones. We review bottom-up and top-down approaches that have been used in the literature to derive Green functions for nonlinear diffusion equations and, in doing so, show that the first claim needs to be revised. We show that the second claim as well needs to be revised. To this end, we point out similarities and differences between non-autonomous linear Fokker-Planck equations and autonomous nonlinear Fokker-Planck equations. In this context, we raise the question whether Bassler et al.’s approach to financial markets is physically plausible because it necessitates the introduction of external traders and causes. Such external entities can easily be eliminated when taking self-organization principles and concepts of nonextensive thermostatistics into account and modeling financial processes by means of nonlinear Fokker-Planck equations.
Perfect and near perfect adaptation in a model of bacterial chemotaxis
Mello, Bernardo A.; Tu, Yuhai
2002-01-01
The signaling apparatus mediating bacterial chemotaxis can adapt to a wide range of persistent external stimuli. In many cases, the bacterial activity returns to its pre-stimulus level exactly and this "perfect adaptability" is robust against variations in various chemotaxis protein concentrations. We model the bacterial chemotaxis signaling pathway, from ligand binding to CheY phosphorylation. By solving the steady-state equations of the model analytically, we derive a full set of conditions...
In this work we report an analytical solution for the monoenergetic neutron diffusion kinetic equation in cartesian geometry. Bearing in mind that the equation for the delayed neutron precursor concentration is a first order linear differential equation in the time variable, to make possible the application of the GITT approach to the kinetic equation, we introduce a fictitious diffusion term multiplied by a positive small value ε. By this procedure, we are able to solve this set of equations. Indeed, applying the GITT technique to the modified diffusion kinetic equation, we come out with a matrix differential equation which has a well known analytical solution when ε goes to zero. We report numerical simulations as well study of numerical convergence of the results attained. (author)
Asymptotic Speed of Wave Propagation for A Discrete Reaction-Diffusion Equation
Xiu-xiang Liu; Pei-xuan Weng
2006-01-01
We deal with asymptotic speed of wave propagation for a discrete reaction-diffusion equation. We find the minimal wave speed c* from the characteristic equation and show that c* is just the asymptotic speed of wave propagation. The isotropic property and the existence of solution of the initial value problem for the given equation are also discussed.
Darboux transformations for (1+2)-dimensional Fokker-Planck equations with constant diffusion matrix
We construct a Darboux transformation for (1+2)-dimensional Fokker-Planck equations with constant diffusion matrix. Our transformation is based on the two-dimensional supersymmetry formalism for the Schrödinger equation. The transformed Fokker-Planck equation and its solutions are obtained in explicit form.
The second boundary value problem for equations of viscoelastic diffusion in polymers
Vorotnikov, Dmitry A.
2009-01-01
The classical approach to diffusion processes is based on Fick's law that the flux is proportional to the concentration gradient. Various phenomena occurring during propagation of penetrating liquids in polymers show that this type of diffusion exhibits anomalous behavior and contradicts the just mentioned law. However, they can be explained in the framework of non-Fickian diffusion theories based on viscoelasticity of polymers. Initial-boundary value problems for viscoelastic diffusion equat...
On the integrability of the generalized Fisher-type nonlinear diffusion equations
In this paper, the geometric integrability and Lax integrability of the generalized Fisher-type nonlinear diffusion equations with modified diffusion in (1+1) and (2+1) dimensions are studied by the pseudo-spherical surface geometry method and prolongation technique. It is shown that the (1+1)-dimensional Fisher-type nonlinear diffusion equation is geometrically integrable in the sense of describing a pseudo-spherical surface of constant curvature -1 only for m = 2, and the generalized Fisher-type nonlinear diffusion equations in (1+1) and (2+1) dimensions are Lax integrable only for m = 2. This paper extends the results in Bindu et al 2001 (J. Phys. A: Math. Gen. 34 L689) and further provides the integrability information of (1+1)- and (2+1)-dimensional Fisher-type nonlinear diffusion equations for m = 2
Finite element approximation of the Fokker-Planck equation for diffuse optical tomography
In diffuse optical tomography, light transport theory is used to describe photon propagation inside turbid medium. A commonly used simplification for the radiative transport equation is the diffusion approximation due to computational feasibility. However, it is known that the diffusion approximation is not valid close to the sources and boundary and in low-scattering regions. Fokker-Planck equation describes light propagation when scattering is forward-peaked. In this article a numerical solution of the Fokker-Planck equation using finite element method is developed. Approach is validated against Monte Carlo simulation and compared with the diffusion approximation. The results show that the Fokker-Planck equation gives equal or better results than the diffusion approximation on the boundary of a homogeneous medium and in turbid medium containing a low-scattering region when scattering is forward-peaked.
Transformed Fourier and Fick equations for the control of heat and mass diffusion
We review recent advances in the control of diffusion processes in thermodynamics and life sciences through geometric transforms in the Fourier and Fick equations, which govern heat and mass diffusion, respectively. We propose to further encompass transport properties in the transformed equations, whereby the temperature is governed by a three-dimensional, time-dependent, anisotropic heterogeneous convection-diffusion equation, which is a parabolic partial differential equation combining the diffusion equation and the advection equation. We perform two dimensional finite element computations for cloaks, concentrators and rotators of a complex shape in the transient regime. We precise that in contrast to invisibility cloaks for waves, the temperature (or mass concentration) inside a diffusion cloak crucially depends upon time, its distance from the source, and the diffusivity of the invisibility region. However, heat (or mass) diffusion outside cloaks, concentrators and rotators is unaffected by their presence, whatever their shape or position. Finally, we propose simplified designs of layered cylindrical and spherical diffusion cloaks that might foster experimental efforts in thermal and biochemical metamaterials
A semi-analytical finite element method for a class of time-fractional diffusion equations
Sun, HongGuang; Sze, K Y
2011-01-01
As fractional diffusion equations can describe the early breakthrough and the heavy-tail decay features observed in anomalous transport of contaminants in groundwater and porous soil, they have been commonly employed in the related mathematical descriptions. These models usually involve long-time range computation, which is a critical obstacle for its application, improvement of the computational efficiency is of great significance. In this paper, a semi-analytical method is presented for solving a class of time-fractional diffusion equations which overcomes the critical long-time range computation problem of time fractional differential equations. In the procedure, the spatial domain is discretized by the finite element method which reduces the fractional diffusion equations into approximate fractional relaxation equations. As analytical solutions exist for the latter equations, the burden arising from long-time range computation can effectively be minimized. To illustrate its efficiency and simplicity, four...
Innovation diffusion equations on correlated scale-free networks
Bertotti, M L; Modanese, G
2016-01-01
We introduce a heterogeneous network structure into the Bass diffusion model, in order to study the diffusion times of innovation or information in networks with a scale-free structure, typical of regions where diffusion is sensitive to geographic and logistic influences (like for instance Alpine regions). We consider both the diffusion peak times of the total population and of the link classes. In the familiar trickle-down processes the adoption curve of the hubs is found to anticipate the total adoption in a predictable way. In a major departure from the standard model, we model a trickle-up process by introducing heterogeneous publicity coefficients (which can also be negative for the hubs, thus turning them into stiflers) and a stochastic term which represents the erratic generation of innovation at the periphery of the network. The results confirm the robustness of the Bass model and expand considerably its range of applicability.
Innovation diffusion equations on correlated scale-free networks
Bertotti, M. L.; Brunner, J.; Modanese, G.
2016-07-01
We introduce a heterogeneous network structure into the Bass diffusion model, in order to study the diffusion times of innovation or information in networks with a scale-free structure, typical of regions where diffusion is sensitive to geographic and logistic influences (like for instance Alpine regions). We consider both the diffusion peak times of the total population and of the link classes. In the familiar trickle-down processes the adoption curve of the hubs is found to anticipate the total adoption in a predictable way. In a major departure from the standard model, we model a trickle-up process by introducing heterogeneous publicity coefficients (which can also be negative for the hubs, thus turning them into stiflers) and a stochastic term which represents the erratic generation of innovation at the periphery of the network. The results confirm the robustness of the Bass model and expand considerably its range of applicability.
Transition probabilities for diffusion equations by means of path integrals
Goovaerts, Marc; DE SCHEPPER, Ann; Decamps, Marc
2002-01-01
In this paper, we investigate the transition probabilities for diffusion processes. In a first part, we show how transition probabilities for rather general diffusion processes can always be expressed by means of a path integral. For several classical models, an exact calculation is possible, leading to analytical expressions for the transition probabilities and for the maximum probability paths. A second part consists of the derivation of an analytical approximation for the transition probab...
Transition probabilities for diffusion equations by means of path integrals.
Goovaerts, Marc; De Schepper, A; Decamps, M.
2002-01-01
In this paper, we investigate the transition probabilities for diffusion processes. In a first part, we show how transition probabilities for rather general diffusion processes can always be expressed by means of a path integral. For several classical models, an exact calculation is possible, leading to analytical expressions for the transition probabilities and for the maximum probability paths. A second part consists of the derivation of an analytical approximation for the transition probab...
Application of The Full-Sweep AOR Iteration Concept for Space-Fractional Diffusion Equation
Sunarto, A.; Sulaiman, J.; Saudi, A.
2016-04-01
The aim of this paper is to investigate the effectiveness of the Full-Sweep AOR Iterative method by using Full-Sweep Caputo’s approximation equation to solve space-fractional diffusion equations. The governing space-fractional diffusion equations were discretized by using Full-Sweep Caputo’s implicit finite difference scheme to generate a system of linear equations. Then, the Full-Sweep AOR iterative method is applied to solve the generated linear system To examine the application of FSAOR method two numerical tests are conducted to show that the FSAOR method is superior to the FSSOR and FSGS methods.
Crouseilles, Nicolas; Hivert, Hélène; Lemou, Mohammed
2015-01-01
In this work, we propose some numerical schemes for linear kinetic equations in the diffusion and anomalous diffusion limit. When the equilibrium distribution function is a Maxwellian distribution, it is well known that for an appropriate time scale, the small mean free path limit gives rise to a diffusion type equation. However, when a heavy-tailed distribution is considered, another time scale is required and the small mean free path limit leads to a fractional anomalous diffusion equation....
Non-probabilistic solutions of imprecisely defined fractional-order diffusion equations
The fractional diffusion equation is one of the most important partial differential equations (PDEs) to model problems in mathematical physics. These PDEs are more practical when those are combined with uncertainties. Accordingly, this paper investigates the numerical solution of a non-probabilistic viz. fuzzy fractional-order diffusion equation subjected to various external forces. A fuzzy diffusion equation having fractional order 0 < α ≤ 1 with fuzzy initial condition is taken into consideration. Fuzziness appearing in the initial conditions is modelled through convex normalized triangular and Gaussian fuzzy numbers. A new computational technique is proposed based on double parametric form of fuzzy numbers to handle the fuzzy fractional diffusion equation. Using the single parametric form of fuzzy numbers, the original fuzzy diffusion equation is converted first into an interval-based fuzzy differential equation. Next, this equation is transformed into crisp form by using the proposed double parametric form of fuzzy numbers. Finally, the same is solved by Adomian decomposition method (ADM) symbolically to obtain the uncertain bounds of the solution. Computed results are depicted in terms of plots. Results obtained by the proposed method are compared with the existing results in special cases. (general)
A CLASS OF SINGULARLY PERTURBED INITIAL BOUNDARY PROBLEM FOR REACTION DIFFUSION EQUATION
Xie Feng
2003-01-01
The singularly perturbed initial boundary value problem for a class of reaction diffusion equation isconsidered. Under appropriate conditions, the existence-uniqueness and the asymptotic behavior of the solu-tion are showed by using the fixed-point theorem.
Jia-qi Mo; Wan-tao Lin
2006-01-01
A class of nonlinear singularly perturbed problems for reaction diffusion equations with boundary perturbation are considered. Under suitable conditions, the asymptotic behavior of solution for the initial boundary value problems is studied using the theory of differential inequalities.
OuyangCheng; MoJiaqi
2005-01-01
The nonlinear singularly perturbed problems for reaction diffusion equations with boundary perturbation are considered. Under suitable conditions, using the theory of differential inequalities the asymptotic behavior of solution for the initial boundary value problems is studied.
THE CORNER LAYER SOLUTION TO ROBIN PROBLEM FOR REACTION DIFFUSION EQUATION
无
2012-01-01
A class of Robin boundary value problem for reaction diffusion equation is considered. Under suitable conditions, using the theory of differential inequalities the existence and asymptotic behavior of the corner layer solution to the initial boundary value problem are studied.
THE NONLINEAR SINGULARLY PERTURBEDPROBLEMS FOR REACTION DIFFUSION EQUATIONS WITH TIME DELAY
莫嘉琪; 冯茂春
2001-01-01
A class of nonlinear for singularly perturbed problems for reaction diffusion equations with time delays are considered. Under suitable conditions, using theory of differential inequalities the asymptotic behavior of solution for the initial boundary value problems are studied.
ASYMPTOTIC BEHAVIOR OF SOLUTION FOR A CLASS OF REACTION DIFFUSION EQUATIONS
MoJiaqi; LinWantao; ZhuJiang
2004-01-01
A class of initial boundary value problems for the reaction diffusion equations are considered. The asymptotic behavior of solution for the problem is obtained using the theory of differential inequality.
Exact solutions to nonlinear delay reaction-diffusion equations of hyperbolic type
Authors present periodic and antiperiodic solutions, composite solutions resulting from a nonlinear superposition of generalized separable and traveling wave solutions, and others. Some results are extended to nonlinear delay reaction-diffusion equations with time-varying delay
Mixed, Nonsplit, Extended Stability, Stiff Integration of Reaction Diffusion Equations
Alzahrani, Hasnaa H.
2016-07-26
A tailored integration scheme is developed to treat stiff reaction-diffusion prob- lems. The construction adapts a stiff solver, namely VODE, to treat reaction im- plicitly together with explicit treatment of diffusion. The second-order Runge-Kutta- Chebyshev (RKC) scheme is adjusted to integrate diffusion. Spatial operator is de- scretised by second-order finite differences on a uniform grid. The overall solution is advanced over S fractional stiff integrations, where S corresponds to the number of RKC stages. The behavior of the scheme is analyzed by applying it to three simple problems. The results show that it achieves second-order accuracy, thus, preserving the formal accuracy of the original RKC. The presented development sets the stage for future extensions, particularly, to multidimensional reacting flows with detailed chemistry.
A modified diffusion equation for room-acoustic predication.
Jing, Yun; Xiang, Ning
2007-06-01
This letter presents a modified diffusion model using an Eyring absorption coefficient to predict the reverberation time and sound pressure distributions in enclosures. While the original diffusion model [Ollendorff, Acustica 21, 236-245 (1969); J. Picaut et al., Acustica 83, 614-621 (1997); Valeau et al., J. Acoust. Soc. Am. 119, 1504-1513 (2006)] usually has good performance for low absorption, the modified diffusion model yields more satisfactory results for both low and high absorption. Comparisons among the modified model, the original model, a geometrical-acoustics model, and several well-established theories in terms of reverberation times and sound pressure level distributions, indicate significantly improved prediction accuracy by the modification. PMID:17552680
On Nonlinear Nonlocal Systems of Reaction Diffusion Equations
B. Ahmad
2014-01-01
Full Text Available The reaction diffusion system with anomalous diffusion and a balance law ut+-Δα/2u=-fu,v, vt+-∆β/2v=fu,v, 0<α, β<2, is con sidered. The existence of global solutions is proved in two situations: (i a polynomial growth condition is imposed on the reaction term f when 0<α≤β≤2; (ii no growth condition is imposed on the reaction term f when 0<β≤α≤2.
Non-negative mixed finite element formulations for a tensorial diffusion equation
Nakshatrala, K. B.; Valocchi, A. J.
2008-01-01
We consider the tensorial diffusion equation, and address the discrete maximum-minimum principle of mixed finite element formulations. In particular, we address non-negative solutions (which is a special case of the maximum-minimum principle) of mixed finite element formulations. The discrete maximum-minimum principle is the discrete version of the maximum-minimum principle. In this paper we present two non-negative mixed finite element formulations for tensorial diffusion equations based on ...
On the stability of some exact solutions to the generalized convection-reaction-diffusion equation
Highlights: → We investigate stability properties of a set of traveling wave solutions to the convection-reaction-diffusion equation taking into account the effects of memory. → We examine the analytical solution previously obtained within the direct method. → Our analytical studies are backed by the numerical simulations. - Abstract: Stability of a set of traveling wave solutions to the convection-reaction-diffusion equation taking into account the effects of memory is studied by means of the qualitative methods and numerical simulation.
Numerical methods for advection-diffusion-reaction equations and medical applications
Montecinos, Gino Ignacio
2014-01-01
The purpose of this thesis is twofold, firstly, the study of a relaxation procedure for numerically solving advection-diffusion-reaction equations, and secondly, a medical application. Concerning the first topic, we extend the applicability of the Cattaneo relaxation approach to reformulate time-dependent advection-diffusion-reaction equations, that may include stiff reactive terms, as hyperbolic balance laws with stiff source terms. The resulting systems of hyperbolic balance laws are solved...
Sound Field Modeling in Architectural Acoustics using a Diffusion Equation Based Model
Fortin, Nicolas; Picaut, Judicaël; Billon, Alexis; Valeau, Vincent; SAKOUT, Anas
2009-01-01
In this paper, an implementation of a model for room-acoustic predictions in COMSOL Multiphysics is presented. The model (called diffusion model) is based on the solving of diffusion equations instead of classical wave equations and allows simulating the sound propagation in complex geometries at high frequency. Instead of using COMSOL Multiplysics to solve directly the problem, a specific tool has been developed. It is composed of a user-friendly interface (I-Simpa) which manipulates all the...
Sound Field Modeling in Architectural Acoustics using a Diffusion Equation Based Model
Fortin, Nicolas; Picaut, Judicaël; Billon, Alexis; Valeau, Vincent; SAKOUT, Anas
2009-01-01
In this paper, an implementation of a model for room-acoustic predictions in COMSOL Multiphysics is presented. The model (called diffusion model) is based on the solving of diffusion equations instead of classical wave equations and allows simulating the sound propagation in complex geometries at high frequency. Instead of using COMSOL Multiplysics to solve directly the problem, a specific tool has been developed. It is composed of a user-friendly interface (I-Simpa) which manipulates a...
Painleve analysis and exact solutions for the predator-prey system of equations with diffusion
A system of reaction-diffusion equations describing predator-prey relations is considered. The mathematical model corresponds to the modified Lotka-Volterra system with the logistic growth of the prey and with both predator and prey dispersing by diffusion. The Painleve analysis of the system of equations is performed. Exact traveling wave solutions are found by means of the logistic-function method
Umarov, Sabir; Steinberg, Stanly
2009-01-01
In this paper diffusion processes with changing modes are studied involving the variable order partial differential equations. We prove the existence and uniqueness theorem of a solution of the Cauchy problem for fractional variable order (with respect to the time derivative) pseudo-differential equations. Depending on the parameters of variable order derivatives short or long range memories may appear when diffusion modes change. These memory effects are classified and studied in detail. Pro...
Highlights: • Uncertain neutron diffusion equation of bare square homogeneous reactor is studied. • Proposed interval arithmetic is extended for fuzzy numbers. • The developed fuzzy arithmetic is used to handle uncertain parameters. • Governing differential equation is modelled by modified fuzzy finite element method. • Fuzzy critical eigenvalues and effective multiplication factors are investigated. - Abstract: The scattering of neutron collision inside a reactor depends upon geometry of the reactor, diffusion coefficient and absorption coefficient etc. In general these parameters are not crisp and hence we get uncertain neutron diffusion equation. In this paper we have investigated the above equation for a bare square homogeneous reactor. Here the uncertain governing differential equation is modelled by a modified fuzzy finite element method. Using modified fuzzy finite element method, obtained eigenvalues and effective multiplication factors are studied. Corresponding results are compared with the classical finite element method in special cases and various uncertain results have been discussed
The stochastic thermodynamics of a dilute, well-stirred mixture of chemically reacting species is built on the stochastic trajectories of reaction events obtained from the chemical master equation. However, when the molecular populations are large, the discrete chemical master equation can be approximated with a continuous diffusion process, like the chemical Langevin equation or low noise approximation. In this paper, we investigate to what extent these diffusion approximations inherit the stochastic thermodynamics of the chemical master equation. We find that a stochastic-thermodynamic description is only valid at a detailed-balanced, equilibrium steady state. Away from equilibrium, where there is no consistent stochastic thermodynamics, we show that one can still use the diffusive solutions to approximate the underlying thermodynamics of the chemical master equation
熊岳山; 韦永康
2001-01-01
The sediment reaction and diffusion equation with generalized initial and boundary condition is studied. By using Laplace transform and Jordan lemma , an analytical solution is got, which is an extension of analytical solution provided by Cheng Kwokming James ( only diffusion was considered in analytical solution of Cheng ). Some problems arisen in the computation of analytical solution formula are also analysed.
A discrete-ordinate discontinuous-streamline diffusion method for the radiative transfer equation
Wang, Cheng; Sheng, Qiwei; Han, Weimin
2016-01-01
The radiative transfer equation (RTE) arises in many different areas of science and engineering. In this paper, we propose and investigate a discrete-ordinate discontinuous-streamline diffusion (DODSD) method for solving the RTE, which is a combination of the discrete-ordinate technique and the discontinuous-streamline diffusion method. Different from the discrete-ordinate discontinuous Galerkin (DODG) method for the RTE, an artificial diffusion parameter is added to the test functions in the...
An equation describing diffusivity of liquid atoms by magnetic confinement
In this work, we report an obvious low field-induced magnetic confinement effect on the diffusivity in binary metallic melts under a weak magnetic field. A quantitative description of this nontrivial dynamic behavior is given by a physical analytical model based on the Hall effect, which is in agreement with our experimental results. Meanwhile, a quadratic B dependence of the dynamic viscosity obtained in the same confined environment is observed. Our results show that one can effectively control the atomic diffusion process of metallic melts by the application of magnetic field. Meanwhile, this magnetic confinement effect at atomic scale should provide an important new ingredient to deeply understand the condensed matter physics under the external magnetic field. (paper)
On the numerical solution of the one-dimensional convection-diffusion equation
Mehdi Dehghan
2005-01-01
also used to develop new methods of high accuracy. This approach allows simple comparison of the errors associated with the partial differential equation. Various difference approximations are derived for the one-dimensional constant coefficient convection-diffusion equation. The results of a numerical experiment are provided, to verify the efficiency of the designed new algorithms. The paper ends with a concluding remark.
L1 Error Estimates for Difference Approximations Of Degenerate Convection-Diffusion Equations
Karlsen, Kenneth H.; Risebro, Nils Henrik; Storrøsten, Erlend B.
2012-01-01
We analyze monotone difference schemes for strongly degenerate convection-diffusion equations in one spatial dimension. These nonlinear equations are well-posed within a class of (discontinuous) entropy solutions. We prove that the L1 difference between the approximate solutions and the unique entropy solution is bounded above by a constant times the cube root of the spatial discretization parameter.
ZHANG Shun-Li; QU Chang-Zheng
2006-01-01
@@ The concept of approximate generalized conditional symmetry (AGCS) as a generalization to both approximate Lie point symmetry and generalized conditional symmetry is introduced, and it is applied to study the perturbed nonlinear diffusion-convection equations. Complete classification of those perturbed equations which admit cer tain types of AGCSs is derived.
Global Null Controllability of the 1-Dimensional Nonlinear Slow Diffusion Equation
Jean-Michel CORON; Jesús Ildefonso D（I）AZ; Abdelmalek DRICI; Tommaso MINGAZZINI
2013-01-01
The authors prove the global null controllability for the 1-dimensional nonlinear slow diffusion equation by using both a boundary and an internal control.They assume that the internal control is only time dependent.The proof relies on the return method in combination with some local controllability results for nondegenerate equations and rescaling techniques.
Strang-type preconditioners for solving fractional diffusion equations by boundary value methods
Gu, Xian-Ming; Huang, Ting-Zhu; Zhao, Xi-Le; Li, Hou-Biao; Li, Liang
2015-01-01
The finite difference scheme with the shifted Grünwarld formula is employed to semi-discrete the fractional diffusion equations. This spatial discretization can reduce to the large system of ordinary differential equations (ODEs) with initial values. Recently, boundary value method (BVM) was develop
Wave front propagation for a reaction–diffusion equation in narrow random channels
We consider a reaction–diffusion equation in narrow random channels. We approximate the generalized solution to this equation by the corresponding one on a random graph. By making use of large deviation analysis we study the asymptotic wave front propagation. (paper)
Ranjit Kumar; R S Kaushal; Awadhesh Prasad
2010-10-01
An auto-Bäcklund transformation derived in the homogeneous balance method is employed to obtain several new exact solutions of certain kinds of nonlinear diffusion-reaction (D-R) equations. These equations arise in a variety of problems in physical, chemical, biological, social and ecological sciences.
A class of Fokker-Planck equations with logarithmic factors in diffusion and drift terms
The x2 ln (x) dependence of the diffusion coefficient in the Fokker-Planck equation is retrieved by means of symmetry arguments. Exact solutions of the equation with logarithmic factors in coefficients are presented. Algebraic and log-algebraic solutions are found. For some values of exponents they seem to be analogues (on an interval) of log-normal distributions. (author)
Ivanova, N M; Sophocleous, C
2007-01-01
This is the second part of the series of papers on symmetry properties of a class of variable coefficient (1+1)-dimensional nonlinear diffusion-convection equations of general form $f(x)u_t=(g(x)A(u)u_x)_x+h(x)B(u)u_x$. At first, we review the results of Part 1 of the series on equivalence transformations and group classification of the class under consideration. Investigation of non-trivial limits of parameterized subclasses of equations from the given class, which generate contractions of the corresponding maximal Lie invariance algebras, leads to the natural notion of contractions of systems of differential equations. After a brief discussion on contractions of symmetries, equations and solutions in general case, such types of contractions are studied for diffusion--convection equations. A detailed symmetry analysis of an interesting equation from the class under consideration is performed. Exact solutions of some subclasses of the considered class are also given.
Reaction-diffusion equation for quark-hadron transition in heavy-ion collisions
Bagchi, Partha; Sengupta, Srikumar; Srivastava, Ajit M
2015-01-01
Reaction-diffusion equations with suitable boundary conditions have special propagating solutions which very closely resemble the moving interfaces in a first order transition. We show that the dynamics of chiral order parameter for chiral symmetry breaking transition in heavy-ion collisions, with dissipative dynamics, is governed by one such equation, specifically, the Newell-Whitehead equation. Further, required boundary conditions are automatically satisfied due to the geometry of the collision. The chiral transition is, therefore, completed by a propagating interface, exactly as for a first order transition, even though the transition actually is a crossover for relativistic heavy-ion collisions. Same thing also happens when we consider the initial confinement-deconfinement transition with Polyakov loop order parameter. The resulting equation, again with dissipative dynamics, can then be identified with the reaction-diffusion equation known as the Fitzhugh-Nagumo equation which is used in population genet...
Fractional Fokker-Planck Equation and Black-Scholes Formula in Composite-Diffusive Regime
Liang, Jin-Rong; Wang, Jun; Lǔ, Long-Jin; Gu, Hui; Qiu, Wei-Yuan; Ren, Fu-Yao
2012-01-01
In statistical physics, anomalous diffusion plays an important role, whose applications have been found in many areas. In this paper, we introduce a composite-diffusive fractional Brownian motion X α, H ( t)= X H ( S α ( t)), 0Black-Scholes formula. We obtain the fractional Fokker-Planck equation governing the dynamics of the probability density function of the composite-diffusive fractional Brownian motion and find the Black-Scholes differential equation driven by the stock asset X α, H ( t) and the corresponding Black-Scholes formula for the fair prices of European option.
A robust numerical solution of the nonlinear Poisson–Boltzmann equation for asymmetric polyelectrolyte solutions in discrete pore geometries is presented. Comparisons to the linearized approximation of the Poisson–Boltzmann equation reveal that the assumptions leading to linearization may not be appropriate for the electrochemical regime in many cementitious materials. Implications of the electric double layer on both partitioning of species and on diffusive release are discussed. The influence of the electric double layer on anion diffusion relative to cation diffusion is examined.
Lattice fractional diffusion equation in terms of a Riesz-Caputo difference
Wu, Guo-Cheng; Baleanu, Dumitru; Deng, Zhen-Guo; Zeng, Sheng-Da
2015-11-01
A fractional difference is defined by the use of the right and the left Caputo fractional differences. The definition is a two-sided operator of Riesz type and introduces back and forward memory effects in space difference. Then, a fractional difference equation method is suggested for anomalous diffusion in discrete finite domains. A lattice fractional diffusion equation is proposed and the numerical simulation of the diffusion process is discussed for various difference orders. The result shows that the Riesz difference model is particularly suitable for modeling complicated dynamical behaviors on discrete media.
Boundary element method for the solution of the diffusion equation in cylindrical symmetry
Equations for the solution of the diffusion equation in plane Cartesian geometry with the Boundary Element method was derived. The equation for the axi-symmetric case were set and included in the computer program. The results were compared to those obtained by the Finite Difference method. Comparing the results some advantages of the proposed method can be observed, with implications on the multidimensional problems. (author)
Conservation Laws and Exact Solutions for a Reaction-Diffusion Equation with a Variable Coefficient
Zhijie Cao
2014-01-01
Full Text Available In this paper a variable-coefficient reaction-diffusion equation is studied. We classify the equation into three kinds by different restraints imposed on the variable coefficient b(x in the process of solving the determining equations of Lie groups. Then, for each kind, the conservation laws corresponding to the symmetries obtained are considered. Finally, some exact solutions are constructed.
Enhanced Group Analysis of Variable Coefficient Semilinear Diffusion Equations with a Power Source
Vaneeva, O O; Sophocleous, C
2007-01-01
A new approach to group classification problems and more general investigations on transformational properties of classes of differential equations are proposed. It is based on mappings between classes of differential equations, generated by families of point transformations. A class of variable coefficient (1+1)-dimensional semilinear reaction--diffusion equations of the general form $f(x)u_t=(g(x)u_x)_x+h(x)u^m$ ($m\
Inci Cilingir Sungu; Huseyin Demir
2015-01-01
A new application of the hybrid generalized differential transform and finite difference method is proposed by solving time fractional nonlinear reaction-diffusion equations. This method is a combination of the multi-time-stepping temporal generalized differential transform and the spatial finite difference methods. The procedure first converts the time-evolutionary equations into Poisson equations which are then solved using the central difference method. The temporal differential transform ...
Gurhan Gurarslan; Halil Karahan; Devrim Alkaya; Murat Sari; Mutlu Yasar
2013-01-01
This study aims to produce numerical solutions of one-dimensional advection-diffusion equation using a sixth-order compact difference scheme in space and a fourth-order Runge-Kutta scheme in time. The suggested scheme here has been seen to be very accurate and a relatively flexible solution approach in solving the contaminant transport equation for Pe≤5. For the solution of the present equation, the combined technique has been used instead of conventional solution techniques. The accuracy and...
Multi-dimensional traveling fronts in bistable reaction-diffusion equations
谷口 雅治; Masaharu Taniguchi
2012-01-01
Multi-dimensional traveling fronts have been studied in the Allen-Cahn equation (Nagumo equation) and also in multistable reaction-diffusion equations recently. Two-dimensional V-form fronts are studied by Ninomiya and myself (2005) and also by Hamel, Monneau and Roquejoffre (2005). Rotationally symmetric traveling fronts are studied by several authors.In this talk I will give a brief survey on multi-dimensional traveling fronts, and explain what is different and what gives the difficulties c...
Traveling Wave Solutions of Reaction-Diffusion Equations Arising in Atherosclerosis Models
Narcisa Apreutesei
2014-05-01
Full Text Available In this short review article, two atherosclerosis models are presented, one as a scalar equation and the other one as a system of two equations. They are given in terms of reaction-diffusion equations in an infinite strip with nonlinear boundary conditions. The existence of traveling wave solutions is studied for these models. The monostable and bistable cases are introduced and analyzed.
Utilization of Weibull equation to obtain soil-water diffusivity in horizontal infiltration
Water movement was studied in horizontal infiltration experiments using laboratory columns of air-dry and homogeneous soil to obtain a simple and suitable equation for soil-water diffusivity. Many water content profiles for each one of the ten soil columns utilized were obtained through gamma-ray attenuation technique using a 137Cs source. During the measurement of a particular water content profile, the soil column was held in the same position in order to measure changes in time and so to reduce the errors in water content determination. The Weibull equation utilized was excellent in fitting water content profiles experimental data. The use of an analytical function for ν, the Boltzmann variable, according to Weibull model, allowed to obtain a simple equation for soil water diffusivity. Comparisons among the equation here obtained for diffusivity and others solutions found in literature were made, and the unsuitability of a simple exponential variation of diffusivity with water content for the full range of the latter was shown. The necessity of admitting the time dependency for diffusivity was confirmed and also the possibility fixing that dependency on a well known value extended to generalized soil water infiltration studies was found. Finally, it was shown that the soil water diffusivity function given by the equation here proposed can be obtained just by the analysis of the wetting front advance as a function of time. (Author)
Vaneeva, O O; Sophocleous, C
2011-01-01
The group classification of a class of variable coefficient reaction-diffusion equations with exponential nonlinearities is carried out up to both the equivalence generated by the corresponding generalized equivalence group and the general point equivalence. The set of admissible transformations of this class is exhaustively described via finding the complete family of maximal normalized subclasses and associated conditional equivalence groups. Limit processes between variable coefficient reaction-diffusion equations with power nonlinearities and those with exponential nonlinearities are simultaneously studied with limit processes between objects related to these equations (including Lie symmetries, exact solutions and conservation laws).
Fractional diffusion equation for an n -dimensional correlated Lévy walk
Taylor-King, Jake P.; Klages, Rainer; Fedotov, Sergei; Van Gorder, Robert A.
2016-07-01
Lévy walks define a fundamental concept in random walk theory that allows one to model diffusive spreading faster than Brownian motion. They have many applications across different disciplines. However, so far the derivation of a diffusion equation for an n -dimensional correlated Lévy walk remained elusive. Starting from a fractional Klein-Kramers equation here we use a moment method combined with a Cattaneo approximation to derive a fractional diffusion equation for superdiffusive short-range auto-correlated Lévy walks in the large time limit, and we solve it. Our derivation discloses different dynamical mechanisms leading to correlated Lévy walk diffusion in terms of quantities that can be measured experimentally.
Fokker-Planck equations with memory: the cross over from ballistic to diffusive processes
The unified description of diffusion processes that crosses over from a ballistic behavior at short times to a fractional diffusion (sub- or superdiffusion), as well as to the ordinary diffusion at longer times, is proposed on the basis of a non-Markovian generalization of the Fokker-Planck equation. The relations between the non-Markovian kinetic coefficients and observable quantities (mean- and mean square displacements) are established. The problem of calculations of the kinetic coefficients using the Langevin equations is discussed. Solutions of the non-Markovian equation describing diffusive processes in the real (coordinate) space are obtained. For long times, such a solution agrees with results obtained within the continuous random walk theory but is much superior to this solution at shorter times, where the effect of the ballistic region is crucial.
Diffusion equation and spin drag in spin-polarized transport
Flensberg, Karsten; Jensen, Thomas Stibius; Mortensen, Asger
2001-01-01
We study the role of electron-electron interactions for spin-polarized transport using the Boltzmann equation, and derive a set of coupled transport equations. For spin-polarized transport the electron-electron interactions are important, because they tend to equilibrate the momentum of the two-spin...... species. This "spin drag" effect enhances the resistivity of the system. The enhancement is stronger the lower the dimension is, and should be measurable in, for example, a two-dimensional electron gas with ferromagnetic contacts. We also include spin-flip scattering, which has two effects: it...... equilibrates the spin density imbalance and, provided it has a non-s-wave component, also a current imbalance....
Fa, Kwok Sau, E-mail: kwok@dfi.uem.br
2015-02-15
An integro-differential diffusion equation with linear force, based on the continuous time random walk model, is considered. The equation generalizes the ordinary and fractional diffusion equations, which includes short, intermediate and long-time memory effects described by the waiting time probability density function. Analytical expression for the correlation function is obtained and analyzed, which can be used to describe, for instance, internal motions of proteins. The result shows that the generalized diffusion equation has a broad application and it may be used to describe different kinds of systems. - Highlights: • Calculation of the correlation function. • The correlation function is connected to the survival probability. • The model can be applied to the internal dynamics of proteins.
An integro-differential diffusion equation with linear force, based on the continuous time random walk model, is considered. The equation generalizes the ordinary and fractional diffusion equations, which includes short, intermediate and long-time memory effects described by the waiting time probability density function. Analytical expression for the correlation function is obtained and analyzed, which can be used to describe, for instance, internal motions of proteins. The result shows that the generalized diffusion equation has a broad application and it may be used to describe different kinds of systems. - Highlights: • Calculation of the correlation function. • The correlation function is connected to the survival probability. • The model can be applied to the internal dynamics of proteins
Møller, Jan Kloppenborg; Madsen, Henrik
Lamperti transform. This note gives an example driven introduction to the Lamperti transform. The general applicability of the Lamperti transform is limited to univariate diffusion processes, but for a restricted class of multivariate diffusion processes Lamperti type transformations are available and the...
FAYEZ MOUSTAFA MOAWAD, RAGAB
2016-01-01
[EN] The neutron diffusion equation is an approximation of the neutron transport equation that describes the neutron population in a nuclear reactor core. In particular, we will consider here VVER-type reactors which use the neutron diffusion equation discretized on hexagonal meshes. Most of the simulation codes of a nuclear power reactor use the multigroup neutron diffusion equation to describe the neutron distribution inside the reactor core.To study the stationary state of a reactor, the r...
Cargo, P.; Samba, G
2007-07-01
We consider the P{sub n} model to approximate the transport equation in one dimension of space. In a diffusive regime, the solution of this system is solution of a diffusion equation. We are looking for a numerical scheme having the diffusion limit property: in a diffusive regime, it gives the solution of the limiting diffusion equation on a mesh at the diffusion scale. The numerical scheme proposed is an extension of the Godunov type scheme proposed by L. Gosse to solve the P{sub 1} model without absorption term. Moreover, it has the well-balanced property: it preserves the steady solutions of the system. (authors)
Valentyn Tychynin
2015-09-01
Full Text Available Additional nonlocal symmetries of diffusion-convection equations and the Burgers equation are obtained. It is shown that these equations are connected via a generalized hodograph transformation and appropriate nonlocal symmetries arise from additional Lie symmetries of intermediate equations. Two entirely different techniques are used to search nonlocal symmetry of a given equation: the first is based on usage of the characteristic equations generated by additional operators, another technique assumes the reconstruction of a parametrical Lie group transformation from such operator. Some of them are based on the nonlocal transformations that contain new independent variable determined by an auxiliary differential equation and allow the interpretation as a nonlocal transformation with additional variables. The formulae derived for construction of exact solutions are used.
A nonlinear equation for ionic diffusion in a strong binary electrolyte
Ghosal, Sandip; 10.1098/rspa.2010.0028
2012-01-01
The problem of the one dimensional electro-diffusion of ions in a strong binary electrolyte is considered. In such a system the solute dissociates completely into two species of ions with unlike charges. The mathematical description consists of a diffusion equation for each species augmented by transport due to a self consistent electrostatic field determined by the Poisson equation. This mathematical framework also describes other important problems in physics such as electron and hole diffusion across semi-conductor junctions and the diffusion of ions in plasmas. If concentrations do not vary appreciably over distances of the order of the Debye length, the Poisson equation can be replaced by the condition of local charge neutrality first introduced by Planck. It can then be shown that both species diffuse at the same rate with a common diffusivity that is intermediate between that of the slow and fast species (ambipolar diffusion). Here we derive a more general theory by exploiting the ratio of Debye length...
Simulate-HEX - The multi-group diffusion equation in hexagonal-z geometry
The multigroup diffusion equation is solved for the hexagonal-z geometry by dividing each hexagon into 6 triangles. In each triangle, the Fourier solution of the wave equation is approximated by 8 plane waves to describe the intra-nodal flux accurately. In the end an efficient Finite Difference like equation is obtained. The coefficients of this equation depend on the flux solution itself and they are updated once per power/void iteration. A numerical example demonstrates the high accuracy of the method. (authors)
Non-Fickian delay reaction-diffusion equations: theoretical and numerical study
Ferreira, J. A.; Branco, J. R.; Silva, P. da
2007-01-01
The Fisher’s equation is established combining the Fick’s law for the flux and the mass conservation law. Assuming that the reaction term depends on the solution at some past time, a delay parameter is introduced and the delay Fisher’s equation is obtained. Modifying the Fick’s law for the flux considering a temporal memory term, integro-differential equations of Volterra type were introduced in the literature. In these paper we study reaction-diffusion equations obtained co...
Conservation Laws of a Family of Reaction-Diffusion-Convection Equations
Bruzón, M. S.; Gandarias, M. L.; de la Rosa, R.
Ibragimov introduced the concept of nonlinear self-adjoint equations. This definition generalizes the concept of self-adjoint and quasi-self-adjoint equations. Gandarias defined the concept of weak self-adjoint. In this paper, we found a class of nonlinear self-adjoint nonlinear reaction-diffusion-convection equations which are neither self-adjoint nor quasi-self-adjoint nor weak self-adjoint. From a general theorem on conservation laws proved by Ibragimov we obtain conservation laws for these equations.
On the sharp front-type solution of the Nagumo equation with nonlinear diffusion and convection
M B A Mansour
2013-03-01
This paper is concerned with the Nagumo equation with nonlinear degenerate diffusion and convection which arises in several problems of population dynamics, chemical reactions and others. A sharp front-type solution with a minimum speed to this model equation is analysed using different methods. One of the methods is to solve the travelling wave equations and compute an exact solution which describes the sharp travelling wavefront. The second method is to solve numerically an initial-moving boundary-value problem for the partial differential equation and obtain an approximation for this sharp front-type solution.
An Efficient Implicit FEM Scheme for Fractional-in-Space Reaction-Diffusion Equations
Burrage, Kevin
2012-01-01
Fractional differential equations are becoming increasingly used as a modelling tool for processes associated with anomalous diffusion or spatial heterogeneity. However, the presence of a fractional differential operator causes memory (time fractional) or nonlocality (space fractional) issues that impose a number of computational constraints. In this paper we develop efficient, scalable techniques for solving fractional-in-space reaction diffusion equations using the finite element method on both structured and unstructured grids via robust techniques for computing the fractional power of a matrix times a vector. Our approach is show-cased by solving the fractional Fisher and fractional Allen-Cahn reaction-diffusion equations in two and three spatial dimensions, and analyzing the speed of the traveling wave and size of the interface in terms of the fractional power of the underlying Laplacian operator. © 2012 Society for Industrial and Applied Mathematics.
Diffusive approximation of a time-fractional Burgers equation in nonlinear acoustics
Lombard, Bruno
2016-01-01
A fractional time derivative is introduced into the Burgers equation to model losses of nonlinear waves. This term amounts to a time convolution product, which greatly penalizes the numerical modeling. A diffusive representation of the fractional derivative is adopted here, replacing this nonlocal operator by a continuum of memory variables that satisfy local-in-time ordinary differential equations. Then a quadrature formula yields a system of local partial differential equations, well-suited to numerical integration. The determination of the quadrature coefficients is crucial to ensure both the well-posedness of the system and the computational efficiency of the diffusive approximation. For this purpose, optimization with constraint is shown to be a very efficient strategy. Strang splitting is used to solve successively the hyperbolic part by a shock-capturing scheme, and the diffusive part exactly. Numerical experiments are proposed to assess the efficiency of the numerical modeling, and to illustrate the e...
Second order time evolution of the multigroup diffusion and P1 equations for radiation transport
Highlights: → An existing multigroup transport algorithm is extended to be second-order in time. → A new algorithm is presented that does not require a grey acceleration solution. → The two algorithms are tested with 2D, multi-material problems. → The two algorithms have comparable computational requirements. - Abstract: An existing solution method for solving the multigroup radiation equations, linear multifrequency-grey acceleration, is here extended to be second order in time. This method works for simple diffusion and for flux-limited diffusion, with or without material conduction. A new method is developed that does not require the solution of an averaged grey transport equation. It is effective solving both the diffusion and P1 forms of the transport equation. Two dimensional, multi-material test problems are used to compare the solution methods.
Gal, Ciprian G.; Warma, Mahamadi
2016-08-01
We investigate the long term behavior in terms of finite dimensional global and exponential attractors, as time goes to infinity, of solutions to a semilinear reaction-diffusion equation on non-smooth domains subject to nonlocal Robin boundary conditions, characterized by the presence of fractional diffusion on the boundary. Our results are of general character and apply to a large class of irregular domains, including domains whose boundary is Hölder continuous and domains which have fractal-like geometry. In addition to recovering most of the existing results on existence, regularity, uniqueness, stability, attractor existence, and dimension, for the well-known reaction-diffusion equation in smooth domains, the framework we develop also makes possible a number of new results for all diffusion models in other non-smooth settings.
Biomixing by chemotaxis and efficiency of biological reactions: the critical reaction case
Kiselev, Alexander
2012-01-01
Many phenomena in biology involve both reactions and chemotaxis. These processes can clearly influence each other, and chemotaxis can play an important role in sustaining and speeding up the reaction. In continuation of our earlier work, we consider a model with a single density function involving diffusion, advection, chemotaxis, and absorbing reaction. The model is motivated, in particular, by the studies of coral broadcast spawning, where experimental observations of the efficiency of fertilization rates significantly exceed the data obtained from numerical models that do not take chemotaxis (attraction of sperm gametes by a chemical secreted by egg gametes) into account. We consider the case of the weakly coupled quadratic reaction term, which is the most natural from the biological point of view and was left open. The result is that similarly to higher power coupling, the chemotaxis plays a crucial role in ensuring efficiency of reaction. However, mathematically, the picture is quite different in the qua...
A Bloch-Torrey Equation for Diffusion in a Deforming Media
Diffusion Tensor Magnetic Resonance Imaging (DTMRI)technique enables the measurement of diffusion parameters and therefore, informs on the structure of the biological tissue. This technique is applied with success to the static organs such as brain. However, the diffusion measurement on the dynamically deformable organs such as the in-vivo heart is a complex problem that has however a great potential in the measurement of cardiac health. In order to understand the behavior of the Magnetic Resonance (MR)signal in a deforming media, the Bloch-Torrey equation that leads the MR behavior is expressed in general curvilinear coordinates. These coordinates enable to follow the heart geometry and deformations through time. The equation is finally discredited and presented in a numerical formulation using implicit methods, in order to get a stable scheme that can be applied to any smooth deformations. Diffusion process enables the link between the macroscopic behavior of molecules and the microscopic structure in which they evolve. The measurement of diffusion in biological tissues is therefore of major importance in understanding the complex underlying structure that cannot be studied directly. The Diffusion Tensor Magnetic Resonance Imaging(DTMRI) technique enables the measurement of diffusion parameters and therefore provides information on the structure of the biological tissue. This technique has been applied with success to static organs such as the brain. However, diffusion measurement of dynamically deformable organs such as the in-vivo heart remains a complex problem, which holds great potential in determining cardiac health. In order to understand the behavior of the magnetic resonance (MR) signal in a deforming media, the Bloch-Torrey equation that defines the MR behavior is expressed in general curvilinear coordinates. These coordinates enable us to follow the heart geometry and deformations through time. The equation is finally discredited and presented in a
Non-local reaction-diffusion equations modelling predator-prey coevolution
Calsina, Ángel; Calsina i Ballesta, Àngel; Perelló Valls, Carles
1994-01-01
In this paper we examine a prey-predator system with a characteristic of the predator subject to mutation. The ultimate equilibrium of the system is found by Maynard-Smith et al. by the so called ESS (Evolutionary Stable Strategy). Using a system of reaction-diffusion equations with non local terms, we conclude the ESS result for the diffusion coefficient tending to zero, without resorting to any optimization criterion.
Incremental Unknowns Method for Solving Three-Dimensional Convection-Diffusion Equations
Lunji Song; Yujiang Wu
2007-01-01
We use the incremental unknowns method in conjunction with the iterative methods to approximate the solution of the nonsymmetric and positive-definite linear systems generated from a multilevel discretization of three-dimensional convection-diffusion equations. The condition numbers of incremental unknowns matrices associated with the convection-diffusion equations and the number of iterations needed to attain an acceptable accuracy are estimated. Numerical results are presented with two-level approximations,which demonstrate that the incremental unknowns method when combined with some iterative methods is very efficient.
A Fully Discrete Galerkin Method for a Nonlinear Space-Fractional Diffusion Equation
Yunying Zheng
2011-01-01
Full Text Available The spatial transport process in fractal media is generally anomalous. The space-fractional advection-diffusion equation can be used to characterize such a process. In this paper, a fully discrete scheme is given for a type of nonlinear space-fractional anomalous advection-diffusion equation. In the spatial direction, we use the finite element method, and in the temporal direction, we use the modified Crank-Nicolson approximation. Here the fractional derivative indicates the Caputo derivative. The error estimate for the fully discrete scheme is derived. And the numerical examples are also included which are in line with the theoretical analysis.
An analytical solution for the two-group kinetic neutron diffusion equation in cylindrical geometry
Recently the two-group Kinetic Neutron Diffusion Equation with six groups of delay neutron precursor in a rectangle was solved by the Laplace Transform Technique. In this work, we report on an analytical solution for this sort of problem but in cylindrical geometry, assuming a homogeneous and infinite height cylinder. The solution is obtained applying the Hankel Transform to the Kinetic Diffusion equation and solving the transformed problem by the same procedure used in the rectangle. We also present numerical simulations and comparisons against results available in literature. (author)
Analytical versus discretized solutions of four-group diffusion equations to thermal reactors
This paper presents the application of four-group Diffusion theory to thermal reactor criticality calculation. The four-group diffusion equations are applied to the spherical nucleus and reflector of an example reactor. The neutrons fluxes depend upon the radial coordinate. The simultaneous linear ordinary differential equations are solved given the solutions for the fluxes. The neutron fluxes for the nucleus are functions of the eight functions linearly independent consisting of sin, cos, sinh, cosh, sin sinh, sin cosh, cos sinh, and cos cosh. The analytical and discretized calculations of keff value give excellent agreement, an error around 0,03%. (author)
An analytical solution for the two-group kinetic neutron diffusion equation in cylindrical geometry
Fernandes, Julio Cesar L.; Vilhena, Marco Tullio, E-mail: julio.lombaldo@ufrgs.br, E-mail: vilhena@pq.cnpq.br [Programa de Pos Graduacao em Matematica Aplicada (DMPA/UFRGS), Universidade Federal do Rio Grande do Sul Porto Alegre, RS (Brazil); Bodmann, Bardo Ernst, E-mail: bardo.bodmann@ufrgs.br [Programa de Pos-Graduacao em Engenharia Mecanica (PROMEC/UFRGS), Universidade Federal do Rio Grande do Sul, Porto Alegre, RS (Brazil)
2011-07-01
Recently the two-group Kinetic Neutron Diffusion Equation with six groups of delay neutron precursor in a rectangle was solved by the Laplace Transform Technique. In this work, we report on an analytical solution for this sort of problem but in cylindrical geometry, assuming a homogeneous and infinite height cylinder. The solution is obtained applying the Hankel Transform to the Kinetic Diffusion equation and solving the transformed problem by the same procedure used in the rectangle. We also present numerical simulations and comparisons against results available in literature. (author)
We develop a piecewise linear (PWL) Galerkin finite element spatial discretization for the multi-dimensional radiation diffusion equation. It uses piecewise linear weight and basis functions in the finite element approximation, and it can be applied on arbitrary polygonal (2-dimensional) or polyhedral (3-dimensional) grids. We show that this new PWL method gives solutions comparable to those from Palmer's finite-volume method. However, since the PWL method produces a symmetric positive definite coefficient matrix, it should be substantially more computationally efficient than Palmer's method, which produces an asymmetric matrix. We conclude that the Galerkin PWL method is an attractive option for solving diffusion equations on unstructured grids. (authors)