Code coupling for simulation of flowinduced vibrations Elisabeth

Code coupling for simulation of flow-induced vibrations Elisabeth LONGATTE Fabien HUVELIN Mhamed SOULI 1 2005 EDF R&D

FRAMEWORK ASTER COMPUTATION SATURNE COMPUTATION LOADING DISPLACEMENT 2 2005 END

FRAMEWORK ASTER COMPUTATION SATURNE LOADING COMPUTATION LOADING DISPLACEMENT 3 2005 END

FRAMEWORK ASTER SUPERVISOR COSMETHYC SATURNE LOADING DISPLACEMENT COMPUTATION LOADING COMPUTATION DISPLACEMENT Convergence Test Conditions on time step Conditions on sub-cycling IF SUBCYCLING IF SUB CYCLING LOADING DISPLACEMENT No 4 2005 T=T + DT Si T > Tfinal END No

CODE_SATURNE Code_Saturne CFD Code developed by EDF R&D • Two- and three-dimensional calculations of steady or transient single-phase, incompressible, laminar or turbulent flows · Finite volume approach · Fully co-located arrangement of all variables · Time discretization based on a predictorcorrector scheme · Any kind of mesh (hybrid, any type of cell) · RANS model, LES · ALE formulation (moving boundary) 5 2005

CODE_ASTER Code_Aster CSD Code developed by EDF R&D • Linear, non linear statics • Linear, non linear dynamics · Finite element method · · 6 2005 Time calculation (Newmark…) Modal calculation (Newmark, Euler…) Stochastic approaches Dynamic response under loading

COSMETHYC Cosmethyc CFSD Code developed by EDF R&D • Fluid loading • Structure velocity · Iterative procedure 7 2005

COUPLING PROCESS Initialization Statics computation Fluid solver Coupling Fluid solver Structure solver Transient computation Time step loop Fluid solver 8 2005 Coupling End Structure solver

DATA TRANSFER Data transfer operators Interface • Inlet, outlet data interpolation • Data projection · 2 D, 3 D / 1 D beam, 2 D, 3 D 9 2005

INTERPOLATION Mesh interface Aster mesh Saturne mesh Non-matching interface 10 2005

SUPERVISOR Supervisor Coupling scheme • Governs time iteration • Governs coupling scheme · Prediction – correction on loading (convergence test) · Sub-cycling (implicit) 11 2005

COUPLING SCHEMES Explicit synchronous scheme time step Prediction of the fluid mesh motion Boundary conditions on the fluid-structure interface : Fluid solver force computation (Fn) Structure solver displacement computation (Un+1 ) Explicit staggered scheme time step Prediction of the fluid mesh motion Boundary conditions on the fluid-structure interface : Fluid solver force computation (Fn+1/2) Structure solver displacement computation (Un+1 ) (Farhat et al. , 1995, 1997; Piperno et al. , 2005 12 1995, 1997)

COUPLING SCHEMES Implicit scheme Initialization : Sub-cycling time step Fluid solver force computation : Structure solver displacement computation : Fluid mesh motion : Convergence criterion on the force values : 13 2005

COUPLING SCHEMES 1 D test case f(-) x (-) Explicit synchronous 0, 0 8, 0 10 -4 Explicit asynchronous 0, 0 7, 0 10 -6 Implicit 0, 0 9, 0 10 – 12 Strong coupling 0, 0 7, 0 10 -12 Analytical 0, 0 Error calculation / theory (-) 14 2005

VALIDATIONS Concentric tubes 15 2005

VALIDATIONS Eccentric tubes 16 2005

VALIDATIONS Tube bundles Frequency (Hz) Damping (Hz) 17 2005 Experimental Analytical Numerical - 20. 3 20. 5 0. 037 ± 0. 004 0. 037 0. 036

VALIDATIONS Tube bundles V<Vc V~Vc 18 2005 V>Vc

APPLICATIONS 19 2005

APPLICATIONS 20 2005

APPLICATIONS 21 2005

PERSPECTIVES PARALLEL Parallel code for distributed memory machines 22 2005

PERSPECTIVES PARALLEL Domain partitioning ALE • Nodes Coupling • Data transfer operators Wi 23 2005

PERSPECTIVES CONTACT 24 2005

PERSPECTIVES CONTACT 25 2005

PERSPECTIVES CONTACT 26 2005

PERSPECTIVES CONTACT 27 2005

PERSPECTIVES SALOME SUPERVISOR : To build schemes et control calculations – (with graphics) Node = component (fluid, mechanics, thermics) Port = inlet / outlet parameters to and from nodes Link = connexion between ports Checking data types 28 2005