Model and Experimental Estimates of Vertical Mixing Intensity in the Sea Upper Homogeneous Layer

2021 ◽  
Vol 37 (3) ◽  
Author(s):  
A. M. Chukharev ◽  
M. I. Pavlov ◽  
◽  

Purpose. The study is aimed at qualitative and quantitative analysis (based on the updated previously proposed multiscale model) of the experimental data on turbulence intensity and their comparison with theoretical and semi-empirical relationships for the purpose of describing the contributions of various turbulence sources. Methods and Results. A comparative analysis of experimental data and model calculations of turbulence characteristics near the sea surface was performed. The methods of theoretical assessing generation of turbulence in the near-surface sea layer by various physical processes are considered. The results of calculations by the well-known models of turbulent exchange were compared with the experimental data collected by the scientists of the Turbulence Department of MHI, RAS, using the specialized equipment. The analysis results made it possible to determine the possibility of applying the considered models for calculating turbulence intensity under different hydrometeorological conditions. At light winds, none of the models yielded the results which matched the measurement data. At moderate winds, the simulation results showed quite satisfactory agreement with the experiment data; and for strong winds, the multiscale model results were the best. This model was modified to assess the contributions of two other mechanisms of turbulence generation: the Stokes drift and the Langmuir circulations. Conclusions. Objective assessment of the turbulent exchange intensity requires taking into account of three main mechanisms of turbulence generation, namely flow velocity shear, wave motions and wave breaking. Depending on the hydrometeorological situation, each of these mechanisms can dominate in a certain depth range. The calculations performed using the updated model showed that the Stokes drift added 2–17 % to the total dissipation in the upper 30-meter layer, whereas the contribution of the Langmuir circulations calculated through dependence of the vertical velocity of kinetic energy transfer upon the Langmuir number, can reach 15 % for small Langmuir numbers.

2021 ◽  
Vol 28 (3) ◽  
Author(s):  
A. M. Chukharev ◽  
M. I. Pavlov ◽  
◽  

Purpose. The study is aimed at qualitative and quantitative analysis (based on the updated previously proposed multiscale model) of the experimental data on turbulence intensity and their comparison with theoretical and semi-empirical relationships for the purpose of describing the contributions of various turbulence sources. Methods and Results. A comparative analysis of experimental data and model calculations of turbulence characteristics near the sea surface was performed. The methods of theoretical assessing generation of turbulence in the near-surface sea layer by various physical processes are considered. The results of calculations by the well-known models of turbulent exchange were compared with the experimental data collected by the scientists of the Turbulence Department of MHI, RAS, using the specialized equipment. The analysis results made it possible to determine the possibility of applying the considered models for calculating turbulence intensity under different hydrometeorological conditions. At light winds, none of the models yielded the results which matched the measurement data. At moderate winds, the simulation results showed quite satisfactory agreement with the experiment data; and for strong winds, the multiscale model results were the best. This model was modified to assess the contributions of two other mechanisms of turbulence generation: the Stokes drift and the Langmuir circulations. Conclusions. Objective assessment of the turbulent exchange intensity requires taking into account of three main mechanisms of turbulence generation, namely flow velocity shear, wave motions and wave breaking. Depending on the hydrometeorological situation, each of these mechanisms can dominate in a certain depth range. The calculations performed using the updated model showed that the Stokes drift added 2–17 % to the total dissipation in the upper 30-meter layer, whereas the contribution of the Langmuir circulations calculated through dependence of the vertical velocity of kinetic energy transfer upon the Langmuir number, can reach 15 % for small Langmuir numbers.


2021 ◽  
Vol 28 (4) ◽  
Author(s):  
A. S. Samodurov ◽  
A. M. Chukharev ◽  
D. A. Kazakov ◽  
◽  
◽  
...  

Purpose. The purpose of the study is to assess the coefficient of vertical turbulent exchange for different layers of the Black Sea basin based on the experimental data on microstructure of the physical fields obtained for the period 2004–2019 in the Black Sea and using the semi-empirical models. Methods and Results. For the upper mixed layer, the turbulent energy dissipation rate ɛ and the exchange coefficient were calculated using the velocity fluctuation spectra based on the Kolmogorov hypotheses on the turbulence spectrum inertial range. In the stratified layers, the turbulence coefficient and the dissipation rate were experimentally determined both from the spectra of the velocity horizontal fluctuations’ gradients and the vertical spectra of temperature fluctuations using the concept of the effective scale of turbulent patches. Depending on the features of the hydrological regime and the prevailing energy contributors to turbulence generation, five layers were identified and described (including their characteristic power dependences of the vertical turbulent diffusion coefficients K on the buoyancy frequency N) using the 1.5D-model of vertical turbulent exchange for the basin under study. For the stratified layers, the 1.5D-model results were comparatively analyzed with those of the other semi-empirical and theoretical models describing the most probable hydrophysical processes in each specific layer; the relations for the vertical turbulent exchange coefficient were obtained depending on the buoyancy frequency. Conclusions. Comparison of the experimental data collected under different hydrometeorological conditions with the simulations resulted from the known turbulence models for the sea upper layer showed that the best agreement between the simulation and measurement data was provided by a multiscale model taking into account three basic mechanisms of turbulence generation: current velocity shear, instability of wave motions, and wave breaking. The turbulent exchange coefficient dependencies on depth are conditioned by the effect of the turbulence dominant source at a given level. In the stratified layers, the exchange coefficient dependence on buoyancy frequency is determined by the hydrophysical processes in each layer; the relations obtained for individual layers indicate intensity of the contributions of vertical advection, internal wave breakings, turbulence diffusion and geothermal flux.


2014 ◽  
Vol 44 (10) ◽  
pp. 2739-2752 ◽  
Author(s):  
A. Basovich

Abstract A new mechanism of instability leading to development of Langmuir circulations is proposed and studied based on the hypothesis that the turbulence and, correspondingly, the eddy viscosity are reduced in regions of higher than average contaminant concentration. Here, bubbles are considered as the contaminant, although it is known that surfactants and some particles also are capable of turbulence reduction. The analysis shows that only a very small local decrease in eddy viscosity is needed to initiate the instability. Simplifications to the momentum and bubble turbulence models, as well as neglect of vertical advection, make it possible to analytically solve the perturbation equations and determine the characteristic scale with the maximum growth rate. The scale of the fastest-growing Langmuir circulations is found to be a function of the concentration of bubbles in the near-surface layer, the surface current shear (wind shear), the Stokes drift created by the surface waves, and the eddy viscosity. In contrast to the results of earlier models, the analysis predicts that the maximum change in the current velocity along the direction of the wind (the so-called jets and wakes) is at the surface, not below it. The ratio of perturbation of along-wind surface current and vertical velocity generated by circulations (the pitch) and the aspect ratio of the Langmuir rolls are in reasonable agreement with the experimental data.


2021 ◽  
Vol 37 (4) ◽  
Author(s):  
A. S. Samodurov ◽  
A. M. Chukharev ◽  
D. A. Kazakov ◽  
◽  
◽  
...  

Purpose. The purpose of the study is to assess the coefficient of vertical turbulent exchange for different layers of the Black Sea basin based on the experimental data on microstructure of the physical fields obtained for the period 2004–2019 in the Black Sea and using the semi-empirical models. Methods and Results. For the upper mixed layer, the turbulent energy dissipation rate ɛ and the exchange coefficient were calculated using the velocity fluctuation spectra based on the Kolmogorov hypotheses on the turbulence spectrum inertial range. In the stratified layers, the turbulence coefficient and the dissipation rate were experimentally determined both from the spectra of the velocity horizontal fluctuations’ gradients and the vertical spectra of temperature fluctuations using the concept of the effective scale of turbulent spots. Depending on the features of the hydrological regime and the prevailing energy contributors to turbulence generation, five layers were identified and described (including their characteristic power dependences of the vertical turbulent diffusion coefficients K on the buoyancy frequency N) using the 1.5D-model of vertical turbulent exchange for the basin under study. For the stratified layers, the 1.5D-model results were comparatively analyzed with those of the other semi-empirical and theoretical models describing the most probable hydrophysical processes in each specific layer; the relations for the vertical turbulent exchange coefficient were obtained depending on the buoyancy frequency. Conclusions. Comparison of the experimental data collected under different hydrometeorological conditions with the simulations resulted from the known turbulence models for the sea upper layer showed that the best agreement between the simulation and measurement data was provided by a multiscale model taking into account three basic mechanisms of turbulence generation: current velocity shear, instability of wave motions, and wave breaking. The turbulent exchange coefficient dependencies on depth are conditioned by the affect of the turbulence dominant source at a given level. In the stratified layers, the exchange coefficient dependence on buoyancy frequency is determined by the hydrophysical processes in each layer; the relations obtained for individual layers indicate intensity of the contributions of vertical advection, internal wave breakings, turbulence diffusion and geothermal flux.


Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCB35-WCB46 ◽  
Author(s):  
Hua-wei Zhou ◽  
Peiming Li ◽  
Zhihui Yan ◽  
Hui Liu

Although first-arrival tomography provides an effective way to estimate near-surface velocities and static corrections, the undulation of velocity interfaces such as the base of the weathered zone may not be easily determined by this method. The main reason is that first arrivals are insensitive to small geometric changes in velocity interfaces because their raypaths tend to traverse along those interfaces. To improve the solution of interface geometry, we developed a deformable layer tomostatics method that approximates the near-surface velocity field as several layers of constant velocity and variable thickness that can be inverted for the geometry of the velocity interfaces. We use a multiscale model parameterization in the inversion for interface geometry. Synthetic and field data tests showed that the method can determine the interface geometry. Constraining the depth range of the basal boundary of the weathered zone increases the convergence rate of the iterative inversion process. Tests on field data showed greater reflection coherency in a stacked section based on constrained static corrections than in one from unconstrained static corrections. The method yielded a better match with statics computed from sand-dune curves than does a match obtained by using two commercial grid tomography packages.


2011 ◽  
Vol 133 (5) ◽  
Author(s):  
F. P. P. Tan ◽  
N. B. Wood ◽  
G. Tabor ◽  
X. Y. Xu

In this study, two different turbulence methodologies are investigated to predict transitional flow in a 75% stenosed axisymmetric experimental arterial model and in a slightly modified version of the model with an eccentric stenosis. Large eddy simulation (LES) and Reynolds-averaged Navier–Stokes (RANS) methods were applied; in the LES simulations eddy viscosity subgrid-scale models were employed (basic and dynamic Smagorinsky) while the RANS method involved the correlation-based transitional version of the hybrid k-ε/k-ω flow model. The RANS simulations used 410,000 and 820,000 element meshes for the axisymmetric and eccentric stenoses, respectively, with y+ less than 2 viscous wall units for the boundary elements, while the LES used 1,200,000 elements with y+ less than 1. Implicit filtering was used for LES, giving an overlap between the resolved and modeled eddies, ensuring accurate treatment of near wall turbulence structures. Flow analysis was carried out in terms of vorticity and eddy viscosity magnitudes, velocity, and turbulence intensity profiles and the results were compared both with established experimental data and with available direct numerical simulations (DNSs) from the literature. The simulation results demonstrated that the dynamic Smagorinsky LES and RANS transitional model predicted fairly comparable velocity and turbulence intensity profiles with the experimental data, although the dynamic Smagorinsky model gave the best overall agreement. The present study demonstrated the power of LES methods, although they were computationally more costly, and added further evidence of the promise of the RANS transition model used here, previously tested in pulsatile flow on a similar model. Both dynamic Smagorinsky LES and the RANS model captured the complex transition phenomena under physiological Reynolds numbers in steady flow, including separation and reattachment. In this respect, LES with dynamic Smagorinsky appeared more successful than DNS in replicating the axisymmetric experimental results, although inflow conditions, which are subject to caveats, may have differed. For the eccentric stenosis, LES with Smagorinsky coefficient of 0.13 gave the closest agreement with DNS despite the known shortcomings of fixed coefficients. The relaminarization as the flow escaped the influence of the stenosis was amply demonstrated in the simulations, graphically so in the case of LES.


2021 ◽  
Vol 11 (7) ◽  
Author(s):  
Nyakno Jimmy George

AbstractAVI (Aquifer vulnerability index), GOD (groundwater occurrence, overlying lithology and depth to the aquifer), GLSI (geo-electric layer susceptibility indexing) and S (longitudinal unit conductance) models were used to assess economically exploitable groundwater resource in the coastal environment of Akwa Ibom State, southern Nigeria. The models were employed in order to delineate groundwater into its category of vulnerability to contamination sources using the first- and second-order geo-electric indices as well as hydrogeological inputs. Vertical electrical sounding technique employing Schlumberger electrode configuration was carried out in 16 locations, close to logged boreholes with known aquifer core samples. Primary or first-order geo-electric indices (resistivity, thickness and depth) measured were used to determine S. The estimated aquifer hydraulic conductivity, K, calculated from grain size diameter and water resistivity values were used to calculate hydraulic resistance (C) used to estimate AVI. With the indices assigned to geo-electric parameters on the basis of their influences, GOD and FSLI were calculated using appropriate equations. The geologic sequence in the study area consists of geo-electric layers ranging from motley topsoil, argillites (clayey to fine sands) and arenites (medium to gravelly sands). Geo-electric parametric indices of aquifer overlying layers across the survey area were utilized to weigh the vulnerability of the underlying water-bearing resource to the contaminations from surface and near-surface, using vulnerability maps created. Geo-electrically derived model maps reflecting AVI, BOD, FLSI and S were compared to assess their conformity to the degree of predictability of groundwater vulnerability. The AVI model map shows range of values of log C ( −3.46—0.07) generally less than unity and hence indicating high vulnerability. GOD model tomographic map displays a range of 0.1–0.3, indicating that the aquifer with depth range of 20.5 to 113.1 m or mean depth of 72. 3 m is lowly susceptible to surface and near-surface impurities. Again, the FLSI map displays a range of FLSI index of 1.25 to 2.75, alluding that the aquifer underlying the protective layer has a low to moderate vulnerability. The S model has values ranging from 0.013 to 0.991S. As the map indicates, a fractional portion of the aquifer at the western (Ikot Abasi) part of the study area has moderate to good protection (moderate vulnerability) while weak to poor aquifer protection (high vulnerability) has poor protection. The S model in this analysis seems to overstate the degree of susceptibility to contamination than the AVI, GOD and GLSI models. From the models, the categorization of severity of aquifer vulnerability to contaminations is relatively location-dependent and can be assessed through the model tomographic maps generated.


2018 ◽  
Vol 48 (1) ◽  
pp. 117-122 ◽  
Author(s):  
Brodie Pearson

AbstractThis study shows that the presence of Stokes drift us in the turbulent upper ocean induces a near-surface Eulerian current that opposes the Stokes drift. This current is distinct from previously studied anti-Stokes currents because it does not rely on the presence of planetary rotation or mean lateral gradients. Instead, the anti-Stokes flow arises from an interaction between the Stokes drift and turbulence. The new anti-Stokes flow is antiparallel to us near the ocean surface, is parallel to us at depth, and integrates to zero over the depth of the boundary layer. The presence of Stokes drift in large-eddy simulations (LES) is shown to induce artificial energy production caused by a combination of the new anti-Stokes flow and LES numerics. As a result, care must be taken when designing and interpreting simulations of realistic wave forcing, particularly as rotation becomes weak and/or us becomes perpendicular to the surface wind stress. The mechanism of the artificial energy production is demonstrated for a generalized LES subgrid scheme.


Author(s):  
M. Stripf ◽  
A. Schulz ◽  
H.-J. Bauer ◽  
S. Wittig

Two extended models for the calculation of rough wall transitional boundary layers with heat transfer are presented. Both models comprise a new transition onset correlation, which accounts for the effects of roughness height and density, turbulence intensity and wall curvature. In the transition region, an intermittency equation suitable for rough wall boundary layers is used to blend between the laminar and fully turbulent state. Finally, two different submodels for the fully turbulent boundary layer complete the two models. In the first model, termed KS-TLK-T in this paper, a sand roughness approach from Durbin et al., which builds upon a two-layer k-ε-turbulence model, is used for this purpose. The second model, the so-called DEM-TLV-T model, makes use of the discrete-element roughness approach, which was recently combined with a two-layer k-ε-turbulence model by the present authors. The discrete element model will be formulated in a new way suitable for randomly rough topographies. Part I of the paper will provide detailed model formulations as well as a description of the database used for developing the new transition onset correlation. Part II contains a comprehensive validation of the two models, using a variety of test cases with transitional and fully turbulent boundary layers. The validation focuses on heat transfer calculations on both, the suction and the pressure side of modern turbine airfoils. Test cases include extensive experimental investigations on a high-pressure turbine vane with varying surface roughness and turbulence intensity, recently published by the current authors as well as new experimental data from a low-pressure turbine vane. In the majority of cases, the predictions from both models are in good agreement with the experimental data.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Wei-Jen Huang ◽  
Ming-Ta Lee ◽  
Kuei-Chen Huang ◽  
Kai-Jung Kao ◽  
Ming-An Lee ◽  
...  

AbstractThe release of anthropogenic radiocesium to the North Pacific Ocean (NPO) has occurred in the past 60 years. Factors controlling 137Cs (half-life, 30.2 year) and 134Cs (half-life, 2.06 year) activity concentrations in the Kuroshio east of Taiwan and the Taiwan Strait (latitude 20° N–27° N, longitude 116° E–123° E) remain unclear. This study collected seawater samples throughout this region and analyzed 134Cs and 137Cs activity concentrations between 2018 and 2019. A principal component analysis (PCA) was performed to analyze the controlling factors of radiocesium. Results of all 134Cs activity concentrations were below the detection limit (0.5 Bq m−3). Analyses of water column 137Cs profiles revealed a primary concentration peak (2.1–2.2 Bq m−3) at a depth range of 200–400 m (potential density σθ: 25.3 to 26.1 kg m−3). The PCA result suggests that this primary peak was related to density layers in the water column. A secondary 137Cs peak (1.90 Bq m−3) was observed in the near-surface waters (σθ = 18.8 to 21.4 kg m−3) and was possibly related to upwelling and river-to-sea mixing on the shelf. In the Taiwan Strait, 137Cs activity concentrations in the near-surface waters were higher in the summer than in the winter. We suggest that upwelling facilitates the vertical transport of 137Cs at the shelf break of the western NPO.


Sign in / Sign up

Export Citation Format

Share Document