共查询到17条相似文献，搜索用时 125 毫秒
1.
基于高阶边界元的三维数值波浪港池--波浪破碎的模拟 总被引：5，自引：1，他引：4
在势流理论的框架内,采用高阶边界元方法和混合欧拉-拉格朗日法,实现了对三维波浪破碎过程的数值模拟.数值模型使用可调节时间步长的基于二阶显式泰勒展开的混合欧拉-拉格郎日时间步进来求解自由表面的演化过程.在所使用的边界元方法中,采用16节点三次滑移四边形单元来表示,这种单元在单元内具有高阶的精度同时在单元之间具有良好的连续性.给出了孤立波的传播和周期性非线性波浪沿缓坡传播的计算结果,表明数值模型具有良好的稳定性. 相似文献
2.
波浪漫滩是近岸波浪的小尺度运动,在实际海域的波生流数值计算中通常被忽略。本文基于Boussinesq方程的FUNWAVE模式,分别采用波浪漫滩边界、固壁边界、海绵边界进行Haller波浪港池物理模型实验的数值模拟,比较三种边界计算结果与实验观测数据的误差,检验波浪漫滩边界对波生流数值计算的影响;然后设计了多种周期、波高的波生流数值模拟试验,分析多种波浪入射条件下波浪漫滩边界对近岸波生流数值计算的影响。结果表明,波浪漫滩对邻近区域波生流有明显影响,漫滩边界下的波生流计算结果更接近实验观测值,在近岸波生流数值模型中引入波浪漫滩边界可以提高波生流计算精度。 相似文献
3.
基于势流理论和时域高阶边界元方法,建立了三维完全非线性数值波浪水槽模型.利用源造波法产生入射波浪,应用五阶斯托克斯波理论给定波浪速度;采用混合欧拉-拉格朗日方法追踪流体瞬时水面,将二阶泰勒级数展开法应用于更新下一时间步的波面和速度势;通过加速势的方法准确计算自由水面速度的法向导数和物面速度势的时间导数.对完全非线性波浪进行了模拟,得到了稳定的波形.当波浪非线性较小时,与四阶Runge- Kutta法(RK4)计算结果和五阶斯托克斯波理论解均吻合良好;随着波浪非线性的增大,计算结果误差逐渐增大.通过数值试验分析,在满足精度要求的基础上,本方法计算时间略大于四阶Runge- Kutta法的四分之一,大大减少了计算量. 相似文献
4.
5.
一种新的波浪变形三维数值模式──0-1混合型边界元 总被引：4，自引：0，他引：4
为提高边界元法模拟三维波场波浪变形的数值计算精度,借鉴常数元和线性元剖分方式、波势函数及波势函数法向导数对单元节点设定的各自适应性,提出了一种新的单元剖分模式──0-1混合型边界元,以控制和减缓由于计算误差累计而造成的波浪数值计算上的“横向振动”,借此结合边界元法的分区模式可实现较大范围的波场线性波浪变形计算,并为时域内的波浪非线性变形计算提供时间步长的数值保证. 相似文献
6.
无限水深聚焦波完全非线性数值模拟 总被引：1，自引：1，他引：0
基于势流理论提出一种新的高阶边界元方法对无限水深的聚焦波浪进行完全非线性数值模拟.自由水面满足完全非线性边界条件,模拟波浪的非线性效果可以达到更高阶.利用镜像原理,建立一种全新的格林函数应用到无限水深的数值波浪水槽中,以致于两无限深水槽侧壁的积分可以被排除.为了产生相应的入射波和吸收出流波浪,一个由点源组成的造波装置被布置于计算域内,同时人工阻尼层被用来吸引出流波浪,由波浪聚焦的方法得到极限波浪.通过开展线性和完全非线性聚焦波浪的数值实验及与理论解对比,验证本数值模型可以用来模拟无限深水域的极限波浪,且在出流边界没有反射. 相似文献
7.
非线性波浪波面追踪的一种新模式 总被引：1，自引：0，他引：1
基于Laplace方程的Green积分表达式和波面BemouUi方程所建立的非线性波动数学模型,是一个时域上具有初始值的边值问题,而精确地追踪自由表面的波动位置,给出波面运动瞬时的波面高度和波面势函数,是建立时域内非线性波浪数值模式的基础。本文采用0-1混合型边界元剖分计算域边界并离散Laplace方程的Green积分表达式,采用有限元剖分自由水面并推导满足自由表面非线性边界条件的波面有限元方程,联立计算域内以节点波势函数和波面位置高度的时间增量为未知量的线性方程组,通过时步内的循环迭代,给出每个时步上的波面位置和波面势函数,从而建立了一种新的非线性波浪波面追踪模式。数值造波水槽内的波浪试验表明,其数值模拟结果具有良好的计算精度。 相似文献
8.
应用基于势流理论的时域高阶边界元方法,建立一个完全非线性的三维数值波浪水槽,通过实时模拟推板造波运动的方式产生波浪。通过混合欧拉-拉格朗日方法和四阶Runge-Kutta方法更新自由水面和造波板的瞬时位置。利用所建模型分别模拟了有限水深波和浅水波,与试验结果、相关文献结果和浅水理论结果吻合较好,且波浪能够稳定传播。系统地讨论造波板的运动圆频率、振幅和水深等对波浪传播和波浪特性的影响,并对波浪的非线性特性进行分析,研究发现造波板运动频率、运动振幅以及水深均将对波浪形态和波浪非线性产生显著影响。结果为真实水槽造波机的运动控制以及波浪生成试验提供了依据,便于实验室设置更合理的参数来准确模拟不同条件下的波浪。 相似文献
9.
10.
本文在时域非线性数值波浪水槽中，研究了不同风速条件下极端波浪的特性。采用推板造波的方式生成非线性波浪，基于Jeffrey遮蔽理论将风压项引入自由面动力学边界条件来模拟风压作用，通过高阶边界元法和混合欧拉-拉格朗日时间步进法来求解初边值问题。通过与已发表的聚焦波群实验结果对比验证了该数值模型的准确性，并研究了风压对极端波浪的最大波高、聚焦位置的偏移和波浪谱的演变等波浪性质的影响。本文进一步在数值波浪水槽中引入均匀水流，来模拟风生流对波浪演变的影响。结果表明，风压的存在会少量增大极端波浪的最大波高，波浪的聚焦和解焦过程伴随着明显的能量传递，并且风生流进一步导致了波浪聚焦位置的偏移。 相似文献
11.
Simulation of Fully Nonlinear 3-D Numerical Wave Tank 总被引：6，自引：0，他引：6
A fully nonlinear numerical wave tank (NWT) has been simulated by use of a three-dimensional higher order bouodary element method (HOBEM) in the time domain. Within the frame of potential flow and the adoption of simply Rankine source, the resulting boundary integral equation is repeatedly solved at each time step and the fully nonlinear free surface boundary conditions are integrated with time to update its position and boundary values. A smooth technique is also adopted in order to eliminate the possible saw-tooth numerical instabilities. The incident wave at the uptank is given as theoretical wave in this paper. The outgoing waves are absorbed inside a damping zone by spatially varying artificial damping on the free surface at the wave tank end. The numerical results show that the NWT developed by these approaches has a high accuracy and good numerical stability. 相似文献
12.
13.
Fully nonlinear wave-body interactions for stationary surface-piercing single and double bodies are studied by a potential-theory-based fully nonlinear 2D numerical wave tank (NWT). The NWT was developed in time domain by using boundary element method (BEM) with constant panels. MEL free surface treatment and Runge–Kutta fourth-order time integration with smoothing scheme was used for free-surface time simulation. The acceleration-potential scheme is employed to obtain accurate time derivative of velocity potential. Using the steady part of nonlinear force time histories, mean and a series of higher-harmonic force components are calculated and compared with the experimental and numerical results of other researchers. The slow-decaying second-harmonic vertical forces are investigated with particle velocities and corresponding body pressure. Typical patterns of two-body interactions, shielding effect, and the pumping/sloshing modes of water column in various gap distances are investigated. The pumping mode in low frequencies is demonstrated by the comparison of velocity magnitudes. 相似文献
14.
Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φn-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated. 相似文献
15.
《Ocean Engineering》2006,33(8-9):983-1006
Nonlinear waves and forces induced by a wedge-shape wave maker were simulated in a potential-theory-based fully nonlinear 2D Numerical Wave Tank (NWT). The NWT is developed in a time domain by using Boundary Element Method (BEM) including Mixed Eulerian–Lagrangian method (MEL) and Runge–Kutta 4th-order (RK4) integration as a time marching process. For ensuring accurate nonlinear free surface both material-node and semi-Lagrangian approach are independently developed for crosschecking. The acceleration-potential scheme is used for obtaining accurate time derivative of velocity potential. The developed NWT is utilized to calculate water particle velocity and a series of higher-harmonic force components on the wave maker. The added-mass and radiation-damping coefficients of the wave maker are also obtained from the least-square method. The simulation results are compared with the experimental and numerical results of other researchers. To compare the relative importance of free-surface and body-surface nonlinearities, a body nonlinear formulation is independently developed. Force by body nonlinear method is in good agreement with fully nonlinear result in case of low body-stroke frequency. 相似文献
16.
Gilbert R.W. Zedler E.A. Grilli S.T. Street R.L. 《Oceanic Engineering, IEEE Journal of》2007,32(1):236-248
In this paper, we report on the use of a numerical wave tank (NWT), based on fully nonlinear potential flow (FNPF) equations, in driving simulations of flow and sediment transport around partially buried obstacles. The suspended sediment transport is modeled in the near-field in a Navier-Stokes (NS) model using an immersed-boundary method and an attached sediment transport simulation module. Turbulence is represented by large eddy simulation (LES). The NWT is based on a higher order boundary element method (BEM), with an explicit second-order time stepping. Hence, only the NWT boundary is discretized. The solution for the velocity potential and its derivatives along the boundary is obtained in the BEM, which subsequently provides a solution at any required internal point within the domain. At initial time, the NS-LES model domain is initialized with the 3-D velocity field provided by the NWT and driven for later time by the pressure gradient field obtained in the NWT. Incident wave fields, as specified in the NWT to drive sediment transport, can be arbitrary. Applications are presented here for single frequency waves, such as produced by a harmonic piston wavemaker in the laboratory, and modulated frequency wave groups. The feasibility of coupling the irrotational flow and NS solutions is demonstrated. 相似文献
17.
In this paper, motion response of a moored floating structure interacting with a large amplitude and steep incident wave field is studied using a coupled time domain solution scheme. Solution of the hydrodynamic boundary value problem is achieved using a three-dimensional numerical wave tank (3D NWT) approach based upon a form of Mixed-Eulerian–Lagrangian (MEL) scheme. In the developed method, nonlinearity arising due to incident wave as well as nonlinear hydrostatics is completely captured while the hydrodynamic interactions of radiation and diffraction are determined at every time step based on certain simplifying approximations. Mooring lines are modelled as linear as well as nonlinear springs. The horizontal tension for each individual mooring line is obtained from the nonlinear load-excursion plot of the lines computed using catenary theory, from which the linear and nonlinear line stiffness are determined. Motions of three realistic floating structures with different mooring systems are analyzed considering various combinations of linear and approximate nonlinear hydrodynamic load computations and linear/nonlinear mooring line stiffness. Results are discussed to bring out the influence and need for consideration of nonlinearities in the hydrodynamics and hydrostatics as well as the nonlinear modelling of the line stiffness. 相似文献