Parallel chemistry acceleration algorithm with ISAT table-size control in the application of gaseous detonation

Abstract

In order to improve the computational efficiency of a parallel ISAT (in situ adaptive tabulation)-based chemistry acceleration algorithm in the computations of transient, chemically reacting flows, a control strategy is proposed to maintain the sizes of the data tables in the ISAT computations. The table-size control strategy is then combined with a parallel algorithm to simulate two-dimensional gaseous detonation wave propagation. In the computation of 2H2 + O2 detonation, two sets of tests are conducted to identify the size control strategy. In the first set, the maximum total table size (Mtot) summed over all sub-zones is fixed, while the maximum size of the table on each sub-zone (Msin) is varied. In the second set, a fixed Msin is used for all the tables on the sub-zones while Mtot is varied. A maximum speedup ratio of 4.29 is found in the former tests, while 5.52 is found in the latter. Two parameters, σf and p, are proposed to analyze the load balance and synchronization among table operations in the parallel ISAT computations in the above tests. It is found that both load balance and synchronization have clear influences on the speedup ratio. A parameter pM is defined, and a strategy to choose the optimal maximum table sizes (both Mtot and Msin) based on pM is proposed and is verified to be universal in the computations of both 2H2 + O2 detonation and C2H4 + 3O2 detonation. Finally, the parallel acceleration algorithm enhanced with table-size control is shown to be highly accurate for the detonations in both fuels.