塘沽海区海底地形的SAR影像仿真与反演研究

利用袁业立、金梅兵(1997)提出的海底地形SAR影像仿真与反演模型，对渤海塘沽海区的一张Radarsat SAR影像进行了仿真和水深反演研究。研究结果表明：仿真影像与真实SAR影像基本吻合，反演水深与实际水深也有较好的一致性；进一步证实了袁业立(1997)SAR成像机理的正确性和在中国近海利用SAR影像进行浅海水深探测的可行性。     

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.     

The dimensions of sand ripples in full-scale oscillatory flows

New large-scale experiments have been carried out in two oscillatory flow tunnels to study ripple regime sand suspension and net sand transport processes in full-scale oscillatory flows. The paper focuses on ripple dimensions and the new data are combined with existing data to make a large dataset of ripple heights and lengths for flows with field-scale amplitudes and periods. A feature of the new experiments is a focus on the effect of flow irregularity. The combined dataset is analysed to examine the range of hydraulic conditions under which oscillatory flow ripples occur, to examine the effects of flow irregularity and ripple three-dimensionality on ripple dimensions and to test and improve existing methods for predicting ripple dimensions.The following are the main conclusions. (1) The highest velocities in a flow time-series play an important role in determining the type of bedform occurring in oscillatory flow. Bedform regime is well characterised by mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth peak velocities in the case of irregular flow. (2) For field-scale flows, sand size is the primary factor determining whether equilibrium ripples will be 2D or 3D. 2D ripples occur when the sand D₅₀ ≥ 0.30 mm and 3D ripples occur when D₅₀ ≤ 0.22 mm (except when the flow orbital diameter is low). (3) Ripple type (2D or 3D) is the same for regular and irregular flows and ripple dimensions produced by equivalent regular and irregular flows follow a similar functional dependence on mobility number, with mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth velocities in the case of irregular flow. For much of the ripple regime, ripple dimensions have weak dependency on mobility number and ripple dimensions are similar for regular and irregular flows with the same flow orbital amplitude. However, differences in ripples produced by equivalent regular and irregular flows become significant at the high mobility end of the ripple regime. (4) Ripple dimensions predicted using the Wiberg and Harris formulae are in poor agreement with measured ripple dimensions from the large-scale experiments. Predictions based on the Mogridge et al. and the Nielsen formulae show better overall agreement with the data but also show systematic differences in cases of 3D ripples and ripples generated by irregular flows. (5) Based on the combined large-scale data, modifications to the Nielsen ripple dimension equations are proposed for the heights and lengths of 2D ripples. The same equations apply to regular and irregular flows, but with mobility number appropriately defined. 3D ripples are generally smaller than 2D ripples and estimates of 3D ripple height and length may be obtained by applying multipliers of 0.55 and 0.73 respectively to the 2D formulae. The proposed modified Nielsen formulae provide an improved fit to the large-scale data, accounting for flow irregularity and ripple three-dimensionality.     

Examination of empirical formulas for wave shoaling and breaking on steep slopes

Estimation of the wave height transformation of shoaling and breaking is essential for the nearshore hydrodynamics and the design of coastal structures. Many empirical formulas have been well recognized to the wave height transformation, but most of them were only applicable for gentle slopes. This paper reports the experimental results of wave shoaling and breaking over the steep slopes to examine the applicability of the previous empirical formulas. Two steep bottom slopes of 1/3 and 1/5, and one gentle slope of 1/10 were conducted in the present experiments. Experimental results show that the shoaling distance of steep slopes become short and the surface waves may be partially reflected from the steep bottom, thus the estimation of wave shoaling using the well-known previous formula did not conform completely to the experimental results. The previous empirical formulas for the wave breaking criteria were also examined, and the modified equations to the steep beaches were proposed in this work. A numerical model was finally adopted to calculate the wave height transformation in the surf zone by introducing the modified breaking index.     

The nonlinear hydrodynamic model for simulating a ship steering in waves with autopilot system

In the paper, a hydrodynamic numerical model including wave effect is developed to simulate ship autopilot systems by using the time domain analysis. The PD controller and the sliding mode controller are adopted as the autopilot systems. The differences of simulation results between two controllers are analyzed by cost function composed of heading angle error and rudder deflection, either in calm water or in waves. The results in calm water show that both controllers are tracking well for the desired route with the similar cost function value by tuning the key design parameters. However, the course tracking ability of the controller using sliding mode in waves is generally better even the cost function value is similar.     

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.     

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

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

近岸海浪模式在中国东海台风浪模拟中的应用--数值模拟及物理过程研究

较为详细地介绍了基于能量平衡方程的第三代近岸海浪数值模式SWAN（Simulation Waves Nearshore）及其包含的物理过程（风生浪、底摩擦、白浪耗散、深度诱导波破碎、非线性波－波相作用等），并利用该模式对影响杭州湾－长江口沿岸海域的一次台风浪过程进行了模拟研究：模式所需风场由藤田台风风场模型嵌入对应台风特征等压线，并对相应时段的NCAR／NCEPT资料、单站资料进行同化后提供；利用自嵌套的方式提供波谱边界条件；模式模拟的结果与实际海浪观测资料相符较好，在此基础上，研究了底摩擦、深度诱导波破碎、三波相互作用等物理过程联合对近岸台风浪的影响，初步认识了它们在近岸台风浪生成、传播过程中的重要作用。     

海浪周期关系的研究

本文利用现场观测数据、实验室风浪槽观测数据和由文氏谱及JONSWAP谱生成的模拟数据,研究了有效周期与谱平均周期、谱峰周期与谱平均周期以及有效周期与谱峰周期的关系,通过数据拟合给出了相应的关系式.研究发现,有效周期与由谱的负阶矩计算的平均周期之间的关系更加稳定,并且有效周期与负2阶矩计算的平均周期几乎相等,均代表了海浪主导波对海浪平均周期的贡献.