Internal wave reflections and transmissions arising from a non‐uniform mesh. Part I: An analysis for the Crank–Nicolson linear finite element scheme

Abstract

AbstractThis is the first of a series of three related papers dealing with some of the consequences of non‐uniform meshes in a numerical model. In this paper the accuracy of the Crank–Nicolson linear finite element scheme, which is applied to the linear shallow water equations, is examined in the context of a single abrupt change in nodal spacing. The (in)accuracy is quantified in terms of reflection and transmission coefficients. An incident wave impinging on the interface between two regions with different nodal spacings is shown to give rise to no reflected waves and two transmitted waves. The analysis is verified using three different wavelengths (2Δx, 4Δx 8Δx) in three ‘hot‐start’ numerical experiments with a mesh expansion factor of 2 and three experiments with a mesh contraction factor of 1/2. An energy flux analysis based on the concept of group velocity shows that energy is conserved across the interface.