Numerical Investigation of Vortex-Induced Vibration of a Circular Cylinder Close to a Plane Boundary Subject to Oscillatory Flow

2016 ◽  
Author(s):  
Adnan Munir ◽  
Ming Zhao ◽  
Helen Wu

This paper presents a numerical study of flow around an elastically mounted circular cylinder in close proximity to a plane boundary vibrating in the transverse and inline directions in an oscillatory flow. The Reynolds-Averaged Navier-Stokes (RANS) equations and the SST k-ω turbulent equations are solved using the Arbitrary Langrangian-Eulerian (ALE) scheme and Petrov-Galerkin Finite Element Method for simulating the flow. The equation of motion is solved using the fourth-order Runge-Kutta method to find the displacements of the cylinder in the transverse and inline directions. The numerical model is validated against the previous results of vortex-induced vibration of an isolated circular cylinder in both cross-flow and inline directions. The flow model is further extended to study the vortex-induced vibration of a cylinder near a plane boundary with a very small gap ratio (e/D) of 0.01, with D and e being the diameter and the gap between the cylinder and the plane boundary, respectively. Simulations are carried out for two Keulegan-Carpenter (KC) numbers of 5 and 10 and a wide range of reduced velocities. It is observed that both the KC number and the reduced velocity affect the vibration of the cylinder significantly.

Author(s):  
Ming Zhao ◽  
Liang Cheng ◽  
Tongming Zhou

Vortex-induced vibration (VIV) of a circular cylinder in oscillatory flow is investigated numerically in this study. The incompressible Reynolds-Averaged Navier-Stokes equations governing fluid flow around a circular cylinder are solved using Arbitrary Langrangian-Eulerian (ALE) scheme and Petrov-Galerkin finite element method. The equation of motion is solved for the displacements of the cylinder both in the inline and cross-flow directions. The numerical model is firstly validated against the experimental results of one-degree-of-freedom VIV in cross-flow direction. It is found that both VIV frequency and amplitude vary with reduced velocity for a fixed KC number. In most of the simulated cases the vibration comprises of multiple frequencies of different amplitudes. Each frequency component is multiple times of the frequency of the oscillatory flow. Two-degree-of-freedom VIV is investigated with the same parameters used in the one-degree-of-freedom case. By examining the XY-trajectory of the vibration, it if found that the vibration follows different trajectory for different KC numbers or reduced velocities.


Author(s):  
Bruno C. Ferreira ◽  
Marcelo A. Vitola ◽  
Juan B. V. Wanderley ◽  
Sergio H. Sphaier

The vortex-induced vibration (VIV) is a classical problem in ocean engineering. Intensive research on this field for flow around a circular cylinder has been observed, due to practical application, mainly the design of risers, cables and pipelines with free span. The relevance of this phenomenon is related to the structure failure, consequence of large displacement or fatigue. In the present study the influence of initial condition on the vortex induced vibration (VIV) of a circular cylinder with two degree of freedom is investigated by the numerical solution of the slightly compressible formulation of Reynolds Average Navier-Stokes equations. An upwind and Total Variation Diminishing (TVD) conservative scheme is used to solve the governing equations written in curvilinear coordinates. The k–ε turbulence model is used to simulate the turbulent flow in the wake of the cylinder. Two different initial conditions have been tested, free-stream and continuous reduced velocity increase (using the previous reduced velocity as initial condition for the next value). Results for the phase angle, amplitude, frequency, and lift coefficient are presented. The numerical results have been compared with experimental data of Jauvtis and Williamson [1]. The results indicate that the history of cylinder movement has a important impact in the amplitude oscillation observed in-line and cross-flow, principally in the reduced velocity range associated with the upper branch. Results obtained for the initial and lower branch seems to be independent of the initial condition. Further investigation are necessary to understand the difference observed such as the absence of the jump in the cross-flow oscillation between the initial and upper branch and the absence of in-line oscillation for reduced velocity in the range of 1–4 and the peak of in-line oscillation at reduced velocity 6.0.


2021 ◽  
Author(s):  
Muhannad W. Gamaleldin ◽  
Alexander V. Babanin

Abstract A numerical study is conducted to investigate a sinusoidal oscillatory cross-flow around a smooth circular cylinder. A numerical model based on the Reynolds-averaged Navier-Stokes (RANS) equations with the k-ω Shear Stress Transport (SST) turbulence model is implemented, using the open source C++ package OpenFOAM®. The study focuses on the level of confidence in using two-dimensional numerical simulations for the predictions of the inherent transient behaviour of the highly nonlinear hydrodynamic induced loads. Simulations are conducted for two values of Keulegan-Carpenter (KC) number corresponding to inertia-dominated and drag-inertia regimes, for a fixed value of the frequency parameter β. Boundary layer separation structures and their temporal evolution are investigated and visualized by means of turbulence kinetic energy contours. Moreover, inline forces inertia and drag coefficients of Morison loads decomposition are found and validated against previous experimental results in literature. Finally, transverse forces and their spectral content is investigated and validated. The study shows that the magnitude and frequency of the unsteady hydrodynamic induced forces coefficients are well reproduced and are in a very good agreement with the experimental observations.


2018 ◽  
Author(s):  
Jiajun Chen ◽  
Yue Sun ◽  
Hang Zhang ◽  
Dakui Feng ◽  
Zhiguo Zhang

Mixing in pipe junctions can play an important role in exciting force and distribution of flow in pipe network. This paper investigated the cross pipe junction and proposed an improved plan, Y-shaped pipe junction. The numerical study of a three-dimensional pipe junction was performed for calculation and improved understanding of flow feature in pipe. The filtered Navier–Stokes equations were used to perform the large-eddy simulation of the unsteady incompressible flow in pipe. From the analysis of these results, it clearly appears that the vortex strength and velocity non-uniformity of centerline, can be reduced by Y-shaped junction. The Y-shaped junction not only has better flow characteristic, but also reduces head loss and exciting force. The results of the three-dimensional improvement analysis of junction can be used in the design of pipe network for industry.


2008 ◽  
Vol 596 ◽  
pp. 49-72 ◽  
Author(s):  
HIROSHI HIGUCHI ◽  
HIDEO SAWADA ◽  
HIROYUKI KATO

The flow over cylinders of varying fineness ratio (length to diameter) aligned with the free stream was examined using a magnetic suspension and balance system in order to avoid model support interference. The drag coefficient variation of a right circular cylinder was obtained for a wide range of fineness ratios. Particle image velocimetry (PIV) was used to examine the flow field, particularly the behaviour of the leading-edge separation shear layer and its effect on the wake. Reynolds numbers based on the cylinder diameter ranged from 5×104 to 1.1×105, while the major portion of the experiment was conducted at ReD=1.0×105. For moderately large fineness ratio, the shear layer reattaches with subsequent growth of the boundary layer, whereas over shorter cylinders, the shear layer remains detached. Differences in the wake recirculation region and the immediate wake patterns are clarified in terms of both the mean velocity and turbulent flow fields, including longitudinal vortical structures in the cross-flow plane of the wake. The minimum drag corresponded to the fineness ratio for which the separated shear layer reattached at the trailing edge of the cylinder. The base pressure was obtained with a telemetry technique. Pressure fields and aerodynamic force fluctuations are also discussed.


1998 ◽  
Vol 120 (1) ◽  
pp. 72-75 ◽  
Author(s):  
V. N. Kurdyumov ◽  
E. Ferna´ndez

A correlation formula, Nu = W0(Re)Pr1/3 + W1(Re), that is valid in a wide range of Reynolds and Prandtl numbers has been developed based on the asymptotic expansion for Pr → ∞ for the forced heat convection from a circular cylinder. For large Prandtl numbers, the boundary layer theory for the energy equation is applied and compared with the numerical solutions of the full Navier Stokes equations for the flow field and energy equation. It is shown that the two-terms asymptotic approximation can be used to calculate the Nusselt number even for Prandtl numbers of order unity to a high degree of accuracy. The formulas for coefficients W0 and W1, are provided.


Author(s):  
Jean-Franc¸ois Sigrist ◽  
Cyrille Allery ◽  
Claudine Beghein

The present paper is the sequel of a previously published study which is concerned with the numerical simulation of vortex-induced-vibration (VIV) on an elastically supported rigid circular cylinder in a fluid cross-flow (A. Placzek, J.F. Sigrist, A. Hamdouni; Numerical Simulation of Vortex Shedding Past a Circular Cylinder at Low Reynolds Number with Finite Volume Technique. Part I: Forced Oscillations, Part II: Flow Induced Vibrations; Pressure Vessel and Piping, San Antonio, 22–26 July 2007). Such a problem has been thoroughly studied over the past years, both from the experimental and numerical points of view, because of its theoretical and practical interest in the understanding on flow-induced vibration problems. In this context, the present paper aims at exposing a numerical study based on a fully coupled fluid-structure simulation. The numerical technique is based on a finite volume discretisation of the fluid flow equations together with i) a re-meshing algorithm to account for the cylinder motion ii) a projection subroutine to compute the forces induced by the fluid on the cylinder and iii) a coupling procedure to describe the energy exchanges between the fluid flow and solid motion. The study is restricted to moderate Reynolds numbers (Re∼2.000–10.000) and is performed with an industrial CFD code. Numerical results are compared with existing literature on the subject, both in terms of cylinder amplitude motion and fluid vortex shedding modes. Ongoing numerical studies with different numerical techniques, such as ROM (Reduced Order Models)-based methods, will complete the approach and will be published in next PVP conference. These numerical simulations are proposed for code validation purposes prior to industrial applications in tube bundle configuration.


Sign in / Sign up

Export Citation Format

Share Document