The analysis of the importance of the forces that act over an ensemble of particles in a turbulent field has been carried out by using direct numerical simulation for a wide range of density ratios (2.65<ρ<2650). It has been observed that, compared to the Stokes drag, the added mass is always negligible, the pressure drag is relevant for density ratios O(1), and the Basset force is appreciable for the whole range investigated. However, the effect of these forces on the particle dispersion is about 1% for ρ∼1 as well as for large density ratios.

Wall roughness produces a downward shift of the mean streamwise velocity profile in the log region, known as the roughness function. The dependence of the roughness function on the height and arrangement of roughness elements has been confirmed in several studies where regular rough walls were analysed; less attention has been paid to non-regular rough walls. Here, a numerical analysis of turbulent flows over irregularly shaped rough walls is performed, clearly identifying the importance of a parameter, called the effective slope (ES) of the wall corrugations, in characterizing the geometry of non-smooth irregular walls. The effective slope proves to be one of the fundamental geometric parameters for scaling the roughness function. Specifically, for a moderate range of roughness heights, both in the transitionally and in the fully rough regime, ES appears to scale the roughness function for a wide range of irregular rough geometric configurations. The effective slope determines the relative importance of friction drag and pressure drag. For ES ~ 0.15 we find that the friction contribution to the total wall stress is nearly in balance with the pressure-drag contribution. This value separates the region where the roughness function ΔU⁺ = f(ES) is linear from that where a smooth nonlinear behaviour is observed. In the cases investigated, value ES ~ 0.15 also separates the transitionally rough regime from the fully rough regime.

Boundary-forced stratified turbulence is studied in the prototypical case of turbulent channel flow subject to stable stratification. The large-eddy simulation approach is used with a mixed subgrid model that involves a dynamic eddy viscosity component and a scale-similarity component. After an initial transient, the flow reaches a new balanced state corresponding to active wall-bounded turbulence with reduced vertical transport which, for the cases in our study with moderate-to-large levels of stratification, coexists with internal wave activity in the core of the channel. A systematic reduction of turbulence levels, density fluctuations and associated vertical transport with increasing stratification is observed. Countergradient buoyancy flux is observed in the outer region for sufficiently high stratification.Mixing of the density field in stratified channel flow results from turbulent events generated near the boundaries that couple with the outer, more stable flow. The vertical density structure is thus of interest for analogous boundary-forced mixing situations in geophysical flows. It is found that, with increasing stratification, the mean density profile becomes sharper in the central region between the two turbulent layers at the upper and lower walls, similar to observations in field measurements as well as laboratory experiments with analogous density-mixing situations.Channel flow is strongly inhomogeneous with alternative choices for the Richardson number. In spite of these complications, the gradient Richardson number, Ri_g, appears to be the important local determinant of buoyancy effects. All simulated cases show that correlation coefficients associated with vertical transport collapse from their nominal unstratified values over a narrow range, 0.15 < Ri_g < 0.25. The vertical turbulent Froude number, Fr_w, has an O(1) value across most of the channel. It is remarkable that stratified channel flow, with such a large variation of overall density difference (factor of 26) between cases, shows a relatively universal behaviour of correlation coefficients and vertical Froude number when plotted as a function of Ri_g. The visualizations show wavy motion in the core region where the gradient Richardson number, Ri_g, is large and low-speed streaks in the near-wall region, typical of unstratified channel flow, where Ri_g is small. It appears from the visualizations that, with increasing stratification, the region with wavy motion progressively encroaches into the zone with active turbulence; the location of Ri_g ≃ 0.2 roughly corresponds to the boundary between the two zones.

In this study, the effects of small-scale velocity fluctuations on the motion of tracer particles is investigated by releasing particles in a turbulent channel flow at Reτ=175, and following their motion in time. Two types of numerical experiments were carried out: first, the Eulerian velocity field was evaluated by the direct numerical simulation (DNS) and the particles were advanced in time using the resolved and several filtered velocity fields to study the effect of the subgrid-scale velocity fluctuations on particle motion without the influence of modeling errors. In the second stage, the particle-motion study was performed using independent DNS and large-eddy simulations (LES), thus including the effect of interpolation and subgrid-scale stress modeling errors on the dispersion statistics. At this Reynolds number the small scales were found to have a limited effect on the statistics examined (one-particle dispersion, one-particle velocity autocorrelation, Lagrangian integral time scale, turbulent diffusivity, and two-particles rms dispersion). Only when a significant percentage of the energy was removed from the velocity field some differences were observed between filtered and unfiltered data (especially near the wall). It was found that when the dynamic eddy-viscosity model was employed, modeling errors did not affect the results as much as the filtering itself; the use of the Smagorinsky model, on the other hand, gave inaccurate results. Additional computations for finite-inertia particles have shown that these results represent a conservative estimate, in the sense that actual particles with inertia are less sensitive than the tracer particles examined in the first part of the investigation, and that LES provides improved results when particles with inertia are used.

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.