风生紊流导致微囊藻群体破碎和形态变化

微囊藻群体大小和形态决定其垂向迁移能力,从而影响着水华的形成.为了探讨湖泊中风生紊流对微囊藻群体大小和形态的影响,本研究于2012年8月26日至9月7日在太湖梅梁湾的围隔内进行了12 d的昼夜不间断的高频采样（采样间隔每2小时一次）.研究期间,水面微囊藻密度呈现4次周期性消涨,藻密度变化范围为4×10⁴~2671×10⁴ cells/mL.而整个水柱中的藻密度变化范围仅为3×10⁴~18×10⁴ cells/mL.皮尔逊相关性分析表明微囊藻的原位生长速率与表面藻密度呈负相关而与风速呈正相关.强风速使微囊藻在水柱中均匀分散,增强了透光性,促进了微囊藻的生长.微囊藻群体粒径随着风速的增大逐渐减小,反之亦然.其中值粒径（D₅₀）变化范围为66.2~768.0 μm.在此期间微囊藻群体形态主要以鱼害微囊藻、不规则的惠氏微囊藻、球状的惠氏微囊藻和铜绿微囊藻群体形态为主,其占比也呈现出波动状态.皮尔逊相关分析结果显示微囊藻群体大小与风速呈负相关,说明湖泊中风生紊流会影响微囊藻群体大小.当紊流强度为2.33×10^-5 m²/s³时,微囊藻群体会发生破碎现象,该紊流强度相当于5 m/s的风在30 m深的水库或湖泊中所产生的紊流强度.微囊藻群体被风生紊流破碎后最大粒径与该风速下紊流的最小涡旋尺度相近,表明紊流的最小涡旋尺度决定了微囊藻所能形成群体的最终大小.监测期间,整水柱中不同群体形态的微囊藻占比发生了明显变化,在监测初期以鱼害微囊藻群体形态为主,随后不规则的惠氏微囊藻和铜绿微囊藻群体形态的比例不断增加,最后鱼害微囊藻群体形态又占据主导地位.球状的惠氏微囊藻群体形态在整个监测期中的比例随时间的增加而逐渐降低.不同群体形态微囊藻之间比例的大幅变化无法用微囊藻生长演替来解释.而皮尔逊相关分析结果显示鱼害微囊藻与惠氏微囊藻（不规则的和球状的惠氏微囊藻之和）群体形态之间存在负相关,且惠氏微囊藻与铜绿微囊藻群体形态呈负相关.但在今后研究中需进一步关注在微囊藻群体形态的动态变化过程中细胞大小、胶被、产毒特性和基因序列等特征,从而验证不同种微囊藻群体是否存在形态转换这一猜想.总而言之,普通强度的风生紊流能够破碎微囊藻群体,而气候变化导致的内陆湖泊周边风速下降会促使微囊藻形成更大的群体,从而有利于水华的形成.     

Modeling the effect of saltation on surface layer turbulence

Particle–turbulence interaction has been a research focus in the field of pneumatic transport, especially in aeolian environments. However, knowledge regarding the effect of saltating particles on the turbulence characteristics is very limited. In this article, a process of sand-laden flow from forming sand streamers to stability is investigated via a coupled mathematical model of wind-blown sand that includes the spatiotemporal development. The variations in the turbulence characteristics, such as the mean velocity and turbulence intensity in clean air or sand-laden flow field, are analyzed. The results show that the splash process of sand grains near the wall decrease the wind speed in the saltation layer and destroy the low-speed streaks. Moreover, the profiles of streamwise turbulence intensity exhibit a transition from ‘decreasing’ to ‘increasing’ and approximately intersect at an ‘intensity focus’, which is presented for the first time. Furthermore, it was found that saltating particles could enhance the Reynolds stress. Meanwhile, it was also noticed that the shear stress at the wall surface is greater than the impact threshold and that there is a tendency towards the impact threshold. Therefore, saltation makes the particle Reynolds number of sand-laden flow higher than that under non-saltation conditions, thus changing the particles’ effect on the turbulence intensity. Gravity-dominated saltation is probably the most essential difference between wind-blown sand and other traditional two-phase flows. © 2020 John Wiley & Sons, Ltd.     

Broken ergodicity in magnetohydrodynamic turbulence

Turbulent magnetofluids appear in various geophysical and astrophysical contexts, in phenomena associated with planets, stars, galaxies and the universe itself. In many cases, large-scale magnetic fields are observed, though a better knowledge of magnetofluid turbulence is needed to more fully understand the dynamo processes that produce them. One approach is to develop the statistical mechanics of ideal (i.e. non-dissipative), incompressible, homogeneous magnetohydrodynamic (MHD) turbulence, known as “absolute equilibrium ensemble” theory, as far as possible by studying model systems with the goal of finding those aspects that survive the introduction of viscosity and resistivity. Here, we review the progress that has been made in this direction. We examine both three-dimensional (3-D) and two-dimensional (2-D) model systems based on discrete Fourier representations. The basic equations are those of incompressible MHD and may include the effects of rotation and/or a mean magnetic field B _o. Statistical predictions are that Fourier coefficients of the velocity and magnetic field are zero-mean random variables. However, this is not the case, in general, for we observe non-ergodic behavior in very long time computer simulations of ideal turbulence: low wavenumber Fourier modes that have relatively large means and small standard deviations, i.e. coherent structure. In particular, ergodicity appears strongly broken when B _o?=?0 and weakly broken when B _o?≠?0. Broken ergodicity in MHD turbulence is explained by an eigenanalysis of modal covariance matrices. This produces a set of modal eigenvalues inversely proportional to the expected energy of their associated eigenvariables. A large disparity in eigenvalues within the same mode (identified by wavevector k ) can occur at low values of wavenumber k?=?| k |, especially when B _o?=?0. This disparity breaks the ergodicity of eigenvariables with smallest eigenvalues (largest energies). This leads to coherent structure in models of ideal homogeneous MHD turbulence, which can occur at lowest values of wavenumber k for 3-D cases, and at either lowest or highest k for ideal 2-D magnetofluids. These ideal results appear relevant for unforced, decaying MHD turbulence, so that broken ergodicity effects in MHD turbulence survive dissipation. In comparison, we will also examine ideal hydrodynamic (HD) turbulence, which, in the 3-D case, will be seen to differ fundamentally from ideal MHD turbulence in that coherent structure due to broken ergodicity can only occur at maximum k in numerical simulations. However, a nonzero viscosity eliminates this ideal 3-D HD structure, so that unforced, decaying 3-D HD turbulence is expected to be ergodic. In summary, broken ergodicity in MHD turbulence leads to energetic, large-scale, quasistationary magnetic fields (coherent structures) in numerical models of bounded, turbulent magnetofluids. Thus, broken ergodicity provides a large-scale dynamo mechanism within computer models of homogeneous MHD turbulence. These results may help us to better understand the origin of global magnetic fields in astrophysical and geophysical objects.     

A one-dimensional turbulence model; Enstrophy cascades and small-scale intermittency in decaying turbulence

Abstract

Severe unidirectional Fourier truncation of the equations for 2-D incompressible flow leads to a system of three coupled PDEs in one space dimension with the same quadratic invariants as the original set (i.e. energy and enstrophy). Numerically generated equilibria for inviscid, truncated versions of the reduced system are well approximated by Kraichnan's energy-enstrophy equipartition spectra. Viscous calculations for decaying turbulence at moderate resolution (1024 degrees of freedom) also appear to be consistent with a direct, k ^?3, enstrophy cascading inertial range when the dissipation is small. Dissipation range intermittency in the form of spatially intermittent enstrophy dissipation with occasional strong bursts producing linear phase locking is also observed. In contrast to full 2-D simulations, no tendency towards the emergence of isolated, coherent vorticity structures is observed. The model consequently mimics some, but not all, of the properties of the full 2-D set.     

高层结构群风压力分布与干扰效应

为研究群风环境中高层结构的风压分布特性与群风干扰效应,采用1:300缩尺刚性模型风洞试验,在24种风向角下分析了高层结构风荷载特性的群风效应,研究了结构的压分布特性与群风干扰效应。结果表明,现有《建筑结构荷载规范》的结构风压设计值不宜用于群风环境中高层结构抗风设计,而模型风洞试验可有效地反映群风环境中的结构风压特性;迎风面与背风面脉动压力系数较侧风面低,最大脉动压力出现在涡旋运动最强烈处,也是最大负压产生处;结构极值风压可以通过脉动风压的概率特性识别,进而对结构采取相应的加强措施。试验研究技术可靠,结果可用于该工程及类似工程的抗风设计参考。     

Mixing efficiency in stably-stratified decaying turbulence

Abstract

The behavior of the flux Richardson number R _f, as a function of the overall Richardson number Ri ₀, was investigated for a stably stratified, grid-generated, turbulent flow evolving in a closed-loop water channel. The turbulent dissipation rate ε, the buoyancy or vertical mass flux p wbar; and the rms density fluctuation ρ′ were obtained from simultaneous single-point measurements of the horizontal and vertical velocity components and density fluctuations. From these, R _f and Ri ₀ were calculated at each point in the spatially evolving flow. The resulting curves of R _f vs. Ri ₀ exhibit the full range of behavior found in the very different case studied by Linden (1980). The length scale arguments of Gibson (1980) and Stillinger et al. (1983b) provide an underlying mechanism which successfully accounts for the shape of the R _f vs. Ri ₀ curve.     

Local turbulence in the Earth's core

Abstract

It is shown that in the Earth's core, where the geodynamo is at work (and is supplied with energy by the prevailing unstable density stratification), a buoyancy instability of a local character exists which is highly supercritical. This instability results in fully developed turbulence dominated by small scale vortices. The influence of the Earth's rotation and of the magnetic field produced by the geodynamo makes this small scale turbulence highly anisotropic. A qualitative picture of this local anisotropic turbulence is devised and the main parameters characterizing it are estimated. Expressions for the turbulent diffusivity are developed and discussed.     

Structure of triads in geostrophic turbulence

The interaction between Fourier waves, which results in energy transfer over the spectrum, has been considered using the Boussinesq model in a plane layer during rapid rotation as an example. It has been indicated that the wave triangle spectrum strongly differs on the scales shorter (longer) than the scale of a cyclonic turbulence leading mode.     

Modelling of convective turbulence with a two-equation k-ɛ turbulence closure scheme

 In this study we aim at comparing turbulence parameters from field observations and model simulations under convective conditions. The comparison is focused on the depth dependence and temporal dynamics of viscous and diffusive dissipation rates ɛ and χ. The near-surface observations were obtained by using a quasi-free-rising profiler which measured small-scale shear and temperature fluctuations to within the vicinity of the water surface. Convective conditions during the experiment are characterized by low wind speeds (between 0 and 4 m s⁻¹) and a typical heat loss of about 150 Wm⁻². We applied a state-of-the-art two-equation k-ɛ turbulence model with an algebraic second-moment closure scheme. The qualitative agreement of the turbulence quantities resulting from observations and simulations is rather good. The temporal dynamics of the temperature field is simulated correctly, whereas in the spatial dynamics some deficiencies of the model due to its local character can be seen. It is concluded that such models realistically reproduce convective turbulence and therefore represent a reasonable compromise between complexity and simplicity, so that they can be used with acceptable costs in large-scale models. Received: 31 October 2001 / Accepted: 20 September 2002 Acknowledgements The whole project was intitiated by Peter Schlittenhardt, who strongly supported the development of the observational technique in uprising mode and encouraged us to undertake the experiment in Lake Maggiore; for which we will thank him most. Performing the measurement campaign was only possible with the help of several colleagues from the Marine Environment Unit at the Joint Research Centre. Thanks to all of them, but especially to Dirk van d. Linde, Ulisse Devisioni, Bjarke Rasmussen and Hartmut Prandke. The Istituto Italiano di Idrobiologia, Verbania Pallanza, provided the boat Livia for the installation and removal of the measurement system; thanks to the master Stefano Maurizio for his engagement. Ute Tschesche provided the data evaluation and presentation software and Judith Challis helped to polish our English. Part of the study was␣supported by the PROVESS project (MAS3-CT97-0159); thanks to our sponsors at the European Commission. We are further grateful to two anonymous referees for their constructive comments. Responsible Editor: Charitha Pattiaratchi     

The effect of wave-induced turbulence on intertidal mudflats: Impact of boat traffic and wind

Semi-diurnal and fortnightly surveys were carried out to quantify the effects of wind- and navigation-induced high-energy events on bed sediments above intertidal mudflats. The mudflats are located in the upper fluvial part (Oissel mudflat) and at the mouth (Vasière Nord mudflat) of the macrotidal Seine estuary. Instantaneous flow velocities and mudflat bed elevation were measured at a high frequency and high resolution with an acoustic doppler velocimeter (ADV) and an ALTUS altimeter, respectively. Suspended particulate matter concentrations were estimated by calibrating the ADV acoustic backscattered intensity with bed sediments collected at the study sites. Turbulent bed shear stress values were estimated by the turbulent kinetic energy method, using velocity variances filtered from the wave contribution. Wave shear stress and maximum wave–current shear stress values were calculated with the wave–current interaction (WCI) model, which is based on the bed roughness length, wave orbital velocities and the wave period (T_S). In the fluvial part of the estuary, boat passages occurred unevenly during the surveys and were characterized by long waves (T_S>50 s) induced by the drawdown effect and by short boat-waves (T_S<10 s). Boat waves generated large bottom shear stress values of 0.5 N m⁻² for 2–5 min periods and, in burst of several seconds, larger bottom shear stress values up to 1 N m⁻². At the mouth of the estuary, west south-west wind events generated short waves (T_S<10 s) of H_S values ranging from 0.1 to 0.3 m. In shallow-water environment (water depth <1.5 m), these waves produced bottom shear stress values between 1 and 2 N m⁻². Wave–current shear stress values are one order of magnitude larger than the current-induced shear stress and indicate that navigation and wind are the dominant hydrodynamic forcing parameters above the two mudflats. Bed elevation and SPM concentration time series showed that these high energy events induced erosion processes of up to several centimetres. Critical erosion shear stress (τ_ce) values were determined from the SPM concentration and bed elevation measurements. Rough τ_ce values were found above 0.2 N m⁻² for the Oissel mudflat and about 1 N m⁻² for the Vasière Nord mudflat.     

Sand settling through bedform-generated turbulence in rivers

Fluvial bedforms generate a turbulent wake that can impact suspended-sediment settling in the passing flow. This impact has implications for local suspended-sediment transport, bedform stability, and channel evolution; however, it is typically not well-considered in geomorphologic models. Our study uses a three-dimensional OpenFOAM hydrodynamic and particle-tracking model to investigate how turbulence generated from bedforms and the channel bed influences medium sand-sized particle settling, in terms of the distribution of suspended particles within the flow field and particle-settling velocities. The model resolved the effect of an engineered bedform, which altered the flow field in a manner similar to a natural dune. The modelling scenarios alternated bed morphology and the simulation of turbulence, using detached eddy simulation (DES), to differentiate the influence of bedform-generated turbulence relative to that of turbulence generated from the channel bed. The bedform generated a turbulent wake that was composed of eddies with significant anisotropic properties. The eddies and, to a lesser degree, turbulence arising from velocity shear at the bed substantially reduced settling velocities relative to the settling velocities predicted in the absence of turbulence. The eddies tended to advect sediment particles in their primary direction, diffuse particles throughout the flow column, and reduced settling likely due to production of a positively skewed vertical-velocity fluctuation distribution. Study results suggest that the bedform wake has a significant impact on particle-settling behaviour (up to a 50% reduction in settling velocity) at a scale capable of modulating local suspended transport rates and bedform dynamics. © 2020 John Wiley & Sons, Ltd.     

Topographically generated steady currents in barotropic turbulence

Abstract

Steady currents develop in oceanic turbulence above topography even in the absence of steady forcing. Mesoscale steady currents are correlated with mesoscale topography with anticyclonic eddies above topographic bumps, and large scale westward flows develop when β is non-zero. The relationship between those two kinds of steady currents, as well as their dependence on various parameters, is studied using a barotropic quasi-geostrophic channel model. The percentage of steady energy is found to depend on the forcing, friction and topography in a non-monotonic fashion. For example, the percentage of steady currents grows with the energy level in the linear regime (low energies) and decreases when the energy level increases in the nonlinear regime (high energies). Mesoscale steady currents are the energy source for the steady westward flow U, and therefore U is the maximum when large scale and mesoscale currents are of the same order of magnitude. This happens when the ratio S of the large scale slope βH/f ₀ and the mesoscale rms topographic slope α is of order one. U decreases for both small and large values of S.     

Wavelet analysis of turbulence in cirrus clouds

Two flights of the UK Meteorological Offices Hercules aircraft through daytime frontal cirrus around Scotland have been analysed using wavelet analysis on the vertical velocity time-series from the horizontal runs. It is shown that wavelet analysis is a useful tool for analysing the turbulence data in cirrus clouds. It finds the largest scales involved in producing turbulence, as does Fourier analysis, such as the 2-km spectral peaks corresponding to convective activity during flight A283. Wavelet spectra have the added advantage that the position is shown, and so they identify smaller-scale, highly localised processes such as the production of turbulent kinetic energy by the breaking of Kelvin-Helmholtz waves due to the vertical shear in the horizontal wind. These may be lost in Fourier spectra obtained for long time-series, though they contribute something to the average spectral density at the appropriate scale. The main disadvantage of this technique is that only octave frequency bands are resolved.     

On the influence of microphysics parametrization on the rainfall rates in numerical models of clouds

Accuracy of formulas for growth by accretion and evaporation of rain in bulk parametrization of these processes for the case of light and moderate precipitation is investigated. It is done by comparison of results from two simple models: with bulk approach and with exact calculations of growth or evaporation of drops in each size bin separately. Growth by accretion is accurately represented in bulk parametrization but rain evaporation is overpredicted. Corrected formula for rain evaporation is suggested.     

Comparing parametrization of the coefficient of molecular thermal conductivity of the thermosphere

At present, different parametrizations of the molecular thermal conductivity coefficient are used in the models of the dynamics and neutral and ion composition of the thermosphere. These parametrizations have been compared by solving the energy balance equation for medium-global conditions. The spread in estimates reaches ～300 K for temperature of the thermosphere and ～40% for density at an altitude of 300 km. A new parametrization of the molecular thermal conductivity coefficient has been proposed. This parametrization is based on the renewed temperature dependences of molecular conductivity [Uribe et al., 1990] according to measurements in pure N₂ and O₂ gases and on the last theoretical estimates of O atom elastic collision sections [Kharchenko et al., 2000].     

On the large-scale atmospheric turbulence

Summary The dynamical equations of the velocity correlation and spectrum functions for three-dimensional anisotropic turbulence in a stratified, rotating atmosphere are derived and discussed. The characteristics of the horizontal two-dimensional isotropic turbulence are analyzed.National Center for Atmospheric Research, sponsored by the National Science Foundation. On leave from the University of Utah, Salt LakeCity, Utah.     

On Wind Profiles over Developing Wind Wave in the Water-Air Interface

The results of studying the vertical profile of wind velocity and the vortex-formation mechanism in the atmospheric layer immediately adjacent to the water surface under the condition of developing wind wave are considered. Materials of field and laboratory experiments are used.__________Translated from Vodnye Resursy, Vol. 32, No. 3, 2005, pp. 295–300.Original Russian Text Copyright © 2005 by Anisimova, Pokazeev, Soboleva, A. Speranskaya, O. Speranskaya.     

High-resolution numerical simulation of turbulence in natural waterways

We develop an efficient and versatile numerical model for carrying out high-resolution simulations of turbulent flows in natural meandering streams with arbitrarily complex bathymetry. The numerical model solves the 3D, unsteady, incompressible Navier-Stokes and continuity equations in generalized curvilinear coordinates. The method can handle the arbitrary geometrical complexity of natural streams using the sharp-interface curvilinear immersed boundary (CURVIB) method of Ge and Sotiropoulos (2007) [1]. The governing equations are discretized with three-point, central, second-order accurate finite-difference formulas and integrated in time using an efficient, second-order accurate fractional step method. To enable efficient simulations on grids with tens of millions of grid nodes in long and shallow domains typical of natural streams, the algebraic multigrid (AMG) method is used to solve the Poisson equation for the pressure coupled with a matrix-free Krylov solver for the momentum equations. Depending on the desired level of resolution and available computational resources, the numerical model can either simulate, via direct numerical simulation (DNS), large-eddy simulation (LES), or unsteady Reynolds-averaged Navier-Stokes (URANS) modeling. The potential of the model as a powerful tool for simulating energetic coherent structures in turbulent flows in natural river reaches is demonstrated by applying it to carry out LES and URANS in a 50-m long natural meandering stream at resolution sufficiently fine to capture vortex shedding from centimeter-scale roughness elements on the bed. The accuracy of the simulations is demonstrated by comparisons with experimental data and the relative performance of the LES and URANS models is also discussed.     

Large-scale vortices in helical turbulence of incompressible fluid

Abstract

The interaction of a mean flow with a random fluctuation field is considered. This interaction is described by the averaged Navier-Stokes equation in which terms nonlinear in the fluctuation field are expressed in terms of the mean flow and the statistical properties of the fluctuation field, which is assumed to be homogeneous, isotropic, and helical. Averaged equations are derived using a functional technique. These equations are solved for a mean background flow that depends linearly on the position vector. The solutions show that large-scale vortices may arise in this system.