线性薛定谔方程实现界面推移

线性薛定谔方程实现界面推移

考虑如下线性薛定谔方程
$$r^2{\alpha }^2{\nabla }^{2}T-T=0$$式中$r$是界面推移的速率，单位$\mathrm{m}/\mathrm{s}$。$\alpha$是扩散项系数，单位为$\mathrm{s}$。$T$是无量纲温度，单位为1
令
$$T=e^{-\frac{\tau }{\alpha }}$$式中$\tau$是时间场。代入线性薛定谔方程中有
$$r^2{\alpha }^2{\nabla }^2{e}^{-\frac{\tau }{\alpha }}-e^{-\frac{\tau }{\alpha }}=0$$$$r^2{\alpha }^2\nabla \cdot \left(-\frac{1}{\alpha }{e}^{-\frac{\tau }{\alpha }}\nabla \tau \right)-e^{-\frac{\tau }{\alpha }}=0$$$$r^2{\alpha }^2\left(\frac{1}{{\alpha }^2}{e}^{-\frac{\tau }{\alpha }}{\left|\nabla \tau \right|}^2-\frac{1}{\alpha }{e}^{-\frac{\tau }{\alpha }}{\nabla }^2\tau \right)-e^{-\frac{\tau }{\alpha }}=0$$$$e^{-\frac{\tau }{\alpha }}\left(r^2{\left|\nabla \tau \right|}^2-r^2\alpha {\nabla }^2\tau -1\right)=0$$$${\left|\nabla \tau \right|}^2-\alpha {\nabla }^2\tau -\frac{1}{r^2}=0$$$${\left|\nabla \tau \right|}^2-\frac{1}{r^2}=\alpha {\nabla }^2\tau$$
可见求解了线性薛定谔方程就等价地求解了程函方程，也就获得了界面推移的结果。
$$\tau =-\alpha \ln T$$
参考自Numerical Solving a Boundary Value Problem for the Eikonal Equation

求解案例

一个正方形可燃烧区域，左侧为高燃速材料，右侧为低燃速材料。上边界、左边界和下边界被点燃，右边界不燃烧。
以下是燃速分布

以下是无量纲温度分布

以下是点燃时间场分布

此方法虽然简单（只用使用有限元法求解一个对称线性方程组），但是界面间的精度还需进一步验证

openFoam实现https://gitee.com/jedi-knight/sef-foam