兴化湾及周边海域潮流动力特征与数值模拟研究

基于2011年兴化湾及周边海域的流速、流向的观测资料,应用短期资料的潮流准调和分析方法,对海流周日连续观测数据进行分析,探讨海域潮流特征。最后,针对兴化湾岸线复杂、岛屿众多的特点,基于FVCOM数值模式,采用非结构三角形网格,建立兴化湾及周边海域潮流动力数值模型,经实测站潮汐、潮流数据对模拟结果验证,模型能够较好的刻画各测站潮汐、潮流过程特点,反映研究海域流场特征规律,可为该海域海洋环境保护和海洋开发利用提供科学论证依据。     

Viscous free surface effect on coastal embankment hydrodynamics

A finite-differnece method was used to calculate the nonlinear hydrodynamic pressures acting on the coastal embankment faces by seismic-wave actions. The nonlinearity of free surface flow, convective acceleration, viscosity and surface tension of fluid are included in the analysis. The kinematic and dynamic free surface boundary conditions are employed for calculating the horizontal fluid velocity, pressure at the free surface and the surface profile of the fluid. The time-dependent water surface is transformed to the horizontal plane, and the flow field is mapped onto a rectangular, making it convenient to model the complex sea bottom geometry and the wavy water surface by the finite-difference method. Fully nonlinear and weakly nonlinear dynamic free surface conditions are used and compared. The effects of surface tension of fluid are also discussed. The nonslip boundary condition is applied on the most part of the interface between fluid and solid face, except the region near the intersection between free surface and wall face. The numerical results are presented for various water depths and ground motion intensities, and their associate viscous effects on coastal embankment hydrodynamics are discussed.     

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.     

Field investigation on the results of non-linear shallow water equations in the Swash Zone

For the study of the cross-shore wave-induced hydrodynamics in the swash zone, a numerical model is developed based on the one-dimensional non-linear shallow water (NSW) equations for prediction of hydrodynamic parameters in the swash zone. In order to evaluate the accuracy of the outputs of the numerical model, the model's predictions in terms of water surface elevations and cross-shore velocities, are compared to field data from full-scale experiments conducted on three sites with different beach slope; mild and steep, several bed particle sizes and under various incident wave conditions. The quantitative and qualitative comparison of the results of the numerical model and the full-scale data reveals that the model can generally predict many aspects of the flow in the surf and swash zone on both types of beach. The accuracy is adequate for application in a sediment transport study. Considering the time-history and probability distribution of water surface elevation, the model is generally more accurate on steep beaches than on the mild beach. The model can adequately simulate the dominant frequency across the beach and saturation of higher frequencies on both mild and steep beaches for various incident wave energy characteristics. With regard to the horizontal (cross-shore) velocity, the sawtooth shape of time-history and negative acceleration of water are well predicted by the model for both mild and steep beaches. Due to the uncertainties in maximum and minimum values of velocity data, clear judgement about the accuracy of the numerical model in this matter was not possible. However, the comparison of the minimum velocities (offshore direction) revealed that the application of friction factors below the range which is suggested by literature best match the data.     

Simulation of the Fate of Faecal Bacteria in Estuarine and Coastal Waters Based on A Fractionated Sediment Transport Model

A two-dimensional depth-integrated numerical model is refined in this paper to simulate the hydrodynamics, graded sediment transport process and the fate of faecal bacteria in estuarine and coastal waters. The sediment mixture is divided into several fractions according to the grain size. A bed evolution model is adopted to simulate the processes of the bed elevation change and sediment grain size sorting. The faecal bacteria transport equation includes enhanced source and sink terms to represent bacterial kinetic transformation and disappearance or reappearance due to sediment deposition or re-suspension. A novel partition ratio and dynamic decay rates of faecal bacteria are adopted in the numerical model. The model has been applied to the turbid water environment in the Bristol Channel and Severn estuary, UK. The predictions by the present model are compared with field data and those by non-fractionated model.     

河口海岸三维双σ坐标斜压水流数值模型研究

在海底地形陡变、垂向密度分层明显的水域,三维σ坐标模式中会出现一种"伪"水平斜压梯度力,并会引起"伪"密度流,以至于影响模拟的精度。垂向上引入双σ坐标变换,建立河口海岸水域三维斜压水流数值模型。数值试验结果表明,在海底地形陡变水域,双σ坐标模式可以减小水平斜压梯度力处理引起的误差。     

Three-dimensional tidal current numerical model of the Oujiang Estuary

The characteristics of three-dimensional (3-D) tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, a 3-D numerical model on an unstructured triangular grid has been developed. The σ coordinate transformation, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simu...     

基于半分析方法的黄河口悬浮物浓度遥感反演

河口区的悬浮物浓度受陆源输入和水动力等其他因素的影响严重,具有重要研究意义。本文利用黄河口及其邻近区域采集的145个站位光谱数据和悬浮物(Suspended Particulate Matter,SPM)浓度数据,检验了Nechad模型、多波段准分析算法(Quasi-Analytical Algorithm,QAA)、最优化模型和半经验辐射传输模型(Semi-Empirical Radiative Transfer,SERT)在黄河口水域的适用性。结果表明,QAA561模型和Nechad561模型在低浓度水域反演结果较好,QAA655模型反演结果较差,Nechad665模型适合于低浓度水域,但反演精度低。QAA865模型和Nechad865模型在中高浓度水域的反演结果较好。SERT655模型反演精度较低。最优化模型和SERT滑动模型的反演精度较高。由于黄河口水域悬浮物浓度变化范围较大,因此,本文建立了分段QAA模型和Nechad模型。分段模型的反演结果均优于单一模型。根据误差敏感性分析可知,Nechad分段模型和QAA分段模型对于遥感反射率的50%以内的误差不敏感,稳定可靠。将分段模型应用于Landsat-8 OLI数据,获取了悬浮物浓度时空分布图。结果表明,反演结果与已有研究具有相似的分布特征。两种模型得到的悬浮物浓度在数值上存在差异,然而分布趋势在中高浓度季节有非常好的一致性。     

长江口南槽最大浑浊带枯季大小潮悬沙峰特征及其动力机制

长江口最大浑浊带是陆海交汇的核心区域,其航槽是扼海-河联运的咽喉,悬沙峰的涨落潮周期变化深刻影响航槽的稳定性。本文利用长江口南槽上、中、下段3个站点枯季小潮和大潮的流速、盐度、悬沙平均粒径和悬沙浓度的实测资料,分析最大浑浊带悬沙峰特征及其动力机制。发现:流速和滩槽交换增强导致大潮平均悬沙浓度比小潮增加了0.78—1.97倍,絮凝也导致憩流底层悬沙浓度增加8%左右,但流速和絮凝与悬沙浓度的关系均非线性。大小潮盐度梯度与底层悬沙浓度关系呈现高线性相关关系,表明盐度梯度强化或突变是泥沙再悬浮形成悬沙峰的主要动力。     

无结构网格二维河口海岸水动力数值模式的建立及其应用

为完全拟合河口近海复杂岸线和工程结构以及有效局部加密,设计并建立了一个无结构三角形网格二维河口海岸水动力数值模式。空间离散主要基于有限体积法以保证守恒性,时间积分采用预估修正法以提高精度。水位在三角形网格中心通过连续方程求解;水平x方向和y方向的流速U和V均在网格边中点上通过动量方程求解。流速平流项的求解中采用了TVD格式。TVD流速平流通量为一个一阶迎风格式通量和一个二阶格式通量的组合,一阶格式通量和二阶格式通量根据流速的局部分布情况得出配比,最终组合得到TVD通量。TVD格式具有低耗散和无频散的优点,提高了模式的稳定性。应用实测资料验证建立的模式,结果显示水位、流速和流向的计算值与实测值均符合良好。     

Conjunction of 2D and 3D modified flow solvers for simulating spatio-temporal wind induced hydrodynamics in the Caspian Sea

The main objective of this study is the simulation of flow dynamics in the deep parts of the Caspian Sea, in which the southern and middle deep regions are surrounded by considerable areas of shallow zones. To simulate spatio-temporal wind induced hydrodynamics in deep waters, a conjunctive numerical model consisting of a 2D depth average model and a 3D pseudo compressible model is proposed. The 2D model is applied to determine time dependent free surface oscillations as well as the surface velocity patterns and is conjunct to the 3D flow solver for computing three-dimensional velocity and pressure fields which coverage to steady state for the top boundary condition. The modified 2D and 3D sets of equations are conjunct considering interface shear stresses. Both sets of 2D and 3D equations are solved on unstructured triangular and tetrahedral meshes using the Galerkin Finite Volume Method. The conjunctive model is utilized to investigate the deep currents affected by wind, Coriolis forces and the river inflow conditions of the Caspian Sea. In this study, the simulation of flow field due to major winds as well as transient winds in the Caspian Sea during a period of 6 hours in the winter season has been conducted and the numerical results for water surface level are then compared to the 2D numerical results.     

湖泊水库湖流模拟系统

在系统研究湖泊、水库各项水动力影响因素的基础上,建立了湖流模拟系统,其中包括水流叠加波浪后波流的共同作用、波浪剩余动量流、风作用剪应力、风吹流以及悬移质扩散分布等。数值模拟中还采用了几种处理技术,并对ADI法进行了双向双步扫描改进。在引滦入津输水工程中的于桥水库湖流计算预报的应用表明,计算结果与现场实测吻合良好,具有计算精度高、收敛快、稳定性好、通用性强等优点。     

A Vertical Two-Dimensional Model to Simulate Tidal Hydrodynamics in A Branched Estuary

A vertical (laterally averaged) two-dimensional hydrodynamic model is developed for tides, tidal current, and salinity in a branched estuarine system. The goveming equations are solved with the hydrostatic pressure distribution assumption and the Boussinesq approximation. An explicit scheme is employed to solve the continuity equations. The momentum and mass balance equations are solved implicitly in the Cartesian coordinate system. The tributaries are govemed by the same dynamic equations. A control volume at the junctions is designed to conserve mass and volume transport in the finite difference schemes, based on the physical principle of continuum medium of fluid. Predictions by the developed model are compared with the analytic solutions of steady wind-driven circulatory flow and tidal flow. The model results for the velocities and water surface elevations coincide with analytic results. The model is then applied to the Tanshui River estuarine system. Detailed model calibration and verification have been conducted with measured water surface elevations,tidal current, and salinity distributions. The overall performance of the model is in qualitative agreement with the available field data. The calibrated and verified numerical model has been used to quantify the tidal prism and flushing rate in the Tanshui River-Tahan Stream, Hsintien Stream, and Keelung River.     

河口大型围垦工程中龙口水动力特点

河口大型围垦工程围区吞吐潮量大,河口滩势条件复杂,龙口位置选择和极值流速确定是工程设计中的关键环节。以两个典型河口大型围垦工程龙口为例,从水动力分析角度提出合理的龙口平面选址和龙口渡汛时因地制宜的结构布置。在青草沙水库工程中,龙口位置的确定充分利用围区原有深槽的过流能力,以保证库区内外及时水交换,并兼顾到围堤实施顺序与工程区整体河势环境相协调,避免口门进出水引起工程河段滩地的大冲大淤,也确保了堤基安全。龙口渡汛需要面临长时间大潮汛过流考验,龙口流速大小决定了龙口结构的保护方案和后期的合龙方案,常规使用的堰流计算方法能较为准确地计算龙口流量和断面平均流速过程,数值模型能对大型龙口流速空间分布情况进行很好的模拟,是常规计算方法的重要补充。计算结果表明口门横向上中心流速大于口门两侧流速,在纵向上底坡内外两侧顶角处水流受重力作用加强,垂向断面收缩,在涨、落急时刻分别形成大流速区,是龙口结构布置时重点抗冲保护区域,需要设置抗冲性较强、自重大且联接牢固的护面材料。此外,同一潮周期内涨急流速通常大于落急流速也是其重要水动力特点。     

河口海岸泥沙输移的分组数学模型研究

根据河口海岸水沙输移的特点,建立了一个新的二维分组数学模型,用来预测该区域的水沙输移过程。该模型耦合了水动力模块、泥沙输移模块和床面演变模块。其中水动力模块基于浅水方程组,综合考虑了柯氏力、床面切应力以及表面风应力的影响,引入干湿判断法处理动边界。泥沙输移模块首先将泥沙按照粒径分组,针对不同泥砂性质,对各组泥沙分别进行建模求解。床面演变模块基于质量守恒方程,实时更新床面高程以及床沙级配变化,并传递给水动力模块,更新底部边界。该模型被应用在了英国塞汶(Severn)河口,其预测的泥沙浓度和实测数据以及不分组的模型的预测结果进行了比较,结果显示,文中建立的分组模型预测的结果要明显好于不分组模型。     

Dynamic Analysis for the Global Performance of An SPM-Feeder-Cage System Under Waves and Currents

In the present study, the dynamic response of a coupled SPM-feeder-cage system under irregular waves and shear currents is analyzed. A numerical model is developed by using the commercial software Orca Flex. Hydrodynamics coefficients of the vessel are calculated by using a 3D diffraction/radiation panel program. First- and second-order wave forces are included in the calculations. Morison equation is used to compute the drag force on line elements representing the net. Drag coefficients are determined at every time step in the simulation considering the relative normal velocity between the structural elements and the fluid flow. The dynamic response of the coupled system is analyzed for various environments and net materials. The results of the study show the effects of solidity ratio of the net and vertical positions of the cage on the overall dynamic response of the system, confirming the viability of this type of configuration for future development of offshore aquaculture in deep waters.     

秦皇岛河流-海岸水动力和水质耦合模拟分析

基于MIKE FLOOD建立秦皇岛入海河流及其近岸海域一维、二维耦合河流-海岸水动力和水质模型,研究了秦皇岛海域和入海河流的水动力和污染物分布规律。结果表明:(1)海域潮流为半日潮,其中石河口南北支入海口流速变化过程差别较大,但水位变化过程基本一致;(2)整体上河流入海口化学需氧量(COD)与流速呈正相关关系;(3)COD在近岸海域的扩散方向与涨落潮潮流方向一致,各入海河口海域落潮流期间的COD高于涨潮流。如上结论为赤潮爆发机理和应急处置研究提供科学依据和技术支持。     

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.     

Adaptive numerical analysis of current velocities in the shelf zone within the framework of a quasiisopycnic model

We propose a method for the adaptive numerical analysis of current velocities in a basin including a shelf zone and having a lengthy open boundary with the sea, at which the current velocity is not given. This method is based on a multilayer quasiisopycnic model. We use an approach according to which one introduces additional terms of the source type into the equations of heat and salt transfer and the equation of evolution of the thicknesses of layers. We describe the results of calculations of the current velocities in the northwest part of the Black Sea in June on the basis of climatic data. We also analyse specific features of thermohaline and dynamic structure of waters. Translated by Peter V. Malyshev and Dmitry V. Malyshev     

Hydrodynamic Study on Energy Capturing Performance of Horizontal Axis Blades Under Sub-Low Speed Tidal Current

The research on the hydrodynamics of blades is mainly focused on sea areas with high-speed current. However, the average velocity in most territorial waters of China is smaller than 1 m/s, and the lift type of airfoil blades has limited application in most of these conditions. Therefore, it is of great significance to study the tidal current energy capture of blades in sub-low speed sea areas. The effect of flow impact resistance on the blade at sub-low current speed is considered and a new type of thin-walled blade based on the lift type of blade is proposed, and then the lift-impact combined hydrodynamic model of horizontal axis blade is established. Based on this model, and considering the characteristics of tidal current and velocity in the sea area of Yushan Islands, simulation and optimization of blade design are carried out. Additionally, the horizontal axis thin-walled blade and the NACA airfoil contrast blade under the same conditions are developed. By using a synthetical experimental test system, the power, torque, rotational speed and load characteristics of these two blades are tested. The performance of the thin-walled blade and the design theory are verified. It shows that this type of blade has much better energy capture efficiency in the sub-low speed sea area. This research will promote the study and development of turbines that can be used in low-speed current sea areas in the future.