RANS Simulation Versus Experiments of Flow Induced Motion of Circular Cylinder With Passive Turbulence Control at 35,000 < RE < 130,000

Author(s):  
Wei Wu ◽  
Michael M. Bernitsas ◽  
Kevin Maki

Two-dimensional (2D) Unsteady Reynolds-Averaged Navier–Stokes equations (URANS) equations with the Spalart–Allmaras turbulence model are used to simulate the flow and body kinematics of the transverse motion of spring-mounted circular cylinder. The flow is in the high-lift TrSL3 regime of a Reynolds number in the range 35,000 < Re < 130,000. Passive turbulence control (PTC) in the form of selectively distributed surface roughness is used to alter the cylinder flow induced motion (FIM). Simulation is performed using a solver based on the open source Computational Fluid Dynamics (CFD) tool OpenFOAM, which solves continuum mechanics problems with a finite-volume discretization method. Roughness parameters of PTC are chosen based on tests conducted in the Marine Renewable Energy Lab (MRELab) of the University of Michigan. The numerical tool is first tested on smooth cylinder in vortex-induced vibration (VIV) and results are compared with available experimental measurements and URANS simulations. For the cylinder with PTC cases, the sandpaper grit on the cylinder wall is modeled as a rough-wall boundary condition. Two sets of cases with different system parameters (spring, damping) are simulated and the results are compared with experimental data measured in the MRELab. The amplitude ratio curve shows clearly three different branches, including the VIV initial and upper branches, and a galloping branch. The numerical branches are similar to those observed experimentally. Frequency ratio, vortex patterns, transitional behavior, and lift are also predicted well for PTC cylinders at such high Reynolds numbers.

Author(s):  
Wei Wu ◽  
Michael M. Bernitsas ◽  
Kevin Maki

Two-dimensional RANS equations with the Spalart-Allmaras turbulence model are used to simulate the flow and body kinematics of a rigid circular cylinder mounted on springs, transversely to a steady uniform flow in the high-lift, TrSL3 regime with 35,000<Re<130,000. Passive Turbulence Control (PTC) in the form of selectively distributed surface roughness is used to alter the cylinder Flow Induced Motion (FIM). Simulation is performed by using a solver based on the open source CFD tool OpenFOAM, which solves continuum mechanics problems with a finite volume discretization method. Roughness parameters of PTC are simulated modeling tests conducted in the Marine Renewable Energy Lab (MRELab) of the University of Michigan. The numerical tool is first tested on smooth cylinder in VIV and results are compared with available experimental measurements and RANS simulations. For the cylinder with PTC cases, the sandpaper grit (k) on the cylinder wall is modeled as a rough-wall boundary condition. Two sets of cases with different system parameters (spring constant, damping) are simulated and the results are compared with experimental data measured in the MRELab. The amplitude-ratio curve shows clearly three different branches, including the VIV initial and upper branches and a galloping branch, similar to those observed experimentally. Frequency ratio, vortex patterns, transitional behavior, and lift are also predicted well for PTC cylinders at such high Reynolds numbers.


Author(s):  
Hongrae Park ◽  
Michael M. Bernitsas ◽  
Che-Chun Chang

Passive turbulence control (PTC) in the form of two straight roughness strips with variable width, and thickness about equal to the boundary layer thickness, is used to modify the flow-induced motions (FIM) of a rigid circular cylinder. The cylinder is supported by two end-springs and the flow is in the TrSL3, high-lift, regime. The PTC-to-FIM Map, developed in previous work, revealed zones of weak suppression, strong suppression, hard galloping, and soft galloping. In this paper the sensitivity of the PTC-to-FIM Map to: (a) the width of PTC covering, (b) PTC covering a single or multiple zones, (c) PTC being straight or staggered is studied experimentally. Experiments are conducted in the Low Turbulence Free Surface Water Channel of the University of Michigan. Fixed parameters are: cylinder diameter D = 8.89cm, m* = 1.725, spring stiffness K = 763N/m, aspect ratio l/D = 10.29, and damping ratio ζ = 0.019. Variable parameters are: circumferential PTC location αPTC ∈ [0°−180°], Reynolds number Re ∈ [30,000–120,000], flow velocity U ∈ [0.36m/s–1.45m/s]. Measured quantities are: amplitude ratio A/D, frequency ratio fosc/fn,w, and synchronization range. As long as the roughness distribution is limited to remain within a zone, the width of the strips does not affect the FIM response. When multiple zones are covered, the strong suppression zone dominates the FIM.


1995 ◽  
Vol 303 ◽  
pp. 215-232 ◽  
Author(s):  
H. M. Badr ◽  
S. C. R. Dennis ◽  
S. Kocabiyik ◽  
P. Nguyen

The transient flow field caused by an infinitely long circular cylinder placed in an unbounded viscous fluid oscillating in a direction normal to the cylinder axis, which is at rest, is considered. The flow is assumed to be started suddenly from rest and to remain symmetrical about the direction of motion. The method of solution is based on an accurate procedure for integrating the unsteady Navier–Stokes equations numerically. The numerical method has been carried out for large values of time for both moderate and high Reynolds numbers. The effects of the Reynolds number and of the Strouhal number on the laminar symmetric wake evolution are studied and compared with previous numerical and experimental results. The time variation of the drag coefficients is also presented and compared with an inviscid flow solution for the same problem. The comparison between viscous and inviscid flow results shows a better agreement for higher values of Reynolds and a Strouhal numbers. The mean flow for large times is calculated and is found to be in good agreement with previous predictions based on boundary-layer theory.


Author(s):  
Lin Ding ◽  
Li Zhang ◽  
Chunmei Wu ◽  
EunSoo Kim ◽  
Michael M. Bernitsas

The effect of tandem spacing on the flow induced motions (FIM) of two circular cylinders with passive turbulence control is investigated using two-dimensional Unsteady Reynolds-Average Navier-Stokes equations with the Spalart-Allmaras turbulence model. Results are compared to experiments in the range of Reynolds number of 30,000<Re<100,000. The center-to-center spacing between the two cylinders is varied from 2 to 6 diameters. Simulation results predict well all ranges of FIM including VIV and galloping and match well with experimental measurements. For the upstream cylinder, the amplitude and frequency responses are not considerably influenced by the downstream cylinder when the spacing is greater than 2D. For the downstream cylinder, a rising amplitude trend in the VIV upper branch can be observed in all cases as is typical of flows in the TrSL3 regime. The galloping branch merges with the VIV upper branch for spacing greater than 3D. Vortex structures show significant variation in different flow regimes in accordance with experimental observations. High-resolution post-processing shows that the interaction between the wakes of cylinders result in various types of FIM.


Author(s):  
Hongrae Park ◽  
Eun Soo Kim ◽  
Michael M. Bernitsas

Passive turbulence control (PTC) in the form of two straight roughness strips with variable width, and thickness about equal to the boundary layer thickness, is used to modify the flow-induced motions (FIM) of a rigid circular cylinder. The cylinder is supported by two end springs and the flow is in the TrSL3, high-lift, regime. The PTC-to-FIM Map, developed in the previous work, revealed zones of weak suppression (WS), strong suppression (SS), hard galloping (HG), and soft galloping (SG). In this paper, the sensitivity of the PTC-to-FIM map to: (a) the width of PTC covering, (b) PTC covering a single or multiple zones, and (c) PTC being straight or staggered is studied experimentally. Experiments are conducted in the low turbulence free surface water channel of the University of Michigan, Ann Arbor, MI. Fixed parameters are: cylinder diameter D = 8.89 cm, m* = 1.725, spring stiffness K = 763 N/m, aspect ratio l/D = 10.29, and damping ratio ζ = 0.019. Variable parameters are circumferential PTC location αPTC∈ (0–180 deg), Reynolds number Re ∈ (30,000–120,000), flow velocity U∈ (0.36–1.45 m/s). Measured quantities are amplitude ratio A/D, frequency ratio fosc/fn,w, and synchronization range. As long as the roughness distribution is limited to remain within a zone, the width of the strips does not affect the FIM response. When multiple zones are covered, the strong suppression zone dominates the FIM.


Author(s):  
Lin Ding ◽  
Li Zhang ◽  
Chunmei Wu ◽  
Eun Soo Kim ◽  
Michael M. Bernitsas

The effect of tandem spacing on the flow-induced motions (FIM) of two circular cylinders with passive turbulence control is investigated using two-dimensional (2D) unsteady Reynolds-averaged Navier–Stokes equations with the Spalart–Allmaras turbulence model. Results are compared to experiments in the range of Reynolds number of 30,000 < Re < 100,000. The center-to-center spacing between the two cylinders is varied from 2 to 6 diameters. Simulation results predict well all the ranges of FIM including vortex-induced vibrations (VIV) and galloping and match well with experimental measurements. For the upstream cylinder, the amplitude and frequency responses are not considerably influenced by the downstream cylinder when the spacing is greater than 2D. For the downstream cylinder, a rising amplitude trend in the VIV upper-branch can be observed in all the cases as is typical of flows in the TrSL3 flow regime (transition in shear layer 3; 2 × 104 < Re < 3 × 105). The galloping branch merges with the VIV upper-branch for spacing greater than three-dimensional (3D). Vortex structures show significant variation in different flow regimes in accordance with experimental observations. High-resolution postprocessing shows that the interaction between the wakes of cylinders results in various types of FIM.


2020 ◽  
Author(s):  
Ningyu Li ◽  
Hongrae Park ◽  
Hai Sun ◽  
Michael M. Bernitsas

Abstract Passive turbulence control (PTC) is being used in the Marine Renewable Energy Laboratory (MRELab) of the University of Michigan to enhance flow induced oscillations (FIO) of cylinders in the VIVACE (Vortex Induced Vibration for Aquatic Clean Energy) Converter. Large PTC triggers VIV and galloping at lower flow speeds for energy harvesting. Currently, FIO of cylinders with large PTC for high Re has received limited attention and, particularly, the effect of variable PTC height on FIO of cylinders. The vast majority of ocean currents, rivers, and tides are too slow for Marine Hydro Kinetic (MHK) energy technologies to harness it. In order to enhance FIO and to initiate galloping earlier, a circular cylinder is geometrically modified using straight strips placed on the cylinder surface symmetrically PTC strips on the cylinder effectively change the flow properties. In the present study, the FIO of a single-cylinder with large PTC, on end linear-springs, is modelled and simulated using a Fluid-Structure Interaction (FSI) code. Results are verified by corresponding experimental data. Results show that VIV onset occurs at lower Re for large-PTC cylinder in comparison with lower-PTC cylinder. Contrary to smooth cylinders for which the amplitude ratio is small in the transition region between VIV and galloping, application of large PTC leads to high amplitude response in the transition region. The mechanism behind this observation is the further departure of the geometry from the smooth circular cylinder. The latter does not exhibit galloping due to flow and geometric symmetry in all directions. Moreover, in the galloping region, the amplitude ratio increases with the height of PTC. Earlier onset of galloping and enhancement of geometric asymmetry support this observation as well.


2000 ◽  
Vol 407 ◽  
pp. 123-144 ◽  
Author(s):  
S. C. R. DENNIS ◽  
P. NGUYEN ◽  
SERPIL KOCABIYIK

The temporal development of two-dimensional viscous incompressible flow induced by an impulsively started circular cylinder which performs time-dependent rotational oscillations about its axis and translates at right angles to this axis is investigated. The investigation is based on the solutions of the unsteady Navier–Stokes equations. A series expansion for small times is developed. The Navier–Stokes equations are also integrated by a spectral–finite difference method for moderate values of time for both moderate and high Reynolds numbers. The numerical method is checked with the results of the analytical solution. The effects of the Reynolds number and of the forcing Strouhal number S on the laminar asymmetric flow structure in the near-wake region are studied. The lift and drag coefficients are also extracted from numerical results. An interesting phenomenon has been observed both in the flow patterns and in the behaviour of drag coefficients for S = π/2 at Reynolds number R = 500 and is discussed. For comparison purposes the start-up flow is determined numerically at a low Reynolds number and is found to be in good agreement with previous experimental predictions.


Author(s):  
V. Tamimi ◽  
M. Zeinoddini ◽  
A. Bakhtiari ◽  
M. Golestani

In this paper results from simulating the vortex shedding phenomena behind a fixed tapered circular cylinder, at relatively high Reynolds numbers, are reported. Ansys-CFX computational fluid dynamics model, based on solving three-dimensional (3D) incompressible transient Navier Stokes equations, is employed for this purpose. The geometries applied in the models resemble those used in wind tunnel experiments by other researchers. The taper slope along the cylinder span is uniform with a tangent of 24:1. The diameter at mid-span of the cylinder equals to 0.0389 m. The Reynolds number (based on the mid-span diameter) is around 29,000. The computational model has first been calibrated against experiments for uniform 3D cylinders as well as results from a Direct Numerical Simulation of turbulent wake with vortex shedding past a uniform circular cylinder, as obtained by other researchers. The main flow characteristics for tapered cylinders such as vortex dislocations and splitting, cellular vortex shedding, oblique vortex shedding and the variation of the vorticity patterns along the tapered cylinder could be obtained from the simulations.


2013 ◽  
Vol 135 (7) ◽  
Author(s):  
A. Ghidoni ◽  
A. Colombo ◽  
S. Rebay ◽  
F. Bassi

In the last decade, discontinuous Galerkin (DG) methods have been the subject of extensive research efforts because of their excellent performance in the high-order accurate discretization of advection-diffusion problems on general unstructured grids, and are nowadays finding use in several different applications. In this paper, the potential offered by a high-order accurate DG space discretization method with implicit time integration for the solution of the Reynolds-averaged Navier–Stokes equations coupled with the k-ω turbulence model is investigated in the numerical simulation of the turbulent flow through the well-known T106A turbine cascade. The numerical results demonstrate that, by exploiting high order accurate DG schemes, it is possible to compute accurate simulations of this flow on very coarse grids, with both the high-Reynolds and low-Reynolds number versions of the k-ω turbulence model.


Sign in / Sign up

Export Citation Format

Share Document