Three-dimensional Layer-integrated Modelling of Estuarine Flows with Flooding and Drying

Details are given of the refinement and application of a thee-dimensional (3-D) layer-integrated numerical model of tidal circulation, with the aim of simulating severe tidal conditions for practical engineering applications. The mode splitting strategy has been used in the model. A set of depth-integrated 2-D equations are first solved to give the pressure gradient, and the layer-integrated 3-D equations are then solved to obtain the vertical distributions of the flow velocities. Attention has been given to maintaining consistency of the physical quantities derived from the 2-D and 3-D equations. A TWO=layer mixing length turbulence model for the vertical shear stress distribution has been included in the model. Emphasis has been focused on applying the model to a real estuary, which is geometrically complicated and has large tidal ranges giving rise to extensive flooding and drying. The model has been applied to three examples, including: wind-driven flow in a rectangular lake, tidal circulation in a model rectangular harbour, and tidal circulation in a large estuary. Favourable results have been obtained for both the simple and complex flow beds.     

分潮波模型多资料数值计算方法研究 Ⅱ.有限差分最小二乘方法

对含多资料的分潮波模型的超定问题进行探讨，建立一种有限差分—最小二乘方法，该方法使得岸边资料直接进入了计算方程，并保证岸界法向流速为零条件。通过对渤海的四大分潮的试算，证实本方法可行。     

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.     

黄、东海环流的数值研究Ⅲ——正压环流的数值模拟

依据黄、东海环流的的动力学模型 ,运用“流速分解法”对黄、东海正压环流进行了数值模拟。计算结果表明冬季黄海正压环流主要受风应力影响 ,基本形态为黄海暖流由济州岛西南进入南黄海中部 ,其东西两侧分别为两支向南流动的沿岸流 ;夏季主要受到潮致体力的影响 ,为一逆时针涡旋。东海环流主要是边界力作用驱动的结果 ,东海黑潮、台湾暖流和对马暖流较稳定。冬季风应力对东海环流表层流场有消弱作用 ,在夏季则有一定增强作用。     

0010号台风"碧利斯"对厦门港湾表层海水中溶解无机砷含量的影响

本文研究了0010号台风“碧利斯”影响期间及其随后的恢复期中,厦门港湾中九龙江口和西港海域站表层海水总溶解砷（TDIAs)、As^3 含量和As^3 /As^5 值的变化特征。结果表明,观测期间两站TDIAs和As^5 含量均比非台风时期呈数倍增加,经10d后仍未恢复到非台风时的水平。As^3 含量也有较大幅度的增加,但波动较大。文中还对观测期间影响溶解无机砷含量的因素进行了初步的探讨。     

New Numerical Scheme for Simulation of Hyperbolic Mild-Slope Equation

The original hyperbolic mild-slope equation can effectively take into account the combined effects of wave shoaling, refraction, diffraction and reflection, but does not consider the nonlinear effect of waves, and the existing numerical schemes for it show some deficiencies. Based on the original hyperbolic mild-slope equation, a nonlinear dispersion relation is introduced in present paper to effectively take the nonlinear effect of waves into account and a new numerical scheme is proposed. The weakly nonlinear dispersion relation and the improved numerical scheme are applied to the simulation of wave transformation over an elliptic shoal. Numerical tests show that the improvement of the numerical scheme makes efficient the solution to the hyperbolic mild-slope equation. A comparison of numerical results with experimental data indicates that the results obtained by use of the new scheme are satisfactory.     

闽江河口潮流数值模拟

应用边界拟合坐标系统，解决了闽江河口河道曲折汊道多造成的计算方向与河底走向交角过大的难点；针对闽江河口长宽比相差较大的特点，移植了通常用于河道计算的沿纵向求解大尺度矩阵问题转化为沿横向小尺度矩阵求逆的水流方程求解方法，采用全隐差分格式，使闽江口数模网格步长最小控制在50m左右，时步长达3min，潮位、流场与实测拟合良好。     

Numerical Simulation of Multi-Directional Random Seas

Multi-directional irregular waves are simulated on the basis of the given directional spectrum using a double summation model, a single direction per frequency model and a single summation model. Their results are compared. It is shown that the single direction per frequency model proposed in this paper can generate a realistic wave field. The effects of the model parameters on the simulated results are also studied in this paper and corresponding suggestions are given.     

分形分析方法在海冰趋势预测中的应用

根据渤海北部1966——2001海冰冰级年际变化的实际资料，利用R／S分析的原理和方法计算了H指数，建立了R(i)/S(i)与i的关系式，对渤海北部海冰趋势进行了预测分析。结果表明，利用分形分析方法进行海冰趋势预测不仅方法简便，而且结论具有客观性，可以将其作为海冰趋势预测中较为有效的一种方法。     

东中国海环流及其季节变化的数值模拟

关于东中国海环流的研究，国内外学者已做了大量的工作。早期科学家们主要依赖于对温盐资料和少数测流资料的分析研究对渤、黄、东海的环流结构有了较系统和深入的认识。东中国海环流是由一个气旋式的“流涡”组成，东侧主要是北上的黑潮-对马暖流-黄海暖流及其延伸部分；西侧为南下的沿岸流系。黑潮对东中国海环流的影响是如此之大，以致于除了某些局部区域外，上述海域主要流系的冬、夏季分布形式比较相似而无本质上的差异(胡敦欣等，1993)。但本文所研究海域正处于世界上最显著的季风区，冬、夏季盛行风向基本相反，过渡季节(春、秋季)风向多变，风力减弱；海洋热盐结构季节变化明显(如冬季混合强，而夏季层化明显等)，这些因素都使得东中国海环流存在着较明显的季节变化。 自20世纪80年代以来，东中国海环流的数值模拟工作逐步展开，并已成为研究环流结构及其形成机制的强有力工具。但由于数值模式本身以及计算方案的缺陷(如有些学者用固定的风场、温盐场对东中国海环流进行诊断模拟等)和观测资料的不足，数值模拟的结果难以得到验证，渤、黄、东海的环流研究中仍有大量的问题存在争议，以待澄清。例如，台湾暖流的来源、流径;对马暖流的来源;夏季黄海暖流的流径以及黄海冷水团环流等均有不同的论述。对黄、东海环流季节变化的数值模拟工作也较少，多用冬、夏典型月份的风场强迫积分至稳定态，给出冬、夏季环流，这种做法值得商榷。三维环流模式很难在1个月内达到稳定态，尤其是夏季层化明显、风力减弱的情况下，非常定风场的影响更应引起人们的重视。 本文采用比较符合实际的计算方案，用年循环风场和海面热通量场为外强迫，对渤、黄、东海的环流及其季节变化进行了模拟，并对一些争议问题进行了探讨。