A mass, energy, vorticity, and potential enstrophy conserving lateral fluid–land boundary scheme for the shallow water equations

A numerical scheme for treating fluid–land boundaries in inviscid shallow water flows is derived that conserves the domain-summed mass, energy, vorticity, and potential enstrophy in domains with arbitrarily shaped boundaries. The boundary scheme is derived from a previous scheme that conserves all four domain-summed quantities only in periodic domains without boundaries. It consists of a method for including land in the model along with evolution equations for the vorticity and extrapolation formulas for the depth at fluid–land boundaries. Proofs of mass, energy, vorticity, and potential enstrophy conservation are given. Numerical simulations are carried out demonstrating the conservation properties and accuracy of the boundary scheme for inviscid flows and comparing its performance with that of four alternative boundary schemes. The first of these alternatives extrapolates or finite-differences the velocity to obtain the vorticity at boundaries; the second enforces the free-slip boundary condition; the third enforces the super-slip condition; and the fourth enforces the no-slip condition. Comparisons of the conservation properties demonstrate that the new scheme is the only one of the five that conserves all four domain-summed quantities, and it is the only one that both prevents a spurious energy cascade to the smallest resolved scales and maintains the correct flow orientation with respect to an external forcing. Comparisons of the accuracy demonstrate that the new scheme generates vorticity fields that have smaller errors than those generated by any of the alternative schemes, and it generates depth and velocity fields that have errors about equal to those in the fields generated by the most accurate alternative scheme.