Numerical investigations on model order reduction to SEM based on POD-DEIM to linear/nonlinear heat transfer problems

Abstract

The motivation in this article is to foster and conduct an in-depth investigation on reduced-order modeling (ROM) techniques via proper orthogonal decomposition (POD) and discrete empirical interpolation method (DEIM) in conjunction with high-order Legendre spectral element method (SEM) applied for time-dependent linear and nonlinear heat conduction problems. This approach significantly relieves the CPU burden, yet preserves the numerical accuracy of high-order schemes; this is especially pronounced for linear transient problems. For nonlinear cases, we use DEIM to save CPU time for nonlinear counterparts. For this purpose, we first obtain training data by solving the system of full-order models (FOM) using snapshot techniques to construct the POD basis. The well-known generalized single step single solve (GSSSS-1) framework is next employed for the first-order system for the time discretization. Numerical results for both linear and nonlinear heat conduction problems such as in a nuclear reactor, etc., show excellent agreement between FOM solutions and ROM solutions, and the CPU advantages via POD ROM techniques are also prominent for high-order Legendre SEM.