Analytical testcase for 2D/3D anisotropic Diffusion (Heat Kernel)

Question

I want to verify and compare different Discretizations of the anisotropic diffusion equation in 2D / 3D image of my testsetting.

I want to verify and compare different Discretizations of the anisotropic diffusion equation in 2D / 3D.

enter image description here

In order to both test the timestepping and the spatial discretisations I had a look at using the heat kernels as an analytical solution to diffusion equations. The one dimensional heat kernel (wiki) looks like this:

$$ u(x,t) = \frac{1}{(4\pi t)^{d/2}} e^{-x^2/4t}\ $$ and solves following diffusion equation of dimension $d$ analytically:

$$ \frac{\partial u}{\partial t} = -\Delta u $$ The idea is that for an analytical solution I can verify my convergence both spatially and temporal. It is important in my case, that the diffusion is anisotropic, and might not be aligned with the underlying grids.

My anisotropic diffusion equation reads: $$ \frac{\partial u}{\partial t} = \nabla (D \nabla u) $$ In order to expand the analytical solution to my case, where the diffusion is not scalar and unity, but a full (symmetric positive definite) tensor D, I had a look at multivariate gaußian distributions (wikipage).

My question is the following:

1. Is there somewhere in the literature a heat kernel for anisotropic (but homogeneous) diffusion which I can use to verify my numerics?

By analogy I would come up with something like: $$ u(\vec{x},\vec{\mu},t) = \frac{1}{(4\pi~ |\mathbf{D}| t)^{d/2}} \exp \large(\frac{-(\vec{x}-\vec{\mu})^{T} \mathbf{D}^{-1}(\vec{x}-\vec{\mu})}{4~|\mathbf{D}|~t}\large). $$

2. Is the above equation a proper analytical solution of the anisotropic diffusion equation?

I am not sure if the exponent $d/2$ is correct for the $|D|$. I am aware that the question might touch both physics, and numerics, but I thought this might be useful to other people trying to verify their codes.

Any hints or help is appreciated!

Answer 1

Thanks to @Andrew I dug out the following solution: $$ u(\vec{x},\vec{\mu},t) = \frac{1}{(4\pi~ t)^{d/2}~\sqrt{|\mathbf{D}|}} \exp \large(\frac{-(\vec{x}-\vec{\mu})^{T} \mathbf{D}^{-1}(\vec{x}-\vec{\mu})}{4~t}\large). $$ It seems the root of the determinant is not dependend on the problem dim. Take it with a grain of salt.