Two alternatives for solving hyperbolic boundary value problems of geophysical fluid dynamics

Linear internal waves in inviscid bounded fluids generally give a mathematically ill-posed problem since hyperbolic equations are combined with elliptic boundary conditions. Such problems are difficult to solve. Two new approaches are added to the existing methods: the first solves the two-dimensional spatial wave equation by iteratively adjusting Cauchy data such that boundary conditions are satisfied along a predefined boundary. After specifying the data in this way, solutions can be computed using the d'Alembert formula. The second new approach can numerically solve a wider class of two dimensional linear hyperbolic boundary value problems by using a ‘boundary collocation’ technique. This method gives solutions in the form of a partial sum of analytic functions that are, from a practical point of view, more easy to handle than solutions obtained from characteristics. Collocation points have to be prescribed along certain segments of the boundary but also in the so-called fundamental intervals, regions along the boundary where Cauchy data can be given arbitrarily without over-or under-determining the problem. Three prototypical hyperbolic boundary value problems are solved with this method: the Poincaré, the Telegraph, and the Tricomi boundary value problem. All solutions show boundary-detached internal shear layers, typical for hyperbolic boundary value problems. For the Tricomi problem it is found that the matrix that has to be inverted to find solutions from the collocation approach is ill-conditioned; thus solutions depend on the distribution of the collocation points and need to be regularized.