  The high computational cost associated with solving for detailed chemistry
poses a significant challenge for predictive computational fluid dynamics (CFD)
simulations of turbulent reacting flows. These models often require solving a
system of coupled stiff ordinary differential equations (ODEs). While deep
learning techniques have been experimented with to develop faster surrogate
models, they often fail to integrate reliably with CFD solvers. This
instability arises because deep learning methods optimize for training error
without ensuring compatibility with ODE solvers, leading to accumulation of
errors over time. Recently, NeuralODE-based techniques have offered a promising
solution by effectively modeling chemical kinetics. In this study, we extend
the NeuralODE framework for stiff chemical kinetics by incorporating mass
conservation constraints directly into the loss function during training. This
ensures that the total mass and the elemental mass are conserved, a critical
requirement for reliable downstream integration with CFD solvers.
Proof-of-concept studies are performed with physics-constrained neuralODE
(PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a
range of composition and thermodynamic conditions. Our results demonstrate that
this enhancement not only improves the physical consistency with respect to
mass conservation criteria but also ensures better robustness. Lastly, a
posteriori studies are performed wherein the trained PC-NODE model is coupled
with a 3D CFD solver for computing the chemical source terms. PC-NODE is shown
to be more accurate relative to the purely data-driven neuralODE approach.
Moreover, PC-NODE also exhibits robustness and generalizability to unseen
initial conditions from within (interpolative capability) as well as outside
(extrapolative capability) the training regime.
