共查询到18条相似文献，搜索用时 484 毫秒
1.
本文就波浪与结构物相互作用问题，提出了一个适用于高阶边界元应用的新的积分方程，并利用修改积分区域的方法得了适用于本积分方程的不规则频率消除方法。最后，通过数值计算对附加区域的选择、单元的离散做了研究 相似文献
2.
3.
潜型水平板水动力特性的数值研究 总被引：4，自引：0，他引：4
运用数值解析方法研究一个位于自由水面下的水平、刚性薄板的水动力特性。这种水平板可作为一种固定式浮型防波堤使用。假定流体为势流，波浪为线性时间谐波，利用格林定理，可将拉普拉斯控制方程简化为第二类Fredholm积分方程。采用线性分布的不连续单元，通过边界元法，可建立求解薄板源强分布的一组线性代数方程组，引入三维格林函数计算各个基本单元，并由积分方程直接求得附加质量，阻尼力和散射力。结果发现当板非常接近自由表面时，在一定频率范围内，附加质量为负值。本文另一主要特点是引入了GMRES快速迭代算法，这种算法可大量节省计算机内存，计算单元数高达数千个并且具有良好的收敛特性。 相似文献
4.
基于线性势流理论,利用高阶边界元法研究了规则波在三维局部渗透海床上的传播。根据Darcy渗透定律推导出渗透海床的控制方程,利用渗透海床顶部和海底处法向速度和压强连续条件得到渗透海床顶部满足的边界条件。根据绕射理论,利用满足自由水面条件的格林函数建立了求解渗透海床绕射势的边界积分方程,采用高阶边界元方法求解边界积分方程进而得到自由水面的绕射势和波浪在局部渗透海床上传播过程中幅值的变化情况。通过与已发表的波浪对圆柱形暗礁的时域全绕射结果对比,证明了本文建立的频域方法计算波幅的正确性和有效性。利用这一模型研究了三维矩形渗透海床区域上波浪的传播特性,并分析了入射波波长、海床渗透特性系数等参数对波浪传播的影响。 相似文献
5.
基于高阶边界元的三维数值波浪港池--波浪破碎的模拟 总被引：5，自引：1，他引：4
在势流理论的框架内,采用高阶边界元方法和混合欧拉-拉格朗日法,实现了对三维波浪破碎过程的数值模拟.数值模型使用可调节时间步长的基于二阶显式泰勒展开的混合欧拉-拉格郎日时间步进来求解自由表面的演化过程.在所使用的边界元方法中,采用16节点三次滑移四边形单元来表示,这种单元在单元内具有高阶的精度同时在单元之间具有良好的连续性.给出了孤立波的传播和周期性非线性波浪沿缓坡传播的计算结果,表明数值模型具有良好的稳定性. 相似文献
6.
深水中波浪与弱流对结构物的作用 总被引：4，自引：0，他引：4
对于深水中波浪和弱流对三维结构物的作用问题,本文提出了一个适用于高阶边界元应用的新的积分方程.基于小流速下格林函数和速度势的摄动展开,本方法避免了移动脉动源的计算,并将未知量限制在物体表面上,使计算速度大为提高,高阶边界元中的不同类型的柯西主值积分,分别采用间接和直接方法加以计算. 相似文献
7.
8.
9.
以方形立柱作为研究对象，对其在聚焦波作用下的波浪力特性进行研究。研究主要基于计算流体力学数值方法，对聚焦波作用下方形立柱在固定、单自由度纵荡、单自由度垂荡三种运动状态下所受波浪力进行研究。利用phase-inversion逆相位分解方法对其高阶力特性进行分析。首先对聚焦波生成方法、数值计算方法与参数设置进行简要介绍，其次展示了三种运动状态下立柱在聚焦波作用下所受波浪力，并对其进行分析，最后利用phase-inversion逆相位分解方法获得高阶波浪力成分并对其载荷特性进行分析讨论。研究发现逆相位分解方法可以有效分离波浪力高阶成分，不同运动状态下立柱所受波浪力会有显著不同。 相似文献
10.
一种新的波浪变形三维数值模式──0-1混合型边界元 总被引：4，自引：0，他引：4
为提高边界元法模拟三维波场波浪变形的数值计算精度,借鉴常数元和线性元剖分方式、波势函数及波势函数法向导数对单元节点设定的各自适应性,提出了一种新的单元剖分模式──0-1混合型边界元,以控制和减缓由于计算误差累计而造成的波浪数值计算上的“横向振动”,借此结合边界元法的分区模式可实现较大范围的波场线性波浪变形计算,并为时域内的波浪非线性变形计算提供时间步长的数值保证. 相似文献
11.
Removing irregular frequencies by a partial discontinuous higher order boundary element method 总被引：1，自引：0，他引：1
The phenomenon of irregular frequencies is a puzzle in the course of calculating the interaction of waves and structures by the boundary element method. To remove the irregular frequencies, the modified integral domain method is adopted, and continuous higher order elements and partial discontinuous higher order elements are used for discretization. By these means, the effects of the irregular frequencies are effectively removed. Effective strategies have been adopted to deal with singular integrals and nearly singular integrals at different situations. The numerical results of the horizontal wave force on a uniform cylinder in the first order and second order diffraction problems show that the present method has a good validity. At the same time, the influence of collocation parameter on accuracy of numerical results is examined in detail. 相似文献
12.
Teng Bin Li Yucheng
Associate Professor Dept. of Civil Engineering Dalian University of Technology Dalian Professor Dept. of Civil Engineering Dalian University of Technology Dalian 《中国海洋工程》1996,(3)
- For the discretization of higher order elements, this paper presents a modified integral domain method to remove the irregular frequencies inherited in the integral equation of wave diffraction and radiation from a surface-piercing body. The set of over-determined linear equations obtained from the method is modified into a normal set of linear equations by superposing a set of linear equations with zero solutions. Numerical experiments have also been carried out to find the optimum choice of the size of the auxiliary domain and the discretization on it. 相似文献
13.
A Fast Multipole Method (FMM) is developed as a numerical approach to the reduction of the computational cost and requirement memory capacity for a large in solving large-scale problems. In this paper it is applied to the boundary integral equation method (BIEM) for current diffraction from arbitrary 3D bodies. The boundary integral equation is discretized by higher order elements, the FMM is applied to avoid the matrix/vector product, and the resulting algebraic equation is solved by the Generalized Conjugate Residual method (GCR). Numerical examination shows that the FMM is more efficient than the direct evaluation method in computational cost and storage of computers. 相似文献
14.
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. 相似文献
15.
SUN Dapeng LI Yucheng TENG Bin
Associate Professor the State Key Laboratory of Coastal Offshore Engineering Dalian University of Technology Dalian
Professor the State Key Laboratory of Coastal Offshore Engineering Dalian University of Technology Dalian 《中国海洋工程》1999,(2)
For higher accuracy in simulating the transformation of three dimensional waves,in consid-eration of the advantages of constant panels and linear elements,a combined boundary elements is appliedin this research.The method can be used to remove the transverse vibration due to the accumulation ofcomputational errors.A combined boundary condition of sponge layer and Sommerfeld radiation condi-tion is used to remove the reflected waves from the computing domain.By following the water particle onthe water surface,the third order Stokes wave transform is simulated by the numerical wave flume tech-nique.The computed results are in good agreement with theoretical ones. 相似文献
16.
基于高阶边界元的三维数值波浪港池 总被引：9，自引：1，他引：8
初步建立了一个基于高阶边界元的三维数值波浪港池,港池具有造波和消波功能。采用高阶边界元16节点四边形单元和基于二阶显式泰勒展开的混合欧拉-拉格朗日时间步进求解带自由表面的完全非线性势流方程。模型中对于影响数值精度的问题作了细致的处理。数值计算结果表明本港池可以用来模拟非线性波浪的传播,具有很高的数值精度和稳定性。 相似文献
17.
Moo-Hyun Kim 《Applied Ocean Research》1992,14(6):353-370
The second-order difference-frequency wave forces on a large three-dimensional body in multi-directional waves are computed by the boundary integral equation method and the so-called FML formulation (assisting radiation potential method). Semi-analytic solutions for a bottom-mounted vertical circular cylinder are also developed to validate the numerical method. Difference-frequency wave loads on a bottom-mounted vertical cylinder and stationary four legs of the ISSC tension-leg platform (TLP) are presented for various combinations of incident wave frequencies and headings. These force quadratic transfer functions (QTF) can directly be used in studying slowly varying wave loads in irregular short-crested seas described by a particular directional spectrum. From our numerical results, it is seen that the slowly varying wave loads are in general very sensitive to the directional spreading function of the sea, and therefore wave directionality needs to be taken into account in relevant ocean engineering applications. It is also pointed out that the uni-directionality of the sea is not necessarily a conservative assumption when the second-order effects are concerned. 相似文献
18.
A fast time-domain method is developed in this paper for the real-time prediction of the six degree of freedom motions of a vessel traveling in an irregular seaway in infinitely deep water. The fully coupled unsteady ship motion problem is solved by time-stepping the linearized boundary conditions on both the free surface and body surface. A velocity-based boundary integral method is then used to solve the Laplace equation at every time step for the fluid kinematics, while a scalar integral equation is solved for the total fluid pressure. The boundary integral equations are applied to both the physical fluid domain outside the body and a fictitious fluid region inside the body, enabling use of the fast Fourier transform method to evaluate the free surface integrals. The computational efficiency of the scheme is further improved through use of the method of images to eliminate source singularities on the free surface while retaining vortex/dipole singularities that decay more rapidly in space. The resulting numerical algorithm runs 2–3 times faster than real time on a standard desktop computer. Numerical predictions are compared to prior published results for the transient motions of a hemisphere and laboratory measurements of the motions of a free running vessel in oblique waves with good agreement. 相似文献