Search NASA⌕ Search

Engineering topics

Barwey, Shivam

Publications and source records attributed to Barwey, Shivam.

Jacobian-scaled K-means clustering for physics-informed segmentation of reacting flows

This work introduces Jacobian-scaled K-means (JSK-means) clustering, which is a physicsinformed clustering strategy centered on the K-means framework. The method allows for the injection of underlying physical knowledge into the clustering procedure through a distance function modification: instead of leveraging conventional Euclidean distance vectors, the JSKmeans procedure operates on distance vectors scaled by matrices obtained from dynamical system Jacobians evaluated at the cluster centroids. The goal of this work is to show how the JSKmeans algorithm - without modifying the input dataset - produces clusters that capture regions of dynamical similarity, in that the clusters are redistributed towards high-sensitivity regions in phase space and are described by similarity in the source terms of samples instead of the samples themselves. The algorithm is demonstrated on a complex reacting flow simulation dataset (a channel detonation configuration), where the dynamics in the thermochemical composition space are known through the highly nonlinear and stiff Arrhenius-based chemical source terms. Interpretations of cluster partitions in both physical space and composition space reveal how JSK-means shifts clusters produced by standard K-means towards regions of high chemical sensitivity (e.g., towards regions of peak heat release rate near the detonation reaction zone). Furthermore, the findings presented here illustrate the benefits of utilizing Jacobian-scaled distances in clustering techniques, and the JSK-means method in particular displays promising potential for improving former partition-based modeling strategies in reacting flow (and other multi-physics) applications.

Clustering↗

Detonation structure in the presence of mixture stratification using reaction-resolved simulations

Many investigations of detonation-based combustors have identified reactant mixture inhomogeneity as having a leading-order impact on wave dynamics and combustion efficiency. To examine this phenomenon in a simplified context, an array of two- and three-dimensional channel detonation simulations are conducted in the present work. The reactant mixture consists of stratified fuel and air, wherein the randomly distributed equivalence ratio field features a characteristic stratification length scale. Detailed chemical kinetics are implemented in an adaptive mesh refinement solution framework where the region near the shock front is resolved with Ο (100) cells per representative ZND induction length. The results show that in comparison to baseline cases with uniform reactant mixtures, reactant stratification has a marked impact on the detonation structure. Increasing the stratification length scale increases the size and irregularity of the detonation cells, yielding larger variations in wave speed. Triple point collisions in fuel-rich regions lead to local wave speeds above the notional mean CJ speed, but wave passage through inert regions causes the local wave speed and strength to diminish. Further, conditional statistics show that increasing the stratification length scale increases the variance in pressure and temperature in the primary reaction zone, as well as the variance in heat release over a range of mixture conditions. In addition to the reactant mixture, the impact of the boundary condition behind the detonation is also investigated. The results show that an inflow boundary condition acts to over-drive the wave, leading to higher peak pressures, smaller detonation cells, and increased reactant consumption. On the other hand, cases with a wall behind the wave exhibit weaker waves with lower peak pressures and heat release rates, as well as greater variance in conditional quantities. Comparisons between complementary two- and three-dimensional simulations show reasonable qualitative agreement in wave structure, speed, and conditional statistics.

42 ENGINEERING↗

Multiscale graph neural network autoencoders for interpretable scientific machine learning

The goal of this work is to address two limitations in autoencoder-based models: latent space interpretability and compatibility with unstructured meshes. This is accomplished here with the development of a novel graph neural network (GNN) autoencoding architecture with demonstrations on complex fluid flow applications. To address the first goal of interpretability, the GNN autoencoder achieves reduction in the number nodes in the encoding stage through an adaptive graph reduction procedure. Further, this reduction procedure essentially amounts to flowfieldconditioned node sampling and sensor identification, and produces interpretable latent graph representations tailored to the flowfield reconstruction task in the form of so-called masked fields. These masked fields allow the user to (a) visualize where in physical space a given latent graph is active, and (b) interpret the time-evolution of the latent graph connectivity in accordance with the time-evolution of unsteady flow features (e.g. recirculation zones, shear layers) in the domain. To address the goal of unstructured mesh compatibility, the autoencoding architecture utilizes a series of multi-scale message passing (MMP) layers, each of which models information exchange among node neighborhoods at various lengthscales. The MMP layer, which augments standard single-scale message passing with learnable coarsening operations, allows the decoder to more efficiently reconstruct the flowfield from the identified regions in the masked fields. Analysis of latent graphs produced by the autoencoder for various model settings are conducted using unstructured snapshot data sourced from large-eddy simulations in a backward-facing step (BFS) flow configuration with an OpenFOAM-based flow solver at high Reynolds numbers.

97 MATHEMATICS AND COMPUTING↗