使用COMSOL提取指定时刻固体火箭发动机内的流体域（观点）

为何要提取流体域？

侵蚀燃烧影响下的复杂三维装药准静态内弹道性能预示的迭代法https://blog.csdn.net/weixin_43325228/article/details/123142114
为了完成含有侵蚀燃烧的内弹道性能预示，必须使用CFD方法计算燃面推移过程中的流场。
但是随着燃面推移，流体域形状不断变化，如何提取指定时刻流体域是一个难点问题

COMSOL要点

1 导入x_t几何体

几何体由初始流体域（灰色）和初始装药（紫色）构成，右上方是喷管的收敛段，未建模喷管的扩张段

2 定义isSolid变量

下图是初始流体域的变量，isSolid=0。对于初始固体域，isSolid=1

3 定义PEF法的源项和扩散项，设置边界条件

PEF法可参考https://blog.csdn.net/weixin_43325228/article/details/128478544

4 对整体划分网格

5 开展PEF法的计算，得到燃面推移结果

PEF法可参考https://blog.csdn.net/weixin_43325228/article/details/128478544
由于我们在整个求解域（初始固体域和初始流体域）上开展计算，所以得到的解也是全局分布的。
初始流体域的$T$场恒为0，初始固体域的$T$场符合燃面推移规律

6 使用过滤器得到实际流体区域

以下是2s时刻的结果

7 将该流体区域导出为STL文件

如果明确质量通量边界条件后，后续可以使用COMSOL、Fluent或OpenFOAM完成CFD计算