• International Journal of Aeronautical and Space Sciences
• /
• 제18권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.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제10권2호
• /
• pp.212-224
• /
• 2018

The characteristics of the relationship between the evolution of propeller trailing vortex wake and skew angle are numerically examined based on four different five-blade David Taylor Model Basin (DTMB) model propellers with different skew angles. Numerical simulations are based on Reynolds-averaged Naviere-Stokes (RANS) equations combined with SST $k-{\omega}$ turbulence model. Results show that the contraction of propeller trailing vortex wake can be restrained by increasing skew angle and loading conditions, and root vortices fade away when the propeller skew angle increases. With the increase of the propeller's skew angle, the deformation of the hub vortex and destabilization of the tip vortices are weakening gradually because the blade-to-blade interaction becomes weaker. The transition trailing vortex wake from stability to instability is restrained when the skew increases. Furthermore, analyses of tip vortice trajectories show that the increasing skew can reduce the difference in trailing vortex wake contraction under different loading conditions.

• 대한기계학회논문집B
• /
• 제22권6호
• /
• pp.771-787
• /
• 1998

A Navier-Stokes code with a modified low Reynolds number k-.epsilon. turbulence model was used to study the unsteady transitional boundary layer flow due to rotor-stator interaction. The modification, proposed by Launder, to improve prediction of stagnation flows was incorporated to the low Reynolds number k-.epsilon. turbulence model by Fan-Lakshminarayana-Barnett. Numerical solution is shown to capture well the calmed laminar flow as well as the wake induced transitional strip due to rotor-stator interaction and shows improvement, in terms of onset of transition and its length, over previous Euler/boundary layer solution. The turbulent kinetic energy shows local maximum along the upstream rotor wake in the wake induced transitional strip and this characteristics is observed untill the end of transition. The wake induced strip also shown apparent even in the laminar sublayer as the upstream rotor wake penetrates inside the boundary layer.

• 대한조선학회논문집
• /
• 제53권6호
• /
• pp.494-502
• /
• 2016

The accurate assessment of hull-appendage interaction in the early design stage is important to control the inflow to the propeller plane, which can cause undesirable hydrodynamic effects in terms of cavitation phenomenon. This paper describes a numerical analysis for the flow around a fully appended surface ship model for which KRISO has carried out a model test in the Large Cavitation Tunnel(LCT). This numerical study was performed with the LCT model test in a complementary manner for a good reproduction of the wake distribution of surface ships. A second order accurate finite volume method provided by a commercial computational fluid dynamics(CFD) program was used to solve the governing Reynolds Averaged Navier-Stokes(RANS) equations, where the SST $k-{\omega}$ model was used for turbulence closure. The numerical results were compared to available LCT experimental data for validation. The calculations gave good predictions for the boundary layer profiles on the walls of the empty cavitation tunnel and the wake at the propeller plane of the fully appended hull model in the LCT.

• Ocean Systems Engineering
• /
• 제12권2호
• /
• pp.247-265
• /
• 2022

In the present study, we focus on the CFD simulations for the performance and the rotor-generated wake of a model-scale wind turbine which was designed for wave tank experiments. The CFD simulations with fully resolved rotor geometry are performed using MARIN's community-based open-source CFD code ReFRESCO. The absolute formulation method (AFM) is leveraged to model the rotating wind turbine. The k - ω SST turbulence model is adopted in the incompressible Reynolds Averaged Navier-Stokes (RANS) simulations. First, the thrust and torque coefficients, C_T and C_P, are calculated at different Tip Speed Ratios (TSR), and the results are compared against the experimental data and previous numerical results. The pressure distribution of the turbine blades at the 70% span is obtained and compared to the results obtained by other tools. Then, a verification study aiming at quantifying the discretization uncertainty of the turbine performance with respect to the grid resolution in the wake region is performed. Last, the rotor-generated wake at the TSR of 7 is presented and discussed.

• Wind and Structures
• /
• 제17권5호
• /
• pp.479-494
• /
• 2013

The vehicle-induced aerodynamic loads bring vibrations to some of the highway sound barriers, for they are designed in consideration of natural wind loads only. As references to the previous field experiment, the vehicle-induced aerodynamic loads is investigated by numerical and theoretical methodologies. The numerical results are compared to the experimental one and proved to be available. By analyzing the flow field achieved in the numerical simulation, the potential flow is proved to be the main source of both head and wake impact, so the theoretical model is also validated. The results from the two methodologies show that the shorter vehicle length would produce larger negative pressure peak as the head impact and wake impact overlapping with each other, and together with the fast speed, it would lead to a wake without vortex shedding, which makes the potential hypothesis more accurate. It also proves the expectation in vehicle-induced aerodynamic loads on Highway Sound Barriers Part1: Field Experiment, that max/min pressure is proportional to the square of vehicle speed and inverse square of separation distance.

• Wind and Structures
• /
• 제35권1호
• /
• pp.35-54
• /
• 2022

A comparative study of aeroelastic vibrations of spring-mass cylinder and chimneys, with the help of a few wake oscillator models available in the literature, is presented. The models include those proposed by Facchinetti, Farshidian and Dolatabadi method-I, Farshidian and Dolatabadi method-II, de Langre, Skop and Griffin. Besides, the linear model proposed by Simiu and Scanlan is also incorporated in the study. For chimneys, the first mode oscillation is considered, and the top displacements of the chimneys are evaluated using the considered models. The results of the analytical model are compared with those obtained from the numerical solution of the wake-oscillator coupled equations. The response behavior of the cylinder and three chimneys of different heights are studied and compared with respect to critical parametric variations. The results of the study indicate that the numerical analysis is essential to capture the effect of non-linear aeroelastic phenomena in the solutions, especially for small damping. Further, except for the models proposed by Farshidian and Dolatabadi, other models predict nearly the same responses. The non-linear model predicts a much higher response as compared to the linear model.

• Wind and Structures
• /
• 제8권6호
• /
• pp.427-442
• /
• 2005

Numerical and wind tunnel simulations of pollutant dispersion around rectangular obstacles with five aspect ratios have been conducted in order to identify the effects of flow patterns induced by buildings on plume dispersion in the near wake of buildings. An emission from a low source located upwind of obstacles was used in this simulation. The local flow patterns and concentrations around a cubical obstacle were initially investigated using three RANS turbulence models, (the standard $k-{\varepsilon}$, Shear Stress Transport (SST), Reynolds-Stress RSM turbulence model) and also using Large-eddy simulation (LES). The computed concentrations were compared with those measured in the wind tunnel. Among the three turbulence models, the SST model offered the best performance and thus was used in further investigations. The results show, for normal aspect ratios of width to height, that concentrations in the near wake are appreciably affected because of plume capture by the horseshoe vortex and convection by the vertical vortex pairs. These effects are less important for high aspect ratios. Vertical vortex pairs present a strong ability to exchange mass vertically and acts efficiently to reduce ground-level concentrations in the near wake.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제7권3호
• /
• pp.595-609
• /
• 2015

The wake characteristics of Contra-Rotating Propeller (CRP) were investigated using numerical simulation and flow measurement. The numerical simulation was carried out with a commercial CFD code based on a Reynolds Averaged Navier-Stokes (RANS) equations solver, and the flow measurement was performed with Stereoscopic Particle Image Velocimetry (SPIV) system. The simulation results were validated through the comparison with the experiment results measured around the leading edge of rudder to investigate the effect of propeller operation under the conditions without propeller, with forward propeller alone, and with both forward and aft propellers. The evolution of CRP wake was analyzed through velocity and vorticity contours on three transverse planes and one longitudinal plane based on CFD results. The trajectories of propeller tip vortex core in the cases with and without aft propeller were also compared, and larger wake contraction with CRP was confirmed.

• 대한기계학회논문집B
• /
• 제22권6호
• /
• pp.757-770
• /
• 1998

A Navier-Stokes code with a low Reynolds number k-.epsilon. turbulence model was tested to investigate its predictability for the unsteady transitional boundary layer flow due to rotor-stator interaction. A preliminary calculation with three different numbers of time steps 300, 600, and 1000 for a rotor wake passing period was carried out to see the effects of time steps on the unsteady flow and pressure fields due to rotor-stator interaction. Numerical solutions showed that unsteady pressure was much more sensitive to the number of time steps and over 600 time steps should be used to get a numerical solution independent of the number of time steps for a rotor wake passing period. The original low Reynolds number k-.epsilon. turbulence model showed very poor prediction of the unsteady transitional boundary layer flow due to rotor-stator interaction. This was due to the excessive production of turbulent kinetic energy near the leading edge. A modification suggested by Launder was incorporated and the modified model captured well the wake induced transitional strip. Present solutions also showed improved prediction over previous Euler/boundary layer solution in terms of the onset of unsteady transition and its extent.