非线性波传播的新型数值模拟模型及其实验验证

以一种新型的Boussinesq型方程为控制方程组,采用五阶Runge-Kutta-England格式离散时间积分,采用七点差分格式离散空间导数,并通过采用恰当的出流边界条件,从而建立了非线性波传播的新型数值模拟模型.通过对均匀水深水域内波浪传播的数值模拟说明,模型能较好地模拟大水深水域和强非线性波的传播.通过设置不同的入射波参数来进行潜堤地形上波浪传播的物理模型实验,并将数值解与物理模型实验结果进行了比较.     

复杂地质体粘滞弹性波正演模拟的有限差分法

运用时空域的二阶中心差分法 ,由运动微分方程、自由边界条件及人工边界上的吸收边界条件 ,导出了二维非均匀粘滞弹性介质中波传播的二阶显示差分格式。该法可以方便且精确地处理自由表面、人工边界及交接层面。实现了在微机上用有限差分法对地震波在传播过程中的波型转换及波形变化的全波场模拟     

非线性波传播的新型数值模拟模型及其实验验证引入变换速度变量

以一种新型的含变换速度变量的Boussinesq型方程为控制方程组,采用五阶Runge-Kutta-England格式离散时间积分,采用七点差分格式离散空间导数,并采用恰当的出流边界条件,从而建立了非线性波传播的新型数值模拟模型.对均匀水深水域内波浪传播的数值模拟,说明在引入变换速度后进一步增大了模型的水深适用范围.对潜堤地形上波浪传播的数值模拟说明,在引入变换速度后进一步提高了模型的数值模拟精度.     

完全非线性孤立波的直墙反射

报道了应用边界积分方法模拟完全非线性孤立波的传播与直墙反射,给出了波形演变过程。结果表明,本模型对计算孤立波的传播与直墙反射是有效的。三阶Boussinesq方程的孤立波解比低阶方程的孤立波解更接近完全非线性的数值解.当来波波高增大时,孤立波直墙反射的相位滞后变小。若考虑大波高孤立波的直墙反射或波——波相互作用,一阶理论预报的相位滞后往往低估实际情况。     

港域波浪数学模型的改进与验证

通过物理模型对改进的港内波浪传播变形数学模型进行验证。该数学模型以推广的时变缓坡方程为控制方程,采用含松弛因子的ADI法求解,并对波浪反射和透射边界模拟方法进行改进。先通过物理模型试验确定斜向浪入射条件下抛石防波堤前的波浪反射系数,作为数学模型中部分反射边界模拟的依据。然后进行了一个典型港口内波浪折射、绕射和反射的模型试验,测量港内波浪分布。对比模型试验和数学模型计算的结果表明,数学模型可较好地模拟港内复杂地形和边界条件下规则波和不规则波的传播变形。     

一维全非线性Green-Naghdi水波方程数值模型

建立了求解一维全非线性Green-Naghdi水波方程的中心有限体积/有限差分混合数值格式。采用结构化网格对守恒形式的控制方程进行离散和积分,界面数值通量采用有限体积法计算,剩余项则采用中心有限差分格式求解。其中,采用中心迎风有限体积格式计算控制体界面数值通量,并结合界面变量的线性重构方法,使其在空间上具有四阶精度,通过引入静压重构技术和波浪破碎指标使模型具备处理海岸水-陆动边界及波浪破碎的能力。时间积分则采用具有总时间变差减小(Total Variation Diminishing,TVD)性质的三阶龙格-库塔法进行。应用该模型对孤立波在常水深和斜坡海岸上的传播过程及规则波跨越潜堤传播的实验进行了数值模型研究,数值计算同解析解及实验数据吻合良好。     

海堤坡角变化对海啸波传播过程影响数值分析

针对近岸海域极端海况的防灾减灾问题,准确模拟追踪海啸波传播过程,再现液面局部射流、崩破波等波面湍动现象。建立SPH数值水槽,边界条件基于固壁粒子法,减少海堤坡角改变对计算域精度的影响,讨论粒子间距设置对模拟精度的影响。模拟7种不同海堤坡角下波浪的爬高与衰减情况,讨论了海堤坡角变化对消波系数的影响。当粒子间距设置为0.002 m时,模型能准确地捕捉海啸波的强非线性现象。随着海堤坡角增大,海啸波峰值爬升率增大,波浪越堤后因崩破波的产生,波能衰减进一步加剧。当坡角较小时,消波系数随坡角的增大提升明显,而后趋于缓慢增长。     

管式透空防波堤反射系数的计算

本文从规则波的研究入手,基于反射波与入射波势函数线性叠加的概念,利用积分形式的流量连续性条件来处理透空堤面上的间断边界条件,确定了规则波与管式透空堤相互作用时的反射系数。在此基础上,利用线性海浪谱的概念,分析了海浪与管式透空堤的相互作用问题。通过对模型试验结果的分析表明,本文提出的理论与实验结果符合良好。     

完全非线性孤立波的稳态解

报道了应用边界积分方法模拟完全非线性孤立波的传播,给急剧记解的波形,流束和数值计算结果。结果表明,本模型对计算孤立波的传播是有铲的。当ａ／ｈ〉０．３时,自由液面上的水平流速、底部流速和垂向平均流速之间的差别是明显的。三阶Ｂｏｕｓｓｉｎｅｓｑ方程的孤立波解比低阶方程的孤立波解更接近本文完全非线性的数值解。     

斜向波对台阶地形上有限吃水半无限长弹性板的散射作用

基于线性势流理论和欧拉-伯努利梁理论,采用匹配特征函数展开法建立斜向波对台阶地形上有限吃水半无限长弹性板散射的解析解。与以往研究不同,本文同时考虑波浪入射角度、弹性板吃水、三种不同板端条件（自由板端、简支板端、固定板端）以及海底地形的影响。本文还直接将板端条件引入匹配边界条件,建立线性方程组求解速度势中的展开系数,计算过程简单、高效。计算结果表明：本文方法收敛性好,并满足能量守恒方程。将本文计算结果与文献中的特例结果进行对比,吻合良好。分析了波浪入射角度、弹性板吃水、不同板端条件以及海底地形对弹性板水动力参数的影响,研究结果可以为工程设计提供参考。     

On the Interaction of Surface Waves with An Elastic Plate of Finite Length in Head Seas

An eigen-function expansion method based on a new orthogonal inner product is proposed by Sahoo et al. (2000) for the study of the hydroelastic response of mat-type VLFS in head seas. However, their main emphasis is on the effect of edge conditions and they assume that the plate is of a semi-infinite length. In reality, the plate is of finile length. For consider-alion of the finite length effect, the reflection and transmission from the other end must be considered. The effect of this reflection and transmission on the hydroelaslic response of VLFS is of interest for praclical application. Furthermore, the physi-cal meaning of the new inner producl was not given in their paper. In this paper, it is shown that the new inner product can be derived from the governing equation and the bottom boundary conditions. Then the same eigen-function expansion method is adopted for the study of the hydroelastic response of an elastic plate of finite length in surface waves. Delailed comparisons are made between the     

Analytic Solution of Interaction Between Surface Waves and Floating Elastic Plates with Elastic End-Restraints

Based on the dynamic theories of water waves and Mindlin plates,the analytic solution of interaction between surface waves and two-dimensional floating elastic plates with edge-restraint is constructed by use of the Wiener-Hopf technique.Firstly,without regard for elastic edge restraint,the wave-induced responses of elastic floating plate analyzed by the present method are in good agreement with the results from literature and experimental results.Therefore,it can be shown that the present method is valid.Secondly,three end-restraint cases(i.e.,the left-end elastic restraints,the both-end elastic restraints,and the right-end elastic restraints)are proposed to reduce the vibration of floating plates,in which the spring is used to connect the sea bottom and the floating plate's left(or right)edge.The relations between the spring stiffness and the parameters of wave-induced responses of floating plates are discussed.Moreover,the effective method to reduce the vibration of floating elastic plates can be obtained through comparison.     

目标回波中的弹性信息提取方法研究

在主动声纳探雷识别中,弹性波提取方法一直是研究的难点。研究表明回波中,几何亮点主要表现为波形包络的强起伏的慢变化,在频域上属于低频,而弹性信息则隐含在包络中,其变化较快且能量微弱,相对于包络来说是比较弱的高频成份。运用中心差分的方法,可使高频成份得到较大的增强,从而实现弹性信息提取的目的。仿真和实验数据均取得了良好的效果。     

海底沉积物弹性参数计算方法及误差分析

由海底沉积物的纵波速、横波速及湿密度,根据岩土力学原理,整理出一组计算海底沉积物弹性参数的公式;以该公式计算了南海南部海域浅层沉积物的弹性参数,计算结果与已有文献资料数据做了比较,分析了两者误差的原因.     

正压冲固平台结构分析中基础边界条件的影响

正压冲固平台是一种采用短桩加固基础的新型海洋采油平台,对于这种新型的平台结构,在结构分析和构件强度校核中必须考虑其有限元模型的基础边界条件处理问题。本文提出了正压冲固平台有限元计算模型中基础边界条件的一种简化方法,将两个水平方向的扭转自由度简化为扭转弹簧边界元,其余自由度简化为固定约束。通过计算分析得到了不同的边界约束刚度系数的取值对平台总体位移和强度校核应力的影响及变化趋势。结论是,平台结构对约束刚度系数K的反应在10^4～10^4之间时比较明显,对K的敏感度最为强烈：在此范围之外,平台反应分别接近于简支约束情况和刚性约束情况。尤其对于接近约束边界的单元,其应力变化最敏感。     

Elastic scattering of electrons from water molecule

This paper uses the momentum--space optical potential method to calculate the e--H₂O scattering elastic cross sections at the energy range from 6 eV to 50 eV, and the differential cross sections in the angle from 0° to 180° at 40 eV and 50 eV. The polarisation is taken into account via an emphab initio equivalent-local potential. The cross sections are compared with experimental measurements and other theoretical calculations.     

An active wave generating–absorbing boundary condition for VOF type numerical model

The objective of the present work is to discuss the implementation of an active wave generating–absorbing boundary condition for a numerical model based on the Volume Of Fluid (VOF) method for tracking free surfaces. First an overview of the development of VOF type models with special emphasis in the field of coastal engineering is given. A new type of numerical boundary condition for combined wave generation and absorption in the numerical model VOFbreak² is presented. The numerical boundary condition is based on an active wave absorption system that was first developed in the context of physical wave flume experiments, using a wave paddle. The method applies to regular and irregular waves. Velocities are measured at one location inside the computational domain. The reflected wave train is separated from the incident wave field in front of a structure by means of digital filtering and subsequent superposition of the measured velocity signals. The incident wave signal is corrected, so that the reflected wave is effectively absorbed at the boundary. The digital filters are derived theoretically and their practical design is discussed. The practical use of this numerical boundary condition is compared to the use of the absorption system in a physical wave flume. The effectiveness of the active wave generating–absorbing boundary condition finally is proved using analytical tests and numerical simulations with VOFbreak².     

Transparent boundary condition for simulating nonlinear water waves by a particle method

Open boundaries are important when simulating water waves. In this study, a transparent boundary condition at an open boundary was developed for simulating nonlinear water waves propagating to a distant area using the Moving Particle Semi-implicit method. The novelty of this study is that the technique of wave analysis used in the experiment was introduced into the particle simulation to absorb incident waves; the simulation cost was reduced by employing inflow and outflow regions instead of a long dissipation region. Incident waves in front of the boundary were evaluated using Fourier analysis, and the particles on the transparent boundary were forced to move at the velocity of the analytical solution for Stokes waves in order to absorb the incident waves. The analysis was restricted to periodic waves. Wave propagation was simulated for two wave periods using the developed transparent boundary condition. The results showed that this transparent boundary transmitted the incident waves with small reflection and the simulation cost was lower than that for wave damping by a conventional highly viscous region.     

二维水波辐射问题的二阶时域解

用二阶理论在时域范围内研究二维二阶非线性水波的辐射问题,一阶、二阶问题分别满足各自的自由表面条件和物面条件,采用人工阻尼来吸收反射波.计算采用有限元法,对水平半圆柱体在水面的振荡问题进行了计算,计算了波高时间历程和水动力,计算结果与有关文献中的数据相比符合较好.     

Hydrodynamics of a submerged hydrofoil advancing in waves

The hydrodynamic problem of a hydrofoil travelling at constant speed in water waves has been investigated through velocity potential theory. The boundary conditions on the free surface have been linearized, and the effects are accounted for through the Green function. The overall problem is decomposed into the steady forward speed problem and periodic wave radiation and diffraction problems. Each of these problems is solved using the boundary integral equation over the hydrofoil surface together with a vortex sheet behind the trailing edge. The body surface boundary condition is imposed on its mean position. As a result the steady potential will contribute a well-known m_j term to the body surface boundary condition on the radiation problem. The numerical difficulty in dealing with this term is effectively resolved through a difference method. The effects of the thickness on the wave radiation and diffraction are investigated. The applicability of various reciprocity relationships in this problem is discussed.