Analysis Of The Numerical Solution Of The Shallow Water Equations

This thesis is concerned with the analysis of various methods for the numerical solution of the shallow water equations along with the stability of these methods. Most of the thesis is concerned with the background and formaulation of the shallow water equations. The derivation of the basic equations will be given, in the primative variable and vorticity divergence formulation. Also the shallow water equations will be written in spherical coordinates. Two main types of methods used in approximating differential equations of this nature will be discussed. The two schemes are finite difference method (FDM) and the finite element method (FEM). After presenting the shallow water equations in several formulations, some examples will be presented. The use of the Fourier transform to find the solution of a semidiscrete analog of the shallow water equations is also demonstrated.

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.

Within this monograph a comprehensive and systematic knowledge on shallow-water hydrodynamics is presented. A two-dimensional system of shallow-water equations is analyzed, including the mathematical and mechanical backgrounds, the properties of the system and its solution. Also featured is a new mathematical simulation of shallow-water flows by compressible plane flows of a special virtual perfect gas, as well as practical algorithms such as FDM, FEM, and FVM. Some of these algorithms have been utilized in solving the system, while others have been utilized in various applied fields. An emphasis has been placed on several classes of high-performance difference schemes and boundary procedures which have found wide uses recently for solving the Euler equations of gas dynamics in aeronautical and aerospatial engineering. This book is constructed so that it may serve as a handbook for practicians. It will be of interest to scientists, designers, teachers, postgraduates and professionals in hydraulic, marine, and environmental engineering; especially those involved in the mathematical modelling of shallow-water bodies.

1. 1 AREAS OF APPLICATION FOR THE SHALLOW WATER EQUATIONS The shallow water equations describe conservation of mass and mo mentum in a fluid. They may be expressed in the primitive equation form Continuity Equation _ a, + V. (Hv) = 0 L(l;,v;h) at (1. 1) Non-Conservative Momentum Equations a M("vjt,f,g,h,A) = at(v) + (v. V)v + tv - fkxv + gV, - AIH = 0 (1. 2) 2 where is elevation above a datum (L) ~ h is bathymetry (L) H = h + C is total fluid depth (L) v is vertically averaged fluid velocity in eastward direction (x) and northward direction (y) (LIT) t is the non-linear friction coefficient (liT) f is the Coriolis parameter (liT) is acceleration due to gravity (L/T2) g A is atmospheric (wind) forcing in eastward direction (x) and northward direction (y) (L2/T2) v is the gradient operator (IlL) k is a unit vector in the vertical direction (1) x is positive eastward (L) is positive northward (L) Y t is time (T) These Non-Conservative Momentum Equations may be compared to the Conservative Momentum Equations (2. 4). The latter originate directly from a vertical integration of a momentum balance over a fluid ele ment. The former are obtained indirectly, through subtraction of the continuity equation from the latter. Equations (1. 1) and (1. 2) are valid under the following assumptions: 1. The fluid is well-mixed vertically with a hydrostatic pressure gradient. 2. The density of the fluid is constant.

This book presents the theory and computation of open channel flows, using detailed analytical, numerical and experimental results. The fundamental equations of open channel flows are derived by means of a rigorous vertical integration of the RANS equations for turbulent flow. In turn, the hydrostatic pressure hypothesis, which forms the core of many shallow water hydraulic models, is scrutinized by analyzing its underlying assumptions. The book’s main focus is on one-dimensional models, including detailed treatments of unsteady and steady flows. The use of modern shock capturing finite difference and finite volume methods is described in detail, and the quality of solutions is carefully assessed on the basis of analytical and experimental results. The book’s unique features include: • Rigorous derivation of the hydrostatic-based shallow water hydraulic models • Detailed treatment of steady open channel flows, including the computation of transcritical flow profiles • General analysis of gate maneuvers as the solution of a Riemann problem • Presents modern shock capturing finite volume methods for the computation of unsteady free surface flows • Introduces readers to movable bed and sediment transport in shallow water models • Includes numerical solutions of shallow water hydraulic models for non-hydrostatic steady and unsteady free surface flows This book is suitable for both undergraduate and graduate level students, given that the theory and numerical methods are progressively introduced starting with the basics. As supporting material, a collection of source codes written in Visual Basic and inserted as macros in Microsoft Excel® is available. The theory is implemented step-by-step in the codes, and the resulting programs are used throughout the book to produce the respective solutions.

The objective of this book is to analyze within reasonable limits (it is not a treatise) the basic mathematical aspects of the finite element method. The book should also serve as an introduction to current research on this subject. On the one hand, it is also intended to be a working textbook for advanced courses in Numerical Analysis, as typically taught in graduate courses in American and French universities. For example, it is the author’s experience that a one-semester course (on a three-hour per week basis) can be taught from Chapters 1, 2 and 3 (with the exception of Section 3.3), while another one-semester course can be taught from Chapters 4 and 6. On the other hand, it is hoped that this book will prove to be useful for researchers interested in advanced aspects of the numerical analysis of the finite element method. In this respect, Section 3.3, Chapters 5, 7 and 8, and the sections on “Additional Bibliography and Comments should provide many suggestions for conducting seminars.