A two-dimensional coastal ocean model based on unstructured C-grid is built, in which the momentum equation is discretized on the faces of each cell, and the continuity equation is discretized on the cell. The model is discretized by semi-implicit finite volume method, in that the free surface is semi-implicit and the bottom friction is implicit, thereby removing stability limitations associated with the surface gravity wave and friction. The remaining terms in the momentum equations are discretized explicitly by integral finite volume method and second-order Adams-Bashforth method. Tidal flow in the polar quadrant with known analytic solution is employed to test the proposed model. Finally, the performance of the present model to simulate tidal flow in a geometrically complex domain is examined by simulation of tidal currents in the Pearl River Estuary.
In this study, porosity was introduced into two-dimensional shallow water equations to reflect the effects of obstructions, leading to the modification of the expressions for the flux and source terms. An extra porosity source term appears in the momentum equation. The numerical model of the shallow water equations with porosity is presented with the finite volume method on unstructured grids and the modified Roe-type approximate Riemann solver. The source terms of the bed slope and porosity are both decomposed in the characteristic direction so that the numerical scheme can exactly satisfy the conservative property. The present model was tested with a dam break with discontinuous porosity and a flash flood in the Toce River Valley. The results show that the model can simulate the influence of obstructions, and the numerical scheme can maintain the flux balance at the interface with high efficiency and resolution.
In this paper, the parametric tropical cyclone models for storm surge modeling are further developed. Instead of tangential wind speed via cyclostrophic balance and radial wind speed using a simple formulation of defection angle, the analyrical expressions of tangential and radial wind speed distribution are derived from the governing momentum equations based on the general symmetric pressure distribution of Holland and Fujita. The radius of the maximum wind is estimated by tropical cyclone wind structure which is characterized by the radial extent of special wind speed. The shape parameter in the pressure model is estimated by the data of several tropical cyclones that occurred in the East China Sea. Finally, the Fred cyclone (typhoon 199417) is calculated, and comparisons of the measured and calculated air pressures and wind speed are presented.