Implementation of baroclinic terms in a three-dimensional hydrodynamic model and its application to Anzali Port, Iran

Baroclinic terms have been implemented in a three-dimensional fully hydrodynamic model developed by Badiei et al. [2008. A three-dimensional non-hydrostatic boundary fitted model for free surface flows. International Journal for Numerical Methods in Fluids, 56(6), 607-627] modifying its momentum equations to account for density gradients and utilizing the scalar (salinity, temperature, etc.) conservation equation (SCE) and a state equation for the calculation of density. In the solution of advection-diffusion terms of the governing Navier-Stokes equations (NSE) and SCE, a symmetric splitting method was applied to ensure the long-term stability of simulations. Correction terms proposed by Ruddic et al. (1995) were applied to SCE to ensure the conservation of the scalar quantity. In the presence of baroclinic terms, the zero gradient pressure in the vertical direction in the vicinity of surface and bottom boundaries assumed by Badiei et al. [2008. A three-dimensional non-hydrostatic boundary fitted model for free surface flows. International Journal for Numerical Methods in Fluids, 56(6), 607-627] created spurious currents. This problem was solved by assuming a hydrostatic pressure variation at those boundaries. The ability of extended model was validated by comparing its results with an experimental test case. The simulation of hydrodynamic and salt intrusion at Anzali Port located at the southern coasts of Caspian Sea in Iran was carried out by the model with both barotropic and baroclinic modes. The simulated results with baroclinic mode show a better agreement with measured data as compared to the results of barotropic mode that clearly demonstrate the significance of baroclinic terms in the simulation of cyclic intrusion of salt wedge into the Port Basin.