Galerkin Methods for Parabolic and SCHR{Ö}DINGER Equations with Dynamical Boundary Conditions and Applications to Underwater Acoustics
For researchers in numerical analysis and underwater acoustics, this work provides rigorous error estimates for a class of boundary conditions, but the alternative models are shown to be superior, limiting the practical impact.
The paper develops Galerkin-finite element methods for parabolic and Schrödinger equations with dynamical boundary conditions, proving optimal error estimates in L2 and H1. It applies these methods to underwater acoustics, showing that alternative models without dynamical boundary conditions yield qualitatively better approximations.
In this paper we consider Galerkin-finite element methods that approximate the solutions of initial-boundary-value problems in one space dimension for parabolic and Schrödinger evolution equations with dynamical boundary conditions. Error estimates of optimal rates of convergence in $L^2$ and $H^1$ are proved for the accociated semidiscrete and fully discrete Crank-Nicolson-Galerkin approximations. The problem involving the Schrödinger equation is motivated by considering the standard `parabolic' (paraxial) approximation to the Helmholtz equation, used in underwater acoustics to model long-range sound propagation in the sea, in the specific case of a domain with a rigid bottom of variable topography. This model is contrasted with alternative ones that avoid the dynamical bottom boundary condition and are shown to yield qualitatively better approximations. In the (real) parabolic case, numerical approximations are considered for dynamical boundary conditions of reactive and dissipative type.