An adaptive mesh refinement algorithm for shallow water equations is presented. The algorithm uses upwind scheme that is Godunov type and which approximately solves the Riemann problem using Roe's technique. A highly accurate solution is achieved by using the adaptive mesh refinement technique of Berger and Oliger for mesh refinement algorithm. The numerical method is second-order accurate and approximately max-min preserving by using van Leer limited-slope technique. One-dimensional nesting algorithm has been implemented successfully. Numerical results on a test problem verify the second order accuracy of the algorithm. The nested grid results yield the equivalent solution to that of the corresponding fine grid solution.

acase@tulane.edu / Shallow water models are widely used to describe and study fluid dynamics phenomena where the horizontal length scale is much greater than the vertical length scale, for example, in the atmosphere and oceans. Since analytical solutions of the shallow water models are typically out of reach, development of accurate and efficient numerical methods is crucial to understand many mechanisms of atmospheric and oceanic phenomena. In this dissertation, we are interested in developing simple, accurate, efficient and robust numerical methods for two shallow water models --- the Saint-Venant system of shallow water equations and the two-mode shallow water equations. We first construct a new second-order moving-water equilibria preserving central-upwind scheme for the Saint-Venant system of shallow water equations. Special reconstruction procedure and source term discretization are the key components that guarantee the resulting scheme is capable of exactly preserving smooth moving-water steady-state solutions and a draining time-step technique ensures positivity of the water depth. Several numerical experiments are performed to verify the well-balanced and positivity preserving properties as well as the ability of the proposed scheme to accurately capture small perturbations of moving-water steady states. We also demonstrate the advantage and importance of utilizing the new method over its still-water equilibria preserving counterpart. We then develop and study numerical methods for the two-mode shallow water equations in a systematic way. Designing a reliable numerical method for this system is a challenging task due to its conditional hyperbolicity and the presence of nonconservative terms. We present several numerical approaches---two operator splitting methods (based on either Roe-type upwind or central-upwind scheme), a central-upwind scheme and a path-conservative central-upwind scheme---and test their performance in a number of numerical experiments. The obtained results demonstrate that a careful numerical treatment of nonconservative terms is crucial for designing a robust and highly accurate numerical method for this system. / 1 / Yuanzhen Cheng

No description available.

This thesis presents different control design approaches for stabilizing networks of quasi-linear hyperbolic partial differential equations. These equations are usually conservative which gives them interesting properties to design stabilizing control laws. Two main design approaches are developed: a methodology based on entropies and Lyapunov functions and a methodology based on the Riemann invariants. The stability theorems are illustrated using numerical simulations. Two practical applications of these methodologies are presented. Netword of navigation channels are modelled using Saint-Venant equations (also known as the Shallow Water Equations). The stabilization problem of such system has an industrial importance in order to satisfy the navigation constraints and to optimize the production of electricity in hydroelectric plants, usually located at each hydraulic gates. A second application deals with the regulation of water waves in moving tanks. This problem is also modelled by a modified version of the shallow water equations and appears in a number of industrial fields which deal with liquid moving parts.

Two kinds of tsunami models are used in this thesis to simulate tsunami propagation in the ocean. One is the linear dispersion tsunami model developed by Port and Airport Research Institute (PARI), Japan. The other is COrnell Multigrid COupled Tsunami model (COMCOT) developed by the School of Civil and Environmental Engineering, Cornell University, that carries on the tsunami run-up computation to the nearshore region. Two kinds of tsunami models have the same mechanism of initial wave profile, which is the vertical seabed displacement as the initial tsunami profile proposed by Mansinha and Smylie (1971). Both models describe the tsunami by the same shallow water equations. At first, the feasibility of the PARI model is established by comparing with the record in Maldives during the South Asia tsunami in December 2004. Then, the COMCOT model in applied to the Pingtong earthquake in December 2006 and is validated by comparing with the tidal station records. Possible submarine fault activities around Taiwan and the Western Pacific ring is simulated by the The PARI model based on moment magnitude scale (¢Ûw) close to the South Asia tsunami. Seven sources are chosen: the Hokkaido, the East Japan, the Ryukyu Islands, the Guishan Island in Taiwan, the Fukien of mainland China, the Luzon Island and the New Guinea. The results suggest the northeast and southwest part of Taiwan have potential tsunami risk. Finally, we simulate the fault activity between Taiwan and Luzon islands by the COMCOT model. The inundation area extends northward to the Tso-Ying and San-Min districts, eastward to the Siao-Gang district and Fengshan city. The Kaohsiung harbor can resist tsunami hazard for moment magnitude scale (¢Ûw) up to 7.58 with maximum wave height of 5.5 meters.

In this thesis a high order finite difference scheme is derived and implemented solving the shallow water equations using the SBP-SAT method. This method was tested against various benchmark problems were convergence was verified. The shallow water equations were also solved on a multi-block setup representing a tsunami approaching a shoreline from the ocean. Experiments show that a bottom topography with many spikes provides a dispersing effect on the incoming tsunami wave. Higher order convergence is not guaranteed for the multi-block simulations and could be investigated further in a future study.

Flood propagation over urban areas can cause an interaction between the free-surface flow and large underground pipe networks used for storm drainage and sewage, causing outflows and inflows at the bed. The associated waves may collide with each other and the surface waves. In this thesis the shallow water equations are used to model this type of wave interaction over dry or wet beds with bathymetry gradients and friction terms. The proposed shallow water scheme is solved based on finite volume high-resolution Godunov-type methods. The solver is well-balanced and can accurately balance the source terms and flux-gradients for the steady-state solutions. The solver also utilises a new type of Riemann wave speed to provide depth-positive results over nearly dry beds and dry states. Additionally a new type of source term is introduced in the continuity equation to model pipe inflow and outflow conditions at bed connections. For the standard one-dimensional shallow water equations the numerical results are validated with analytical solutions or other reference solutions provided in the literature. This includes the incipient Riemann problems for nearly dry and dry-states, steady flow over a hump in a rectangular channel and the wave propagation problem. Eventually, the generation of dry bed in the middle, over discontinuous topography is considered. Close agreement is achieved between the shallow water scheme and analytical or reference solutions for the above test cases. For the shallow water problems with influx/efflux source terms comparisons are made with STAR-CD, a commercial Navier-Stokes solver for general fluid flow prediction. The shallow water model is first used to simulate vertical flows through finite gaps in the bed. Next, the interaction of the vertical flows with a dam-break flow is considered for both dry and wet beds. An efflux number, En, is defined based on the vertical efflux velocity and the gap length. A parameter study is undertaken to investigate the effect of the one-dimensional approximation of the present model, for a range of non-dimensional efflux numbers. It is found that the shallow flow model gives sensible predictions at all times provided En<0.5, and for long durations for En>0.5. Dam break flow over an underground connecting pipe is also considered for the one-dimensional efflux problems. To solve two-dimensional problems the shallow water scheme uses the dimensional-splitting method which solves each one-dimensional Riemann problem in the x- and y-directions separately. The cross-derivative terms for second-order accuracy are incorporated by solving another Riemann problem in the orthogonal direction. For two-dimensional problems first the dam-break problems are considered over wet and dry beds. Then, flood propagation over complex terrain is demonstrated. Next, efflux discharge is modelled in isolation over a dry bed and then with dam-break interaction, comparing with STAR-CD results. Again very good agreement is shown between the two-dimensional shallow water model and STAR-CD for the efflux numbers of En<0.5. For modelling the inundation problem over an underground pipe network the solver is coupled with the general underground pipe network solver to calculate the efflux discharge as the flood waves pass through the pipe network. For analysing the pipe network with unknown effluxes an additional set of equations is incorporated into the solution of a general pipe network solver. The shallow water solver coupled to an underground pipe network is then used to simulate dam-break interaction with pipe networks with 9 and 25 nodes to demonstrate the versatility of the method.

Accurate prediction of tsunami waves requires complete boundary and initial condition data, coupled with the appropriate mathematical model. However, necessary data is often missing or inaccurate, and may not have sufficient resolution to capture the dynamics of such nonlinear waves accurately. In this thesis we demonstrate that variational data assimilation for the continuous shallow water equations (SWE) is a feasible approach for recovering both initial conditions and bathymetry data from sparse observations. Using a Sadourny finite-difference finite volume discretisation for our numerical implementation, we show that convergence to true initial conditions can be achieved for sparse observations arranged in multiple configurations, for both isotropic and anisotropic initial conditions, and with realistic bathymetry data in two dimensions. We demonstrate that for the 1-D SWE, convergence to exact bathymetry is improved by including a low-pass filter in the data assimilation algorithm designed to remove scale-scale noise, and with a larger number of observations. A necessary condition for a relative L2 error less than 10% in bathymetry reconstruction is that the amplitude of the initial conditions be less than 1% of the bathymetry height. We perform Second Order Adjoint Sensitivity Analysis and Global Sensitivity Analysis to comprehensively assess the sensitivity of the surface wave to errors in the bathymetry and perturbations in the observations. By demonstrating low sensitivity of the surface wave to the reconstruction error, we found that reconstructing the bathymetry with a relative error of about 10% is sufficiently accurate for surface wave modelling in most cases. These idealised results with simplified 2-D and 1-D geometry are intended to be a first step towards more physically realistic settings, and can be used in tsunami modelling to (i) maximise accuracy of tsunami prediction through sufficiently accurate reconstruction of the necessary data, (ii) attain a priori knowledge of how different bathymetry and initial conditions can affect the surface wave error, and (iii) provide insight on how these can be mitigated through optimal configuration of the observations. / Thesis / Candidate in Philosophy

The research detailed in this study focuses on the determination of an appropriate objective function to aid parameter estimation when simulating areas influenced by tidally varying flows. Three objective functions that are measures of how well the model results match field data at several locations and times were tested. A set of test cases is developed to represent tidally influenced systems and allow for the testing of the objective functions. These objective functions were tested by computing their values and comparing them for the various estimated parameters. Based on results of the first method of testing a further analysis was performed using PEST, an automatic parameter estimation tool. A weighted least squares of the velocity and water surface values with a weight function on the velocity term based on the shallow water equations is found to be a reasonable objective function at this point in the research.

Riemann problems are of prime importance in computational fluid dynamics simulations using finite elements or finite volumes discretizations. In some applications, billions of Riemann problems might need to be solved in a single simulation, therefore it is important to have reliable and computationally efficient algorithms to do so. Given the nonlinearity of the flux function in most systems considered in practice, to obtain an exact solution for the Riemann problem explicitly is often not possible, and iterative solvers are required. However, because of issues found with existing iterative solvers like lack of convergence and high computational cost, their use is avoided and approximate solvers are preferred. In this thesis work, motivated by the advances in computer hardware and algorithms in the last years, we revisit the possibility of using iterative solvers to compute the exact solution for Riemann problems. In particular, we focus on the development, implementation, and performance comparison of iterative Riemann solvers for the shallow water and Euler equations. In a one-dimensional homogeneous framework for these systems, we consider several initial guesses and iterative methods for the computation of the Riemann solution. We find that efficient and reliable iterative solvers can be obtained by using recent estimates on the Riemann solution to modify and combine well-known methods. Finally, we consider the application of these solvers in finite volume simulations using the wave propagation algorithms implemented in Clawpack.