• 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.

• Journal of Institute of Control, Robotics and Systems
• /
• v.18 no.2
• /
• pp.76-80
• /
• 2012

This paper demonstrates how to implement inherent pitching stability in an insect-mimicking flapping-wing system for vertical takeoff. Design and fabrication of the insect-mimicking flapping-wing system is briefly described focusing on the recent modification. Force produced by the flapping-wing systems is estimated using the UBET (Unsteady Blade Element Theory) developed in the previous work. The estimation shows that the wing twist placed in the modified system can improve thrust production for about 10 %. The estimated thrust is compared with the measured thrust, which proves that the UBET provides fairly good estimations for the thrust produced by the flapping-wing systems. The vertical takeoff test shows that inherent pitching stability can be implemented in an insect-mimicking flapping-wing system by aligning the aerodynamic force center and center of gravity.

• 한국신재생에너지학회:학술대회논문집
• /
• 2008.05a
• /
• pp.315-318
• /
• 2008

Numerical calculation for the 1MW class horizontal axis wind turbine blade has been carried out to estimate the magnitude between discrete noise and random noise. Farassat formula 1A was adopted to get the discrete noise signal, and blade element momentum theory was used to obtain the distribution of the aerodynamic data along the blade span. Fukano's approach was also adopted to calculate the unsteady aerodynamic random noise due to the Karman vortex generation at the trailing edge of the wind turbine blade. From the noise prediction for the 1MW class horizontal axis wind turbine, the frequency band of the discrete noise lies in the infrasound region, and that of the random noise lies in the audible band region.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.6-15
• /
• 2008

This paper deals with the aeroelastic instability of vibrating multiple blade rows under aerodynamic coupling with each other. A model composed of three blade rows, e.g., rotor-stator-rotor, where blades of the two rotor cascades are simultaneously vibrating, is considered. The displacement of a blade vibrating under aerodynamic force is expanded in a modal series with the natural mode shape functions, and the modal amplitudes are treated as the generalized coordinates. The generalized mass matrix and the generalized stiffness matrix are formulated on the basis of the finite element concept. The generalized aerodynamic force on a vibrating blade consists of the component induced by the motion of the blade itself and those induced not only by vibrations of other blades of the same cascade but also vibrations of blades in another cascade. To evaluate the aerodynamic forces, the unsteady lifting surface theory for the model of three blade rows is applied. The so-called k method is applied to determine the critical flutter conditions. A numerical study has been conducted. The flutter boundaries are compared with those for a single blade row. It is shown that the effect of the aerodynamic blade row coupling substantially modifies the critical flutter conditions.

• Structural Engineering and Mechanics
• /
• v.85 no.1
• /
• pp.105-118
• /
• 2023

The present work is an attempt to develop a simple and accurate finite element formulation for the assessment of thermal shock/thermally induced vibrations in pretwisted and tapered functionally graded material thin (FGM) blades obtained from Voigt and local representative volume elements (LRVE) homogenization models, based on neutral surface approach. The neutral surface of the FGM blade does not coincide with its mid-surface. A finite element model (FEM) is developed using first-order shear deformation theory (FSDT) and the FGM turbine blade is modelled according to the shallow shell theory. The top and the bottom layers of the FGM blade are made of pure ceramic and pure metal, respectively and temperature-dependent material properties are functionally graded in the thickness direction, the position of the neutral surface also depends on the temperature. The material properties are estimated according to two different homogenization models viz., Voigt or LRVE. The top layer of the FGM blade is subjected to high temperature and the bottom surface is either thermally insulated or kept at room temperature. The solution of the nonlinear profile of the temperature in the thickness direction is obtained from the Fourier law of heat conduction in the unsteady state. The results obtained from the present FEM are compared with the benchmark examples. Next, the effect of angle of twist, intensity of thermal shock, variable chord and span and volume fraction index on the transient response due to thermal shock obtained from the two homogenization models viz., Voigt and LRVE scheme is investigated. It is shown that there can be a significant difference in the transient response calculated by the two homogenization models for a particular set of material and geometric parameters.

• International Journal of Aeronautical and Space Sciences
• /
• v.15 no.4
• /
• pp.345-355
• /
• 2014

This paper focuses on aerodynamic and inertial modeling of the propeller for its applications in flight dynamics analyses of a propeller-driven airplane. Unsteady aerodynamic and inertial loads generated by the propeller are formulated using the blade element method, where the local velocity and acceleration vectors for each blade element are obtained from exact kinematic relations for general maneuvering conditions. Vortex theory is applied to obtain the flow velocities induced by the propeller wake, which are used in the computation of the aerodynamic forces and moments generated by the propeller and other aerodynamic surfaces. The vortex lattice method is adopted to obtain the induced velocity over the wing and empennage components and the related influence coefficients are computed, taking into account the propeller induced velocities by tracing the wake trajectory trailing from each of the propeller blades. Aerodynamic forces and moments of the fuselage and other aerodynamic surfaces are computed by using the wind tunnel database and applying strip theory to incorporate viscous flow effects. The propeller models proposed in this paper are applied to predict isolated propeller performances under steady flight conditions. Trimmed level forward and turn flights are analyzed to investigate the effects of the propeller on the flight characteristics of a propeller-driven light-sports airplane. Flight test results for a series of maneuvering flights using a scaled model are employed to run the flight dynamic analysis program for the proposed propeller models. The simulations are compared with the flight test results to validate the usefulness of the approach. The resultant good correlations between the two data sets shows the propeller models proposed in this paper can predict flight characteristics with good accuracy.

• The KSFM Journal of Fluid Machinery
• /
• v.18 no.2
• /
• pp.14-21
• /
• 2015

Design lifetime of a wind turbine is required to be at least 20 years. The most important step to ensure the deign is to evaluate the loads on the wind turbine as accurately as possible. In this study, extreme design load of a offshore wind turbine using Garrad Hassan (GH) Bladed and National Renewable Energy Laboratory (NREL) FAST codes are calculated considering structural dynamic loads. These wind turbine aeroelastic analysis codes are high efficiency for the rapid numerical analysis scheme. But, these codes are mainly based on the mathematical and semi-empirical theories such as unsteady blade element momentum (UBEM) theory, generalized dynamic wake (GDW), dynamic inflow model, dynamic stall model, and tower influence model. Thus, advanced CFD-dynamic coupling method is also applied to conduct cross verification with FAST and GH Bladed codes. If the unsteady characteristics of wind condition are strong, such as extreme design wind condition, it is possible to occur the error in analysis results. The NREL 5 MW offshore wind turbine model as a benchmark case is practically considered for the comparison of calculated designed loads. Computational analyses for typical design load conditions such as normal turbulence model (NTM), normal wind profile (NWP), extreme operation gust (EOG), and extreme direction change (EDC) have been conducted and those results are quantitatively compared with each other. It is importantly shown that there are somewhat differences as maximum amount of 18% among numerical tools depending on the design load cases.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.4
• /
• pp.297-305
• /
• 2011

In this study, the aeroelastic analysis of rotorcraft in forward flight has been performed using dynamic inflow model to handle unsteady aerodynamics. The quasi-steady airload model based on the blade element method has been coupled with dynamic inflow model developed by Peters and He. The nonlinear steady response to periodic motion is obtained by integrating the full finite element equation in time through a coupled trim procedure with a vehicle trim for stability analysis. The aerodynamic and structural characteristics of dynamic inflow model are validated against other numerical analysis results by comparing induced inflow and blade tip deflections(flap, lag). In order to validate aeroelastic stability of dynamic inflow model, lag damping are also compared with those of linear inflow model.

• International Journal of Aeronautical and Space Sciences
• /
• v.18 no.2
• /
• pp.186-196
• /
• 2017

In order to resolve the trimming difficulty in rotor CFD calculations, a high-efficiency and improved "delta trim method" is established to compute the blade control settings that are necessary to identify the blade motion. In this method, a simplified model which combines the blade element theory and different inflow models is employed to calculate the control settings according to the target aerodynamic forces, then it is coupled into a CFD solver with unsteady Navier-Stokes equations by the delta methodology, which makes the control settings and aerodynamics calculated and updated in the meantime at every trim cycle. Different from the previous work, the current research combines the inflow model based on prescribed wake theory. Using the method established, the control settings and aerodynamic characteristics of Helishape 7A, AH-1G and Caradonna-Tung rotors are calculated. The influence of different inflow models on trimming calculations is analyzed and the computational efficiency of the current "delta trim method" is compared with that of the "CFD-based trim method". Furthermore, for the sake of improving the calculation efficiency, a novel acceleration factor method is introduced to accelerate the trimming process. From the numerical cases, it is demonstrated that the current "delta trim method" has higher computational efficiency than "CFD-based trim method" in both hover and forward flight, and up to 70% of the amount of calculation can be saved by current "delta trim method" which turns out to be satisfactory for engineering applications. In addition, the proposed acceleration factor shows a good ability to accelerate the trim procedure, and the prescribed wake inflow model is always of better stability than other simple inflow models whether the acceleration factor is utilized in trimming calculations.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.46 no.5
• /
• pp.368-375
• /
• 2018

In the present study, aerodynamic load analysis for a floating off-shore wind turbine was conducted to examine the effect of periodic platform motion in the direction of 6-DOF on rotor aerodynamic performance. Blade-element momentum method(BEM) was used for a numerical simulation, the unsteady airload effects due to the flow separation and the shed wake were considered by adopting a dynamic stall model based on the indicial response method. Rotor induced downwash was estimated using the momentum theory, coupled with empirical corrections for the turbulent wake states. The periodic platform motions including the translational motion in the heave, sway and surge directions and the rotational motion in the roll, pitch and yaw directions were considered, and each platform motion was applied as a sinusoidal function. For the numerical simulation, NREL 5MW reference wind turbine was used as the target wind turbine. The results showed that among the translation modes, the surge motion has the largest influence on changing the rotor airloads, while the effect of pitch motion is predominant for the rotations.