海岸沿岸流垂向分布实验研究

进行海岸沿岸流模型实验,利用ADV(acoustic doppler velocimeter)测量了沿岸流流速沿垂直岸线方向多个断面的垂向分布。采用Faria等的方法对实验结果进行了分析,将沿岸流沿水深分为上层和下层两部分:下层采用对数分布来表达沿岸流的垂向分布;上层考虑了波浪引起的自由表面的波动,流速分布为修正的对数分布。最后对沿岸流垂向分布特征做了分析,并与Visser和Hamilton的实验结果做了对比。     

根据流速剖面估计海底粗糙长度的研究

由于底床的摩擦作用,河口海岸近底层的水体流动与远离海底上层的水体流动产生很大的差别,形成重要的底边界层。在近底层由于水层与床面有频繁的物质交换,所以研究河口海岸边界层对水流、泥沙、温度、盐度垂向混合、底沙运动等物质过程有着重要的作用。底床粗糙度是研究河口海岸沉积物运移和水流结构的重要因素,它由颗粒粗糙度、底形粗糙度和推移质输沙粗糙度组成。目前还没有一个经验公式来计算底床粗糙度,但是在无波浪的情况下粗糙度可以通过拟合流速的对数分布(卡门-普朗特公式)来求得。利用对数拟合方法必须满足两个条件:一是必须测量近海底三个层面以上的流速,二是测得的流速剖面必须满足对数分布。     

不同湍流模式下钱塘江涌潮水流三维模拟

钱塘江涌潮具有动力作用强和流速变化快等特点。涌潮水流紊动复杂,流速的垂向分布和紊动强度息息相关。通过涌潮水流实测资料的分析可以发现,涌潮作用下流速垂向分布在底部和上层存在差异。为研究涌潮作用下流速垂向分布的特征,应用基于非结构网格下有限体积法模型FVCOM对钱塘江涌潮河段水流运动进行三维数值模拟。考虑到涌潮紊动作用复杂且对流速的垂向分布起着重要影响,采用不同的湍流模式对涌潮传播过程中水流的运动特征开展研究。通过与涌潮河段实测资料的验证,复演涌潮到达前后水流运动特征,给出涌潮水流湍动能的变化过程。研究成果有助于深入认识涌潮水流紊动特征和流速的分布规律,为涌潮作用下物质输运的研究提供基础。     

近岸海域流速与含沙量垂线分布分析

根据日照、舟山等近岸海域实测流速与含沙量资料,对流速采用构造的方式得到结构简单、适用于工程计算的流速垂线分布规律;基于工程应用的便利性原则,采用理论分析方法改进了Van Rijn含沙量垂向分布公式,并利用实测资料对其进行验证。验证结果表明,推导得到的流速垂线分布公式和改进的含沙量垂向分布公式与实测资料吻合程度较好,能够客观反映流速垂线和含沙量垂线分布基本规律。     

垂向变化的波浪辐射应力影响下流场数值研究

采用浪流沙耦合模型COHEREN-SED模型进行黄河三角洲的海流三维运动数值模拟。文中将垂向平均的波浪辐射应力与垂向变化的波浪辐射应力引入至COHERENS-SED,通过该模型进行波浪辐射应力影响下的海流场演变数值研究。研究结果显示波浪辐射应力的存在显著改变了水深小于5m的流场。垂向变化的辐射应力则使得其表层波生沿岸流速大于传统的垂向均匀辐射应力引起的波生沿岸流速,底层则刚好相反。当波浪强度很大时,流场中的波生沿岸流对水深较浅的海岸区海流场影响较大。     

海州湾湾顶浅海区ADCP测流中垂向流速分量分析

根据2012年6月在海州湾湾顶浅海区利用ADCP定点测流以及CTD等观测资料,采用经验模态分解(EMD)和Hilbert变换方法,对该海区海流的垂向流速分量及其影响因素进行了分析,研究了垂向流速分量的频率结构及其随时间的变化。结果表明,海州湾湾顶浅海区垂向流速分量变化不规则,以团块状分布为主。经验模态分解显示垂向流速分量存在不同的模态,各模态对整个波动的贡献差异明显,其中心频率构成具有垂向变频特征。垂向流速分量与水平流速分量在低频模态和趋势项具有较好的相关性,表明垂向流速分量与水平流速分量具有一定的联系。     

内孤立波波致流场数值模拟研究

基于KdV、mKdV理论,利用Fluent计算软件,采用"平板拍击"造波方法,进行内孤立波数值模拟,并与物理实验结果进行对比验证。利用数值模拟结果,分析内孤立波波致流场变化,结果表明:上下层流体中波致水平流速方向相反,均呈现先增加后减小的变化趋势,且波谷经过时刻流速最大;在波谷经过断面处,波致水平流速在上层流体中沿垂向分布无明显变化,在波面以下的下层流体中有衰减趋势,但衰减很小;两层流体界面与波谷之间存在过渡水深范围,水平流速在该水深范围内沿垂向衰减明显,且随内孤立波振幅的增大,过渡水深范围有所增大。     

源区黑潮季节变异及其动力机制的数值研究

应用POM2K模式对中国海黑潮区气候态平均环流进行了数值模拟。采用正交曲线网格, 模式区域为太平洋海盆, 特别的在中国海区域进行加密并较好的拟合了岸线; 垂向分为21层, 并在海表9层以上采用对数网格分布; 采用COADS气候态月平均的风应力, 并将模式的温度结果和MODIS月平均的SST数据进行同化, 然后将模式模拟出的流量、海表高度异常同实测数据和卫星观测数据进行了对比验证, 结果均显示模拟结果可信度较高。接着本文探讨了北赤道流分岔位置季节性的变化对源区黑潮流量的影响, 结果表明, 秋冬季节北赤道流分岔位置较靠北, 源区黑潮流量较大, 而春夏季节北赤道流分岔位置较靠南, 黑潮流量较小。在此基础上, 针对源区黑潮的动力机制进行了数值实验。实验中主要考虑了以下动力因子对源区黑潮季节性变化的影响: (1)风应力; (2)非线性; (3)黑潮的斜压敏感性, 然后通过与控制实验的对比, 讨论了不同的动力因素对吕宋海峡净流量和吕宋海峡上层环流场的影响。     

浙江近海夏季流场特征分析

为了揭示浙江近岸流场特征及沿程变化规律,于2006年和2009年夏季在浙江岸外3个固定点利用ADCP潜标进行了多个潮周期分层海流流速、流向观测。研究结果表明:(1)浙江沿岸流在中北部海域(A和B站位)为旋转流,流向呈顺时针方向旋转,在南部(C站位)涨潮流方向基本为北向,落潮方向为东偏北向;各站位海流在垂向上流向较一致。(2)3个站位垂线平均流速相近(44.4~51.1 cm/s),但平均流速的垂向分布差异明显;各站的最大流速均大于110 cm/s,且均出现在大潮涨急时刻。(3)观测期间,A(北部)、B(中北部)和C(南部)站位平均余流的大小分别为21.9,12.3和22.3 cm/s;受长江冲淡水影响,A和B站位中上层余流为西南向,从中层向底层流向呈逆时针方向偏转,下层流向呈东南向,可能为台湾暖流牵引所致,C站位余流流向在垂向较为一致,均为东北向,主要受季风影响。(4)夏季浙江沿岸流在沿浙江沿岸北上的过程中,在浙江中部(B和C站位中间)逐渐向东偏转(可能受台湾暖流的牵引),流经海域水深变大。(5)在夏季长江径流量偏小时段,浙江中北部近岸海域也存在向南的沿岸流(同冬季),其范围从长江口以南一直至浙江中北部。浙江近岸海流受季风、长江冲淡水和台湾暖流共同制约,但各区域的主要受控因素不同。     

台湾海峡西侧水动力环境的季节变化特征

本文利用2006 ～ 2007年“中国近海海洋综合调查与评价”项目在台湾海峡西侧获取的大面水温、盐度走航资料,结合座底式海床基系统获取的海流剖面以及近海底水温观测结果,系统地阐述了台湾海峡西侧水动力环境在春、夏和晚秋初冬3个季节的变化特征.研究结果表明:春季台湾海峡西侧海水垂向混合比较均一,等温线和等盐线分布大致与岸线平行,盐度和水温呈现从北往南逐渐升高、自近岸向远岸递增的特点.福建沿岸海流具有明显的正压海流特征,且为逆风的东北向海流,其最大流速(约0.25 m/s)出现在泉州外海;与此相对应,浙闽沿岸水仅局限于泉州以北近岸海域的中上水层,其盐度小于32.4.夏季台湾海峡西侧局部海域呈现典型的上升流特征,汕头、东山和南日岛外海近海底存在着明显的低温(＜25℃)高盐水(＞33.8)涌升中心.福建沿岸海流流向与盛行的西南季风一致,除汕头外海10 m层流速(0.19 m/s)比25 m层流速(约0.15 m/s)略大,其他海域仍具有正压流特点,最大流速(约0.40 m/s)出现在平潭外海附近.晚秋初冬台湾海峡西侧海域盐度和水温的空间分布比较接近,整体呈近岸低外海高,并随纬度降低而增加.调查期间浙闽沿岸水在2007年12月26日左右开始出现在平潭岛外海附近,随后从2008年1月15日左右整个台湾海峡西侧被南下的浙闽沿岸冷水(＜17℃)控制;厦门以南近海余流流速明显小于厦门以北海域,并存在着显著的垂向流剪切.     

潮滩水流加减速阶段边界层流速剖面的比较

江苏大丰地区潮滩由于水深较浅,潮流、波浪等动力较强,整个水层可视为边界层,其主体部分是对数层,即水流流速在垂向上呈对数分布。在潮流的加减速阶段,流速剖面将可能偏离对数分布,从而使对数剖面法计算出来的边界层参数造成误差。使用MIDAS-400用户化数据采集系统,在大丰潮滩获得了多层流速、浊度等同步高频观测数据,基于修正后的von Karman-Prandtl模型对u-lnz进行回归分析、数据内部一致性分析来定义流速对数剖面并与未修正前经典理论得到的边界层参数进行比较。分析结果表明,修正后的流速剖面更符合实际情况,边界层参数除了受水流加速度的影响外,还和沙纹等因素有关。另外,边界层参数的变化量与特征加速度负相关。     

Study on the double-logarithmic profile of tidal flow velocity in the near-bed layers

Tidal current velocity profile in the near-bed layers has been widely studied. The results showed that velocity profile in the near-bed layer obviously departure from the traditional logarithmic profile, due to the acceleration or deceleration. Although the logarithmic linear profile can reduce the rate of deviation from this, only it is a lower-order approximate solution. In this paper, considering the unsteady and non-linear features of tidal motion, the double logarithmic profile near-bed layers in estuarine and coastal waters is established on the assumption that the turbulent shear stress along the water depth was parabolic distribution, and on the basis of Prandtl''s mixing length theory and von Karman''s self-similar theory. Having been verified the data observed at the West Solent in the south of England, and comparison of the logarithmic linear profile, it found that the double logarithmic profile is more precious than the latter. At last, the discussed results showed that:(1) The parabolic distribution of the tidal shear stresses verified good by the field data and experimental data, can be better reflected the basic features of the tidal shear stress deviating from linear distribution that is downward when to accelerate, upward when to decelerate. (2) The traditional logarithmic velocity profile is the zero-order approximation solution of the double logarithmic profile, the logarithmic linear profile is the first order, and the logarithmic parabolic profile is the second order. (3) Ignoring the conditions of diffusion and convection in the tidal movement, the double logarithmic profile can reflect the tidal properties of acceleration or deceleration, so that the calculation of the friction velocity and roughness length are more reasonable. When the acceleration or the deceleration is about zero, the double logarithmic profile becomes the logarithmic profile.     

Study on the Log-Linear Velocity Profile of Near-Bed Tidal Current in Estuarine and Coastal Waters

Many observed data show that the near-bed tidal velocity profile deviates from the usual logarithmic law. The amount of deviation may not be large, but it results in large errors when the logarithmic velocity profile is used to calculate the bed roughness height and friction velocity （or shear stress）. Based on their investigation, Kuo et al. （1996） indicate that the deviation amplitude may exceed 100%. On the basis of fluid dynamic principle, the profile of the near-bed tidal velocity in estuarine and coastal waters is established by introducing Prandtl＇ s mixing length theory and Von Kannan selfsimilarity theory. By the fitting and calculation of the near-bed velocity profde data observed in the west Solent, England, the results are compared with those of the usual logarithmic model, and it is shown that the present near-bed tidal velocity profile model has such advantages as higher fitting precision, and better inner consistency between the roughness height and friction velocity. The calculated roughness height and friction velocity are closer to reality. The conclusions are validated that the logarithmic model underestimates the roughness height and friction velocity during tidal acceleration and overestimates them during tidal deceleration.     

Linear Approximations of the Second Turbulent Moments of the Atmospheric Convective Surface Layer in a Forced-Convection Sublayer

The Monin–Obukhov similarity theory for the convective surface layer distinguishes two limiting cases: a dynamic limit and a free-convection limit. The dynamic limit for the convective surface layer is defined as a flow with a logarithmic profile of wind and a zero buoyancy flux at the underlying surface. The free-convection limit is characterized by a zero wind speed and a positive buoyancy flux at the underlying surface. The limits of the generalized Monin–Obukhov similarity theory are able to describe the higher order turbulent moments. In this paper, it is assumed that the convective surface layer consists of two sublayers: the lower dynamic sublayer adjacent to the surface and the upper forced-convection sublayer. The turbulent moments can be approximated separately for each sublayer. Linear approximations are suggested for the turbulent moments of the vertical velocity and the potential temperature variance in the forced-convection sublayer. The first-order expansion terms of them correspond to the free-convection limits of the Monin–Obukhov theory under no-wind conditions. The second-order expansion terms describe profiles of the turbulent moments in under convective conditions with a moderate wind. A comparison between the proposed approximations and experimental data strongly suggests that the linear approximation is correct within a forced-convection range.     

Modeling a neutrally stratified turbulent air flow over a horizontal rough surface

The neutrally stratified boundary layer over a smooth rough surface is consider. The turbulent flow is simulated using a finite-difference eddy-resolving model of the atmospheric boundary layer (ABL). The model includes different turbulence closure schemes and numerical approximations for advection components of the momentum balance equation. We investigate the quality of reproduction of spectral characteristics of the turbulent flow and the model’s capabilities to reproduce the observed profile of mean wind velocity near the rough surface. It is shown that the best result is obtained by coupling a numerical scheme of higher order of accuracy with a mixed closure scheme based on an adaptive estimation of the mixing length for subgrid-scale fluctuations. Here, we are able to reproduce the asymptotics of the fluctuation spectrum of the longitudinal component of wind velocity near the surface and within the boundary layer as well as the logarithmic profile of mean velocity near the surface.     

台风期间厦门湾潮流边界层特征研究

本文基于2007年8月圣帕台风影响台湾海峡期间的厦门湾垂向海流流速剖面观测资料,运用Karman-Prandtl模型统计和分析了流速对数剖面的出现频率和边界层参数（摩阻流速与粗糙长度）．结果表明流速剖面满足对数分布出现频率都超过50％,但进行内部一致性原则分析后,流速剖面满足对数分布情况均不足30％．造成流速剖面偏离对数分布的因素主要有：台风、径流、流速加减速、水体分层、底床结构、海底悬沙浓度等等．台风期间的摩阻流速和粗糙长度都大于非台风期间的数值．     

Experimental study on internal waves in a stratified shear flow

This paper describes experiments on interfacial phenomena in a stratified shear flow having a sharp velocity shear at a density interface. The interface was visualized in vertical cross-section using dye, and the flow pattern was traced using aluminum powder. Two kinds of internal waves with different phase velocities and wave profiles were observed. They are here named p(positive)-waves and n(negative)-waves, respectively. By means of a two-dimensional visualization technique, the following facts have been confirmed regarding these waves. (1) The two kinds of waves propagate in the opposite direction relative to a system moving with the mean velocity at the interface, and their dispersion relations approximately agree with the two solutions of interfacial waves in a two-layer system of a linear basic shear flow. (2) The p-wave has sharp crests and flat troughs, and the n-wave has the reverse of this. This difference in wave profile is due to the finite amplitude effect. (3) Phase velocity of each wave lies within the range of the mean velocity profile, so that a critical layer exists and each wave has a “cat's eye” flow pattern in the vicinity of the critical layer, when observed in a system moving with the phase velocity. Consequently, these two waves are symmetrical with respect to the interface. The mechanisms of generation of these waves, and the entrainment process are discussed. It is inferred that when the “cat's eye” flow pattern is distorted and a stagnation point approaches the interface, entrainment in the form of a stretched wisp from the lower to the upper layer occurs for the p-wave, and from the upper to the lower layer for the n-wave.     

江苏大丰潮间带粉砂滩的潮流边界层特征

对江苏大丰潮间带粉砂滩的潮流观测数据，运用Karman-parndtl模型分析了流速对数剖面的出现频率和边界层参数（摩阻流速与粗糙长度）。结果表明，不同时间尺度（10^1-10^3s）的平均流速对计算u-lnz线性关系出现频率的影响不大，而不同层位组合的影响较为显著，影响流速对数剖面分布的自然因素为波浪、风、潮流旋转、表底层流向差异、滩面坡度和水流加减速等。相邻层位的流速在时间序列上存在很高的线性关系，回归直线的斜率与水层距底床的高度、粗糙长度和摩阻流速有关。涨潮阶段床面粗糙长度大于落潮阶段，涨潮阶段摩阻流速小于落潮阶段。     

Effects of time dependence in unstratified tidal boundary layers: results from large eddy simulations

A three-dimensional Large Eddy Simulation (LES) model is used to simulate oscillating tidal boundary layers and test previous results obtained from one-dimensional boundary layer models and turbulence measurements in tidal channels. The LES model produces low-order turbulence statistics in agreement with the semi-analytic theory and observations. It shows a logarithmic layer in the mean velocity profile and a linear distribution of Reynolds stress with water depth. However, the eddy viscosity profile predicted by the LES model is not parabolic but better matches a parabolic profile modified by wake effect observed in the outer part of depth-limited steady boundary layers. Low-order turbulence statistics can be scaled by the instantaneous friction velocity at the bottom boundary. Although turbulence intensities in three directions fluctuate over a tidal cycle, their normalized values are in good agreement with those determined from laboratory experiments of steady open-channel flows. The LES model confirms that tidal turbulence is in quasi-equilibrium. However, it also demonstrates the importance of flow acceleration/deceleration term in the depth-integrated momentum balance for the mean flow. Phase differences are found between flows at different heights above the bottom boundary.     

Effects of sea roughness and atmospheric stability on wind wave growth

The paper considers the effects of sea roughness and atmospheric stability on the wind wave growth by using the logarithmic boundary layer profile including a stability function, as well as adopting Toba et al.'s [J. Phys. Ocean. 34 (1990) 705] significant wave height formula combined with some commonly used sea surface roughness formulations. The wind wave growth is represented by the non-dimensional total wave energy relative to that for neutral stability used by Young [Coast. Engng 34 (1998) 23]. For a given velocity at the 10 m elevation, spectral peak period and stability parameter, the wind wave growth is determined.