Publications

Abstract not provided.

This work will be presented at the combustion symposium in Milan, Italy.

Abstract not provided.

The ability to accurately predict the structure and dynamics of pool fires using computational simulations is of great interest in a wide variety of applications, including accidental and wildland fires. However, the presence of physical processes spanning a broad range of spatial and temporal scales poses a significant challenge for simulations of such fires, particularly at conditions near the transition between laminar and turbulent flow. In this study, we examine the transition to turbulence in methane pool fires using high-resolution simulations with multi-step finite rate chemistry, where adaptive mesh refinement (AMR) is used to directly resolve small-scale flow phenomena. We perform three simulations of methane pool fires, each with increasing diameter, corresponding to increasing inlet Reynolds and Richardson numbers. As the diameter increases, the flow transitions from organized vortex roll-up via the puffing instability to much more chaotic mixing associated with finger formation along the shear layer and core collapse near the inlet. These effects combine to create additional mixing close to the inlet, thereby enhancing fuel consumption and causing more rapid acceleration of the fluid above the pool. We also make comparisons between the transition to turbulence and core collapse in the present pool fires and in inert helium plumes, which are often used as surrogates for the study of buoyant reacting flows.

Tabulated chemistry models are widely used to simulate large-scale turbulent fires in applications including energy generation and fire safety. Tabulation via piecewise Cartesian interpolation suffers from the curse-of-dimensionality, leading to a prohibitive exponential growth in parameters and memory usage as more dimensions are considered. Artificial neural networks (ANNs) have attracted attention for constructing surrogates for chemistry models due to their ability to perform high-dimensional approximation. However, due to well-known pathologies regarding the realization of suboptimal local minima during training, in practice they do not converge and provide unreliable accuracy. Partition of unity networks (POUnets) are a recently introduced family of ANNs which preserve notions of convergence while performing high-dimensional approximation, discovering a mesh-free partition of space which may be used to perform optimal polynomial approximation. We assess their performance with respect to accuracy and model complexity in reconstructing unstructured flamelet data representative of nonadiabatic pool fire models. Our results show that POUnets can provide the desirable accuracy of classical spline-based interpolants with the low memory footprint of traditional ANNs while converging faster to significantly lower errors than ANNs. For example, we observe POUnets obtaining target accuracies in two dimensions with 40 to 50 times less memory and roughly double the compression in three dimensions. We also address the practical matter of efficiently training accurate POUnets by studying convergence over key hyperparameters, the impact of partition/basis formulation, and the sensitivity to initialization.

Numerical simulations were performed in 3D Cartesian coordinates to examine the post-detonation processes produced by the detonation of a 12 mm-diameter hemispherical PETN explosive charge in air. The simulations captured air dissociation by the Mach 20+ shock, chemical equilibration, and afterburning using finite-rate chemical kinetics with a skeletal chemical reaction mechanism. The Becker-Kistiakowsky-Wilson real-gas equation of state is used for the gas-phase. A simplified programmed burn model is used to seamlessly couple the detonation propagation through the explosive charge to the post-detonation reaction processes inside the fireball. Four charge sizes were considered, including diameters of 12 mm, 38 mm, 120 mm, and 1200 mm. The computed blast, shock structures, and chemical composition within the fireball agree with literature. The evolution of the flow at early times is shown to be gas dynamic driven and nearly self-similar when the time and space was scaled. The flow fields were azimuthally averaged and a mixing layer analysis was performed. The results show differences in the temperature and chemical composition with increasing charge size, implying a transition from a chemical kinetic-limited to a mixing-limited regime.

As large systems of Li-ion batteries are being increasingly deployed, the safety of such systems must be assessed. Due to the high cost of testing large systems, it is important to extract key safety information from any available experiments. Developing validated predictive models that can be exercised at larger scales offers an opportunity to augment experimental data In this work, experiments were conducted on packs of three Li-ion pouch cells with different heating rates and states of charge (SOC) to assess the propagation behavior of a module undergoing thermal runaway. The variable heating rates represent slow or fast heating that a module may experience in a system. As the SOC decreases, propagation slows down and eventually becomes mitigated. It was found that the SOC boundary between propagation and mitigation was higher at a heating rate of 50 °C/min than at 10 °C/min for these cells. However, due to increased pre-heating at the lower heating rate, the propagation speed increased. Simulations were conducted with a new intra-particle diffusion-limited reaction model for a range of anode particle sizes. Propagation speeds and onset times were generally well predicted, and the variability in the propagation/mitigation boundary highlighted the need for greater uncertainty quantification of the predictions.

This report documents the generation of a mechanism to predict the inclusion of carbon soot particles in a high explosive flow. The mechanism includes gasification and oxidation reactions, formation, sublimation, radiation, and agglomeration. Each part of the mechanism is derived from properties in the literature. The influence of each part of the mechanism is explored using simple, example simulations consisting of a 12 mm diameter 2,4,6-Trinitrotoluene charge detonated in ambient air. The mechanism has not been quantitatively compared to experiments. Additional efforts will be required to tune and validate it, which will require continued advancements in experimental diagnostics and simulation techniques.

This report documents the generation of a skeletal chemical reaction mechanism for use with hemispherical pentaerythritol tetranitrate charges. Skeletal mechanisms can substantially reduce computation time while maintaining accuracy. The methodology within uses faster running sample simulations to build a representative thermodynamic state space. These thermodynamic states are used with a constant-volume reactor analysis and a reaction flow analysis to remove unimportant species and reactions from a full chemical reaction mechanism. For the given test case, this results in a 6x speedup in computation time for directly comparable simulations in 2D axisymmetric simulations. We see a 30x speedup in simulations in 3D Cartesian coordainates when compared to a prior full kinetics simulation. There is strong agreement between temperature and species mass fraction profiles between the full and skeletal chemical reaction mechanisms. These methodologies can be applied to any explosive, given the availability of sample simulations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A numerical simulation study was performed to examine the post-detonation reaction processes produced by the detonation of a 12 mm diameter hemispherical pentaerythritol tetranitrate (PETN) explosive charge. The simulations used a finite rate detailed chemical reaction model consisting of 59 species and 368 reactions to capture post-detonation reaction processes including air dissociation from Mach 19+ shock waves that initially break out of the PETN charge, reactions within the detonation products during expansion, and afterburning when the detonation products mix with the shock heated air. The multi-species and thermodynamically complete Becker-Kistiakowsky-Wilson real-gas equation of state is used for the gaseous phase to allow for the mixing of reactive species. A recent simplified reactive burn model is used to propagate the detonation through the charge and allow for detailed post-detonation reaction processes. The computed blast, shock structures, and mole fractions of species within the detonation products agree well with experimental measurements. A comparison of the simulation results to equilibrium calculations indicates that the assumption of a local equilibrium is fairly accurate until the detonation products rapidly cool to temperatures in the range of 1500-1900 K by expansion waves. Below this range, the computed results show mole fractions that are nearly chemically frozen within the detonation products for a significant portion of expansion. These results are consistent with the freeze out approximation used in the blast modeling community.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Chemistry tabulation is a common approach in practical simulations of turbulent combustion at engineering scales. Linear interpolants have traditionally been used for accessing precomputed multidimensional tables but suffer from large memory requirements and discontinuous derivatives. Higher-degree interpolants address some of these restrictions but are similarly limited to relatively low-dimensional tabulation. Artificial neural networks (ANNs) can be used to overcome these limitations but cannot guarantee the same accuracy as interpolants and introduce challenges in reproducibility and reliable training. These challenges are enhanced as the physics complexity to be represented within the tabulation increases. In this manuscript, we assess the efficiency, accuracy, and memory requirements of Lagrange polynomials, tensor product B-splines, and ANNs as tabulation strategies. We analyze results in the context of nonadiabatic flamelet modeling where higher dimension counts are necessary. While ANNs do not require structuring of data, providing benefits for complex physics representation, interpolation approaches often rely on some structuring of the table. Interpolation using structured table inputs that are not directly related to the variables transported in a simulation can incur additional query costs. This is demonstrated in the present implementation of heat losses. We show that ANNs, despite being difficult to train and reproduce, can be advantageous for high-dimensional, unstructured datasets relevant to nonadiabatic flamelet models. We also demonstrate that Lagrange polynomials show significant speedup for similar accuracy compared to B-splines.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.