Stabilizing PDE--ML coupled systems

Partial differential equations (PDEs) are an essential modeling tool in engineering and physical sciences. The numerical methods used for solving the more descriptive and sophisticated of these models comprise many computationally expensive modules. Machine learning (ML) provides a way of replacing some of these modules by surrogates that are much more efficient at the time of inference. The resulting PDE-ML coupled systems, however, can be highly susceptible to instabilities [1-3]. Efforts towards ameliorating these have mostly concentrated on improving the accuracy of the surrogates, imbuing them with additional structure, or introducing problem-specific stabilizers, and have garnered limited success [4-7]. In this article, we study a prototype problem to understand the mathematical subtleties involved in PDE-ML coupling, and draw insights that can help with more complex systems.