胶州湾三维变动边界的潮流数值模拟

基于Blumberg等（1986）的河口、陆架和海洋模式，引入变边界处理技术，建立了胶州湾三维变动边界的潮流模型，模型以σ-正交曲线坐标下三维非线性潮波方程为基本方程，引入2.5阶瑞封闭方程组，采用分裂算子法数值求解方程组，利用湍封闭方程求解垂直紊动粘滞系数，采用干湿网格方法模拟潮流漫滩过程，三维变动边界潮流模型计算结果与实测值吻合良好。     

乐清湾的动力沉积过程

主要采用二维潮波、泥沙输运途径的数值模拟，定量地展示了泥沙在乐清湾内的输运途径和沉降过程。模拟结果表明，乐清湾大致以10m等深线为界，其东部的涨潮流主要来自玉环岛与横趾山之间的深槽；水流在一周日内的轨迹为一宽约1km的狭带；此流态有利于泥沙向南的净输运，从面明利于此海域深水航道的维持。10m等深线以西水流轨迹为一较宽的带状结构，有利于泥沙向西、西北的净输运。华歧潮滩附近水流运动轨迹得于泥沙向北、向西的净输运，从而有利于华歧潮滩的发育。本文还定性讨论了楚门潮滩特征的成因,主要是漩门口堵截前丰富的湾外泥沙以及驻波型潮波与漏斗状地形。     

闽江口水动力和污染扩散数值模拟

本文建立了闽江口二维水动力、污染扩散数值模型，模型采用ＡＤＩ双向隐式求解。在水动力模拟中，对岸界的丁坝有用线性边界进行拟合，并考虑到河口区侧向摩擦作用；在污染扩散模拟中，用动态拟合的方法对浓度水边界进行了较为合理的改进。计算结果表明，本模型是成功的。     

Numerical simulation of random wave slamming on structures in the splash zone

The numerical investigation of random wave slamming on superstructures of marine structures in the splash zone is presented in this paper. The impact pressures on the underside of the structure are computed based on the improved volume of fluid method (VOF). The governing equations are Reynolds time-averaged equations and the two equation k– model. The third order upwind difference scheme is applied to the convection term to reduce the effect of numerical viscosity. The numerical wave flume with random wave-maker suitable for VOF is established. Appropriate moving contact-line boundary conditions are introduced to the model wave in contact with and separated from the underside of structure. Parametric studies have been carried out for different incident waves, structure dimensions and structure clearance. The numerical results are verified by the experimental results.     

质点-网格海冰模式在渤海的数值预报试验

在海冰动力学和热力学的研究基础上,开发应用于渤海的质点-网格海冰模式.该模式采用质点-网格法,有效地避免了传统海冰模式的数值扩散问题.该模式采用冰厚分布函数取代传统渤海业务海冰预报模式中平整冰、堆积冰和开阔水的三类海冰分型.应用质点-网格海冰模式对2003年～2004年冬季渤海冰情进行了业务化逐日数值预报试验,并对预报结果进行统计检验和分析比较.     

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.     

Dynamics and control of a towed underwater vehicle system, part I: model development

Autonomous vehicles are being developed to replace the conventional, manned surface vehicles that tow mine hunting towed platforms. While a wide body of work exists that describes numerical models of towed systems, they usually include relatively simple models of the towed bodies and neglect the dynamics of the towing vehicle. For systems in which the mass of the towing vehicle is comparable to that of the towed vehicle, it becomes important to consider the dynamics of both vehicles. In this work, we describe the development of a numerical model that accurately captures the dynamics of these new mine hunting systems. We use a lumped mass approximation for the towcable and couple this model to non-linear numerical models of an autonomous surface vehicle and an actively controlled towfish. Within the dynamics models of the two vehicles, we include non-linear controllers to allow accurate maneuvering of the towed system.     

用样条函数计算结果研究冲绳海槽热源机制

对冲绳海槽１９８４－１９９０年的实测热流值，利用数字滤波方法计算，得到反映该区深部地壳热状态的区域热流值为８０－１６０ｍＷ／ｍ＾２，用Ｂ样条函数法对其作数值模拟，计算出海槽地温场模式，并以此研究冲绳海槽海底高热流的形成机制，得出冲绳海槽热流值的区域性变化与綦中地质构造格局呈正相关，高热流是该区现代构造活动所致，其明显特征是地壳减薄，地幔上拱等结论。     

Influence of large-diameter pipe pile driving on surrounding marine soils and adjacent submarine pipelines

Abstract

With the large-scale development and utilization of ocean resources and space, it is inevitable to encounter existing submarine facilities in pile driving areas, which necessitates a safety assessment. In this article, by referring to a wharf renovation project as a reference, the surrounding soil response and buried pipe deformation during pile driving in a near-shore submarine environment are investigated by three-dimensional (3D) numerical models that consider the pore water effect. Numerical studies are carried out in two different series: one is a case of a single pile focusing on the effect of the minimum plane distance of the pile–pipe, and the other is a case of double piles focusing on the effect of the pile spacing.     

Wave reflection from partially perforated-wall caisson breakwater

In 1995, Suh and Park developed a numerical model that computes the reflection of regular waves from a fully perforated-wall caisson breakwater. This paper describes how to apply this model to a partially perforated-wall caisson and irregular waves. To examine the performance of the model, existing experimental data are used for regular waves, while a laboratory experiment is conducted in this study for irregular waves. The numerical model based on a linear wave theory tends to over-predict the reflection coefficient of regular waves as the wave nonlinearity increases, but such an over-prediction is not observed in the case of irregular waves. For both regular and irregular waves, the numerical model slightly over- and under-predicts the reflection coefficients at larger and smaller values, respectively, because the model neglects the evanescent waves near the breakwater.