3D Variable-Boundary Numerical Modeling of Tidal Current in the Xiangshangang Bay
吴晓燕,管卫兵
DOI: https://doi.org/10.3969/j.issn.1001-909x.2009.02.004
2009-01-01
Journal of Marine Sciences
Abstract:The Xiangshangang Bay is a typical narrow semi-enclosed bay in the middle of the Zhejiang Province,located at the south edge of the Changjiang River Delta.Along its axis from southwest to northeast,the Xiangshangang Bay extends about 60 km with a water area of 563 km2,in which 171 km2 is tidal flat(above the base level of theoretical depth).The tidal range in the bay is large.It increases gradually from the entrance inwards with a mean value of above 3 m.At the innermost part it reaches a maximum of 6 m with a mean range of 4 m or so.The water area of the mudflats changes markedly with tide.When the water level is high,the whole bay is merged by water,and when the water level falls,the shallow areas are exposed.Therefore,the influence of drying and flooding of the mudflats should be considered when a numerical model is applied for such a region with wide mudflats and large tidal range.Various drying and of flooding algorithms have been developed very well for 2D models.However,if those algorithms used in 2D models are directly introduced into 3D models,the numerical oscillations or instability could occur,and thus it is necessary to find an easy and feasible mobile boundary technique for 3D models.The original POM(Princeton Ocean Model) with a fixed boundary cannot deal with the overbank flow.Base on the work of Burchard et al in 2004,the fixed boundary in original POM was replaced with a naturally variable boundary following a strategy to let the physics do the best to meet the requirement of flow reduction.Some simplifications of the model physics will be applied for very shallow water in order to describe the physical process in a proper way.In this method,a dimensionless damping factor,which varies linearly between unity(deep water) and zero(very shallow water),was introduced to the horizontal advection,horizontal diffusion,Coriolis force and baroclinicity terms of the momentum equations.When the whole water depth is larger than the set critical water depth,usual momentum equations are resulted,and when the whole water depth trends approaching the minimum value,the momentum equations are simplified to be the balance among the tendency,friction and external pressure gradient.The friction dominates in such a way that for water depths approaching zero,the current velocities converge to zero as well,and as a result,the unphysical negative water depths can be avoided.In addition to the simplification of the momentum equations at very shallow water depths,a correction to the barotropic pressure gradient should be applied where the sea level cell is below the sea bed at the adjacent grid cell,so as to avoid unreasonable pressure gradient.Furthermore,the expression of the frictional coefficient at the bottom was modified so that the model becomes more stable to deal with drying and flooding of mudflats.The improved POM was applied to simulate 3D tidal current field of the Xiangshangang Bay.The computed results matched very well the observed data from several survey stations including 2 tidal level stations,1 mooring site and 2 anchor stations.The model has reasonably reproduced the 3D structure and variation of the tidal current in the whole Xiangshangang Bay.Horizontal and vertical distributions of the tidal velocity in the innermost part of the bay at the beginning of the ebb tide,the peak ebb tidal currents,the beginning of the flood tide and the peak flood tidal currents show that the model is able to present the correct processes of drying and flooding of mudflats.It is hence indicated that the modification of the origin POM is very reasonable and feasible,and will enable POM to be used more extensively.