深海垂向流速剖面数值模拟

运用CFD软件对深海各种不同的垂向流速剖面进行数值模拟,分析研究试验区域内水流在空间上的流速分布情况以及湍流强度的变化,并与目标流速进行了比较.结合模型试验,表明数值模拟的结果是可靠的,造流系统能够较为准确地模拟出不同垂向流速剖面的海流,满足深水试验池对深海海流的模拟要求.     

Generation of internal waves by moving non-stationary disturbances in a real stratified ocean

The plane problem on the generation of linear internal waves by a moving area of time-harmonic surface pressures in a continuously-stratified ocean of constant depth is considered. An analytical relation has been derived for forced internal waves off the site of their generation in the form of an internal wave field superposition corresponding to individual vertical modes. The possible wave regimes are determined. For the Brunt-Väisälä frequency distribution in the North Atlantic, the generation conditions and amplitudes of diverse radiated waves are numerically determined.Translated by Vladimir A. Puchkin.     

Turbulence in the ocean: Development of investigations, evolution of ideas, and prospects

We present an analytic survey of a series of works devoted to the direct methods for the investigation of turbulence in natural basins originated under the guidance of A. G. Kolesnikov in the 50s. We describe the principal results of the early stage which formed the basis for future investigations. The main scientific directions of the latest period are also considered from the viewpoint of the development of Kolesnikov's ideas: investigation of the processes of generation and energy supply of small-scale turbulence, relationship between turbulence and the fine vertical structure and internal waves, and the mechanisms of vertical turbulent exchange in the ocean. A model of vertical exchange under the conditions of intermittent “source” turbulence is distinguished. We also discuss the directions of development and prospects of the investigation of turbulence. Translated by Peter V. Malyshev and Dmitry V. Malyshev     

Generation of long internal waves in the ocean by a moving pressure zone

The generation of plane internal waves in a continuously stratified ocean by a steadily moving anomaly of the atmospheric pressure is studied in a linear statement using a quasi-static approximation. Quantitative estimates for the parameters of internal waves in the wake are obtained for the average distribution of the Väisälä-Brunt frequency in the Caribbean Sea.Translated by Mikhail M. Trufanov.     

The energy cycle of wind waves on the sea surface

The problems of wind-induced waves on the sea surface are considered. To this end, the empirical fetch laws that determine variations in the basic periods and heights of waves in relation to their fetch are used. The relation between the fetch and the physical time is found, as are the dependences of the basic characteristics of waves on the time of wind forcing. It is found that about 5% of wind energy dissipated in the near-water air layer contributes to the growth of wave heights, i.e. wave energy, although this quantity depends on the age of waves and the exponent in the fetch laws. With consideration for estimates of the probability distribution functions for the wind over the world ocean [11], it is found that the rate of wind-energy dissipation in the near-water air layer is on the order of 1 W/m². The calculations of wind waves [19] for the world ocean for 2007 have made it possible to assess the mean characteristics of the cycle of wave development and their seasonal variations. An analysis of these calculations [19] shows that about 20% of wind energy is transferred to the water surface. The remaining amount (80%) of wind energy is spent on the generation of turbulence in the near-water air layer. About 2%, i.e., one tenth of the energy transferred to water, is spent on turbulence generation due to the instability of the vertical velocity profile of the Stokes drift current and on energy dissipation in the surf zones. Of the remaining 18%, 5% is spent directly on wave growth and 13% is spent on the generation of turbulence during wave breaking and on a small-scale spectral region. These annually and globally mean estimates have a seasonal cycle with an amplitude on the order of 20% in absolute values but with a smaller amplitude in relative values. According to [19] and to the results of this study, the annually mean height of waves is estimated as 2.7 m and their age is estimated as 1.17.     

Experimental analysis of the flow field around horizontal axis tidal turbines by use of scale mesh disk rotor simulators

Understanding the flow field around horizontal axis marine current turbines is important if this new energy generation technology is to advance. The aim of this work is to identify and provide an understanding of the principal parameters that govern the downstream wake structure and its recovery to the free-stream velocity profile. This will allow large farms or arrays of devices to be installed whilst maximising device and array efficiency. Wake characteristics of small-scale mesh disk rotor simulators have been measured in a 21 m tilting flume at the University of Southampton. The results indicate that wake velocities are reduced in the near wake region (close behind the rotor disk) for increasing levels of disk thrust. Further downstream all normalised wake velocity values converge, enforcing that, as for wind turbines, far wake recovery is a function of the ambient flow turbulence. Varying the disk proximity to the water surface/bed introduces differential mass flow rates above and below the rotor disk that can cause the wake to persist much further downstream. Finally, the introduction of increased sea bed roughness whilst increasing the depth-averaged ambient turbulence actually decreases downstream wake velocities. Results presented demonstrate that there are a number of interdependent variables that affect the rate of wake recovery and will have a significant impact on the spacing of marine current turbines within an array.     

Investigation of the sensitivity of the Mellor–Yamada parametrization to the choice of finite-difference analogs in a numerical three-dimensional model of the operative prediction of currents in the black sea

Within the framework of the Mellor–Yamada approach, we realize a numerical scheme for the calculation of the coefficients of turbulent viscosity and diffusion in the z -system of coordinates for the three-dimensional model of operative prediction of currents in the Black Sea. Some discrete analogs of the equations for turbulent kinetic energy and turbulence macroscale are studied. Their high sensitivity to the choice of finite-difference approximations is demonstrated. On the basis of the comparison of the results of prognostic experiments with the data of observations, we choose the best approximation of the term used to describe the generation pf turbulence energy.     

长江口湍流剖面的观测与分析

小尺度湍流过程对河口物质输运与能量交换至关重要。受传统观测方法的限制, 河口浅水区域的剖面观测资料至今较为匮乏, 进而限制了湍流过程的研究。为此, 采用新型5波束声学多普勒流速剖面仪(Nortek Signature 1 000 kHz AD2CP)在长江口开展湍流剖面观测, 获取高频、低噪的高质量湍流剖面数据, 并与声学多普勒点式流速仪(acoustic doppler velocimeters, ADV)同步观测的数据进行对比。结果表明, 通过AD2CP与ADV获得的近底部边界层摩阻流速u_*、拖曳系数C_d、雷诺应力S_R等特征参数基本一致, 底摩擦与波浪能量为河口区域湍动能的主要输入源。湍流垂向结构存在显著的非局地平衡, 即温盐等斜压作用引起的浮力通量、对流项以及强波浪作用影响的脉动压力做功、黏性输运等因素可能为长江口湍流非局地平衡的主因。     

Turbulence and zooplankton production: insights from PROVESS

Zooplankton are directly influenced by turbulence in both a passive and an active manner. Passively, zooplankton are at the mercy of turbulence in how it affects their vertical mixing, encounter rate, detection abilities and feeding current efficiency. Many zooplankton species, however, are actively able to mitigate the effects of turbulence by modifying their behaviour, e.g. vertical migration, prey switching and habituation to hydromechanical stimuli. Both theoretical treatments of these processes and field observations from the northern North Sea are examined. Field observations show that some copepod species actively migrate to avoid high turbulence levels in surface waters. Furthermore, observations show a negative relationship between turbulence and zooplankton ingestion rates. This supports the paradigm of a dome-shaped response for zooplankton production with environmental turbulence. A theoretical treatment shows that the reaction distance, R, for an ambush-feeding copepod feeding on swimming organisms follows R ∝ ^−1/6 where is the turbulent dissipation rate, a result that shows close agreement with previously reported experimental results.     

Some properties of tidal currents estimated from analytical and LES simulation studies

Some features of tidal current ellipses near the ocean floor are estimated from the bottom-limit solutions of the tidal equations with constant eddy viscosity. It is verified that a clockwise rotating ellipse becomes broad descending through the bottom boundary, whereas a counterclockwise ellipse becomes narrow except near the critical latitude. The major axis of an ellipse near the floor is directed 45° leftward (rightward) against that of the overlying tidal flow ellipse poleward from the critical latitude in the northern (southern) hemisphere. Equatorward from the critical latitude, the major axis near the floor lines up with that of the overlying tidal current ellipse in both hemispheres. The retrograde rotating hodograph against the overlying clockwise tidal hodograph appears near the floor around the critical latitude. These features are confirmed from the numerical studies using the large eddy simulation model. The deflection of the major axis poleward from the critical latitude is a little smaller, about 19°. Furthermore, the effects of nontraditional Coriolis force resulting from the horizontal component of the earth-rotating vector on the development of turbulence were investigated. When the diurnal tidal current is directed westward, the effective vertical eddy viscosity becomes maximum with developing turbulence. While in the case of a semidiurnal tide, a northward tidal current is preferred for turbulence development.     

Turbulence modeling in ocean circulation problems

A physical formulation of the problem is considered. A mathematical model and a numerical algorithm of the turbulence model as part of the ocean circulation model for simulations for decades are formulated. The model is based on the evolution equations for turbulent kinetic energy (TKE) and the frequency of its viscous dissipation. A numerical solution algorithm for both the circulation model and the turbulence model is based on implicit schemes of splitting with respect to physical processes and geometric coordinates. For the turbulence model, this provided analytical solutions at a splitting step related to TKE generation and dissipation. Numerical experiments have been performed with a model of the joint circulation of the North Atlantic, the Arctic Ocean, and the Bering Sea to reproduce the annual cycle and synoptic disturbances of ocean characteristics. The model has a resolution of 0.25° in latitude and longitude and 40 levels in the vertical, which are compressed toward the surface to reproduce the process of developed turbulence better. The results are compared with observations and with the results of simulations using traditional parameterizations of the upper ocean mixing. It is shown that the model reproduces ocean characteristics correctly, only slightly increasing the computation time in comparison with simple parameterizations. Spatial and temporal characteristics of turbulence are analyzed.     

A method of submersion in the problem on the determination of dispersion curves of internal waves in the presence of a current with vertical velocity shear

The previously proposed method [1] of submerging aimed at determining internal wave dispersion characteristics is generalized to consider the case of a background horizontal flow with a vertical velocity shear. The results of calculating, dispersion curves of discrete spectrum modes for the profiles of the Brunt-Väisälä frequency and the mean current velocity are given; they agree well with the results of studies carried out in ref. [5].Translated by Mikhail M. Trufanov.     

黄河三角洲堤前泥沙起动现场观测研究

使用波潮仪、流速计以及浊度计现场测量记录海堤前的波高、潮高、流速和浊度值,与室内测定的浊度和水中悬浮颗粒物含量对应值比较,分析计算浊度与现场测量波高、潮高、流速水动力因素间的相关性,得出堤前泥沙起动的决定性水动力因素。在悬移质运移为主情况下,由于海堤的阻挡,波浪破碎产生的旋涡和紊动是泥沙起动的主要动力因素,泥沙起动量与波高值呈现较高相关性;海堤前复杂水动力条件中,波浪破碎前的流速值与悬移质泥沙起动量相关系数较小,启动流速的主要组成部分是破波产生的旋涡和紊动的流速。     

潮流波浪联合输沙及海床冲淤演变的理论体系与其数学模拟

根据近岸带及河口区潮流、波浪、湍流各自物理尺度的不同,从Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程和质量传输方程出发,利用Ｒｅｙｎｏｌｄｓ分解的方法,建立了模拟波浪 流联合输沙及海床冲淤演变的理论体系,给出了潮流作用下近岸波浪传播方程、波浪作用下潮流运动方程并通过利用波流合成底部切应力、底层湍流脉动随机特性,得出了波流联合作用下不平衡沙计算中泥沙起悬与沉降量的确定方法。本文模型应用于“广西合浦围垦工程潮流波浪     

珠江河口底边界层憩流动力特征

通过在珠江黄茅海河口进行的25 h定点连续观测,分析了潮流底边界层在憩流时刻的动力特征。结果表明,憩流期非恒定性增强,改变了平均流和湍流结构。主要表现:1)平均流结构出现三种不同类型,分别是对数分布(摩擦力主导)、垂向均匀分布(惯性主导)和过渡状态(摩擦力和惯性力平衡);2)湍流各向异性增强,湍谱惯性子区带宽减小。各湍流特征参数在转流时刻达到最小值,表明惯性作用抑制了湍流的发育。     

Laboratory study of the dynamics of a layer of turbulence excited by standing surface waves

The evolution of a turbulent layer excited by standing surface parametric waves is studied experimentally. The time dependence of the depth of the turbulized region is obtained. It is found that the thikness of the layer is described by a nonlinear dependence as distinguished from the previously obtained theoretical results. The dynamics of turbulent-layer embedding is described on the basis of a semiempirical theory and the assumption that the outer scale of turbulence is constant. Satisfactory agreement with experimental data is obtained. In the context of the proposed model, a relationship is found between the quantities characterizing the turbulence and surface waves. The dependence of the steepness of surface waves on the acceleration of the bottom is obtained experimentally. The data obtained in this study can be used for estimation of turbulence generation during a submarine earthquake.     

悬对沙Mellor-Yamada海洋边界层模型计算结果的影响

为了分析海洋水体垂向水流紊动及紊动交换情况而采用了一维的海洋边界层模型(Mellor-Yamada)并利用数值实验的方法对悬沙、盐度、温度等数据进行分析。原模型未将悬沙考虑在内，本文试将它引入进去探讨由于它的存在对紊动混合特性的影响。2000年4月，Mellor将最初的模型引入了依赖于Richardson数的紊动动能耗散率。本文通过比较具有悬沙和不具有悬沙两种情况下的速度、温度和盐度垂向分布随时间的变化，分析讨论由于悬沙的存在所引发的密度层化对紊动混合作用的影响，并发现悬浮泥沙抑制了部分模拟时间的紊动混合作用。     

Second-order wavemaker theory for multidirectional waves

The present paper develops the complete second-order wavemaker theory for the generation of multidirectional waves in a semi-infinite basin. The theory includes superharmonics and subharmonics and is valid for a rotational as well as a translatory serpent-type wave-board motion. The primary goal is to obtain the second-order motion of the wave paddles required to get a prescribed multidirectional irregular wave field correct to second order, i.e. to suppress spurious free-wave generation. The wavemaker theory is a 3D extension of the full second-order wavemaker theory for wave flumes by Schäffer (1996).