Effects of the bottom boundary condition in numerical investigations of dense water cascading on a slope

The flow of dense water along continental slopes is considered. There is a large literature on the topic based on observations and laboratory experiments. In addition, there are many analytical and numerical studies of dense water flows. In particular, there is a sequence of numerical investigations using the dynamics of overflow mixing and entrainment (DOME) setup. In these papers, the sensitivity of the solutions to numerical parameters such as grid size and numerical viscosity coefficients and to the choices of methods and models is investigated. In earlier DOME studies, three different bottom boundary conditions and a range of vertical grid sizes are applied. In other parts of the literature on numerical studies of oceanic gravity currents, there are statements that appear to contradict choices made on bottom boundary conditions in some of the DOME papers. In the present study, we therefore address the effects of the bottom boundary condition and vertical resolution in numerical investigations of dense water cascading on a slope. The main finding of the present paper is that it is feasible to capture the bottom Ekman layer dynamics adequately and cost efficiently by using a terrain-following model system using a quadratic drag law with a drag coefficient computed to give near-bottom velocity profiles in agreement with the logarithmic law of the wall. Many studies of dense water flows are performed with a quadratic bottom drag law and a constant drag coefficient. It is shown that when using this bottom boundary condition, Ekman drainage will not be adequately represented. In other studies of gravity flow, a no-slip bottom boundary condition is applied. With no-slip and a very fine resolution near the seabed, the solutions are essentially equal to the solutions obtained with a quadratic drag law and a drag coefficient computed to produce velocity profiles matching the logarithmic law of the wall. However, with coarser resolution near the seabed, there may be a substantial artificial blocking effect when using no-slip.     

The role of eddies on pathways,transports, and entrainment in dense water flows along a slope

Ocean Dynamics - The flow of dense water along slopes has been investigated in several numerical investigations based on the Dynamics of Overflow Mixing and Entrainment (DOME) setup. In the present...     

Gravity currents down canyons: effects of rotation

The flow of dense water in a V-shaped laboratory-scale canyon is investigated by using a non-hydrostatic numerical ocean model with focus on the effects of rotation. By using a high-resolution model, a more detailed analysis of plumes investigated in the laboratory (Deep-Sea Res I 55:1021–1034 2008) for laminar flow is facilitated. The inflow rates are also increased to investigate plume structure for higher Reynolds numbers. With rotation, the plumes will lean to the side of the canyon, and there will be cross-canyon geostrophic currents and Ekman transports. In the present study, it is found that the cross-canyon velocities are approximately 5 % of the down-canyon velocities over the main body of the plume for the rotational case. With rotation, the flow of dense water through the body of the plume and into the plume head is reduced. The plume head becomes less developed, and the speed of advance of the head is reduced. Fluid parcels near the top of the plume will to a larger extent be left behind the faster flowing dense core of the plume in a rotating system. Near the top of the plume, the cross-canyon velocities change direction. Inside the plume, the cross-flow is up the side of the canyon, and above the interface to the ambient there is a compensating cross-flow down the side of the canyon. This means that parcels of fluid around the interface become separated. Parcels of fluid around the interface with small down-canyon velocity components and relative large cross-canyon components will follow a long helix-like path down the canyon. It is found that the entrainment coefficients often are larger in the rotational experiments than in corresponding experiments without rotation. The effects of rotation and higher inflow rates on the areal patterns of entrainment velocities are demonstrated. In particular, there are bands of higher entrainment velocities along the lateral edges of the plumes in the rotational cases.