The two-dimensional problem of gravity-driven laminar flow of a thin layer of fluid down a heated wavy inclined surface is discussed. The coupled effect of bottom topography, variable surface tension and heating has been investigated both analytically and numerically. A stability analysis is conducted while nonlinear simulations are used to validate the stability predictions and also to study thermocapillary effects. The governing equations are based on the Navier–Stokes equations for a thin fluid layer with the cross-stream dependence eliminated by means of a weighted residual technique. Comparisons with experimental data and direct numerical simulations have been carried out and the agreement is good. New interesting results regarding the combined role of surface tension and sinusoidal topography on the stability of the flow are presented. The influence of heating and the Marangoni effect are also deduced.

The method of weighted residuals for thin film flow down an inclined plane is extended to include the effects of bottom waviness, heating, and permeability in this study. A bottom slip condition is used to account for permeability and a constant temperature bottom boundary condition is applied. A weighted residual model (WRM) is derived and used to predict the combined effects of bottom waviness, heating, and permeability on the stability of the flow. In the absence of bottom topography, the results are compared to theoretical predictions from the corresponding Benney equation and also to existing Orr-Sommerfeld predictions. The excellent agreement found indicates that the model does faithfully predict the theoretical critical Reynolds number, which accounts for heating and permeability, and these effects are found to destabilize the flow. Floquet theory is used to investigate how bottom waviness influences the stability of the flow. Finally, numerical simulations of the model equations are also conducted and compared with numerical solutions of the full Navier-Stokes equations for the case with bottom permeability. These results are also found to agree well, which suggests that the WRM remains valid even when permeability is included. V

Internal hydraulic jumps in flows with upstream shear are investigated using two-layer shock-joining theories and numerical solutions of the Navier-Stokes equations. The role of upstream shear has not previously been thoroughly investigated, although it is important in many oceanographic situations, including exchange flows. The full solution spaces of several two-layer theories, distinguished by how dissipation is distributed between the layers, with upstream shear are found, and the physically allowable solution space is identified. These two-layer theories are then evaluated using more realistic numerical simulations that have continuous density and velocity profiles and permit turbulence and mixing. Two-dimensional numerical simulations show that none of the two-layer theories reliably predicts the relation between jump height and speed over the full range of allowable solutions. The numerical simulations also show that different qualitative types of jumps can occur, including undular bores, energy-conserving conjugate state transitions, smooth front jumps with trailing turbulence, and overturning turbulent jumps. Simulation results are used to investigate mixing, which increases with jump height and upstream shear. A few three-dimensional simulations results were undertaken and are in quantitative agreement with the two-dimensional simulations.

The abyssal overturning circulation is thought to be primarily driven by small-scale turbulent mixing. Diagnosed watermass transformations are dominated by rough topography “hotspots”, where the bottom-enhancement of mixing causes the diffusive buoyancy flux to diverge, driving widespread downwelling in the interior—only to be overwhelmed by an even stronger up-welling in a thin Bottom Boundary Layer (BBL). These watermass transformations are significantly underestimated by one-dimensional (1D) sloping boundary layer solutions, suggesting the importance of three-dimensional physics. Here, we use a hierarchy of models to generalize this 1D boundary layer approach to three-dimensional eddying flows over realistically rough topography. When applied to the Mid-Atlantic Ridge in the Brazil Basin, the idealized simulation results are roughly consistent with available observations. Integral buoyancy budgets isolate the physical processes that contribute to realistically strong BBL upwelling. The downwards diffusion of buoyancy is primarily balanced by upwelling along the sloping canyon sidewalls and the surrounding abyssal hills. These flows are strengthened by the restratifying effects of submesoscale baroclinic eddies and by the blocking of along-ridge thermal wind within the canyon. Major topographic sills block along-thalweg flows from restratifying the canyon trough, resulting in the continual erosion of the trough’s stratification. We propose simple modifications to the 1D boundary layer model which approximate each of these three-dimensional effects. These results provide local dynamical insights into mixing-driven abyssal overturning, but a complete theory will also require the non-local coupling to the basin-scale circulation.