共查询到20条相似文献，搜索用时 31 毫秒
1.
A numerical study of vortex-induced rotations(VIRs) of an equivalent triangular cylinder, which is free to rotate in the azimuthal direction in a uniform flow, is presented. Based on an immersed boundary method, the numerical model is established, and is verified through the benchmark problem of flow past a freely rotating rectangular body.The computation is performed for a fixed reduced mass of m~*=2.0 and the structural stiffness and damping ratio are set to zero. The effects of Reynolds number(Re=25-180) on the characteristics of VIR are studied. It is found that the dynamic response of the triangular cylinder exhibits four distinct modes with increasing Re: a rest position,periodic rotational oscillation, random rotation and autorotation. For the rotational oscillation mode, the cylinder undergoes a periodic vibration around an equilibrium position with one side facing the incoming flow. Since the rotation effect, the outset of vortex shedding from cylinder shifts to a much lower Reynolds number. Further increase in Re leads to 2 P and P+S vortex shedding modes besides the typical 2 S pattern. Our simulation results also elucidate that the free rotation significantly changes the drag and lift forces. Inspired by these facts, the effect of free rotation on flow-induced vibration of a triangular cylinder in the in-line and transverse directions is investigated. The results show that when the translational vibration is coupled with rotation, the triangular cylinder presents a galloping response instead of vortex-induced vibration(VIV). 相似文献
2.
When fluid flow passes a cylinder, the drag crisis phenomenon occurs between the sub-critical and the super-critical Reynolds numbers. The focus of the present studies was on the numerical prediction of the drag crisis based on CFD methods. In this work, block structured meshes with refined grids near the cylinder surface and in the downstream were employed. Both 2D and 3D simulations were performed using various turbulence models, including the SST k − ω model, the k − ϵ model, the SST with LCTM, the DES model, and the LES model. In the convergence studies, the effects of the grid size, the time step, the first grid size and the aspect ratio (for 3D simulations) on the solutions were examined. The errors due to spatial and time discretizations were quantified according to a V&V procedure. Validation studies were carried out for various Reynolds numbers between Re = 6.31 × 104 and 7.57 × 105. The averaged drag force, the RMS of lift force and the Strouhal number were compared with experimental data. The studies indicated that standard 2D and 3D RANS methods were inadequate to capture the drag crisis phenomenon. The LES method however has the potential to address the problem. 相似文献
3.
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. 相似文献
4.
两层粘性流体中圆柱体受迫振荡数值模拟 总被引：1，自引：1，他引：1
研究两层粘性流体中无限长水平圆柱体的受迫振荡问题。在湍流模式下,采用VOF方法追踪两层流体的内界面,基于动网格技术模拟圆柱体的运动边界,对均匀流中横向振荡圆柱体的绕流场进行了数值模拟。计算受迫振荡圆柱体的升力系数、阻力系数随时间的演化曲线和圆柱体的尾涡分布,以及圆柱体的受迫振荡激发两层流体内界面的扰动,并与均匀流体的情况进行了比较分析。研究表明,流体的两层分层效应对受迫振荡圆柱体的升阻力系数和尾涡分布特性都有显著影响,在水下输油气管道涡激振动特性的工程评估中,应考虑流体的密度分层效应。 相似文献
5.
利用计算流体力学软件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%。 相似文献
6.
《Coastal Engineering》2001,42(2):163-172
Applying the Green's theorem and a moving oscillating source as the Green's function, an integral equation method is developed for predicting wave height over a shoal in a weak current. The integral equation is discretized by a higher order element method and a numerical code is implemented. To validate the numerical code, comparison is made on wave run-up around a vertical cylinder with McCamy and Fuchs' analytic solution in the still water and the results in currents from numerical wave tanks. Computation has also been carried out for the wave diffraction over a parabolic shoal in a current, and numerical results are given at many sections. The numerical results from the calculation can be used for benchmark test for other methods on wave diffraction in current. 相似文献
7.
Data on the acceleration of an oscillating cylinder in water are analysed by a Kalman extended filter. A hydrodynamical forcing model based upon Morison's equation appears to be reasonable and drag estimation appears simple using simulated data. Estimation based upon real data gives an average drag coefficient decreasing with the oscillation amplitude. An assymetric drag coefficient variation over the oscillation period is estimated. 相似文献
8.
When a high-speed body with cavity passes through water-air free surface and exits water, its mechanical environment and dynamic characteristics change significantly due to the great difference in density and viscosity between water and air. With focusing on this problem, the Computational Fluid Dynamics (CFD) method is applied to perform numerical calculation on the process of this vapor-liquid-gas flow during the water exit of a high-speed cylinder, with the Volume of Fraction (VOF) multiphase flow interface-capturing techniques and the overset grid technology. After the verification and validation of the CFD model through mesh convergence study and a water-entry experiment, cavity evolution and flow characteristics including pressure and velocity distribution during the water exit are analyzed. The effects of different initial velocities on the pressure distribution and drag characteristics of the cylinder are investigated. Calculated results show that the cavity collapse during water exit causes strong pressure fluctuation on the cylinder; when the cylinder exits water enveloped in a supercavity, the pressure distribution on its wall surface and surrounding water region is relatively uniform, and the drag changes gently, and thus the cylinder has good motion stability.
相似文献9.
基于光滑粒子流体动力学(SPH)方法,开发了能够准确描述水流作用下圆柱强迫振动特性的数学模型。通过引入适合于无网格粒子法的开边界算法,来模拟出入流边界条件,建立了具有造流功能的SPH数值水槽。圆柱及计算域的上下边界均采用修正的动力边界条件进行模拟。借助于粒子位移矫正和压力修正算法,避免了圆柱周围流体粒子压力大幅震荡以及结构下游区域出现空腔等非物理性现象。使用典型的圆柱绕流数据,验证了所建SPH模型的计算性能,研究了固定圆柱在低雷诺数情况下的尾涡脱落模式和升阻力变化规律。明确了低雷诺数下强迫振动圆柱在频率锁定以及非锁定区间内的升力变化规律,量化了升力与外界激励频率之间的关系。 相似文献
10.
尝试建立1种基于无单元法的涡激振动数值模拟算法,并给出适用于涡激振动分析的前处理自动布点方法.基于无单元法理论,使用动最小二乘法构造形函数,利用无单元伽辽金法,采用速度和压力分离模式,及手工布点和自动布点2种方法所得出的节点布置形式,对流场控制方程进行空间离散,模拟二维涡激振动的流场形态,并对VIV相关参数进行分析,计算不同节点布置情况下的升力系数(Cl)、曳力系数(Cd)及斯特罗哈数(St),并与物理模型实验结果进行对比.计算结果表明,无单元伽辽金法应用于立管VIV分析是可行的,且文中采用2种布点方法均能较好的模拟流场中泻涡脱落的形态,计算结果与传统方法和物理模型实验结果吻合良好.证明文中2种布点方法都能用于固定圆柱的二维VIV分析,但自动布点法能够更好地适用于复杂问题的计算及圆柱体在流场中的VIV动力响应分析. 相似文献
11.
12.
A parametric study was carried out to investigate the hydrodynamics of a cylindrical wave energy absorber. Established methods of hydrodynamic analysis were applied to the case of a damped vertically oriented cylinder pivoted near the sea floor in intermediate depth water. The simple geometry provides a canonical reference for more complex structure shapes and configurations that may be considered for either wave energy conversion or wave energy absorption. The study makes use of the relative velocity Morison equation, with force coefficients derived from radiation and diffraction theory. Viscous effects were accounted for by including a drag term with an empirically derived coefficient, CD. A non-linear first-order formulation was used to calculate the cylinder motion response in regular waves. It was found that the non-linear drag term, which is often neglected in studies on wave energy conversion, has a large effect on performance. Results from the study suggest a set of design criteria based on Keulegan–Carpenter (KC) number, ratio of cylinder radius to water depth (a/h), and ratio of water depth to wavelength (h/L). Respectively, these parameters account for viscous, wave radiation, and water depth effects, and optimal ranges are provided. 相似文献
13.
Drag and inertia coefficients for a circular cylinder in steady plus low-amplitude oscillatory flows
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. 相似文献
14.
We present a numerical study on the hydrodynamic performance of undulation NACA0012 foil in the near wake of D-section cylinder. Computations are conducted using unsteady incompressible Navier-Stokes equations with a moving adaptive mesh based on laminar flow. Investigations are focused on the effect of distance ratio between foil tip and centre of cylinder (L/D≤2.0) on the thrust/drag performance of foil and cylinder at various foil undulation frequency (St). We found that, foil thrust coefficient (Ct) increases considerably with the appearance of cylinder and an optimal distance exists at which Ct reaches maxima. The maximum increment is about eleven times that of its counterpart of single foil, which is obtained at St=0.23 and L/D=0.5. Our results for the cylinder drag coefficient (Cd) observed the existence of optimal parametric map, combined with various gap ratios and foil frequencies. With these parameters, insertion of an undulation foil can significantly lead to the drag reduction indicating that undulating foil could work efficiently as a passive vortex control device for cylinder drag reduction. 相似文献
15.
W. Koterayama 《Ocean Engineering》1980,7(3):399-412
This paper presents a method of estimating wave forces acting on a submerged horizontal circular cylinder fixed in oblique waves.The experiments show that drag and inertia coefficients in beam sea are available for calculating the wave forces in oblique waves.Wave forces exerted on a vertical circular cylinder in deep waves are also investigated.The experimental results show that wave forces acting on the vertical cylinder coincide approximately with hydrodynamic forces acting on a submerged circular cylinder in an oscillating fluid. 相似文献
16.
A parameterization of subgridscale surface fluxes over the marginal sea ice zone which has been used earlier in several studies is modified and applied to a nonhydrostatic mesoscale model. The new scheme accounts for the form drag of ice floes and is combined with a so-called flux averaging method for the determination of surface fluxes over inhomogeneous terrain. Individual fluxes over ice and water are calculated as a function of the blending height. It is shown by comparison with observations that the drag coefficients calculated with the new parameterization agree well with data. The original scheme strongly overestimates the form drag effect. An improvement is mainly obtained by an inclusion of stratification and by use of a more adequate coefficient of resistance for individual ice floes. The mesoscale model is applied to off-ice flows over the polar marginal sea ice zone. The model results show that under certain meteorological conditions the form drag can have a strong influence on the near-surface wind velocity and on the turbulent fluxes of momentum. Four case studies are carried out. The maximum influence of form drag occurs in the case with strong unstable stratification and with wind oblique to the ice edge. Under these conditions the wind stress on sea ice is modified by at least 100% for ice concentrations less than 50% if form drag is taken into account. 相似文献
17.
The incompressible viscous uniform and shear flow past a circular cylinder is studied. The two-dimensional Navier-Stokes equations are solved by a finite element method. The governing equations are discretized by a weighted residual method in space. The stable three-step scheme is applied to the momentum equations in the time integration. The numerical model is firstly applied to the computation of the lid-driven cavity flow for its validation. The computed results agree well with the measured data and other numerical results. Then, it is used to simulate the viscous uniform and shear flow over a circular cylinder for Reynolds numbers from lO0 to lO00. The transient time interval before the vortex shedding occurs is shortened considerably by introduction of artificial perturbation. The computed Strouhal number, drag and lift coefficients agree well with the experimental data. The computation shows that the finite element model can be successfully applied to the viscous flow problem. 相似文献
18.
This paper presents the results of an experimental investigation on the variation in the tension and the distribution of drag force coefficients along flexible risers under vortex-induced vibration (VIV) in a uniform flow for Reynolds numbers (Re) up to 2.2×105. The results show that the mean tension is proportional to the square of the incoming current speed, and the tension coefficient of a flexible riser undergoing VIV can be up to 12. The mean drag force is uniformly and symmetrically distributed along the axes of the risers undergoing VIV. The corresponding drag coefficient can vary between 1.6 and 2.4 but is not a constant value of 1.2, as it is for a fixed cylinder in the absence of VIV. These experimental results are used to develop a new empirical prediction model to estimate the drag force coefficient for flexible risers undergoing VIV for Reynolds number on the order of 105, which accounts for the effects of the incoming current speed, the VIV dominant modal number and the frequency. 相似文献
19.
Measurements are presented of the effects of corner radius on the hydrodynamic forces experienced by cylindrical bluff bodies in oscillatory flow over the Keulegan Carpenter number range from 1 to 100. The range of Reynolds number covered is from 200 to 2 × 104 and two basic cross-sectional shapes, namely a square and a diamond, were used in the investigation. Drag coefficients in oscillatory flow are found to be more sensitive to corner radius than in steady flow. At a fixed Keulegan Carpenter number the drag coefficient of a diamond section decreases with increasing corner radius and the reduction is more marked at low Keulegan Carpenter numbers. The drag of a square section does not follow a consistent trend with corner radius and for high Keulegan Carpenter numbers the rounded square with a corner radius of 0.265D (where D is the width of the model normal to the flow) can have a drag coefficient equal to half the value for a circular cylinder. 相似文献
20.
The phenomenon of drag reduction by the injection of micro-bubbles into turbulent boundary layer has been investigated using an Eulerian-Eulerian two-fluid model. Multiple-size group (MUSIG) based on population balance models, which resolve a wide range of bubble sizes taking into account the bubble break-up and coalescence have been used for this purpose. The simulated results are compared against the experimental findings of Madavan et al. [1984. Reduction of turbulent skin friction by micro-bubbles. Physics of Fluids 27, 356-363] and also other numerical studies explaining the sophisticated phenomena of drag reduction. For the two Reynolds number cases considered, the buoyancy with the plate on the bottom configuration is investigated, as from the experiments it is seen that buoyancy seem to play a role in the drag reduction. Numerical model employed in the investigation comprises of a micro-bubble laden flow wherein two independent sets of Reynolds averaged Navier-Stokes (RANS) transport equations were used to describe both the phases of the flow. The shear stress transport (SST) turbulence model is used as the turbulent closure for the primary phase and a zero equation turbulence model is used for the micro-bubbles. Change in the mean streamwise velocity profiles, void fraction, turbulence modification and other results are presented and discussed with corresponding change in the gas injection rates. The complex mechanism of drag reduction are scrutinised and explained in context to our numerical findings. Special attentions have been also devoted to divulge the effect of bubble coalescence and break-up caused by random collision and turbulent impact. Numerical results showed good agreement for the skin-friction coefficients against experimental data throughout various air injection rates. The MUSIG model was found to be one of the best candidates to resolve the bubble dynamics in micro-bubble-induced drag reduction problems. 相似文献