使用PEF相场的方法实现不规则形状燃面推移（观点）

不规则形状的定义

对于$n$维几何图形，可以使用$n+1$维的某个函数的等值线定义。这里称为相场$\phi$。

定义$\phi$取1的位置是固体装药，$\phi$取-1的位置是气体。那么气-固边界自然而然地就定义好了（虽然没有解析表达式，但是边界是可以绘制出来的）

不规则形状燃面推移

根据相场$\phi$，提出两个变量：虚拟燃速$r$和源项因子$isSolid$。它们的定义为
$$r=\{\begin{array}{c}1,\text{如果}\phi >0\\ 1+5000{\phi }^2,\text{其他情况}\end{array}$$$$isSolid=\{\begin{array}{c}1,\text{如果}\phi >0\\ 0,\text{其他情况}\end{array}$$
之所以虚拟燃速使用二次罚函数的形式，是为了使得在有限元框架下交界面更加光滑连续。
那么PEF求解的控制方程为
$$\alpha r{\nabla }^{2}W=\left(r^2\nabla W\cdot \nabla W-1\right)\cdot isSolid$$

一个实例

相场$\phi$如下

燃面推移结果$W$场如下

有限元法燃面推移的收敛性和收敛速度（因为使用一阶单元）有明显提升，且气-固边界清晰光滑。

进一步的性能还有待考察。