双层局部开孔板沉箱对波浪反射的理论研究

提出了一种用于研究由双层开孔板和一个不透水后板的开孔结构对斜向波反射率的理论分析方法。整个流域被分成三个子域，在每个子域内应用特征函数展开法以得到该域内包含未知展开系数的势函数的表达式，在速度势的展开中，考虑了非传播模态波浪的影响。通过匹配开孔板处的边界条件可以求解待定的展开系数，继而求解双层开孔板防波堤结构对斜向波的反射率。数值计算结果与试验结果进行了比较，符合较好。并进一步讨论了几个重要因素对反射系数的影响。     

A polar method for obtaining wave resonating quadruplets in finite depths

A polar method for obtaining wave resonating quadruplets {K₁, K₂, K₃, K₄} in the computation of nonlinear wave–wave interaction source term of the wave model is presented with results for both deep and finite water depths. The method first determines the end radial points of the locus equation for K₂, for each set of input wave vectors (K₁, K₃) on the symmetry. The locus of K₂ (and hence K₄) is then traced in the anti-clockwise direction starting with the maximum radial point on the line of symmetry. It is shown that when k₃>k₁, the number of points on the locus varies when the orientations of the input wave vectors are changed and reduces when the difference in the magnitude of the input wave vectors is increased. A significant advantage in this method is that the angular increment on the locus for K₂ can be kept constant.     

An Application of Direct Ritz Method for Random Seismic Responses of Jacket Platforms

1 .IntroductionNowadays tremendous efforts have been devotedtothe analysis of randomseismic responses .Forexample ,American PetroleumInstitute (API) has declared that the analysis of randomseismic re-sponsesis a useful tool for safety inspection. Owing to the complexity of platformand randomness ofload,it is difficult to carry out the randomresponse analysis of a jacket platform. The conventionalrandomvibration methods ,such as square root of the sumof squares (SRSS) and complete quadratic…     

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

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

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.     

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.     

Baroclinic Buoyancy-Inertia Joint Stability Parameter

At present, the barotropic buoyant stability parameter has been derived from a vertical virtual displacement of a water parcel. The barotropic inertial stability parameter in the eccentrically cyclogeostrophic, basic current field was derived in 2003 from a horizontal cross-stream virtual displacement of a parcel. By expressing acceleration of a parcel due to a virtual displacement, which is arbitrarily sloping within a vertical section across the basic current, in terms of natural coordinates, we derived the vertical component of baroclinic buoyant stability parameter B ₂ ², the horizontal component of baroclinic inertial stability parameter I ₂ ², the baroclinic joint stability parameter J ², its buoyant component B ² and its inertial component I ². B ² is far greater than I ₂ ², and when neglecting relative vorticity except for vertical shear, a downward convex curve of J ² plotted against the slope of a virtual displacement follows a trend of B ² curve. If a parcel displaces along a horizontal surface or an isopycnal surface, however, B ² vanishes, and J ² becomes equal to I ². Actual parcel is apt to displace not only along the bottom slope, but also along the sea surface and an isopycnal interfacial surface, which is approximately equivalent to an isentropic surface, preferred by lateral mixing and exchange of momentum. Such actual displacement makes B ² vanishing, and grants I ² an important role. The present analysis of I ² examining effects due to curvature and horizontal and vertical shear vorticities are useful in deepening our understanding of baroclinic instability in actual oceanic streams.