Numerical solution of the two-dimensional Poincaré equation

This paper deals with numerical approximation of the two-dimensional Poincaré equation that arises as a model for internal wave motion in enclosed containers. Inspired by the hyperbolicity of the equation we propose a discretisation particularly suited for this problem, which results in matrices whose size varies linearly with the number of grid points along the coordinate axes. Exact solutions are obtained, defined on a perturbed boundary. Furthermore, the problem is seen to be ill-posed and there is need for a regularisation scheme, which we base on a minimal-energy approach.