A Three-Dimensional Environmental Fluid Dynamics Computer Code : Theoretical and computational aspects

This report describes and documents the theoretical and computational aspects of a three-dimensional computer code for environmental fluid flows. The code solves the three-dimensional primitive variable v1ertically hydrostatic equations of motion for turbulent flow in a coordinate system which is curvilinear and orthogonal in the horizontal plane and stretched to follow bottom topography and free surface displacement in the vertical direction which is aligned with the gravitational vector. A second moment turbulence closure scheme relates turbulent viscosity and diffusivity to the turbulence intensity and a turbulence length scale. Transport equations for the turbulence intensity and length scale as well as transport equations for and a dye tracer are also to pressure, salinity, concentration. salinity, temperature, suspended sediment solved. An equation of state relates density temperature and suspended sediment The computational scheme utilizes an external-internal mode splitting to solve the horizontal momentum equations and the continuity equation on a staggered grid. The external mode, associated with barotropic long wave motion, is solved using a semiimplicit three time level scheme with a periodic two time level correction. A multi-c:olor successive over relaxation scheme is used to solve the resulting system of equations for the free surface displacement. The internal mode, associated with vertical shear of the horizontal velocity components is solved using a fractional step scheme combining an implicit step for the vertical shear terms, with an explicit step for all other terms. The transport equations for the turbulence intensity, turbulence length scale, salinity, temperature, suspended sediment and dye tracer are also solved using a fractional step scheme with implicit vertical diffusion and explicit advection and horizontal diffusion. A number of alternate advection schemes are implemented in th1e code.

[1]  J. E. Glynn,et al.  Numerical Recipes: The Art of Scientific Computing , 1989 .

[2]  T. J. Simons Verification of Numerical Models of Lake Ontario: Part I. Circulation in Spring and Early Summer , 1974 .

[3]  G. Mellor An Equation of State for Numerical Models of Oceans and Estuaries , 1991 .

[4]  John M. Hamrick,et al.  Long-term dispersion in unsteady skewed free surface flow , 1986 .

[5]  G. Mellor,et al.  Development of a turbulence closure model for geophysical fluid problems , 1982 .

[6]  E. D Chikhliwala,et al.  Application of orthogonal mapping to some two-dimensional domains , 1985 .

[7]  Akio Arakawa,et al.  Computational Design of the Basic Dynamical Processes of the UCLA General Circulation Model , 1977 .

[8]  W. Grabowski,et al.  The multidimensional positive definite advection transport algorithm: nonoscillatory option , 1990 .

[9]  M. Foreman,et al.  An analysis of two-step time discretizations in the solution of the linearized shallow water equations , 1983 .

[10]  Jan J. Leendertse,et al.  Multidimensional Numerical Modeling of Estuaries and Coastal Seas , 1978 .

[11]  Pz,et al.  Mesoscale meteorological modeling , 1986 .

[12]  Ralph T. Cheng,et al.  Tidal and tidally averaged circulation characteristics of Suisun Bay, California , 1987 .

[13]  C. Fletcher Computational techniques for fluid dynamics , 1992 .

[14]  P. Smolarkiewicz A Fully Multidimensional Positive Definite Advection Transport Algorithm with Small Implicit Diffusion , 1984 .

[15]  M. Cullen,et al.  Numerical Prediction and Dynamic Meteorology, 2nd Edn. By G. J. HALTINER and R. T. WILLIAMS. Wiley, 1980. 477 pp. £26.90. , 1984, Journal of Fluid Mechanics.

[16]  Richard B. Rood,et al.  Numerical advection algorithms and their role in atmospheric transport and chemistry models , 1987 .

[17]  P. Smolarkiewicz,et al.  The multidimensional positive definite advection transport algorithm: further development and applications , 1986 .

[18]  Alan F. Blumberg,et al.  Modeling vertical and horizontal diffusivities with the sigma coordinate system , 1985 .

[19]  H. Fischer Mixing in Inland and Coastal Waters , 1979 .

[20]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .