鲁苏榴辉岩的地球化学

以微量元素、稀土元素、Ｓｒ和Ｎｄ同位素变异特征为依据，确定鲁苏榴辉岩为多成因、多来源和多阶段，指出主要是在印支期扬子陆块与华北陆块碰撞造山作用过程中，挤入的上地幔碎片以及不同原岩类型的壳内高压变质岩碎块。燕山晚期的区域构造热事件使得某些榴辉岩的同位素体系再平衡。     

Strength Analysis of Turret Mooring Productive/Storage Tanker

The longitudinal strength of turret mooring productive/storage tanker is studied.A numeri-cal example has been implemented according to the method presented in this paper to give practical il-lustration.From the results of the numerical example,it is concluded that the turret hole located near theforward of the amidships has small effect on the longitudinal strength of the ship hull.As for design ex-treme value of wave bending moment of storage tanker,statistic method is a more reasonablemethodology,especially with the consideration of the servere environmental conditions.The primary esti-mation of design section modulus of turret storage tanker can be determined by this design bending mo-ment.     

Acoustic plane-wave scattering from a rough surface over a random fluid medium

Acoustic plane-wave scattering from a rough surface overlying a fluid half-space with a sound-speed distribution subject to a small random variation is considered. Under the assumption that the surface roughness and medium randomness are statistically independent, the scattered field may be derived by first solving for the mean field in the medium, and then incorporating with boundary-perturbation method to obtain the total mean field and the power spectral density of the scattered field. The employed algorithm is compatible to the analysis available in the existing literature so that the formulations are conveniently integrated. The results for the power spectral density have shown that the effects of medium inhomogeneities on the rough surface scattering are limited in a spectral regime where the scattered components have shallow grazing angles. The distribution of the power spectral density over the space is primarily governed by the power spectrum and correlation lengths of the rough surface.     

Structural modeling and testing of a concrete canoe

A finite element model is set up and experimental tests are performed to help understand the behavior of a concrete canoe and subsequently optimize its design. First, the performance criteria that must be satisfied to participate at the annual ASCE/Master Builders competition are described. Then, the finite element model and the different loading cases that were studied are presented. Results from these loading cases are discussed and used to optimize the material properties as well as the thickness of the hull and the dimensions of the gunwale and reinforcing ribs. Static and dynamic experimental tests were also conducted to validate the results of the finite element analyses. The results indicate that the main stresses and strains are caused by the static load cases. The additional stresses caused during races are small.     

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.     

Boussinesq-type modelling using an unstructured finite element technique

A model for solving the two-dimensional enhanced Boussinesq equations is presented. The model equations are discretised in space using an unstructured finite element technique. The standard Galerkin method with mixed interpolation is applied. The time discretisation is performed using an explicit three-step Taylor–Galerkin method. The model is extended to the surf and swash zone by inclusion of wave breaking and a moving boundary at the shoreline. Breaking is treated by an existing surface roller model, but a new procedure for the detection of the roller thickness is devised. The model is verified using four test cases and the results are compared with experimental data and results from an existing finite difference Boussinesq model.