均匀土质边坡稳定性计算的一种新方法

根据土力学原理 ,推导出一种在不考虑渗流力的情况下 ,确定斜坡临界滑动面和计算斜坡稳定系数的新方法。介绍了该方法在若干土石方工程中的应用实例。     

西太平洋暖池区海—气通量计算分析

用J.Launiaimen和T.Vihma提出的近地面层湍流通量计算方法,对我国在1992年11月至1993年2月TOGA—COARE—IOP实验中所获资料计算处理。得出所在站位的海一气间显效、潜热及动量通量。指出西大平洋暖池海区游热通量与显效通量之比为10.14:1;风速大于8m/s后各通量随风速的变化率明显增加;动量与热量的块体通量系数Cd和Ce,h随风速变化有相似的规律;Monin—Obukhov大气稳定度参数Z/L与△T/U_(10)之间有较好的统计关系。     

Validation and Improvement of Satellite-Derived Surface Solar Radiation over the Northwestern Pacific Ocean

The purpose of this study is to validate and improve satellite-derived downward surface shortwave radiation (DSSR) over the northwestern Pacific Ocean using abundant in situ data. The DSSR derivation model used here assumes that the reduction of solar radiation by clouds is proportional to the product of satellite-measured albedo and a cloud attenuation coefficient. DSSR is calculated from Geostationary Meteorological Satellite-5/Visible Infrared Spin-Scan Radiometer data in 0.05° × 0.05° grids. The authors first compare the satellite DSSR derived with a cloud attenuation coefficient table determined in past research with in situ values. Although the hourly satellite DSSR agrees well with land in situ values in Japan, it has a bias of +13∼+34 W/m² over the ocean and the bias is especially large in the low latitudes. The authors then improve the coefficient table using the ocean in situ data. Usage of the new table successfully reduces the bias of the satellite DSSR over the ocean. The cloud attenuation coefficient for low-albedo cases over the ocean needs to be larger in the low latitudes than past research has indicated. Daily and hourly DSSR can be evaluated from the satellite data with RMS errors of 11–14% and 30–33%, respectively, over a wide region of the ocean by this model. It is also shown that the cloud attenuation coefficient over land needs to be smaller than over the ocean because the effect of the radiation reflected by the land surface cannot be ignored.     

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.     

两参量的海面阻力系数模式的探讨

从风浪的能量平衡方程出发，引进若干风要素与波要素以及波要素之间的定性关系，经演算可导出海面阻力系数（Ｃｐ）或是风速（Ｕ）和波龄（β）或是Ｕ和波高（Ｈ）的函数，然后沿用最小二乘法，终将得出４组１２个回归方程。当β（或β）或Ｈ为某一给定值，惟有Ｕ为唯一参量时，所提各式均可简化为非线性方程：ＣＤ＝ａ＋ｂ，Ｕ＋ｃ．Ｕ＾２；式中ａ，ｂ和ｃ为三个经验系数，就所检验的例子而言，本文的结果与实际的符合前人的为好。     

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.