共查询到20条相似文献，搜索用时 140 毫秒
1.
实际工程中存在大量的曲边界,因此在曲边界上的计算准确性可以考察出一个数值模型的实用价值。利用Beji的改进型Boussinesq方程建立了一个有限元方法的数值波浪模型。造波方面采用Fenton提出的非线性规则波浪解;在墙边界处,以求解法线方向和切线方向的速度和导数代替求解x、y方向的速度和导数,从而使边界条件直接适用、严格满足,保证了对曲边界计算的准确性。"重开始广义极小残量法"的使用保证了求解方程组的效率和精度,使造波和边界处理方法的有效性和准确性得到了合理地诠释。通过与试验数据、他人数值结果、解析解的比对,显示出该模型计算稳定、结果准确,真正体现出了有限元方法对曲边界适用的优势。 相似文献
2.
波浪对透空式结构物冲击作用的光滑粒子流体动力学数值模拟 总被引：3，自引：0，他引：3
基于光滑粒子流体动力学(SPH)方法的二维数值波浪水槽模拟了规则波对透空式结构物的冲击作用。通过黎曼解和CSPM相结合的方法对连续方程和动量方程进行了修正。对造波边界采用虚粒子法模拟。提出了一种耦合计算方法来消除造波边界附近的压力波动现象,在结构物边界处设置了适合的冲击边界条件。应用修正的SPH模型模拟了规则波对浪溅区结构物的冲击作用,给出了结构物附近流场和压力场的变化特征,并应用物理模型试验结果对数值模型进行了验证。 相似文献
3.
椭圆型缓坡方程的一个有效的有限元解 总被引：7，自引：0，他引：7
将海绵层消波的方法用于有限元方法中,提出了缓坡方程的一种有效的有限元求解方法.在应用有限元法求解椭圆型缓坡方程时,通过在方程中加摩阻项,并在入射边界(波浪由此边界进入计算域)处使用不连续单元,将绕射势从总势中分离出来,在边界上利用海绵层进行消波处理,有效地消除了由于引用放射边界条件引起的误差和数值反射现象. 相似文献
4.
非线性波浪波面追踪的一种新模式 总被引：1，自引：0，他引：1
基于Laplace方程的Green积分表达式和波面BemouUi方程所建立的非线性波动数学模型,是一个时域上具有初始值的边值问题,而精确地追踪自由表面的波动位置,给出波面运动瞬时的波面高度和波面势函数,是建立时域内非线性波浪数值模式的基础。本文采用0-1混合型边界元剖分计算域边界并离散Laplace方程的Green积分表达式,采用有限元剖分自由水面并推导满足自由表面非线性边界条件的波面有限元方程,联立计算域内以节点波势函数和波面位置高度的时间增量为未知量的线性方程组,通过时步内的循环迭代,给出每个时步上的波面位置和波面势函数,从而建立了一种新的非线性波浪波面追踪模式。数值造波水槽内的波浪试验表明,其数值模拟结果具有良好的计算精度。 相似文献
5.
缓坡方程是描述近岸波浪运动较好的数学模型之一。在发展的自适应有限元求解缓坡方程的基础上,采用迭代求解的方法,确定波浪相对于边界的入射方向,从而对边界条件进行改进,建立了求解缓坡方程的数值计算模型。典型算例表明,考虑波浪相对于边界的入射角度后,模型可以更好地模拟吸收波浪边界,同时对多向波对双突堤的绕射进行了模拟研究,与试验结果比较表明,所建立的数值计算模型能够适用于多向不规则波传播过程的模拟研究。 相似文献
6.
在快速模拟波浪运动的谱方法基础上,引入造波边界,建立了模拟波浪产生和运动的二阶计算模型。采用摄动展开方法简化了带有造波边界的水波运动问题,将速度势分解,得到了满足造波边界和自由面边界的速度势的一般解,运用快速Fourier变换和时间积分,建立了模拟波浪产生和运动的数学模型。基于该模型,采用不同的数值造波条件,模拟了波浪的产生问题;考虑了波浪的初始运动问题;通过把数值结果与物理实验的比较,验证了波浪计算模型的有效性。 相似文献
7.
8.
不规则波Boussinesq型方程的造波,消波和反射 总被引：5，自引：1，他引：4
对前人提出的造波、消波和反射边方法分析表明,其方法是极浅水波近似,不适用于任意水深的水域,本文就任意水深变化Ｂｏｕｓｓｉｎｅｓｑ型方程,提出了不规则波新的造波原理、方法和消波边界及部分反射边界波动方程,试验表明,本文提出的造波、消波和反射方程有效而可靠的。 相似文献
9.
基于高阶边界元的三维数值波浪港池 总被引：8，自引：1，他引：8
初步建立了一个基于高阶边界元的三维数值波浪港池，港池具有造波和消波功能。采用高阶边界元16节点四边形单元和基于二阶显式泰勒展开的混合欧拉-拉格朗日时间步进求解带自由表面的完全非线性势流方程。模型中对于影响数值精度的问题作了细致的处理。数值计算结果表明本港池可以用来模拟非线性波浪的传播，具有很高的数值精度和稳定性。 相似文献
10.
非静压波浪模型是一种近年来得到迅速发展且有望在工程计算中获得广泛应用的波浪模型。为了在非静压波浪模型中实现无反射造波,在连续性方程中添加了质量源项进行域内源造波,并采用阻尼消波方法进行数值消波。数值算例结果表明,所建立的无反射造波方法可以较好地避免非静压波浪数学模型中的二次反射,实现长时间稳定波浪模拟。 相似文献
11.
The behavior of a highly deformable membrane to ocean waves was studied by coupling a nonlinear boundary element model of the fluid domain to a nonlinear finite element model of the membrane. The hydrodynamic loadings induced by water waves are computed assuming large body hydrodynamics and ideal fluid flow and then solving the transient diffraction/radiation problem. Either linear waves or finite amplitude waves can be assumed in the model and thus the nonlinear kinematic and dynamic free surface boundary conditions are solved iteratively. The nonlinear nature of the boundary condition requires a time domain solution. To implicitly include time in the governing field equation, Volterra's method was used. The approach is the same as the typical boundary element method for a fluid domain where the governing field equation is the starting point. The difference is that in Volterra's method the time derivative of the governing field equation becomes the starting point.The boundary element model was then coupled through an iterative process to a finite element model of membrane structures. The coupled model predicts the nonlinear interaction of nonlinear water waves with highly deformable bodies. To verify the coupled model a large scale test was conducted in the OH Hinsdale wave Research Laboratory at Oregon State University on a 3-ft-diameter fabric cylinder submerged in the wave tank. The model data verified the numerical prediction of the structure displacements and of the changes in the wave field.The boundary element model is an ideal modeling technique for modeling the fluid domain when the governing field equations is the Laplace equation. In this case the nonlinear boundary element model was coupled with a finite element model of membrane structures, but the model could have been coupled with other finite element models of more rigid structures, such as a pontoon floating breakwater. 相似文献
12.
Lin Mingchung Hsiao Sungshan Hsu Yungcheng
Professor Dept. of Naval Architecture Ocean Engineering National Taiwan University Taiwan. Doctor Course Student Dept. of Naval Architecture Ocean Engineering National Taiwan University Taiwan 《中国海洋工程》1994,(3)
-Wave refraction-diffraction due to a large ocean structure and topography in the presence of a 'current are studied numerically. The mathematical model is the mild-slope equation developed by Kirby (1984). This equation is solved using a finite and boundary element method. The physical domain is devid-ed into two regions: a slowly varying topography region and a constant water depth region. For waves propagating in the constant water depth region, without current interfering, the mild- slope equation is then reduced to the Helmholtz equation which is solved by boundary element method. In varying topography region, this equation will be solved by finite element method. Conservation of mass and energy flux of the fluid between these two regions is required for composition of these two numerical methods. The numerical scheme proposed here is capable of dealing with water wave problems of different water depths with the main characters of these two methods. 相似文献
13.
二维各向异性介质中三分量波动方程有限元法模拟 总被引：1，自引：0，他引：1
利用六面体单元和三线性插值函数 ,推导出在任意弹性各向异性介质中三维三分量波动方程所满足的有限元方程 ,同时给出也适用于矩形单元和双线性插值的二维三分量各向异性介质中的波动方程模拟的一般性公式 ,最后讨论两个数值模拟的结果。 相似文献
14.
修正型缓坡方程的有限元模型 总被引：1，自引：1，他引：0
与缓坡方程相比,修正型缓坡方程增加了地形曲率项和坡度平方项,从而提高了数值求解的复杂性。本文将计算域划分为内域和外域,内域为水深变化区域,使用修正型缓坡方程,其中的地形曲率项和坡度平方项可用有限单元各节点的水深信息和单元插值函数表示,外域为水深恒定区,速度势满足Helmholtz方程,通过内外域的边界匹配建立有限元方程,并用高斯消去法求解。进而分别模拟了波浪传过Homma岛和圆形浅滩的变形,其结果与相关的解析解和实验数据吻合良好,证明了本文有限元模型的正确性。同时,通过与实验数据的对比也明显看出,在地形坡度较陡的情况下,修正型缓坡方程较缓坡方程具有更高的计算精度。 相似文献
15.
The finite element method(FEM) is employed to analyze the resonant oscillations of the liquid confined within multiple or an array of floating bodies with fully nonlinear boundary conditions on the free surface and the body surface in two dimensions.The velocity potentials at each time step are obtained through the FEM with 8-node quadratic shape functions.The finite element linear system is solved by the conjugate gradient(CG) method with a symmetric successive overelaxlation(SSOR) preconditioner.The waves at the open boundary are absorbed by the combination of the damping zone method and the Sommerfeld-Orlanski equation.Numerical examples are given by an array of floating wedgeshaped cylinders and rectangular cylinders.Results are provided for heave motions including wave elevations,profiles and hydrodynamic forces.Comparisons are made in several cases with the results obtained from the second order solution in the time domain.It is found that the wave amplitude in the middle region of the array is larger than those in other places,and the hydrodynamic force on a cylinder increases with the cylinder closing to the middle of the array. 相似文献
16.
In this paper, based on the linear wave theory, the interaction of short-crested waves with a concentric dual cylindrical system with a partially porous outer cylinder is studied by using the scaled boundary finite element method (SBFEM), which is a novel semi-analytical method with the advantages of combining the finite element method (FEM) with the boundary element method (BEM). The whole solution domain is divided into one unbounded sub-domain and one bounded sub-domain by the exterior cylinder. By weakening the governing differential equation in the circumferential direction, the SBFEM equations for both domains can be solved analytically in the radial direction. Only the boundary on the circumference of the exterior porous cylinder is discretized with curved surface finite elements. Meanwhile, by introducing a variable porous-effect parameter G, non-homogeneous materials caused by the complex configuration of the exterior cylinder are modeled without additional efforts. Comparisons clearly demonstrate the excellent accuracy and computational efficiency associated with the present SBFEM. The effects of the wide range wave parameters and the structure configuration are examined. This parametric study will help determine the various hydrodynamic effects of the concentric porous cylindrical structure. 相似文献
17.
The finite element method (FEM) is employed to analyze the resonant oscillations of the liquid confined within multiple or an array of floating bodies with fully nonlinear boundary conditions on the free surface and the body surface in two dimensions. The velocity potentials at each time step are obtained through the FEM with 8-node quadratic shape functions. The finite element linear system is solved by the conjugate gradient (CG) method with a symmetric successive overelaxlation (SSOR) preconditioner. The waves at the open boundary are absorbed by the combination of the damping zone method and the Sommerfeld-Orlanski equation. Numerical examples are given by an array of floating wedge- shaped cylinders and rectangular cylinders. Results are provided for heave motions including wave elevations, profiles and hydrodynamic forces. Comparisons are made in several cases with the results obtained from the second order solution in the time domain. It is found that the wave amplitude in the middle region of the array is larger than those in other places, and the hydrodynamic force on a cylinder increases with the cylinder closing to the middle of the array. 相似文献
18.
A numerical model to compute wave field is developed. It is based on the Berkhoff diffraction-refraction equation, in which an energy dissipation term is added, to take into account the breaking and the bottom friction phenomena. The energy dissipation function, by breaking and by bottom friction, is introduced in the Berkhoff equation to obtain a new equation of propagation.The resolution is done with the hybrid finite element method, where lagrangians elements are used. 相似文献
19.
M. Arockiasamy 《Ocean Engineering》1983,10(5):303-312
Response of a compliant platform to irregular waves is determined using finite element method. The tower is idealized by 2-D beam elements with an elastic support at the guy lines location. The flexural characteristics of the beam correspond to the four corner members of the trusses. The guying system is modelled by an axial element with linear load deformation characteristics. A computer program based on the linearized Morison's equation and the linear (Airy) wave theory, is developed to calculate the total force based on the storm wave height data at different levels of the structure. The response of the structure to random waves is based on the spectral approach. The direct and cross spectral densities of the generalized wave forces are determined and used to obtain the spectral densities of the generalized modal coordinates and mean square response at each level. Possible extension of the method is indicated to compute the evolutionary response to nonstationary wave forces. 相似文献
20.
A numerical simulation model is developed to investigate wave profiles and velocities in harbor channels. The Galerkin Technique is coupled with the finite element method in the flow model. An important feature is the addition of terms in the momentum equation to account for the vertical accelerations. The model can be used by ocean engineers to analyze the effects of varying channel geometry on wave oscillations and velocities in a channel connecting a harbor to the ocean. 相似文献