Vertical Distribution of Tidal Flow Reynolds Stress in Shallow Sea

Based on the results of the tidal flow Reynolds stresses of the field observations,indoor experiments,and numerical models,the parabolic distribution of the tidal flow Reynolds stress is proposed and its coefficients are determined theoretically in this paper.Having been well verified with the field data and experimental data,the proposed distribution of Reynolds stress is also compared with numerical model results,and a good agreement is obtained,showing that this distribution can well reflect the basic features of Reynolds stress deviating from the linear distribution that is downward when the tidal flow is of acceleration,upward when the tidal flow is of deceleration.Its dynamics cause is also discussed preliminarily and the influence of the water depth is pointed out from the definition of Reynolds stress,turbulent generation,transmission,and so on.The established expression for the vertical distribution of the tidal flow Reynolds stress is not only simple and explicit,but can also well reflect the features of the tidal flow acceleration and deceleration for further study on the velocity profile of tidal flow.     

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.     

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.     

江苏近海潮流湍流切应力研究

潮流的流速分布和湍流切应力的分布密切相关。为了研究潮流湍流切应力的特征和变化规律,本文从潮流运动方程推导出随潮流呈周期性变化的湍流切应力随深度的分布表达式。从江苏近海现场实测流速剖面中计算得到一个潮周期内的湍流切应力的变化过程。通过将实测值与本文理论解的比较发现,计算值能很好地反映出实测值的变化特征,尤其是能正确地反映潮流在加速和减速过程中,湍流切应力偏离线性分布呈现上凹和下凹的现象,以及只有在加减速转换的时候切应力才呈现线性分布的特点。最后,通过分析湍流切应力振幅及相位沿水深的变化情况,认为反映潮流周期、涡粘性和水深之间关系的参数Ri是决定湍流切应力弯曲程度及分布形态的重要参数。对于近海潮流,水深是影响湍流切应力偏离线性程度的主要因素。     

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.     

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

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

Water exchange and material exchange through a strait due to tidal flow

The process of material transport through a strait due to tidal flow is modeled, and then the differences between various concepts of tidal exchange which have been used hitherto are pointed out using this model. In particular, the exchange of water itself and the exchange of material should be distinguished even in the case where the material of interest is carried by the water,i.e., the material and water move as one body. Further, the physical meaning of “tidal trapping” (Fischeret al., 1979) is discussed by using the model in this paper. The relationship between the exchange ratio for the water itself (r) and the phase lag (δ) of material concentration to the tidal stream in a section of the strait, which is an important factor in tidal trapping, is obtained as follows: $$\delta = \tan ^{ - 1} \left( {\tfrac{1}{r} - 1} \right)$$ Observational results at Lake Hamana (Shizuoka Pref) and at Kabira Cove (Okinawa Pref.) support the validity of the above relationship.     

ECOMSED模式在杭州湾海域流场模拟中的应用

针对杭州湾独特的喇叭型强潮河口湾的特点,基于Blumberg等(1996)的ECOMSED模式,引入动边界技术,建立杭州湾三维动边界的潮流模型.模型以正交曲线坐标下三维非线性水动力方程为基本方程,应用Mellor和Yamada的2.5阶湍流闭合模型计算紊动黏滞系数,嵌入Grant和Madsen的底边界层模型考虑波浪对底部应力的作用,采用干湿网格法模拟潮流漫滩过程;综合考虑径流,风应力,密度流和M₂,S₂,K₁,O₁四个主要分潮和M₄,S₄,MS₄三个浅水分潮的作用,从而提高杭州湾潮流模拟的精度.通过验潮站调和常数和多次海流连续观测资料的验证,表明该文建立的模型可以更好的用于杭州湾流场的预报模拟.     

Distribution of bottom stress and tidal energy dissipation in a well-mixed estuary

Estimates of area-averaged tidal bottom stress are made for four channel segments of the Great Bay Estuary, N.H. Current and sealevel measurements are used to estimate acceleration and pressure gradient terms in the equation of motion, while the equation of motion itself is used to infer the remaining stress term. Dynamic terms, bottom stress values, friction coefficients and energy dissipation rates are estimated for each site. The analysis shows that while throughout the estuary the principal force balance is between the frictional stress and the pressure gradient forcing, RMS values of total bottom stress range from 2·67 to 10·38 Nm^?2 and friction coefficients vary from 0·015 to 0·054. Both stress and energy dissipation are largest in the seaward portion of the estuary with an order of magnitude decrease in dissipation at the most inland site.These distributions of stress and energy dissipation are consistent with cotidal charts of the principal semi-diurnal tidal constituent (M₂) which indicate that the estuary is composed of a highly dissipative more progressive tidal wave regime seaward and a less dissipative standing wave regime landward.     

Seasonal variation and tidal influences on estuarine use by bottlenose dolphins (Tursiops aduncus)

In order to show that dolphins use estuary habitats differently depending on the season and tidal state, possibly in response to prey distribution, temperature, risk of stranding and accessibility, Indo-Pacific bottlenose dolphins (Tursiops aduncus) were observed year-round during a 3-year study in the Clarence River estuary (CR) and Richmond River estuary (RR) in northern New South Wales, Australia. Peak dolphin sightings occurred during the spring season and one or 2 h prior to high tide. The spatial distribution of the dolphins in each estuary was analysed using the distance in kilometres that the dolphins travelled upstream with seasons and tidal phase as determinants. A General Linear Model showed that in the CR the dolphin spatial distribution in the estuary was not determined by season (F = 0.434, df = 3, P = 0.729) but was by tidal phase (F = 9.943, df = 3, P < 0.001) and the interaction between season and tidal phase (F = 3.398, df = 9, P < 0.002). However, in the RR the spatial distribution of the dolphin use of the estuary was not determined by either season (F = 1.647, df = 3, P = 0.194) or tidal phase (F = 0.302, df = 3, P = 0.824). In the CR, the spatial distribution of the dolphins was largest on high and flood tides. This pattern of spatial distribution may occur because the CR is a relatively shallow estuary and this increased spatial distribution may reflect a lower stranding risk and an increase in accessibility of shallow areas during periods of higher tide. These areas could also provide access to their preferred prey items of sea mullet (Mugil cephalus) and sand whiting (Sillago ciliata).     

福建省罗源湾互花米草扩展过程及其特征分析

为探讨外来引种的互花米草（Spartina alterniflora Loisel）在我国海湾的扩展过程,根据1989-2010年遥感影像和野外观测数据,对罗源湾互花米草的分布特征进行分析,结果显示,对该区互花米草近20 a的扩展过程划分为3个阶段:前5 a其沿着高潮带平行于海岸线呈带状扩展;之后5 a在中潮带向湾内快速拓殖;后10 a则扩展缓慢,种群内部的空隙地逐渐被米草覆盖。截至2010年,罗源湾互花米草的总面积约为741.13 hm2。受潮滩环境分带性影响,罗源湾互花米草在高潮带平行于海岸线的扩展速度明显大于垂直于海岸线的扩展速度;在中潮带的快速扩展过程中互花米草倾向于先在潮沟两侧分布,这可能与归槽水和滩面摩阻作用有关,相比较,互花米草在浅洼地的扩展较滞后,该区域随着滩面的淤高过程逐渐被米草覆盖。在其整个扩展过程中,由于互花米草对潮流及悬浮颗粒的阻滞作用,使得其扩展前缘米草斑块和连续种群之间的光滩区成为悬浮颗粒及米草种子容易沉降的区域,互花米草易朝向该区域扩展。     

Vertical Distribution of Momentum Exchange Coefficient and Sediment Concentration in Estuarine and Coastal Waters

Vertical structures of momentum exchange coefficient and sediment concentration are the keys in the research on estuarine and coastal suspended sediment transport. Based on the parabolic mixing length distribution pattern,the distribution pattern of vertical momentum exchange coefficient which is suitable for estuarine and coastal waters is constructed. A comparison with steady flow and measured momentum exchange coefficient during one tidal cycle in the Menai Strait of England shows that the result of this...     

挡潮闸矩形孔口式鱼道紊流结构及放鱼试验研究

对挡潮闸枢纽中矩形中孔、底孔鱼道中紊流结构进行了较为系统的试验研究,并做了放鱼试验。选择了一种鱼类偏爱流速所对应的流量作为典型流量,考虑了不同的孔口位置(中孔和底孔),用声学多普勒测速仪(ADV)量测了测点的三维瞬时流速及流向,分析了矩形孔口鱼道的三维时均流速分布特征、断面最大流速沿程变化规律、流速矢量场、紊动强度分布及雷诺应力分布。此外,还通过放鱼试验,利用在鱼体植入T形标签和高速摄影机观察了过鱼对象对中孔、底孔的反应情况,分析了过鱼对象与矩形孔口鱼道紊流结构的关系。试验结果表明:水流经中孔形成三维紊动自由射流,经底孔形成三维壁面射流,中孔纵向流速呈高斯分布,而底孔纵向流速则近似为高斯分布,流速由孔口向两侧逐渐减小;中孔和底孔横向流速在位于孔口范围内的纵剖面上沿程减小,孔口之外则变化较小;中孔和底孔垂向流速分布特征表现为在铅垂方向上均存在旋涡;在中孔水平面和纵剖面上,纵向最大流速均沿程衰减;中孔和底孔情形孔口处紊动强度和雷诺应力比非孔口处大得多,而非孔口处不同水深平面上紊动强度和雷诺应力变化趋于平缓;过鱼对象喜爱在紊动强度分布的峰值区和雷诺应力较大变幅区溯游。     

On salinity regimes and the vertical structure of residual flows in narrow tidal estuaries

An examination is made of the circulation in narrow estuaries subject to a predominant tidal forcing. Velocity structures are derived separately for residual flow components associated with (a) river flow, (b) wind stress, (c) a well-mixed longitudinal density gradient and (d) a fully stratified saline wedge. Dimensionless parameters are introduced to indicate the magnitude of each component and these parameters are evaluated for 9 major estuaries, thereby revealing their sensitivity to each component.For a channel of constant breadth and depth, formulae are deduced for the length of saline intrusion, L. Comparisons with observed data show that such formulae may be used with confidence to predict changes in L arising from variations in river flow, tidal range or channel depths.The level of stratification is shown to be related to a product of two parameters, one associated with velocity structure and a second involving the square of the ‘flow ratio’ $\text{u}\text{u}\text{?}$ (i.e. residual velocity/amplitude of the tidal velocity). This relationship provides a simple classification system for estuarine stratification which can be used to indicate the sensitivity of any particular estuary to changing conditions.     

Variations in the boundary-drag coefficient in the tidal entrance to Chesapeake Bay,Virginia

Use of the quadratic shear-stress law for estimating boundary drag requires specific knowledge of the magnitude of a drag coefficient, C_D, and sectional mean velocity, $\text{u}\text{?}$. In previous attempts to adapt the relationship for use in studies of marine-sediment transport, the flow measurement has been standardized at a level 100 cm above the bed. The particularized value of the drag coefficient has been designated as C₁₀₀.In the entrance area to Chesapeake Bay, Virginia, C₁₀₀ has been found to range through unacceptably wide limits. Two-thirds of the values obtained are between 3.5 · 10^?3 and 5.4 · 10^?2. Mean C₁₀₀ for the area is 1.3 · 10^?2 as compared to 3 · 10^?3 for tidal channels within Puget Sound, Washington.Present data suggest that, given a moveable bed, a size hierarchy of mobile bed forms, time-varying flow, and a lack of equilibrium between flow and bed, C₁₀₀ changes continuously with boundary shear stress.Accurate evaluation of boundary shear stress in tidal entrances with high flow rates and mobile beds presently requires measurement of velocity profiles.     

Effect of wave amplitude on turbulent flow in a wavy channel by direct numerical simulation

The present study numerically investigates the characteristics of three-dimensional turbulent flow in a wavy channel. For the purpose of a careful observation of the effect of the wave amplitude on the turbulent flow, numerical simulations are performed at a various range of the wave amplitude to wavelength ratio (0.01?α/λ?0.05), where the wavelength is fixed with the same value of the mean channel height (H). The immersed boundary method is used to handle the wavy surface in a rectangular grid system, using the finite volume method. The Reynolds number (Re=U_bH/ν) based on the bulk velocity (U_b) is fixed at 6760. The present computational results for a wavy surface are well compared with those of references. When α/λ=0.02, the small recirculating flow occurs near the trough at the instant, but the mean reverse flow is not observed. In the mean flow field, the reverse flow appears from α/λ=0.03 among the wave amplitude considered in this study. The domain of the mean reverse flow defined by the locations of separation and reattachment depends strongly on the wave amplitude. The pressure drag coefficient augments with increasing the wave amplitude. The friction drag coefficient shows the increase and decrease behavior according to the wave amplitude. The quantitative information about the flow variables such as the distribution of pressure and shear stress on the wavy surface is highlighted.     

清水沟北汊流路入海泥沙对东营港影响的数值分析

本文用平面二维潮流泥沙数学模型,模拟了清水沟流路北汊入海口海域的潮流泥沙扩散和河口淤积延伸方向等问题。在此基础上从流速场分布,泥沙扩散浓度分布及淤积厚度零米线的范围及河口淤积延伸方向等方面探讨了入海泥沙对东营港的影响。结论认为如果按此口门入海,入海泥沙不会对黄河海港产生直接淤积影响。     

Laboratory observation of boundary layer flow under spilling breakers in surf zone using particle image velocimetry

A boundary layer flow under spilling breakers in a laboratory surf zone with a smooth bottom is investigated using a high resolution particle image velocimetry (PIV) technique. By cross-correlating the images, oscillatory velocity profiles within a viscous boundary layer of O(1) mm in thickness are resolved over ten points. Using PIV measurements taken for an earlier study and the present study, flow properties in the wave bottom boundary layer (WBBL) over the laboratory surf zone are obtained, including the mean velocities, turbulence intensity, Reynolds stresses, and intermittency of coherent events. The data are then used to estimate the boundary layer thickness, phase variation, and bottom shear stress. It is found that while the time averaged mass transport inside the WBBL is onshore in the outer surf zone, it changes to offshore in the inner surf zone. The zero Eulerian mass transport occurs at h/h_b ≈ 0.92 in the outer surf zone. The maximum overshoot of the streamwise velocity and boundary layer thickness are not constant across the surf zone. The bottom shear stress is mainly contributed by the viscous stress through mean velocity gradient while the Reynolds stress is small and negligible. The turbulence level is higher in the inner surf zone than that in the outer surf zone, although only a slight increase of turbulent intensity is observed inside the WBBL from the outer surf zone to the inner surf zone. The variation of phase inside and outside the WBBL was examined through the spatial velocity distribution. It is found the phase lead is not constant and its value is significantly smaller than previous thought. By analyzing instantaneous velocity and vorticity fields, a remarkable number of intermittent turbulent eddies are observed to penetrate into the WBBL in the inner surf zone. The size of the observed large eddies is about 0.11 to 0.16 times the local water depth. Its energy spectra follow the − 5/3 slope in the inertial subrange and decay exponentially in the dissipation subrange.