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 -