潜堤对波浪在直立式防浪墙上爬高的影响

对直立式防浪墙前潜堤的空间布置对波浪在直墙上爬高的影响进行了研究.建立简化模型,即直立式防浪墙前海床设计为平底,潜堤设计为直立薄板. 改变潜堤高度及其与防浪墙之间距,研究不同的潜堤布置对波浪在直墙上爬高的影响. 针对线性波浪场,利用数值波浪水槽模拟了潜堤作用下直墙上波浪的爬高现象. 同时建立了理论模型,系统地分析了潜堤布置形式对直墙上波浪爬高的影响. 研究结果表明潜堤相对于水深的高度越大,对波浪在直墙上爬高的影响越大,而在一定的相对高度条件下,潜堤与直立墙之间距对波浪爬高的影响呈现出一定的规律性.      

非周期波浪与直墙作用的非线性数值研究

基于时域高阶边界元方法,建立了完全非线性二维数值波浪水槽,对非周期波浪与直墙的相互作用问题进行了模拟和研究.自由表面满足完全非线性自由水面运动学和动力学边界条件,采用混合欧拉-拉格朗日方法追踪瞬时自由面流体质点,采用四阶Runge-Kutta法对下一时间步的波面和自由面速度势进行更新.采用加速度式法求解直墙表面速度势的时间导数,对瞬时物体湿表面上的水动力压强积分,得到作用在物体上的瞬时波浪力.首先,将全非线性与Serre-Green-Naghdi(SGN)模型的结果进行了对比分析,发现对于大幅值双入射波问题,仅满足弱色散关系的SGN模型大大低估了最大波浪爬高;其次,研究了双入射波与直墙的非线性作用问题,发现线性预报对波浪最大爬高有较大低估,而波浪的非线性成分不只导致了自由面爬高的异常增大,也引起了局部自由面的高频振荡,该物理过程中,直墙所受的波浪载荷,也展示出了与波浪爬高相似的非线性特性;最后,对波浪爬升和波浪力的时间历程进行了频谱分析,发现入射主频波的部分能量传递给了更高频的波浪成分,反映出该问题具有典型的非线性特性.     

A deforming finite element mesh for use in moving one-dimensional boundary wave problems

The finite element method is developed to solve the problem of wave run-up on a mild, plane slope. A novel approach to implementing a deforming mesh of one-dimensional, three-node, isoparametric elements is described and demonstrated. The discrete time interval (DTI), arbitrary Lagrangian–Eulerian (ALE) and space–time element (STE) methods are used to solve the unsteady one-dimensional shallow water wave equations. The boundary condition required is simply the seaward water surface elevation, and although the method has only been tested for monochromatic waves, it should be equally valid for any sea state which can be described as a water surface elevation as a function of time. All three solution methods are shown to given good results. Time histories of the terms of the governing equations are calculated and used to demonstrate how the ALE and STE methods account for mesh deformation. The model could be extended to two dimensions, which would have practical application to the run-up of obliquely incident waves. © 1997 John Wiley & Sons, Ltd.     

Numerical analysis of a moving boundary problem in coastal hydrodynamics

A finite element model to tackle the moving boundary problem of wave run-up on moderately steep slopes is developed. The special aspects considered in this study are (1) the modification of shallow water equations to accommodate the effect of vertical accelerations and (2) the use of Lagrangian acceleration coupled with an element that adapts itself to the moving boundary closely. The pressure term in the one-dimensional momentum equation is derived using the Eulerian equation in the vertical direction. This takes care of the vertical accelerations which are significant during the motion of a wave on moderately steep slopes. The element near the boundary is allowed to change its dimension so that the fluid boundary is closely followed. Such a flexible element precludes the need for approximation of the variables with regard to the indefinite position of the boundary. This element is split into two when its dimension becomes unduly large compared to the unchanging elements. The need for such a splitting is shown by an examination of the entries in the global matrix. Results of water profile as a wave runs up a structure are given. A brief history of the work on similar problems is outlined.     

Estimation of Wave Propagation using a Kalman Filter

An estimation method for gravity water wave propagation combining the finite element method and the Kalman filtering technique is presented. The method is further simplified to reduce the computational effort. According to the results obtained with this method, it is possible to estimate the water surface profile by taking into account the measured data, and to filter out the noise included in the measured data and the finite element approximation.     

An efficient RPIM for simulating wave motions in saturated porous media

An efficient meshless technique for simulating wave motions in saturated porous media is introduced in this paper. Using radial point interpolation method (RPIM) and considering the precision requirement in wave motion simulation, a lumped-mass RPIM for saturated porous media with clear physical concepts is derived. Combining explicit time integration with the lumped-mass RPIM leads to a decoupled radial point interpolation method (decoupled RPIM) for dynamic analysis without solving algebra equation set. A two-dimensional problem is calculated by the proposed decoupled RPIM, the ordinary RPIM, and the finite element method (FEM). The accuracy and efficiency of these methods are compared, which demonstrate that the proposed decoupled RPIM has comparable accuracy with the ordinary RPIM and the finite element method, but has high efficiency in simulating wave motions in saturated porous media.     

Uncertainty analysis of deflagration-to-detonation run-up distance

Three methods were adopted to estimate the deflagration-to-detonation run-up distance in a smooth tube, which comprise (1) the measurement of the propagation speed of pressure or combustion waves compared with the C–J detonation speed, (2) the time of the onset of detonation by the emission of visible light and the trajectory of pressure wave or combustion waves, and (3) the trajectory intersection with the presence of retonation wave. A nonstationary cross-correlation technique was applied to evaluate the uncertainty in estimating the run-up distance. Evaluation of the pressure wave (pressure wave speed or the pressure wave trajectory) appears to be more suitable to determine the deflagration-to-detonation run-up distance. PACS 47.40.-x Communicated by Z. Jiang     

有限元畸变单元离散模型中的SH波动

本文采用晶格动力学分析方法,研究了有限元畸变单元离散模型中的SH波的传播规律。结果表明单元畸变对波动规律影响显著,当单元畸变较小时,单元畸变对波动传播并不产生明显的不利影响,当单元畸变较大时,这一影响不可忽略。通过详细分析和对比集中质量、混合质量和一致质量模型中的波动规律,评价了不同有限元模型的抗畸变能力。文中还讨论了时域离散化与单元畸变联合作用对波动传播的影响。     

Spectrally formulated wavelet finite element for wave propagation and impact force identification in connected 1-D waveguides

In this paper, a novel wavelet based spectral finite element is developed for studying elastic wave propagation in 1-D connected waveguides. First the partial differential wave equation is converted to simultaneous ordinary differential equations (ODEs) using Daubechies wavelet approximation in time. These ODEs are then solved using finite element (FE) technique by deriving the exact interpolating function in the transformed domain. Spectral element captures the exact mass distribution and thus the system size required is very much smaller then conventional FE. The localized nature of the compactly supported Daubechies wavelet allows easy imposition of initial-boundary values. This circumvents several disadvantages of the conventional spectral element formulation using Fast Fourier Transforms (FFT) particularly in the study of transient dynamics. The proposed method is used to study longitudinal and flexural wave propagation in rods, beams and frame structures. Numerical experiments are performed to show the advantages over FFT-based spectral element methods. The efficiency of the spectral formulation for impact force identification is also demonstrated.     

等波高双孤立波直墙爬高的数值模拟

基于RANS方程、VOF方法以及修正的Goring造波方法建立了模拟活塞式推波板运动的二维数值波浪水槽,实现了双孤立波直墙爬高的数值模拟．利用动边界技术模拟造波机推波板的运动,有效地实现了不同波峰间距双孤立波的造波方法．在验证单孤立波直墙爬高的基础上,模拟了不同相对波高、相对波峰间距的等波高双孤立波的直墙爬高过程,给出了波面、速度场及波动能量的变化规律．数值模拟结果表明：对于等波高的双孤立波,当入射波波高较大及两个波峰间距相对较小时,跟随在后孤立波的爬高放大系数小于先导孤立波的爬高放大系数;双孤立波在直墙爬高过程中,波动场的势能时间过程线呈现三峰形态,其中居中的最大势能峰值出现在第二个孤立波与经直墙反射后反向传播的第一个孤立波完全对撞的时刻．     

成层半空间出平面自由波场的一维化时域算法

提出了一种计算出平面SH波斜入射时弹性水平成层半空间中自由波场时域计算的一维化有 限元方法. 在进行有限元网格划分时,竖向单元取满足有限元模拟精度的任意尺寸,水平向 网格尺寸由时间离散步长和水平视波速确定,并自动进行虚拟网格划分. 基底设置人工边界, 并将波动输入转化为等效荷载施加在边界节点上. 然后将集中质量有限元法和中心差分法相 结合建立节点运动方程,并将水平方向相邻节点的运动用该节点相邻时刻的运动表示,从而 将求解节点运动的二维方程组转化为一维方程组. 求解此方程组,即得到自由场中竖向一列 节点的运动. 最后根据行波传播的特点,可方便地确定全部自由波场. 理论分析和数值算例 表明,该方法具有较高的精度和良好的稳定性.     

Wave propagation in double helical rods

The propagation of waves in helical rods has been studied extensively. However, studying the wave propagation in double helical rods have received less attention although this can be useful in multiple fields of science and engineering. Obtaining an analytical model for a double helical rod is challenging since the curvature and tortuosity are not constant. Thus, resolving the wave behaviour analytically is nearly impossible. In this paper, wave propagation in a double helical rod will be studied using the wave and finite element method which is a technique that can be used to model homogeneous and periodic one and two dimensional structures based on the periodic structure theory. For modelling a double helical rod, the finite element model of a single turn is processed using Bloch waves. The dispersion curves and wavemodes are obtained and the similarities and differences of waves in helical and double helical rods are highlighted.     

A particle finite element method applied to long wave run‐up

This paper presents a Lagrangian–Eulerian finite element formulation for solving fluid dynamics problems with moving boundaries and employs the method to long wave run‐up. The method is based on a set of Lagrangian particles which serve as moving nodes for the finite element mesh. Nodes at the moving shoreline are identified by the alpha shape concept which utilizes the distance from neighbouring nodes in different directions. An efficient triangulation technique is then used for the mesh generation at each time step. In order to validate the numerical method the code has been compared with analytical solutions and a preexisting finite difference model. The main focus of our investigation is to assess the numerical method through simulations of three‐dimensional dam break and long wave run‐up on curved beaches. Particularly the method is put to test for cases where different shoreline segments connect and produce a computational domain surrounding dry regions. Copyright © 2006 John Wiley & Sons, Ltd.     

波浪与外壁开孔双筒柱群的相互作用

应用速度势的特征函数展开和透空壁内两壁间压力差和流体速度成正比的线性模型,建立 了波浪与外壁开孔同轴双筒桩柱群相互作用的线性解析解. 应用这一模型进行了数值计算, 用以检验孔隙系数对双筒柱上的波浪力和波面高度的影响. 结果表明,外壁孔隙系数的增加 对减小波浪力和柱外波面高度有很大影响.     

水下声波传播过程的时域间断Galerkin有限元方法

本文构建了声压波动方程的改进时域间断Galerkin有限元方法.传统时域连续有限元方法在计算高梯度、强间断特征水中声波传播问题时往往会出现虚假数值振荡现象,这些数值振荡会影响正常波动的计算精度.为了解决这一问题,本文通过引入人工阻尼的方式构建了改进的时域间断Galerkin有限元方法,并针对具有高梯度、强间断特征的多障...     

A non-hydrostatic numerical scheme for dispersive waves generated by bottom motion

A numerical scheme based on the staggered finite volume method is presented at the aim of studying surface waves generated by a bottom motion. We address the 2D Euler equations in which the vertical domain is resolved only by one layer. The resulting non-hydrostatic scheme is used to simulate surface waves generated by bottom motion in a water tank. Here we mimic Hammack experiments numerically, in which a bed section is moved upwards or downwards, resulting in transient dispersive waves. For an impulsive downward bottom thrust, free surface responds in terms of a negative leading wave, followed with dispersive train of waves. For an upward bottom thrust, amplitude of the leading wave decays as the wave propagates, and no wave of permanent form evolves— instead, there appears a train of solitons. In this article, we show that our numerical scheme can produce the correct wave profiles, comparable with the analytical and experimental results of Hammack. Simulations using intermediate and slow bottom motions are also presented. In addition, we perform a simulation of a wave generated by submerged landslide, that compares well against previous numerical simulations. Via this simulation, we demonstrate that our scheme can incorporate a moving wet–dry boundary algorithm in the run-up simulation.     

基于非结构网格的TTGC有限元格式的实现及在超声速流动中的应用

基于非结构网格系统,实现了时空三阶精度的TTGC有限元格式,并在三阶TTGC格式上发展了基于人工粘性的激波捕捉技术。在非结构网格下,采用这种方法对若干典型的超声速流动问题(SOD激波管、马赫数为3的前台阶流动以及马赫数为8的高超声速圆柱流动)进行了验证计算。结果表明,TTGC格式分辨率高,在粗糙网格下能够准确的模拟超声速流场中的激波、接触间断等复杂流动现象,并且能有效的控制间断附近的数值色散现象。与传统的有限体积方法相比,本文实现的TTGC有限元格式在模拟超声速流动问题方面具有格式精度高、数值耗散小等优点。     

Fast Fourier homogenization for elastic wave propagation in complex media

In the context of acoustic or elastic wave propagation, the non-periodic asymptotic homogenization method allows one to determine a smooth effective medium and equations associated with the wave propagation in a given complex elastic or acoustic medium down to a given minimum wavelength. By smoothing all discontinuities and fine scales of the original medium, the homogenization technique considerably reduces meshing difficulties as well as the numerical cost associated with the wave equation solver, while producing the same waveform as for the original medium (up to the desired accuracy). Nevertheless, finding the effective medium requires one to solve the so-called “cell problem”, which corresponds to an elasto-static equation with a finite set of distinct loadings. For general elastic or acoustic media, the cell problem is a large problem that has to be solved on the whole domain and its resolution implies the use of a finite element solver and a mesh of the fine scale medium. Even if solving the cell problem is simpler than solving the wave equation in the original medium (because it is time and source independent, based on simple tetrahedral meshes and embarrassingly parallel) it is still a challenge. In this work, we present an alternative method to the finite element approach for solving the cell problem. It is based on a well-known method designed by H. Moulinec and P. Suquet in 1998 in structural mechanics. This iterative technique relies on Green functions of a simple reference medium and extensively uses Fast Fourier Transforms. It is easy to implement, very efficient and relies on a simple regular gridding of the medium. Through examples we show that the method gives excellent results, even, under some conditions, for discontinuous media.     

PARTITION OF UNITY FINITE ELEMENT METHOD FOR SHORT WAVE PROPAGATION IN SOLIDS

IntroductionIt is known that standard finite element procedure is unable to simulate the wavepropagation with high oscillations or gradients in space in the media with reasonableefficiency and accuracy due to the nature of polynomial interpolation approxi…     

Analysis of the pressure at a vertical barrier due to extreme wave run-up over variable bathymetry

The pressure load at a vertical barrier caused by extreme wave run-up is analysed numerically, using the conformal mapping method to solve the two-dimensional free surface Euler equations in a pseudospectral model. Previously this problem has been examined in the case of a flat-bottomed geometry. Here,the model is extended to consider a varying bathymetry. Numerical experiments show that an increasing step-like bottom profile may enhance the extreme run-up of long waves but result in a reduced pressure load.