波浪运动在底边界层的湍流结构数值研究

本文基于$ k $-$ \varepsilon $ 模型研究了波流边界层内湍流结构特征。研究结果表明,时均流速分布数值解与实验结果高度吻合。一个波周期内湍流结构特征（如：涡量、湍动能、湍动能耗散率等）呈周期性变化规律,波浪作用引起涡量、湍动能及湍动能耗散率均在减速阶段减小,在波谷处达到最低值,而后在加速阶段增大,并在波峰处达到最大值。近壁面处湍流结构变化幅值较大（湍动能耗散率变化可达53%）,远离壁面处变化幅值较平均值较小（仅3%）。波流边界层厚度在减速阶段增加,在加速阶段减小。本文所建立的数值模型克服了现有模型因采用“高雷诺数方法”引起的近壁区精度不高问题,可较好地描述波浪作用下湍流结构演变过程的物理机制,为河口海岸地区泥沙运动、岸滩演变及海洋可再生能源的开发利用提供一些指导意义。     

Research on the Effects of in-Line Oscillatory Flow on the Vortex-Induced Motions of A Deep Draft Semi-Submersible in Currents

A Deep Draft Semi-submersible (DDS) under certain flow conditions could be subjected to Vortex-Induced Motions (VIM), which significantly influences the loads on and life fatigue of the moorings and the risers. To investigate the VIM of a DDS with four rectangular section columns in waves coupled with a uniform current, a numerical study using the computational fluid dynamics (CFD) method was conducted. The issues of the VIM of multi-column floaters can be conveniently converted to the issues of oscillating cylinders in fluid cross flows. This paper looks into the CFD numerical simulation of infinite cylinders having rectangular sections in a two-dimensional sinusoidal time-dependent flow field coupled with a uniform current. The resulted hydrodynamic forces and motion responses in different oscillatory flows plus currents both aligned in the same direction for the incidence of 135° of the DDS relative to the flow are compared with the ones in current only cases. The results show that the VIM response of this geometric arrangement of a DDS with four rectangular columns in a current combined with oscillatory flows is more evident than that in the current only case. The oscillatory flows and waves have the significant influence on the VIM response, forces and trajectory, in-plane motions of the DDS.     

近地层风廓线与垂直湍流强度关系研究

通过改进Prandtl 混合长理论,得出近地层风廓线形式     

两种湍流模型在潮流能水轮机数值模拟中的适用性研究

在对特定湍流条件下的水平轴潮流能水轮机进行数值模拟时,不同的湍流模型对于计算结果具有一定的影响。为选择合适的湍流模型对所设计水轮机的三维流场进行描述,将采用常用的两方程模型(SSTk-ω与Realizable k-ε)进行数值模拟,得到不同工况下水轮机的性能参数以及水轮机后1.5D处流场的流速、湍流强度分布情况,并进行模型水槽试验,试验值与采用两种湍流模型进行数值模拟得到的仿真值对比表明:在低尖速比工况下,采用SSTk-ω进行数值模拟得到的仿真值比Realizable k-ε更接近试验值,在较高尖速比时,采用Realizable k-ε进行模拟得到的数值更接近试验值;采用两种湍流模型进行数值模拟得到的速度分布具有相似性,但对于湍流强度的预测,SSTk-ω更具有优势。考虑此水轮机的实际工况,选择SSTk-ω湍流模型进行描述更加合理。     

A Study of the Performances of Different Turbulence Schemes in Numerical Simulation of Hydrodynamics of a Semi-Closed Sea (Persian Gulf)

Complex process of turbulent mixing in Persian Gulf that is a semi-closed sea makes it a good media to test the performance of different turbulence schemes. In this research, we used the 3D ocean model COHERENS (COupled Hydrodynamical Ecological model for REgioNal Shelf seas) for the Persian Gulf with the open boundary in the Hormuz Strait. Of the turbulence schemes for the vertical diffusion available in the COHERENS, we tested four models to investigate the hydrodynamic characteristics of the Persian Gulf. The results show that all of the schemes presented the sea surface salinity (SSS) distribution rather accurately but the k-l and flow-dependent models results have better agreements with observations. The most noticeable difference between the results of four schemes is the differences found in the simulation of turbulent parameters. The turbulent closure schemes generally provide better results, but the algebraic schemes show turbulent parameters far from reality and they do not show substantial changes with time. Generally, the vertical structures of turbulence in the water basins and parameterization of turbulence in water column is very sensitive to the selection of the type of the turbulence scheme. However, large-scale structures that take place within the inflow and outflow area are approximately quasi-horizontal, and the vertical small-scale turbulence does not affect them as much. As a result, they show less sensitivity to the performance of various turbulence schemes.     

珠江河口底边界层湍流积分尺度研究

介绍湍流积分尺度的几种常用计算方法.基于三组珠江口崖门的湍流观测资料,分别计算了一个潮周期内珠江口底边界层湍流积分尺度的变化,并对比不同计算方法的稳定性.计算结果表明,基于Taylor假定的自相关函数积分法的计算结果稳定可靠.湍流积分尺度在一个潮周期内的变化很大,积分尺度的变化与平均流的流速变化有着很好的对应关系,三个测次纵向、横向和垂向积分尺度的平均值分别约为4.9～10.4 cm、3.3～6.2 cm和4.3～5.2 cm.     

Turbulent and Macro-turbulent Structures Developed in the Benthic Boundary Layer Downstream of Topographic Features

The characteristics and effects of large-scale flow structures developed in the benthic boundary layer downstream from large topographic features were analysed throughout a tidal cycle. The observed signature of the macro-turbulent features consisted of streamwise modules of low horizontal velocity and high suspended sediment concentration (SSC), alternating with modules of high horizontal velocity and low SSC. These modules extended 10 to 20 m streamwise and exceeded 1 m vertically, and are believed to be related to flow separation effects over large bedforms upstream of the deployment site. The macroscale flow modules intensified the ‘ burst-like ’ turbulent events and favoured sediment transport. ‘ Ejection-like ’ events were magnified during modules of decreasing horizontal velocity and increasing turbidity, whereas ‘ sweep-like ’ events were magnified during modules of increasing horizontal velocity and decreasing SSC. The enhanced turbidity of the macroscale modules may be the result of enhanced upward diffusion of sediment by ejection events, whereas the low-turbidity modules may be induced by increased downward transport of suspended sediment by sweep events. These hypotheses were supported by cross-spectral analysis performed on velocity and suspended sediment concentration time-series recorded at the site. An enhanced (negative) contribution of outward and inward interaction events to the Reynolds stress, compared to those reported in uniform BBLs, resulted in ‘ abnormally ’ low stress values.     

黄海、渤海底层温度数值预报方法的研究

根据文献［ｌ］建立的底层温度（ＴＨ）与其水柱垂向平均温度（）的经验关系，结合流体动力学方程和（垂向平均）热传导方程，发展了以水气温差和风速为已知量的底层温度二维数值预报模式。该模式避开了海面热量和动量输入在垂直水柱中分配的复杂物理过程而直接报出底层水温场，具有较好的实用性；此外，从试报结果看，效果令人满意。     

A Study on Residual Flow in the Gulf of Tongking

On the basis of observations of water temperature, salinity and wind during winter (December to following February) and summer (June to August) in the Gulf of Tongking, a robust diagnostic three-dimensional model has established that reveals the seasonal variation in residual flow, including wind-driven current, density-driven current and tide-induced residual current. It is shown that in the Gulf of Tongking the wind-driven current plays the most important role in the seasonal variation of residual flow field. Due to strong NE monsoon (9 m/s) during winter the residual flow (reaches 30 cm/s) is stronger than that in summer. At the surface the residual flow direction during winter is basically westward or southwestward, while during summer it is almost reversed.     

Numerical Study on Instantaneous Discharge of Unsorted Particle Cloud in Cross Flow

The mixing characteristics of particles such as dredged sediment of variable size discharged into cross flow are studied by a 3D numerical model, which is developed to model the particle-fluid two-phase flow. The Eulerian method with the modified k-ε parameterization of turbulence for the fluid phase is used to solve fluid phase, while a Lagrangian method for the solid phase (particles), both the processes are coupled through the momentum sources. In the model the wake turbulence induced by particles has been included as additional source term in the k-ε model; and the variable drift velocities of the particles are treated efficiently by the Lagrangian method in which the particles are tracked explicitly and the diffusion process is approximated by a random walk model. The hydrodynamic behavior of dumping a cloud of particles is governed by the total buoyancy of the cloud, the drag force on each particle and the velocity of cross-flow. The computed results show a roughly linear relationship between the displacement of the frontal position and the longitudinal width of the particle cloud. The particle size in the cloud and the velocity of cross flow dominate the flow behavior. The computed results are compared with the results of laboratory experiments and satisfactory agreement is obtained.     

黄海夏季海雾的边界层结构特征及其与春季海雾的对比

利用海上浮标站、高分辨率数字式探空仪等多种观测手段和三维中尺度模式,对2008年7月7~11日1次夏季黄海海雾过程的边界层结构特征进行了观测分析与数值模拟,并将结果与春季的黄海海雾个例进行对比。结果表明:(1)夏季海洋大气边界层(MABL)中无强逆温层,静力稳定度较春季下降,有利于湍流的发展。加之水汽量较大,容易形成比较厚的雾(500 m);春季低空有明显逆温层,水汽供应量较少,但强稳定的层结可以使水汽局限于比较低的空中,形成比较薄的雾层(200 m)。(2)在夏季风控制下,青岛近海海洋大气边界层(MABL)中高、低层气块均来自海洋上空,温湿属性差异不大,使得温度垂直差异较小;春季MABL中高、低层气块分别来自陆地和海洋,来自陆地的暖、干气流和海洋的冷却效应导致强逆温层和雾区上方干层的出现。(3)夏季海雾的含水量大,水汽在凝结成雾的过程中放出更多的凝结潜热,雾中的海表面气温(SAT)明显高于海表面水温(SST);春季雾顶强烈的长波辐射冷却和湍流混合使雾中气温明显下降,雾中SAT与SST更加接近,甚至出现SAT小于SST。以上结果有助于对海雾形成机制的认识。     

Numerical Study on the 3-D Complex Characteristics of Flow Around the Hull Structure of TLP

Vortex-induced motion is based on the complex characteristics of the flow around the tension leg platform (TLP) hull. By considering the flow field of the South China Sea and the configuration of the platform, three typical flow velocities and three flow directions are chosen to study the numerical simulation of the flow field characteristics around the TLP hull. Reynolds-averaged Navier–Stokes equations combined with the detached eddy simulation turbulence model are employed in the numerical study. The hydrodynamic coefficients of columns and pontoons, the total drag and lift coefficients of the TLP, the formation and development of the wake, and the vorticity iso-surfaces for different inlet velocities and current directions are discussed in this paper. The average value of the drag coefficient of the upstream columns is considerably larger than that of the downstream columns in the inlet direction of 0°. Although the time history of the lift coefficient demonstrates a “beating” behavior, the plot shows regularity in general. The Strouhal number decreases as the inlet velocity increases from the power spectral density plot at different flow velocities. The mean root values of the lift and drag coefficients of the front column decrease as the current direction increases. Under the symmetrical configuration of 45°, the streamwise force on C4 is the smallest, whereas the transverse force is the largest. The broken vortex conditions in current directions of 22.5° and 45° are more serious than that in the current direction of 0°. In addition, turbulence at the bottom of the TLP becomes stronger when the current direction changes from 0° to 45°. However, a high inlet velocity indicates a large region influenced by the broken vortex and shows the emergence of the wake behind the TLP under the same current angle.     

瓯江口外洞头岛北岙后涂围垦工程潮流数值模拟研究

通过建立基于不规则三角形差分法的洞头岛附近海区二维潮流数学模型,在对现场实测资料验证的基础上,对北岙后涂围垦工程前后的潮流场进行了数值模拟,就北岙后涂围垦工程对周边海区的影响及围塘根部的冲刷情况进行了分析研究。分析研究结果表明,围垦工程不会对客运和货运码头及航道通航条件产生明显影响,对洞头峡的涨落潮水流、三盘大桥桥墩及其通航孔基本没有影响;单纯潮流对围塘根部冲刷的可能性很小,因此,围垦工程是可行的。     

A Numerical Modeling of the Upper and the Intermediate Layer Circulation in the East Sea

Circulation in the upper and the intermediate layer of the East Sea is investigated by using a fine resolution, ocean general circulation model. Proper separation of the East Korean Warm Current from the coast is achieved by adopting the isopycnal mixing, and using the observed heat flux (Hirose et al., 1996) and the realistic wind stress (Na et al., 1992). The simulated surface circulation exhibits a remarkable seasonal variation in the flow patterns of the Nearshore Branch, the East Korean Warm Current and the Cold Currents. East of the Oki Bank, the Nearshore Branch follows the isobath of shelf topography from late winter to spring, while in summer and autumn it meanders offshore. The Nearshore Branch is accompanied by cyclonic and anticyclonic eddies in a fully developed meandering phase. The meandering and the eddy formation of the Nearshore Branch control the interior circulation in the Tsushima Current area. A recirculation gyre is developed in the region of the East Korean Warm Current in spring and grown up to an Ulleung Basin scale in summer. A subsurface water is mixed with the fresh surface water by winter convection in the northeastern coastal region of Korea. The well-mixed low salinity water is transported to the south by the Cold Currents, forming the salinity minimum layer (Intermediate Water) beneath the East Korean Warm Current water. The recirculation gyre redistributes the core water of the salinity minimum layer in the Ulleung Basin. This revised version was published online in August 2006 with corrections to the Cover Date.     

春季南海南部上混合层数值模拟与数值实验

采用一维湍动能模式对南海南部的 SST及混合层进行数值模拟和数值试验。结果表明 :TKE模式能够模拟南海南部的海表面温度 SST以及除南海南部 5月中旬以外的上混合层深度随时间变化基本特征。在 5～ 6月 ,SST的日振荡主要依赖于短波辐射的日变化 ,风的混合作用抑制了 SST的日周期振荡。春季夏季风爆发期间 ,南海海面潜热通量和感热通量与短波辐射和风应力相比较 ,是一个对 SST和混合层影响较小的量。在春季南海南部 ,短波辐射作用能使 SST升高的最大值约为 4℃ ;潜热和感热通量能使 SST的下降的最大值为 3℃。风应力对南海混合层深度随时间变化趋势起着决定的作用 ,并能使其深度加深 2 0～ 30 m,而短波辐射则使混合层的深度变浅2～ 3m,潜热和感热通量会使混合层的深度加深 1～ 2 m。在春季南海南部 ,热通量对混合层深度的影响与风应力相比要小得多     

A Numerical Simulation on the Mesoscale Dynamics of the Spring Bloom in the Sea of Japan

The northward migration of spring bloom was observed in the Sea of Japan from April to May 1997 by the Ocean Color and Temperature Scanner (OCTS) on board the Advanced Earth Observing Satellite (ADEOS). This phenomenon is well simulated with a numerical ecosystem model coupled with a hydrodynamic model. The hydrodynamic model is the Geophysical Fluid Dynamics Laboratory (GFDL) Modular Ocean Model (MOM). The ecosystem model consists of five components: dissolved inorganic nitrogen (DIN), dissolved inorganic phosphorus (DIP), phytoplankton, zooplankton and detritus. Results of the numerical ecosystem model suggest that the mesoscale development of the spring bloom in the Sea of Japan is related to that of sea water temperature, and that the bloom is limited by the depletion of DIN. This revised version was published online in August 2006 with corrections to the Cover Date.     

内蒙孪井灌区地下水数值模拟及土壤盐渍化分析

系统分析了内蒙孪井灌区的水文地质条件和包气带水文地质特征，概化出该区水文地质概念模型。在典型地段对不同作物进行了大量的野外灌溉入渗试验，确定了灌溉渗漏量；通过抽水试验，确定了该区的水文地质参数；利用长观资料。建立了该区地下水的数值模拟模型并利用该模型预报了现状、规划和下游排水3种条件下地下水的变化趋势。研究表明：该区随灌溉面积的扩大，地下水位上升迅速，灌区下游水位将超过临界埋深，从而引起土壤次生盐渍化。     

A Preliminary Study on the Calculating Method for the Evolution of Siltation and Erosion of Tidal Flat Based on GIS

The estuary and coast is an area where the land and the sea interact and a place in which human beings frequently move about so that understanding and controlling the change and development modes of the coastal landform plays a vital part in exploiting and protecting coastal resources. A model is the generalization and abstraction of objective things. This paper summarizes four methods for the landform development of the tidal shore and underwater delta, mainly discusses the model‘s structural elements, and presents their specific application on the basis of the authors‘ case study. With the application of the profile model, the dynamic change of coastal landform can be clearly seen by contrasting the different profiles of different years. Through the shrinking, expanding and transformation of the isobath, plane model is used to study the macro-change of the shoal and the coastal landform. Speed model is an efficient means to analyze the trend of erosion and deposition and the local change in a great area of the sea. Statistical survey model is a static analysis, which can be used to establish the relationship between the erosion and deposition of the shoal and the altitude and slope of survey spot.     

Numerical Simulation and Experimental Investigation on Performance of the Wells Turbine in Irregular Oscillating Flow

The Wells turbine is an axial-flow air-turbine designed to extract energy from ocean waves. An important consideration is the self-starting capability of the Wells turbine, a phenomenon encountered where the turbine accelerate by itself up to a certain speed for the best turbine performance. In order to clarify the self-starting characteristic and running performance of the Wells turbine in an irregular oscillating flow, a numerical simulation process is established in this paper on the rational assumption of quasi-steady flow conditions. Both self-starting characteristics and running performance are obtained through the numerical simulation and subsequently compared with the experimental data achieved on a computer-controlled oscillating flow test rig which could realize some irregular oscillating flow according to the specified spectrum. Results show that the self-starting time decreases with the increase of the significant wave height and the mean frequency of the irregular oscillating flow. Therefor     

江苏如东西太阳沙及烂沙洋海域潮流泥沙数值模拟

基于不规则三角形网格有限差分法并考虑波浪及其破碎作用,建立了平面二维潮流场和泥沙场数学模型.该模型对有望建设成深水码头和深水航道的江苏如东西太阳沙和烂沙洋海域的潮流场和泥沙场进行了细化数值模拟.数值模拟结果表明：（1）本海区潮流基本上是顺深槽流动的往复流,潮流流速大,烂沙洋北水道和西太阳沙附近大潮涨落潮最大流速分别在2 m/s和1 m/s以上;（2）本海区的潮平均水体含沙量在0.5 kg/m^3以下,落潮含沙量大于涨潮含沙量;（3）小浪对水体含沙量影响很小,大浪作用下水体含沙量明显增加.