A 2-D Numerical Model for Linear Long Wave Propagation In Boundary-Fitted Curvilinear Grids (original) (raw)
Related papers
Fully dispersive nonlinear water wave model in curvilinear coordinates
Journal of Computational Physics, 2004
A vertically integrated fully dispersive nonlinear wave model is expressed in curvilinear coordinates with non-orthogonal grids for the simulation of broad-banded nonlinear random water waves in regions of arbitrary geometry. The transformation is performed for both dependent and independent variables, hence an irregular physical domain is converted into a rectangular computational domain with contravariant velocities. Use of contravariant velocity components as dependent variables ensures easy and accurate satisfaction of the wall condition for lateral enclosures surrounding a physical domain, such as a coastal area, channel, or harbor. The numerical scheme is based on finitedifference approximations with staggered grids which results in implicit formulations for the momentum equations and a semi-explicit formulation for the continuity equation. Linear long wave propagation in a channel of varying crosssection and linear random wave propagation in a circular channel are presented as test cases for comparisons with the corresponding analytical solutions. Cnoidal and Stokes waves in a circular channel are also simulated as examples to nonlinear wave propagation within curved walls.
Numerical wave propagation on the hexagonal C-grid
Journal of Computational Physics, 2008
Inertio-gravity mode and Rossby mode dispersion properties are examined for discretizations of the linearized rotating shallow water equations on a regular hexagonal C-grid in planar geometry. It is shown that spurious non-zero Rossby mode frequencies found by previous authors in the f-plane case can be avoided by an appropriate discretization of the Coriolis terms. Three generalizations of this discretization that conserve energy even for non-constant Coriolis parameter are presented. A quasigeostrophic β-plane analysis is carried out to investigate the Rossby mode dispersion properties of these three schemes. The Rossby mode dispersion relation is found to have two branches. The primary branch modes are good approximations, in terms of both structure and frequency, to corresponding modes of the continuous governing equations, and offer some improvements over a quadrilateral C-grid scheme. The secondary branch modes have vorticity structures approximating those of smallscale modes of the continuous governing equations, suggesting that the hexagonal C-grid might have an advantage in terms of resolving extra Rossby modes; however, the frequencies of the secondary branch Rossby modes are much smaller than those of the corresponding continuous modes, so this potential advantage is not fully realized.
2012
A new approach is developed to simulate the effect of far field tsunami in a limited area model domain where the coastal and Island boundaries are curvilinear in nature and the bending is high. The model is designed in a boundary fitted curvilinear grid system. To simulate the effect of far field tsunami, it is considered that the tsunami source is located far away from the model domain. The coastal and island boundaries and the other open boundaries of the model domain are represented by some functions so as to generate the boundary fitted grids. To use the regular finite difference scheme a transformation is used so that the physical domain is transformed into a rectangular one. The transformed shallow water equations are then solved in the transformed domain. The response of the tsunami source due to 26 December 2004 Indonesian tsunami is computed along the western open boundary of the model domain. Based on the response of the tsunami source, an appropriate boundary condition is...
Non-linear Shallow Water Equations numerical integration on curvilinear boundary-conforming grids
2015
An Upwind Weighted Essentially Non-Oscillatory scheme for the solution of the Shallow Water Equations on generalized curvilinear coordinate systems is proposed. The Shallow Water Equations are expressed in a contravariant formulation in which Christoffel symbols are avoided. The equations are solved by using a high-resolution finite-volume method incorporated with an exact Riemann Solver. A procedure developed in order to correct errors related to the difficulties of numerically satisfying the metric identities on generalized boundary-conforming grids is presented; this procedure allows the numerical scheme to satisfy the freestream preservation property on highly-distorted grids. The capacity of the proposed model is verified against test cases present in literature. The results obtained are compared with analytical solutions and alternative numerical solutions. Key-Words: 2D Shallow Water Equations, Upwind WENO scheme, Contravariant formulation, Christoffel Symbols, Freestream pre...
Numerical Simulation of Propagation and Run-Up of Long Waves in U-Shaped Bays
Fluids, 2021
Wave propagation and run-up in U-shaped channel bays are studied here in the framework of the quasi-1D Saint-Venant equations. Our approach is numerical, using the momentum conserving staggered-grid (MCS) scheme, as a consistent approximation of the Saint-Venant equations. We carried out simulations regarding wave focusing and run-ups in U-shaped bays. We obtained good agreement with the existing analytical results on several aspects: the moving shoreline, wave shoaling, and run-up heights. Our findings also confirm that the run-up height is significantly higher in the parabolic bay than on a plane beach. This assessment shows the merit of the MCS scheme in describing wave focusing and run-up in U-shaped bays. Moreover, the MCS scheme is also efficient because it is based on the quasi-1D Saint-Venant equations.
Modeling of wave run-up by using staggered grid scheme implementation in 1D Boussinesq model
A new numerical method is presented to study free surface waves in coastal areas. The method is based on the phase resolving variational Boussinesq model (VBM) which is solved on a computational staggered grid domain. In this model, the nonhydrostatic pressure term has been incorporated in order to correctly described short wave dynamics. In simulating run-up phenomena, a special treatment, so-called thin layer method, is needed for solving the elliptic equation of the Boussinesq model. As a result, the proposed scheme is capable of simulating various run-up phenomena with great accuracy. Several benchmark tests were conducted, i.e., run-up experiments by Synolakis (Int. J. Numer. Methods Fluids 43(12), 1329-1354 1987) for non-breaking and breaking case, a run-up case proposed by Carrier and Greenspan (J. Fluid Mech. 4(1), 97-109 1958) and a dam-break with shock wave Aureli et al. (J. Hydraul. Res. 38(3), 197-206 2000). Moreover, the ability of the numerical scheme in simulating dispersion and nonlinearity effects were shown via simulation of the broad band waves propagation, i.e., focusing wave and irregular wave. The propagation of regular wave above a submerged trapezoidal bar was shown to confirm Beji-Batjes experiment (Coast. , 151-162 1993). Moreover, the numerical model is tested for simulating regular wave breaking on a plane beach of Ting and Kirby (Coast. , 51-80 1995), and for simulating random wave over a barred beach of Boer (1996).
WSEAS transactions on fluid mechanics, 2022
We present a new three-dimensional numerical model for the simulation of breaking waves. In the proposed model, the integral contravariant form of the Navier-Stokes equations is expressed in a curvilinear moving coordinate system and are integrated by a predictor-corrector method. In the predictor step of the method, the equations of motion are discretized by a shock-capturing scheme that is based on an original highorder scheme for the reconstruction of the point values of the conserved variables on the faces of the computational grid. On the cell faces, the updating of the point values of the conserved variables is carried out by an exact Riemann solver. The final flow velocity field is obtained by a corrector step which is based exclusively on conserved variables, without the need of calculating an intermediate field of primitive variables. The new three-dimensional model significantly reduces the kinetic energy numerical dissipation introduced by the scheme. The proposed model is validated against experimental tests of breaking waves and is applied to the three-dimensional simulation of the local vortices produced by the interaction between the wave motion and an emerged barrier.
Journal of Computational Physics, 2003
In the present study, an existing two-dimensional boundary-fitted model [J. Hydraul. Eng.-ASCE 122 (9) (1996) 512] is used to study the effect of grid non-orthogonality on the solution of shallow water equations using boundary-fitted grids. The linearized two-dimensional shallow water equations are expressed in terms of the grid angle and aspect ratio. The truncation errors of the finite difference approximations used in the solution of the governing equations are shown to be dependent on the grid angle and the aspect ratio. The coefficient of the truncation error was shown to increase, with the decrease in the grid angle. The RMS errors in model predicted surface elevations and velocities for the case of seiching in a rectangular basin are found to increase gradually, as the grid resolution decreases from 174 to 80 gridpoints per wavelength or as the grid angle decreases from 90°to 50°and increases rather sharply for a grid angle of 30°a t grid resolutions less than 80 gridpoints per wavelength. The model predicted surface elevations for the case of tidal forcing in a rectangular basin are found to be insensitive to the grid angle at grid resolutions higher than 600 gridpoints per wavelength. The RMS error in the model predicted velocities is found to increase gradually as the grid angle decreases from 90°to 30°or as the grid resolution decreases from 1400 gridpoints per wavelength to 400 gridpoints per wavelength and increases sharply as the grid resolution decreases from 400 to 150 gridpoints per wavelength. Twodimensional depth averaged hydrodynamic modeling of tidal circulation in Narragansett Bay, using three different boundary-fitted grids showed that the model predicted surface elevations are insensitive to the grid angle at grid resolutions as low as 200 gridpoints per wavelength. However, the model predicted velocities were found to increase as the grid resolution decreases from 600 to 200 gridpoints per wavelength. We conclude from this study that grid angle and grid resolution affects the accuracy of the model predicted currents and the numerical dispersion increases with the decrease in grid angle or grid resolution and these are in agreement with that reached by Sankaranarayanan and Spaulding [Dispersion and Stability Analyses of Shallow Water Equations in Boundary-fitted Coordinates,
2017
A numerical technique in a boundary-fitted curvilinear grid model is developed to simulate the extent of inland inundation along the coastal belts of Peninsular Malaysia and Southern Thailand due to 2004 Indian ocean tsunami. Tsunami propagation and run-up are also studied in this paper. The vertically integrated shallow water equations are solved by using the method of lines (MOL). For this purpose the boundary-fitted grids are generated along the coastal and island boundaries and the other open boundaries of the model domain. A transformation is used to the governing equations so that the transformed physical domain is converted into a rectangular one. The MOL technique is applied to the transformed shallow water equations and the boundary conditions so that the equations are converted into ordinary differential equations initial value problem. Finally the 4<sup>th</sup> order Runge-Kutta method is used to solve these ordinary differential equations. The moving boundar...