Two-Dimensional Boundary-Fitted Circulation Model in Spherical Coordinates
Publication: Journal of Hydraulic Engineering
Volume 122, Issue 9
Abstract
A two-dimensional, vertically averaged, unsteady circulation model, using a nonorthogonal boundary-fitted technique, was developed in spherical coordinates for predicting sea level and currents in estuarine and shelf waters. Both the dependent and independent variables are transformed into a curvilinear coordinate system. The governing equations are solved by a semiimplicit method in which the elevations are solved implicitly and the vertically averaged velocities are solved explicitly. The model employs a space-staggered grid system and a three-level time discretization. Truncation errors are second order both in space and time. The model was tested against analytic solutions for a standing wave in a closed basin, tidal circulation in a simple rectangular channel with an irregular grid system and various degrees of rotation, and tidal flow in an annular section channel with quadratic bottom topography. The model was also tested against steady-state wind-induced setup in a closed irregular basin with constant depth represented by an irregular grid system. Comparison of the model predictions with the corresponding analytical solutions were very good. The model was applied to simulate tidal circulation in the Providence River. The agreement with available observations is very good. The model predicts that the tide exhibits a cooscillating wave pattern with tidal currents leading tidal elevation by 2.8 to 3.8 h depending on location for the M2 tidal constituent. The M4 and M6 tidal components are significantly amplified because their frequency is close to the resonant frequency of the bay.
Get full access to this article
View all available purchase options and get full access to this article.
References
1.
Arakawa, A., and Lamb, V. R.(1977). “Computational design of the basic dynamical process of the UCLA General Circulation Model.”Method in Computational Physics, 17, 173–265.
2.
Barber, R. W.(1992). “Solving the shallow water equations using a non-orthogonal curvilinear coordinate system.”Proc. 2nd Int. Conf. on Hydr. and Envir. Modelling of Coast., Estuarine and River Waters, U.K., 1, 469–480.
3.
Blumberg, A. F., and Herring, H. J. (1987). “Circulation modelling using orthogonal curvilinear coordinates.”Three Dimensional Model of Marine and Estuarine Dynamics, J. C. J. Nihoul and B. M. Jamart, eds., Univ. of Liège, Liège Belgium, 55–75.
4.
Chandler-Wilde, S. N., and Lin, B.(1992). “A finite-difference method for the shallow water equations with conformal boundary-fitted mesh generation.”Proc. 2nd Int. Conf. on Hydr. and Envir. Modelling of Coast., Estuarine and River Waters, U.K., 1, 507–518.
5.
Gordon, R. B., and Spaulding, M. L.(1987). “Numerical simulation of the tidal and wind driven circulation in Narragansett Bay.”Estuarine, Coast., and Shelf Sci., 24, 611–636.
6.
Hauser, J., Paap, H. G., Eppel, D., and Mueller, A.(1985). “Solution of shallow water equations for complex flow domains via boundary-fitted coordinates.”Int. J. for Numer. Methods in Fluids, 5, 727–744.
7.
Johnson, B. H. (1980). “VAHM-A vertically averaged hydrodynamic model using boundary-fitted coordinates.”MP HL-80-3, U.S. Army Engrs. Wtrwy. Experiment St., Vicksburg, Miss.
8.
Kinmark, I. (1985). The shallow water wave equations: formulation, analysis and application . Springer Verlag, New York, N.Y.
9.
Leendertse, J. J. (1967). “Aspects of computational model for long period water-wave propagation.”Rep. RM-5294-PR, Rand Corp., Santa Monica, Calif.
10.
Lynch, D. R., and Gray, W. G.(1978). “Analytic Solutions for Computer Flow Model testing.”J. Hydr. Div., ASCE, 104(6), 1409–1428.
11.
Madala, R. V., and Piacsek, S. A.(1977). “A semi implicit numerical model for baroclinic oceans.”J. Computation Physics, 23, 167–178.
12.
Mendelsohn, D., and Swanson, C. (1990). “A diagnostic simulation of the density induced flow in the Providence River.”Report prepared for the Narragansett Bay Project, Dept. of Envir. Mgmt., Providence, R.I.
13.
Muin, M. (1993). “A three dimensional boundary-fitted circulation model in spherical coordinates, PhD dissertation, Dept. of Oc. Engrg., Univ. of Rhode Island, Narragansett Bay Campus, R.I.
14.
Nielson, P., and Skovgaard, O.(1990). “The effect of using non-orthogonal boundary fitted grids for solving the shallow water equations.”Int. J. for Numer. Methods in Fluids, 11, 177–188.
15.
Oey, L. Y., Mellor, G., and Hires, R. I.(1985). “A three-dimensional simulation of the Hudson-Raritan Estuary. part II: comparison with observation.”J. Phys. Oceanography, 15(12), 1693–1709.
16.
Pearson, R. V., and Barber, R. W. (1993). “Mathematical simulation of the Humber Estuary using a depth-averaged boundary-fitted tidal model.”Numer. Methods in Laminar and Turbulent Flow, VIII, Part 2, 1244–1255.
17.
Sheng, Y. P. (1986). “Numerical modeling of coastal and estuarine processes using boundary-fitted grids.” in 3rd Int. Symp. on River Sedimentation, 1426–1442.
18.
Sheng, Y. P. (1990). “Evolution of a three-dimensional curvilinear-grid model for estuaries, lakes and coastal waters: CH3D. Estuarine and Coast. Modeling, Malcolm L. Spaulding, ed., ASCE, New York, N.Y., 40–49.
19.
Spaulding, M. L., Bedford, K., Blumberg, A., Cheng, R., and Swanson, C., eds., (1994). Estuarine and Coast. Modelling III, Proc. 3rd Int. Conf., ASCE Wtrwy., Port, Coast., and Oc. Div., New York, N.Y.
20.
Spaulding, M. L., Bedford, K., Blumberg, A., Cheng, R., and Swanson, C., eds., (1992). Estuarine and Coast. Modeling II, Proc., 2nd Int. Conf., ASCE Wtrwy., Port Coast. and Oc. Div., New York, N.Y.
21.
Spaulding, M. L., Swanson, J. C., and Turner, C. (1990). “The new tide and tidal current of Narragansett Bay.”Marine Tech. Rep., Univ. of Rhode Island. Kingston, R.I.
22.
Spaulding, M. L.(1984). “A vertically averaged circulation model using boundary-fitted coordinates.”J. Physical Oceanography, 14, 973–982.
23.
Swanson, J. C. (1986). “A three dimensional numerical model system of coastal circulation and water quality,” PhD dissertation, Dept. of Oc. Engrg., Univ. of Rhode Island, Kingston, R.I.
24.
Thompson, J. F., Thames, F. C., and Mastin, C. W. (1977). “Boundary fitted curvilinear coordinate system for solution of partial differential equations on fields containing any number of arbitrary two dimensional bodies.”NASA CR-2729, Nat. Aeronautics and Space Admin., Washington, D.C.
25.
Thompson, J. F., Warsi, Z. U. A., and Martin, C. W. (1985). Numerical grid generation, foundation and applications . North-Holland, Amsterdam.
26.
Turner, A. C. (1984). “Tidal and subtidal circulation in the Providence River,” MS thesis, Dept. of Oc. Engrg., Univ. of Rhode Island, Kingston, R.I.
27.
Turner, A. C., Puckett, S., Isaji, T., and Feng, S. S. (1989). “Providence River wet weather current meter study data report,”Rep. prepared for Narragansett Bay Project, Dept. of Envir. Mgmt., Providence, R.I.
28.
Willemse, J. B. T. M., Stelling, G. S., and Verboom, G. K. (1985). “Solving the shallow water equations with an orthogonal coordinate transformation.”Int. Symp. on Computational Fluid Dyn., Delft Hydr. Communication No. 356, Jan. 86, Delft Hydraulics Laboratory, Delft, The Netherlands.
Information & Authors
Information
Published In
Journal of Hydraulic Engineering
Volume 122 • Issue 9 • September 1996
Pages: 512 - 521
Copyright
Copyright © 1996 American Society of Civil Engineers.
History
Published online: Sep 1, 1996
Published in print: Sep 1996
Authors
Metrics & Citations
Metrics
Citations
Download citation
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.