A non-intrusive model order reduction approach for parameterized time-domain Maxwell's equations

International audienceWe present a non-intrusive model order reduction (NIMOR) method with an offline-online decoupling for the solution of parameterized time-domain Maxwell's equations. During the offline stage, the training parameters are chosen by using Smolyak sparse grid method with an approximation level L (L ≥ 1) over a target parameterized space. This method can deal with the so-called curse of dimensionality in high dimensional space. For each selected parameter, the snapshot vectors are first produced by a high order discontinuous Galerkin time-domain (DGTD) solver formulated on an unstructured simplicial mesh. In order to minimize the overall computational cost in the offline stage and to improve the accuracy of the NIMOR method, a radial basis function (RBF) interpolation method is then used to construct more snapshot vectors at the sparse grid with approximation level L + 1, which includes the sparse grids from approximation level L. A nested proper orthogonal decomposition (POD) method is employed to extract time-and parameter-independent POD basis functions. By using the singular value decomposition (SVD) method, the principal components of the reduced coefficient matrices of the high-fidelity solutions onto the reduced-order subspace spaned by the POD basis functions are extracted. Moreover, a Gaussian process regression (GPR) method is proposed to approximate the dominating time-and parameter-modes of the reduced coefficient matrices. During the online stage, the reduced-order solutions for new time and parameter values can be rapidly recovered via outputs from the regression models without using the DGTD method. The performance of this NIMOR method is illustrated numerically by considering two classical test cases: the scattering of a plane wave by a 2-D dielectric disk and the scattering of a plane wave by a multi-layer heterogeneous medium. The prediction capabilities of the NIMOR method are evaluated by varying the relative permittivity. Numerical results indicate that the NIMOR method is a promising approach for simulating accurately and in fast way parameterized timedomain electromagnetic problems

A deep learning enabler for non-intrusive reduced order modeling of fluid flows

In this paper, we introduce a modular deep neural network (DNN) framework for data-driven reduced order modeling of dynamical systems relevant to fluid flows. We propose various deep neural network architectures which numerically predict evolution of dynamical systems by learning from either using discrete state or slope information of the system. Our approach has been demonstrated using both residual formula and backward difference scheme formulas. However, it can be easily generalized into many different numerical schemes as well. We give a demonstration of our framework for three examples: (i) Kraichnan-Orszag system, an illustrative coupled nonlinear ordinary differential equations, (ii) Lorenz system exhibiting chaotic behavior, and (iii) a non-intrusive model order reduction framework for the two-dimensional Boussinesq equations with a differentially heated cavity flow setup at various Rayleigh numbers. Using only snapshots of state variables at discrete time instances, our data-driven approach can be considered truly non-intrusive, since any prior information about the underlying governing equations is not required for generating the reduced order model. Our \textit{a posteriori} analysis shows that the proposed data-driven approach is remarkably accurate, and can be used as a robust predictive tool for non-intrusive model order reduction of complex fluid flows.Comment: 36 pages, 21 figure

Non-intrusive reduced-order modeling of parameterized electromagnetic scattering problems using cubic spline interpolation

International audienceThis paper presents a non-intrusive model order reduction (MOR) for the solution of parameterized electromagnetic scattering problems, which needs to prepare a database offline of full-order solution samples (snapshots) at some different parameter locations. The snapshot vectors are produced by a high order discontinuous Galerkin time-domain (DGTD) solver formulated on an unstructured simplicial mesh. Because the second dimension of snapshots matrix is large, a two-step or nested proper orthogonal decomposition (POD) method is employed to extract time- and parameter-independent POD basis functions. By using the singular value decomposition (SVD) method, the principal components of the projection coefficient matrices (also referred to as the reduced coefficient matrices) of full-order solutions onto the RB subspace are extracted. A cubic spline interpolation-based (CSI) approach is proposed to approximate the dominating time- and parameter-modes of the reduced coefficient matrices without resorting to Galerkin projection. The generation of snapshot vectors, the construction of POD basis functions and the approximation of reduced coefficient matrices based on the CSI method are completed during the offline stage. The RB solutions for new time and parameter values can be rapidly recovered via outputs from the interpolation models in the online stage. In particular, the offline and online stages of the proposed RB method, termed as the POD-CSI method, are completely decoupled, which ensures the computational validity of the method. Moreover, a surrogate error model is constructed as an efficient error estimator for the POD-CSI method. Numerical experiments for the scattering of plane wave by a 2-D dielectric cylinder and a multi-layer heterogeneous medium nicely illustrate the performance of POD-CSI method

Shape optimization through proper orthogonal decomposition with interpolation and dynamic mode decomposition enhanced by active subspaces

We propose a numerical pipeline for shape optimization in naval engineering involving two different non-intrusive reduced order method (ROM) techniques. Such methods are proper orthogonal decomposition with interpolation (PODI) and dynamic mode decomposition (DMD). The ROM proposed will be enhanced by active subspaces (AS) as a pre-processing tool that reduce the parameter space dimension and suggest better sampling of the input space. We will focus on geometrical parameters describing the perturbation of a reference bulbous bow through the free form deformation (FFD) technique. The ROM are based on a finite volume method (FV) to simulate the multi-phase incompressible flow around the deformed hulls. In previous works we studied the reduction of the parameter space in naval engineering through AS [38, 10] focusing on different parts of the hull. PODI and DMD have been employed for the study of fast and reliable shape optimization cycles on a bulbous bow in [9]. The novelty of this work is the simultaneous reduction of both the input parameter space and the output fields of interest. In particular AS will be trained computing the total drag resistance of a hull advancing in calm water and its gradients with respect to the input parameters. DMD will improve the performance of each simulation of the campaign using only few snapshots of the solution fields in order to predict the regime state of the system. Finally PODI will interpolate the coefficients of the POD decomposition of the output fields for a fast approximation of all the fields at new untried parameters given by the optimization algorithm. This will result in a non-intrusive data-driven numerical optimization pipeline completely independent with respect to the full order solver used and it can be easily incorporated into existing numerical pipelines, from the reference CAD to the optimal shape