Simulation of flows around cylinders by a Lagrangian vortex scheme

A method is presented for calculating flows around cylinders of arbitrary shape, which may be in an array. Vorticity is discretised and is moved in a Lagrangian time-stepping procedure by the vortex-in-cell method. The zero-velocity surface condition is maintained at each time step by representing the surface by a vortex sheet through a boundary integral calculation. Vorticity diffusion is simulated through random walks. Results are presented for an isolated circular cylinder in steady and oscillatory incident flows.     

Wave forces on a large, horizontal submerged cylinder

A numerical boundary integral equation method combined with a non-linear time stepping procedure is used for the calculation of wave forces on a large, submerged, horizontal circular cylinder. As the method is based on potential theory, all computations are performed in the inertia dominated domain, that is, for small Keulegan-Carpenter numbers. Computations are carried out for the Eulerian mean current under wave trough level equal to zero. When the cylinder is moved towards the sea bed the computations show that the inertia coefficients increase significantly, which is associated with a blockage effect. Furthermore, the effect of the wave steepness is reduced when the submergence of the cylinder is increased. In the vicinity of the free water surface the vertical inertia coefficient is highly dependent upon the wave steepness, which tends to reduce it, whereas the horizontal inertia coefficient is only slightly dependent on the wave steepness. Computations are also carried out for cylinder diameters comparable with the wave length. Finally, inertia coefficients computed by the present method are compared with some analytical results by Ogilvie [(1963), First and second order forces on a cylinder submerged under a free surface. J. Fluid Mech. 16, 451–472]. As long as the assumptions leading to Ogilvie's theory are fulfilled (cylinder radius small compared to the wave length), the results are quite similar.     

Numerical Investigation of Run-ups on Cylinder in Steep Regular Wave

The run-up on offshore structures induced by the steep regular wave is a highly nonlinear flow with a free surface.This article focuses on the investigation of the steep regular wave run-up on a single vertical cylinder by solving the Navier-Stokes equations. A numerical wave tank is established based on the open-source package to simulate the wave scattering induced by a vertical cylinder. The VOF method is applied to capture the large deformation and breaking of the free surface. The numerical model is validated by experimental results. The relative wave run-ups on the front face and the back face along the centerline of a cylinder are analyzed. The changes of the relative run-ups with the wave steepness, the relative diameter and the relative depth are studied. It is found that the relative run-ups on the front face and the back face of the cylinder depend mainly on the wave steepness and the relative diameter,while the dependence on the relative depth is weak. The empirical formulae are proposed to calculate the relative run-ups in terms of the wave steepness of incident regular waves and the relative diameter of a cylinder.     

Third order wave force on axisymmetric bodies

The third order triple-frequency wave load on fixed axisymmetric bodies by monochromatic waves is considered within the frame of potential theory. Waves are assumed to be weak non-linearity and a perturbation method is used to expand velocity potentials and wave loadings into series according to a wave steepness of kA. Integral equation method is used to compute velocity potentials up to second order in wave steepness. The third order triple-frequency wave loads are computed by an indirect method and an efficient method is applied to form the third order forcing term on the free surface quickly. The method can be used to compute third order triple-frequency surge force, heave force and pitch moment on any revolution bodies with vertical axes. The comparison with Malenica and Molin's results is made on surge force on a uniform cylinder, and comparison with experimental results is made on third order surge force, heave force and pitch moment on a truncated cylinder. More numerical computations are carried out for third order forces and moments on a uniform cylinder, truncated cylinders and a hemisphere.     

A Finite Element Solution of Wave Forces on Submerged Horizontal Circular Cylinders

In this paper, a numerical model is established for estimating the wave forces on a submerged horizontal circular cylinder. For predicting the wave motion, a set of two-dimensional Navier-Stokes equations is solved numerically with a finite element method. In order to track the moving non-linear wave surface boundary, the Navier-Stokes equations are discretized in a moving mesh system. After each computational time step, the mesh is modified according to the changed wave surface boundary. In order to stabilize the numerical procedure, a three-step finite element method is applied in the time integration. The water sloshing in a tank and wave propagation over a submerged bar are simulated for the first time to validate the present model. The computational results agree well with the analytical solution and the experimental data.Finally, the model is applied to the simulation of interaction between waves and a submerged horizontal circular cylinder.The effects of the KC number and the cylinder depth on the wave forces are studied.     

双矩形浮子波能装置辐射问题的解析方法

在有限水深下1个漂浮在水中的矩形浮子和1个淹没在水下的浮子构成双矩形波能装置模型。基于特征函数展开法求解了线性入射波作用下双矩形浮子波能装置的辐射问题,得出了双矩形浮子辐射速度势的1种新解析式,然后根据Haskind关系由入射势和辐射势来计算波浪激励力,并且采用数值方法对相同算例进行了计算,得到了完全一致的结果,从而证明这种方法是正确的。研究了在不同工况下的波浪激励力和系统的水动力学系数变化的规律。     

Interaction of uniform current with a circular cylinder submerged below an ice sheet

The problem of a uniform current passing through a circular cylinder submerged below an ice sheet is considered. The fluid flow is described by the linearized velocity potential theory, while the ice sheet is modelled through a thin elastic plate floating on the water surface. The Green function due to a source is first derived, which satisfies all the boundary conditions apart from that on the body surface. Through differentiating the Green function with respect to the source position, the multipoles are obtained. This allows the disturbed velocity potential to be constructed in the form of an infinite series with unknown coefficients which are obtained from the boundary condition. The result shows that there is a critical Froude number which depends on the physical properties of the ice sheet. Below this number there will be no flexural waves propagating to infinity and above this number there will be two waves, one on each side of the body. When the depth based Froude number is larger than 1, there will always be a wave at far upstream of the body. This is similar to those noticed in the related problem and is different from that in the free surface problem without ice sheet. Various results are provided, including the properties of the dispersion equation, resistance and lift, ice sheet deflection, and their physical features are discussed.     

The Nonlinear Wave Force on a Circular Cylinder in Shallow Water

Based on the 1st order cnoidal wave theory, the nonlinear wave diffraction around a circular cylinder in shallow water is studied in this paper. The equation of the wave surface around the cylinder is formulated and by using this formula the wave surface elevation on the cylinder surface can be obtained. In this paper, the formula for calculating the cnoidal wave force on a circular cylinder is also derived. For the wave conditions which are often encountered in practical engineering designs, the ratios of the nonlinear wave forces to the linear wave forces are calculated, and the results are plotted in this paper for design purposes. In order to verify the theoretical results, model tests are conducted. After comparing the test results with the theoretical ones, it is concluded that, in shallow water, for the case of T g / d~(1/2) > 8-10 and H / d > 0.3, the cnoidal wave theory should be used to calculate the wave action on a cylindrical pier.     

三维二阶水波绕射问题的有限元时域计算

本文用有限元法配合时步处理来求解三维非线性水波的绕射问题,自由表面条件和物面条件都满足到二阶,采用人工阻尼区来吸收反射波,流场内的速度势通过求解有限元方程得到。对垂直圆柱体的绕射问题进行了计算,得到了自由表面波高时间历程和圆柱所受到的波浪力,计算结果和有关文献的理论计算结果进行了比较。     

Benchmark computations of wave run-up on single cylinder and four cylinders by naoe-FOAM-SJTU solver

The benchmark simulations of wave run-up on a fixed single truncated circular cylinder and four circular cylinders are presented in this paper. Our in-house CFD solver naoe-FOAM-SJTU is adopted which is an unsteady two-phase CFD code based on the open source package OpenFOAM. The Navier-Stokes equations are employed as the governing equations, and the volume of fluid (VOF) method is applied for capturing the free surface. Monochromatic incident waves with the specified wave period and wave height are simulated and wave run-up heights around the cylinder are computed and recorded with numerical virtual wave probes. The relationship between the wave run-up heights and the incident wave parameters are analyzed. The numerical results indicate that the presented naoe-FOAM-SJTU solver can provide accurate predictions for the wave run-up on one fixed cylinder and four cylinders, which has been proved by the comparison of simulated results with experimental data.     

马蹄波波形特征和对圆柱作用实验研究

为了研究真实海洋表面马蹄波特性以及对建筑物的作用,通过物理模型实验研究了马蹄波的波形特征参数以及马蹄波对圆柱体的作用。实验中通过对浪高仪采集的波面升高时间历程曲线进行分析得出了不同水深情况下马蹄波的垂向几何特征,并通过快速傅里叶变换得出了马蹄波波幅谱的特征,研究了马蹄波各组成波波幅沿空间的变化从而得出了圆柱存在对马蹄波演化的影响,同时给出马蹄波绕射形成的波面分布和不同频率谐波在圆柱周围的分布,讨论了马蹄波不同于Stokes波对圆柱作用的特征。结果表明,马蹄波波形受水深影响较大,水深越浅,马蹄波的波面形状越接近椭圆余弦波。圆柱体的存在干扰了马蹄波不稳定的增长,使其在接近圆柱时呈下降趋势,导致不稳定幅值最大值的位置提前并且出现在偏离圆柱迎浪点的侧表面,从而使圆柱受到侧向力的作用。     

Second-Order Wave Diffraction Around 3-D Bodies by A Time-Domain Method

A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.     

水下平板与波浪相互作用完全非线性数值模拟(英文)

利用完全非线性数值波浪水槽技术研究水下平板与波浪的相互作用。假定水下平板厚度极薄、刚性,位于有限水深并且非常接近自由水面。应用四阶龙格库塔方法追踪每一时刻的波面形状,采用阻尼层来吸收反射波以保证算法的稳定性,同时引入平滑和重组的方法抑制自由表面控制点的较高梯度。通过对波浪与浮动圆柱相互作用的数值模拟证实了数值波浪水槽方法的有效性,计算结果与线性理论吻合良好。在波浪数值水槽方法中引入造波板模拟波浪产生并与水下平板发生相互作用,应用傅立叶解析方法对波面变形、波浪力作了分析。结果表明在板非常接近自由水面的情况下会表现出现很强的非线性,揭示了线性理论的局限性。     

Numerical simulation of solitary wave scattering by a circular cylinder array

The interaction of a solitary wave with an array of surface-piercing vertical circular cylinders is investigated numerically. The wave motion is modeled by a set of generalized Boussinesq equations. The governing equations are discretized using a finite element method. The numerical model is validated against the experimental data of solitary wave reflection from a vertical wall and solitary wave scattering by a vertical circular cylinder respectively. The predicted wave surface elevation and the wave forces on the cylinder agree well with the experimental data. The numerical model is then employed to study solitary wave scattering by arrays of two circular cylinders and four circular cylinders respectively. The effect of wave direction on the wave forces and the wave runup on the cylinders is quantified.     

Diffraction of Water Waves by A Vertically Floating Cylinder in A Two-Layer Fluid

In this paper, the diffraction of water waves by a vertically floating cylinder in a two-layer fluid of a finite depth is studied. Analytical expressions for the hydrodynamic loads on the vertically floating cylinder are obtained by use of the method of eigenfunction expansions. The hydrodynamic loads on the vertically floating cylinder in a two-layer fluid inelude not only the surge, heave and pitch exciting forces due to the incident wave of the surface-wave mode, but also those due to the incident wave of the internal-wave mode. This is different from the case of a homogenous fluid. Some given examples show that, for a two-layer fluid system with a small density difference, the hydrodynamic loads for the surface-wave mode do not differ significantly from those due to surface waves in a single-layer fluid, but the hydrodynamic loads for the internal-wave mode are important over a wide range of frequencies. Moreover, also considered are the free surface and interface elevations generated by the diffraction wave due to the incident wave of the surface-wave and interhal-wave modes, and transfer of energy between modes.     

大尺度海工圆柱的非线性波浪载荷

引入Stokes二阶水波理论,建立了对海洋工程中大尺度综合型圆柱式结构波浪载荷的二阶修正统一积分算式。通过对波势特征解法的推广及引入干扰波叠加法、最小二乘法,使原仅适合于轴对称单柱的二阶散射波浪载荷快速积分算法可完全适用于类型广泛的任意型实用圆柱的计算。几则实例(圆台,单、双柱)的计算检验了该方法的有效性。     

透水cosine-type同心圆柱波浪绕射研究

在势能流及微小振幅波理论假设下,采用复合边界元素法(CBEM)数值解析规则波通过单根外壁透水cosine-type型同心圆柱结构物绕射。将cosine-type外壁不透水条件改为透水外,并在内部设有一半径为b的不透水内圆柱共同组成同心圆柱,为了验证CBEM数值模式正确性与可行性,将其退化成圆柱,均得到合理结果。数值计算条件包含波浪正向入射凸端、透水参数及绕射参数,探讨波浪通过单根cosine-type同心圆柱体四周水面波动变化情况,并与双重圆筒同心圆柱作比较。计算结果显示结构物外壁采用透水型式,可以大幅降低cosine-type同心圆柱体四周整个波浪绕射场的水面波动。     

圆筒形水工建筑物波浪荷工的试验研究

在圆筒形水工建筑物的波压力测定试验基础上 ,通过多种波压力理论计算方法进行计算、比较和分析 ,采用半经验半理论方法 ,探讨性地提出一种适合圆筒形水工建筑物波浪荷载的分布及计算方法 ,该法表现形式简单 ,计算结果接近实测值 ,可为相关的工程设计提供参考依据     

基于粘性流模型的筒型基础防波堤波浪力数值分析

筒型基础防波堤是一种新型港口海岸工程结构,其基础上部是由连续排列的圆筒构成的直立防浪墙.采用粘性流数值模型,研究连续圆筒防波堤上波浪力竖向分布、水平(沿圆筒环向)分布和波浪力合力特性,并对粘性流数值模型计算的平面直墙波浪力与海港水文规范方法计算结果;粘性流数值模型计算的连续圆筒墙面波浪力与平面直墙波浪力;无限长连续圆筒墙面波浪力与有限长连续圆筒墙面波浪力进行比较分析.针对所选工程算例,建议按<海港水文规范>中平面直墙波浪力计算方法确定连续圆筒防波堤上的波浪力时,波峰时考虑0.90左右的折减系数,波谷时考虑0.95左右的折减系数.     

Short-crested waves interaction with a concentric porous cylinder system with partially porous outer cylinder

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.