|
Escapil-Inchauspe, P., & Jerez-Hanckes, C. (2020). Helmholtz Scattering by Random Domains: First-Order Sparse Boundary Elements Approximation. SIAM J. Sci. Comput., 42(5), A2561–A2592.
Abstract: We consider the numerical solution of time-harmonic acoustic scattering by obstacles with uncertain geometries for Dirichlet, Neumann, impedance, and transmission boundary conditions. In particular, we aim to quantify diffracted fields originated by small stochastic perturbations of a given relatively smooth nominal shape. Using first-order shape Taylor expansions, we derive tensor deterministic first-kind boundary integral equations for the statistical moments of the scattering problems considered. These are then approximated by sparse tensor Galerkin discretizations via the combination technique [M. Griebel, M. Schneider, and C. Zenger, A combination technique for the solution of sparse grid problems, in Iterative Methods in Linear Algebra, P. de Groen and P. Beauwens, eds., Elsevier, Amsterdam, 1992, pp. 263-281; H. Harbrecht, M. Peters, and M. Siebenmorgen, J. Comput. Phys., 252 (2013), pp. 128-141]. We supply extensive numerical experiments confirming the predicted error convergence rates with polylogarithmic growth in the number of degrees of freedom and accuracy in approximation of the moments. Moreover, we discuss implementation details such as preconditioning to finally point out further research avenues.
|
|
|
Escapil-Inchauspé, P., & Ruz, G. A. (2023). Hyper-parameter tuning of physics-informed neural networks: Application to Helmholtz problems. Neurocomputing, 561, 126826.
Abstract: We consider physics-informed neural networks (PINNs) (Raissiet al., 2019) for forward physical problems. In order to find optimal PINNs configuration, we introduce a hyper-parameter optimization (HPO) procedure via Gaussian processes-based Bayesian optimization. We apply the HPO to Helmholtz equation for bounded domains and conduct a thorough study, focusing on: (i) performance, (ii) the collocation points density r and (iii) the frequency kappa, confirming the applicability and necessity of the method. Numerical experiments are performed in two and three dimensions, including comparison to finite element methods.
|
|
|
Fierro, I., & Jerez-Hanckes, C. (2020). Fast Calderon preconditioning for Helmholtz boundary integral equations. J. Comput. Phys., 409, 22 pp.
Abstract: Calderon multiplicative preconditioners are an effective way to improve the condition number of first kind boundary integral equations yielding provable mesh independent bounds. However, when discretizing by local low-order basis functions as in standard Galerkin boundary element methods, their computational performance worsens as meshes are refined. This stems from the barycentric mesh refinement used to construct dual basis functions that guarantee the discrete stability of L-2-pairings. Based on coarser quadrature rules over dual cells and H-matrix compression, we propose a family of fast preconditioners that significantly reduce assembly and computation times when compared to standard versions of Calderon preconditioning for the three-dimensional Helmholtz weakly and hyper-singular boundary integral operators. Several numerical experiments validate our claims and point towards further enhancements. (C) 2020 Elsevier Inc. All rights reserved.
|
|