Elliptic harbor wave model with perfectly matched layer and exterior bathymetry effects

Standard strategies for dealing with the Sommerfeld condition in elliptic mild-slope models require strong assumptions on the wave field in the region exterior to the computational domain. More precisely, constant bathymetry along (and beyond) the open boundary, and parabolic approximations-based boundary conditions are usually imposed. Generally, these restrictions require large computational domains, implying higher costs for the numerical solver. An alternative method for coastal/harbor applications is proposed here. This approach is based on a perfectly matched layer (PML) that incorporates the effects of the exterior bathymetry. The model only requires constant exterior depth in the alongshore direction, a common approach used for idealizing the exterior bathymetry in elliptic models. In opposition to standard open boundary conditions for mild-slope models, the features of the proposed PML approach include (1) completely noncollinear coastlines, (2) better representation of the real unbounded domain using two different lateral sections to define the exterior bathymetry, and (3) the generation of reliable solutions for any incoming wave direction in a small computational domain. Numerical results of synthetic tests demonstrate that solutions are not significantly perturbed when open boundaries are placed close to the area of interest. In more complex problems, this provides important performance improvements in computational time, as shown for a real application of harbor agitation.