A discontinuous finite element method for the hydroelastic analysis of submerged structures
Abstract
This paper develops a finite element method for wave-structure interaction problems arising in the hydroelastic modelling of submerged elastic plates. The approach is formulated for both two-dimensional settings and three-dimensional channels. For the fluid part of the problem, we use a symmetric discontinuous Galerkin scheme that captures the discontinuity in the potential across the plate. To impose the Sommerfeld radiation conditions, we take inspiration from the Dirichlet-to-Neumann map approach and use the analytic solution in the semi-infinite domains to construct appropriate non-local boundary conditions. For the plate, we use a continuous/discontinuous Galerkin method to resolve the 4th order operator without using continuously differentiable finite elements. We prove well-posedness of the method, showing that the resulting sesquilinear form is bounded and satisfies the Gårding inequality, leading to a unique and stable solution. We implement our formulation as an open-source tool using the Gridap package in Julia. We validate our method against a hypersingular boundary integral method, finding excellent agreement. The flexibility of the methodology is demonstrated through simulations of a two-dimensional slanted plate and three-dimensional plate geometries including centred rectangular plates, offset rectangular plates, and annulus-shaped plates. The finite element method presented here can be readily applied to model wave energy conversion of piezoelectric bimorphs. Furthermore, the method could easily be extended to consider wave scattering from fixed or rigid structures or from variable bottom topography, and is not limited to hydroelastic applications.
Related
- Anisotropic error estimates for the L2-projection onto discontinuous finite element spaces
- A Hybridized Staggered Discontinuous Galerkin--Mixed Finite Element Method for Strain Gradient Elasticity
- A discontinuous Petrov-Galerkin finite-element framework for the simulation of microwave-heated flows
- High-Order Discontinuous Cut Finite Element Methods for Scalar Hyperbolic Conservation Laws
- Least-Squares Finite Element Methods with Weakly Imposed Boundary Conditions: Norm and Error Equivalences, Over-Penalization, and Fully Discontinuous Discretizations
BibTeX
@misc{wegert2026discontinuous,
title = {{A discontinuous finite element method for the hydroelastic analysis of submerged structures}},
author = {Zachary J. Wegert and Oriol Colomés and Michael H. Meylan and Ngamta Thamwattana and Vivien J. Challis},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.08313},
eprint = {2609.08313},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.08313},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - A discontinuous finite element method for the hydroelastic analysis of submerged structures
AU - Zachary J. Wegert
AU - Oriol Colomés
AU - Michael H. Meylan
AU - Ngamta Thamwattana
AU - Vivien J. Challis
PY - 2026
DA - 2026/09/10
PB - arXiv
DO - 10.48550/arxiv.2609.08313
UR - https://arxiv.org/abs/2609.08313
AB - This paper develops a finite element method for wave-structure interaction problems arising in the hydroelastic modelling of submerged elastic plates. The approach is formulated for both two-dimensional settings and three-dimensional channels. For the fluid part of the problem, we use a symmetric discontinuous Galerkin scheme that captures the discontinuity in the potential across the plate. To impose the Sommerfeld radiation conditions, we take inspiration from the Dirichlet-to-Neumann map approach and use the analytic solution in the semi-infinite domains to construct appropriate non-local boundary conditions. For the plate, we use a continuous/discontinuous Galerkin method to resolve the 4^th order operator without using continuously differentiable finite elements. We prove well-posedness of the method, showing that the resulting sesquilinear form is bounded and satisfies the Gårding inequality, leading to a unique and stable solution. We implement our formulation as an open-source tool using the Gridap package in Julia. We validate our method against a hypersingular boundary integral method, finding excellent agreement. The flexibility of the methodology is demonstrated through simulations of a two-dimensional slanted plate and three-dimensional plate geometries including centred rectangular plates, offset rectangular plates, and annulus-shaped plates. The finite element method presented here can be readily applied to model wave energy conversion of piezoelectric bimorphs. Furthermore, the method could easily be extended to consider wave scattering from fixed or rigid structures or from variable bottom topography, and is not limited to hydroelastic applications.
ER -