聚类分析法在浅海水团分析中的应用及黄、东海变性水团的分析

本文引入变性水团概念,用聚类分析确定黄、东海浅海海域变性水团的边界及其变化。根据给定的标准进行聚类,确定了水团的个数、边界及混合区。 由聚类方法所得结果看出,在该海区中有八个变性水团。根据温、盐相对指标,将这些水团分为九种不同的特征。也可以把这些水团划分为三种盐度类型。在温—盐点聚图上,不同变性水团的温盐点,分布于一条曲线附近,它体现出逐级变性的特点。各变性水团的变化范围都很大,这说明该海区中水团变性强烈。在暖季中的增温降盐和冷季中的降温增盐,可认为是整个黄海和东海水体变性的特征。所有季节性水团都经历一个形成与消亡的过程。 文中讨论了变性水团和海洋环境的相互关系。水团变性是热力因素和动力因素共同作用而产生的。通过变性水团的分析,作者认为,在该海区水团的变性中,热力因素起着主要作用。变性水团界线的舌状分布与流向之间的关系是很明显的,可作为分析环流的旁证。最后描述了变性水团和底层中心渔场之间的关系。     

A finite element model of interaction between viscous free surface waves and submerged cylinders

This paper describes the simulation of the flow of a viscous incompressible Newtonian liquid with a free surface. The Navier–Stokes equations are formulated using a streamline upwind Petrov–Galerkin scheme, and solved on a Q-tree-based finite element mesh that adapts to the moving free surface of the liquid. Special attention is given to fitting the mesh correctly to the free surface and solid wall boundaries. Fully non-linear free surface boundary conditions are implemented. Test cases include sloshing free surface motions in a rectangular tank and progressive waves over submerged cylinders.     

A dynamic response analysis of tension leg platforms including hydrodynamic interaction in regular waves

A numerical procedure is described for predicting the motion and structural responses of tension leg platforms (TLPs) in waves. The developed numerical approach, in a TLP is assumed to be flexible instead of rigid, is based on a combination of the three dimensional source distribution method and the finite-element method. The hydrodynamic interactions among TLP members, such as columns and pontoons, are included in the motion and structural response analysis. Numerical results are compared with the experimental and numerical ones. The results of comparison confirmed the validity of the proposed approach.     

A B-spline based higher order panel method for analysis of steady flow around marine propellers

A higher order panel method based on B-spline representation for both the geometry and the solution is developed for the analysis of steady flow around marine propellers. The self-influence functions due to the normal dipole and the source are desingularized through the quadratic transformation, and then shown to be evaluated using conventional numerical quadrature. By selecting a proper order for numerical quadrature, the accuracy of the present method can be increased to the machine limit. The far- and near-field influences are shown to be evaluated based on the same far-field approximation, but the near-field solution requires subdividing the panels into smaller subpanels continuously, which can be effectively implemented due to the B-spline representation of the geometry. A null pressure jump Kutta condition at the trailing edge is found to be effective in stabilizing the solution process and in predicting the correct solution. Numerical experiments indicate that the present method is robust and predicts the pressure distribution on the blade surface, including very close to the tip and trailing edge regions, with far fewer panels than existing low-order panel methods.     

长江口及济州岛附近海域变性水团的初步分析

基于变性水团的概念,本文把聚类分析用于确定长江口及济州岛附近海域变性水团的边界。由聚类分析得到的结果表明,在该海区有十个水团。对它们的特征、分布与变化进行了初步分析。作者得出的结论是:1、在十个水团中,有四个大洋性水团,六个变性水团。2、该海域水团变性的特点为暧季增温、降盐、降氧,而冷季则相反。3、水团变性是由海区内、外因素综合作用而发生的,而后者在浅水区域起主要作用。4、水团边界的舌状分布与流向之间有明显关系。因之,海流的方向及强度,大致可依水团舌状分布而判断。5、底层中心渔场基本上位于各变性水团之间的混合区或其附近。     

A Reconstruction of Observed Profiles in the Sea East of Japan Using Vertical Coupled Temperature-Salinity EOF Modes

It is important to estimate hard-to-observe parameters in the ocean interior from easy-to-observe parameters. This study therefore demostrates a reconstruction of observed temperature and salinity profiles of the sea east of Japan (30°≈40°N, 140°≈150°E). The reconstruction was done by estimating suboptimal state from several values of the observed profiles and/or sea surface dynamic height (SDH) calculated from the profiles. The estimation used a variational method with vertical coupled temperature-salinity empirical orthogonal function (EOF) modes. Profiles of temperature and salinity in the subtropical region are effectively reconstructed from in situ temperature profile data, or sea surface temperature (SST) and SDH. For example, the analyzed temperature field from SST and SDH has an accuracy to within 1°C in the subtropical region. Salinity in the sea north of Kuroshio, however, is difficult to estimate because of its complex variability which is less correlated with temperature than in the subtropical region. Sea surface salinity is useful to estimate the subsurface structure. We also show the possibility that the estimation is improved by considering nonlinearity in the equation calculating SDH from temperature and salinity analysis values in order to examine the misfit between analysis and observation. Analysis using TOPEX/POSEIDON altimetry data instead of SDH was also performed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

P矢量方法在南海夏季环流诊断计算中的应用

基于1998年6～7月南海调查航次的CTD资料,对南海环流采用最近发展的P矢量方法进行诊断计算.计算结果:黑潮向西入侵南海,然后做反气旋弯曲向东北方向流动,最终有通过巴士海峡流出南海的趋势.在南海北部存在一个气旋性环流,这个环流的强度和范围随深度增加而减小.该环流的冷中心位置随深度增加稍向南移.南海中部、越南以东海域存在一个明显的气旋涡和反气旋涡,尤其在200m及其以上水层均相当稳定,反气旋涡位于越南以东,其中心位置在11°53'N,111°50'E,气旋涡的中心位置在13°17'N,112°55'E,两者的尺度皆约为250km.吕宋岛西侧存在一个反气旋涡.在计算海区南部、巴拉望岛西南海域,100m以上层存在一个反气旋式涡.从各层流场分布均可以显示海流在西部强化的现象.     

Modified mild slope equation and open boundary conditions

A numerical model to compute wave field is developed. It is based on the Berkhoff diffraction-refraction equation, in which an energy dissipation term is added, to take into account the breaking and the bottom friction phenomena. The energy dissipation function, by breaking and by bottom friction, is introduced in the Berkhoff equation to obtain a new equation of propagation.The resolution is done with the hybrid finite element method, where lagrangians elements are used.     

Analysis of seabed instability using element free Galerkin method

Wave-induced seabed instability, either momentary liquefaction or shear failure, is an important topic in ocean and coastal engineering. Many factors, such as seabed properties and wave parameters, affect the seabed instability. A non-dimensional parameter is proposed in this paper to evaluate the occurrence of momentary liquefaction. This parameter includes the properties of the soil and the wave. The determination of the wave-induced liquefaction depth is also suggested based on this non-dimensional parameter. As an example, a two-dimensional seabed with finite thickness is numerically treated with the EFGM meshless method developed early for wave-induced seabed responses. Parametric study is carried out to investigate the effect of wavelength, compressibility of pore fluid, permeability and stiffness of porous media, and variable stiffness with depth on the seabed response with three criteria for liquefaction. It is found that this non-dimensional parameter is a good index for identifying the momentary liquefaction qualitatively, and the criterion of liquefaction with seepage force can be used to predict the deepest liquefaction depth.