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
Fick's laws of diffusion describe diffusion and were derived by Adolf Fick in 1855. They can be used to solve for the diffusion coefficient, . Fick's first law can be used to derive his second law which in turn is identical to the diffusion ...
). In mathematics, it is related to
Markov processes, such as
random walks, and applied in many other fields, such as
materials science,
information theory
Information theory is the scientific study of the quantification, storage, and communication of information. The field was originally established by the works of Harry Nyquist and Ralph Hartley, in the 1920s, and Claude Shannon in the 1940s. ...
, 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 and convection (advection) equations, and describes physical phenomena where particles, energy, or other physical quantities are transferred inside a physical system due to two ...
, 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. Mathematicall ...
of the diffusing material at location and time and is the collective
diffusion coefficient for density at location ; and represents the vector
differential operator del. 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; in the case of anisotropic diffusion, is a symmetric
positive definite matrix, and the equation is written (for three dimensional diffusion) as:
If is constant, then the equation reduces to the following
linear differential equation:
:
which is identical to the
heat equation.
Historical origin
The
particle diffusion equation was originally derived by
Adolf Fick in 1855.
Derivation
The diffusion equation can be trivially derived from the
continuity equation
A continuity equation or transport equation is an equation that describes the transport of some quantity. It is particularly simple and powerful when applied to a conserved quantity, but it can be generalized to apply to any extensive quantity. ...
, 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:
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:
If drift must be taken into account, the
Fokker–Planck equation
In statistical mechanics, the Fokker–Planck equation is a partial differential equation that describes the time evolution of the probability density function of the velocity of a particle under the influence of drag forces and random forces, a ...
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 becomes the
discrete Gaussian kernel, rather than the continuous
Gaussian kernel. 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 ...
.
Discretization (Image)
The
product rule 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:
where "tr" denotes the
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 tens ...
, 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 of the image
structure tensors. 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 mathematical operation on two functions ( and ) that produces a third function (f*g) that expresses how the shape of one is modified by the other. The term ''convolution' ...
with a varying kernel (stencil) of size 3 × 3 in 2D and 3 × 3 × 3 in 3D.
See also
*
Continuity equation
A continuity equation or transport equation is an equation that describes the transport of some quantity. It is particularly simple and powerful when applied to a conserved quantity, but it can be generalized to apply to any extensive quantity. ...
*
Heat equation
*
Fokker–Planck equation
In statistical mechanics, the Fokker–Planck equation is a partial differential equation that describes the time evolution of the probability density function of the velocity of a particle under the influence of drag forces and random forces, a ...
*
Fick's laws of diffusion
Fick's laws of diffusion describe diffusion and were derived by Adolf Fick in 1855. They can be used to solve for the diffusion coefficient, . Fick's first law can be used to derive his second law which in turn is identical to the diffusion ...
*
Maxwell–Stefan equation
*
Radiative transfer equation and diffusion theory for photon transport in biological tissue
*
Streamline diffusion
*
Numerical solution of the convection–diffusion equation
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