The finite volume method (FVM) has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Ro...The finite volume method (FVM) has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Roe scheme is applied in the calculation of the intercell numerical flux, and the MUSCL method is introduced to improve its accuracy. The time integral is a two-step scheme of forecast and revision. For the verification of the present method, the Stoker's problem is calculated and the result is compared with the mathematically analytic solutions. The comparison indicates that the method is feasible. A sea area of a port is used as an example to test the method established here. The result shows that the present computational method is satisfactory, and it could be applied to the engineering fields.展开更多
A depth-averaged 2-D numerical model for unsteady tidal flow in estuaries is established by use of the finite volume WENO scheme which maintains both uniform high order accuracy and an essentially non-oscillatory shoc...A depth-averaged 2-D numerical model for unsteady tidal flow in estuaries is established by use of the finite volume WENO scheme which maintains both uniform high order accuracy and an essentially non-oscillatory shock transition on unstructured triangular grid. The third order TVD Range-Kutta method is used for time discretization. The model has been firstly tested against four cases: 1) tidal forcing, 2) seiche oscillation, 3) wind setup in a closed bay, and 4) onedimensional dam-break water flow. The results obtained in the present study compare well with those obtained from the corresponding analytic solutions idealized for the above four cases. The model is then applied to the simulation of tidal circulation in the Yangpu Bay, and detailed model calibration and verification have been conducted with measured tidal current in the spring tide, middle tide, and neap tide. The overall performance of the model is in qualitative agreement with the data observed in 2005, and it can be used to calculate the flow in estuaries and coastal waters.展开更多
The characteristics of three-dimensional (3-D) tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, ma...The characteristics of three-dimensional (3-D) tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, σ 3-D numerical model oil an unstructured triangular grid has been degeloped. The σ coordinate transforination, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simulations in the Oujing Estuary with submerged dyke and moveable boundary problems. The model is evaluated with the in situ data, and the results show that the calculated water elevations at 19 stations and currents at 19 profiler stations are in good agreement with measured data both in magnitude and phase. This numerical model is applied to the 3-D tidal circulation simulations of experiments in stopping flow transport through the South Branch of the Oujiang Estuary, and the feasibility to cutoff the flow in the South Branch of the Oujiang Estuary is demonstrated by numerical simulation experiments. The developed numerical model simulated the 3-D tidal current circulations in complicated coastal and estuarine waters very well.展开更多
Constructing Bernstein-Bezier triangular interpolating curve surface interpolating a series of arbitrary disordered data points is of considerable importance for the design and modeling of surfaces with a variety of c...Constructing Bernstein-Bezier triangular interpolating curve surface interpolating a series of arbitrary disordered data points is of considerable importance for the design and modeling of surfaces with a variety of continuity information. In this article. a kind of simple and reliable algorithm that can process complex field triangular grid generating is presented, and a group of formulae for determining triangular curved surface with wholly C1 continuity are given. It can process arbitrary non-convex boundary and can be used to construct surfaces inner holes.展开更多
基金This paper was supported bythe Natural Science Foundation of Shandong Province (Grant No.y2004f13)
文摘The finite volume method (FVM) has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Roe scheme is applied in the calculation of the intercell numerical flux, and the MUSCL method is introduced to improve its accuracy. The time integral is a two-step scheme of forecast and revision. For the verification of the present method, the Stoker's problem is calculated and the result is compared with the mathematically analytic solutions. The comparison indicates that the method is feasible. A sea area of a port is used as an example to test the method established here. The result shows that the present computational method is satisfactory, and it could be applied to the engineering fields.
基金This work was supported by Open Research Fund Programof State Key Laboratory of Water Resources and Hydropow-er Engineering Science ( Grant No. 2005C011)National Natural Science Foundation of China ( Grant No.50479038)
文摘A depth-averaged 2-D numerical model for unsteady tidal flow in estuaries is established by use of the finite volume WENO scheme which maintains both uniform high order accuracy and an essentially non-oscillatory shock transition on unstructured triangular grid. The third order TVD Range-Kutta method is used for time discretization. The model has been firstly tested against four cases: 1) tidal forcing, 2) seiche oscillation, 3) wind setup in a closed bay, and 4) onedimensional dam-break water flow. The results obtained in the present study compare well with those obtained from the corresponding analytic solutions idealized for the above four cases. The model is then applied to the simulation of tidal circulation in the Yangpu Bay, and detailed model calibration and verification have been conducted with measured tidal current in the spring tide, middle tide, and neap tide. The overall performance of the model is in qualitative agreement with the data observed in 2005, and it can be used to calculate the flow in estuaries and coastal waters.
基金The Natural Science Foundation of Tianjin, China under contract No.08JCZDZT00200
文摘The characteristics of three-dimensional (3-D) tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, σ 3-D numerical model oil an unstructured triangular grid has been degeloped. The σ coordinate transforination, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simulations in the Oujing Estuary with submerged dyke and moveable boundary problems. The model is evaluated with the in situ data, and the results show that the calculated water elevations at 19 stations and currents at 19 profiler stations are in good agreement with measured data both in magnitude and phase. This numerical model is applied to the 3-D tidal circulation simulations of experiments in stopping flow transport through the South Branch of the Oujiang Estuary, and the feasibility to cutoff the flow in the South Branch of the Oujiang Estuary is demonstrated by numerical simulation experiments. The developed numerical model simulated the 3-D tidal current circulations in complicated coastal and estuarine waters very well.
文摘Constructing Bernstein-Bezier triangular interpolating curve surface interpolating a series of arbitrary disordered data points is of considerable importance for the design and modeling of surfaces with a variety of continuity information. In this article. a kind of simple and reliable algorithm that can process complex field triangular grid generating is presented, and a group of formulae for determining triangular curved surface with wholly C1 continuity are given. It can process arbitrary non-convex boundary and can be used to construct surfaces inner holes.