Shallow-water flow solver with non-hydrostatic pressure: 2D vertical plane problems

A numerical solution for shallow-water flow is developed based on the unsteady Reynolds-averaged Navier-Stokes equations without the conventional assumption of hydrostatic pressure. Instead, the non-hydrostatic pressure component may be added in regions where its influence is significant, notably where bed slope is not small and separation in a vertical plane may occur or where the free-surface slope is not small. The equations are solved in the sigma-co-ordinate system with semi-implicit time stepping and the eddy viscosity is calculated using the standard k-epsilon turbulence model. Conventionally, boundary conditions at the bed for shallow-water models only include vertical diffusion terms using wall functions, but here they are extended to include horizontal diffusion terms which can be significant when bed slope is not small. This is consistent with the inclusion of non-hydrostatic pressure. The model is applied to the 2D vertical plane flow of a current over a trench for which experimental data and other numerical results are available for comparison. Computations with and without non-hydrostatic pressure are compared for the same trench and for trenches with smaller side slopes, to test the range of validity of the conventional hydrostatic pressure assumption. The model is then applied to flow over a 2D mound and again the slope of the mound is reduced to assess the validity of the hydrostatic pressure assumption. (C) 1998 John Wiley & Sons, Ltd.