Entrainment and mixing in lock-exchange gravity currents are investigated by large eddy simulations. Nine cases are analysed, varying the initial excess density driving the motion and the aspect ratio r of the initial water depth to the lock length. Laboratory experiments are also performed and a fair agreement between numerical simulations and measurements is found. Mixing between the gravity current and the ambient fluid, in both the slumping and self-similar phases, is studied for a range of entrainment parameters, gravity current fractional area and using an energy budget method. The entrainment is found to increase as r decreases. The occurrence of irreversible mixing is detected during the entire development of the flow, i.e. both in the slumping and self-similar phases. A higher amount of mixing is observed as r decreases and the initial excess density increases

Lock-exchange gravity currents propagating up a slope are investigated by large eddy simulations, focusing on the entrainment and mixing processes occurring between the dense current and the ambient fluid. Relevant parameters, such as the aspect ratio of the initial volume of dense fluid in the lock R, the angle between the bottom boundary and the horizontal direction θ and the depth aspect ratio φ, are varied. The numerical results are compared with laboratory experiments and a good agreement is found. Entrainment and mixing in a lock-release gravity current are studied using different entrainment parameters and an energy budget method. The entrainment is found to depend on both Froude, Fr, and Reynolds, Re, numbers. In addition, the dependence of both entrainment and mixing on the parameters varied is discussed. The entrainment decreases with increasing steepness of the bottom and R. Irreversible mixing is not affected by the varied parameters during the slumping phase, while during the successive phases of motion, it is found to decrease with the increase of θ and R. Low entrainment and mixing occur for φ < 1

No abstract

Dense oceanic overflows descend over the rough topography of the continental slope entraining and mixing with surrounding waters. The associated dilution dictates the fate of these currents and thus is of fundamental importance to the formation of deep water masses. The entrainment in a dense current flowing down a sloping bottom in a rotating homogeneous fluid is investigated using laboratory experiments, focusing on the influence of the bottom roughness on the flow dynamics. The roughness is idealized by an array of vertical rigid cylinders and both their spacing and height are varied as well as the inclination of the sloping bottom. The presence of the roughness is generally observed to decelerate the dense current, with a consequent reduction of the Froude number, when compared to the smooth bottom configuration. However, the dilution of the dense current due to mixing with the ambient fluid is enhanced by the roughness elements, especially for low Froude numbers. When the entrainment due to shear instability at the interface between the dense current and the ambient fluid is low, the additional turbulence and mixing arising at the bottom of the dense current due to the roughness elements strongly affects the dilution of the current. Finally, a strong dependence of the entrainment parameter on the Reynolds number is observed.

Triggering and evolution of internal solitary waves (ISWs) generated by intrusive gravity currents (IGCs) propagating into a stratified ambient fluid is analyzed by laboratory experiments. After the release of a fluid of uniform density, intermediate with respect to the upper (lower-density) and lower (higher-density) layers in the channel, the IGC develops and flows downstream, intruding into the pycnocline. Near the IGC leading front, the compression of the upper layer generates ISWs: they gradually separate from the current that propagates slower. Shoaling downstream over a uniform sloping boundary, solitons break and partially reflect. We investigate the dynamics of the interaction between the reflected ISWs and the incoming IGC. During the engage, an increase in the ISW celerity occurs, leading the celerity of the reflected waves to be even larger than the incident wave. Our analysis shows how both ISWs and IGCs can significantly change their features as they experience a change of the density structure in the water column. This is expected to occur, for example, in stratified small-scale basins, where river plumes intrude the seasonal thermocline. The radial ISWs, originated by IGCs, can then be reflected by the adjacent bottom bathymetry, spreading against the intrusive current from which they are generated.

The dynamics of lock-release Intrusive Gravity Currents (IGCs) generating Internal Solitary Waves (ISWs) are investigated by threedimensional large eddy simulations. We set the numerical, laboratory-scale domain in order to release a uniform fluid in multi-layer, stratified ambient, exciting pycnocline displacements. By adopting different initial settings, we analyzed the influence of the ambient stratification on both IGCs and ISWs features. We present the main flow dynamics and the time evolution of IGC and ISW front and trough positions, respectively. During the simulations, the ISW is allowed to reach the vertical wall at the end of the domain, and it undergoes reflection. We then analyzed the interaction between the IGC and the reflected ISW: the wave is observed to accelerate as it is pushed upwards by the intrusion, which, in turns, flows below the ISW, decelerating. By analyzing instantaneous velocity fields and flow rates, we found that during this interaction, the ISW increases its celerity in response of the reduced area available for its propagation, partially occupied by the intrusion, and because the velocity field in the IGC interface surroundings acts to facilitate the ISW passage.

No abstract

Abstract. The aim of this study was to determine the dispersion of passive pollutants associated with the Tiber discharge into the Tyrrhenian Sea using numerical marine dispersion models and satellite data. Numerical results obtained in the simulation of realistic discharge episodes were compared with the corresponding evolution of the spatial distributions of MODIS diffuse light attenuation coefficient at 490 nm (K490), and the results were discussed with reference to the local climate and the seasonal sub-regional circulation regime. The numerical model used for the simulation of the sub-tidal circulation was a Mediterranean sub-regional scale implementation of the Princeton Ocean Model (POM), nested in the large-scale Mediterranean Forecasting System. The nesting method enabled the model to be applied to almost every area in the Mediterranean Sea and also to be used in seasons for which imposing climatological boundary conditions would have been questionable. Dynamical effects on coastal circulation and on water density due to the Tiber discharge were additionally accounted for in the oceanographic model by implementing the river estuary as a point source of a buoyant jet. A Lagrangian particle dispersion model fed with the POM current fields was then run in order to reproduce the effect of the turbulent transport of passive tracers mixed in the plume with the coastal flow. Two significant episodes of river discharge in both winter and summer conditions were discussed in this paper. It was found that the winter regime was characterized by the presence of a strong coastal jet flowing with the ambient current. In summer the prevailing wind regime induced coastal downwelling conditions, which tended to confine the riverine waters close to the shore. In such conditions sudden wind reversals due to local weather perturbations, causing moderate local upwelling, proved to be the only effective way to disperse the tracers offshore, moving the plume from the coast and detaching large pools of freshwater.