A Time Splitting Method For The Three Dimensional Shallow Water Equations

The resulting method is mass conservative and efficient on vector and parallel computers. The accuracy, stability, and computational efficiency of this method is demonstrated for wind induced problems ina rectangular basis."

Holl. Zusammenfass.

A wide variety of problems are associated with the flow of shallow water, such as atmospheric flows, tides, storm surges, river and coastal flows, lake flows, tsunamis. Numerical simulation is an effective tool in solving them and a great variety of numerical methods are available. The first part of the book summarizes the basic physics of shallow-water flow needed to use numerical methods under various conditions. The second part gives an overview of possible numerical methods, together with their stability and accuracy properties as well as with an assessment of their performance under various conditions. This enables the reader to select a method for particular applications. Correct treatment of boundary conditions (often neglected) is emphasized. The major part of the book is about two-dimensional shallow-water equations but a discussion of the 3-D form is included. The book is intended for researchers and users of shallow-water models in oceanographic and meteorological institutes, hydraulic engineering and consulting. It also provides a major source of information for applied and numerical mathematicians.

This book, first published in 2002, contains an introduction to hyperbolic partial differential equations and a powerful class of numerical methods for approximating their solution, including both linear problems and nonlinear conservation laws. These equations describe a wide range of wave propagation and transport phenomena arising in nearly every scientific and engineering discipline. Several applications are described in a self-contained manner, along with much of the mathematical theory of hyperbolic problems. High-resolution versions of Godunov's method are developed, in which Riemann problems are solved to determine the local wave structure and limiters are then applied to eliminate numerical oscillations. These methods were originally designed to capture shock waves accurately, but are also useful tools for studying linear wave-propagation problems, particularly in heterogenous material. The methods studied are implemented in the CLAWPACK software package and source code for all the examples presented can be found on the web, along with animations of many of the simulations. This provides an excellent learning environment for understanding wave propagation phenomena and finite volume methods.

The International Workshop on "The Use of Supercomputers in Theoretical Science" took place on January 24 and 25, 1991, at the University of Antwerp (UIA), Antwerpen, Belgium. It was the sixth in a series of workshops, the fIrst of which took place in 1984. The principal aim of these workshops is to present the state of the art in scientific large-scale and high speed-computation. Computational science has developed into a third methodology equally important now as its theoretical and experimental companions. Gradually academic researchers acquired access to a variety of supercomputers and as a consequence computational science has become a major tool for their work. It is a pleasure to thank the Belgian National Science Foundation (NFWO-FNRS) and the Ministry of ScientifIc Affairs for sponsoring the workshop. It was organized both in the framework of the Third Cycle "Vectorization, Parallel Processing and Supercomputers" and the "Governemental Program in Information Technology". We also very much would like to thank the University of Antwerp (Universitaire Instelling Antwerpen -VIA) for fInancial and material support. Special thanks are due to Mrs. H. Evans for the typing and editing of the manuscripts and for the preparation of the author and subject indexes. J.T. Devreese P.E. Van Camp University of Antwerp July 1991 v CONlENTS High Perfonnance Numerically Intensive Applications on Distributed Memory Parallel Computers .................... . F.W. Wray Abstract ......................................... .

This thesis concerns the use of time-split methods for the numerical solution of time-dependent partial differential equations. Frequently the differential operator splits additively into two or more pieces such that the corresponding subproblems are each easier to solve than the original equation, or are best handled by different techniques. In the time-split method the solution to the original equation is advanced by alternately solving the subproblems. In this thesis a unified approach to splitting methods is developed which simplifies their analysis. Particular emphasis is given to splittings of hyperbolic problems into subproblems with disparate wave speeds. Three main aspects of the method are considered. The first is the accuracy and efficiency of the time-split method relative to unsplit methods. The second topic is stability for split methods. The final topic is the proper specification of boundary data for the intermediate solutions, e.g., the solution obtained after solving only one of the subproblems. The main emphasis is on hyperbolic problems, and the one-dimensional shallow water equations are used as a specific example throughout. The final chapter is devoted to some other applications or the theory. Two-dimensional hyperbolic problems, convection-diffusion equations, and the Peaceman-Rachford ADI method for the heat equation are considered.