HOME

TheInfoList



OR:

The diffusion equation is a parabolic partial differential equation. In physics, it describes the macroscopic behavior of many micro-particles in
Brownian motion Brownian motion, or pedesis (from grc, πήδησις "leaping"), is the random motion of particles suspended in a medium (a liquid or a gas). This pattern of motion typically consists of random fluctuations in a particle's position insi ...
, resulting from the random movements and collisions of the particles (see Fick's laws of diffusion). In mathematics, it is related to
Markov processes Markov ( Bulgarian, russian: Марков), Markova, and Markoff are common surnames used in Russia and Bulgaria. Notable people with the name include: Academics *Ivana Markova (born 1938), Czechoslovak-British emeritus professor of psychology at ...
, such as
random walks In mathematics, a random walk is a random process that describes a path that consists of a succession of random steps on some mathematical space. An elementary example of a random walk is the random walk on the integer number line \mathbb Z ...
, and applied in many other fields, such as materials science,
information theory Information theory is the scientific study of the quantification (science), quantification, computer data storage, storage, and telecommunication, communication of information. The field was originally established by the works of Harry Nyquist a ...
, and
biophysics Biophysics is an interdisciplinary science that applies approaches and methods traditionally used in physics to study biological phenomena. Biophysics covers all scales of biological organization, from molecular to organismic and populations. ...
. The diffusion equation is a special case of the
convection–diffusion equation The convection–diffusion equation is a combination of the diffusion equation, diffusion and convection (advection equation, advection) equations, and describes physical phenomena where particles, energy, or other physical quantities are transferr ...
, when bulk velocity is zero. It is equivalent to the heat equation under some circumstances.


Statement

The equation is usually written as: where is the
density Density (volumetric mass density or specific mass) is the substance's mass per unit of volume. The symbol most often used for density is ''ρ'' (the lower case Greek letter rho), although the Latin letter ''D'' can also be used. Mathematical ...
of the diffusing material at location and time and is the collective
diffusion coefficient Diffusivity, mass diffusivity or diffusion coefficient is a proportionality constant between the molar flux due to molecular diffusion and the gradient in the concentration of the species (or the driving force for diffusion). Diffusivity is enco ...
for density at location ; and represents the vector differential operator
del Del, or nabla, is an operator used in mathematics (particularly in vector calculus) as a vector differential operator, usually represented by the nabla symbol ∇. When applied to a function defined on a one-dimensional domain, it denotes ...
. If the diffusion coefficient depends on the density then the equation is nonlinear, otherwise it is linear. The equation above applies when the diffusion coefficient is
isotropic Isotropy is uniformity in all orientations; it is derived . Precise definitions depend on the subject area. Exceptions, or inequalities, are frequently indicated by the prefix ' or ', hence ''anisotropy''. ''Anisotropy'' is also used to describe ...
; in the case of anisotropic diffusion, is a symmetric
positive definite matrix In mathematics, a symmetric matrix M with real entries is positive-definite if the real number z^\textsfMz is positive for every nonzero real column vector z, where z^\textsf is the transpose of More generally, a Hermitian matrix (that is, a c ...
, and the equation is written (for three dimensional diffusion) as: If is constant, then the equation reduces to the following
linear differential equation In mathematics, a linear differential equation is a differential equation that is defined by a linear polynomial in the unknown function and its derivatives, that is an equation of the form :a_0(x)y + a_1(x)y' + a_2(x)y'' \cdots + a_n(x)y^ = b ...
: :\frac = D\nabla^2\phi(\mathbf,t), which is identical to the heat equation.


Historical origin

The particle diffusion equation was originally derived by
Adolf Fick Adolf Eugen Fick (3 September 1829 – 21 August 1901) was a German-born physician and physiologist. Early life and education Fick began his work in the formal study of mathematics and physics before realising an aptitude for medicine. He th ...
in 1855.


Derivation

The diffusion equation can be trivially derived from the continuity equation, which states that a change in density in any part of the system is due to inflow and outflow of material into and out of that part of the system. Effectively, no material is created or destroyed: \frac+\nabla\cdot\mathbf=0, where j is the flux of the diffusing material. The diffusion equation can be obtained easily from this when combined with the phenomenological Fick's first law, which states that the flux of the diffusing material in any part of the system is proportional to the local density gradient: \mathbf=-D(\phi,\mathbf)\,\nabla\phi(\mathbf,t). If drift must be taken into account, the Fokker–Planck equation provides an appropriate generalization.


Discretization

The diffusion equation is continuous in both space and time. One may discretize space, time, or both space and time, which arise in application. Discretizing time alone just corresponds to taking time slices of the continuous system, and no new phenomena arise. In discretizing space alone, the
Green's function In mathematics, a Green's function is the impulse response of an inhomogeneous linear differential operator defined on a domain with specified initial conditions or boundary conditions. This means that if \operatorname is the linear differenti ...
becomes the
discrete Gaussian kernel In the areas of computer vision, image analysis and signal processing, the notion of scale-space representation is used for processing measurement data at multiple scales, and specifically enhance or suppress image features over different ranges ...
, rather than the continuous
Gaussian kernel In mathematics, a Gaussian function, often simply referred to as a Gaussian, is a function of the base form f(x) = \exp (-x^2) and with parametric extension f(x) = a \exp\left( -\frac \right) for arbitrary real constants , and non-zero . It is ...
. In discretizing both time and space, one obtains the
random walk In mathematics, a random walk is a random process that describes a path that consists of a succession of random steps on some mathematical space. An elementary example of a random walk is the random walk on the integer number line \mathbb Z ...
.


Discretization (Image)

The
product rule In calculus, the product rule (or Leibniz rule or Leibniz product rule) is a formula used to find the derivatives of products of two or more functions. For two functions, it may be stated in Lagrange's notation as (u \cdot v)' = u ' \cdot v ...
is used to rewrite the anisotropic tensor diffusion equation, in standard discretization schemes, because direct discretization of the diffusion equation with only first order spatial central differences leads to checkerboard artifacts. The rewritten diffusion equation used in image filtering: \frac = \nabla\cdot \left (\phi,\mathbf)\right\nabla \phi(\mathbf,t) + \Big D(\phi,\mathbf)\big(\nabla\nabla^T \phi(\mathbf,t)\big)\Big where "tr" denotes the
trace Trace may refer to: Arts and entertainment Music * ''Trace'' (Son Volt album), 1995 * ''Trace'' (Died Pretty album), 1993 * Trace (band), a Dutch progressive rock band * ''The Trace'' (album) Other uses in arts and entertainment * ''Trace'' ...
of the 2nd rank
tensor In mathematics, a tensor is an algebraic object that describes a multilinear relationship between sets of algebraic objects related to a vector space. Tensors may map between different objects such as vectors, scalars, and even other tenso ...
, and superscript "''T''" denotes
transpose In linear algebra, the transpose of a matrix is an operator which flips a matrix over its diagonal; that is, it switches the row and column indices of the matrix by producing another matrix, often denoted by (among other notations). The tr ...
, in which in image filtering ''D''(''ϕ'', r) are symmetric matrices constructed from the
eigenvectors In linear algebra, an eigenvector () or characteristic vector of a linear transformation is a nonzero vector that changes at most by a scalar factor when that linear transformation is applied to it. The corresponding eigenvalue, often denoted ...
of the image
structure tensor In mathematics, the structure tensor, also referred to as the second-moment matrix, is a matrix derived from the gradient of a function. It describes the distribution of the gradient in a specified neighborhood around a point and makes the inf ...
s. The spatial derivatives can then be approximated by two first order and a second order central
finite difference A finite difference is a mathematical expression of the form . If a finite difference is divided by , one gets a difference quotient. The approximation of derivatives by finite differences plays a central role in finite difference methods for t ...
s. The resulting diffusion algorithm can be written as an image
convolution In mathematics (in particular, functional analysis), convolution is a operation (mathematics), mathematical operation on two function (mathematics), functions ( and ) that produces a third function (f*g) that expresses how the shape of one is ...
with a varying kernel (stencil) of size 3 × 3 in 2D and 3 × 3 × 3 in 3D.


See also

* Continuity equation * Heat equation * Fokker–Planck equation * Fick's laws of diffusion * Maxwell–Stefan equation *
Radiative transfer equation and diffusion theory for photon transport in biological tissue Photon transport in biological tissue can be equivalently modeled numerically with Monte Carlo simulations or analytically by the radiative transfer equation (RTE). However, the RTE is difficult to solve without introducing approximations. A common ...
* Streamline diffusion *
Numerical solution of the convection–diffusion equation The convection–diffusion equation describes the flow of heat, particles, or other physical quantities in situations where there is both diffusion and convection or advection. For information about the equation, its derivation, and its conceptual i ...


References


Further reading

*Carslaw, H. S. and Jaeger, J. C. (1959). ''Conduction of Heat in Solids''. Oxford: Clarendon Press *Crank, J. (1956). ''The Mathematics of Diffusion''. Oxford: Clarendon Press *Mathews, Jon; Walker, Robert L. (1970). ''Mathematical methods of physics'' (2nd ed.), New York: W. A. Benjamin, *Thambynayagam, R. K. M (2011). ''The Diffusion Handbook: Applied Solutions for Engineers''. McGraw-Hill


External links


Diffusion Calculator for Impurities & Dopants in Silicon

A tutorial on the theory behind and solution of the Diffusion Equation.


{{DEFAULTSORT:Diffusion Equation Diffusion Partial differential equations Parabolic partial differential equations Functions of space and time it:Leggi di Fick