Vortex shedding from a row of square bars

Interactions of wakes in a flow past a row of square bars, which is placed across a uniform flow, are investigated by numerical simulations and experiments on the tassumption that the flow is two-dimensional and incompressible. At small Reynolds numbers the flow is steady and symmetric with respect not only to streamwise lines through the center of each square bar but also to streamwise centerlines between adjacent square bars. However, the steady symmetric flow becomes unstable at larger Reynolds numbers and make a transition to a steady asymmetric flow with respect to the centerlines between adjacent square bars in some cases or to an oscillatory flow in other cases. It is found that vortices are shed synchronously from adjacent square bars in the same phase or in anti-phase depending upon the distance between the bars when the flow is oscillatory. The origin of the transition to the steady asymmetric flow is identified as a pitchfork bifurcation, while the oscillatory flows with synchronous shedding of vortices are clarified to originate from a Hopf bifurcation. The critical Reynolds numbers of the transitions are evaluated numerically and the bifurcation diagram of the flow is obtained.

Vortex Shedding Dynamics in Long Aspect-Ratio Aerodynamics Bodies

The focus of the current dissertation is to study qualitatively the underlying physics of vortex-shedding and wake dynamics in long aspect-ratio aerodynamics in incompressible viscous flow through the use of the KLE method. We carried out a long series of numerical experiments in the cases of flow around the cylinder at low Reynolds numbers. The study of flow at low Reynolds numbers provides an insight in the fluid physics and also plays a critical role when applying to stalled turbine rotors. Many of the conclusions about the qualitative nature of the physical mechanisms characterizing vortex formation, shedding and further interaction analyzed here at low Re could be extended to other Re regimes and help to understand the separation of the boundary layers in airfoils and other aerodynamic surfaces. In the long run, it aims to provide a better understanding of the complex multi-physics problems involving fluid-structure-control interaction through improved mathematical computational models of the multi-physics process. Besides the scientific conclusions produced, the research work on streamlined and bluff-body condition will also serve as a valuable guide for the future design of blade aerodynamics and the placement of wind turbines and hydrakinetic turbines, increasing the efficiency in the use of expensive workforce, supplies, and infrastructure. After the introductory section describing the main fields of application of wind power and hydrokinetic turbines, we describe the main features and theoretical background of the numerical method used here. Then, we present the analysis of the numerical experimentation results for the oscillatory regime right before the onset of vortex shedding for circular cylinders. We verified the wake length of the closed near-wake behind the cylinder and analysed the decay of the wake at the wake formation region, and then studied the St-Re relationship at the Reynolds numbers before the wake sheds compared to the experimental data. We found a theoretical model that describes the time evolution of the amplitude of fluctuations in the vorticity field on the twin vortex wake, which accurately matches the numerical results in terms of the frequency of the oscillation and rate of decay. We also proposed a model based on an analog circuit that is able to interpret the concerning flow by reducing the number of degrees of freedom. It follows the idea of the non-linear oscillator and resembles the dynamics mechanism of the closed near-wake with a common configured sine wave oscillator. This low-dimensional circuital model may also help to understand the underlying physical mechanisms, related to vorticity transport, that give origin to those oscillations.

3D CFD Simulation of Vortex-induced Vibration of Cylinder

The 3D flow around a circular cylinder free to oscillate transversely to the free stream was simulated using Computational Fluid Dynamics (CFD) and the Spalart-Allmaras Detached Eddy Simulation (DES) turbulence model for a Reynolds number Re = 10(4). Simulations were carried out for a small mass-damping parameter m*zeta = 0.00858, where m* = 3.3 and zeta = 0.0026. We found good agreement between the numerical results and experimental data. The simulations predicted the high observed amplitudes of the upper branch of vortex-induced vibrations for low mass-damping parameters.

Volumetric reconstruction of the mean flow around circular cylinders fitted with strakes

The volumetric reconstruction technique presented in this paper employs a two-camera stereoscopic particle image velocimetry (SPIV) system in order to reconstruct the mean flow behind a fixed cylinder fitted with helical strakes, which are commonly used to suppress vortex-induced vibrations (VIV). The technique is based on the measurement of velocity fields at equivalent adjacent planes that results in pseudo volumetric fields. The main advantage over proper volumetric techniques is the avoidance of additional equipment and complexity. The averaged velocity fields behind the straked cylinders and the geometrical periodicity of the three-start configuration are used to further simplify the reconstruction process. Two straked cylindrical models with the same pitch (p = 10d) and two different heights (h = 0.1 and 0.2d) are tested. The reconstructed flow shows that the strakes introduce in the wake flow a well-defined wavelength of one-third of the pitch. Measurements of hydrodynamic forces, fluctuating velocity, vortex formation length, and vortex shedding frequency show the interdependence of the wake parameters. The vortex formation length is increased by the strakes, which is an important effect for the suppression of vortex-induced vibrations. The results presented complement previous investigations concerning the effectiveness of strakes as VIV suppressors and provide a basis of comparison to numerical simulations.

Experimental investigation of flow-induced vibration on isolated and tandem circular cylinders fitted with strakes

The effect of varying the geometric parameters of helical strakes on vortex-induced vibration (VIV) is investigated in this paper. The degree of oscillation attenuation or even suppression is analysed for isolated circular cylinder cases. How a cylinder fitted with strakes behaves when immersed in the wake of another cylinder in tandem arrangement is also investigated and these results are compared to those with a single straked cylinder. The experimental tests are conducted at a circulating water channel facility and the cylindrical models are mounted on a low-damping air bearing elastic base with one degree-of-freedom, restricted to oscillate in the transverse direction to the channel flow. Three strake pitches (p) and heights (h) are tested: p = 5, 10, 15d, and h = 0.1, 0.2, 0.25d. The mass ratio is 1.8 for all models. The Reynolds number range is from 1000 to 10000, and the reduced velocity varies up to 21. The cases with h = 0.1d strakes reduce the amplitude response when compared to the isolated plain cylinder, however the oscillation still persists. On the other hand, the cases with h = 0.2, 0.25d strakes almost completely suppress VIV. Spanwise vorticity fields, obtained through stereoscopic digital particle image velocimetry (SDPIV), show an alternating vortex wake for the p = 10d and h = 0.1d straked cylinder. The p = 10d and h = 0.2d cylinder wake has separated shear layers with constant width and no roll-up close to the body. The strakes do not increase the magnitude of the out-of-plane velocity compared to the isolated plain cylinder. However, they deflect the flow in the out-of-plane direction in a controlled way, which can prevent the vortex shedding correlation along the span. In order to investigate the wake interference effect on the strake efficiency, an experimental arrangement with two cylinders in tandem is employed. The centre-to-centre distance for the tandem arrangement varies from 2 to 6. When the downstream p = 10d and h = 0.2d cylinder is immersed in the wake of an upstream fixed plain cylinder, it loses its effectiveness compared with the isolated case. Although the oscillations have significant amplitude, they are limited, which is a different behaviour from that of a tandem configuration with two plain cylinders. For this particular case, the amplitude response monotonically increases for all gaps, except one, a trait usually found in galloping-like oscillations. SDPIV results for the tandem arrangements show alternating vortex shedding and oscillatory wake. (C) 2010 Elsevier Ltd. All rights reserved.

Wind Loads on Dynamic Message Cabinets and Behavior of Supporting Trusses, TR-612, 2013

Large Dynamic Message Signs (DMSs) have been increasingly used on freeways, expressways and major arterials to better manage the traffic flow by providing accurate and timely information to drivers. Overhead truss structures are typically employed to support those DMSs allowing them to provide wider display to more lanes. In recent years, there is increasing evidence that the truss structures supporting these large and heavy signs are subjected to much more complex loadings than are typically accounted for in the codified design procedures. Consequently, some of these structures have required frequent inspections, retrofitting, and even premature replacement. Two manufacturing processes are primarily utilized on truss structures - welding and bolting. Recently, cracks at welding toes were reported for the structures employed in some states. Extremely large loads (e.g., due to high winds) could cause brittle fractures, and cyclic vibration (e.g., due to diurnal variation in temperature or due to oscillations in the wind force induced by vortex shedding behind the DMS) may lead to fatigue damage, as these are two major failures for the metallic material. Wind and strain resulting from temperature changes are the main loads that affect the structures during their lifetime. The American Association of State Highway and Transportation Officials (AASHTO) Specification defines the limit loads in dead load, wind load, ice load, and fatigue design for natural wind gust and truck-induced gust. The objectives of this study are to investigate wind and thermal effects in the bridge type overhead DMS truss structures and improve the current design specifications (e.g., for thermal design). In order to accomplish the objective, it is necessary to study structural behavior and detailed strain-stress of the truss structures caused by wind load on the DMS cabinet and thermal load on the truss supporting the DMS cabinet. The study is divided into two parts. The Computational Fluid Dynamics (CFD) component and part of the structural analysis component of the study were conducted at the University of Iowa while the field study and related structural analysis computations were conducted at the Iowa State University. The CFD simulations were used to determine the air-induced forces (wind loads) on the DMS cabinets and the finite element analysis was used to determine the response of the supporting trusses to these pressure forces. The field observation portion consisted of short-term monitoring of several DMS Cabinet/Trusses and long-term monitoring of one DMS Cabinet/Truss. The short-term monitoring was a single (or two) day event in which several message sign panel/trusses were tested. The long-term monitoring field study extended over several months. Analysis of the data focused on trying to identify important behaviors under both ambient and truck induced winds and the effect of daily temperature changes. Results of the CFD investigation, field experiments and structural analysis of the wind induced forces on the DMS cabinets and their effect on the supporting trusses showed that the passage of trucks cannot be responsible for the problems observed to develop at trusses supporting DMS cabinets. Rather the data pointed toward the important effect of the thermal load induced by cyclic (diurnal) variations of the temperature. Thermal influence is not discussed in the specification, either in limit load or fatigue design. Although the frequency of the thermal load is low, results showed that when temperature range is large the restress range would be significant to the structure, especially near welding areas where stress concentrations may occur. Moreover stress amplitude and range are the primary parameters for brittle fracture and fatigue life estimation. Long-term field monitoring of one of the overhead truss structures in Iowa was used as the research baseline to estimate the effects of diurnal temperature changes to fatigue damage. The evaluation of the collected data is an important approach for understanding the structural behavior and for the advancement of future code provisions. Finite element modeling was developed to estimate the strain and stress magnitudes, which were compared with the field monitoring data. Fatigue life of the truss structures was also estimated based on AASHTO specifications and the numerical modeling. The main conclusion of the study is that thermal induced fatigue damage of the truss structures supporting DMS cabinets is likely a significant contributing cause for the cracks observed to develop at such structures. Other probable causes for fatigue damage not investigated in this study are the cyclic oscillations of the total wind load associated with the vortex shedding behind the DMS cabinet at high wind conditions and fabrication tolerances and induced stresses due to fitting of tube to tube connections.

Numerical Investigation of Flow Past a Circular Cylinder and in a Staggered Tube Bundle Using Various Turbulence Models

Transitional flow past a three-dimensional circular cylinder is a widely studied phenomenon since this problem is of interest with respect to many technical applications. In the present work, the numerical simulation of flow past a circular cylinder, performed by using a commercial CFD code (ANSYS Fluent 12.1) with large eddy simulation (LES) and RANS (κ - ε and Shear-Stress Transport (SST) κ - ω! model) approaches. The turbulent flow for ReD = 1000 & 3900 is simulated to investigate the force coefficient, Strouhal number, flow separation angle, pressure distribution on cylinder and the complex three dimensional vortex shedding of the cylinder wake region. The numerical results extracted from these simulations have good agreement with the experimental data (Zdravkovich, 1997). Moreover, grid refinement and time-step influence have been examined. Numerical calculations of turbulent cross-flow in a staggered tube bundle continues to attract interest due to its importance in the engineering application as well as the fact that this complex flow represents a challenging problem for CFD. In the present work a time dependent simulation using κ – ε, κ - ω! and SST models are performed in two dimensional for a subcritical flow through a staggered tube bundle. The predicted turbulence statistics (mean and r.m.s velocities) have good agreement with the experimental data (S. Balabani, 1996). Turbulent quantities such as turbulent kinetic energy and dissipation rate are predicted using RANS models and compared with each other. The sensitivity of grid and time-step size have been analyzed. Model constants sensitivity study have been carried out by adopting κ – ε model. It has been observed that model constants are very sensitive to turbulence statistics and turbulent quantities.

Stochastic field processes in the mathematical modelling of damped transmission line vibrations

Wind-excited vibrations in the frequency range of 10 to 50 Hz due to vortex shedding often cause fatigue failures in the cables of overhead transmission lines. Damping devices, such as the Stockbridge dampers, have been in use for a long time for supressing these vibrations. The dampers are conveniently modelled by means of their driving point impedance, measured in the lab over the frequency range under consideration. The cables can be modelled as strings with additional small bending stiffness. The main problem in modelling the vibrations does however lay in the aerodynamic forces, which usually are approximated by the forces acting on a rigid cylinder in planar flow. In the present paper, the wind forces are represented by stochastic processes with arbitrary crosscorrelation in space; the case of a Kármán vortex street on a rigid cylinder in planar flow is contained as a limit case in this approach. The authors believe that this new view of the problem may yield useful results, particularly also concerning the reliability of the lines and the probability of fatigue damages. © 1987.

Estudo do escoamento ao redor de cilindros elípticos

Pós-graduação em Engenharia Mecânica - FEIS

Estudo experimental do escoamento sobre placas planas munidas de protuberâncias quadradas e onduladas

Pós-graduação em Engenharia Mecânica - FEIS

Estudo experimental do escoamento em torno de cilindros circulares em movimento de rotação

Pós-graduação em Engenharia Mecânica - FEIS

FLOW AROUND A SQUARE CYLINDER POSITIONED NEAR A SMOOTH FLAT PLATE

Vortex shedding phenomenon produced by a square cylinder placed close to a smooth flat plate is experimentally studied by means of flow visualization techniques and hot-film anemometry. Qualitative and quantitative information about the flow field has been obtained for Reynolds numbers up to 1,000. Vortex shedding images in several Reynolds number have been captured and the non dimensional vortex shedding frequency has been obtained as a function of the Reynolds number.

Simulación URANS 2D del flujo alrededor de un cilindro para 104 ≤ Re ≤ 106 con tratamiento de paredes

Máster Universitario en Sistemas Inteligentes y Aplicaciones Numéricas en Ingeniería (SIANI)

Experimental analysis of coaxial jets: instability, flow and mixing characterization

The velocity and mixing field of two turbulent jets configurations have been experimentally characterized by means of cold- and hot-wire anemometry in order to investigate the effects of the initial conditions on the flow development. In particular, experiments have been focused on the effect of the separation wall between the two streams on the flow field. The results of the experiments have pointed out that the wake behind a thick wall separating wall has a strong influence on the flow field evolution. For instance, for nearly unitary velocity ratios, a clear vortex shedding from the wall is observable. This phenomenon enhances the mixing between the inner and outer shear layer. This enhancement in the fluctuating activity is a consequence of a local absolute instability of the flow which, for a small range of velocity ratios, behaves as an hydrodynamic oscillator with no sensibility to external perturbations. It has been suggested indeed that this absolute instability can be used as a passive method to control the flow evolution. Finally, acoustic excitation has been applied to the near field in order to verify whether or not the observed vortex shedding behind the separating wall is due to a global oscillating mode as predicted by the theory. A new scaling relationship has been also proposed to determine the preferred frequency for nearly unitary velocity ratios. The proposed law takes into account both the Reynolds number and the velocity ratio dependence of this frequency and, therefore, improves all the previously proposed relationships.

Uso de DES y RANS no estacionario para la simulación del flujo alrededor de una carena Wigley - Use of Unsteady RANS and DES to simulate flows around a Wigley Hull)

Computer Fluid Dynamics tools have already become a valuable instrument for Naval Architects during the ship design process, thanks to their accuracy and the available computer power. Unfortunately, the development of RANSE codes, generally used when viscous effects play a major role in the flow, has not reached a mature stage, being the accuracy of the turbulence models and the free surface representation the most important sources of uncertainty. Another level of uncertainty is added when the simulations are carried out for unsteady flows, as those generally studied in seakeeping and maneuvering analysis and URANS equations solvers are used. Present work shows the applicability and the benefits derived from the use of new approaches for the turbulence modeling (Detached Eddy Simulation) and the free surface representation (Level Set) on the URANS equations solver CFDSHIP-Iowa. Compared to URANS, DES is expected to predict much broader frequency contents and behave better in flows where boundary layer separation plays a major role. Level Set methods are able to capture very complex free surface geometries, including breaking and overturning waves. The performance of these improvements is tested in set of fairly complex flows, generated by a Wigley hull at pure drift motion, with drift angle ranging from 10 to 60 degrees and at several Froude numbers to study the impact of its variation. Quantitative verification and validation are performed with the obtained results to guarantee their accuracy. The results show the capability of the CFDSHIP-Iowa code to carry out time-accurate simulations of complex flows of extreme unsteady ship maneuvers. The Level Set method is able to capture very complex geometries of the free surface and the use of DES in unsteady simulations highly improves the results obtained. Vortical structures and instabilities as a function of the drift angle and Fr are qualitatively identified. Overall analysis of the flow pattern shows a strong correlation between the vortical structures and free surface wave pattern. Karman-like vortex shedding is identified and the scaled St agrees well with the universal St value. Tip vortices are identified and the associated helical instabilities are analyzed. St using the hull length decreases with the increase of the distance along the vortex core (x), which is similar to results from other simulations. However, St scaled using distance along the vortex cores shows strong oscillations compared to almost constants for those previous simulations. The difference may be caused by the effect of the free-surface, grid resolution, and interaction between the tip vortex and other vortical structures, which needs further investigations. This study is exploratory in the sense that finer grids are desirable and experimental data is lacking for large α, especially for the local flow. More recently, high performance computational capability of CFDSHIP-Iowa V4 has been improved such that large scale computations are possible. DES for DTMB 5415 with bilge keels at α = 20º were conducted using three grids with 10M, 48M and 250M points. DES analysis for flows around KVLCC2 at α = 30º is analyzed using a 13M grid and compared with the results of DES on the 1.6M grid by. Both studies are consistent with what was concluded on grid resolution herein since dominant frequencies for shear-layer, Karman-like, horse-shoe and helical instabilities only show marginal variation on grid refinement. The penalties of using coarse grids are smaller frequency amplitude and less resolved TKE. Therefore finer grids should be used to improve V&V for resolving most of the active turbulent scales for all different Fr and α, which hopefully can be compared with additional EFD data for large α when it becomes available.