• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.6
• /
• pp.786-798
• /
• 2001

This paper describes the phenomena of wake-induced transition of the boundary layers on a NACA0012 airfoil using measured phase-averaged data. Especially, the phase-averaged wall shear stresses are reasonably evaluated using the principle of Computational Preston Tube Method. Due to the passing wake, the turbulent patch is generated in the laminar boundary layer on the airfoil and the boundary layer becomes temporarily transitional. The patches propagate downstream with less speed than free-stream velocity and merge with each other at further down stream station, and the boundary layer becomes more transitional. The generation of turbulent patch at the leading edge of the airfoil mainly depends on velocity defects and turbulent intensity profiles of passing wakes. However, the growth and merging of turbulent patches depend on local streamwise pressure gradients as well as characteristics of turbulent patches. In this transition process, the present experimental data show very similar features to the previous numerical and experimental studies. It is confirmed that the two phase-averaged mean velocity dips appear in the outer region of transitional boundary layer for each passing cycle. Relatively high values of the phase-averaged turbulent fluctuations in the outer region indicate the possibility that breakdown occurs in the outer layer not near the wall.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.10
• /
• pp.1417-1426
• /
• 2001

Turbulent flow characteristics in the near wake of a square cylinder have been studied experimentally by using a Digital PIV method. Experiments are performed at the Reynolds numbers of 1600 and 3900 based on the free-stream velocity and the square height. The ensemble averaged turbulence statistics are acquired from 2030 realizations of instantaneous fluctuating velocity field after the conventional Reynolds decomposition. The differences in turbulent intensity and Reynolds shear stress profiles fur both oases indicate that the effect of Reynolds number seems to be descernible mainly due to the occurrence of transition in the separated shear layer. Because of the periodic nature of vortex shedding process, transverse velocity fluctuations contribute dominantly , to turbulent kinetic energy distribution. A comparison with previous LDV data obtained at much higher Reynolds number shows a fairly good agreement each other. It turns out that the effect of Reynolds number diminishes as increasing Reynolds number, which is a well-known feature of a sharp-edged bluff body wake. The streamwise variation of turbulence intensities are compared with those from a circular cylinder along the centerline at the same Reynolds number. The overall magnitudes and the decay rates of turbulence intensities are quite similar, but some differences are noticeble especially in the transverse intensity variation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.12
• /
• pp.1821-1830
• /
• 2001

The flow around free end of a finite circular cylinder (FC) embedded in an atmospheric boundary layer has been investigated experimentally. The experiments were carried out in a closed-return type subsonic wind tunnel with varying aspect ratio of the finite cylinder mounted vertically on a flat plate. The wakes behind a 2-D cylinder and a finite cylinder located in a uniform flow were measured for comparison. Reynolds number based on the cylinder diameter was about Re=20,000. A hot-wire anemometer was employed to measure the wake velocity and the mean pressure distributions on the cylinder surface were also measured. The flow past the FC free end shows a complicated three-dimensional wake structure and flow phenomenon is quite different from that of 2-D cylinder. The three-dimensional flow structure was attributed to the downwashing counter rotating vortices separated from the FC free end. As the FC aspect ratio decreases, the vortex shedding frequency decreases and the vortex formation length increases compared to that of 2-D cylinder. Due to the descending counter-rotating twin-vortex, near the FC free end, regular vortex shedding from the cylinder is suppressed and the vortex formation region is hardly distinguished. Around the center of the wake, the mean velocity for the FC located in atmospheric boundary layer has large velocity deficit compared to that of uniform flow.

• Journal of the Society of Naval Architects of Korea
• /
• v.56 no.6
• /
• pp.507-514
• /
• 2019

Computational studies of accelerating flow around 2D Circular Cylinder was performed to investigate characteristics of wake field and drag forces. Previous studies had revealed that drag on the cylindrical body in accelerating flow is much greater than that in the flow with constant velocity; however, the underlying physics on the drag increase has not been clearly investigated. In order to investigate the drag increase and its relationship with wake development, this study employed a finite-volume based CFD code, Fluent 13.0 with k-ω SST model for turbulence effects. Inflows are modeled with varied accelerations from 0.4905 to 9.81m/s². The drag computed in the present study is in good agreement with previous studies, and clearly shows the increase compared to the drag on the body in the flow with constant velocity. The results also show that drag crisis observed at high Reynolds number in the case of the flow with constant velocity is also found in the case of accelerating flow. The analysis for wake and recirculation length shows that conventional vortex shedding does not occur even at high Reynolds number and the drag increase is larger at higher acceleration.

• Wind and Structures
• /
• v.30 no.5
• /
• pp.525-532
• /
• 2020

This work numerically investigates the effects of Reynolds number Re_D (= 100 - 150), cross-sectional aspect ratio AR = ( 0.25 -1.0), and attack angle α (= 0° - 90°) on the forces, Strouhal number, and wake of an elliptical cylinder, where Re_D is based on the freestream velocity and cylinder cross-section height normal to the freestream flow, AR is the ratio of the minor axis to the major axis of the elliptical cylinder, and α is the angle between the cylinder major axis and the incoming flow. At Re_D = 100, two distinct wake structures are identified, namely 'Steady wake' (pattern I) and 'Karman wake followed by a steady wake (pattern II)' when AR and α are varied in the ranges specified. When Re_D is increased to 150, an additional wake pattern, 'Karman wake followed by secondary wake (pattern III)' materializes. Pattern I is characterized by two steady bubbles forming behind the cylinder. Pattern II features Karman vortex street immediately behind the cylinder, with the vortex street transmuting to two steady shear layers downstream. Inflection angle α_i = 32°, 37.5° and 45° are identified for AR = 0.25, 0.5 and 0.75, respectively, where the wake asymmetry is the greatest. The α_i effectively distinguishes the dependence on α and AR of force and vortex shedding frequency at either Re_D. In Pattern III, the Karman street forming behind the cylinder is modified to a secondary vortex street. At a given AR and α, Re_D = 150 renders higher fluctuating lift and Strouhal number than Re_D = 100.

• International Journal of Aeronautical and Space Sciences
• /
• v.18 no.1
• /
• pp.17-27
• /
• 2017

The treatment of rotor wake has been a critical issue in the field of the rotor aerodynamics. This paper presents a new free wake model for the unsteady analysis for a wind turbine. A blade-wake-tower interaction is major source of unsteady aerodynamic loading and noise on the wind turbine. However, this interaction can not be considered in conventional free wake model. Thus, the free wake model named Finite Vortex Element (FVE hereafter) was devised in order to consider the interaction effects. In this new free wake model, the wake-tower interaction was described by dividing one vortex filament into two vortex filaments, when the vortex filament collided with a tower. Each divided vortex filaments were remodeled to make vortex ring and horseshoe vortex to satisfy Kelvin's circulation theorem and Helmholtz's vortex theorem. This model was then used to predict aerodynamic load and wake geometry for the horizontal axis wind turbine. The results of the FVE model were compared with those of the conventional free wake model and the experimental results of SNU wind tunnel test and NREL wind tunnel test under various inflow velocity and yaw condition. The result of the FVE model showed better correlation with experimental data. It was certain that the tower interaction has a strong effect on the unsteady aerodynamic load of blades. Thus, the tower interaction needs to be taken into account for the unsteady load prediction. As a result, this research shows a potential of the FVE for an efficient and versatile numerical tool for unsteady loading analysis of a wind turbine.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.1
• /
• pp.20-29
• /
• 2010

The complicated flow characteristics of upper propeller wake influenced by hull wake are investigated in detail in the present study. A two-frame PIV (particle image velocimetry) technique was employed to visualize the upper propeller wake region. As the upper hull wake affects strongly propeller inflow, upper propeller wake shows much unstable vortical behavior, especially in the tip vortices. Velocity field measurements were conducted in a cavitation tunnel with a simulated hull wake. Generally, the hull wake generated by the hull of a marine ship may cause different loading distributions on the propeller blade in both upper and lower propeller planes. The unstable upper propeller wake caused by the ship's hull is expressed in terms of turbulent kinetic energy (TKE) and is identified by using the proper orthogonal decomposition (POD) method to characterize the coherent flow structure in it. Instabilities appeared in the eigen functions higher than the second one, giving unsteadiness to the downstream flow characteristics. The first eigen mode would be useful to find out the tip vortex positions immersed in the unstable downstream region.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.4
• /
• pp.42-47
• /
• 2011

In this paper, wake equalizing duct (WED) form optimization was carried out using computational fluid dynamics (CFD) techniques. A WED is a ring-shaped flow vane with a foil-type cross-section fitted to a hull in front of the upper propeller area. The main advantage of a WED is the power savings resulting from the uniformity of the velocity distribution on the propeller plane, a reduction in the flow separation at the aft-body, and lift generation with a forward force component on the foil section. This paper intends to evaluate these functions and find an optimized WED form for minimizing the viscous resistance and equalizing the wake distribution. In the optimization process, the study uses four WED parameters: the angle of the section, longitudinal location, and angles of the axes for the half rings against the longitudinal and transverse planes of the ship. KRISO 300K VLCC2 (KVLCC2) is chosen as an example ship to demonstrate the WED optimization. The optimization procedure uses genetic algorithms (GAs), a gradient-based optimizer for the refinement of the solution, and Non-dominated Sorting GA-II(NSGA-II) for Multiobjective Optimization. The results show that the optimized WED can reduce the viscous resistance at the expense of the uniformity of the wake distribution.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.181-186
• /
• 1999

The effects of the frequency of upstream gust on the unsteady boundary characteristics on a downstream blade was simulated by using a Navier-Stokes code. The Navier-Stokes code is based on an unstructured finite volume method and uses a low Reynolds k-e turbulence model to close the momentum equations. The MIT flapping foil experiment set-up is used to simulate the interaction between the upstream wake and a blade. The frequency of the upstream wake is simulated by varying rate of pitching motion of the flapping airfoils. Three reduced frequencies. 3.62. 7.24. and 10.86. are simulated. As the frequency increases, the unsteady fluctuation on the surfaces of the downstream hydrofoil is shown to decrease while the upstream flapper wake has larger first harmonics of y-velocity component. The unsteady vortices are shown to interact with each other and. as a result. the upstream wake becomes undiscernible inside the inner layer. The turbulence kinetic energy shows a similar behavior. Limiting streamlines around the trailing edge of the flapper are shown to conform with the unsteady Kutta condition for a round trailing edge. while limiting streamlines around the trailing edge of the hydrofoil conforms with the unsteady Kutta condition for a sharp edge.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.12 no.1
• /
• pp.501-517
• /
• 2020

Three-dimensional numerical simulations were performed to study the effects of flow direction and flow velocity on the flow regime behind a curved pipe represented by a curved circular cylinder. The cylinder is based on a previous study and consists of a quarter segment of a ring and a horizontal part at the end of the ring. The cylinder was rotated in the computational domain to examine five incident flow angles of 0-180° with 45° intervals at Reynolds numbers of 100 and 500. The detailed wake topologies represented by λ2 criterion were captured using a Large Eddy Simulation (LES). The curved cylinder leads to different flow regimes along the span, which shows the three-dimensionality of the wake field. At a Reynolds number of 100, the shedding was suppressed after flow angle of 135°, and oblique flow was observed at 90°. At a Reynolds number of 500, vortex dislocation was detected at 90° and 135°. These observations are in good agreement with the three-dimensionality of the wake field that arose due to the curved shape.