TY - JOUR

T1 - An adaptive time-integration scheme for stiff chemistry based on Computational Singular Perturbation and Artificial Neural Networks

AU - Malpica Galassi, Riccardo

AU - Ciottoli, Pietro P.

AU - Valorani, Mauro

AU - Im, Hong G.

N1 - KAUST Repository Item: Exported on 2021-12-14
Acknowledgements: We acknowledge the fruitful discussions and the technical support kindly offered by Dr. Shivam Barwey and Professor Venkat Raman at the University of Michigan, Ann Arbor, MI, USA, and Dr. Mattia Soldan at KAUST, Saudi Arabia. R. Malpica Galassi acknowledges the financial support of the Fédération Wallonie-Bruxelles (FWB), Cellule Europe.

PY - 2021/11/29

Y1 - 2021/11/29

N2 - We leverage the computational singular perturbation (CSP) theory to develop an adaptive time-integration scheme for stiff chemistry based on a local, projection-based, reduced order model (ROM) freed of the fast time-scales. Its construction is such that artificial neural networks (ANN) can be plugged-in as cheap surrogates of the local projection basis, which is a state function, to alleviate the computational cost, without sacrificing the geometrical and physical foundation of the method. In fact, the solver relies on the synthetic basis in place of the more expensive on-the-fly calculated basis, i.e. the eigenvectors of the Jacobian matrix of the chemical source term, to define the local slow invariant manifold (SIM) and the projection matrix, then integrates explicitly the projected, i.e., non-stiff, chemical source term.

AB - We leverage the computational singular perturbation (CSP) theory to develop an adaptive time-integration scheme for stiff chemistry based on a local, projection-based, reduced order model (ROM) freed of the fast time-scales. Its construction is such that artificial neural networks (ANN) can be plugged-in as cheap surrogates of the local projection basis, which is a state function, to alleviate the computational cost, without sacrificing the geometrical and physical foundation of the method. In fact, the solver relies on the synthetic basis in place of the more expensive on-the-fly calculated basis, i.e. the eigenvectors of the Jacobian matrix of the chemical source term, to define the local slow invariant manifold (SIM) and the projection matrix, then integrates explicitly the projected, i.e., non-stiff, chemical source term.

UR - http://hdl.handle.net/10754/673852

UR - https://linkinghub.elsevier.com/retrieve/pii/S0021999121007701

U2 - 10.1016/j.jcp.2021.110875

DO - 10.1016/j.jcp.2021.110875

M3 - Article

SN - 0021-9991

SP - 110875

JO - Journal of Computational Physics

JF - Journal of Computational Physics

ER -