SWASH: An operational public domain code for simulating wave fields and rapidly varied flows in coas

Abstract A computational procedure has been developed for simulating non-hydrostatic, free-surface, rotational flows in one and two horizontal dimensions. Its implementation in the publicly available SWASH (an acronym of Simulating WAves till SHore) is intended to be used for predicting transformation of surface waves and rapidly varied shallow water flows in coastal waters. This open source code ( http://swash.sourceforge.net ) has been developed based on the work of Stelling and Zijlema, 2003 , Stelling and Duinmeijer, 2003 , Zijlema and Stelling, 2005 , Zijlema and Stelling, 2008 ). The governing equations are the nonlinear shallow water equations including non-hydrostatic pressure and provide a general basis for describing complex changes to rapidly varied flows typically found in coastal flooding resulting from e.g. dike breaks and tsunamis, and wave transformation in both surf and swash zones due to nonlinear wave–wave interactions, interaction of waves with currents, and wave breaking as well as runup at the shoreline. The present paper provides a complete description of the numerical algorithms currently used in the code. The code is benchmarked using some analytical problems. Moreover, the numerical results are validated with various cases of laboratory data with the principal aim to convey the capabilities of the SWASH code. In particular, emphasis is put on an analysis of model performance and associated physical implications. Serial and parallel performance scalings are also presented.

[1]  R. B. Simpson,et al.  Centered differencing and the box scheme for diffusion convection problems , 1976 .

[2]  Guus S. Stelling,et al.  NUMERICAL SIMULATION OF 3D QUASI-HYDROSTATIC, FREE-SURFACE FLOWS , 1998 .

[3]  A. C. Radder,et al.  Verification of numerical wave propagation models for simple harmonic linear water waves , 1982 .

[4]  Michel Benoit,et al.  Non-linear propagation of unidirectional wave fields over varying topography , 1999 .

[5]  I. A. Svendsen,et al.  Quasi-three-dimensional modeling of rip current systems , 2003 .

[6]  Harry Yeh,et al.  Tsunami run-up and draw-down on a plane beach , 2003, Journal of Fluid Mechanics.

[7]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[8]  Franklin Liu,et al.  Modeling wave runup with depth-integrated equations , 2002 .

[9]  J. Kirby,et al.  Boussinesq modeling of a rip current system , 1999 .

[10]  N. Booij,et al.  A third-generation wave model for coastal regions-1 , 1999 .

[11]  Kwok Fai Cheung,et al.  Depth‐integrated, non‐hydrostatic model with grid nesting for tsunami generation, propagation, and run‐up , 2010 .

[12]  Michael J. Briggs,et al.  Laboratory experiments of tsunami runup on a circular island , 1995 .

[13]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[14]  Kwok Fai Cheung,et al.  Well-Balanced Finite-Volume Model for Long-Wave Runup , 2006 .

[15]  Peter Stansby,et al.  Shallow‐water flow solver with non‐hydrostatic pressure: 2D vertical plane problems , 1998 .

[16]  J. Kan A second-order accurate pressure correction scheme for viscous incompressible flow , 1986 .

[17]  S. Elgar Waves in Oceanic and Coastal Waters , 2007 .

[18]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[19]  Robert A. Dalrymple,et al.  Experimental study of nearshore dynamics on a barred beach with rip channels , 2002 .

[20]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[21]  Benedict D. Rogers,et al.  Numerical Modeling of Water Waves with the SPH Method , 2006 .

[22]  Marcel Zijlema,et al.  An accurate and efficient finite‐difference algorithm for non‐hydrostatic free‐surface flow with application to wave propagation , 2003 .

[23]  H. L. Stone ITERATIVE SOLUTION OF IMPLICIT APPROXIMATIONS OF MULTIDIMENSIONAL PARTIAL DIFFERENTIAL EQUATIONS , 1968 .

[24]  J. A. Battjes,et al.  ENERGY LOSS AND SET-UP DUE TO BREAKING OF RANDOM WAVES , 1978 .

[25]  D. Causon,et al.  The surface gradient method for the treatment of source terms in the shallow-water equations , 2001 .

[26]  G. Stelling,et al.  On the construction of computational methods for shallow water flow problems , 1983 .

[27]  Guus S. Stelling,et al.  A staggered conservative scheme for every Froude number in rapidly varied shallow water flows , 2003 .

[28]  Peter Stansby,et al.  2D shallow water flow model for the hydraulic jump , 1999 .

[29]  E. Blayo,et al.  Revisiting open boundary conditions from the point of view of characteristic variables , 2005 .

[30]  M. Boers,et al.  Simulation of a surf zone with a barred beach. Report 1. Wave heights and wave breaking , 1996 .

[31]  W. Thacker Some exact solutions to the nonlinear shallow-water wave equations , 1981, Journal of Fluid Mechanics.

[32]  Per A. Madsen,et al.  Simulation of nonlinear wave run-up with a high-order Boussinesq model , 2008 .

[33]  Marcel Zijlema,et al.  Further experiences with computing non‐hydrostatic free‐surface flows involving water waves , 2005 .

[34]  Akihide Tada,et al.  Applicability of numerical models to nonlinear dispersive waves , 1995 .

[35]  Guus S. Stelling,et al.  A finite volume analogue of the P1NC-P1 finite element: With accurate flooding and drying , 2010 .

[36]  Marcel Zijlema,et al.  Efficient computation of surf zone waves using the nonlinear shallow water equations with non-hydrostatic pressure , 2008 .

[37]  Kwok Fai Cheung,et al.  Depth‐integrated, non‐hydrostatic model for wave breaking and run‐up , 2009 .

[38]  Chih-Chieh Young,et al.  A σ‐coordinate non‐hydrostatic model with embedded Boussinesq‐type‐like equations for modeling deep‐water waves , 2009 .

[39]  Jurjen A. Battjes,et al.  Experimental investigation of wave propagation over a bar , 1993 .

[40]  Alfredo Bermúdez,et al.  Upwind methods for hyperbolic conservation laws with source terms , 1994 .

[41]  W. Hansen Theorie zur Errechnung des Wasserstandes und der Strömungen in Randmeeren nebst Anwendungen , 1956 .

[42]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[43]  Congfang Ai,et al.  A new fully non‐hydrostatic 3D free surface flow model for water wave motions , 2011 .

[44]  Stephan C. Kramer,et al.  A conservative unstructured scheme for rapidly varied flows , 2008 .