Gregory P. Chini

Langmuir–Submesoscale Interactions: Descriptive Analysis of Multiscale Frontal Spindown Simulations

AbstractThe interactions between boundary layer turbulence, including Langmuir turbulence, and submesoscale processes in the oceanic mixed layer are described using large-eddy simulations of the spindown of a temperature front in the presence of submesoscale eddies, winds, and waves. The simulations solve the surface-wave-averaged Boussinesq equations with Stokes drift wave forcing at a resolution that is sufficiently fine to capture small-scale Langmuir turbulence. A simulation without Stokes drift forcing is also performed for comparison. Spatial and spectral properties of temperature, velocity, and vorticity fields are described, and these fields are scale decomposed in order to examine multiscale fluxes of momentum and buoyancy. Buoyancy flux results indicate that Langmuir turbulence counters the restratifying effects of submesoscale eddies, leading to small-scale vertical transport and mixing that is 4 times greater than in the simulations without Stokes drift forcing. The observed fluxes are also sh...

Reduced description of exact coherent states in parallel shear flows.

A reduced description of exact coherent structures in the transition regime of plane parallel shear flows is developed, based on the Reynolds number scaling of streamwise-averaged (mean) and streamwise-varying (fluctuation) velocities observed in numerical simulations. The resulting system is characterized by an effective unit Reynolds number mean equation coupled to linear equations for the fluctuations, regularized by formally higher-order diffusion. Stationary coherent states are computed by solving the resulting equations simultaneously using a robust numerical algorithm developed for this purpose. The algorithm determines self-consistently the amplitude of the fluctuations for which the associated mean flow is just such that the fluctuations neither grow nor decay. The procedure is used to compute exact coherent states of a flow introduced by Drazin and Reid [Hydrodynamic Stability (Cambridge University Press, Cambridge, UK, 1981)] and studied by Waleffe [Phys. Fluids 9, 883 (1997)]: a linearly stable, plane parallel shear flow confined between stationary stress-free walls and driven by a sinusoidal body force. Numerical continuation of the lower-branch states to lower Reynolds numbers reveals the presence of a saddle node; the saddle node allows access to upper-branch states that are, like the lower-branch states, self-consistently described by the reduced equations. Both lower- and upper-branch states are characterized in detail.

An asymptotically reduced model of turbulent Langmuir circulation

An asymptotically exact reduced model describing Langmuir circulation (LC), an upper ocean convective flow, in the strong surface wave forcing limit is constructed. Linear and secondary stability analyses suggest that the model captures correctly the dominant instability modes arising in strongly anisotropic, turbulent LC. The model is conceptually simpler, and may be more amenable to numerical simulation, than the three-dimensional Craik–Leibovich equations from which the reduced equation set is derived.

Large Rayleigh number thermal convection: heat flux predictions and strongly nonlinear solutions

We investigate the structure of strongly nonlinear Rayleigh–Benard convection cells in the asymptotic limit of large Rayleigh number and fixed, moderate Prandtl number. Unlike the flows analyzed in prior theoretical studies of infinite Prandtl number convection, our cellular solutions exhibit dynamically inviscid constant-vorticity cores. By solving an integral equation for the cell-edge temperature distribution, we are able to predict, as a function of cell aspect ratio, the value of the core vorticity, details of the flow within the thin boundary layers and rising/falling plumes adjacent to the edges of the convection cell, and, in particular, the bulk heat flux through the layer. The results of our asymptotic analysis are corroborated using full pseudospectral numerical simulations and confirm that the heat flux is maximized for convection cells that are roughly square in cross section.

Prospectus: towards the development of high-fidelity models of wall turbulence at large Reynolds number

Recent and on-going advances in mathematical methods and analysis techniques, coupled with the experimental and computational capacity to capture detailed flow structure at increasingly large Reynolds numbers, afford an unprecedented opportunity to develop realistic models of high Reynolds number turbulent wall-flow dynamics. A distinctive attribute of this new generation of models is their grounding in the Navier–Stokes equations. By adhering to this challenging constraint, high-fidelity models ultimately can be developed that not only predict flow properties at high Reynolds numbers, but that possess a mathematical structure that faithfully captures the underlying flow physics. These first-principles models are needed, for example, to reliably manipulate flow behaviours at extreme Reynolds numbers. This theme issue of Philosophical Transactions of the Royal Society A provides a selection of contributions from the community of researchers who are working towards the development of such models. Broadly speaking, the research topics represented herein report on dynamical structure, mechanisms and transport; scale interactions and self-similarity; model reductions that restrict nonlinear interactions; and modern asymptotic theories. In this prospectus, the challenges associated with modelling turbulent wall-flows at large Reynolds numbers are briefly outlined, and the connections between the contributing papers are highlighted. This article is part of the themed issue ‘Toward the development of high-fidelity models of wall turbulence at large Reynolds number’.

Two-dimensional streaming flows in high-intensity discharge lamps

High-intensity discharge (HID) lamps embody a practical application in which acoustically generated streaming flows are used to significantly improve energy efficiency. Streaming in these lamps is examined using finite-element simulations in conjunction with available experimental results on the basis of the assumption that the streaming motion is excited by two-dimensional acoustic standing waves. Neither the magnitude nor the direction of the time-averaged flows is adequately explained by existing theory. Consequently, a modified streaming analysis is proposed in which the fluctuating flow is driven by an oscillating pressure field rather by a moving boundary and convective terms in both the instantaneous and streaming flows are included. Density variations are also shown to be important to the generation of the observed and simulated streaming. This analysis highlights the differences between streaming flows in HID lamps and those described in canonical problems appearing elsewhere in the literature.

A self-sustaining process model of inertial layer dynamics in high Reynolds number turbulent wall flows

Field observations and laboratory experiments suggest that at high Reynolds numbers Re the outer region of turbulent boundary layers self-organizes into quasi-uniform momentum zones (UMZs) separated by internal shear layers termed ‘vortical fissures’ (VFs). Motivated by this emergent structure, a conceptual model is proposed with dynamical components that collectively have the potential to generate a self-sustaining interaction between a single VF and adjacent UMZs. A large-Re asymptotic analysis of the governing incompressible Navier–Stokes equation is performed to derive reduced equation sets for the streamwise-averaged and streamwise-fluctuating flow within the VF and UMZs. The simplified equations reveal the dominant physics within—and isolate possible coupling mechanisms among—these different regions of the flow. This article is part of the themed issue ‘Toward the development of high-fidelity models of wall turbulence at large Reynolds number’.

Modulated patterns in a reduced model of a transitional shear flow

We consider a close relative of plane Couette flow called Waleffe flow in which the fluid is confined between two free-slip walls and the flow driven by a sinusoidal force. We use a reduced model of such flows constructed elsewhere to compute stationary exact coherent structures of Waleffe flow in periodic domains with a large spanwise period. The computations reveal the emergence of stationary states exhibiting strong amplitude and wavelength modulation in the spanwise direction. These modulated states lie on branches exhibiting complex dependence on the Reynolds number but no homoclinic snaking.

Exact coherent structures in an asymptotically reduced description of parallel shear flows

A reduced description of shear flows motivated by the Reynolds number scaling of lower-branch exact coherent states in plane Couette flow (Wang J, Gibson J and Waleffe F 2007 Phys. Rev. Lett. 98 204501) is constructed. Exact time-independent nonlinear solutions of the reduced equations corresponding to both lower and upper branch states are found for a sinusoidal, body-forced shear flow. The lower branch solution is characterized by fluctuations that vary slowly along the critical layer while the upper branch solutions display a bimodal structure and are more strongly focused on the critical layer. The reduced equations provide a rational framework for investigations of subcritical spatiotemporal patterns in parallel shear flows.

Reduced modeling of porous media convection in a minimal flow unit at large Rayleigh number

Abstract Direct numerical simulations (DNS) indicate that at large values of the Rayleigh number (Ra) convection in porous media self-organizes into narrowly-spaced columnar flows, with more complex spatiotemporal features being confined to boundary layers near the top and bottom walls. In this investigation of high-Ra porous media convection in a minimal flow unit, two reduced modeling strategies are proposed that exploit these specific flow characteristics. Both approaches utilize the idea of decomposition since the flow exhibits different dynamics in different regions of the domain: small-scale cellular motions generally are localized within the thermal and vorticity boundary layers near the upper and lower walls, while in the interior, the flow exhibits persistent large-scale structures and only a few low (horizontal) wavenumber Fourier modes are active. Accordingly, in the first strategy, the domain is decomposed into two near-wall regions and one interior region. Our results confirm that suppressing the interior high-wavenumber modes has negligible impact on the essential structural features and transport properties of the flow. In the second strategy, a hybrid reduced model is constructed by using Galerkin projection onto a fully a priori eigenbasis drawn from energy stability and upper bound theory, thereby extending the model reduction strategy developed by Chini et al. (2011) [45] to large Ra. The results indicate that the near-wall upper-bound eigenmodes can economically represent the small-scale rolls within the exquisitely-thin thermal boundary layers. Relative to DNS, the hybrid algorithm enables over an order-of-magnitude increase in computational efficiency with only a modest loss of accuracy.