Abstract
A new numerical method for Nwogu’s (ASCE Journal of Waterway, Port, Coastal and Ocean En-gineering1993;119:618) two-dimensional extended Boussinesq equations is presented using a lineartriangular nite element spatial discretization coupled with a sophisticated adaptive time integrationpackage. The authors have previously presented a nite element method for the one-dimensional formof these equations (M. Walkley and M. Berzins (International Journal for Numerical Methods inFluids1999;29(2):143)) and this paper describes the extension of these ideas to the two-dimensionalequations and the application of the method to complex geometries using unstructured triangular grids.Computational results are presented for two standard test problems and a realistic harbour model. Copy-right?2002 John Wiley & Sons, Ltd