skip to main content


Title: Verification of DES for Flow over Rigidly and Elastically-Mounted Circular Cylinders in Normal and Yawed Flow
A computational approach based on a k-ω delayed detached eddy simulation model for predicting aerodynamic loads on a smooth circular cylinder is verified against experiments. Comparisons with experiments are performed for flow over a rigidly mounted (static) cylinder and for an elastically-mounted rigid cylinder oscillating in the transverse direction due to vortex-induced vibration (VIV). For the static cases, measurement data from the literature is used to validate the predictions for normally incident flow. New experiments are conducted as a part of this study for yawed flow, where the cylinder axis is inclined with respect to the inflow velocity at the desired yaw angle, β = 30◦. Good agreement is observed between the predictions and measurements for mean and rms surface pressure. Three yawed flow cases (β = 15◦, 30◦, & 45◦) are simulated and the results are found to be independent of β (independence principle) when the flow speed normal to the cylinder axis is selected as the reference velocity scale. Dynamic (VIV) simulations for an elastically-mounted rigid cylinder are performed by coupling the flow solver with a solid dynamics solver where the cylinder motion is modeled as a mass–spring–damper system. The simulations accurately predict the displacement amplitude and unsteady loading over a wide range of reduced velocity, including the region where ‘‘lock-in’’ (synchronization) occurs. VIV simulations are performed at two yaw angles, β = 0◦ and 45◦ and the independence principle is found to be valid over the range of reduced velocities tested with a slightly higher discrepancy when the vortex shedding frequency is close to the natural frequency of the system.  more » « less
Award ID(s):
1537917
NSF-PAR ID:
10203863
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of fluids and structures
Volume:
94
ISSN:
0889-9746
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Cables of suspension, cable-stayed and tied-arch bridges, suspended roofs, and power transmission lines are prone to moderate to large-amplitude vibrations in wind because of their low inherent damping. Structural or fatigue failure of a cable, due to these vibrations, pose a significant threat to the safety and serviceability of these structures. Over the past few decades, many studies have investigated the mechanisms that cause different types of flow-induced vibrations in cables such as rain-wind induced vibration (RWIV), vortex-induced vibration (VIV), iced cable galloping, wake galloping, and dry-cable galloping that have resulted in an improved understanding of the cause of these vibrations. In this study, the parameters governing the turbulence-induced (buffeting) and motion-induced wind loads (self-excited) for inclined and yawed dry cables have been identified. These parameters facilitate the prediction of their response in turbulent wind and estimate the incipient condition for onset of dry-cable galloping. Wind tunnel experiments were performed to measure the parameters governing the aerodynamic and aeroelastic forces on a yawed dry cable. This study mainly focuses on the prediction of critical reduced velocity 〖(RV〗_cr) as a function of equivalent yaw angle (*) and Scruton number (Sc) through measurement of aerodynamic- damping and stiffness. Wind tunnel tests using a section model of a smooth cable were performed under uniform and smooth/gusty flow conditions in the AABL Wind and Gust Tunnel located at Iowa State University. Static model tests for equivalent yaw angles of 0º to 45º indicate that the mean drag coefficient 〖(C〗_D) and Strouhal number (St) of a yawed cable decreases with the yaw angle, while the mean lift coefficient 〖(C〗_L) remains zero in the subcritical Reynolds number (Re) regime. Dynamic one degree-of-freedom model tests in across-wind and along-wind directions resulted in the identification of buffeting indicial derivative functions and flutter derivatives of a yawed cable for a range of equivalent yaw angles. Empirical equations for mean drag coefficient, Strouhal number, buffeting indicial derivative functions and critical reduced velocity for dry-cable galloping are proposed for yawed cables. The results indicate a critical equivalent yaw angle of 45° for dry-cable galloping. A simplified design procedure is introduced to estimate the minimum damping required to arrest dry-cable galloping from occurring below the design wind speed of the cable structure. Furthermore, the results from this study can be applied to predict the wind load and response of a dry cable at a specified wind speed for a given yaw angle. 
    more » « less
  2. Motivated by the need for compact descriptions of the evolution of non-classical wakes behind yawed wind turbines, we develop an analytical model to predict the shape of curled wakes. Interest in such modelling arises due to the potential of wake steering as a strategy for mitigating power reduction and unsteady loading of downstream turbines in wind farms. We first estimate the distribution of the shed vorticity at the wake edge due to both yaw offset and rotating blades. By considering the wake edge as an ideally thin vortex sheet, we describe its evolution in time moving with the flow. Vortex sheet equations are solved using a power series expansion method, and an approximate solution for the wake shape is obtained. The vortex sheet time evolution is then mapped into a spatial evolution by using a convection velocity. Apart from the wake shape, the lateral deflection of the wake including ground effects is modelled. Our results show that there exists a universal solution for the shape of curled wakes if suitable dimensionless variables are employed. For the case of turbulent boundary layer inflow, the decay of vortex sheet circulation due to turbulent diffusion is included. Finally, we modify the Gaussian wake model by incorporating the predicted shape and deflection of the curled wake, so that we can calculate the wake profiles behind yawed turbines. Model predictions are validated against large-eddy simulations and laboratory experiments for turbines with various operating conditions. 
    more » « less
  3. Effect of airfoil thickness on onset of dynamic stall is investigated using large eddy simulations at chord-based Reynolds number of 200 000. Four symmetric NACA airfoils of thickness-to-chord ratios of 9 %, 12 %, 15 % and 18 % are studied. The three-dimensional Navier–Stokes solver, FDL3DI is used with a sixth-order compact finite difference scheme for spatial discretization, second-order implicit time integration and discriminating filters to remove unresolved wavenumbers. A constant-rate pitch-up manoeuver is studied with the pitching axis located at the airfoil quarter chord. Simulations are performed in two steps. In the first step, the airfoil is kept static at a prescribed angle of attack ( $=4^{\circ }$ ). In the second step, a ramp function is used to smoothly increase the pitch rate from zero to the selected value and then the pitch rate is held constant until the angle of attack goes past the lift-stall point. The solver is verified against experiments for flow over a static NACA 0012 airfoil. Static simulation results of all airfoil geometries are also compared against XFOIL predictions with a generally favourable agreement. FDL3DI predicts two-stage transition for thin airfoils (9 % and 12 %), which is not observed in the XFOIL results. The dynamic simulations show that the onset of dynamic stall is marked by the bursting of the laminar separation bubble (LSB) in all the cases. However, for the thickest airfoil tested, the reverse flow region spreads over most of the airfoil and reaches the LSB location immediately before the LSB bursts and dynamic stall begins, suggesting that the stall could be triggered by the separated turbulent boundary layer. The results suggest that the boundary between different classifications of dynamic stall, particularly leading edge stall versus trailing edge stall, is blurred. The dynamic-stall onset mechanism changes gradually from one to the other with a gradual change in some parameters, in this case, airfoil thickness. 
    more » « less
  4. The hydrodynamic forces on an oscillating circular cylinder are predicted using neural networks under flow conditions where Vortex-Induced Vibrations (VIV) are known to occur. The derived neural network approximators are then incorporated in a dynamical model that allows prediction of the cylinder motion given flow conditions and initial conditions. Using experimental data, a minimum-least-squares compensator is tuned that includes linear stiffness and damping su-perimposed with a constant force offset. The compensator is decoupled, i.e., with equations in-dependent for each degree of freedom. By applying the neural network approximators and the derived compensator simulated experiments can be performed. These simulated experiments show that the compensator which cancels the linear components and any bias in the hydrody-namic forces effectively stabilizes the VIV motion. To support this time-domain analysis is per-formed along with phase-plane investigations. Maximum Lyapunov exponent analysis is also shown. 
    more » « less
  5. High-resolution large eddy simulations and complementary laboratory experiments using particle image velocimetry were performed to provide a detailed quantitative assessment of flow response to gaps in cylinder arrays. The base canopy consists of a dense array of emergent rigid cylinders placed in a regular staggered pattern. The gaps varied in length from [Formula: see text] to 24, in intervals of 4 d, where d is the diameter of the cylinders. The analysis was performed under subcritical conditions with Froude numbers [Formula: see text] and bulk Reynolds numbers [Formula: see text]. Results show that the gaps affect the flow statistics at the upstream and downstream proximity of the canopy. The affected zone was [Formula: see text] for the mean flow and [Formula: see text] for the second-order statistics. Dimensionless time-averaged streamwise velocity within the gap exhibited minor variability with gap spacing; however, in-plane turbulent kinetic energy, k, showed a consistent decay rate when normalized with that at [Formula: see text] from the beginning of the gap. The emergent canopy acts as a passive turbulence generator for the gap flow for practical purposes. The streamwise dependence of k follows an exponential trend within [Formula: see text] and transitions to a power-law at [Formula: see text]. The substantially lower maximum values of k within the gap compared to k within the canopy evidence a limitation of gap measurements representative of canopy flow statistics. We present a base framework for estimating representative in-canopy statistics from measurements in the gap. 
    more » « less