Diffusion tensor 2026-11-05
A diffusion tensor describes direction-dependent diffusive flux. In the local, linear Fickian model,where is concentration and the entries of have the units of a diffusion coefficient. A scalar coefficient is recovered when . Material models can instead specify a diagonal matrix or symmetric matrix to represent anisotropic transport.[ref]
Diffusion in solids in the COMSOL 6.4 documentation. Solid, Diffusion: isotropic, diagonal and symmetric diffusion coefficients
For an original two-dimensional example, chooseA concentration gradient directed along , with components , produces diffusive flux in consistent units. The flux need not be antiparallel to the gradient. Nevertheless,for every nonzero gradient in this positive-definite matrix example.
The eigenvectors and have eigenvalues 3 and 1. They identify the principal directions and coefficients. Rotating the coordinate basis changes the matrix entries according to ; it does not change these principal coefficients. In an unbounded, homogeneous, drift-free Gaussian model, displacement covariance is , so spreading also differs along those directions.
Symmetry and positive definiteness are assumptions of this ordinary reciprocal diffusion model, not a definition covering every transport law called a tensor. Boundaries also matter: a reflecting boundary condition constrains , which need not be proportional to the normal derivative alone.
Restricted diffusion 2026-11-05
Restricted diffusion is diffusion whose displacement statistics are constrained by boundaries that particles cannot cross. A reflecting boundary condition models an impermeable wall: particles continue moving, but there is no normal diffusive flux through it. Finite permeability would require a different boundary model.
A physical demonstration uses a liquid confined between solid lamellae. In his first-person account, John E. Tanner described constructing 99 mica layers separated by aluminium-foil strips, removing wrinkles and burrs to obtain regular spacing; Edward O. Stejskal inspected the stack with a stereomicroscope. Pulsed magnetic-field gradients encode displacements in the nuclear magnetic resonance signal, allowing motion normal to the barriers to be compared with less confined directions.[ref]
John Tanner, the return to Germelshausen, and the work that nearly never happened by Derek Jones. Printed p.6, mica-stack account and Figure3
This is a material confinement experiment, not a change in the liquid's chemistry.For an original one-dimensional model, let particles diffuse in with constant diffusion coefficient :The wall conditions are homogeneous Neumann boundary conditions for the diffusion equation. Start in equilibrium, so has the uniform distribution on . Its variance isAt long time separation, the positions decorrelate. Consequently the mean squared displacement approachesAt short times relative to , most displacements have not encountered a wall, and the leading scale is the free-diffusion result . For a hypothetical gap , the long-time plateau is and its root is about .
Equilibrium restricted diffusion saturates in a reflecting interval
. Original Neumann-eigenmode calculation with uniform initial positions; time is and mean squared displacement is divided by . Perfectly reflecting walls give a plateau , while free diffusion continues as . These are model curves, not measurements of the mica stack. Details.The initial condition matters: fixing instead gives the long-time value . Likewise, echo attenuation in a confined sample must not automatically be converted into a free-diffusion coefficient. Confinement geometry, timing and the selected boundary model all affect the inference.
