Proper generalized decomposition for parameterized Helmholtz problems in heterogeneous and unbounded domains: Application to harbor agitation

Abstract

Solving the Helmholtz equation for a large number of input data in an heterogeneous media and unbounded domain still represents a challenge. This is due to the particular nature of the Helmholtz operator and the sensibility of the solution to small variations of the data. Here a reduced order model is used to determine the scattered solution everywhere in the domain for any incoming wave direction and frequency. Moreover, this is applied to a real engineering problem: water agitation inside real harbors for low to mid-high frequencies.The Proper Generalized Decomposition (PGD) model reduction approach is used to obtain a separable representation of the solution at any point and for any incoming wave direction and frequency. Here, its applicability to such a problem is discussed and demonstrated. More precisely, the contributions of the paper include the PGD implementation into a Perfectly Matched Layer framework to model the unbounded domain, and the separability of the operator which is addressed here using an e cient higher-order projection scheme.Then, the performance of the PGD in this framework is discussed and improved using the higher-order projection and a Petrov-Galerkin approach to construct the separated basis. Moreover, the e ciency of the higherorder projection scheme is demonstrated and compared with the higher-order singular value decomposition.

An efficient quasi-optimal space-time PGD application to frictional contact mechanics

“…This extension could be useful for PGD multiparametric studies and may converge to recent works concerning the High Order SVD (HOSVD) or similar tensor decomposition [28][29][30][31][32].…”

Application of PGD on Parametric Modeling of a Piezoelectric Energy Harvester

“…electrical) test and trial functions can be constructed using the Galerkin method, as in Eq. (9) . In the implementation, it is more convenient to arrange u jh ′s and φ h ′s in vectors.…”

“…For instance, in [8] operator separating in both linear and nonlinear cases are discussed. In [9], the same issue in high-dimensional cases are investigated. After separating, PGD formulations can be constructed using the Galerkin orthogonality criteria or residual minimization, depending on whether the operator is self-adjoint [10].…”

An error estimator for real-time simulators based on model order reduction

“…An additional compression of the modes with the so-called PGD-projection, see [38], provides with a very restricted number of modes. In this case, the modes could be compressed so as to consider less than 12 modes for the whole loading process without further increase in the error in the approximation.…”