 Award ID(s):
 1005363
 NSFPAR ID:
 10218845
 Date Published:
 Journal Name:
 Journal of Fluid Mechanics
 Volume:
 801
 ISSN:
 00221120
 Page Range / eLocation ID:
 459 to 474
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

We perform direct numerical simulations of a gas bubble dissolving in a surrounding liquid. The bubble volume is reduced due to dissolution of the gas, with the numerical implementation of an immersed boundary method, coupling the gas diffusion and the Navier–Stokes equations. The methods are validated against planar and spherical geometries’ analytical moving boundary problems, including the classic Epstein–Plesset problem. Considering a bubble rising in a quiescent liquid, we show that the mass transfer coefficient $k_L$ can be described by the classic Levich formula $k_L = (2/\sqrt {{\rm \pi} })\sqrt {\mathscr {D}_l\,U(t)/d(t)}$ , with $d(t)$ and $U(t)$ the timevarying bubble size and rise velocity, and $\mathscr {D}_l$ the gas diffusivity in the liquid. Next, we investigate the dissolution and gas transfer of a bubble in homogeneous and isotropic turbulence flow, extending Farsoiya et al. ( J. Fluid Mech. , vol. 920, 2021, A34). We show that with a bubble size initially within the turbulent inertial subrange, the mass transfer coefficient in turbulence $k_L$ is controlled by the smallest scales of the flow, the Kolmogorov $\eta$ and Batchelor $\eta _B$ microscales, and is independent of the bubble size. This leads to the nondimensional transfer rate ${Sh}=k_L L^\star /\mathscr {D}_l$ scaling as ${Sh}/{Sc}^{1/2} \propto {Re}^{3/4}$ , where ${Re}$ is the macroscale Reynolds number ${Re} = u_{rms}L^\star /\nu _l$ , with $u_{rms}$ the velocity fluctuations, $L^*$ the integral length scale, $\nu _l$ the liquid viscosity, and ${Sc}=\nu _l/\mathscr {D}_l$ the Schmidt number. This scaling can be expressed in terms of the turbulence dissipation rate $\epsilon$ as ${k_L}\propto {Sc}^{1/2} (\epsilon \nu _l)^{1/4}$ , in agreement with the model proposed by Lamont & Scott ( AIChE J. , vol. 16, issue 4, 1970, pp. 513–519) and corresponding to the high $Re$ regime from Theofanous et al. ( Intl J. Heat Mass Transfer , vol. 19, issue 6, 1976, pp. 613–624).more » « less

We discuss the Onsager theory of wallbounded turbulence, analysing the momentum dissipation anomaly hypothesized by Taylor. Turbulent drag laws observed with both smooth and rough walls imply ultraviolet divergences of velocity gradients. These are eliminated by a coarsegraining operation, filtering out smallscale eddies and windowing out nearwall eddies, thus introducing two arbitrary regularization lengthscales. The regularized equations for resolved eddies correspond to the weak formulation of the Navier–Stokes equation and contain, in addition to the usual turbulent stress, also an inertial drag force modelling momentum exchange with unresolved nearwall eddies. Using an Onsagertype argument based on the principle of renormalization group invariance, we derive an upper bound on wall friction by a function of Reynolds number determined by the modulus of continuity of the velocity at the wall. Our main result is a deterministic version of Prandtl’s relation between the Blasius − 1 / 4 drag law and the 1/7 powerlaw profile of the mean streamwise velocity. At higher Reynolds, the von Kármán–Prandtl drag law requires instead a slow logarithmic approach of velocity to zero at the wall. We discuss briefly also the largeeddy simulation of wallbounded flows and use of iterative renormalization group methods to establish universal statistics in the inertial sublayer. This article is part of the theme issue ‘Scaling the turbulence edifice (part 1)’.more » « less

Abstract We quantify the lake dynamics, near‐bottom turbulence, flux of dissolved oxygen (DO) across the sediment‐water interface (SWI) and their interactions during oxygenation in two lakes. Field observations show that the lake dynamics were modified by the bubble plumes, showing enhanced mixing in the near‐field of the plumes. The interaction of the bubble‐induced flow with the internal density structure resulted in downwelling of warm water into the hypolimnion in the far‐field of the plumes. Within the bottom boundary layer (BBL), both lakes show weak oscillating flows primarily induced by seiching. The vertical profile of mean velocity within 0.4 m above the bed follows a logarithmic scaling. One lake shows a larger drag coefficient than those in stationary BBLs, where the classic law‐of‐the‐wall is valid. The injection of oxygen elevated the water column DO and hence, altered the DO flux across the SWI. The gas transfer velocity is driven by turbulence and is correlated with the bottom shear velocity. The thickness of the diffusive boundary layer was found to be consistent with the Batchelor length scale. The dynamics of the surface renewal time follow a log‐normal distribution, and the turbulent integral time scale is comparable to the surface renewal time. The analyses suggest that the effect of bubble plumes on the BBL turbulence is limited and that the canonical scales of turbulence emerge for the time‐average statistics, validating the turbulence scaling of gas transfer velocity in low‐energy lakes.

We employ numerically implicit subgridscale modeling provided by the wellknown streamlined upwind/Petrov–Galerkin stabilization for the finite element discretization of advection–diffusion problems in a Large Eddy Simulation (LES) approach. Whereas its original purpose was to provide sufficient algorithmic dissipation for a stable and convergent numerical method, more recently, it has been utilized as a subgridscale (SGS) model to account for the effect of small scales, unresolvable by the discretization. The freestream Mach number is 2.5, and direct comparison with a DNS database from our research group, as well as with experiments from the literature of adiabatic supersonic spatially turbulent boundary layers, is performed. Turbulent inflow conditions are generated via our dynamic rescaling–recycling approach, recently extended to highspeed flows. Focus is given to the assessment of the resolved Reynolds stresses. In addition, flow visualization is performed to obtain a much better insight into the physics of the flow. A weak compressibility effect is observed on thermal turbulent structures based on twopoint correlations (IC vs. supersonic). The Reynolds analogy (u′ vs. t′) approximately holds for the supersonic regime, but to a lesser extent than previously observed in incompressible (IC) turbulent boundary layers, where temperature was assumed as a passive scalar. A much longer power law behavior of the mean streamwise velocity is computed in the outer region when compared to the log law at Mach 2.5. Implicit LES has shown very good performance in Mach 2.5 adiabatic flat plates in terms of the mean flow (i.e., Cf and UVD+). iLES significantly overpredicts the peak values of u′, and consequently Reynolds shear stress peaks, in the buffer layer. However, excellent agreement between the turbulence intensities and Reynolds shear stresses is accomplished in the outer region by the present iLES with respect to the external DNS database at similar Reynolds numbers.more » « less

The inertial subrange of turbulent scales is commonly reflected by a power law signature in ensemble statistics such as the energy spectrum and structure functions – both in theory and from observations. Despite promising findings on the topic of fractal geometries in turbulence, there is no accepted image for the physical flow features corresponding to this statistical signature in the inertial subrange. The present study uses boundary layer turbulence measurements to evaluate the selfsimilar geometric properties of velocity isosurfaces and investigate their influence on statistics for the velocity signal. The fractal dimension of streamwise velocity isosurfaces, indicating statistical selfsimilarity in the size of ‘wrinkles’ along each isosurface, is shown to be constant only within the inertial subrange of scales. For the transition between the inertial subrange and production range, it is inferred that the largest wrinkles become increasingly confined by the overall size of largescale coherent velocity regions such as uniform momentum zones. The selfsimilarity of isosurfaces yields powerlaw trends in subsequent onedimensional statistics. For instance, the theoretical 2/3 powerlaw exponent for the structure function can be recovered by considering the collective behaviour of numerous isosurface level sets. The results suggest that the physical presence of inertial subrange eddies is manifested in the selfsimilar wrinkles of isosurfaces.more » « less