象山港水交换数值研究──Ⅰ.对流-扩散型的水交换模式

以溶解态的保守性物质作为湾内水的示踪剂,建立了对流－扩散型的海湾水交换数值模型。数值模型使用参数化的方法把重力环流和潮振荡的垂向剪切作用的水平混合效应包纳在水平二维的示踪剂对流－扩散方程中。在空间网距较小时,模型的稳定性和守恒性均可满足海湾水交换研究的需要。     

象山港水交换特性研究

在验证良好的三维斜压潮流数学模型的基础上,以溶 解态的保守性物质为示踪剂,建立对流-扩散型的海湾水交换数值模型,计算了象山港水体半交换时间和平均滞留时间,并研究了斜压动力对湾内外水交换的贡献。研究结果表明,象山港水交换速度的区域性变化较大,水体半交换时间和平均滞留时间由象山港口门向湾顶逐渐增加,口门附近半交换时间在5d以内,平均滞留时间为5~10 d;湾顶水交换速度缓慢,水体半交换时间为30~35 d,平均滞留时间为35~40d。斜压动力对狭湾外段水交换影响较弱,对狭湾内段有较大的影响。     

Numerical Simulation of Residual Circulation due to Bottom Roughness Variability Under Tidal Flows in A Semi-Enclosed Bay

1 .IntroductionWater pollution has a large influence on water environment and our life ;for example ,it some-times induces oxygen depression,blue tide andredtide in a bay.For water purificationin a semi-en-closed sea ,some new methods to activate the tidal exchange between an inner bay and an outer seaarea by control of a tidal residual current have been proposed (Awayaet al .,1995) . Huthnance(1973) ,Zimmerman (1978 ,1981) , Loder (1980) ,and Robinson (1983) considered cases wherethe fricti…     

港口内靠码头系泊船运动的计算

本文以英国南海岸处Ｓｈｏｒｅｈａｍ港内系泊船为例，研究了港口内靠码头系泊船运动的数值计算问题。该船的实船实验和模型实验已经完成，为验证本文计算结果提供了依据。运动方程在时域内求解。在计算船体脉冲响应函数时，引入了船体阻尼系数在频率很低时的渐近表达式和一迭代算法。分析了港内共振波浪和其对船体的作用力。讨论了靠码头系泊船运动的非线性特征。计算结果与实验结果符合很好。     

Effects of a marine ridge to western boundary current in a three-dimensional source-sink flow model

A numerical experiment is made with a three-dimensional source-sink flow model in order to study effects of a marine ridge to the western boundary current. The model represents the Kuroshio current crossing the Izu Ridge which has a submarine narrow passage. It is indicated that most of the surface water of the western boundary current crosses the ridge through the narrow passage. A part of the surface water circulates around a warm eddy formed in an upstream region of the ridge. A water below the thermocline makes a long detour around the ridge and forms a subsurface northward current along the eastern slope of the ridge.     

Numerical Simulation of Local Scour Around A Large Circular Cylinder Under Wave Action

A horizontal two- dimensional numerical model is developed for estimation of sediment transport and sea bed change around a large circular cylinder under wave action. The wave model is based on an elliptic mild slope equation. The wave-induced current by the gradient of radiation stress is considered and a depth integrated shallow water equation is applied to the calculation of the current. The mass transport velocity and the bed shear stress due to streaming are considered, which are important factors affecting the sediment transport around a structure due to waves, especially in reflective areas. Wave-current interaction is taken into account in the model for computing the bed shear stress. The model is implemented by a finite element method. The results of this model are compared with those from other methods and agree well with experimental data.     

日照港工程对海水环境影响的数值预测

日照港位于山东省南部日照市石臼湾畔,濒临南黄海,港域开阔水深,常年不冻，又无淤积之患，,陆域西通中原和欧亚大陆,是一个不可多得的天然良港。石臼湾水深由西向东逐渐加深,离岸3km处的最大水深即达10m以上,平均水深也在6-8m之间。湾内无大径流注入,几条季节性小河兼作工业废水和生活污水的排放沟,流量为500t/d，其中COD含量约为100mg/L。 日照港始建于1979年,一期工程由万平口向南,建成开敝式深水栈桥煤炭专用码头。栈桥墩对海水动力条件改变很小,所以对湾内海水物理自净能力的影响甚微。二期工程以后,码头全选用重力式结构。正在建设中的东港区码头又位于石臼湾中间,将海湾西部海岸与外海隔开,破坏了原来湾内的海水动力环境、使石臼镇的工业废水和生活污水先排入港区,再外排黄海,污水在港池内滞留,恶化了港区水质,并造成湾西部水质质量下降。本文通过数值预测,给出了日照港二期工程以后各建设阶段工程对湾内海水物理自净能力的影响和水质恶化程度。     

电厂温排水中余氯浓度预测

利用有限元分步杂交方法,对胶州湾海域潮流场进行了数值模拟。考虑余氯在水体中的化学反应,利用前人大量实验成果选取余氯衰减系数,建立胶州湾余氯的二维输运-扩散模型,并将此模型应用于青岛市黄岛发电厂三期温排水工程,预测了温排水中余氯的浓度分布及影响范围,为电厂冷却水排放口邻近水域的生态环境保护提供依据。     

环抱式港区水体交换能力数值研究——以闸坡渔港为例

环抱式港池的水体交换能力强弱对港池的水质有重要影响,水体交换通道是改善港池水质的有效措施。本研究以闸坡渔港防波堤改造工程为例,采用不规则三角网格和有限体积方法,在潮流数模研究成果的基础上,建立平面二维对流-扩散数学模型,开展闸坡渔港扩建前后水体交换能力数值模拟研究,评估工程前后港内水体半交换周期和水体交换率。研究结果表明,闸坡渔港的水体交换能力较强,工程前水体半交换时间为129 h,工程后闸坡渔港水体半交换时间为104 h。工程后洲仔峡防波堤被改为涵洞式桥梁,使渔港与外海形成了水体交换通道,增加了港池内外水体相互作用和水体交换量,相应的水体交换能力增强。而港池外的两个水产养殖区域由于与外海相连,水体交换能力较强,工程前后该区域水体交换能力变化不大。     

Water levels in a dual-basin harbour in response to infragravity and edge waves

Seiche modes in a compound harbour (an “Outer Harbour” connected both to the sea and to an “Inner Harbour”) were studied using water level data and a numerical model. A variety of harbour oscillations are present, with periods up to 67 min. Periods longer than 25 min exceed resonant modes of the harbour. This paper addresses the characteristics and causes of the open-basin modes. The dual harbour open-basin mode is modified by constriction at the connection between harbours, by partial reflection at the antinode, and by the geometry of the entrance. The single-harbour open-basin mode excites the dual harbour closed-basin mode, which has nearly the same period. This forcing moves the closed-basin antinode and slightly changes the modal period, but the coupling permits the amplitude to increase through the closed-basin resonance. The water level response to wind stress is weak, but significant residual currents can occur, which take the form of clockwise gyres in each basin. Energetic peaks in the water level spectrum at 26, 35, and 67 min are shown to correspond to possible edge waves on the local shelf. The work has practical implications to port design, e.g. towards minimisation of ship ranging while at anchor.     

Numerical Experiment on Two-Dimensional Line Thermal

The time evolution of a two-dimensional line thermal-a turbulent flow produced by an initial element with signifi-cant buoyancy released in a large water body, is numerically studied with the two-equation k - e model for turbulence closure. The numerical results show that the thermal is characterized by a vortex pair flow and a kidney shaped concentra-tion structure with double peak maxima; the computed flow details and scalar mixing characteristics can be described by self-similar relations beyond a dimensionless time around 10. There are two regions in the flow field of a line thermal: a mixing region where the concentration of tracer fluid is high and the flow is turbulent and rotational with a pair of vortex eyes, and an ambient region where the concentration is zero and the flow is potential and well-described by a model of doublet with strength very close to those given by early experimental and analytical studies. The added virtual mass coeffi-cient of the thermal motion is found to be approximat     

Numerical models for evolution of extreme wave groups

The paper considers the application of two numerical models to simulate the evolution of steep breaking waves. The first one is a Lagrangian wave model based on equations of motion of an inviscid fluid in Lagrangian coordinates. A method for treating spilling breaking is introduced and includes dissipative suppression of the breaker and correction of crest shape to improve the post breaking behaviour. The model is used to create a Lagrangian numerical wave tank, to reproduce experimental results of wave group evolution. The same set of experiments is modelled using a novel VoF numerical wave tank created using OpenFOAM. Lagrangian numerical results are validated against experiments and VoF computations and good agreement is demonstrated. Differences are observed only for a small region around the breaking crest.     

Verification of A Numerical Harbour Wave Model

A numerical model for wave propagation in a harbour is verified by use of physical models.The extended time-dependent mild slope equation is employed as the governing equation,and the model is solved by use of ADI method containing the relaxation factor.Firstly,the reflection coefficient of waves in front of rubble-mound breakwaters under oblique incident waves is determined through physical model tests,and it is regarded as the basis for simulating partial reflection boundaries of the numerical model.Then model tests on refraction,diffraction and reflection of waves in a harbour are performed to measure wave height distribution.Comparative results between physical and numerical model tests show that the present numerical model can satisfactorily simulate the propagation of regular and irregular waves in a harbour with complex topography and boundary conditions.     

Validation of the three-wave quasi-kinetic approximation for the spectral evolution in shallow water

This paper aims at validating the three-wave quasi-kinetic approximation for the spectral evolution of weakly nonlinear gravity waves in shallow water. The problem is investigated using a one-dimensional numerical wave propagation model, formulated in the spectral representation. This model includes both a nonlinear triad interactions term and a wave breaking dissipation term. Some numerical tests were carried out in order to show the importance of using the triad nonlinear term in wave propagation spectral models, particularly to describe both behavior of the spectral integral parameters and of the spectral shape evolution in shallow water depth. Furthermore; a comparison against different set of experimental observations was carried out. Comparing the numerical results with the experimental observations made it possible to show the modeling efficiency of the three-wave quasi-kinetic approximation.     

象山港分区水交换数值研究

在验证良好的三维斜压潮流数学模型的基础上,将象山港划分成7个海区,以7种不同的溶解态保守性物质为示踪剂,建立对流-扩散型的海湾水交换数值模型,采用关联矩阵法来描述各区之间以及各区与外海之间的水交换特性。通过不同示踪剂的浓度分布,计算了15、30、45、60和80d的水交换矩阵。由各时刻水交换矩阵可知,在象山港的7个子区中与湾外水的水交换速度最快的为Ⅰ区,最慢的为Ⅶ区。Ⅲ区(西沪港)、Ⅵ区(黄墩港)和Ⅶ区(铁港)内水体与湾外海水之间的水交换速度缓慢,但与区域外水体之间的交换速度要明显快于其他几个海区。从水交换矩阵可以看出西沪港、黄墩港和3个内港之间相互交换的海水量很小。     

Three-dimensional numerical experiments on tidal exchange through a narrow strait in a homogeneous and a stratified sea

We have investigated the three-dimensional Lagrangian motion of water particles related with tidal exchange between two basins with a constant depth connected through a narrow strait and the effects of density stratification on the exchange processes by tracking a number of the labeled particles. Tide-induced transient eddies (TITEs), which are similar to those in two-dimensional basin, are generated behind the headlands. Upwelling appears around the center of the eddy and sinking around the boundary. When the basins are filled with homogeneous water, a pair of vortices are produced in the vertical cross section of the strait due to bottom stress, with upwellings along the side walls of the strait and sinking in the center of the strait. These circulations form the horizontally convergent field in the cross-strait direction in the upper layers while the horizontal divergence takes place in the bottom layer. These vertical water-motions produce the three-dimensional distribution of velocity shear and phase lag of the tidal current around the strait, and the Lagrangian drifts of water particles become large. As a result, water exchange through the strait is greatly enhanced: The water exchange rate reaches 94.1% which is much larger than that obtained in the vertically integrated two-dimensional model. When the basins are stratified, the stable stratification suppresses the vertical motion so that a pair of vertical vortices are confined in the lower layers. This leads to a decrease in the exchange rate, down to 88.6%. Our numerical results show that the three-dimensional structure of tidal currents should be taken into account in tidal exchange through a narrow strait.     

一个改进的二阶Boussinesq方程模型在线性波浪反射问题中的适用性研究

基于改进型的二阶Boussinesq方程,在交错网络下建立数值模型.利用模型模拟波浪在常水深情况下的传播,波浪反射系数均低于2％.利用该模型模拟波浪在平斜坡前的反射,并将数值结果与解析解进行对比.结果表明,对于相对水深较大情况,坡度较陡时模拟结果明显偏大;对 于相对水深较小情况,坡度超过1:1时,数值结果仍与解析解有....     

用活有孔虫分布指示三亚汊道外港潮流格局

由因子分析得知,海南岛三亚湾汊道外港海底分布有四个主因子活有孔虫群,它们分别代表海湾内水体、湖冲出半盐水水体、湾外内陆架水体及内陆架与海湾水相混的过渡水体.这些水体在潮汐作用下的运动、相互混合,可由各水体的有孔虫指示.由每个主因子因子得分最高的三个属种百分含量之和等值线图所指示的外港潮流循环格局,与海洋学研究所得的结果相当一致.利用有孔虫指示潮汐汊道潮流格局尚属首次尝试,这方法对海岸开发分析泥沙运动、污染物排放等有实践意义.     

海南岛博鳌港洪水过程模拟

采用一个包括陆地淡水径流和风暴潮影响的模型计算了海南省博鳌港的洪水水位。模型中的外海潮汐考虑K^1和O^1分潮，风暴增水和河流流量的变化过程采用峰值函数的形式模拟，取得了与实测结果相近的计算结果。模拟实验表明，洪水水位峰值与河流洪峰流量，风暴增水峰值，系统的过水断面形状以及洪水历时的长短等因素有关。对于博鳌港口门通航条件整治，采用复合的梯形过水断面最为合适，因为相对于天然的口门过水断面，复合梯形断面能使湾内的最高洪水水位下降约20％。     

锦州港潮流悬沙淤积模型试验研究

本文首先分析了锦州湾的自然条件，明确了风浪掀沙和潮流输沙是引起港池和航道淤积的主要原因并得到了海水含沙量与风、浪、潮的定量关系，为确定模型中应施放的特征含沙值提供了依据。文中首次将窦国仁全沙模型相似理论应用于海湾潮流泥沙模型并取得了成功。验证试验表明，模型中潮位、潮流速、流向、含沙最分布，淤积强度分布，淤积量和淤沙粒径级配等均与原型基本一致。利用此泥沙模型对锦州港各期工程方案的淤积量和淤积强度分布进行了预报并对今后港口扩建工程布置提出了建议。