Comparison of Finite Elements and Meshless Methods for Modeling of Solid Phase Chemical Reactions under Explosive Loading

Abstract

The paper represents a comparison of the results of modeling the reaction of solid-phase synthesis in a mixture of aluminum and sulfur under impact loading by meshless and finite element methods. The model of solid-phase chemical reaction of aluminum with sulfur with zero order kinetics and pressure-dependent acceleration is considered. It is established that the divergent form of the equations of motion of elastic plastic medium with solid-phase reactions is poorly suited for SPH calculations. It has been established that the weak variation formulation of SPH is more preferable due to its greater stability, conservative properties, and stability of calculations. In numerical calculations, good agreement with experiments on explosion initiation of solid-phase chemical reactions was obtained. It has been shown that the zero-order pressure-accelerated reaction model qualitatively describes the process of explosive initiation and simulates the formation of zones of high and low reaction degrees. It was also found that this model is very sensitive to the value of the coefficient describing the reaction acceleration at high pressure.