非均匀水流中非线性波传播的数值模拟

以一种考虑波流相互作用的新型{Boussinesq}型方程为控制方程组， 采用五阶{Runge}-{Kutta}-{England}格式离散时间积分，采用七点 差分格式离散空间导数，并通过采用恰当的出流边界条件，从而建立了非均匀水流中非线性 波传播的数值模拟模型. 通过对均匀水流与水深水域内和潜堤地形上存在弱流或强流时波浪 传播的数值模拟，说明模型能有效地反映水流对波浪传播的影响.     

随机波浪作用下近岸波流场的数值模拟

结合近岸波浪抛物型缓坡模型和近岸波流场模型，对近岸不规则波浪及其破碎后所产生的流 场进行了数值模拟. 在不规则波浪场的模拟中，采用JONSWAP波浪谱对入射单向不规则波浪 要素按等分频率法进行离散，应用考虑波浪不规则性和破碎效应的抛物型缓坡方程对波浪场进行 数值模拟，并基于抛物型缓坡方程中的波浪势函数等参数计算波浪辐射应力，以波浪辐射应 力为主要动力因素基于近岸流数学模型对近岸波浪破碎所产生的近岸流场进行数值模拟，并 对数值模拟结果进行了验证. 模拟结果表明该模型可有效地用于研究波浪破碎产生的近岸波 流场.     

黄茅海波-流共同作用下的三维悬沙数值模拟

针对黄茅海河口区波流相互作用显著、三维空间结构明显及泥沙运动复杂等特点,联合第3代海浪模式SWAN,建立了波流共同作用下的三维悬沙数学模型.在水流计算中,考虑了波浪产生的辐射应力影响,波浪依赖的表面风应力影响及波浪影响下的底部剪切应力,并采用k-kl紊流闭合模型提供垂向涡粘系数和垂向扩散系数.悬沙扩散方程中的源函数利用切应力法确定,泥沙沉降速度考虑絮凝作用,从而提高黄茅海悬沙场数值模拟精度.通过含沙量实测资料验证,模拟值与实测值符合较好.在模型验证的基础上,讨论了不考虑波浪与考虑波浪两种情况下黄茅海的悬沙分布特征,为工程实际研究提供了依据.     

考虑表面效应纳米压电双晶中波的频散分析

基于非局部应变梯度理论探究了考虑表面弹性和表面残余应力的纳米压电双晶中波的频散特性,压电双晶的上下压电层暴露在电场之中并且整体沉积在粘弹性基底之上.利用哈密顿原理和正弦剪切理论推导了控制方程,利用含非局部参数和长度尺度参数的尺度依赖本构关系得到了运动方程,带入谐波解求解相应的特征方程.数值揭示了表面弹性和表面残余应力、尺度参数和波数以及粘弹性基底对压电双晶的作用规律.研究表明,表面效应的存在对压电纳米双晶频率特性的研究至关重要,尺度参数和波数对频散特性具有耦合作用,弹性系数、阻尼系数和压电层厚度对频率的作用表现出区域性.     

一种求解复波数域波动弥散方程的通用数值方法

波动弥散方程的求解是分析弹性波或声波问题时的一种比较常见的过程,求解所得的曲线即为弥散曲线。其方程形式一般是关于波数与频率的一个二元超越方程,当波数为复数时,方程实际为三元超越方程。这类方程没有解析解,由于方程的变量和参数涉及复数域,一般的数值求解方法也难以处理。本文提出了一种基于模值收敛性的数值方法,求解一般波动弥散方程复波数域内的解。首先论述了本方法的理论依据,然后介绍了其求解步骤,最后将本方法应用于求解复波数域波动弥散方程的具体算例。通过对比可以发现:本方法广泛有效地解决了复波数域内波动弥散曲线的求解问题,适用于不同形式波动弥散方程的求解。     

功能梯度压电压磁材料中裂纹对SH波的散射

研究无限大功能梯度压电/压磁复合材料中裂纹对SH波的散射问题.为了便于分析,假设材料性质沿着裂纹的法线方向是指数变化.利用Fourier余弦变化,将问题转化为对偶积分方程的求解,此对偶积分方程采用Copson方法求解.然后求得应力强度因子、电位移强度因子、磁通量强度因子的解析表达式,最后数值算例给出了材料参数、入射角及波数对标准动应力强度因子的影响.     

移动载荷作用下浮冰的非线性动力学响应

本文考虑非线性、惯性和阻尼的影响, 研究了任意深度二维理想流体顶部浮冰的振动. 对相关的拟微分算子进行展开并将非线性项保留至三阶后, 完全非线性问题被简化为仅与自由面上的变量相关的三阶截断模型. 为了验证简化模型的准确性, 重点关注了自由孤立波解. 在不考虑阻尼的情况下, 采用多重尺度方法推导了三阶非线性薛定谔方程(NLS), 利用该方程预测了任意水深下原始欧拉方程中自由波包型孤立波解的存在性及三阶截断模型的准确性. 相比于Dinvay等所提出的二阶模型, 三阶截断模型的优势在于其对应的三阶NLS具有准确的非线性项系数, 能够在最小相速度附近更好地模拟冰层的动力学响应. 进一步地对自由孤立波解进行数值计算, 数值结果表明三阶截断模型在分岔曲线和孤立波波形上均与完全欧拉方程吻合良好, 准确性高于二阶截断模型. 基于三阶截断模型, 探究了匀速局域化载荷作用下的浮冰非线性动力学响应并将时间依赖解与实验测量数据进行比较, 数值计算结果与实验记录吻合良好.      

?????????е?????????????????????

悬浮液中纤维的运动和旋转行为在许多现代工业领域都非常重要. Jeffery解析得到 了牛顿流中椭球型纤维的动力学演化方程，但方程只有在平面简单剪切和拉伸流的条 件下才有解析解. 采用数值求解时，由于方程中$\cot \theta $的存在会出现奇异性. 采用随体坐标的方法，消除了方程的奇异性，算法的有效性通过与解析解的比较得到了 验证，并为复杂流场下纤维的动力学数值计算提供了一个可行的方法.     

波-流相互作用的缓坡方程及其波作用量守恒

当表面波从开阔海域传播至近岸水域时，普遍的波一流相互作用经受着海底的强烈影响．运用水波Hamilton变分原理，建立了近岸水域任意水深变化海底上波一流相互作用的缓坡方程．它可包含波、流和水深一般变化的二阶效应，约化为某些典型的缓坡型方程．据此得出广义程函方程，并且证明该缓坡方程的波作用量守恒．     

移动荷载作用下饱和土地基中的波动特性分析

基于Biot波动方程,经过Fourier变换和逆Fourier变换后可获得波数-频率域以及时间-空间域的解析解。通过数值分析的手段研究了移动荷载作用下饱和多孔弹性地基中波的传播特性。重点就弥散曲线、多谱勒效应、波的成分和动力响应频率等几个特性进行了分析,发现饱和土地基由于比弹性地基多了一项流体介质,波动特性明显差异于弹性介质。     

A compact numerical algorithm for solving the time‐dependent mild slope equation

The mild slope equation has been widely used to describe combined wave refraction and diffraction. In this study, a new numerical algorithm is developed to solve the time‐dependent mild slope equation in a second‐order hyperbolic form. The numerical algorithm is based on a compact and explicit finite difference method that is second‐order accurate in both time and space. The algorithm has the similar structure to the leap‐frog method but is constructed on three time levels for the second‐order time derivative term. The numerical model has the capability of simulating transient wave motion by correctly predicting the speed of wave energy propagation, which is important for the real‐time forecast of the arrival time of storm waves generated in the far field. The model is validated against analytical solution for wave shoaling and experimental data for combined wave refraction and diffraction over a submerged elliptic shoal on a slope (Coastal Eng. 1982; 6 :255). Lastly, the realistic scale Homma's island (Geophys. Mag. 1950; 21 :199) is studied with the use of various wave periods of T = 720s, T = 120 s, and T = 24 s. These wave periods correspond to long, intermediate, and short waves for the given topography, respectively. Comparisons are made between numerical results and existing analytical solutions in terms of the wave amplification around the island, which serves as the indicator for the potential wave runup. Excellent agreements are obtained. The model runs on a PC (Pentium IV 1.8GHz) and the computer capacity allows the computation of a mesh system up to 3000 × 3000, which is equivalent to about 150 × 150 waves or a large area of 540km × 540km for a wave train with the period of T = 60 s. Copyright 2004 John Wiley & Sons, Ltd.     

基于相平均方法的折射绕射联合波浪模型

近岸带波浪运动的研究具有很重要的工程意义，近年来已获得了较丰硕的研 究成果并发展了许多波浪模型，而基于不同理论的波浪模型往往具有特定的适用性. 在海岸 工程中应用比较广泛的一类波浪模型以波能(波作用量)守恒为基本依据，如SWAN模型. 该 类模型在实际工程中已经得到了大量的应用，但该类模型未计及波浪绕射效应，成为其突出 的缺陷之一. 如何对模型做适当的改进，使之适用于波浪绕射的模拟，从而在原有基础上拓 广模型的应用范围是一项具有实际意义的研究工作. 该文采用波能(波作用量)守恒方程描 述近岸带波浪运动，通过引入绕射因子，得到折射、绕射联合波浪模型，从而拓广了模型的 应用范围. 通过实际算例验证，表明所建立的模型计及了波浪折射、绕射作用，对相平 均波浪模型在波浪绕射效应模拟方面的改进具有一定的意义.     

Numerical simulation of pollutant transport acted by wave for a shallow water sea bay

A numerical model of pollutant transport acted by water waves on a shallow‐water mild‐slope beach is established in this study. The numerical model is combined with a wave propagation model, a multiple wave‐breaking model, a wave‐induced current model and a pollutant convection–dispersion model. The wave propagation model is based on the higher‐order approximation of parabolic mild‐slope equation which can be used to simulate the wave refraction, diffraction and breaking in a large area of near‐shore zone combined with the wave‐breaking model. The wave‐induced current model is established using the concept of the radiation stress and considering the effect of bottom resistance caused by waves. The numerical model is verified by laboratory experiment results of regular and irregular waves over two mild beaches with different slopes. The numerical results agree well with experimental results. The numerical model has been applied in the near‐shore zone of Bohai bay in China. It is concluded that pollutant transport parallel to the shoreline due to the action of waves, which will induce serious pollution on the beach. Copyright © 2005 John Wiley & Sons, Ltd.     

基于VPM-THINC/QQ模型的波浪高保真模拟

为实现波浪传播的高保真数值模拟,采用包含单元均值和点值(volume-average/point-value method,VPM)的有限体积法求解纳维-斯托克斯方程和具有二次曲面性质和高斯积分的双曲正切函数(THINC method with quadratic surface representation and Gaussian quadrature,THINC/QQ)方法来重构自由面,建立以开源求解库OpenFOAM底层函数库为基础的VPM-THINC/QQ模型. 在本模型中添加推板造波法实现波浪的产生功能,采用松弛法实现消波功能,构建高精度黏性流数值波浪水槽. 分别采用VPM-THINC/QQ模型和InterFoam求解器(OpenFOAM软件包中广泛使用的多相流求解器)开展规则波的数值模拟,重点探究网格大小和时间步长等因素对波浪传播过程的影响,定量地分析波高衰减程度;为验证本模型的适应性,对长短波进行模拟. 结果表明,在相同网格大小或时间步长条件下,VPM-THINC/QQ模型的预测结果与参考值吻合较好,波高衰减较少,且无相位差,在波浪传播过程的模拟中呈现出良好的保真性. 本文工作 为波浪传播的模拟研究提供了一种高精度的黏性数值波浪水槽模型.      

具有精确色散性的非线性波浪数学模型

以完全非线性的自由表面边界条件为基础,以波面升高\eta和自由表面速度势\phi _\eta为待求变量,建立了新的波浪方程．方程在色散性上是完全精确的,非线性近似至三阶．与缓坡方程相比较,两者都具有精确的色散性,但该方程属于非线性模型,可模拟波浪的非线性效应,且适用于不规则波．方程的特点是属于微分-积分方程,对如何处理方程中积分项进行了讨论,并数值模拟了不同周期的线性波和二阶Stokes波,也模拟了波群的非线性演化,以对模型进行验证．      

Improvement of moving particle semi‐implicit method for simulation of progressive water waves

Precise simulation of the propagation of surface water waves, especially when involving breaking wave, takes a significant place in computational fluid dynamics. Because of the strong nonlinear properties, the treatment of large surface deformation of free surface flow has always been a challenging work in the development of numerical models. In this paper, the moving particle semi‐implicit (MPS) method, an entirely Lagrangian method, is modified to simulate wave motion in a 2‐D numerical wave flume preferably. In terms of consecutive pressure distribution, a new and simple free surface detection criterion is proposed to enhance the free surface recognition in the MPS method. In addition, a revised gradient model is deduced to diminish the effect of nonuniform particle distribution and then to reduce the numerical wave attenuation occurring in the original MPS model. The applicability and stability of the improved MPS method are firstly demonstrated by the calculation of hydrostatic problem. It is revealed that these modifications are effective to suppress the pressure oscillation, weaken the local particle clustering, and boost the stability of numerical algorithm. It is then applied to investigate the propagation of progressive waves on a flat bed and the wave breaking on a mild slope. Comparisons with the analytical solutions and experimental results indicate that the improved MPS model can give better results about the profiles and heights of surface waves in contrast with the previous MPS models. Copyright © 2017 John Wiley & Sons, Ltd.     

A meshless numerical wave tank for simulation of nonlinear irregular waves in shallow water

Time domain simulation of the interaction between offshore structures and irregular waves in shallow water becomes a focus due to significant increase of liquefied natural gas (LNG) terminals. To obtain the time series of irregular waves in shallow water, a numerical wave tank is developed by using the meshless method for simulation of 2D nonlinear irregular waves propagating from deep water to shallow water. Using the fundamental solution of Laplace equation as the radial basis function (RBF) and locating the source points outside the computational domain, the problem of water wave propagation is solved by collocation of boundary points. In order to improve the computation stability, both the incident wave elevation and velocity potential are applied to the wave generation. A sponge damping layer combined with the Sommerfeld radiation condition is used on the radiation boundary. The present model is applied to simulate the propagation of regular and irregular waves. The numerical results are validated by analytical solutions and experimental data and good agreements are observed. Copyright © 2008 John Wiley & Sons, Ltd.