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.     

A New Fall-Rate Equation for T-5 Expendable Bathythermograph (XBT) by TSK

The accuracy of the manufacturer’s fall-rate equation for the T-5 Model of expendable bathythermograph (XBT) has been investigated based on about 300 collocated pairs of XBT-CTD (Conductivity-Temperature-Depth profiler) measurements in various climatological regions. We found that the equation systematically overestimates depth by about 5% for the T-5 produced by Tsurumi Seiki, Co. Ltd. (TSK), but almost no bias is associated with the T-5 produced by Sippican, Inc., in USA. The cause of this difference is not clear, because the two manufacturers’ T-5 probes are reported to have identical shape and weight in water. We propose a new fall-rate equation for the TSK T-5: z(t) = 6.54071t - 0.0018691t ², where z(t) is depth in meters at time, t, in seconds.     

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.     

Geodynamic Evolution of the Eastern Segment of the Azores-Gibraltar Zone: The Gorringe Bank and the Gulf of Cadiz Region

Detailed structural interpretation of the recently acquired deep seismic multichannel profiles along the Iberian Atlantic Margins (IAM Project) provides new results on the geodynamic evolution of the eastern part of the Azores-Gibraltar plate boundary. Thrusting and folding of the oceanic basement and of Mesozoic and Cenozoic sedimentary cover of the Gorringe Bank region are consistent with the N–S convergence of Iberia and Africa. Compressive structures in the Gorringe Bank region are spread over a wide area. Deformation under compression took place mainly in Tertiary times, as is evidenced by a basal unconformity and several discontinuities in Tertiary sediments, although some deformation has also been recorded in Quaternary sediments. The compressive structures in the Gulf of Cadiz are E–W oriented thrusts, folds and related diapiric structures. N–S oriented transpressive deformation is likely to occur in the western Portuguese platform. There is no continuity of structures from the oceanic to the continental domain, suggesting that deformation transfers from one side to the other through a transcurrent fault zone. The fault contact between the two domains is located in the ocean-continent transition zone.     

一维扩散方程整体观测资料下的初值变分同化反演

本文考虑一维扩散方程的反问题,利用变分同化方法通过观测资料来确定方程中的未知初值,通过分析观测误差对于初值误差的影响,证明变分同化初值收敛于原问题的真实参数,并得到了参数的收敛精度。同时将得到的初值代入预报模式中,得到预报解,并分析了预报解的收敛性和预报误差。     

Modelling linear wave transformation induced by dissipative structures—Random waves

The relevant theory is presented and numerical results are compared with the analytical solution for the interaction of non-breaking waves with an array of vertical porous circular cylinders on a horizontal bed. The extension to the cases of unidirectional and multidirectional waves is obtained by means of a transfer function. The influence of the mechanical properties of porous structures and wave irregularity on wave transformation is analysed. Results for unidirectional and multidirectional wave spectra are compared to those obtained for regular waves. The model presented reproduces well the analytical results and provides a tool for analysing several engineering problems.     

赤道西太平洋大气边界层的诊断分析

利用１９９３年ＥＮＳＯ事件爆发（４月）前酝酿时期“热带大洋与全球大气－海洋耦合响应试验”强化观测阶段“向阳红五号”科学考察船１５５°Ｅ，２°Ｓ定点海洋气象和高空大气探测资料，分析了赤道西太平洋大气边界层特征。结果表明：赤道西太平洋上空辐合对流区边界层内也有等温或逆温层存在，大气边界层物理参数变化与大尺度海－气变异有密切关系，赤道太平洋沃克环流加强、东移，边界层上部逆温层出现概率明显增大，边界层内高湿层湿度减小。另外，边界层内实测风速、风向廓线随高度的变化基本上不服从Ｅｋｍａｎ规律，但由其平均，分量合成得到的平均风随高度的变化却基本符合Ｅｋｍａｎ规律。     

Study on the characteristics of the marine boundary layer in the Equatorial Pacific

ＳｔｕｄｙｏｎｔｈｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｔｈｅｍａｒｉｎｅｂｏｕｎｄａｒｙｌａｙｅｒｉｎｔｈｅＥｑｕａｔｏｒｉａｌＰａｃｉｆｉｃ￥ＺｈａｎｇＺｉｙｕａｎｄＺｈｏｕＭｉｎｇｙｕ（ＲｅｃｅｉｖｅｄＡｕｇｕｓｔ２１，１９９３；ａｃｃｅｐｔｅｄ...     

Second-Order Wave Diffraction Around 3-D Bodies by A Time-Domain Method

A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.