Near shore hydrodynamics has been an important research area dealing with coastal processes. The nearshore coastal region is the region between the shoreline and a fictive offshore limit which usually is defined as the limit where the depth becomes so large that it no longer influences the waves. This spatially limited but highly energetic zone is where water waves shoal, break and transmit energy to the shoreline and are governed by highly dispersive and non-linear effects. An accurate understanding of this phenomena is extremely useful, especially in emergency situations during hurricanes and storms. While the shallow water assumption is valid in regions where the characteristic wavelength exceeds a typical depth by orders of magnitude, Boussinesq-type equations have been used to model near-shore wave motion. Unfortunately these equations are complex system of coupled non-linear and dispersive differential equations that have made the developement of numerical approximations extremely challenging. In this dissertation, a local discontinuous Galerkin method for Boussinesq-Green Naghdi Equations is presented and validated against experimental results. Currently Green-Naghdi equations have many variants. We develop a numerical method in one horizontal dimension for the Green-Naghdi equations based on rotational characteristics in the velocity field. Stability criterion is also established for the linearized Green-Naghdi equations and a careful proof of linear stability of the numerical method is carried out. Verification is done against a linearized standing wave problem in flat bathymetry and h,p (denoted by K in this thesis) error rates are plotted. The numerical method is validated with experimental data from dispersive and non-linear test cases. / text

Storm surge, the pileup of seawater occurring as a result of high surface stresses and strong currents generated by extreme storm events such as hurricanes, is known to cause greater loss of life than these storms' associated winds. For example, inland flooding from the storm surge along the Gulf Coast during Hurricane Katrina killed hundreds of people. Previous storms produced even larger death tolls. Simultaneously, dune, barrier island, and channel erosion taking place during a hurricane leads to the removal of major flow controls, which significantly affects inland inundation. Also, excessive sea bed scouring around pilings can compromise the structural integrity of bridges, levees, piers, and buildings. Modeling these processes requires tightly coupling a bed morphology equation to the shallow water equations (SWE). Discontinuous Galerkin finite element methods (DGFEMs) are a natural choice for modeling this coupled system, given the need to solve these problems on large, complicated, unstructured computational meshes, as well as the desire to implement hp-adaptivity for capturing the dynamic features of the solution. Comprehensive modeling of these processes in the coastal zone presents several challenges and open questions. Most existing hydrodynamic models use a fixed-bed approach; the bottom is not allowed to evolve in response to the fluid motion. With respect to movable-bed models, there is no single, generally accepted mathematical model in use. Numerical challenges include coupling models of processes that exhibit disparate time scales during fair weather, but possibly similar time scales during intense storms. The main goals of this dissertation include implementing a robust, efficient, tightly-coupled morphological model using the local discontinuous Galerkin (LDG) method within the existing Advanced Circulation (ADCIRC) modeling framework, performing systematic code and model verification (using test cases with known solutions, proven convergence rates, or well-documented physical behavior), analyzing the stability and accuracy of the implemented numerical scheme by way of a priori error estimates, and ultimately laying some of the necessary groundwork needed to simultaneously model storm surges and bed morphodynamics during extreme storm events. / text