2D/3D 各向异性扩散的分析测试用例(热核)

计算科学 收敛 扩散 传播热量
2021-12-11 08:45:10

我想验证和比较 2D / 3D 中各向异性扩散方程的不同离散化。

在此处输入图像描述

为了测试时间步长和空间离散,我研究了使用热核作为扩散方程的解析解。一维热核(wiki-link)如下所示:

u(x,t)=1(4πt)d/2ex2/4t 
并求解以下维数扩散方程d分析:
ut=Δu

这个想法是,对于分析解决方案,我可以在空间和时间上验证我的收敛性。就我而言,重要的是扩散是各向异性的,并且可能与底层网格不对齐。

我的各向异性扩散方程如下:

ut=(Du)

为了将分析解决方案扩展到我的情况,其中扩散不是标量和统一,而是完整(对称正定)张量 D,我查看了多元高斯分布(wiki-link)。

我的问题如下:

  1. 文献中是否有用于各向异性(但均匀)扩散的热核,我可以用它来验证我的数值?

以此类推,我会想出类似的东西:

u(x,μ,t)=1(4π |D|t)d/2exp((xμ)TD1(xμ)4 |D| t).

  1. 上述方程是各向异性扩散方程的适当解析解吗?

我不确定指数是否d/2是正确的|D|. 我知道这个问题可能涉及物理和数字,但我认为这可能对其他试图验证他们的代码的人有用。

任何提示或帮助表示赞赏!

1个回答

由于拉普拉斯算子是一个椭圆算子,您正在寻找假设的常数扩散矩阵的 Cholesky 分解D

D=LLT

因此抛物线方程可以写成:

tu=T(Du)=(LT)T(LTu)=~T(~u)=~2u

命名~=LT因此x~=L1x. 最后一个 PDE 完全解耦(使用傅里叶变换和变量分离)。

非标准化热核的解决方案将是(使用x~=L1x):

Φ=n=1dexp(x~n2/(4t))4πt=exp(x~Tx~/(4t))(4πt)d/2=exp(xTD1x/(4t))(4πt)d/2