首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 125 毫秒
1.
给出了一种联合瞬态格林函数和Rankine源进行有航速浮体时域水动力分析的混合——边界元方法。在三维混合边界元方法中,通过一个匹配面将流体域划分为内域和外域,在内域中使用Rankine源以模拟直壁或非直壁船体及线性或非线性自由面条件,在外域中使用瞬态格林函数以满足自由面条件和远方辐射条件。使用该方法计算了一个有航速潜没圆球的波浪力,和解析结果的比较证明了该方法的正确性。进一步给出了一个有航速Wigley船的水动力结果,计算结果稳定,没有外传波向内反射的现象发生。  相似文献   

2.
在快速模拟波浪运动的谱方法基础上,引入造波边界,建立了模拟波浪产生和运动的二阶计算模型。采用摄动展开方法简化了带有造波边界的水波运动问题,将速度势分解,得到了满足造波边界和自由面边界的速度势的一般解,运用快速Fourier变换和时间积分,建立了模拟波浪产生和运动的数学模型。基于该模型,采用不同的数值造波条件,模拟了波浪的产生问题;考虑了波浪的初始运动问题;通过把数值结果与物理实验的比较,验证了波浪计算模型的有效性。  相似文献   

3.
采用边界条件造波法与人工阻尼消波法开展波浪数值模拟时,存在自由面抬升现象,导致波浪模拟精度略有损失。为解决上述问题,提出了一种基于单方向阻尼源项的数值消波方法,通过采用VOF模型在大、小两种波陡条件下开展波浪数值模拟进行了验证,并进一步研究了阻尼系数与密度依赖性对波浪模拟精度的影响。结果表明,单方向阻尼源项能够显著降低平均自由面高度,有效抑制平均自由面抬升,进而提高波浪模拟精度。阻尼系数及密度依赖性的优取值与波陡条件有关,且阻尼系数的影响大于密度依赖性。  相似文献   

4.
船舶螺旋桨尾流场的数值分析   总被引:16,自引:1,他引:16  
利用基于速度势的低阶面元法计算船舶螺旋桨的尾流场。采用计算较为简捷的关于扰动速度势的基本积分微分方程,并采用双曲面形状的面凶以消除面元间的缝隙。Newton-Raphson迭代过程被用来在桨叶随边满足压力Kutta条件,使桨叶面上表面的压力在随边有良好的一致性。在计算面元的影响系数时,应用了Morino导出的解析计算公式,加快了数值计算的速度。从解面元法的基本积分方程得到的偶极强度和源汇强度,直接求得尾流场的速度分布。  相似文献   

5.
本文提出了一个生成自由面的三维船体贴体网格的数值方法。文中以静水面以上的船体型线作为描述非定常自由面的计算网格域,进一步完善了我们以往开发的无法考虑兴波与静水面以上的船体型线相互作用对船舶粘性流动影响的网格生成方法。静水面以上的初始自由面网格分别沿船面法线η及切线ζ方向的曲线长的指数进行分布,即可极方便地调节自由面在船面及静水面附近的网格间距,也保证了在数值离散三维Poisson方程时具有相当好的收敛性。这种方法的特点是计算量相当小,尤其能极方便地控制网格在船面簿粘性层中及静水面附近的分布。以Series60船模作为计算算例,带自由面的网格分布是相当满意的  相似文献   

6.
消除"不规则频率"的非连续高阶元方法   总被引:2,自引:0,他引:2  
针对使用边界元法计算波浪与结构物相互作用时所出现的“不规则频率”现象,采用连续高阶元和部分非连续高阶元对通过修改积分区域所获得的边界积分方程进行离散,有效地消除了“不规则频率”现象的发生。波浪作用下的截断圆柱所受到的水平波浪力和垂向波浪力的数值计算结果验证了该方法的有效性,同时考虑了非连续单元配置点的选择及单元划分数目对消除效果的影响。  相似文献   

7.
鉴于钱塘江涌潮水流复杂多变和带有自由表面的情况,采用VOF方法追踪自由表面,考虑各向异性,用修正的K吨模型封闭Reynolds方程,用施主一受主法求解流体体积函数控制方程,计算了涌潮冲击丁坝时的压力场和自由面状况,并与现场实测值进行了比较,取得了良好的效果,为涌潮的小尺度数值模拟提供了借鉴的经验。  相似文献   

8.
引入Stokes有限振幅水波的多色波作用机制 ,采用二层层化海洋模型 ,分析了层化海洋中非线性波对结构的散射作用问题。提出了二层海洋中二阶水波散射条件的数学表达式 ,进而推导了二层海洋中二阶波对一般型圆柱散射作用的格林积分解式 ,并给出了对解式中自由面及分层界面无限积分的递推算式。  相似文献   

9.
三维自由面流动模拟中GPU并行计算技术   总被引:1,自引:0,他引:1       下载免费PDF全文
MPS(Moving Particle Semi-implicit)法能够有效地处理溃坝、晃荡等自由面大变形流动问题。在三维MPS方法中,粒子数量的急剧增加会导致其计算效率的降低并限制其在大规模流动问题中的应用。基于自主开发的MPS求解器MLParticleSJTU,本文对求解过程中耗时最多的邻居粒子搜寻和泊松方程求解两个模块采用了GPU并行加速,详细探讨了CPU+GPU策略。以三维晃荡和三维溃坝这两种典型的自由面大变形流动为例,比较了CPU+GPU相对于MLParticle-SJTU串行求解时的加速情况,结果表明CPU+GPU在邻居粒子和泊松方程这两个模块中的加速比最高能达到十倍左右。此外,采用CPU+GPU并行能够较准确地模拟溃坝、晃荡等自由面大变形问题。  相似文献   

10.
配置积分方法在GPS卫星轨道确定中的应用   总被引:1,自引:0,他引:1  
重点介绍了配置(Collocation)积分方法的原理和步骤。比较了其他轨道积分方法和Collocation积分方法在积分卫星轨道时的区别,阐明了Collocation方法在计算卫星速度、加速度等高阶导数方面的优越性,并通过采用固定步长、变化多项式阶数和固定多项式阶数、变化步长的方法,定性和定量地分析了不同参数的选择对GPS卫星轨道积分的影响。  相似文献   

11.
Based on the Rankine source, this paper proposed a time-domain method for analyzing the three-dimensional wave–structure interaction problem in irregular wave. A stable integral form of the free-surface boundary condition (IFBC) is employed to update the velocity potential on the free surface. A multi-transmitting formula, with an artificial wave speed, is used to eliminate the wave reflection for radiation condition on the artificial boundary. An effective multi-transmitting formula, coupled with damping zone method, is further used to analyze the irregular wave diffraction at the artificial boundary. We investigate hydrodynamic forces on floating structure and compare our solution to the frequency-domain solution. It is shown that long time simulation can be done with high stability and the numerical results agree well with the solution obtained under the frequency domain. The efficiency of the proposed multi-transmitting formula and the coupled methods for radiation boundary make them promising candidates in studying the irregular water wave problem in time domain.  相似文献   

12.
在悬沙输运的数值模拟中,初始场的准确给定至关重要。目前诸多确定初始场的方案均存在一定的缺陷,初始场的准确性有待进一步提高。本文基于一个三维悬沙输运伴随同化模型,通过孪生实验和实际实验,对模型初始场进行了伴随法反演研究。在孪生实验中,首先验证了初始场的相对重要性;其次,探讨了初始场的反演结果对优化算法、初始猜测值、卫星遥感数据数量、同化时间窗口宽度和背景流场误差的敏感性;最后,比较了伴随法和插值法重构初始场的能力。孪生实验结果表明:最速下降法对初始场的优化反演效果要优于三种共轭梯度法和有限记忆BFGS法;初始场的反演效果对初始猜测值、卫星遥感数据数量和背景流场误差不敏感,而对同化窗口宽度较为敏感;与插值法相比,伴随法是重构模型初始场更有效的手段。实际实验中,在杭州湾海域同化典型的小潮时期和大潮时期的GOCI卫星遥感资料所得表层悬沙浓度数据,优化反演了初始场。实际实验结果表明:数据同化后,得到了更符合实际的最优初始场,表明伴随法是实现初始场优化反演的有效手段。该研究对进一步改进悬沙输运模型的初始化方案具有一定的参考价值,也对其他数值模型的初始化方案具有一定的借鉴价值。  相似文献   

13.
A potential based panel method for the hydrodynamic analysis of 2-D hydrofoils moving beneath the free surface with constant speed without considering cavitation is described. By applying Green's theorem and the Green function method, an integral equation for the perturbation velocity potential is obtained under the potential flow theory. Dirichlet type boundary condition is used instead of Neumann type boundary condition. The 2-D hydrofoil is approximated by line panels which have constant source strength and constant doublet strength distributions. The free surface condition is linearized and the method of images is used for satisfying this free surface condition. All the terms in fundamental solution (Green function) of perturbation potential are integrated over a line panel. Pressure distribution, lift, residual drag and free surface deformations are calculated for NACA4412, symmetric Joukowski and van de Vooren profile types of hydrofoil. The results of this method show good agreement with both experimental and numerical methods in the literature for the NACA4412 and symmetric Joukowski profile types. The lift and residual drag values of the van de Vooren profile are also presented. The effect of free surface is examined by a parametric variation of Froude number and depth of submergence.  相似文献   

14.
《Ocean Engineering》1999,26(4):343-361
A potential based panel method for the hydrodynamic analysis of 2-D hydrofoils moving beneath the free surface with constant speed without considering cavitation is described. By applying Green's theorem and the Green function method, an integral equation for the perturbation velocity potential is obtained under the potential flow theory. Dirichlet type boundary condition is used instead of Neumann type boundary condition. The 2-D hydrofoil is approximated by line panels which have constant source strength and constant doublet strength distributions. The free surface condition is linearized and the method of images is used for satisfying this free surface condition. All the terms in fundamental solution (Green function) of perturbation potential are integrated over a line panel. Pressure distribution, lift, residual drag and free surface deformations are calculated for NACA4412, symmetric Joukowski and van de Vooren profile types of hydrofoil. The results of this method show good agreement with both experimental and numerical methods in the literature for the NACA4412 and symmetric Joukowski profile types. The lift and residual drag values of the van de Vooren profile are also presented. The effect of free surface is examined by a parametric variation of Froude number and depth of submergence.  相似文献   

15.
The objective of the present work is to discuss the implementation of an active wave generating–absorbing boundary condition for a numerical model based on the Volume Of Fluid (VOF) method for tracking free surfaces. First an overview of the development of VOF type models with special emphasis in the field of coastal engineering is given. A new type of numerical boundary condition for combined wave generation and absorption in the numerical model VOFbreak2 is presented. The numerical boundary condition is based on an active wave absorption system that was first developed in the context of physical wave flume experiments, using a wave paddle. The method applies to regular and irregular waves. Velocities are measured at one location inside the computational domain. The reflected wave train is separated from the incident wave field in front of a structure by means of digital filtering and subsequent superposition of the measured velocity signals. The incident wave signal is corrected, so that the reflected wave is effectively absorbed at the boundary. The digital filters are derived theoretically and their practical design is discussed. The practical use of this numerical boundary condition is compared to the use of the absorption system in a physical wave flume. The effectiveness of the active wave generating–absorbing boundary condition finally is proved using analytical tests and numerical simulations with VOFbreak2.  相似文献   

16.
Ship floating condition in regular waves is calculated. New equations controlling any ship's floating condition are proposed by use of the vector operation. This form is a nonlinear optimization problem which can be solved using the penalty function method with constant coefficients. And the solving process is accelerated by dichotomy. During the solving process, the ship's displacement and buoyant centre have been calculated by the integration of the ship surface according to the waterline. The ship surface is described using an accumulative chord length theory in order to determine the displacement, the buoyancy center and the waterline. The draught forming the waterline at each station can be found out by calculating the intersection of the ship surface and the wave surface. The results of an example indicate that this method is exact and efficient. It can calculate the ship floating condition in regular waves as well as simplify the calculation and improve the computational efficiency and the precision of results.  相似文献   

17.
Nan Xie  Dracos Vassalos   《Ocean Engineering》2007,34(8-9):1257-1264
The purpose of the present paper is to develop a potential-based panel method for determining the steady potential flow about three-dimensional hydrofoil under free surface. The method uses constant-strength doublets and source density distribution over the foil body surface and thereby Dirichlet-type boundary condition is used instead of Neumann-type condition. On the undisturbed free surface source density is used to meet the free surface condition that is linearised in terms of double-body model approach and is discretised by a one-side, upstream, four-point finite difference operator. After solving the doublets on the foil and sources on the free surface, the numerical results of pressure, lift and resistance coefficients and also wave profiles can then be calculated for different Froude number and depth of submergence to demonstrate the influence of free surface and aspect ratio effects on performance of the hydrofoil.  相似文献   

18.
The purpose of this study is to investigate the phenomenon of free surface flow impacting on elastic structures, which is a research topic of great interest in ship and ocean engineering. A series of quasi two-dimensional experiments on dam-break with an elastic plate are conducted. The main features of free surface flow impacting on elastic structures including large impacting force, structural vibration, violent free surface flow, are investigated. The coupled FDM–FEM method developed by the authors is applied for numerical simulation of such dam-break problem. Extensive analysis and discussion based on the comparisons between experimental data and numerical results are made and presented in this paper.  相似文献   

19.
A nonlinear analysis is made for determining the two-dimensional unsteady potential-flow characteristics about a wing subject to wing-in-ground effect (WIG) operating above progressive water waves. The dynamic boundary condition requiring the constant pressure and the kinematic boundary condition prescribing the continuity in the vertical velocity are satisfied on the undisturbed free surface. The boundary conditions imposed on the free surface are linear, but the kinematic boundary condition satisfied on the foil surface is nonlinear. Through the derivation and evaluation of the time-domain Green's functions for two-dimensional singularities above a free surface, the influence of water waves on the lift performance of the two-dimensional WIG is addressed using the discrete vortex method. Furthermore, the roll-up of the wake vorticity is considered. The comparison of present work with other numerical results available in the literature shows the validation of the present approach.  相似文献   

20.
In this paper, the three-dimensional water exit of a sphere with different vertical velocities is investigated numerically using the lattice Boltzmann method (LBM). In this method, the liquid-gas two-phase flow is simplified as a single-phase free surface flow. To capture the free surface, a mass tracking algorithm is incorporated into the LBM. The gravity as a body force is introduced in the form of calculating the equilibrium distribution with an altered velocity, while the surface tension is neglected. Besides, the employed bounce-back boundary conditions are used for a moving sphere. What’s more, the Wall-Adapting Local Eddy (WALE) viscosity model is employed to capture the turbulent structures of the flow and stabilize the simulation. The accuracy of the numerical results is demonstrated through comparisons with the previous numerical and experimental results in the literature. The results show that the spike height is significantly influenced under the Froude number (Fr) below 4.12 and slightly affected under the Fr varying from 4.12 to 8.24. After the sphere exits water totally, the evolution of the free surface waterfall can be described as two phases and becomes more intense with the Froude number increasing. The non-uniform distribution of velocity results in the breaking of the free surface after the sphere completely exits the water. Moreover, the Reynolds number greatly affects the wake dynamics and hydrodynamics acting on the sphere when it moves beneath the water surface.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号