基于等值面(Level Set)函数方法的复杂装药燃面算法研究

通过使用等值面（Level Set）函数界面追踪方法，将固体火箭发动机燃面看作不同物质界面，建立了新的燃面算法，通过追踪固体装药燃烧界面计算出燃面变化规律。采用高精度Weighted—ENO和TVD Runge—Kutta进行数值求解，该方法可以计算复杂几何构形的发动机装药、含缺陷装药和变燃速装药的构型变化和燃烧面积，具有良好的通用性和准确性。新的燃面算法具有了界面追踪的能力，为固体火箭发动机内弹道计算以及装药燃烧与发动机内流场数值模拟耦合计算提供了良好的接口条件。

A New Algorithm for Burning Surface Calculation of Solid Rocket Motor with Complicated Grain Based on Level Set Method

We apply level set method 1,2] to create a new burning surface calculation algorithm, viewing burning surface as interface of different materials. New algorithm attains burning surface variation regulation of arbitrary configuration grain by tracing grain burning surface recession. In the full paper, we explain our new algorithm in detail, here we give only a briefing. Our explanation centers around three topics: (1) calculation of burning surface area of complicated 3D grain; (2) calculation of burning surface area of grain with flaw; (3) calculation of burning surface area of grain with changing burning rate. We used explicit weighted essentially nonoscillatory(WENO) discretization scheme and total variation diminishing(TVD) Runge-Kutta scheme in numerical simulation. Our numerical simulation used grains of assumed configurations, not grains used in actual SRM(solid rocket motor). Simulation results show preliminarily that this new method has high calculation precision and good adaptability. This paper can be used to couple grain burning surface recession process to the numerical simulation of SRM chamber flowfield.