Abstract:We present an energy/entropy stable and high order accurate finite difference (FD) method for solving the nonlinear (rotating) shallow water equations (SWEs) in vector invariant form using the newly developed dual-pairing and dispersion-relation preserving summation by parts (SBP) FD operators. We derive new well-posed boundary conditions (BCs) for the SWE in one space dimension, formulated in terms of fluxes and applicable to linear and nonlinear SWEs. For the nonlinear vector invariant SWE in the subcritical regime, where energy is an entropy functional, we find that energy/entropy stability ensures the boundedness of numerical solution but does not guarantee convergence. Adequate amount of numerical dissipation is necessary to control high frequency errors which could negatively impact accuracy in the numerical simulations. Using the dual-pairing SBP framework, we derive high order accurate and nonlinear hyper-viscosity operator which dissipates entropy and enstrophy. The hyper-viscosity operator effectively minimises oscillations from shocks and discontinuities, and eliminates high frequency grid-scale errors. The numerical method is most suitable for the simulations of subcritical flows typically observed in atmospheric and geostrophic flow problems. We prove both nonlinear and local linear stability results, as well as a priori error estimates for the semi-discrete approximations of both linear and nonlinear SWEs. Convergence, accuracy, and well-balanced properties are verified via the method of manufactured solutions and canonical test problems such as the dam break and lake at rest. Numerical simulations in two-dimensions are presented which include the rotating and merging vortex problem and barotropic shear instability, with fully developed turbulence.

A higher order time domain rankine panel method for linear and weakly non-linear computation

Finite Element Analysis of Second Order Wave Radiation by a Group of Cylinders in the Time Domain

Numerical and Experimental Study on A 2-D Floating Body under Extreme Wave Conditions

Development of time-domain Rankine source method for multiple floating bodies with large plane displacement

A unified steady and unsteady formulation for hydrodynamic potential flow simulations with fully nonlinear free surface boundary conditions

Three dimensional numerical modelling for wave radiation problem under arbitrary seabed condition

A High-Order Hybrid-Spectral Incompressible Navier-Stokes Model For Nonlinear Water Waves

Towards High-Performance Linear Potential Flow BEM Solver with Low-Rank Compressions

High order finite volume schemes with IMEX time stepping for the Boltzmann model on unstructured meshes

Immersed boundary method for high-order flux reconstruction based on volume penalization

A fictitious domain approach based on a viscosity penalty method to simulate wave/structure interaction

An efficient hybrid Laplace-time domain method for dynamic analysis of nonlinear floating systems

High-Resolution Viscous Terms Discretization and ILW Solid Wall Boundary Treatment for the Navier–Stokes Equations

Very high order treatment of embedded curved boundaries in compressible flows: ADER discontinuous Galerkin with a space-time Reconstruction for Off-site data

3D WC-MPS coupled with geometrically nonlinear shell for hydro-elastic free-surface flows

Strongly stable dual-pairing summation by parts finite difference schemes for the vector invariant nonlinear shallow water equations -- I: Numerical scheme and validation on the plane

Numerical Simulation of Floating Bodies in Extreme Free Surface Waves

An Open-Source Python-Based Boundary-Element Method Code for the Three-Dimensional, Zero-Froude, Infinite-Depth, Water-Wave Diffraction-Radiation Problem

On the robustness of high-order upwind summation-by-parts methods for nonlinear conservation laws

Wave interaction with multiple adjacent floating solar panels with arbitrary constraints