Using a dissolved conservative material as a tracer in the water, a vertically integrated 2-D advection-diffusion water exchange numerical model was developed for the water exchange study in Xiangshangang Bay. The model consists of a hydrodynamic model and an advection-diffusion tracer model. The governing equation of the tracer model is formed by integrating the 3-dimensional tracer mass conservation equation. The effects of the gravity residual circulation and vertical tidal current shear on the water exchange are parameterized with the dispersion coefficients as follows:
where, ui is the velocities; c is the tracer concentration; kir and kis are the horizontal dispersion coefficients caused by the gravity residual circulation and the vertical shear effect of the tidal current oscillation in the tidal period respectively; kit is the vertically integrated turbulent coefficient. In the numerical model, the tracer flux normal to the solid boundary is set to zero. The tracer concentration is also set to be zero on the open boundary when there is inflow. Providing the initial tracer concentration is c(t0) in the bay, the water exchange rate R(x, y, t) can be calculated in terms of the tracer concentration c(x, y, t) as follow:
The TRIM model is accepted as the hydrodynamic model and its numerical scheme is applied in the numerical tracer model. Generally, the stability and the conservation of tracer mass of the numerical model can satisfy the requirement of the water exchange study in the estuary and the embayment. Because the water exchange model is 2-dimensional, a lot of time is saved in calculation. As the main vertical mixing mechanisms are included, the model simulates well the water exchange process. |