• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.19-19
• /
• 2009

Dynamic Stall is a flow phenomenon which occurs on the retreating side of helicopter rotor blades during forward flight. It also occurs on blades of stall regulated wind turbines under yawing conditions as well as during gust loads. Time scales occurring during this process are comparable on both helicopter and wind turbine blades. Dynamic Stall limits the speed of the helicopter and its manoeuvrability and limits the amount of power production of wind turbines. Extensive numerical as well as experimental investigations have been carried out recently to get detailed insight into the very complex flow structures of the Dynamic Stall process. Numerical codes have to be based on the full equations, i.e. the Navier-Stokes equations to cover the scope of the problems involved: Time dependent flow, unsteady flow separation, vortex development and shedding, compressibility effects, turbulence, transition and 3D-effects, etc. have to be taken into account. In addition to the numerical treatment of the Dynamic Stall problem suitable wind tunnel experiments are inevitable. Comparisons of experimental data with calculated results show us the state of the art and validity of the CFD-codes and the necessity to further improve calculation procedures. In the present paper the phenomenon of Dynamic Stall will be discussed first. This discussion is followed by comparisons of some recently obtained experimental and numerical results for an oscillating helicopter airfoil under Dynamic Stall conditions. From the knowledge base of the Dynamic Stall Problems, the next step can be envisaged: to control Dynamic Stall. The present discussion will address two different Dynamic Stall control methodologies: the Nose-Droop concept and the application of Leading Edge Vortex Generators (LEVoG's) as examples of active and passive control devices. It will be shown that experimental results are available but CFD-data are only of limited comparison. A lot of future work has to be done in CFD-code development to fill this gap. Here mainly 3D-effects as well as improvements of both turbulence and transition modelling are of major concern.

• Journal of computational fluids engineering
• /
• v.22 no.1
• /
• pp.59-66
• /
• 2017

The guideline of selecting the number of snapshot dataset, $N_s$ in proper orthogonal decomposition(POD) was presented via the analysis of Eigen values based on the singular value decomposition(SVD). In POD, snapshot datasets from the solutions of Euler or Navier-Stokes equations are utilized to SVD and a reduced order model(ROM) is constructed as the combination of Eigen vectors. The ROM is subsequently applied to reconstruct the flowfield data with new set of flow conditions, thereby enhancing the computational efficiency. The overall computational efficiency and accuracy of POD is dependent on the number of snapshot dataset; however, there is no reliable guideline of determining $N_s$. In order to resolve this problem, the order of maximum to minimum Eigen value ratio, O(R) from SVD was analyzed and presented for the decision of $N_s$; in case of steady flow, $N_s$ should be determined to make O(R) be $10^9$. For unsteady flow, $N_s$ should be increased to make O(R) be $10^{11\sim12}$. This strategy of selecting the snapshot dataset was applied to two dimensional NACA0012 airfoil and vortex flow problems including steady and unsteady cases and the numerical accuracies according to $N_s$ and O(R) were discussed.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.23 no.4
• /
• pp.28-34
• /
• 2019

In this study, pin-fin arrays, which are widely used for cooling turbine blades, were studied. The vortex generator in pin-fin arrays is located in front of the circular tube. The cross-section of the vortex generator is NACA-9410. The purpose of this study is to analyze heat transfer performance and flow characteristics of pin-fin arrays. The position of vortex generator is changed with the vertical flow direction on the bottom wall. Pin-fin arrays were calculated with 6000, 10000 and 15000 Reynolds number. The commercial program ANSYS v18.0 CFX and the turbulence model $k-{\omega}$ SST were used. As a result, the heat transfer performance increased up to 5.8% and pressure loss increased less than 1%.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.8
• /
• pp.809-815
• /
• 2009

The previous work(Cheong et al., 2006) where the characteristics of acoustic emissions of wind turbines has been investigated according to the methods of power regulation, has showed that the acoustic power of wind turbine using the stall control for power regulation is more correlated with the wind speed than that using the pitch control. In this paper, basically extending this work, the noise generation characteristics of large modern upwind wind turbines are experimentally indentified according to the power regulation methods. To investigate the noise generation mechanisms, the distribution of noise sources in the rotor plane is measured by using the beam-forming measurement system(B&K 7768, 7752, WA0890) consisting of 48 microphones. The array results for the 660 kW wind turbine show that all noise is produced during the downward movement of the blades. This result show good agreement with the theoretical result using the empirical formula with the parameters: the convective amplification; trailing edge noise directivity; flow-speed dependence. This agreement implies that the trailing edge noise is dominant over the whole frequency range of the noise from the 660 kW wind turbine using the pitch control for power regulation.

• Wind and Structures
• /
• v.25 no.5
• /
• pp.415-432
• /
• 2017

The nonlinear aerodynamic instability of a tensioned plane orthotropic membrane structure is theoretically investigated in this paper. The interaction governing equation of wind-structure coupling is established by the Von $K\acute{a}rm\acute{a}n's$ large amplitude theory and the D'Alembert's principle. The aerodynamic force is determined by the potential flow theory of fluid mechanics and the thin airfoil theory of aerodynamics. Then the interaction governing equation is transformed into a second order nonlinear differential equation with constant coefficients by the Bubnov-Galerkin method. The critical wind velocity is obtained by judging the stability of the second order nonlinear differential equation. From the analysis of examples, we can conclude that it's of great significance to consider the orthotropy and geometrical nonlinearity to prevent the aerodynamic instability of plane membrane structures; we should comprehensively consider the effects of various factors on the design of plane membrane structures; and the formula of critical wind velocity obtained in this paper provides a more accurate theoretical solution for the aerodynamic stability of the plane membrane structures than the previous studies.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.3
• /
• pp.27-36
• /
• 2000

A higher order panel method based on B-spline representation for both the geometry and the velocity potential is developed for the solution of the flow around two-dimensional lifting bodies. Unlike Lee/Kerwin, who placed multiple control points on each panel and solved the overdetermined system of equation by the least square approach, the present method places only as many number of control points as required by the unknowns of the problem. Especially, a null pressure jump Kutta condition at the trailing edge is found to be effective in stabilizing the solution process and in predicting the correct solution. The new approach, is validated to be accurate through comparison with the analytic solution for a 2-D airfoil and to be less time-consuming due to fewer number of panels required than that used in Lee/Kerwin.

• The KSFM Journal of Fluid Machinery
• /
• v.5 no.3 s.16
• /
• pp.40-45
• /
• 2002

Since the previous cut-and-try design algorithm requires much cost and time, the automated design technique with the CFD and optimum design algorithm has recently been concerned. In this work, the Navier-Stokes equation was solved to gain more detailed viscous flow information of cascade with rotor-stator interactions. The H-grid embedded by O-grid was generated to obtain more accurate solution by eliminating the branch cut of H-grid near airfoil surface. To handle the relative motion of the rotor to the stationary stator, the sliding multiblock method was applied and the cubic-spline interpolation was used on the block interface boundary. To validate present procedure, the time-averaged aerodynamic loads were compared with experimeatal data. A good agreement was obtained. The Modified Method of Feasible Direction (MMFD) was used to carry out the sensitivity analysis of the change of aerodynamic performance by the changes of the cascade geometry. The present optimization of the cascade gave a dramatic reduction of the drag while the lift maintains at the value within the user-specified tolerance.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.4
• /
• pp.18-27
• /
• 2002

Aerodynamic shape optimization of two-dimensional airfoils in inviscid compressible flows is performed using a continuous adjoint formulation on unstructured meshes. Accurate evaluation of the gradient is achieved by using a reconstruction scheme based on the Laplacian averaging. A least-square method with extended stencil is used for flow gradient calculations. Proper convergence criterion is studied on Euler and adjoint equations for efficient design. The present method has been applied to RAE2822 and NACA0012 airfoils such that wave drag can be minimized by removing the shock wave. An inverse design is also performed to recover the shock wave on the designed RAE2822 airfoil.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.8
• /
• pp.10-20
• /
• 2002

The present paper investigates methods to control dynamic stall using an optimal approach. An unsteady aerodynamic sensitivity analysis code is developed by a direct differentiation method from a two-dimensional unsteady compressible Navier-Stokes solver including a two-equation turbulence model. Dynamic stall control is conducted by minimizing an objective function defined at an instant instead of integrating for a period of time. Unsteady sensitivity derivatives of the objective function are calculated by the sensitivity code, and optimization is carried out using a linear line search method at every physical time step. Numerous examples of dynamic stall control using control parameters such as nose radius, maximum thickness of airfoil, or suction show satisfactory results.

• Advances in aircraft and spacecraft science
• /
• v.7 no.4
• /
• pp.309-333
• /
• 2020

The proposed study aims to numerically investigate the performance of hydrofoils in the context of amphibious aircraft application. In particular, we also study the effectiveness of a slotted hydrofoil in minimizing the cavitation phenomenon, to improve the overall water take-off performance of an amphibious aircraft. We use the ICON A5 as a base model for this study. First, we propose an approach to estimate the required hydrofoil surface area and to select the most suitable airfoil shape that can minimize cavitation, thus improving the hydrodynamic efficiency. Once the hydrofoil is selected, we perform 2D numerical studies of the hydrodynamic and cavitating characteristics of a non-slotted hydrofoil on ANSYS Fluent. In this work, we also propose to use a slotted hydrofoil to be a passive method to control the cavitation performance through the boundary layer control. Numerical results of several slotted configurations demonstrate notable improvement on the cavitation performance. We then perform a multiobjective optimization with a response surface model to simultaneously minimize the cavitation and maximize the hydrodynamic efficiency of the hydrofoil. The optimization takes the slot geometry, including the slot angle and lengths, as the design variables. In addition, a global sensitivity study has been carried and it shows that the slot widths are the more dominant factors.