首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
赵宇蒙  温鸿杰  任冰  王超 《海洋工程》2021,39(4):134-143
基于光滑粒子流体动力学(SPH)方法,开发了能够准确描述水流作用下圆柱强迫振动特性的数学模型。通过引入适合于无网格粒子法的开边界算法,来模拟出入流边界条件,建立了具有造流功能的SPH数值水槽。圆柱及计算域的上下边界均采用修正的动力边界条件进行模拟。借助于粒子位移矫正和压力修正算法,避免了圆柱周围流体粒子压力大幅震荡以及结构下游区域出现空腔等非物理性现象。使用典型的圆柱绕流数据,验证了所建SPH模型的计算性能,研究了固定圆柱在低雷诺数情况下的尾涡脱落模式和升阻力变化规律。明确了低雷诺数下强迫振动圆柱在频率锁定以及非锁定区间内的升力变化规律,量化了升力与外界激励频率之间的关系。  相似文献   

2.
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.  相似文献   

3.
In this paper, theoretical models are developed and numerical methods are used to analyze the loads, motions and cavity dynamics for freefall wedges with different deadrise angles vertically entering the water surface at Froude numbers: 1  Fn < 9. The time evolutions of the penetration depth, the velocity and the acceleration are analyzed and expressed explicitly. The maximum and average accelerations are predicted. The theoretical results are compared with numerical data obtained through a single-fluid BEM model with globally satisfactory agreement. The evolution of the pressures on the impact side is investigated. Before flow separation, gravity and the acceleration of the wedge have negligible influence on the pressure on the impact side for large Froude numbers or small deadrise angles; with increasing the deadrise angle or decreasing Froude number, the effects of gravity and the acceleration of the wedge tend to become more important. Global loads, with the main emphasis on the drag coefficient, are also studied. It is found that for the light wedge, the transient drag coefficient has slow variation in the first half of the collapse stage and rapid variation in the last half of the collapse stage. For the heavy wedge, the transient drag coefficients vary slowly during the whole collapse stage and can be treated as constant. The characteristics of the transient cavity during its formation are investigated. The non-dimensional pinch-off time, pinch-off depth and submergence depth at pinch-off scale roughly linearly as the Froude number.  相似文献   

4.
The problem of oscillating inviscid flow over an elliptic cylinder is studied for various angles of attack. The flow is incompressible and two-dimensional and the oscillations are harmonic. The flow direction is always normal to the cylinder axis and oscillations are only allowed in the magnitude of the free stream velocity. The study focuses on the hydrodynamic forces acting on the cylinder as well as the surface pressure distribution and their time variation. The parameters involved are the cylinder axis ratio, the angle of attack and the Strouhal number. The variables are normalized in such a way that the solutions for the problems of oscillating flow over flat plates and circular cylinders can be easily obtained as special cases. Analytical expressions are given for the drag coefficient, the lift coefficient, and the surface pressure distribution and their variation with time. Similar analytical expressions are given for the special cases of circular cylinders and inclined flat plates in oscillating inviscid flows.  相似文献   

5.
The successful designs of hulls for ships employing drag reduction by air bottom cavitation have been based on solutions of inverse problems of the theory of ideal incompressible fluid. However, prediction of the drag reduction ratio, the air demand by ventilated cavities and the cavity impact on the hull–propeller interaction is impossible in the framework of this theory because all mentioned characteristics depend on interaction of air cavities with the ship boundary layers. Because the known CFD tools are not fitted to ventilated cavitation at low Froude numbers, an analysis of this interaction requires a novel flow model. This model includes the incompressible air flow in the ventilated cavity, the compressible flow of a water–air mixture in the boundary layer on cavities and downstream of them and the curl-free incompressible outer water flow. The provided 2D computations employing this model allows for explanations of the earlier observed effects and for prediction of the air demand by ventilated cavities. The computed velocity profiles downstream of cavities are in the accordance with the available experimental data.  相似文献   

6.
This paper presents a numerical study on the drag reduction mechanism created by a ventilated partial cavity and its associated effects by the downstream dispersed microbubbles. A semi-empirical approach is introduced to model the discrete interface of the ventilated cavity and its complex gas leakage rate induced by the local turbulent shear stress. Based on the Eulerian–Eulerian two-fluid modeling framework, a population balance approach based on MUltiple-SIze-Group (MUSIG) model is incorporated to simulate the dynamical effects of bubbly flow along the test body. Particular attention is also directed to grasp a better understanding of the size evolution of microbubble and its associated effects on drag reduction. Model predictions are validated against three experimental measurements carried out in a high-speed water tunnel by Schauer (2003) and Wosnik et al. (2005). Close examination of the flow structures, gas void fraction distributions and its resultant density ratio provides valuable insights on the complex physical phenomenon, helping to consolidate idea to maximize the drag reduction for ventilated cavitating vehicles.  相似文献   

7.
In the present study, flow characteristics were investigated experimentally using particle image velocimetry technique (PIV) in a gap between a solid cylinder and a shroud to reveal the effect of shroud diameter (Ds) and porosity (β) on the vortex shedding mechanism of the cylinder. Porosity (varied from β = 0.3 to 0.7) and diameter ratio (D/Ds = 0.4, 0.5 and 0.6) were main parameters examined at a Reynolds number of Re = 5000. For the porosity values of β ≤ 0.5, it is observed that vortex formation of the cylinder occurs only in the gap and shroud produces its own wake flow patterns. Penetrating flow through the shroud extends the shear layers on the both sides of the shroud through the downstream direction and prevents the interaction of shear layers in the near wake region. The diameter ratio and the porosity are impactful on the wake flow patterns in outer region of the shroud since they are determinant of the penetrating flow rate. Force measurements were also performed in the air tunnel in order to reveal the effect of shroud on the drag coefficient of cylinder. It is found that the drag coefficient of the cylinders are reduced significantly by shrouds when compared with that obtained from the bare cylinder case. However, the drag coefficient of the cylinder together with the shroud is higher than the bare cylinder for all cases since the shrouds enlarge the area exposed to the flow.  相似文献   

8.
Direct numerical simulation was conducted to investigate the flow past a slotted cylinder at low Reynolds number (Re) of 100. The slotting of cylinder affects the boundary layer separation, vortex formation position, recirculation region length and wake width, which are determined by the type of slit. The streamwise slit (SS1), T-shaped slit (SS3) and Y-shaped slit (SS4) act as passive jets, while the transverse slit (SS2) achieves an alternate self-organized boundary layer suction and blowing. The flow rate in slits fluctuates over time due to the alternate vortex shedding and fluctuating pressure distribution around the cylinder surface. One fluctuation cycle of flow rate is caused by a pair of vortices shedding for SS2, SS3 and SS4, while it is created by each vortex shedding for SS1. The wall shear stress and flow impact on the slit wall partly contribute to the hydrodynamic forces acting on the slotted cylinder. Taking into account the internal wall of slit, the transverse slit plays the best role in suppressing the fluid forces with drag reduction of 1.7% and lift reduction of 17%.  相似文献   

9.
The wave induced dynamic pressures around a circular cylinder of diameter 0.2 m due to regular waves were measured in a wave flume in a water depth of 1 m and in a wave basin in a water depth of 3 m. The experimental investigations were carried out with the cylinder inclined along and against the direction of wave propagation. The least-squares technique was employed to evaluate the coefficients of drag (CD) and inertia (CM) from the sectional force time histories obtained by integrating the measured circumferential pressure distribution. The variation of drag and inertia coefficients are presented as a function of Keulegan-Carpenter number (KC) for different inclinations of the cylinder. The comparison between the measured and the theoretical force derived from the evaluated hydrodynamic coefficients is found to be good.  相似文献   

10.
The main purpose of this study is to establish a better understanding of the relationship between drag reduction and surface roughness. Experiments were conducted to measure the force and flow characteristics of a circular cylinder with different types of artificial surface roughness over the range 6 × 103 < Re < 8 × 104 (Re is based on the cylinder diameter D). The roughness cylinder was formed by covering the exterior surface of the cylinder with uniformly distributed (1) sandpaper, (2) netting, and (3) dimples. The roughness coefficient ranged from k/D = 0.0028 to 0.025 (k is the roughness height). A detailed quantitative measurement of the flow field around the cylinder using Particle Imaging Velocimetry (PIV) was carried out. The hydrodynamic force coefficients (drag and lift) of the rough cylinders are compared against those of a smooth cylinder measured under the same flow conditions. It is found that certain configuration of surface roughness significantly reduces the mean drag coefficient of the cylinder, particularly at large Reynolds numbers. In addition, the root-mean-square (r.m.s.) lift coefficient of the rough cylinders is considerably lower than that of a smooth cylinder.  相似文献   

11.
利用计算流体力学软件Fluent开展了三维动波浪壁圆柱绕流的数值模拟,建立了三维运动波浪壁圆柱模型,通过C语言自编程序实现波浪壁面的运动控制,并保证壁面变形时网格的高质量。在来流速度u=0.125 m/s、雷诺数Re=12 500的情况下,开展了动波浪壁波动速度w=0、0.062 5、0.125、0.187 5 m/s四个工况的计算分析,并比较了不同波动速度对流场结构、升力、阻力特性的影响。结果表明:动波浪壁圆柱能有效抑制流动的分离,消除交替脱落的尾涡,从而消除周期振荡的升力;在消除卡门涡街的同时,圆柱后驻点处的涡量值随波动速度增加而增加,其原因在于波形移动加大了壁面流体的速度,从而减小了圆柱前后的压力差,减小了阻力;随着波动速度的增大,平均阻力系数呈明显下降趋势,当波动速度为来流速度的1.5倍时,平均阻力系数相对于光滑圆柱下降了53.76%。  相似文献   

12.
Numerical benchmark studies on drag and lift coefficients of a marine riser have been organized by the 27th ITTC Ocean Engineering Committee. The purpose of the studies was to benchmark the capabilities of CFD methods through quantitative comparisons and validation studies against the model test results of a circular cylinder by MARIN. Studies were focused on the drag crisis phenomenon for the stationary smooth cylinder in the critical Reynolds number regime. Eight organizations have participated in the studies by using RANS, DES and LES methods. An overview of the model test results, test cases, submissions and comparison results are presented in this paper. Conclusions and recommendations are made for future studies.  相似文献   

13.
A horizontal, circular cylinder fitted with one bilge keel is forced to rotate harmonically around its axis. The bilge keel load and hull pressure distribution are investigated. A fully submerged condition (infinite fluid), and three partly-submerged conditions are considered. A two-dimensional numerical study is performed, and the results are validated against recently published experimental data by van’t Veer et al. [30]. In addition, comparisons for mass and drag coefficients are also made with experimental data for plate in infinite fluid (Keulegan and Carpenter [8]), and wall-mounted plate (Sarpkaya and O’Keefe [9]) in oscillatory flow.A Navier–Stokes solver based on the Finite Volume Method is adopted for solving laminar flow of incompressible water. The free-surface condition is linearized by neglecting the nonlinear free-surface terms and the influence of viscous stresses in the free surface zone, while the body-boundary condition is exact. This simplified modeling of the problem required the mesh to be fine only around the bilge keels, leading to a total number of cells around N  1 ×104, which reduced computational cost significantly.The influence of draft and amplitude of oscillations on the bilge keel force and hull pressure distribution are considered. The bilge keel force is presented in terms of non-dimensional drag and mass coefficients including higher harmonic components. The numerical results are also compared with the industry standard empirical method for calculation of roll damping proposed by Ikeda et al. [4]. In general, a good agreement between the results of the present numerical method and the experimental data is obtained and the differences with those predicted by the empirical method are addressed.  相似文献   

14.
The flow around a cylinder, surface mounted or partly buried has been examined at Reynolds number of 1.3×104 and 2.6×104 by using the flow visualization method in the case of steady current. The lengths of the upstream and downstream separation regions were determined and the characteristics of the flow area were discussed for the cylinder having different burial depth to the diameter ratio (G/D). For comparsion, some numerical results obtained from the application of the FLUENT, computational fluid dynamics (CFD) software package have also been presented. Both experimental and numerical results indicated that the lengths of the separation regions near the upstream and downstream of the cylinder decreased with the increasing burial ratio (G/D=0, 0.20, 0.50). The results have significant implications both for the illustration of the potential for using CFD software and highlighting the need for data of two dimensional flow around the buried cylinder.  相似文献   

15.
Computer simulations of steady plus low-amplitude oscillatory flow about a circular cylinder are reported at a fixed Reynolds number of 150 based on the steady component. The conventional Keleugan–Carpenter number based on the oscillatory component is fixed at π/5. The oscillation frequency is varied so as to study a wide spectrum of flows where inertial forces dominate at one end and viscous drag forces at the other as a function of the modified Keleugan–Carpenter number. The hydrodynamic force on the cylinder in-line with the flow direction is represented by Morison's equation and an extended version with three terms. The drag and inertia coefficients in Morison's equation are determined by least-squares fits to data directly computed from integration of skin friction and pressure distributions around the periphery of the cylinder. The root-mean-square value of the residue of reconstructed minus directly-computed forces varies between 2 and 41% depending on the flow parameters. Comparable results can be obtained with a semi-theoretical approach using inviscid inertia and quasi-steady viscous drag terms. Physical explanations for the variation of the force coefficients are provided and implications for pertinent flow–structure interactions are discussed.  相似文献   

16.
The present study extends the investigations of the hydrodynamic forces on a cylinder, laid on, or partly buried in the bed. They were determined by measuring the pressure distribution on the cylinder surface in the case of steady current, waves and coexisting flow. The pressure distribution around the cylinder was measured by using pressure transducers, which were replaced in the cylinder. Force coefficients were obtained for the ranges of Re=0.8×104–1.5×104, for steady current, low KC numbers (KC<5) for wave alone case and, for current-to-wave velocity RATIO=0, 3, 6 and infinity (current) for coexisting flow. The forces were also determined for the various burial-depth-to-diameter ratios between 0 and 0.7 values of the cylinder.  相似文献   

17.
在亚临界区雷诺数下,采用脱体涡模拟方法对弱电解质中电磁力作用下圆柱绕流场及其升阻力特性进行了数值模拟与分析。结果表明,电磁力可以提高圆柱体边界层内的流体动能,延缓圆柱体近壁面流动分离,减弱绕流场中流向和展向大尺度涡的强度,减小圆柱体阻力及其升力脉动幅值;当电磁力作用参数大于某个临界值后,流动分离角消失,在圆柱体尾部产生射流现象,电磁力产生净推力作用,出现负阻力现象,而且升力脉动幅值显著减小且接近于零。  相似文献   

18.
A series of experiments is described in which energy losses are measured in a reversing oscillatory flow through a pipe exit. The experiments provide a comparative study of the losses associated with sinusoidal flow through exits of various shape, and in particular how the curvature of the exit determines the degree of flow loss due to the shedding of vortices. It is shown that geometrical factors relating the amplitude of motion to the curvature of the mouth and radius of the tube have an important bearing on the incidence of energy loss; significant losses are only observed to occur once a particular amplitude of motion is exceeded. In addition photographic visualization is used to provide evidence connecting the generation of eddies with increased losses.  相似文献   

19.
高洋洋  张演明  刘彩  王滨 《海洋工程》2020,38(1):86-100
基于计算流体力学(CFD)开源代码OpenFOAM开展了不同雷诺数(Re=100、1500和3900)和倾斜角度(-60°≤α≤60°)工况下倾斜圆柱绕流流场的三维数值模拟,研究倾斜圆柱绕流的三维瞬时及时均尾流流场、流线拓扑、升阻力系数与旋涡脱落频率随雷诺数及倾斜角度变化的规律,探讨在顺流向及逆流向情况下独立性原则对倾斜圆柱绕流的适用性。研究结果表明:随着圆柱倾角的增大,倾斜圆柱尾流产生较为明显的轴向流,尾流旋涡脱落受到明显抑制,细碎旋涡逐渐消失,尾流宽度随之减小;随着雷诺数的增大,圆柱尾流涡管发生明显的变形,展向掺混使得大量细碎旋涡产生,呈现出明显的三维特性。在不同雷诺数下,阻力系数均值、升力系数均方根及无量纲涡脱频率在一定倾角范围内符合独立性原则。  相似文献   

20.
Experimental studies of the transverse motions induced in a long flexible cable by axial flows have been conducted in a large blowdown water channel. A cable with a diameter of 1.59 cm and a length of 9 m was employed. Experiments were conducted over a fluid velocity range of 4.6-9.1 m/s. Both free and fixed downstream terminations of the cable were utilized, the former simulating a towed cable and the latter a mooring cable. Measurements of the drag coefficient of the cable with the free end were also made. These results were compared with those obtained previously for a flexible cylinder with a smooth exterior surface, and the following conclusions were drawn. 1) The amplitudes of the flow-induced transverse motions in a flexible cable or cylinder with a free downstream end in general decrease with increasing distance from this end, but they exhibit weak dependence on flow velocity and surface roughness near this end. 2) The amplitudes of the flow-induced motions are in general smaller when the cable is fixed at both ends than when the downstream end is free. 3) The drag coefficient of the flexible cylinder or cable increases with increasing surface roughness.  相似文献   

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

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