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

• 제어로봇시스템학회논문지
• /
• 제18권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년도 춘계학술대회 논문집
• /
• 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.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년 영문 학술대회
• /
• 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
• /
• 제85권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
• /
• 제15권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.

• 한국유체기계학회 논문집
• /
• 제18권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.

• 한국항공우주학회지
• /
• 제39권4호
• /
• pp.297-305
• /
• 2011

본 논문에서는 비정상 공기력을 모델링하기위해 Peters-He의 동적유입류모델을 고려한 2차원 준정상 공기력 이론을 적용하여 회전익기 전진비행에 대한 공탄성 해석을 수행하였다. 또한, 공력탄성학적 안정성 해석을 수행하기 위하여, 전진비행 시 주기적인 특성을 갖는 비선형 정적 트림 해를 얻기 위해 동체 평형을 고려한 연계 트림 해석을 통한 완전 유한요소 방정식을 이용하였다. 동적유입류모델의 공력과 구조 특성을 검증하기 위해 유도 유입류와 깃끝에서의 구조변형을 타 수치해석결과와 비교하였다. 또한, 공탄성 안정성을 검증하기 위해 두 모델의 래그 감쇠값을 비교하였다.

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

• 한국항공우주학회지
• /
• 제46권5호
• /
• pp.368-375
• /
• 2018

본 연구에서는 부유식 플랫폼의 6자유도 방향으로의 주기 운동이 로터 공력 성능에 미치는 영향을 확인하기 위해 부유식 해상 풍력터빈에 대한 공력 해석이 수행되었다. 수치 해석을 위해 블레이드 요소 운동량 방법을 이용하였으며, 유동 박리와 후류 영향에 의한 비정상 공력 효과를 포착하기 위해 인디셜 응답 방법에 기반한 동적 실속 모델을 이용하였다. 로터에 의해 유도되는 내리 흐름은 운동량 이론과 난류 후류 상태에 대한 경험적 모델을 연계하여 계산하였다. heave, sway, surge 방향으로의 병진 운동과 roll, pitch, yaw 방향으로의 회전 운동을 포함한 플랫폼 주기 운동을 고려하였으며, 각각의 모션은 사인함수 형태로 적용되었다. 수치해석을 위한 대상 풍력터빈으로는 NREL 5MW 풍력터빈이 사용되었다. 해석 결과로부터 세 방향 병진 운동 모드 중, surge 운동 시 로터 공력 변화가 상대적으로 크게 나타났으며, 회전 운동 모드의 경우, pitch 운동에 의해 로터 공력이 크게 변화됨을 확인할 수 있었다.