Research

Scientific Machine Learning (SciML)

Generative models

Generative models can be useful for manipulating high-dimensional probability density functions which are often critical for closure modeling, state estimation or data curation. I develop methods for generative adversarial networks (GANs) for uncertainty-aware state estimation of atmospheric flows.

Since samples generated with GANs are typically of high quality, they can be amenable to forward simulations that can assist uncertainty quantification tasks that use ensembles of forward simulations (typical in atmospheric modeling). When estimating rare-event probability, the sampling strategy can be critical and can benefit from a well-informed generative procedure.

Internally, some generative models directly learn the PDF that is being manipulated. Crucially, they work well in high dimensions. The PDF estimate can be used for various scientific tasks. For instance, I use them to curate datasets and transform inhomogeneously distributed data into a uniform-in-phase-space dataset.

  1. M. Hassanaly, A. Glaws, K. Stengel, R. N. King, “Adversarial sampling of unknown and high-dimensional conditional distributions” in Journal of Computational Physics, 2022. [PDF], [Code]
  2. M. Hassanaly, A. Glaws, R. N. King, “GANISP: a GAN-assisted Importance SPlitting Probability Estimator” in AAAI-ADAM, 2022. [PDF], [Code]
  3. M. Hassanaly, B. Perry, M. E. Mueller, S. Yellapantula, “Uniform-in-Phase-Space Data Selection with Iterative Normalizing Flows”, Data-Centric Engineering, 2023. [PDF], [Code]
  4. A. Rybchuk, M. Hassanaly, N. Hamilton, P. Doubrawa, M. J. Fulton, L. A. Martinez-Tossas, “Ensemble flow reconstruction in the atmospheric boundary layer from spatially limited measurements through latent diffusion models”, Physics of Fluids, 2023. [PDF], [Code]
  5. A. Rybchuk, L. A Martínez-Tossas, N. Hamilton, P. Doubrawa, G. Vijayakumar, M. Hassanaly, M. B. Kuhn and D. S. Zalkind, “A baseline for ensemble-based, time-resolved inflow reconstruction for a single turbine using large-eddy simulations and latent diffusion models”, Wake Conference, 2023. [PDF]
  6. H.-W. Pang, M. Hassanaly, B. Perry, M. Day, W. H. Green, “Iterative Workflow for Quantification and Minimization of Reduced Chemistry-Induced Uncertainties in Reacting Flow Simulation”, WIPP 39th International Symposium on Combustion, 2022. [Poster]
  7. S. Barwey, M. Hassanaly, V. Raman, A. Steinberg, “Using Machine Learning to Construct Velocity Fields from OH-PLIF Images”, in Combustion Science and Technology, 2019. [PDF]

Information extraction

Despite large numerical and experimental available datasets it is not always straightforward to extract useful information for design purpose. I have found that reducing datasets either via sparse-sensing or data clustering can be a valuable tool. Sparse sensing can allow to pinpoint the cause of a variability in observations. Practically, we have used it to understand the cause being ignition kernels that catch fire or that that blow-out in aircraft engines.

Clustering can be useful to extract an interpretable and statistically significant sequence of states that explain a given phenomenon. I have used clustering to identify ignition mode and blow-out modes for ignition kernels that interact with different fuels. I used it to understand why a flame may attack or detach from a nozzle in a swirl combustor.

  1. M. Hassanaly, Y. Tang, S. Barwey, V. Raman, “Data-driven analysis of relight variability of jet fuels induced by turbulence” in Combustion and Flame, 2021. [PDF]
  2. S. Barwey, M. Hassanaly, Q. An, V. Raman, A. Steinberg, “Experimental data-based reduced-order model for analysis and prediction of flame transition in gas turbine combustors” in Combustion Theory and Modelling, 2019. [PDF]
  3. S. Barwey, H. Ganesh, M. Hassanaly, V. Raman, S. Ceccio, “Data-based analysis of multimodal partial cavity shedding dynamics” in Experiments in Fluids, 2020. [PDF]

Complex fluid flows and High-Performance Computing (HPC)

Reacting flows

Reacting flows modeling is a challenging computational problem that require efficiently modeling microscopic scale phenomena that influence macroscopic scales, without resolving the entire range of scales. I work on population-balance equations for modeling of soot particles. Additionally, reacting flows require resolving stiff systems of ODE that can accurately model complex chemistry. However, chemical composition usually spans a reduced part of composition space. It is therefore possible to represent the chemistry on low-dimensional manifold that can represent the dynamics of the chemistry at a reduced computational cost. I help develop this approach for modeling turbulent ignition of jet fuels. Alternatively, surrogate reduced chemistry models can be developed to replace detailed chemisty approaches. I develop such a strategy for surface chemistry models in deposition reactors.

  1. M. Hassanaly, H. Sitaraman, K. L. Schulte, A. J. Ptak, J. Simon, K. Udwary, J. H. Leach, H. Splawn, “Surface chemistry models for GaAs epitaxial growth and hydride cracking using reacting flow simulations” in Journal of Applied Physics, 2021. [PDF]
  2. M. Hassanaly, H. Sitaraman, K. Udwary, K. L. Schulte, H. Splawn, A. Ptak, J. Simon, “Optimization of Hydride Vapor Phase Epitaxy (HVPE) Deposition Reactor Manufacturing III-V Materials Using Multiple Large-Eddy Simulations (LES)” in AIChE Annual Meeting, 2020. [Abstract]
  3. Y. Tang, M. Hassanaly, V. Raman, B. Sforzo, J. Seitzman, “A comprehensive modeling procedure for estimating statistical properties of forced ignition” in Combustion and Flame, 2019. [PDF]
  4. H. Koo, M. Hassanaly, V. Raman, M. E. Mueller, K.-P. Geigle, “Large-eddy simulation of soot formation in a model gas turbine combustor” in Journal of Engineering for Gas Turbines and Power, 2017. [PDF]
  5. S. T. Chong, M. Hassanaly, H. Koo, M. E. Mueller, V. Raman, K.-P. Geigle, “Large eddy simulation of pressure and dilution-jet effects on soot formation in a model aircraft swirl combustor” in Combustion and Flame, 2018. [PDF]
  6. N. Wimer, L. Esclapez, N. Brunhart-Lupo, M. Henry de Frahan, M. Rahimi, M. Hassanaly, J. Rood, S. Yellapantula, H. Sitaraman, B. Perry, M. Martin, O. Doronina, S. Appukuttan, M. Rieth, M. Day, “Visualizations of a methane/diesel RCCI engine using PeleC and PeleLMeX” in 75th Annual Meeting of the APS Division of Fluid Dynamics, 2022. [Gallery of Fluid Motion Video] 🏆 Milton Van Dyke Award [PDF]
  7. J. Parra-Alvarez, M. Hassanaly, M. J. Rahimi, H. Sitaraman, “Computational Analysis of Different Sparging Systems and Their Influence in the Fluid-Dynamic Behavior of Bubble Column Reactors”, 2023. [Code]

Numerical solver

No matter how accurate are physics models, they intimately interact with numerical solvers which can introduce larger errors than the physics-modeling approach. For instance, numerical discretization may introduce dissipation that outweighs the effect of turbulence models. I have developed a solver that minimize the amount of numerical dissipation of variable density low-Mach solvers in OpenFoam.

In turn, numerical solvers also depend on the computing architecture on which they are deployed. With the rise of GPU/CPU computing architectures, writing efficient code requires appropriate memory management. I have developed an approach where numerical chemistry Jacobian can be symbolically encoded, allowing to control the common subexpressions precomputations. A tradeoff between memory and computational intensity can be done that way. The approach has been deployed on Exascale machines.

  1. M. Hassanaly, H. Koo, C. F. Lietz, S. T. Chong, V. Raman, “A minimally-dissipative low-Mach number solver for complex reacting flows in OpenFOAM” in Computer and Fluids, 2018. [PDF]
  2. [Documentation] on common-subexpression precomputations, and symbolic chemical Jacobian.

Uncertainty quantification (UQ)

Bayesian inference

Inverse Bayesian parametric inference can allow objective system identification given a limited dataset. In the context of Li-ion battery degradation, data availability is usually sparse, and the degradation may occur through different parameter sets. Inverse parametric inference can help identify the degradation mode or decide that more data is needed. However, the cost of parametric inference may be prohibitive when complex forward physics models are used. I develop surrogate approaches based on physics-informed neural networks to replace expensive Li-ion battery models.

  1. M. Hassanaly, P. J. Weddle, R. N. King, S. De, A. Doostan, C. R. Randall, E. J. Dufek, A. M. Colclasure, K. Smith, “PINN surrogate of Li-ion battery models for parameter inference. Part I: Implementation and multi-fidelity hierarchies for the single-particle model”, in Under Review, 2024. [Paper] [Code]
  2. M. Hassanaly, P. J. Weddle, R. N. King, S. De, A. Doostan, C. R. Randall, E. J. Dufek, A. M. Colclasure, K. Smith, “PINN surrogate of Li-ion battery models for parameter inference. Part II: Regularization and application of the pseudo-2D model”, in Under Review, 2024. [Paper]
  3. M. Hassanaly, P. J. Weddle, K. Smith, S. De, A. Doostan, R. N. King, “Physics-Informed Neural Network Modeling of Li-Ion Batteries”, in 242nd ECS Meeting, 2022. [Abstract]

Rare event probability estimation

Efficiently estimating rare event probability requires artificially increasing the number of rare events observed. I develop methods for encouraging the occurrence of rare events in physics simulations.

  1. M. Hassanaly, V. Raman, “A self-similarity principle for the computation of rare event probability” in Journal of Physics A: Mathematical and Theoretical, 2019. [PDF]
  2. M. Hassanaly, A. Glaws, R. N. King, “GANISP: a GAN-assisted Importance SPlitting Probability Estimator” in AAAI-ADAM, 2022. [PDF], [Code]

Uncertainty propagation

Coming Soon!