Rdxfbgrtj

The diffusion equation is a partial differential equation. In physics, it describes the behavior of the collective motion of micro-particles in a material resulting from the random movement of each micro-particle. In mathematics, it is applicable in common to a subject relevant to the Markov process as well as in various other fields, such as the materials sciences, information science, life science, social science, and so on. These subjects described by the diffusion equation are generally called Brown problems.

Statement

The equation is usually written as:

∂ϕ(r,t)∂t=∇⋅[D(ϕ,r) ∇ϕ(r,t)],{displaystyle {frac {partial phi (mathbf {r} ,t)}{partial t}}=nabla cdot {big [}D(phi ,mathbf {r} ) nabla phi (mathbf {r} ,t){big ]},}

where ϕ(r, t) is the density of the diffusing material at location r and time t and D(ϕ, r) is the collective diffusion coefficient for density ϕ at location r; and ∇ represents the vector differential operator del. If the diffusion coefficient depends on the density then the equation is nonlinear, otherwise it is linear.

More generally, when D is a symmetric positive definite matrix, the equation describes anisotropic diffusion, which is written (for three dimensional diffusion) as:

∂ϕ(r,t)∂t=∑i=13∑j=13∂∂xi[Dij(ϕ,r)∂ϕ(r,t)∂xj]{displaystyle {frac {partial phi (mathbf {r} ,t)}{partial t}}=sum _{i=1}^{3}sum _{j=1}^{3}{frac {partial }{partial x_{i}}}left[D_{ij}(phi ,mathbf {r} ){frac {partial phi (mathbf {r} ,t)}{partial x_{j}}}right]}

If D is constant, then the equation reduces to the following linear differential equation:

∂ϕ(r,t)∂t=D∇2ϕ(r,t),{displaystyle {frac {partial phi (mathbf {r} ,t)}{partial t}}=Dnabla ^{2}phi (mathbf {r} ,t),}

also called the heat equation.

Historical origin

The particle diffusion equation was originally derived by Adolf Fick in 1855.^[1]

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:

∂ϕ∂t+∇⋅j=0,{displaystyle {frac {partial phi }{partial t}}+nabla cdot mathbf {j} =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:

j=−D(ϕ,r)∇ϕ(r,t).{displaystyle mathbf {j} =-D(phi ,mathbf {r} ),nabla phi (mathbf {r} ,t).}

If drift must be taken into account, the Smoluchowski 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 becomes the discrete Gaussian kernel, rather than the continuous Gaussian kernel. In discretizing both time and space, one obtains the random walk.

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:

∂ϕ(r,t)∂t=∇⋅[D(ϕ,r)]∇ϕ(r,t)+tr[D(ϕ,r)(∇∇Tϕ(r,t))]{displaystyle {frac {partial phi (mathbf {r} ,t)}{partial t}}=nabla cdot left[D(phi ,mathbf {r} )right]nabla phi (mathbf {r} ,t)+{rm {tr}}{Big [}D(phi ,mathbf {r} ){big (}nabla nabla ^{T}phi (mathbf {r} ,t){big )}{Big ]}}

where "tr" denotes the trace of the 2nd rank tensor, and superscript "T" denotes transpose, 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 differences. The resulting diffusion algorithm can be written as an image convolution 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 law of diffusion: Fick's Second Law
  


• 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 Jager, 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,
  ISBN 0-8053-7002-1
  


• 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.


• Classical and nanoscale diffusion (with figures and animations)