Adaptive unstructured finite element method for two-dimensional detonation simulations

Abstract

A non-equilibrium reacting flow methodology has been added to a conservative, monotonic, compressible flow solver to allow numerical simulations of gas detonations. This flow solver incorporates unstructured dynamically adaptive meshes with the Finite Element Method – Flux Corrected Transport (FEM-FCT) scheme, which has shown excellent predictive capability of various non-reacting compressible flows. A two-step induction parameter model was used to model the combustion of the gas phase coupled with an energy release equation which was simulated with a point implicit finite element scheme. This combustion model was then applied to a two-dimensional detonation test case of a hypothetical fuel:oxygen mixture. The detonation simulation yielded two transverse waves which continued to propagate. This feature and the detonation shock speed mean and fluctuations were found to be grid-independent based on a resolution of about twenty elements within the average induction length. The resolution was efficiently achieved with the unstructured dynamically adaptive finite elements, which were three orders of magnitude less in number then required for uniform discretization.