-
Symbolic construction of the chemical Jacobian of quasi-steady state (QSS) chemistries for Exascale computing platforms
Authors:
Malik Hassanaly,
Nicholas T. Wimer,
Anne Felden,
Lucas Esclapez,
Julia Ream,
Marc T. Henry de Frahan,
Jon Rood,
Marc Day
Abstract:
The Quasi-Steady State Approximation (QSSA) can be an effective tool for reducing the size and stiffness of chemical mechanisms for implementation in computational reacting flow solvers. However, for many applications, stiffness remains, and the resulting model requires implicit methods for efficient time integration. In this paper, we outline an approach to formulating the QSSA reduction that is…
▽ More
The Quasi-Steady State Approximation (QSSA) can be an effective tool for reducing the size and stiffness of chemical mechanisms for implementation in computational reacting flow solvers. However, for many applications, stiffness remains, and the resulting model requires implicit methods for efficient time integration. In this paper, we outline an approach to formulating the QSSA reduction that is coupled with a strategy to generate C++ source code to evaluate the net species production rate, and the chemical Jacobian. The code-generation component employs a symbolic approach enabling a simple and effective strategy to analytically compute the chemical Jacobian. For computational tractability, the symbolic approach needs to be paired with common subexpression elimination which can negatively affect memory usage. Several solutions are outlined and successfully tested on a 3D multipulse ignition problem, thus allowing portable application across a chemical model sizes and GPU capabilities. The implementation of the proposed method is available at https://github.com/AMReX-Combustion/PelePhysics under an open-source license.
△ Less
Submitted 30 April, 2024;
originally announced May 2024.
-
SUNDIALS Time Integrators for Exascale Applications with Many Independent ODE Systems
Authors:
Cody J. Balos,
Marc Day,
Lucas Esclapez,
Anne M. Felden,
David J. Gardner,
Malik Hassanaly,
Daniel R. Reynolds,
Jon Rood,
Jean M. Sexton,
Nicholas T. Wimer,
Carol S. Woodward
Abstract:
Many complex systems can be accurately modeled as a set of coupled time-dependent partial differential equations (PDEs). However, solving such equations can be prohibitively expensive, easily taxing the world's largest supercomputers. One pragmatic strategy for attacking such problems is to split the PDEs into components that can more easily be solved in isolation. This operator splitting approach…
▽ More
Many complex systems can be accurately modeled as a set of coupled time-dependent partial differential equations (PDEs). However, solving such equations can be prohibitively expensive, easily taxing the world's largest supercomputers. One pragmatic strategy for attacking such problems is to split the PDEs into components that can more easily be solved in isolation. This operator splitting approach is used ubiquitously across scientific domains, and in many cases leads to a set of ordinary differential equations (ODEs) that need to be solved as part of a larger "outer-loop" time-stepping approach. The SUNDIALS library provides a plethora of robust time integration algorithms for solving ODEs, and the U.S. Department of Energy Exascale Computing Project (ECP) has supported its extension to applications on exascale-capable computing hardware. In this paper, we highlight some SUNDIALS capabilities and its deployment in combustion and cosmology application codes (Pele and Nyx, respectively) where operator splitting gives rise to numerous, small ODE systems that must be solved concurrently.
△ Less
Submitted 2 May, 2024;
originally announced May 2024.
-
Inverse Modeling of Bubble Size Dynamics for Interphase Mass Transfer and Gas Holdup in CO2 Bubble Column Reactors
Authors:
Malik Hassanaly,
John M. Parra-Alvarez,
Mohammad J. Rahimi,
Hariswaran Sitaraman
Abstract:
The use of microbial gas fermentation for transforming captured CO2 into sustainable fuels and chemicals has been identified as a promising decarbonization pathway. To accelerate the scale-up of gaseous CO2 fermentation reactors, computational models need to predict gas-to-liquid mass transfer which requires capturing the bubble size dynamics, i.e. bubble breakup and coalescence. In this work, an…
▽ More
The use of microbial gas fermentation for transforming captured CO2 into sustainable fuels and chemicals has been identified as a promising decarbonization pathway. To accelerate the scale-up of gaseous CO2 fermentation reactors, computational models need to predict gas-to-liquid mass transfer which requires capturing the bubble size dynamics, i.e. bubble breakup and coalescence. In this work, an inverse modeling approach is used to calibrate the breakup and coalescence closure models, that are used in the Multiple-Size-Group (MUSIG) population balance modeling (PBM). The calibration aims at replicating experimental results obtained in a CO2-air-water-coflowing bubble column reactor. Bayesian inference is used to account for noise in the experimental dataset and bias in the simulation results. The estimated simulation bias also allows identifying the best-performing closure models irrespective of the model parameters used. The calibration results suggest that the breakage rate is underestimated by one order of magnitude in two different breakup modeling approaches.
△ Less
Submitted 30 April, 2024;
originally announced April 2024.
-
A Priori Uncertainty Quantification of Reacting Turbulence Closure Models using Bayesian Neural Networks
Authors:
Graham Pash,
Malik Hassanaly,
Shashank Yellapantula
Abstract:
While many physics-based closure model forms have been posited for the sub-filter scale (SFS) in large eddy simulation (LES), vast amounts of data available from direct numerical simulation (DNS) create opportunities to leverage data-driven modeling techniques. Albeit flexible, data-driven models still depend on the dataset and the functional form of the model chosen. Increased adoption of such mo…
▽ More
While many physics-based closure model forms have been posited for the sub-filter scale (SFS) in large eddy simulation (LES), vast amounts of data available from direct numerical simulation (DNS) create opportunities to leverage data-driven modeling techniques. Albeit flexible, data-driven models still depend on the dataset and the functional form of the model chosen. Increased adoption of such models requires reliable uncertainty estimates both in the data-informed and out-of-distribution regimes. In this work, we employ Bayesian neural networks (BNNs) to capture both epistemic and aleatoric uncertainties in a reacting flow model. In particular, we model the filtered progress variable scalar dissipation rate which plays a key role in the dynamics of turbulent premixed flames. We demonstrate that BNN models can provide unique insights about the structure of uncertainty of the data-driven closure models. We also propose a method for the incorporation of out-of-distribution information in a BNN. The efficacy of the model is demonstrated by a priori evaluation on a dataset consisting of a variety of flame conditions and fuels.
△ Less
Submitted 28 February, 2024;
originally announced February 2024.
-
PINN surrogate of Li-ion battery models for parameter inference. Part II: Regularization and application of the pseudo-2D model
Authors:
Malik Hassanaly,
Peter J. Weddle,
Ryan N. King,
Subhayan De,
Alireza Doostan,
Corey R. Randall,
Eric J. Dufek,
Andrew M. Colclasure,
Kandler Smith
Abstract:
Bayesian parameter inference is useful to improve Li-ion battery diagnostics and can help formulate battery aging models. However, it is computationally intensive and cannot be easily repeated for multiple cycles, multiple operating conditions, or multiple replicate cells. To reduce the computational cost of Bayesian calibration, numerical solvers for physics-based models can be replaced with fast…
▽ More
Bayesian parameter inference is useful to improve Li-ion battery diagnostics and can help formulate battery aging models. However, it is computationally intensive and cannot be easily repeated for multiple cycles, multiple operating conditions, or multiple replicate cells. To reduce the computational cost of Bayesian calibration, numerical solvers for physics-based models can be replaced with faster surrogates. A physics-informed neural network (PINN) is developed as a surrogate for the pseudo-2D (P2D) battery model calibration. For the P2D surrogate, additional training regularization was needed as compared to the PINN single-particle model (SPM) developed in Part I. Both the PINN SPM and P2D surrogate models are exercised for parameter inference and compared to data obtained from a direct numerical solution of the governing equations. A parameter inference study highlights the ability to use these PINNs to calibrate scaling parameters for the cathode Li diffusion and the anode exchange current density. By realizing computational speed-ups of 2250x for the P2D model, as compared to using standard integrating methods, the PINN surrogates enable rapid state-of-health diagnostics. In the low-data availability scenario, the testing error was estimated to 2mV for the SPM surrogate and 10mV for the P2D surrogate which could be mitigated with additional data.
△ Less
Submitted 26 March, 2024; v1 submitted 28 December, 2023;
originally announced December 2023.
-
PINN surrogate of Li-ion battery models for parameter inference. Part I: Implementation and multi-fidelity hierarchies for the single-particle model
Authors:
Malik Hassanaly,
Peter J. Weddle,
Ryan N. King,
Subhayan De,
Alireza Doostan,
Corey R. Randall,
Eric J. Dufek,
Andrew M. Colclasure,
Kandler Smith
Abstract:
To plan and optimize energy storage demands that account for Li-ion battery aging dynamics, techniques need to be developed to diagnose battery internal states accurately and rapidly. This study seeks to reduce the computational resources needed to determine a battery's internal states by replacing physics-based Li-ion battery models -- such as the single-particle model (SPM) and the pseudo-2D (P2…
▽ More
To plan and optimize energy storage demands that account for Li-ion battery aging dynamics, techniques need to be developed to diagnose battery internal states accurately and rapidly. This study seeks to reduce the computational resources needed to determine a battery's internal states by replacing physics-based Li-ion battery models -- such as the single-particle model (SPM) and the pseudo-2D (P2D) model -- with a physics-informed neural network (PINN) surrogate. The surrogate model makes high-throughput techniques, such as Bayesian calibration, tractable to determine battery internal parameters from voltage responses. This manuscript is the first of a two-part series that introduces PINN surrogates of Li-ion battery models for parameter inference (i.e., state-of-health diagnostics). In this first part, a method is presented for constructing a PINN surrogate of the SPM. A multi-fidelity hierarchical training, where several neural nets are trained with multiple physics-loss fidelities is shown to significantly improve the surrogate accuracy when only training on the governing equation residuals. The implementation is made available in a companion repository (https://github.com/NREL/pinnstripes). The techniques used to develop a PINN surrogate of the SPM are extended in Part II for the PINN surrogate for the P2D battery model, and explore the Bayesian calibration capabilities of both surrogates.
△ Less
Submitted 26 March, 2024; v1 submitted 28 December, 2023;
originally announced December 2023.
-
Ensemble flow reconstruction in the atmospheric boundary layer from spatially limited measurements through latent diffusion models
Authors:
Alex Rybchuk,
Malik Hassanaly,
Nicholas Hamilton,
Paula Doubrawa,
Mitchell J. Fulton,
Luis A. Martínez-Tossas
Abstract:
Due to costs and practical constraints, field campaigns in the atmospheric boundary layer typically only measure a fraction of the atmospheric volume of interest. Machine learning techniques have previously successfully reconstructed unobserved regions of flow in canonical fluid mechanics problems and two-dimensional geophysical flows, but these techniques have not yet been demonstrated in the thr…
▽ More
Due to costs and practical constraints, field campaigns in the atmospheric boundary layer typically only measure a fraction of the atmospheric volume of interest. Machine learning techniques have previously successfully reconstructed unobserved regions of flow in canonical fluid mechanics problems and two-dimensional geophysical flows, but these techniques have not yet been demonstrated in the three-dimensional atmospheric boundary layer. Here, we conduct a numerical analogue of a field campaign with spatially limited measurements using large-eddy simulation. We pose flow reconstruction as an inpainting problem, and reconstruct realistic samples of turbulent, three-dimensional flow with the use of a latent diffusion model. The diffusion model generates physically plausible turbulent structures on larger spatial scales, even when input observations cover less than 1% of the volume. Through a combination of qualitative visualization and quantitative assessment, we demonstrate that the diffusion model generates meaningfully diverse samples when conditioned on just one observation. These samples successfully serve as initial conditions for a large-eddy simulation code. We find that diffusion models show promise and potential for other applications for other turbulent flow reconstruction problems.
△ Less
Submitted 11 December, 2023; v1 submitted 1 March, 2023;
originally announced March 2023.
-
Bi-fidelity Modeling of Uncertain and Partially Unknown Systems using DeepONets
Authors:
Subhayan De,
Matthew Reynolds,
Malik Hassanaly,
Ryan N. King,
Alireza Doostan
Abstract:
Recent advances in modeling large-scale complex physical systems have shifted research focuses towards data-driven techniques. However, generating datasets by simulating complex systems can require significant computational resources. Similarly, acquiring experimental datasets can prove difficult as well. For these systems, often computationally inexpensive, but in general inaccurate, models, know…
▽ More
Recent advances in modeling large-scale complex physical systems have shifted research focuses towards data-driven techniques. However, generating datasets by simulating complex systems can require significant computational resources. Similarly, acquiring experimental datasets can prove difficult as well. For these systems, often computationally inexpensive, but in general inaccurate, models, known as the low-fidelity models, are available. In this paper, we propose a bi-fidelity modeling approach for complex physical systems, where we model the discrepancy between the true system's response and low-fidelity response in the presence of a small training dataset from the true system's response using a deep operator network (DeepONet), a neural network architecture suitable for approximating nonlinear operators. We apply the approach to model systems that have parametric uncertainty and are partially unknown. Three numerical examples are used to show the efficacy of the proposed approach to model uncertain and partially unknown complex physical systems.
△ Less
Submitted 18 August, 2022; v1 submitted 3 April, 2022;
originally announced April 2022.
-
Uniform-in-Phase-Space Data Selection with Iterative Normalizing Flows
Authors:
Malik Hassanaly,
Bruce A. Perry,
Michael E. Mueller,
Shashank Yellapantula
Abstract:
Improvements in computational and experimental capabilities are rapidly increasing the amount of scientific data that is routinely generated. In applications that are constrained by memory and computational intensity, excessively large datasets may hinder scientific discovery, making data reduction a critical component of data-driven methods. Datasets are growing in two directions: the number of d…
▽ More
Improvements in computational and experimental capabilities are rapidly increasing the amount of scientific data that is routinely generated. In applications that are constrained by memory and computational intensity, excessively large datasets may hinder scientific discovery, making data reduction a critical component of data-driven methods. Datasets are growing in two directions: the number of data points and their dimensionality. Whereas dimension reduction typically aims at describing each data sample on lower-dimensional space, the focus here is on reducing the number of data points. A strategy is proposed to select data points such that they uniformly span the phase-space of the data. The algorithm proposed relies on estimating the probability map of the data and using it to construct an acceptance probability. An iterative method is used to accurately estimate the probability of the rare data points when only a small subset of the dataset is used to construct the probability map. Instead of binning the phase-space to estimate the probability map, its functional form is approximated with a normalizing flow. Therefore, the method naturally extends to high-dimensional datasets. The proposed framework is demonstrated as a viable pathway to enable data-efficient machine learning when abundant data is available. An implementation of the method is available in a companion repository (https://github.com/NREL/Phase-space-sampling).
△ Less
Submitted 27 February, 2023; v1 submitted 28 December, 2021;
originally announced December 2021.
-
GANISP: a GAN-assisted Importance SPlitting Probability Estimator
Authors:
Malik Hassanaly,
Andrew Glaws,
Ryan N. King
Abstract:
Designing manufacturing processes with high yield and strong reliability relies on effective methods for rare event estimation. Genealogical importance splitting reduces the variance of rare event probability estimators by iteratively selecting and replicating realizations that are headed towards a rare event. The replication step is difficult when applied to deterministic systems where the initia…
▽ More
Designing manufacturing processes with high yield and strong reliability relies on effective methods for rare event estimation. Genealogical importance splitting reduces the variance of rare event probability estimators by iteratively selecting and replicating realizations that are headed towards a rare event. The replication step is difficult when applied to deterministic systems where the initial conditions of the offspring realizations need to be modified. Typically, a random perturbation is applied to the offspring to differentiate their trajectory from the parent realization. However, this random perturbation strategy may be effective for some systems while failing for others, preventing variance reduction in the probability estimate. This work seeks to address this limitation using a generative model such as a Generative Adversarial Network (GAN) to generate perturbations that are consistent with the attractor of the dynamical system. The proposed GAN-assisted Importance SPlitting method (GANISP) improves the variance reduction for the system targeted. An implementation of the method is available in a companion repository (https://github.com/NREL/GANISP).
△ Less
Submitted 28 December, 2021;
originally announced December 2021.
-
Adversarial sampling of unknown and high-dimensional conditional distributions
Authors:
Malik Hassanaly,
Andrew Glaws,
Karen Stengel,
Ryan N. King
Abstract:
Many engineering problems require the prediction of realization-to-realization variability or a refined description of modeled quantities. In that case, it is necessary to sample elements from unknown high-dimensional spaces with possibly millions of degrees of freedom. While there exist methods able to sample elements from probability density functions (PDF) with known shapes, several approximati…
▽ More
Many engineering problems require the prediction of realization-to-realization variability or a refined description of modeled quantities. In that case, it is necessary to sample elements from unknown high-dimensional spaces with possibly millions of degrees of freedom. While there exist methods able to sample elements from probability density functions (PDF) with known shapes, several approximations need to be made when the distribution is unknown. In this paper the sampling method, as well as the inference of the underlying distribution, are both handled with a data-driven method known as generative adversarial networks (GAN), which trains two competing neural networks to produce a network that can effectively generate samples from the training set distribution. In practice, it is often necessary to draw samples from conditional distributions. When the conditional variables are continuous, only one (if any) data point corresponding to a particular value of a conditioning variable may be available, which is not sufficient to estimate the conditional distribution. This work handles this problem using an a priori estimation of the conditional moments of a PDF. Two approaches, stochastic estimation, and an external neural network are compared here for computing these moments; however, any preferred method can be used. The algorithm is demonstrated in the case of the deconvolution of a filtered turbulent flow field. It is shown that all the versions of the proposed algorithm effectively sample the target conditional distribution with minimal impact on the quality of the samples compared to state-of-the-art methods. Additionally, the procedure can be used as a metric for the diversity of samples generated by a conditional GAN (cGAN) conditioned with continuous variables.
△ Less
Submitted 8 November, 2021;
originally announced November 2021.
-
Surface chemistry models for GaAs epitaxial growth and hydride cracking using reacting flow simulations
Authors:
Malik Hassanaly,
Hariswaran Sitaraman,
Kevin L. Schulte,
Aaron J. Ptak,
John Simon,
Kevin Udwary,
Jacob H. Leach,
Heather Splawn
Abstract:
Hydride Vapor Phase Epitaxy (HVPE) is a promising technology that can aid in the cost reduction of III-V materials and devices manufacturing, particularly high-efficiency solar cells for space and terrestrial applications. However, recent demonstrations of ultra fast growth rates ($\sim$ 500 $μ$m/h) via uncracked hydrides are not well described by present models for the growth. Therefore, it is ne…
▽ More
Hydride Vapor Phase Epitaxy (HVPE) is a promising technology that can aid in the cost reduction of III-V materials and devices manufacturing, particularly high-efficiency solar cells for space and terrestrial applications. However, recent demonstrations of ultra fast growth rates ($\sim$ 500 $μ$m/h) via uncracked hydrides are not well described by present models for the growth. Therefore, it is necessary to understand the kinetics of the growth process and its coupling with transport phenomena, so as to enable fast and uniform epitaxial growth. In this work, we derive a kinetic model using experimental data and integrate it into a computational fluid dynamics simulation of an HVPE growth reactor. We also modify an existing hydride cracking model that we validate against numerical simulations and experimental data. We show that the developed growth model and the improved cracking model are able to reproduce experimental growth measurements of \ce{GaAs} in an existing HVPE system.
△ Less
Submitted 22 September, 2021;
originally announced September 2021.
-
Data-driven Analysis of Relight variability of Jet Fuels induced by Turbulence
Authors:
Malik Hassanaly,
Yihao Tang,
Shivam Barwey,
Venkat Raman
Abstract:
For safety purposes, reliable reignition of aircraft engines in the event of flame blow-out is a critical requirement. Typically, an external ignition source in the form of a spark is used to achieve a stable flame in the combustor. However, such forced turbulent ignition may not always successfully relight the combustor, mainly because the state of the combustor cannot be precisely determined. Un…
▽ More
For safety purposes, reliable reignition of aircraft engines in the event of flame blow-out is a critical requirement. Typically, an external ignition source in the form of a spark is used to achieve a stable flame in the combustor. However, such forced turbulent ignition may not always successfully relight the combustor, mainly because the state of the combustor cannot be precisely determined. Uncertainty in the turbulent flow inside the combustor, inflow conditions, and spark discharge characteristics can lead to variability in sparking outcomes even for nominally identical operating conditions. Prior studies have shown that of all the uncertain parameters, turbulence is often dominant and can drastically alter ignition behavior. For instance, even when different fuels have similar ignition delay times, their ignition behavior in practical systems can be completely different. In practical operating conditions, it is challenging to understand why ignition fails and how much variation in outcomes can be expected. The focus of this work is to understand relight variability induced by turbulence for two different aircraft fuels, namely Jet-A and a variant named C1. A detailed, previously developed simulation approach is used to generate a large number of successful and failed ignition events. Using this data, the cause of misfire is evaluated based on a discriminant analysis that delineates the difference between turbulent initial conditions that lead to ignition or failure. It was determined that the cause of ignition failure is different for the two fuels. While it was found that Jet-A is influenced by fuel entrainment, C1 was found to be more sensitive to small scale turbulence features. A larger variability is found in the ignition modes of C1, which can be subject to extreme events induced by kernel breakdown.
△ Less
Submitted 12 November, 2020;
originally announced November 2020.
-
A self-similarity principle for the computation of rare event probability
Authors:
Malik Hassanaly,
Venkat Raman
Abstract:
The probability of rare and extreme events is an important quantity for design purposes. However, computing the probability of rare events can be expensive because only a few events, if any, can be observed. To this end, it is necessary to accelerate the observation of rare events using methods such as the importance splitting technique, which is the main focus here. In this work, it is shown how…
▽ More
The probability of rare and extreme events is an important quantity for design purposes. However, computing the probability of rare events can be expensive because only a few events, if any, can be observed. To this end, it is necessary to accelerate the observation of rare events using methods such as the importance splitting technique, which is the main focus here. In this work, it is shown how a genealogical importance splitting technique can be made more efficient if one knows how the rare event occurs in terms of the mean path followed by the observables. Using Monte Carlo simulations, it is shown that one can estimate this path using less rare paths. A self-similarity model is formulated and tested using an a priori and a posteriori analysis. The self-similarity principle is also tested on more complex systems including a turbulent combustion problem with $10^7$ degrees of freedom. While the self-similarity model is shown to not be strictly valid in general, it can still provide a good approximation of the rare mean paths and is a promising route for obtaining the statistics of rare events in chaotic high-dimensional systems.
△ Less
Submitted 31 October, 2019;
originally announced November 2019.
-
Lyapunov spectrum of forced homogeneous isotropic turbulent flows
Authors:
Malik Hassanaly,
Venkat Raman
Abstract:
In order to better understand deviations from equilibrium in turbulent flows, it is meaningful to characterize the dynamics rather than the statistics of turbulence. To this end, the Lyapunov theory provides a useful description of turbulence through the study of the perturbation dynamics. In this work, the Lyapunov spectrum of forced homogeneous isotropic turbulent flows is computed. Using the Ly…
▽ More
In order to better understand deviations from equilibrium in turbulent flows, it is meaningful to characterize the dynamics rather than the statistics of turbulence. To this end, the Lyapunov theory provides a useful description of turbulence through the study of the perturbation dynamics. In this work, the Lyapunov spectrum of forced homogeneous isotropic turbulent flows is computed. Using the Lyapunov exponents of a flow at different Reynolds numbers, the scaling of the dimension of the chaotic attractor for a three-dimensional homogeneous isotropic flow (HIT) is obtained for the first time through direct computation. The obtained Gram-Schmidt vectors (GSV) are analyzed. For the range of conditions studied, it was found that the chaotic response of the flow coincides with regions of large velocity gradients at lower Reynolds numbers and enstrophy at higher Reynolds numbers, but does not coincide with regions of large kinetic energy. Further, the response of the flow to perturbations is more and more localized as the Reynolds number increases. Finally, the energy spectrum of the GSV is computed and is shown to be almost insensitive to the Lyapunov index.
△ Less
Submitted 25 October, 2019;
originally announced October 2019.
-
Using machine learning to construct velocity fields from OH-PLIF images
Authors:
Shivam Barwey,
Malik Hassanaly,
Venkat Raman,
Adam Steinberg
Abstract:
This work utilizes data-driven methods to morph a series of time-resolved experimental OH-PLIF images into corresponding three-component planar PIV fields in the closed domain of a premixed swirl combustor. The task is carried out with a fully convolutional network, which is a type of convolutional neural network (CNN) used in many applications in machine learning, alongside an existing experiment…
▽ More
This work utilizes data-driven methods to morph a series of time-resolved experimental OH-PLIF images into corresponding three-component planar PIV fields in the closed domain of a premixed swirl combustor. The task is carried out with a fully convolutional network, which is a type of convolutional neural network (CNN) used in many applications in machine learning, alongside an existing experimental dataset which consists of simultaneous OH-PLIF and PIV measurements in both attached and detached flame regimes. Two types of models are compared: 1) a global CNN which is trained using images from the entire domain, and 2) a set of local CNNs, which are trained only on individual sections of the domain. The locally trained models show improvement in creating mappings in the detached regime over the global models. A comparison between model performance in attached and detached regimes shows that the CNNs are much more accurate across the board in creating velocity fields for attached flames. Inclusion of time history in the PLIF input resulted in small noticeable improvement on average, which could imply a greater physical role of instantaneous spatial correlations in the decoding process over temporal dependencies from the perspective of the CNN. Additionally, the performance of local models trained to produce mappings in one section of the domain is tested on other, unexplored sections of the domain. Interestingly, local CNN performance on unseen domain regions revealed the models' ability to utilize symmetry and antisymmetry in the velocity field. Ultimately, this work shows the powerful ability of the CNN to decode the three-dimensional PIV fields from input OH-PLIF images, providing a potential groundwork for a very useful tool for experimental configurations in which accessibility of forms of simultaneous measurements are limited.
△ Less
Submitted 22 September, 2019;
originally announced September 2019.
-
Experimental Data Based Reduced Order Model for Analysis and Prediction of Flame Transition in Gas Turbine Combustors
Authors:
Shivam Barwey,
Malik Hassanaly,
Qiang An,
Venkat Raman,
Adam Steinberg
Abstract:
In lean premixed combustors, flame stabilization is an important operational concern that can affect efficiency, robustness and pollutant formation. The focus of this paper is on flame lift-off and re-attachment to the nozzle of a swirl combustor. Using time-resolved experimental measurements, a data-driven approach known as cluster-based reduced order modeling (CROM) is employed to 1) isolate key…
▽ More
In lean premixed combustors, flame stabilization is an important operational concern that can affect efficiency, robustness and pollutant formation. The focus of this paper is on flame lift-off and re-attachment to the nozzle of a swirl combustor. Using time-resolved experimental measurements, a data-driven approach known as cluster-based reduced order modeling (CROM) is employed to 1) isolate key flow patterns and their sequence during the flame transitions, and 2) formulate a forecasting model to predict the flame instability. The flow patterns isolated by the CROM methodology confirm some of the experimental conclusions about the flame transition mechanism. In particular, CROM highlights the key role of the precessing vortex core (PVC) in the flame detachment process in an unsupervised manner. For the attachment process, strong flow recirculation far from the nozzle appears to drive the flame upstream, thus initiating re-attachment. Different data-types (velocity field, OH concentration) were processed by the modeling tool, and the predictive capabilities of these different models are also compared. It was found that the swirling velocity possesses the best predictive properties, which gives a supplemental argument for the role of the PVC in causing the flame transition. The model is tested against unseen data and successfully predicts the probability of flame transition (both detachment and attachment) when trained with swirling velocity with minimal user input. The model trained with OH-PLIF data was only successful at predicting the flame attachment, which implies that different physical mechanisms are present for different types of flame transition. Overall, these aspects show the great potential of data-driven methods, particularly probabilistic forecasting techniques, in analyzing and predicting large-scale features in complex turbulent combustion problems.
△ Less
Submitted 31 March, 2019;
originally announced April 2019.
-
Emerging Trends in Numerical Simulations of Combustion Systems
Authors:
Venkat Raman,
Malik Hassanaly
Abstract:
Numerical simulations have played a vital role in the design of modern combustion systems. Over the last two decades, the focus of research has been on the development of the large eddy simulation (LES) approach, which leveraged the vast increase in computing power to dramatically improve predictive accuracy. Even with the anticipated increase in supercomputing capabilities, the use of LES in desi…
▽ More
Numerical simulations have played a vital role in the design of modern combustion systems. Over the last two decades, the focus of research has been on the development of the large eddy simulation (LES) approach, which leveraged the vast increase in computing power to dramatically improve predictive accuracy. Even with the anticipated increase in supercomputing capabilities, the use of LES in design is limited by its high computational cost. Moreover, to aid decision making, such LES computations have to be augmented to estimate underlying uncertainties in simulation components. At the same time, other changes are happening across industries that build or use combustion devices. While efficiency and emissions reduction are still the primary design objectives, reducing cost of operation by optimizing maintenance and repair is becoming an important segment of the enterprise. This latter quest is aided by the digitization of combustors, which allows collection and storage of operational data from a host of sensors over a fleet of devices. Moreover, several levels of computing including low-power hardware present on board the combustion systems are becoming available. Such large data sets create unique opportunities for design and maintenance if appropriate numerical tools are made available. As LES revolutionized computing-guided design by leveraging supercomputing, a new generation of numerical approaches is needed to utilize this vast amount of data and the varied nature of computing hardware. In this article, a review of emerging computational approaches for this heterogeneous data-driven environment is provided. A case is made that new but unconventional opportunities for physics-based combustion modeling exist in this realm.
△ Less
Submitted 22 August, 2018;
originally announced August 2018.
-
Ensemble-LES Analysis of Perturbation Response of Turbulent Partially-Premixed Flames
Authors:
Malik Hassanaly,
Venkat Raman
Abstract:
The response to small perturbations of the Sandia CH4/air flame series is studied using a dynamical systems formula- tion. Here, the Sandia D and E turbulent partially-premixed flames are computed using large eddy simulation (LES) with a flamelet/progress variable approach. Using 300 simultaneous LES computations of each flame, the partial Lyapunov spectrum of the flow is obtained, which provides…
▽ More
The response to small perturbations of the Sandia CH4/air flame series is studied using a dynamical systems formula- tion. Here, the Sandia D and E turbulent partially-premixed flames are computed using large eddy simulation (LES) with a flamelet/progress variable approach. Using 300 simultaneous LES computations of each flame, the partial Lyapunov spectrum of the flow is obtained, which provides a set of Lyapunov exponents (LEs) and Lyapunov vectors (LVs). Special numerical procedures to handle such large data sets and perform matrix manipulations on-the-fly are discussed. From the LEs, the dimension of the attractor for the flames was estimated, for the first time, to be at least 5000 (flame D) and 10500 (flame E). This is a significant result indicating that the chaoticity of the turbulent flame is confined to a volume of phase-space that is quantified. The growth of perturbations occurs at time-scales that are much larger than the smallest strain-rate based times scales, but smaller than the integral scales. Further, the Lya- punov vectors that determine the perturbation response show that the strongest perturbations are aligned with the jet breakdown region, while other LVs are aligned with the shear layer. It is found that both the extinction/re-ignition process and the turbulence-induced jet motion strongly contribute to the flame chaoticity but are not equally important for all the LVs. These results provide new insights into the perturbation response of turbulent flames. More broadly, this Lyapunov approach allows a full characterization of the chaotic dynamics associated with turbulent flames, which have applications for both modeling and control of such complex flows.
△ Less
Submitted 8 June, 2018;
originally announced June 2018.
-
Numerical convergence of the Lyapunov spectrum computed using low Mach number solvers
Authors:
Malik Hassanaly,
Venkat Raman
Abstract:
In the dynamical systems approach to describing turbulent or otherwise chaotic flows, an important quantity is the Lyapunov exponents and vectors that characterize the strange attractor of the flow. In particular, knowledge of the Lyapunov exponents and vectors will help identify perturbations that the system is most sensitive to, and quantify the dimension of the attractor. However, reliably comp…
▽ More
In the dynamical systems approach to describing turbulent or otherwise chaotic flows, an important quantity is the Lyapunov exponents and vectors that characterize the strange attractor of the flow. In particular, knowledge of the Lyapunov exponents and vectors will help identify perturbations that the system is most sensitive to, and quantify the dimension of the attractor. However, reliably computing these Lyapunov quantities requires robust nu- merical algorithms. While several perturbation-based techniques are available in literature, their application to commonly used turbulent flow solvers as well the numerical convergence properties have not been studied in detail. The goal of this work is two-fold: a) develop a robust algorithm for obtaining Lyapunov exponents and vectors for low-Mach based sim- ulation of turbulent flows, b) quantify the spatial and temporal convergence properties of this algorithm using a series of progressively complex flow problems. In particular, a manu- factured solutions approach is devised using the Orr-Sommerfeld (OS) perturbation theory to extract Lyapunov exponents and vectors from OS eigen solutions. While individual test cases show interesting results, overall the spatial convergence rates of Lyapunov exponents follow the truncation order for the discretization scheme. However, temporal convergence rates were found to be only a weak function of time step used. Additionally, for some config- urations, the convergence properties are found to depend on the Lyapunov exponent itself. These results indicate that convergence properties do not follow universal rates, and require careful analysis for specific configurations considered.
△ Less
Submitted 1 May, 2018;
originally announced May 2018.
-
A minimally-dissipative low-Mach number solver for complex reacting flows in OpenFOAM
Authors:
Malik Hassanaly,
Heeseok Koo,
Christopher Lietz,
Shao Teng Chong,
Venkat Raman
Abstract:
Large eddy simulation (LES) has become the de-facto computational tool for modeling complex reacting flows, especially in gas turbine applications. However, readily usable general-purpose LES codes for complex geometries are typically academic or proprietary/commercial in nature. The objective of this work is to develop and disseminate an open source LES tool for low-Mach number turbulent combusti…
▽ More
Large eddy simulation (LES) has become the de-facto computational tool for modeling complex reacting flows, especially in gas turbine applications. However, readily usable general-purpose LES codes for complex geometries are typically academic or proprietary/commercial in nature. The objective of this work is to develop and disseminate an open source LES tool for low-Mach number turbulent combustion using the OpenFOAM framework. In particular, a collocated-mesh approach suited for unstructured grid formulation is provided. Unlike other fluid dynamics models, LES accuracy is intricately linked to so-called primary and secondary conservation properties of the numerical discretization schemes. This implies that although the solver only evolves equations for mass, momentum, and energy, the implied discrete equation for kinetic energy (square of velocity) should be minimally-dissipative. Here, a specific spatial and temporal discretization is imposed such that this kinetic energy dissipation is minimized. The method is demonstrated using manufactured solutions approach on regular and skewed meshes, a canonical flow problem, and a turbulent sooting flame in a complex domain relevant to gas turbines applications.
△ Less
Submitted 12 May, 2017;
originally announced May 2017.
-
Classification and Simulation of Anomalous Events in Turbulent Combustion
Authors:
Malik Hassanaly,
Stephen Voelkel,
Venkat Raman
Abstract:
In practical combustion systems, the high nonlinearity of the turbulent combustion process can result in unexpected behavior even for nominal operating conditions. This includes flame flashback in premixed gas turbines, engine unstart in scramjets, and failure of engines to ignite at high altitudes. Many of these events have catastrophic consequences. Hence, understanding the nature of such irregu…
▽ More
In practical combustion systems, the high nonlinearity of the turbulent combustion process can result in unexpected behavior even for nominal operating conditions. This includes flame flashback in premixed gas turbines, engine unstart in scramjets, and failure of engines to ignite at high altitudes. Many of these events have catastrophic consequences. Hence, understanding the nature of such irregular excursions of combustion systems is important to ensure robust operation. The focus of this work is to identify the sources of such rare excursions from a dynamical systems perspective. Specifically, such anomalous behavior can be explained in terms of three different classes based on the trajectories that the system traverses in its phase-space. For each class, computational algorithms and validation exercises are briefly discussed.
△ Less
Submitted 27 April, 2017;
originally announced May 2017.