Numerical model of the Moving Burning Surface in High-energy Systems

Experimental investigation of processes occurring in high-energy systems is complicated and expensive. High temperatures and pressures, transience, and nascent toxic gases make mathematical tools demanding in this domain. Hence, using numerical models is a key component for an efficient design of high-energy systems.
In this study, we consider the coupled problem of tracking the burning surface together with calculating flow parameters in high-energy systems. For this purpose, we develop the numerical model that takes into account both complex grain geometry and multidimensional gas flow pattern. We parallelize this algorithm using CUDA technology. It allows us to gain up to 30x acceleration in comparison with consequent implementation.
We test our numerical model on the configuration without nozzle, that is when the nozzle section is made of propellant as well. We prove approximation convergence for the problem and depict the evolution of the propellant burning surface. Also, we show dependencies of thrust, burning surface area, the average speed on the nozzle edge, and maximum pressure on time during the firing. In the end, we present distributions of flow parameters at initial, intermediate, and final time points.