• International Journal of Aeronautical and Space Sciences
• /
• v.10 no.2
• /
• pp.125-133
• /
• 2009

To numerically simulate aerodynamics of rotor-airframe interaction in a rigorous manner, we need to solve the Navier-Stokes system for a rotor-airframe combination as a whole. This often imposes a serious computational burden since rotating blades and a stationary body have to be simultaneously dealt with. An efficient alternative is to adopt a momentum source method in which the action of rotor is approximated as momentum source over a rotor disc plane in a stationary computational domain. This makes the simulation much simpler. For unsteady simulation, the instantaneous momentum sources are assigned only to a portion of disk plane corresponding to blade passage. The momentum source is obtained by using blade element theory with dynamic inflow model. Computations are carried out for the simple rotor-airframe model (the Georgia Tech model) and the results of the simulation are compared with those of the full Navier-Stokes simulation with moving mesh system for rotor and with experimental data. It is shown that the present simulation yields results as good as those of the full Navier-Stokes simulation.

• 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 Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.2
• /
• pp.11-21
• /
• 2005

A three-dimensional parallel Euler flow solver has been developed for the simulation of unsteady rotor-fuselage interaction aerodynamics on unstructured meshes. In order to handle the relative motion between the rotor and the fuselage, the flow field was divided into two zones, a moving zone rotating with the blades and a stationary zone containing the fuselage. A sliding mesh algorithm was developed for the convection of the flow variables across the cutting boundary between the two zones. A quasi-unsteady mesh adaptation technique was adopted to enhance the spatial accuracy of the solution and to better resolve the wake. A low Mach number pre-conditioning method was implemented to relieve the numerical difficulty associated with the low-speed forward flight. Validations were made by simulating the flows around the Georgia Tech configuration and the ROBIN fuselage. It was shown that the present method is efficient and robust for the prediction of complicated unsteady rotor-fuselage aerodynamic interaction phenomena.

• Journal of Mechanical Science and Technology
• /
• v.20 no.7
• /
• pp.1019-1029
• /
• 2006

The ballistic range has long been employed in a variety of engineering fields such as high-velocity impact engineering, projectile aerodynamics, creation of new materials, etc, since it can create an extremely high-pressure state in very short time. Of many different types of ballistic ranges developed to date, two-stage light gas gun is being employed most extensively. In the present study, a theoretical work has been made to develop a new type of ballistic range which can easily simulate a flying projectile. The present ballistic range consists of high-pressure tube, piston, pump tube, shock tube and launch tube. The effect of adding a shock tube in between the pump tube and launch tube is investigated. This improvement is identified as the reduction in pressures in the high pressure tube and pump tube while maintaining the projectile velocity. Equations of motions of piston and projectile are solved using Runge-Kutta methods. Dependence of projectile velocity on various design factors such as high pressure tube pressure, piston mass, projectile mass, area ratio of pump tube to launch tube and type of driver gas in the pump tube are also analyzed. Effect of various gas combinations is also investigated. Calculations show that projectile velocities of the order 8 km/sec could be achieved with the present ballistic range.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.2
• /
• pp.98-107
• /
• 2002

An experimental and numerical study of aerodynamic charateristics of KSR-III and jet impingement on a deflector during launch has been conducted. To investigate aerodynamic characteristics of KSR-III configuration, wind tunnel tests using 6.4% scale model were performed by 4x4 feet ADD trisonic wind tunnel on the Mach number range of 0.4~3.8. Solutions of Three dimensional Euler equations were also obtained and compared with test result. For the study of KSR-III jet impingement flowfield on a deflector during launch operation, unsteady computation using CFD-FASTRAN was performed.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2005.11a
• /
• pp.86-89
• /
• 2005

A computational work has been performed to investigate the aerodynamics of a projectile which is launched from the two-stage light gas gun. A moving coordinate method for a multi-domain technique is employed to simulate unsteady projectile flows with a moving boundary. The effect of a virtual mass is added to the axisymmetric unsteady Euler equation system. The computed results reasonably capture the major flow characteristics which we generated in launching the projectile supersonically, such as the interaction between the shock wave and the blast wave, the interaction between the vortical flow and the barrel shock, and the steady under-expanded jet. The present computational results properly predict the velocity, acceleration, and drag histories of the projectile.

• Wind and Structures
• /
• v.12 no.1
• /
• pp.1-19
• /
• 2009

In this paper, the airflow around an ideal thin plate (hereafter referred to as ITP) with various ratios of central slot is simulated by using the finite-difference-method (FDM)-based Arbitrary-Lagrangian-Eulerian descriptions for the rigid oscillating body. The numerical procedure employs the second-order projection scheme to decouple the governing equations, and the multigrid algorithm with three levels to improve the computational efficiency in evaluating of the pressure equation. The present CFD method is validated through comparing the computed flutter derivatives of the ITP without slot to Theodorsen analytical solutions. Then, the unsteady aerodynamics of the ITP with and without central slot is investigated. It is found that even a smaller ratio of central slot of the ITP has notable effects on pressure distributions of the downstream section, and the pressure distributions on the downstream section will further be significantly affected by the slot ratio and the reduced wind speeds. Continuous increase of $A_2^*$ with the increase of central slot may be the key feature of the slotted ITP. Finally, flutter analyses based on the flutter derivatives of the slotted ITP are performed, and moreover, flutter instabilities of a scaled sectional model of a twin-deck bridge with various ratios of deck slot are investigated. The results confirm that the central slot is effective to improve bridge flutter stabilities, and that the flutter critical wind speeds increase with the increase of slot ratio.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.31-35
• /
• 2006

The delta-formulation of the Navier-Stokes equations has been popularly used in the aerodynamics area. Implicit algorithm can be easily implemented in that by using Taylor series expansion. This formulation is extended for an unsteady analysis by using a dual-time integration. In the meanwhile, the incompressible flows with heat transfers which occur in the area of thermo-hydraulics have been solved by a segregated algorithm such as the SIMPLE method, where each equation is discretised by using an under-relaxed deferred correction method and solved sequentially. In this study, the dual-time delta formulation is implemented in the segregated Navier-Stokes solver which is based on the collocated cell-centerd scheme with un unstructured mesh FVM. The pressure correction equation is derived by the SIMPLE method. From this study, it was found that the Euler dual-time method in the delta formulation can be combined with the SIMPLE method.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.4
• /
• pp.571-580
• /
• 2015

In this paper, a computational study was conducted in order to investigate the rotor blade sweep effect on the aerodynamics of a small axial supersonic impulse turbine stage. For this purpose, three-dimensional unsteady RANS simulations have been performed with three different rotor blade sweep angles ($-15^{\circ}$, $0^{\circ}$, $+15^{\circ}$) and the results were compared with each other. Both NTG (No tip gap) and WTG (With tip gap) models were applied to examine the effect on tip leakage flow. As a result of the simulation, the positive sweep model ($+15^{\circ}$) showed better performance in relative flow angle, Mach number distribution, entropy rise, and tip leakage mass flow rate compared with no sweep model. With the blade static pressure distribution result, the positive sweep model showed that hub and tip loading was increased and midspan loading was reduced compared with no sweep model while the negative sweep model ($-15^{\circ}$) showed the opposite result. The positive sweep model also showed a good aerodynamic performance around the hub region compared with other models. Overall, the positive sweep angle enhanced the turbine efficiency.

• Wind and Structures
• /
• v.31 no.6
• /
• pp.561-573
• /
• 2020

Unsteady self-excited forces are commonly represented by parametric models such as rational functions. However, this requires complex multiparametric nonlinear fitting, which can be a challenging task that requires know-how. This paper explores the alternative nonparametric modeling of unsteady self-excited forces based on relations between flutter derivatives. By exploiting the properties of the transfer function of linear causal systems, we show that damping and stiffness aerodynamic derivatives are related by the Hilbert transform. This property is utilized to develop exact simplified expressions, where it is only necessary to consider the frequency dependency of either the aeroelastic damping or stiffness terms but not both simultaneously. This approach is useful if the experimental data on aerodynamic derivatives that are related to the damping are deemed more accurate than the data that are related to the stiffness or vice versa. The proposed numerical models are evaluated with numerical examples and with data from wind tunnel experiments. The presented method can evaluate any continuous fitted table of interpolation functions of various types, which are independently fitted to aeroelastic damping and stiffness terms. The results demonstrate that the proposed methodology performs well. The relations between the flutter derivatives can be used to enhance the understanding of experimental modeling of aerodynamic self-excited forces for bridge decks.