Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow...

34
© Copyright Cenaero 2009 – All rights reserved EXCELLENCE IN SIMULATION TECHNOLOGIES Development and Application of Partioned Methods for Large Scale FSI Problems Philippe Geuzaine CFD-MP Group, Cenaero Email: [email protected]

Transcript of Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow...

Page 1: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

© Copyright Cenaero 2009 – All rights reserved

EXCELLENCE IN SIMULATION TECHNOLOGIES

Development and Application of Partioned Methods for Large Scale FSI Problems

Philippe GeuzaineCFD-MP Group, Cenaero

Email: [email protected]

Page 2: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Sample Targeted Applications

Page 3: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Unsteady Aeroelastic Response

Page 4: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Four-field Formulation

Heat Transfer Fluid (ALE)

Structure Fluid Mesh

Temperature(Displacements) Force Mesh Motion

TemperatureHeat flux

DisplacementsVelocities

Page 5: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Aero- & Elasto-Dynamic Data Exchange

• Interpolation across non matching interfaces (accuracy and conservation)

Page 6: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Exchange With Non-Matching Interfaces

Page 7: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Example of Displacement Exchange

Page 8: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Data Exchange & Code Coupling

• Dynamic data communication between CSD and CFD modules

• Available tools (integration platforms) for interpolation and/or data communication– MpCCI (Fraunhofer Institute SCAI)– MDICE (CFD RC)– SimServer (EADS)– PALM (CERFACS)– CALCIUM (EDF)– SALOME (EDF/CEA)– ...

Page 9: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Fluid Mesh Motion

• Spring-type analogy• No node collapse• No face penetration• Mesh recovery

(after one oscillation cycle)• Mesh shearing?

Page 10: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

CFD on Moving Grids

• Navier-Stokes equations on fixed grids

• Arbitrary Lagrangian Eulerian (ALE) formulation on moving grids

• The discrete geometric conservation law (DGCL) can guide the selection among the multiple accuracy preserving ALE extensions

Page 11: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Additional Issues for Nonlinear FSI

• Non specific issues (CAD to mesh generation, turbulence modeling, nonlinear structural effects modeling, ...)

• Modularity & compatibility with existing tools (plug-and-play/off-the-shelf approach)

• Speed & turnaround time– Second-order accurate & robust coupled time-

integration for large time-steps– Parallel computing

• Verification & validation

Page 12: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Monolithic Schemes for Coupled Systems

• Single time-integrator (system of first-order semi-discrete equations)

• Solution at the nonlinear level (single set of equations for fluid, fluid mesh and structure)– Newton type method

• Solution at the linear level (very large system)– Single approach (direct, iterative)– Partioned or staggered approach (block Jacobi, block

Gauss-Seidel)

• Appealing for accuracy and stability, but issues for software modularity

• Demonstrated usually for academic problems

Page 13: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Partioned Schemes for Coupled Systems

• Collocated staggered approach

• Loosely-coupled scheme (maximize efficiency)• Strongly-coupled scheme (similar monolithic)

CSD

CFD

CFD Mesh

Page 14: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Individual Time-integrator Selection

• Second-order accurate time-integrators (necessary condition to achieve second-order time-accuracy for the coupled system)

• Newmark scheme (with the equilibrium written at time level tn+1) for the CSD solver (no mass or stiffness matrices exported from the CSD solver)

• ALE version of the three-point backward difference scheme for the CFD solver (not trivially second-order time-accurate on moving meshes)

Page 15: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Sample Predictors and Correctors

• Structural displacement predictor

• Aerodynamic force corrector

Page 16: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Control Parameters Design

• Control parameters– Structural displacement predictor– Mesh integrator– Aerodynamic force corrector

• Preserve second-order time-accuracy of individual time-integrators (three-point backward difference and Newmark schemes)

Page 17: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

AGARD 445.6 Wing Configuration

• Nonlinear flow model (inviscid Argo model)• Linear structural model (SAMCEF Mecano model)• Steady & unsteady aeroelastic analyses• Flutter analysis

Page 18: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

CFD Meshes

22,014 points

178,938 points

Page 19: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Steady Simulation

Page 20: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Unsteady Simulations

Page 21: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Time-Accuracy Analysis

Page 22: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Unsteady Simulations with P/C

Page 23: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Unsteady Simulations with Large Δt

Page 24: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

AGARD 445.6 Wing Flutter Analysis

• Flow conditions: 0.5 < M < 1.141 (zero incidence)• Density taken from Yates’ experiment• Pressure adjusted to produce flutter• Coupling time-step Δt = 1 ms• Flutter index definition

Page 25: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Response at Mach = 0.957

Page 26: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Flutter Frequency

Page 27: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Flutter Index

Page 28: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

CAD-based Wing MDO

CAD : CATIAV5

AutomaticMesh Generation

Parallel Aeroelastic Computation

Argo + Samcef

Page 29: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Optimization Workflow with Minamo

UserUserSpecifications

Approximate ModelANN, RBF, Kriging, …

OptimizationOptimizationEA, gradient-based, …

Performance Check

DATABASEDATABASE

Accurate Model Accurate Model CFD / Structure / Exp. / ...

END

ONLINE modeling

Page 30: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Simulation Workflow with Argo/Samcef

Design Variables

Inviscid Aeroelastic Computation

CFD CSD

Wing Deformed Shape

Viscous Computation

CFD

Lift & dragStress, buckling,tip displacement

CAD generation& meshing

Meshing with Boundary Layers

Page 31: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Wing Parameterization

Tip chord

Root chord

Span

Sweep angle

twist

Angle ofattack

3 spars

10 ribs

Stringers 1 Stringer

s 2 Stringers 3 Stringer

s 4 Stringers 5

Panels 1 Panels

2Panels

3

Page 32: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Optimized Stringers and Planform

• Objective: wing mass reduction• Constraints: L/D, stress, buckling

Page 33: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Future Work on Supersonic Transport

Page 34: Development and Application of Partioned Methods for … 445.6 Wing Flutter Analysis • Flow conditions: 0.5 < M < 1.141 (zero incidence) • Density taken from Yates’ experiment

GDR IFS – Apr 23, 2009 © Copyright Cenaero 2009

Aerothermoelasticity Challenge at Mach 5