aepw-1 hirenasd test cases · ruag aviation, department aerodynamics operates two subsonic wind...

AePW-1 HIRENASD test cases

A. Gehri and D.Steiling

RUAG Aviation

Who we are

RUAG Aviation, Department Aerodynamics

Operates two subsonic wind tunnels

Large Wind Tunnel Emmen (LWTE), 7x5m, aerospace (powered/unpowered),

full scale automotive rain testing

Automotive Wind Tunnel Emmen (AWTE), 2.45x1.55m, with road simulation

Manufactures wind tunnel balances for other wind tunnels

CFD in collaboration with CFS Engineering at the EPFL (Swiss Federal

Institute of Technology) in Lausanne

Alain Gehri

Experienced CFD engineer, within AePW responsible for meshing and

setup of calculations

Daniel Steiling

Aerodynamic engineer, within AePW responsible for coordination and

post-processing

27.04.2012 RUAG Aviation 2

NSMB flow solver and settings

Multiblock Navier-Stokes solver, hence the name “NSMB”

Developed at the EPFL in Lausanne since 1991, together with other

universities and industrial partners

Settings used for the HIRENASD test cases:

Space discretization: 4th-order central scheme (Jameson)

Time integration: implicit LU-SGS scheme

Unsteady calculations: dual time stepping, w/ time correction procedure

Turbulence model: k-ω Menter Shear Stress (DES, RANS for static)

Particular version of NSMB had a bug in the ALE formulation

Dissipation for the turbulent equations was wrong, grid velocity not

included

Corrected now, and one case recalculated (not part of data submittal, but

shown in this presentation)

27.04.2012 RUAG Aviation 3

Test case specific settings and

assumptions

Cases calculated

ct311-4, ETW 159: static steady (c/m/f) -> static aeroelastic (m) ->

dynamic (m), each with forced transition (12% U, 5%L)

ct311-5, ETW 271: static steady (c/m/f) -> static aeroelastic (m) ->

dynamic (m), each fully turbulent

Dynamic calculations with forced oscillation of the second bending

mode

Modal amplitude of 0.01 reference chord used (wing tip), translates to

an amplitude of 2.04mm at position of accelerometer 15

Six periods have been simulated, with the last two periods used to

determine the FRF, 64 time steps per period

𝐹𝑅𝐹 =𝑓𝑓𝑡 𝑒𝑥𝑖𝑡𝑎𝑡𝑖𝑜𝑛 ⋅𝑓𝑓𝑡(𝑟𝑒𝑠𝑝𝑜𝑛𝑠𝑒)

𝑓𝑓𝑡 𝑒𝑥𝑡𝑎𝑡𝑖𝑜𝑛 ⋅𝑓𝑓𝑡(𝑒𝑥𝑡𝑎𝑡𝑖𝑜𝑛)

27.04.2012 RUAG Aviation 4

Grid overview

27.04.2012 RUAG Aviation 5

Grid convergence static steady

27.04.2012 RUAG Aviation 6

,0.0040

,0.0060

,0.0080

,0.0100

,0.0120

,0.0140

,0.0160

,0.0180

,0.0200

,0.3635

,0.3640

,0.3645

,0.3650

,0.3655

,0.3660

,0.3665

,0.3670

,0.3675

0.0E+00 1.0E-05 2.0E-05 3.0E-05 4.0E-05 5.0E-05

CD

CL

N-2/3

ct311-4 ETW159

CL

CD

,0.0129

,0.0131

,0.0133

,0.0135

,0.0137

,0.0139

,0.0141

,0.0330

,0.0335

,0.0340

,0.0345

,0.0350

,0.0355

,0.0360

0.0E+00 1.0E-05 2.0E-05 3.0E-05 4.0E-05 5.0E-05

CD

CL

N-2/3

ct311-5 ETW271

Global picture

27.04.2012 RUAG Aviation 7

ct311-4

ETW159

ct311-5

ETW271

Global picture (video)

ct311-5 ETW271

27.04.2012 RUAG Aviation 8

Global coefficients ct311-4 ETW159

27.04.2012 RUAG Aviation 9

Global coefficients ct311-5 ETW271

27.04.2012 RUAG Aviation 10

Global coefficients ct311-5 ETW271

ALE fixed

27.04.2012 RUAG Aviation 11

Cp at wing root

ct311-4 ETW159 section 1

27.04.2012 RUAG Aviation 12

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 10

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

ct311-4

ct311-4 ETW159

ct311-4 ETW159 aepw

ct311-4 ETW159 aepw conj

Cp mid wing

ct311-4 ETW159 section 3

27.04.2012 RUAG Aviation 13

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 10

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

ct311-4

ct311-4 ETW159

ct311-4 ETW159 aepw

ct311-4 ETW159 aepw conj

Cp time resolved

ct311-4 ETW159 section 3

27.04.2012 RUAG Aviation 14

Cp wing root

ct311-5 ETW271 section 1

27.04.2012 RUAG Aviation 15

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 10

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

ct311-5

ct311-5 ETW271

ct311-5 ETW271 aepw

ct311-5 ETW271 aepw conj

ct311-5-alefix

Cp mid wing

ct311-5 ETW271 section 3

27.04.2012 RUAG Aviation 16

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 10

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

ct311-5

ct311-5 ETW271

ct311-5 ETW271 aepw

ct311-5 ETW271 aepw conj

ct311-5-alefix

Cp time resolved

ct311-5 ETW271 section 3

27.04.2012 RUAG Aviation 17

Issues encountered & challenges

Flow at wing root, effect of fuselage dummy

Periodicity of unsteady flow for ETW271 in some areas questionable

after only six periods of excitation

Scaling/normalizing of results, plus complex conjugate issue

ALE bug encountered

27.04.2012 RUAG Aviation 18

27.04.2012 RUAG Division 19

Thank you for your attention!

Questions?

27.04.2012 RUAG Division 19

27.04.2012 RUAG Division 20

Backup

27.04.2012 RUAG Division 20

Unsteady pressure coefficients

ct311-4 ETW159 section 1

27.04.2012 RUAG Aviation 21

Unsteady pressure coefficients

ct311-4 ETW159 section 2

27.04.2012 RUAG Aviation 22

Unsteady pressure coefficients

ct311-4 ETW159 section 3

27.04.2012 RUAG Aviation 23

Unsteady pressure coefficients

ct311-4 ETW159 section 4

27.04.2012 RUAG Aviation 24

Unsteady pressure coefficients

ct311-4 ETW159 section 5

27.04.2012 RUAG Aviation 25

Unsteady pressure coefficients

ct311-4 ETW159 section 6

27.04.2012 RUAG Aviation 26

Unsteady pressure coefficients

ct311-4 ETW159 section 7

27.04.2012 RUAG Aviation 27

Unsteady pressure coefficients

ct311-5 ETW271 section 1

27.04.2012 RUAG Aviation 28

Unsteady pressure coefficients

ct311-5 ETW271 section 2

27.04.2012 RUAG Aviation 29

Unsteady pressure coefficients

ct311-5 ETW271 section 3

27.04.2012 RUAG Aviation 30

Unsteady pressure coefficients

ct311-5 ETW271 section 4

27.04.2012 RUAG Aviation 31

Unsteady pressure coefficients

ct311-5 ETW271 section 5

27.04.2012 RUAG Aviation 32

Unsteady pressure coefficients

ct311-5 ETW271 section 6

27.04.2012 RUAG Aviation 33

Unsteady pressure coefficients

ct311-5 ETW271 section 7

27.04.2012 RUAG Aviation 34