GPU-accelerated discontinuous Galerkin methods on hybrid meshes

We present a time-explicit discontinuous Galerkin (DG) solver for the time-domain acoustic wave equation on hybrid meshes containing vertex-mapped hexahedral, wedge, pyramidal and tetrahedral elements. Discretely energy-stable formulations are presented for both Gauss-Legendre and Gauss-Legendre-Lobatto (Spectral Element) nodal bases for the hexahedron. Stable timestep restrictions for hybrid meshes are derived by bounding the spectral radius of the DG operator using order-dependent constants in trace and Markov inequalities. Computational efficiency is achieved under a combination of element-specific kernels (including new quadrature-free operators for the pyramid), multi-rate timestepping, and acceleration using Graphics Processing Units.Comment: Submitted to CMAM

Analysis of time-harmonic electromagnetic problems with elliptic material coefficients

We consider time-harmonic electromagnetic problems with material coefficients represented by elliptic fields, covering a wide range of complex and anisotropic material media. The properties of elliptic fields are analyzed, with particular emphasis on scalar fields and normal tensor fields. Time-harmonic electromagnetic problems with general elliptic material fields are then studied. Well-posedness results for classical variational formulations with different boundary conditions are reviewed, and hypotheses for the coercivity of the corresponding sesquilinear forms are investigated. Finally, the proposed framework is applied to examples of media used in the literature: isotropic lossy media, geometric media, and gyrotropic media

An efficient domain decomposition method with cross-point treatment for Helmholtz problems

National audienceThe parallel finite-element solution of large-scale time-harmonic scattering problems is addressed with a non-overlapping domain decomposition method. The efficiency of this method strongly depends on the transmission condition enforced on the interfaces between the subdomains. Conditions based on high-order local absorbing boundary conditions have proved well suited for configurations without cross points (i.e. points where more than two subdomains meet). In this work, we extend this approach to efficiently deal with cross points. Two-dimensional numerical results are presented

Corner treatments for high-order local absorbing boundary conditions in high-frequency acoustic scattering

International audienceThis paper deals with the design and validation of accurate local absorbing boundary conditions set on convex polygonal and polyhedral computational domains for the finite element solution of high-frequency acoustic scattering problems. While high-order absorbing boundary conditions (HABCs) are accurate for smooth fictitious boundaries, the precision of the solution drops in the presence of corners if no specific treatment is applied. We present and analyze two strategies to preserve the accuracy of Padé-type HABCs at corners: first by using compatibility relations (derived for right angle corners) and second by regularizing the boundary at the corner. Exhaustive numerical results for two- and three-dimensional problems are reported in the paper. They show that using the compatibility relations is optimal for domains with right angles. For the other cases, the error still remains acceptable, but depends on the choice of the corner treatment according to the angle

An efficient domain decomposition method with cross-point treatment for Helmholtz problems

International audienceThe parallel finite-element solution of large-scale time-harmonic scattering problems is addressed with a non-overlapping domain decomposition method (DDM). It is well known that the efficiency of this method strongly depends on the transmission condition enforced on the interfaces between the subdomains. Local conditions based on high-order absorbing boundary conditions (HABCs) are well suited for configurations without cross points (where more than two subdo-mains meet). In this work, we extend this approach to efficiently deal with cross points. Two-dimensional finite-element results are presented

An automatic PML for acoustic finite element simulations in convex domains of general shape

International audienceThis article addresses the efficient finite element solution of exterior acoustic problems with truncated computational domains surrounded by perfectly matched layers (PMLs). The PML is a popular nonreflecting technique that combines accuracy, computational efficiency, and geometric flexibility. Unfortunately, the effective implementation of the PML for convex domains of general shape is tricky because of the geometric parameters that are required to define the PML medium. In this work, a comprehensive implementation strategy is proposed. This approach, which we call the automatically matched layer (AML) implementation, is versatile and fully automatic for the end‐user. With the AML approach, the mesh of the layer is extruded, the required geometric parameters are automatically obtained during the extrusion step, and the practical implementation relies on a simple modification of the Jacobian matrix in the elementwise integrals. The AML implementation is validated and compared with other implementation strategies using numerical benchmarks in two and three dimensions, considering computational domains with regular and nonregular boundaries. A three‐dimensional application with a generally shaped domain generated using a convex hull is proposed to illustrate the interest of the AML approach for realistic industrial cases

An efficient domain decomposition method with cross-point treatment for Helmholtz problems

National audienceThe parallel finite-element solution of large-scale time-harmonic scattering problems is addressed with a non-overlapping domain decomposition method. The efficiency of this method strongly depends on the transmission condition enforced on the interfaces between the subdomains. Conditions based on high-order local absorbing boundary conditions have proved well suited for configurations without cross points (i.e. points where more than two subdomains meet). In this work, we extend this approach to efficiently deal with cross points. Two-dimensional numerical results are presented