Numerical Simulation of Small-Scale Turbulent Diffusion for Injected CO2 in the Deep Ocean

Author(s):  
Shinichiro Hirabayashi ◽  
Toru Sato

Small-scale turbulent diffusion is numerically simulated in relation to its time variability dependent on tidal current. From time series of 3D (3-dimensional) velocity and temperature measured at the depth of 2000 m, main constituents of internal tides are extracted. Wavenumbers of eddies, the scale of which is O (102) m, are determined by Taylor’s frozen eddy hypothesis with time-varying tidal speed taken as advection velocity. Large eddy simulation was applied to generate turbulence numerically by forcing the measured low-wavenumber components in the computational domain. The result shows that the energy dissipation rate varies in time, ranging from O (10−11) to O (10−7) m2 s−3. These values of energy dissipation correspond with vertical diffusivity of 10−6 to 10−4 m2 s−1 in this area, with temporal average value of 6.9 × 10−5 m2 s−1.

Author(s):  
Shinichiro Hirabayashi ◽  
Toru Sato

In this study, vertical diffusivity, the scale of which was O (10 m), at a particular site in the deep ocean was estimated by using numerical simulations with forcing low-wavenumber components, which had been reproduced from measurement data. Spatial information of velocity field was reproduced by spectral analysis of 4 sets of time-series measured simultaneously at different places in the real ocean. In order to estimate finer-scale structures, which are necessary to obtain statistical quantities such as energy dissipation rate, large eddy simulations were carried out with forcing low-wavenumber components of velocity reproduced in the spectral analysis. The low-wavenumber components generated by the nonlinear interaction of forced components and resolved components were successfully removed from the computational domain by introducing a partial spectral filter in place of the conventional FFT filter. Vertical diffusivity was estimated by using the energy dissipation rate of the reproduced flow field, which was 3.3×10−5 m2s−1 on the time average.


2012 ◽  
Vol 7 (1) ◽  
pp. 53-69
Author(s):  
Vladimir Dulin ◽  
Yuriy Kozorezov ◽  
Dmitriy Markovich

The present paper reports PIV (Particle Image Velocimetry) measurements of turbulent velocity fluctuations statistics in development region of an axisymmetric free jet (Re = 28 000). To minimize measurement uncertainty, adaptive calibration, image processing and data post-processing algorithms were utilized. On the basis of theoretical analysis and direct measurements, the paper discusses effect of PIV spatial resolution on measured statistical characteristics of turbulent fluctuations. Underestimation of the second-order moments of velocity derivatives and of the turbulent kinetic energy dissipation rate due to a finite size of PIV interrogation area and finite thickness of laser sheet was analyzed from model spectra of turbulent velocity fluctuations. The results are in a good agreement with the measured experimental data. The paper also describes performance of possible ways to account for unresolved small-scale velocity fluctuations in PIV measurements of the dissipation rate. In particular, a turbulent viscosity model can be efficiently used to account for the unresolved pulsations in a free turbulent flow


2008 ◽  
Vol 38 (6) ◽  
pp. 1296-1312 ◽  
Author(s):  
Johannes R. Gemmrich ◽  
Michael L. Banner ◽  
Chris Garrett

Abstract Video observations of the ocean surface taken from aboard the Research Platform FLIP reveal the distribution of the along-crest length and propagation velocity of breaking wave crests that generate visible whitecaps. The key quantity assessed is Λ(c)dc, the average length of breaking crests per unit area propagating with speeds in the range (c, c + dc). Independent of the wave field development, Λ(c) is found to peak at intermediate wave scales and to drop off sharply at larger and smaller scales. In developing seas breakers occur at a wide range of scales corresponding to phase speeds from about 0.1 cp to cp, where cp is the phase speed of the waves at the spectral peak. However, in developed seas, breaking is hardly observed at scales corresponding to phase speeds greater than 0.5 cp. The phase speed of the most frequent breakers shifts from 0.4 cp to 0.2 cp as the wave field develops. The occurrence of breakers at a particular scale as well as the rate of surface turnover are well correlated with the wave saturation. The fourth and fifth moments of Λ(c) are used to estimate breaking-wave-supported momentum fluxes, energy dissipation rate, and the fraction of momentum flux supported by air-entraining breaking waves. No indication of a Kolmogorov-type wave energy cascade was found; that is, there is no evidence that the wave energy dissipation is dominated by small-scale waves. The proportionality factor b linking breaking crest distributions to the energy dissipation rate is found to be (7 ± 3) × 10−5, much smaller than previous estimates.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Chu-Fang Yang ◽  
Wu-Cheng Chi ◽  
Hans van Haren

AbstractTurbulent mixing in the deep ocean is not well understood. The breaking of internal waves on sloped seafloor topography can generate deep-sea turbulence. However, it is difficult to measure turbulence comprehensively due to its multi-scale processes, in addition to flow–flow and flow–topography interactions. Dense, high-resolution spatiotemporal coverage of observations may help shed light on turbulence evolution. Here, we present turbulence observations from four broadband ocean bottom seismometers (OBSs) and a 200-m vertical thermistor string (T-string) in a footprint of 1 × 1 km to characterize turbulence induced by internal waves at a depth of 3000 m on a Pacific continental slope. Correlating the OBS-calculated time derivative of kinetic energy and the T-string-calculated turbulent kinetic energy dissipation rate, we propose that the OBS-detected signals were induced by near-seafloor turbulence. Strong disturbances were detected during a typhoon period, suggesting large-scale inertial waves breaking with upslope transport speeds of 0.2–0.5 m s−1. Disturbances were mostly excited on the downslope side of the array where the internal waves from the Pacific Ocean broke initially and the turbulence oscillated between < 1 km small-scale ridges. Such small-scale topography caused varying turbulence-induced signals due to localized waves breaking. Arrayed OBSs can provide complementary observations to characterize deep-sea turbulence.


2006 ◽  
Vol 63 (4) ◽  
pp. 1291-1307 ◽  
Author(s):  
B. B. Balsley ◽  
R. G. Frehlich ◽  
M. L. Jensen ◽  
Y. Meillier

Abstract Some 50 separate high-resolution profiles of small-scale turbulence defined by the energy dissipation rate (ɛ), horizontal wind speed, and temperature from near the surface, through the nighttime stable boundary layer (SBL), and well into the residual layer are used to compare the various definitions of SBL height during nighttime stable conditions. These profiles were obtained during postmidnight periods on three separate nights using the Tethered Lifting System (TLS) during the Cooperative Atmosphere–Surface Exchange Study (CASES-99) campaign in east-central Kansas, October 1999. Although the number of profiles is insufficient to make any definitive conclusions, the results suggest that, under most conditions, the boundary layer top can be reasonably estimated in terms of a very significant decrease in the energy dissipation rate (i.e., the mixing height) with height. In the majority of instances this height lies slightly below the height of a pronounced minimum in wind shear and slightly above a maximum in N 2, where N is the Brunt–Väisälä frequency. When combined with flux measurements and vertical velocity variance data obtained from the nearby 55-m tower, the results provide additional insights into SBL processes, even when the boundary layer, by any definition, extends to heights well above the top of the tower. Both the TLS profiles and tower data are then used for preliminary high-resolution studies into various categories of SBL structure, including the so-called upside-down boundary layer.


The multifractal formalism for the eulerian statistics of small-scale dynamics in turbulent flows is reviewed. Theoretical extensions of these results (the statistics of small volume averages of the energy dissipation rate) are used to predict properties of the probability distribution of the local energy dissipation rate at a fixed point. The improved parametrization of the eulerian statistics allows the lagrangian statistics (those for a fixed fluid particle in contrast to the eulerian statistics at a fixed point) to be determined exactly by using results derived as a consequence of incompressibility. Several properties of particle trajectories in a turbulent flow can be predicted with these new lagrangian statistics. In particular, a trajectory is typically smooth and generally unremarkable in its features. This contrasts the often suggested description: that of a highly convoluted and intricately structured ‘fractal’ curve. Some of the traditional dispersion results, which depend on the lagrangian statistics, are shown to be only weakly influenced by the intermittency inherent in the multifractal character of turbulence.


2019 ◽  
Vol 10 (1) ◽  
Author(s):  
Clément Vic ◽  
Alberto C. Naveira Garabato ◽  
J. A. Mattias Green ◽  
Amy F. Waterhouse ◽  
Zhongxiang Zhao ◽  
...  

2015 ◽  
Vol 784 ◽  
pp. 109-129 ◽  
Author(s):  
S. L. Tang ◽  
R. A. Antonia ◽  
L. Djenidi ◽  
Y. Zhou

The transport equation for the isotropic turbulent energy dissipation rate $\overline{{\it\epsilon}}_{iso}$ along the centreline in the far-wake of a circular cylinder is derived by applying the limit at small separations to the two-point energy budget equation. It is found that the imbalance between the production and the destruction of $\overline{{\it\epsilon}}_{iso}$, respectively due to vortex stretching and viscosity, is governed by both the streamwise advection and the lateral turbulent diffusion (the former contributes more to the budget than the latter). This imbalance differs intrinsically from that in other flows, e.g. grid turbulence and the flow along the centreline of a fully developed channel, where either the streamwise advection or the lateral turbulent diffusion of $\overline{{\it\epsilon}}_{iso}$ governs the imbalance. More importantly, the different types of imbalance represent different constraints on the relation between the skewness of the longitudinal velocity derivative $S$ and the destruction coefficient of enstrophy $G$. This results in a non-universal approach of $S$ towards a constant value as the Taylor microscale Reynolds number $R_{{\it\lambda}}$ increases. For the present flow, the magnitude of $S$ decreases initially ($R_{{\it\lambda}}\leqslant 40$) before increasing ($R_{{\it\lambda}}>40$) towards this constant value. The constancy of $S$ at large $R_{{\it\lambda}}$ violates the modified similarity hypothesis introduced by Kolmogorov (J. Fluid Mech., vol. 13, 1962, pp. 82–85) but is consistent with the original similarity hypotheses (Kolmogorov, Dokl. Akad. Nauk SSSR, vol. 30, 1941b, pp. 299–303 (see also 1991 Proc. R. Soc. Lond. A, vol. 434, pp. 9–13)) ($K41$), and, more importantly, with the almost completely self-preserving nature of the plane far-wake.


2021 ◽  
Vol 932 ◽  
Author(s):  
M.A. Gorokhovski ◽  
S.K. Oruganti

In this work we introduce a Lagrangian stochastic model for particle motion and evaporation to be used in large-eddy simulations (LES) of turbulent liquid sprays. Effects of small-scale intermittency, usually under-resolved in LES, are explicitly included via modelling of the energy dissipation rate seen by a droplet along its trajectory. Namely, the dissipation rate is linked to the norm of the droplet sub-filtered acceleration which is included in the droplet motion equation. This norm, along with the direction of the droplet sub-filtered acceleration, is simulated as a stochastic process. With increasing Reynolds number, the distribution of the sub-filtered acceleration develops longer tails, with a slower decay in auto-correlation functions of the norm and direction of this acceleration. The stochastic models are specified for particles larger and smaller the Kolmogorov length scale. The assumption of the droplet evaporation model is similar, i.e. the evaporation rate is strongly enhanced when a droplet is subjected to very localized zones of intense velocity gradients. Thereby, the overall evaporation process is assumed to be a succession of two steady-state sub-processes with equal intensities, i.e. evaporation and vapour mixing. Then the stochastic properties of the overall evaporation rate are also controlled by fluctuations of the energy dissipation rate along the droplet path, and with increasing Reynolds number, the intensity of fluctuations of this rate is also increasing. The assessment of the presented stochastic models in LES of high-speed non-evaporating and evaporating sprays show the accurate prediction of experimental data on relatively coarser grids along with a remarkably weaker sensitivity to the grid spacing. The joint statistics and Voronoi tessellations exhibit strong intermittency of evaporation rate. The intensity of turbulence along the droplet pathway substantially promotes the vaporization rate.


Sign in / Sign up

Export Citation Format

Share Document