Numerical Simulation of Vortex Induced Vibration on an Elastically Mounted Cylinder

Author(s):  
Juan B. V. Wanderley ◽  
Carlos A. Levi

The Vortex-induced vibration on a circular cylinder is investigated by the numerical solution of the unsteady Reynolds Average Navier-Stokes equations and results are compared with experimental measurements obtained by different authors. The Beam and Warming implicit factored scheme is used to solve the governing equations and the Baldwin and Lomax model is used to simulate the turbulent flow in the wake of the cylinder. The cylinder is laterally supported by a spring and a damper and is free to oscillate in the transversal direction in an initially uniform flow for the first flow speed investigated. For the subsequent speeds, the final condition obtained for the previous speed is used as initial condition to reproduce the actual experimental set up. In that case, the measurements are done by progressive increments of the flow speed retaining the fluid memory effect. The complexity and high sensitivity of the flow phenomenon at this configuration requires a very accurate and robust numerical model. Most of the known algorithms failed to duplicate the available experimental measurements. The numerical results for the transversal oscillation amplitude are compared to experimental data showing a fairly precise agreement at the difficult to simulate regime of the lock-in phenomenon.

Author(s):  
Juan B. V. Wanderley ◽  
Carlos A. Levi

The Vortex-induced vibration on a circular cylinder is investigated by the numerical solution of the unsteady Navier-Stokes equations and results are compared with experimental measurements obtained by different authors. The Beam and Warming implicit factored scheme is used to solve the governing equations and Large Eddy Simulation is used together with the Smagorinsky subgrid-scale model (SGS) to simulate the turbulent flow in the wake of the cylinder. The cylinder is laterally supported by a spring and a damper and is free to oscillate in the transversal direction in an initially uniform flow for the first flow speed investigated. For the subsequent speeds, the final condition obtained for the previous speed is used as initial condition to reproduce the actual experimental set up. In that case, the measurements are done by progressive increments of the flow speed retaining the fluid memory effect. The complexity and high sensitivity of the flow phenomenon at this configuration requires a very accurate and robust numerical model. Most of the known algorithms failed to duplicate the available experimental measurements. The proposed numerical solution was able to provide a good picture of the real physics of the phenomenon showing the Ka´rma´n vortex street effects on the lift and drag coefficients. The numerical results for the transversal oscillation amplitude are compared to experimental data showing a fairly precise agreement at the difficult to simulate regime of the lock-in phenomenon.


1976 ◽  
Vol 78 (2) ◽  
pp. 355-383 ◽  
Author(s):  
H. Fasel

The stability of incompressible boundary-layer flows on a semi-infinite flat plate and the growth of disturbances in such flows are investigated by numerical integration of the complete Navier–;Stokes equations for laminar two-dimensional flows. Forced time-dependent disturbances are introduced into the flow field and the reaction of the flow to such disturbances is studied by directly solving the Navier–Stokes equations using a finite-difference method. An implicit finitedifference scheme was developed for the calculation of the extremely unsteady flow fields which arose from the forced time-dependent disturbances. The problem of the numerical stability of the method called for special attention in order to avoid possible distortions of the results caused by the interaction of unstable numerical oscillations with physically meaningful perturbations. A demonstration of the suitability of the numerical method for the investigation of stability and the initial growth of disturbances is presented for small periodic perturbations. For this particular case the numerical results can be compared with linear stability theory and experimental measurements. In this paper a number of numerical calculations for small periodic disturbances are discussed in detail. The results are generally in fairly close agreement with linear stability theory or experimental measurements.


Author(s):  
П.А. Поливанов

In this paper a numerical and experimental study of the effect of blowing/suction through a perforated surface on a turbulent boundary layer at a Mach number M = 1.4 is carried out. Most of the calculations were performed by Reynolds-averaged Navier-Stokes equations with the k-w SST turbulence model. The calculated geometry completely repeated the experimental one including the perforated surface. The numerical data were compared with experimental measurements obtained by the PIV method. Analysis of the data made it possible to find the limits of applicability of the numerical method for this flow.


Author(s):  
Xingwei Zhang ◽  
Chaoying Zhou

Fundamental research on interaction between flow and structure is presented for computation the fluid dynamics of different two-dimensional oscillating models. The Navier-Stokes equations are solved using finite volume method. A multigrid mesh method which was applied to the situation of flow past the stagnating or vibrating cylinder is developed to simulate this type of flow. The interactive results between flow and structure rigid cylinders have been present. The computation fluid dynamic codes mainly with low Reynolds RANS solver are used to solve the impressible viscous Navier-Stokes equations. Finite volume method which is coupled with conformal hybrid mesh method is developed to simulate this type of flow. Numerical investigation focused on the response and the fluid forces on the cylinders and also observed the different shedding model in the wake. The numerical results are compared in detail with recent experimental and computational work. Present numerical comparison also showed that solution using different turbulence model will make the result have a little discrepancy and each turbulence model has respective characteristics in numerical solution on the vortex-induced vibration of the cylinder. In addition, the formation of the 2P vortex shedding model through the lock-in region and the beginning of the shedding model transformation in numerical calculation from 2S model to 2P model has been analyzed.


Author(s):  
Juan B. V. Wanderley ◽  
Gisele H. B. Souza ◽  
Carlos Levi

Numerical simulations of Vortex Induced Vibration have been failing to duplicate accurately experimental data mostly due to the complexity of the physics involved in the real problem. Therefore, a careful and comprehensive investigation on CFD algorithms is still required to indicate the most suitable numerical scheme to handle such a complicate problem. Grid generation, boundary condition implementation, and coupling between the fluid flow governing equations and body motion equation are known to have strong influence on the qualities of the numerical results. This work presents results obtained from a long-term investigation featuring different CFD methods. The investigations enabled the selection of a very effective algorithm that showed an outstanding agreement between experiment and numerical simulation of the VIV phenomenon. Good agreement is obtained in the entire range of reduced velocity covered by the experimental investigations. The successful algorithm discussed here applies the Beam and Warming implicit scheme to solve the two-dimensional slightly compressible Navier–Stokes equations with the K-ε turbulence model to simulate the turbulent flow at the wake of the cylinder.


1996 ◽  
Vol 451 ◽  
Author(s):  
G. Marshall ◽  
P. Mocskos ◽  
F. Molina ◽  
S. Dengra

ABSTRACTRecent work demonstrates the relevant influence of convection during growth pattern formation in thin-layer electrochemical deposition. Convection is driven mainly by coulombic forces due to local charges at the tip of the aggregation and by buoyancy forces due to concentration gradients. Here we study through physical experiments and numerical modeling the regime under which coulombic forces are important. In the experimental measurements fluid motion near the growing tips of the deposit is visualized with neutrally buoyant latex spheres and its speed measured with videomicroscope tracking techniques and image processing software. The numerical modeling consists in the solution of the 2D dimensionless Nernst-Planck equations for ion concentrations, the Poisson equation for the electric field and the Navier-Stokes equations for the fluid flow, and a stochastic growth rule for ion deposition. A new set of dimensionless numbers governing electroconvection dominated flows is introduced. Preliminary experimental measurements and numerical results indicate that in the electroconvection dominated regime coulombic forces increase with the applied voltage, and their influence over growth pattern formation can be assessed with the magnitude of the dimensionless electric Froude number. It is suggested that when this number decreases the deposit morphology changes from fractal to dense branching.


Author(s):  
Xin Guo ◽  
Di Yang ◽  
Yi Liu ◽  
Lian Shen

We perform direct numerical simulations to simulate the interaction between surface waves and the turbulence underneath. The Navier–Stokes equations are simulated using a pseudo-spectral method in horizontal directions and a finite-difference method in vertical direction, with fully nonlinear viscous free-surface kinematic and dynamic boundary conditions at the free surface. We set up the turbulence and the waves by a random forcing method in the bulk flow and a pressure forcing method at the surface, which were recently developed by [1]. It is found that there are surface waves generated on the free surface due to the excitation by the turbulence. The surface elevation is sensitive to the effect of gravity and surface tension. In the presence of progressive waves at the free surface, the turbulent vortical structure is turned, stretched, and compressed periodically by the strain field of waves.


2019 ◽  
Vol 123 (1266) ◽  
pp. 1170-1192 ◽  
Author(s):  
M. M. De ◽  
J. S. Mathur ◽  
S. Vengadesan

ABSTRACTOrnithopters and entomopters should be insensitive to the gusty environment during outdoor operations. Hence, it becomes imperative to understand their behaviour under the influence of gust for ensuring stable flight. In light of this, the present numerical study focused on understanding the aerodynamics of flapping wings with five different planform shapes under the influence of a spatiotemporally varying frontal gust. 3D, unsteady, laminar, and incompressible Navier-Stokes equations were solved using finite volume formulation. A canonical case of asymmetric 1 degree of freedom (DoF) flapping kinematics was considered. Horizontal and vertical force patterns in constant and gusty inflow conditions were numerically computed and compared. Findings were analyzed quantitatively by comparing the differences in the instantaneous force patterns, ordinal scoring approach, and phase space plots. Qualitative comparisons were made based on plots of vortex structures and surface pressure contours for constant and gusty inflow conditions for wings with different planform shapes. Spanwise Lagrangian Coherent Structures (LCS) of all the five wings were also compared. Studies revealed that the elliptical wing exhibited low sensitivity and inverse semi-elliptical wing exhibited high sensitivity to the gusty inflow. Rectangular, triangular and semi-elliptical shaped wings were moderately sensitive to the gusty inflow. This finding, within the limitations of the flapping kinematics and simulation conditions considered for the present study, supported the fact that many natural flyers like forest raptors, non-migratory passerines, pheasants, and partridges have adopted elliptical wing planform for efficient flight.


Author(s):  
John Dunham ◽  
Georges Meauzé

Computer codes which solve the Reynolds-averaged Navier-Stokes equations are now used by manufacturers to design turbomachines, but there is no consensus among experts about which grids and which turbulence models are good enough to provide a reliable basis for design decisions. The AGARD Propulsion and Energetics Panel set up a Working Group to help to clarify these issues, by analysing predictions (using as wide a range of codes as possible) of two representative but difficult single blade row test cases: NASA Rotor 37 and an annular turbine cascade tested by DLR. This paper summarises the Group’s results and conclusions. Recommendations are made about the type and density of grid, which depend on many factors. Although mixing-length turbulence models give good results for quasi-two-dimensional boundary layers, they are essentially unsuitable for turbomachines with their complex end wall flows; it is essential to adopt some kind of turbulent transport model.


1995 ◽  
Vol 117 (2) ◽  
pp. 227-233 ◽  
Author(s):  
Dartzi Pan ◽  
Yu-Chi Chin ◽  
Chih-Hao Chang

The vortex lock-in in the laminar wake behind a circular cylinder induced by the unsteady monopole source is numerically simulated in this paper. The artificial compressibility method is employed to solve the incompressible Navier-Stokes equations. A high-order accurate upwind flux-difference finite-volume scheme is used to discretize the flow field. The unsteady monopole source is simulated by a pulsating volume flux through the cylinder surface at a prescribed forcing frequency and amplitude. The forcing amplitude is set to a fixed value while the frequency is varied to search for the lock-in region. The flow field of the periodic lock-in state is examined in detail. Finally, the effects of a higher amplitude and a different source location are briefly investigated.


Sign in / Sign up

Export Citation Format

Share Document