复杂网络河的数值模拟

依据复杂网络河的特点及其内在联系，本文建立了其特征矩阵及由Newton-Raphson格式表述的隐式圣－维南差分方程组的系数矩阵，经输入边界条件及解方程组，可得到不同各支段上，下游端点上的水位和流量值。     

Modified PIC Method for Sea Ice Dynamics

The sea ice cover displays various dynamical characteristics such as breakup, rafting, and ridging under external forces. To model the ice dynamic process accurately, the effective numerical modeling method should be established. In this paper, a modified particle-in-cell （PIC） method for sea ice dynamics is developed coupling the finite difference （FD） method and smoothed particle hydrodynamics （SPH）. In this method, the ice cover is first discretized into a series of lagrangian ice particles which have their own sizes, thicknesses, concentrations and velocities. The ice thickness and concentration at Eulerian grid positions are obtained by interpolation with the Gaussian function from their surrounding ice particles. The momentum of ice cover is solved with FD approach to obtain the Eulerian cell velocity, which is used to estimate the ice particle velocity with the Gaussian function also. The thickness and concentration of ice particles are adjnsted with particle mass density and smooth length, which are adjusted with the redistribution of ice particles. With the above modified PIC method, numerical simulations for ice motion in an idealized rectangular basin and the ice dynamics in the Bohai Sea are carried out. These simulations show that this modified PIC method is applicable to sea ice dynamics simulation.     

波峰特征量的统计分布

本文在作者另文数值模拟得到的三维海浪基础上 ,进一步分析给出了三维波峰长度、高度及方向角度等特征量的统计分布。发现考虑海浪的方向性质后波峰分布不再是均匀连续的 ;波峰的长度分布受方向函数影响 ,方向分布越宽 ,波峰的平均长度越短 ,波峰的方向角分布越宽 ;波峰高度和长度在波峰高度较小时有很大的相关性 ,而在波峰高度很大时无关     

粤东海岸台风风暴潮数值预报

粤东海岸台风风暴潮数值预报取风暴潮运行的全流方程组及其有限差分方程组，并取Jelesnianski飓风模型，在注意到复杂海岸和海洋边界的同时，采用独特的边界布局、复合开边界及使用简便的固化预报程序处理方法，提出了一组可适用本海岸的台风风暴潮数值预报模式。在闽粤海岸经袭击影响粤东海岸的18场历史台风风暴潮的数值计算试验，有实测增水资料可对比的试验结果中，多数试验的计算与实测台风增水曲线拟合较佳，达到本项目研究的预期目标。试验结果有力地证明了该风暴潮数值模型是可靠的，可加以推广和应用。     

海南东寨港红树林港湾潮汐动力研究

在海南东寨港林市村附近红树林港湾潮汐动力测量的基础上，建立二维可移动边界有限元数值模型，研究红树林港湾浅水湖滩潮流场特征。实测与计算结果表明，潮沟区流速最大，具驻波特征，呈明显时间－速度不对称，东潮沟涨潮优势，百潮沟落潮优势，整个港湾水体呈顺时针循环。白滩与林地周期性淹没，白滩流速较湖沟略小。林地曼宁系数为白滩的１０倍，流速小于１０ｃｍ·ｓ（－１），约为相应潮沟流速的１／１０。林地如去掉林木，流速可增加３倍以上，流向也有变化。     

洋浦近岸海域污染控制规划方案优化——Ⅰ.排污口位置优选

在洋浦近岸海域计算潮流场基础上,计算预选排污口附近海水质点运动轨迹及预测污染物浓度分布,最后,从环保角度出发,推荐了洋浦地区各开发区的排污口位置,为洋浦近岸海域污染控制规划方案的制定及其优化提供了依据。     

环台湾岛海域半日潮波特征的三维模拟

用1997版POM海洋模式,首次应用于环台湾岛海域的潮波数值研究.得到该海域的半日潮波主要为23°N以南西太平洋传来的胁振潮.影响台湾海峡的半日潮波分别由海峡南北口传入的两支潮波,且北支强于南支.福建沿岸湄州湾-兴化湾为最强潮区,其M₂分潮最大振幅可达240cm.最强潮流区位于澎湖水道,M₂分潮最大潮流达196cm/s.环台湾岛海域潮波潮流水平结构上除海峡北部原有一个圆流点外,还发现另外存在4个新的圆流点.潮流垂直结构上主要为右偏,接近底层处为左偏.     

内波频散关系的一种数值解法

基于Ｆｌｉｇｅｌ等（１９７５）采用的Ｔｈｏｍｐｓｏｎ－Ｈａｓｋｅｌｌ算法（Ｈａｓｋｅｌｌ，１９５１），在取消Ｂｏｕｓｓｉｎｅｓｑ近似的条件下求解线性内波方程，以研究内波的频散关系及其各种模态的垂直速度振幅随深度分布的特征。本文以南海（１３°３０′Ｎ，１１８°３０′Ｅ）、（９°３０′Ｎ，１１１°３０′Ｅ）两个测站夏季和冬季的内波特征的计算为例，分别用Ｆｌｉｇｅｌ等原先的方法（引入Ｂｏｕｓｓｉｎｅｓｑ     

渤海潮汐及大风作用下水位的数值模拟

本文以M₂,S₂,K₁及O₁ 4个主要分潮之和为开边界条件,用“ADI”法对渤海的潮位和潮流进行了数值模拟。潮位和潮流的计算值与潮汐表的预报值(或无风情况下的实测值)比较吻合。在潮汐模拟的基础上,还进行了大风影响下水位场和流场的数值模拟,亦获得了比较满意的结果。