Stochastic Operator Learning for Chemistry in Non-Equilibrium Flows

TL;DR

This paper introduces a Bayesian operator learning framework for non-equilibrium chemical kinetics that ensures physical consistency, accounts for uncertainties, and provides interpretable surrogates for hypersonic flow simulations.

Contribution

It develops a novel operator learning methodology that separates time dynamics, incorporates model uncertainty, and enhances interpretability for chemical kinetics in non-equilibrium flows.

Findings

Surrogate model remains stable during time integration.

Provides physically consistent probabilistic predictions.

Achieves maximum relative error below 10%.

Abstract

This work presents a novel framework for physically consistent model error characterization and operator learning for reduced-order models of non-equilibrium chemical kinetics. By leveraging the Bayesian framework, we identify and infer sources of model and parametric uncertainty within the Coarse-Graining Methodology across a range of initial conditions. The model error is embedded into the chemical kinetics model to ensure that its propagation to quantities of interest remains physically consistent. For operator learning, we develop a methodology that separates time dynamics from other input parameters. Karhunen-Loeve Expansion (KLE) is employed to capture time dynamics, yielding temporal modes, while Polynomial Chaos Expansion (PCE) is subsequently used to map model error and input parameters to KLE coefficients. The proposed model offers three significant advantages: i) Separating time dynamics from other inputs ensures stability of chemistry surrogate when coupled with fluid solvers; ii) The framework fully accounts for model and parametric uncertainty, enabling robust probabilistic predictions; iii) The surrogate model is highly interpretable, with visualizable time modes and a PCE component that facilitates analytical calculation of sensitivity indices. We apply this framework to O2-O chemistry system under hypersonic flight conditions, validating it in both a 0D adiabatic reactor and coupled simulations with a fluid solver in a 1D shock case. Results demonstrate that the surrogate is stable during time integration, delivers physically consistent probabilistic predictions accounting for model and parametric uncertainty, and achieves maximum relative error below 10%. This work represents a significant step forward in enabling probabilistic predictions of non-equilibrium chemistry with coupled fluid solvers, offering a physically accurate approach for hypersonic flow predictions.