Frequency Domain Flutter Boundary Computations Using Navier-Stokes Equations on Superclusters

In this paper, an energy method for flutter analysis of wing using one-way fluid structure coupling was developed. To consider the effect of wing vibration, Reynolds-averaged Navier–Stokes equations based on the arbitrary Lagrangian Eulerian coordinates were employed to model the flow. The flow mesh was updated using a fast dynamic mesh technology proposed by our research group. The pressure was calculated by solving the Reynolds-averaged Navier–Stokes equations through the SIMPLE algorithm with the updated flow mesh. The aerodynamic force for the wing was computed using the pressure on the wing surface. Then the aerodynamic damping of the wing vibration was computed. Finally, the flutter stability for the wing was decided according to whether the aerodynamic damping was positive or not. Considering the first four modes, the aerodynamic damping for wing 445.6 was calculated using the present method. The results show that the aerodynamic damping of the first mode is lower than the aerodynamic damping of higher order modes. The aerodynamic damping increases with the increase of the mode order. The flutter boundary for wing 445.6 was computed using the aerodynamic damping of the first mode in this paper. The calculated flutter boundary is consistent well with the experimental data.

Download Full-text

Numerical Simulation of Transonic Fan Flutter With 3D N-S CFD Code

Volume 5: Structures and Dynamics, Parts A and B ◽

10.1115/gt2008-50573 ◽

2008 ◽

Cited By ~ 2

Author(s):

Mizuho Aotsuka ◽

Naoki Tsuchiya ◽

Yasuo Horiguchi ◽

Osamu Nozaki ◽

Kazuomi Yamamoto

Keyword(s):

Shock Wave ◽

Stokes Equations ◽

Navier Stokes ◽

Suction Surface ◽

Navier Stokes Equations ◽

Blade Loading ◽

Flow Solver ◽

Flutter Boundary ◽

Transonic Fan ◽

Volume Method

This paper describes the calculation of transonic stall flutter of a fan. A new CFD code has been developed and validated. The code is an unsteady 3D multi-block flow solver. The Reynolds-Averaged Navier-Stokes equations are solved using a finite volume method with Spallart-Allmaras 1 equation turbulence model. A grid deforming system is applied, so the new code is capable of simulating an oscillating blade row. This grid deforming system produces less grid distortion and the code has robustness for a blade oscillating calculation. The code has validated on an IHI’s research transonic fan rig test, and the result was in good agreement with the test data in the prediction of the flutter boundary. In the rig test at part-speed condition, stall-side flutter was experienced. In that condition, the inlet relative Mach number in the tip region is about unity. The aerodynamic work by the CFD at the near flutter condition is positive, which means that the flutter characteristic is unstable, while at other conditions the aerodynamic work is negative. The aerodynamic work increases rapidly just before the zero damping point with the increase of the blade loading. From the detailed CFD result, the shock wave on the suction surface contributes to the excitement of the blade oscillation, and the aerodynamic work of the shock wave has large value at the flutter condition.

Download Full-text