【文献2005】快速步进法FMM求解固体火箭发动机装药燃面推移

THREE-DIMENSIONAL INTERNAL BALLISTIC ANALYSIS BY FAST MARCHING METHOD APPLIED TO PROPELLANT GRAIN BURN-BACK

作者 Toker K A, Tinaztepe H T, Aksel M H.

简介

使用FMM法燃面推移

In order to obtain the control volume with moving walls, Fast Marching Method on fixed unstructured grid is implemented.

使用界面捕获模型切割cells

The interface capturing modeling is maintained on a fixed tetrahedron grid where
the interfacial cells are sliced at each time interval, at the interface and the mass injection boundary condition is applied to these cut cell surfaces.

使用FVM法计算流场（Euler方程）

The flow solutions are carried out by a three-dimensional finite volume Euler solver and results are compared with test data and previously validated algorithms.

Only steady Euler flow solutions with moving boundaries are presented in this article.

移动界面问题有两种处理办法：界面跟踪和界面捕获

The interface tracking methods on the other hand need elaborate mesh refinement techniques in order to overcome the fact that grid elements will move into each other.

FMM法

程函方程的一阶迎风格式如下

FMM的步骤如下：

1. 根据Known的T值和上式计算Close的T值
2. 选择Close中T最小者为Trial
3. 把Trial的邻居变为新的Close
4. 把Trial变为Known
5. 回到第1步
   这样Known的节点越来越多，Close的节点越来越少，实现燃面推移
   

作者认为燃面推移的有限差分法的复杂度为$O(N^4)$而快速步进法的复杂度为$O(N^{3}logN)$,此结论有待进一步考证

作者还讨论了在非结构化四面体网格中FMM的算法。。。

单元切割

The aim of cut-cell methodology is to create the control volume in which the flow solutions will take place.Only the elements that carry the interface inside their volume are affected by this procedure. In this case, the grid is sliced into two domains, where flow solutions are carried out in one of them.

In order to obtain three-dimensional mesh for the flow solver, the mesh is sliced through the interface. The mesh elements are separated into 3 groups:
i. Mesh elements, which are inside the propellant side.
ii. Mesh elements, which are inside the gas side.
iii. Mesh elements, which are sliced by the interface.

被界面切割的单元(sliced by the interface)有可能变为非四面体单元。但是没有关系，混合网格求解器依然可以完成求解。

If the flow solver is capable of handling hybrid mesh, that is to say tetrahedron, pyramid, wedge and hexahedron elements, then the sliced mesh can be used as it is. In this work, the flow solver is adapted so that hybrid grid may be used and solutions are given only for the hybrid solver.

结果

作者取了4个点来计算CFD，没有考虑侵蚀燃烧

作者还考虑了Hump Effect，即由于组分迁移，粘合剂在金属附近聚集，燃速在金属壳体附近降低。作者给出了经验公式。

因为使用一阶迎风格式，误差较大，精度取决于网格细度，燃面的不连续性明显，有必要进行磨光。

The solutions are simply dependent on mesh intensity because the definition of the surface is discontinuous. As the propagation extends a large distance, the discontinuities become visible. An intermediate smoothing may be used; however, this is against the nature of the solution method.