Velocity potential formulations of highly accurate Boussinesq-type models

The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.     

两个同源地不同路径热带气旋的大尺度场对比分析

本文利用NCEP/NCAR再分析格点资料,着重对8903号"BRENDA"和0601号"珍珠"两个同源地进入南海后不同路径的典型热带气旋的大尺度环境场进行对比分析。结果表明:西风槽的强弱和东移的快慢、西太平洋副高的形状变化和东退西进、冷空气的强弱和南下路径及热带气旋内部的不对称结构是造成二者移动路经差异的主要原因。另,印缅槽的变化也对其路径差异有一定的作用。     

Design of ship's bow form by potential-based panel method

A surface panel method treating a boundary-value problem of the Dirichlet type with the free surface is presented to design a three-dimensional body corresponding to a prescribed pressure distribution. The free surface boundary condition is linearized with respect to the oncoming flow, and computed by four-point finite difference scheme. Sample designs for submerged spheroids and Wigley hull are carried out to demonstrate the stable convergence, the effectiveness and the robustness of the method. The design of a 5500TEU container carrier is performed with respect to reduction of the wave resistance. To reduce the wave resistance, calculated pressure on the hull surface is modified to have the lower fluctuation, and is applied as a Dirichlet type dynamic boundary condition on the hull surface. The designed hull form is verified to have the lower wave resistance than the initial one not only by computation but also by experiment.     

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.     

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.     

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 laboratory study of directional spectra with maximum likelihood method-Ⅰ Developing wind wave

measurements by a circular array consisting of 18 wave gauges in a large wave tank, directional spectra of wind-generated waves in deep water are systematically determined by using maximum likelihood method The investigations reveal that the angular spreading of the wave energy is consistent with cos^2s(θ/2) proposed by Longuet-Higgins et al.(1963, Ocean Wave Spectra, 111-136), if the bimodal distributions of wave energy are not taken into account. Bimodality occurring on higher frequency than peak frequency is too rare to affect our whole resalts. Surprisingly, a much broader directional spreading than that of the field, which is interpreted by the strongly nonlinear energy transfer because of the very young waves in laboratory, is found. The parameter s depends on frequency in the same way as observed by Mitsuyasu et al.(1975, Journal of Physical Oceanography, 5, 750-760) and Hasselmann et al.(1980, Journal of Physical Oceanography, 10, 1264-1280) in the field, and the relationship between the four nondimensional parameters s_m, f_o, b₁ and b₂, determining the directional width, and ū₁₀ (corresponding to the inverse of wave age) are given respectively. The observed distributions are found to agree well with the suggestion of Donelan et al.(1985, Philosophical Transaction of Royal Society of London, A315, 509-562) when applied to field waves.     

Comparison of numerical solution and semi-empirical formulas to predict the effects of important design parameters on porpoising region of a planing vessel

Stability of the marine vessels in different conditions is one of the most important problems in the design of a planing vessel. In this research, the effects of some important design parameters (mass, longitudinal center of mass, deadrise angle, and length) of DTMB 62 model 4667-1 planing hull on the drag and also on the longitudinal dynamic stability (porpoising) are investigated in the velocity range of 2.12–8.486 m/s in calm water. In this paper, both numerical simulation of Reynolds Average Naiver Stokes (RANS) equations and semi-empirical formulas of Savitsky are used to analyze the motion of a 4667 planing vessel in calm water with two degrees of freedom (2DOF). For this purpose a finite volume, ANSYS-FLUENT, code is used to solve the Navier-Stokes equations for the simulation of the flow field around the vessel. In addition, an explicit VOF scheme and SST- Kω model is used with dynamic mesh scheme to capture the interface of a two-phase flow and to model the turbulence respectively, in 2DOF model (heave and pitch). Also, the results of both methods are compared with each other. According to the present results, changing the aspect ratio of the vessel and also the longitudinal center of gravity have the most effect on the porpoising region.     

基于杨赤中插值法的坡度、坡向计算的实现与分析

应用杨赤中滤波推估法,将地形碎部点高程观测值看作一种包含规律性和随机性变化的复合变量,求其二项系数加权游动平均,建立估值数学模型,求得各格网点高程估值。将离散的数据经插值计算转换为DEM数据,并根据其网格顶点来分析、推算和显示网格单元的坡度和坡向。     

基于声波探测海底浅层沉积物特性的方法研究

利用声波在水中的传输特性进行海底浅层沉积物特性的探测。在实践应用中有实验室方法、现场测量方法和海底声学遥测方法,并研制出大量的仪器。由于海上和海底实际条件的限制,声波探测方法都存在一定的不足之处,由此提出了声衰减测试沉积物特性以及建立实验室模拟海底环境探测方法的设想。