• International Journal of Aeronautical and Space Sciences
• /
• v.2 no.2
• /
• pp.28-38
• /
• 2001

A numerical method for the assessment and correction of tunnel wall interference effects on forced-oscillation testing is presented. The method is based on the wall pressure signature method using computed wall pressure distributions. The wall pressure field is computed using unsteady three-dimensional full Navier-Stokes solver for a 70-degree pitching delta wing in a wind tunnel. Approximately-factorized alternate direction implicit (AF-ADI) scheme is advanced in time by solving block tri-diagonal matrices. The algebraic Baldwin-Lomax turbulence, model is included to simulate the turbulent flow effect. Also, dual time sub-iteration with, local, time stepping is implemented to improve the convergence. The computed wall pressure field is then imposed as boundary conditions for Euler re-simulation to obtain the interference flow field. The static computation shows good agreement with experiments. The dynamic computation demonstrates reasonable physical phenomena with a good convergence history. The effects of the tunnel wall in upwash and blockage are analyzed using the computed interference flow field for several reduced frequencies and amplitudes. The corrected results by pressure signature method agree well with the results of free air conditions.

• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.1
• /
• pp.30-45
• /
• 2013

The present work focuses on the unsteady aerodynamics and aeroelastic properties of a small-medium sized wind-turbine blade operating under ideal conditions. A tapered/twisted blade representative of commercial blades used in an experiment setup at the National Renewable Energy Laboratory is considered. The aerodynamic loads are computed using Computational Fluid Dynamics (CFD) techniques. For this purpose, FLUENT$^{(R)}$, a commercial finite-volume code that solves the Navier-Stokes and the Reynolds-Averaged Navier-Stokes (RANS) equations, is used. Turbulence effects in the 2D simulations are modeled using the Wilcox k-w model for validation of the CFD approach. For the 3D aerodynamic simulations, in a first approximation, and considering that the intent is to present a methodology and workflow philosophy more than highly accurate turbulent simulations, the unsteady laminar Navier-Stokes equations were used to determine the unsteady loads acting on the blades. Five different blade pitch angles were considered and their aerodynamic performance compared. The structural dynamics of the flexible wind-turbine blade undergoing significant elastic displacements has been described by a nonlinear flap-lag-torsion slender-beam differential model. The aerodynamic quasi-steady forcing terms needed for the aeroelastic governing equations have been predicted through a strip-theory based on a simple 2D model, and the pertinent aerodynamic coefficients and the distribution over the blade span of the induced velocity derived using CFD. The resulting unsteady hub loads are achieved by a first space integration of the aeroelastic equations by applying the Galerkin's approach and by a time integration using a harmonic balance scheme. Comparison among two- and three- dimensional computations for the unsteady aerodynamic load, the flap, lag and torsional deflections, forces and moments are presented in the paper. Results, discussions and pertinent conclusions are outlined.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.8
• /
• pp.91-100
• /
• 2004

A Green's function approach is adopted for analyzing the thermoelastic deformations and stresses of a plate made of functionally graded materials(FGMs). The solution to the 3-dimensional unsteady temperature is obtained by using the laminate theory. The fundamental equations for thermoelastic problems are derived in terms of out-plane deformation and in-plane force, separately. The thermoelastic deformation and the stress distributions due to the bending and in-plane forces are analyzed by using a Green's function based on the Galerkin method. The eigenfunctions of the Galerkin Green's function for the thermoelastic deformation and the stress distributions are approximated in terms of a series of admissible functions that satisfy the homogeneous boundary conditions of the rectangular plate. Numerical analysis for a simply supported plate is carried out and effects of material properties on unsteady thermoclastic behaviors are discussed.

• 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.39 no.10
• /
• pp.919-926
• /
• 2011

The fin of high speed air vehicle is composed of skins and strong skeletons. In the flutter analysis, the eigenmodes of a fin are used for evaluating the unsteady aerodynamic force and the modal approach is applied for solving the flutter equation in both time and frequency domain. Therefore, the proper eigenmodes used for a modal flutter analysis should be chosen. For the appropriate choice of eigenmodes, when there exist local modes of a skin in the high modes, the effects of those modes on the unsteady aerodynamic force and flutter characteristics are anlalyzed.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.12
• /
• pp.1097-1106
• /
• 2011

Unsteady numerical analysis was performed to simulate air-launched missiles from a complete helicopter in hover by using an unstructured overset mesh flow solver coupled with a module of six degree-of-freedom motion of equations. The unsteady computations have been performed to obtain flow fields around the complete helicopter including main rotor, tail rotor, and fuselage equipped with multiple missiles, and six-DOF simulation has been performed to predict the behavior of the air-launched missile. The effects of the launching position and the missile thrust on the trajectory of the missile were investigated as well as the aerodynamic interference of the air-launched missile under the unsteady downwash produced by main rotor.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.294-298
• /
• 2003

The complex, unsteady, self-sustained pressure oscillations induced by supersonic flow past a rectangular cavity is investigated using numerical simulations. The present numerical study is performed using a parallel, multiblock solver for the two-dimensional, compressible Navier­Stokes equations. Open cavities with length-to-depth (L / D) ratio in the range 0.5 - 3.3 are considered. This paper sheds light on the cavity physics, cavity oscillatory mechanism, and the organisation of vortical structures inside the cavity. The vortex shedding phenomenon, the shear layer impingement event at the aft wall and the movement of the acoustic/compression wave within the cavity are well predicted. The vortical structures· and the source of the acoustic disturbances are found to be located near the aft wall of the cavity. With the increase in the cavity length, strong recompression of the flow near the aft wall leading to a sudden jump in the cavity form drag is observed. The estimated cavity tones are in good agreement with the available semi­empirical relation. Multiple peaks are noticed in deep and long cavities. For the present free­stream Mach number 1.71, it is observed that around L/D=2.0, the cavity oscillatory mechanism changes from the transverse to longitudinal oscillatory mode. The effects of this transition on various fluid dynamics and acoustic properties are also discussed.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.12 no.1
• /
• pp.42-47
• /
• 2002

Sirocco fans are widely used in HVAC and air conditioning systems, and the noise generated by these machines causes one of the most serious problems. In general, the sirocco fan noise is often dominated by tones at BPF(blade passage frequency) and broadband noise. However, only a few researches have been carried out on predicting the aeroacoustic noise because of the difficulty in obtaining detailed information about the flow field and casing effects on noise radiation. The objective of this study is to develop a prediction method for the unsteady flow field and the acoustic pressure field of a sirocco fan. We assume that the impeller rotates with a constant angular velocity and the flow field around the impeller is incompressible and inviscid. So, a discrete vortex method (DVM) is used to model the centrifugal fan and to calculate the flow field. The force of each element on the blade is calculated by the unsteady Bernoulli equation. Lowson\`s method is used to predict the acoustic source. Reasonable results are obtained not only fur the tonal noise but also far the amplitudes of the broadband noise. Acoustic pressure is proportional to (Ω)2.3, which is the similar value with the measured data.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.4
• /
• pp.431-438
• /
• 2009

A Large Eddy Simulation (LES) of the flow in an inducer is carried out under flow rate oscillations. The present study focuses on the dynamic response of the backflow and the unsteady pressure performance to the flow rate fluctuations under non-cavitation conditions. The amplitude of angular momentum fluctuation evaluated by LES is larger than that evaluated by RANS. However, the phase delay of backflow is nearly the same as RANS calculation. The pressure performance curve exhibits a closed curve caused by the inertia effect associated with the flow rate fluctuations. Compared with simplified one dimensional evaluation of the inertia component, the component obtained by LES is smaller. The negative slope of averaged performance curve becomes larger under unsteady conditions. From the conservations of angular momentum and energy, an expression useful for the evaluation of unsteady pressure rise was obtained. The examination of each term of this expression show that the apparent decrease of inertia effects is caused by the response delay of Euler's head and that the increase of negative slope is caused by the delay of inertial term associated with the delay of backflow response. These results are qualitatively confirmed by experiments.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.3
• /
• pp.14-21
• /
• 2006

The unsteady vortex lattice method is used to model twisting and flapping motions of a rectangular flat plate wing. The results for plunging and pitching motions were compared with the limited experimental results available and other numerical methods. They show that the method is capable of simulating many of the features of complex flapping flight. The lift, thrust and propulsive efficiency of a rectangular flat plate wing have been calculated for various twisting angles and reduced frequency with an amplitude of flapping angle($20^{\circ}$). And the effects of the twisting on the aerodynamic characteristics of the flapping wing are discussed by examination of their trends.