Perfectly Matched Layers Methods for Mixed Hyperbolic-Dispersive Equations

Absorbing boundary conditions are important when one simulates the propagation of waves on a bounded numerical domain without creating artificial reflections. In this paper, we consider various hyperbolic-dispersive equations modeling water wave propagation. A typical example is the Korteweg-de Vries equation u(t) + u u(x) + epsilon mu(xxx) = 0, for all x is an element of R, for all t > 0. (1) In the case of linearized equations, some progress was recently done for one dimensional scalar dispersive equations using discrete transparent boundary conditions. However, a generalization of this approach to multi-dimensional setting is not obvious. In this paper, we consider the alternative perfectly matched layer (PML) approach for the linearized Korteweg-de Vries equation: u(t) + U u(x) + epsilon mu(xxx) = 0, for all x is an element of R, for all t > 0, (2) where U is an element of R denotes a reference speed. We first propose a direct perfectly matched layer approach and study the stability of the modified system. These equations are not always stable, the main obstruction being the classical condition v(g)(k)v(phi)(k) >= 0 found in the literature onPML(Becache et al in J Comput Phys 188(2):399-433, 2003) that we recover in our analysis. Then, we introduce a hyperbolic system with a source term that is an approximation of the Korteweg-de Vries equations. In this case, the complete PML equations are not, again, completely stable. However, a version of the PML equations for this system derived without the source term is found to be stable and can absorb outgoing waves although it may create reflections as it is not perfectly matched. Finally, we consider the BBM-Boussinesq system that models bi-directional waves at the surface of an inviscid fluid layer. The dispersive properties for a subclass of physically relevant models are better suited for PML techniques since the condition v(g)(k)v(phi)(k) >= 0 is always satisfied. We show that the PML equations are always stable in this case. We illustrate numerically the absorbing and stability properties of these PML models and provide also KdV type simulation by choosing properly initial data.