波流联合作用下的物体上的波浪绕射力

本文研究了在波流联合作用下的任意三维物体的波浪绕射问题,建立了求解物体上的波浪绕射力的时域解理论模型。本文证明了对于任意三维物体,如果在自由面条件简化时,假设定常兴波势为有限量阶,则在波浪绕射力公式中,不存在含有绕射势的水线积分项。从数学上证明了水线积分项反映了自由面条件中定常兴波势和不定常势之间的耦合效应。在此基础上,本文提出了在波流联合作用情况下,作用在物体上的绕射力与辅助辐射势之间的关系式,从而不须求解绕射势便可得到绕射力。     

用耦合元方法计算二维物体的波浪荷载

本文中采用边界积分方程和有限元耦合的数值方法,计算了二维波浪绕射和辐射问题。给出圆形和长方形浮体、潜体在波场中的水动力系数和所受的波浪力,并详细讨论了计算方法。耦合元方法是近年提出的,用于无限区域中波浪的绕射和辐射问题时单元少,计算方便,经济可行,在海洋工程和波能利用的计算流体中有一定的推广意义。该方法也可用于三维问题。     

非线性波浪绕射问题的数值计算

本文给出非线性波浪绕射问题的计算模型,并用时间步进和边界元法求解。该模型能适用于严重非线性情形,且能满足初始条件和辐射条件。作为对该模型的检验,我们对孤立波绕过穿刺水面的直立圆柱的算例进行了计算,并同Isaacson的解析解、李宝元的差分解进行了比较,结果吻合很好。此外我们还用数值方法同时模拟造波机造波和求该波与直立圆柱的相互作用。结果表明,数值波浪水池能够以尽可能短的暂态过程,模拟波浪与物体的相互作用,提供良好的初始条件。数值波浪水池还能用于任意水深、任意水底形状下的相互作用问题。本文的计算模型和算法具有机时省,计算稳定的优点。     

基于速度势分离法的波物相互作用完全非线性数值模拟

该文基于势流理论和高阶边界元方法,建立了开敞水域波浪与三维物体相互作用的完全非线性数值模型。采用速度势分离技术,将总速度势分离成入射势和散射势,入射势通过理论解给定。采用混合欧拉-拉格朗日方法追踪流体瞬时水面,引进虚拟函数耦合求解波浪力和运动方程,自由水面网格仅在初始时刻生成一次,随后采用弹簧近似法在不改变网格节点排列的情况下对瞬时水面进行网格重构。利用所建立模型模拟完全非线性波浪对固定圆柱的绕射作用,并对规则浪三倍频作用下漂浮圆柱的共振响应进行了研究。     

波浪对圆柱外接圆弧型防波堤绕射的解析计算

在微幅波绕射理论的基础上,使用特征函数展开法,推导了圆柱外接圆弧型防波堤复合结构的水波绕射的解析解,并据此计算了该结构所受到的波浪载荷以及绕射波面分布。计算结果表明:与单一圆柱相比,外接圆弧型防波堤的存在可明显降低作用于圆柱的波浪载荷和圆柱表面的波浪爬高;采用外接堤表面密实或透空而内圆柱表面透空的结构形式可以取得圆柱防浪的较佳效果;结构透空系数、波浪入射角度、外接堤张角、外堤半径与水深比、水深与波长比等因素的变化对绕射波浪作用均存在一定影响。     

大直径圆柱体上的二阶波浪力

在二阶波浪力的计算中,二阶绕射势满足的非齐次自由条件,其非齐次项在整个自由面上不是绝对可积的,给解析或数值求解造成了困难。本文推广了周清甫的工作,对有限水深引进了一个“局部扰动源”,用此局部扰动源在自由面上的分布,构造了满足非齐次自由面条件的二阶势的特解。并进一步证明,此特解在远处收敛得足够快,只剩下外传的自由波,后者在无穷远处满足通常的Som-mevfeld条件。 按这个方法,本文得到了大直径圆柱体在有限水深中的完整的二阶绕射势,计算了二阶波浪力,并与实验数据和现有的计算结果进行了比较。比较表明,在波长稍大的地方,本文的结果有明显改善。     

近水面三维物体运动的时域计算

本文提出了求解近水面三维物体作任意运动时受到的水动力和产生的波浪的压力波方法。运动可以是大振幅的,但自由面条件是线性的。整个问题由自由面上分布偶极子和物面上布置源汇来模拟后,就可以分解为自由面上布有压力的压力波解和无限流场的绕流解。本文导出了水面布有压力为一般函数时的解析解,然后利用特殊函数的各种性质和公式积出了其中的四重无限积分,大大地简化了。三维的计算主要是计算机速度和内存空间的限制,故本文的方法有较大的实用价值。作为例子,本文计算了突然起动的圆球和回转体作定常运动时的水动力以及它幅射的波浪。结果与其它数值和PMM试验结果比较证实了本方法的合理性和计算的准确性。     

二阶波浪力的切片理论与潜体上的垂向定常力

二阶波浪力的三维理论基本上是对无航速的结构物建立起来的，航速的影响不甚清楚，一些理论问题，例如二阶绕射势辐射条件的提法，直到最近才得以澄清，即使如此，用三维方法完整地求取物体上的二阶波浪力仍十分费时，本文中对有航速的细长物体建立了二阶波浪力的切片理论，旨在提供理论上合理，实用上相对简便的二阶力计算方法，作为完整的二阶波浪力的一个组成部分，文中给出了细长潜体在不同潜深和浪向角时所受的垂向常力及变化趋     

椭圆余弦波对V形防波堤绕射波浪力的解析计算

应用椭圆余弦波理论,给出了V形薄壁防波堤的浅水波绕射理论解,对现有Airy微幅波理论进行了有效拓展。通过对V形防波堤的浅水波浪力进行试算,揭示了椭圆余弦波对防波堤的作用规律。计算结果表明:椭圆余弦波理论计算的V形防波堤的最大无量纲波浪力明显大于相同浅水条件下Airy微幅波理论的对应结果,该方法有效拓展了无限长直立薄壁堤的反射波理论。此外,浅水波入射角、V形堤张角和防波堤臂长与水深比等参数的变化均将对V形堤的波浪作用产生相关影响,而V形堤的实际绕射波浪力幅值将随浅水波特征参数值的增大而增大。     

基于时域格林函数的有航速船舶与波浪相互作用的数值模拟

基于深水时域格林函数的数值处理方法及开发的计算程序,采用三维时域数值模拟方法,在线性势流理论框架内对深海波浪中航行船舶的辐射问题和绕射问题进行了求解,研究探讨了航行船舶与波浪相互作用问题.研究中对不同航行速度的Wigley型船水动力系数和波浪力进行了数值计算,通过与试验及文献计算结果的比较,验证了三维时域数值计算程序的有效性.     

FREE-SURFACE WAVES AND FAR WAKES GENERATED BY A FLOATING BODY IN A VISCOUS FLUID

1 .　INTRODUCTIONThewavemotionduetoamovingdisturbancehasbeenoflong standinginterest .Butthecorre spondingwakestructuredidnotreceivemuchat tentionuntilthesyntheticapertureradarcollectedanumberofimagesofV likewakepatternsgener atedbyships (Milgram 1988;Brownetal.1989) .Thewakestructuresarethoughttobeclassifiedin tothreecategories :theKelvinwavewakes ,thetrubulentwakesandtheinternalwavewakes(Griffin 1988;Lydenetal.1988;GuandPhillips1994 ) .TheKelvinwakesarethefar fieldwavepatternsgenerate…     

Iterative Rankine HOBEM analysis of hull-form effects in forward-speed diffraction problem

A time-domain numerical algorithm based on the higher-order boundary element method and the iterative time-marching scheme is proposed for seakeeping analysis. The ship waves generated by a hull advancing at a constant forward speed in incident waves and the resultant diffraction forces acting on the hull are computed to investigate the hull-form effects on the hydrodynamic forces. A rectangular computational domain travelling at ship's speed is considered. An artificial damping beach for satisfying the radiation condition is installed at the outer portion of the free surface except the downstream side. An iterative time-marching scheme is employed for updating both kinematic and dynamic free-surface boundary conditions for numerical accuracy and stability. The boundary integral equation is solved by distributing higher-order boundary elements over the wetted body surface and the free surface. The hull-form effects on the naval hydrodynamics are investigated by comparing three different Wigley models. Finally, the corresponding unsteady wave patterns and the wave profiles around the hulls are illustrated and discussed.     

Tracking methods for free surface and simulation of a liquid droplet impacting on a solid surface based on SPH

With some popular tracking methods for free surface, simulations of several typical examples are carried out under various flow field conditions.The results show that the Smoothed Particle Hydrodynamics (SPH) method is very suitable in simulating the flow problems with a free surface.A viscous liquid droplet with an initial velocity impacting on a solid surface is simulated based on the SPH method, and the surface tension is considered by searching the free surface particles, the initial impact effect is co...     

UNSTEADY FREE-SURFACE WAVES GENERATED BY BODIES IN A VISCOUS FLUID

BRIEFINTRODUCTIONOFTHEPAPER :　Inthefirstpartofthisthesis ,thefundamentalsolutionsforthesingularStokesandOseenflowsinanun boundedfluidarederivedinauniversalformwhichinvolvestheHamiltonian ,Hessian ,andLaplacianoperators ,andelementaryfunctions .Inthesecondpart,theinteractionofunsteadylow Reynolds numberflowswithafreesurfaceisinvestigatedanalytically .Thedisturbedflows ,generatedbysubmergedbodiesmovingverticallydownwardsawayfromthesurfaceofthefluid ,aregovernedbytheunsteadyStokeseq…     

An iterative Rankine boundary element method for wave diffraction of a ship with forward speed

A 3-D time-domain seakeeping analysis tool has been newly developed by using a higher-order boundary element method with the Rankine source as the kernel function. An iterative time-marching scheme for updating both kinematic and dynamic free-surface boundary conditions is adopted for achieving numerical accuracy and stability. A rectangular computational domain moving with the mean speed of ship is introduced. A damping beach at the outer portion of the truncated free surface is installed for satisfying the radiation condition. After numerical convergence checked, the diffraction unsteady problem of a Wigley hull traveling with a constant forward speed in waves is studied. Extensive results including wave exciting forces, wave patterns and pressure distributions on the hull are presented to validate the efficiency and accuracy of the proposed 3-D time-domain iterative Rankine BEM approach. Computed results are compared to be in good agreement with the corresponding experimental data and other published numerical solutions.     

Time domain simulations of radiation and diffraction by a Rankine panel method

The radiation and the diffraction of a ship with a forward speed are studied by using a time domain Rankine panel method. The free surface conditions are linearized onto an undisturbed free surface based on the double body flow. The linearized body boundary condition is applied on the mean wetted hull surface. The fluid domain boundary is discretized by a collection of quadric panels. The unknown quantities, including the free surface elevation, the normal flux over the free surface and the potential on the fluid domain boundary, are determined at each time step. The numerical results are compared with experimental data and other numerical solutions, showing satisfactory agreements.     

Free-surface flow past a submerged cylinder

This study employs a Lagrangian frame numerical method to investigate two-dimensional free-surface flow induced by a submerged moving cylinder. This method combines the advantages of vortex methods and boundary integral methods, and is capable of capturing the complex motion of free surface and vortices. A series of computations are performed to investigate the effects of Froude number, the depth of submergence and still water depth on the flow motion. Free surface deformation, wake vortex and hydrodynamic forces are studied.

     

Surface and interfacial gravity waves due to a disturbance steadily moving in a two-layer inviscid fluid

The generation and interaction of surface and interfacial gravity waves due to an submerged source moving in a two-layer fluid are investigated analytically for two-dimensional cases. The fluid is assumed to be inviscid and incompressible. The density of each of the two layers is constant. Two different boundary conditions are considered for the upper surface. The upper fluid of finite depth is bounded above by a rigid lid or a free surface. Based on the assumption of small-amplitude waves, a linear system is established. The integral solutions for the free-surface and interfacial elevations are obtained by means of the Fourier transform. Then the corresponding asymptotic representations are derived for far-field waves by the residue theorem. The critical Froude numbers for the existence of far-field waves are derived for the two-layer system bounded above by a rigid lid or a free surface. The effect of different upper boundary conditions on the wave generation are discussed.

     

Wavy properties and analytical modeling of free-surface flows in the development of the multi-domain method

The wavy(oscillatory both in space and in time) properties of free-surface flows due to presence of floating bodies are analyzed within the framework of the potential-flow theory by assuming that the fluid is perfect and flow irrotational. A so-called new multi-domain method has been developed based on the fluid domain division by an analytical control surface surrounding bodies and the application of different methods adapted in the external and internal domains. In the analytical domain external to the control surface, the fundamental solution satisfying the linear boundary condition on the free surface associated with a point singularity(often called Green function and referred here as point solution) is applied to capture all wavy features of free-surface flows extending horizontally to infinity. Unlike classical studies in which the control surface is discretized, the unknown velocity potential and its normal derivatives are expressed by expansions of orthogonal elementary functions. The velocity potential associated with each elementary distribution(elementary solutions) on the control surface can be obtained by performing multi-fold integrals in an analytical way. In the domain internal to the control surface containing the bodies, we could apply different methods like the Rankine source method based on the boundary integral equations for which the elementary solutions obtained in the external domain playing the role of Dirichlet-to-Neumann operator close the problem.     

核电厂重水池内自由液面之全非线性运动模拟

本论文利用除奇点方法建立一个数值分析模型。用来模拟核电厂重水池受地震影响，进行受迫运动时，其自由液面的运动情形。文中将模拟地震进行方式分为升降，横荡等运动模型。设定初始值包含振幅，频率，重水池宽度及流体深度等。并将模拟之结果做讨论，计算结果之比较显示使用本法可正确地模拟受到地震时，核能发电厂的重水池自由液面的运动情形，比较计算结果显示，使用除奇点技巧所建立之数值模型相对于有限元素法（FEM）能以较少的计算时间及更简单的程式架构。模拟出相同或更好的结果。