Cast3M www-cast3m.cea.fr

Size: px

Start display at page:

Download "Cast3M www-cast3m.cea.fr"

Reynard O’Brien’
5 years ago
Views:

SHIVA, APPLICATION BASÉE SUR SALOME POUR LE

1 SHIVA, APPLICATION BASÉE SUR SALOME POUR LE DIMENSIONNEMENT DES COUVERTURES TRITIGÈNES DES RÉACTEURS À FUSION : CONCEPT ET PREMIÈRES APPLICATIONS Cast3M www-cast3m.cea.fr Julien Aubert, G. Aiello (CEA/DEN/DM2S/SEMT/BCCR) T. Charras (CEA/DEN/DM2S/SEMT/LM2S) A. Li Puma (CEA/DEN/DM2S/SERMA/LPEC) F. Kloss, N. Gore (CEA/DEN/DM2S/STMF/LGLS) 21 NOVEMBRE 2013

2 CONTENT CONTEXT OBJECTIVE TEST CASE PLATFORM WORKFLOW DEVELOPMENTS & RESULTS SALOME developments CAST3M developments Results CONCLUSION AND PERSPECTIVES 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 2

3 CONTEXT WP12 DTM01 IRM, Garching 30 Oct 2012, SALOME Platform PAGE 3

CONTEXT FUSION REACTOR : ITER - DEMOnstration REACTOR ITER (International Thermonuclear Experimental Reactor) is a large-scale scientific experiment that aims to demonstrate

4 CONTEXT FUSION REACTOR : ITER - DEMOnstration REACTOR ITER (International Thermonuclear Experimental Reactor) is a large-scale scientific experiment that aims to demonstrate that it is possible to produce commercial energy from fusion. 500 MWth ; 400 sec. DEMOnstrator power plant 2000 MWth. ITER 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 4

(self-sufficiency) Power production: Provide high grade heat production and removal ITER DEMO Coil Shielding

5 CONTEXT FUSION REACTOR : ITER - DEMOnstration REACTOR Blanket (first structure surrounding the plasma) 3 goals Shielding coils and surrounding structures against irradiation damage Tritium breeding (self-sufficiency) Power production: Provide high grade heat production and removal ITER DEMO Coil Shielding YES YES Tritium breeding Tests with TBM YES ITER Power production - YES 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 5

6 CONTEXT Present EU DEMO blanket Multi Module Segment chosen for blanket modules Different sizes Different shapes Toroidal axis Poloidal axis Radial axis 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 6

7 CONTEXT Present EU DEMO blanket Multi Module Segment chosen for blanket modules Different sizes Different shapes Toroidal axis Poloidal axis Radial axis Four concepts of blankets investigated in the present EU programme Helium Cooled Lithium Lead (HCLL) Helium Cooled Pebble Bed (HCPB) Water Cooled Lithium Lead (WCLL) Dual Coolant Lithium Lead (DCLL) 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 7

8 CONTEXT Present EU DEMO blanket Multi Module Segment chosen for blanket modules Different sizes Different shapes Toroidal axis Poloidal axis Radial axis Four concepts of blankets investigated in the present EU programme Helium Cooled Lithium Lead (HCLL) Tested in the EU ITER TBM programme* Helium Cooled Pebble Bed (HCPB) Water Cooled Lithium Lead (WCLL) Dual Coolant Lithium Lead (DCLL) 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 8

breeder, neutron multiplier and tritium carrier Pb-Li slowly re-circulating (10/50 rec/day)

9 CONTEXT Present EU DEMO blanket HCLL concept HCLL equatorial outboard module Main features EUROFER as structural He at 8MPa, Tinl/Tout = 300/500 C coolant Pb-Li (Li at 90% in 6Li ) breeder, neutron multiplier and tritium carrier Pb-Li slowly re-circulating (10/50 rec/day) T extraction from Pb-Li and He (coolant) 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 9

CONTEXT Present EU DEMO blanket HCPB concept Main features EUROFER as

Breeder (CB): Li 4 SiO 4 in form of pebble bed.

Beryllium as multiplier in form of a pebble bed T extraction (Be and CB)

10 CONTEXT Present EU DEMO blanket HCPB concept Main features EUROFER as structural material He at 8MPa, T inl /T out = 300/500 C coolant Ceramic Breeder (CB): Li 4 SiO 4 in form of pebble bed. Li 2 TiO 3 is an alternative candidate. Beryllium as multiplier in form of a pebble bed T extraction (Be and CB) through low pressure purge He (few bars) HCPB outboard blanket module & BU details (adaptation to PPP&T DEMO) 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 10

11 CONTEXT Present EU DEMO blanket - Loads Key components since it has to withstand extremely severe operating conditions Thermal loads: Heat flux on the FW [MW.m -2 ] Fth 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 11

CONTEXT Present EU DEMO blanket - Loads Key components since it has to withstand extremely severe operating conditions Thermal loads: Heat flux on the FW [MW.

12 CONTEXT Present EU DEMO blanket - Loads Key components since it has to withstand extremely severe operating conditions Thermal loads: Heat flux on the FW [MW.m -2 ] Nuclear heating: radial and poloidal power density [Wm -3 ] profiles from nuclear analyses 16,00 14,00 FW Radial power density (*10^^3 W.mm^-3) 12,00 10,00 8,00 6,00 4,00 H2O (*10^-3 W.mm^-3) LiPb (*10^-3 W.mm^-3) Eurofer (*10^-3 W.mm^-3) 2,00 0, R (mm) 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 12

13 CONTEXT Present EU DEMO blanket - Loads Key components since it has to withstand extremely severe operating conditions Thermal loads: Heat flux on the FW [MW.m -2 ] Power density distributions in poloidal and radial directions. Mechanical loads Thermal field from thermal analyses Pressure inside coolant channels Pressure on faces of stiffener plates, back plate and first wall structure 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 13

14 CONTEXT Present EU DEMO blanket - Loads Key components since it has to withstand extremely severe operating conditions Thermal loads: Heat flux on the FW [MW.m -2 ] Power density distributions in poloidal and radial directions. Mechanical loads Thermal field from thermal analyses Pressure inside coolant channels Pressure on faces of stiffener plates, back plate and first wall structure Electro-Magnetic loads Magneto-hydrodynamics loads 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 14

15 CONTEXT Present EU DEMO blanket - Loads Key components since it has to withstand extremely severe operating conditions Used in SHIVA Thermal loads: Heat flux on the FW [MW.m -2 ] Power density distributions in poloidal and radial directions. Mechanical loads Thermal field from thermal analyses Pressure inside coolant channels Pressure on faces of stiffener plates, back plate and first wall structure Electro-Magnetic loads Magneto-hydrodynamics loads 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 15

16 OBJECTIVE WP12 DTM01 IRM, Garching 30 Oct 2012, SALOME Platform PAGE 16

17 OBJECTIVE Assessment of the Open Source SALOME software for the use as design platform for In Vessel Components 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 17

OBJECTIVE Assessment of the Open Source SALOME software for the use as design platform for In Vessel Components In 2012, European Fusion Development Agreement (EFDA) task to compare platforms :

18 OBJECTIVE Assessment of the Open Source SALOME software for the use as design platform for In Vessel Components In 2012, European Fusion Development Agreement (EFDA) task to compare platforms : SALOME and ANSYS WB. Build a computation chain similar to the ANSYS WB workflow but built with the Open Source SALOME platform, and using in-house codes. Evaluation performed on a scenario where SALOME is used to chain together CEA codes (for instance, CAST3M, Trio-U, TRIPOLI). ANSYS WB used by KIT 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 18

19 OBJECTIVE Assessment of the Open Source SALOME software for the use as design platform for In Vessel Components Thermal-hydraulic / thermal calculations to assess the performances of t he BU in terms of the heat extraction capability. Results from neutronic calculations as input for heat sources, namely the breakdown of the power deposited in the subcomponents of the TBM as well as the corresponding radial and poloidal power density distributions. Simplified thermal hydraulic modeling: explicit 1D FE model of the He flow; where heat exchange with coolant takes place the FE mesh has been coupled with the 3D FE mesh; 1D FE model allows calculating the temperature evolution along the helium flow path on the basis of the local enthalpy balance between the coolant and the structure; Gnielinsky correlation to estimate HTC as a function of fluid properties, flow parameters, bulk temperature of the fluid and wall temperature Steady state analyses at first followed by transient analysis taking into account the characteristics of a typical ITER D-T plasma pulses.. Thermo-mechanical calculations Resulting temperature field in the box structure Mechanical loads corresponding to the assessed situation (pressure ) Relevant stress/strain intensities in the subcomponents identified and compared to allowable limits defined by design codes (SDC IC, RCC MRx) in case of normal and accidental (LOCA) situations. 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 19

20 OBJECTIVE Assessment of the Open Source SALOME software for the use as design platform for In Vessel Components Several aspects of the computational chain : CAD stage: import of the CAD model from a STEP file; creation of the CAD model directly with SALOME tools ability of SALOME to provide proper input for the different codes of the workflow (either via mesh or CAD transmission); ability of SALOME to provide interface tools for the data transfer between the different codes of the workflow. 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 20

21 PLATFORM WORKFLOW WP12 DTM01 IRM, Garching 30 Oct 2012, SALOME Platform PAGE 21

22 PLATFORM WORKFLOW 22 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 22

23 TEST CASE WP12 DTM01 IRM, Garching 30 Oct 2012, SALOME Platform PAGE 23

TEST CASE Test case: Simplified CAD model of a HCPB Test Blanket

defined with the foreseen tools of the platform (i.e. roughly

computation, model preparation for mechanical computation, and

24 TEST CASE Test case: Simplified CAD model of a HCPB Test Blanket Module Test case performed following the workflow such as defined with the foreseen tools of the platform (i.e. roughly CAD import, model preparation for thermal analysis, thermal computation, model preparation for mechanical computation, and mechanical analysis). 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 24

25 TEST CASE Steel component FW VSG HSG Fluid (He) CAP BP + Manifolds 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 25

26 TEST CASE Steel component FW Fluid (He) 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 26

27 DEVELOPMENTS & RESULTS SALOME DEVELOPMENTS: GEOMETRY AND MESHING WP12 DTM01 IRM, Garching 30 Oct 2012, SALOME Platform PAGE 27

28 DEVELOPMENTS & RESULTS - SALOME Description Input: a STEP file Recommendations for the STEP file Ouput: a MED file Algorithm Scripting Improvments realized in SALOME 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 28

29 DEVELOPMENTS & RESULTS - SALOME Where SALOME is used? Preparing the geometry to be meshed Managing the generation of the mesh Launching the computation with Cast3m Using YACS to build a computation chain Visualizing the result Script built for SHIVA : Preparing the geometry and managing the generation of the mesh : Retrieving the STEP file with some naming rules Creating geometrical groups for boundaries conditions Meshing Creating mesh groups from geometrical groups Creating 1D mesh from 3D fluid mesh Saving in Cast3M format file 19 NOVEMBRE ème journée utilisateurs Salomé, 21 November 2013 PAGE 29

liquid parts - void parts Specific naming rules to respect 19

30 DEVELOPMENTS & RESULTS - SALOME Description of the input STEP file Complete description of the CAD : - metallic parts - liquid parts - void parts Specific naming rules to respect 19 NOVEMBRE ème journée utilisateurs Salomé, 21 November 2013 PAGE 30

DEVELOPMENTS & RESULTS - SALOME Recommendations for the STEP file : Naming rules General naming rules for the sub-solids: FLUID_MAN : Collectors FLUID_xxxxx : Helium BERYLLIUM : Beryllium Li4SiO4 :

31 DEVELOPMENTS & RESULTS - SALOME Recommendations for the STEP file : Naming rules General naming rules for the sub-solids: FLUID_MAN : Collectors FLUID_xxxxx : Helium BERYLLIUM : Beryllium Li4SiO4 : Lithium xxxxx : Metallic part Specific naming rules for the metallic and the helium part : Metallic part : - maincomponent_subarea Helium part : - FLUID_mainComponent_numPipe.numsubPipe (for sub-solid) - FLUID_mainComponent_numPipe.numsubPipe.numFace (for all faces of the sub-solid) o numface = 1 => start of the pipe o numface = nbfaces => end of the pipe 19 NOVEMBRE ème journée utilisateurs Salomé, 21 November 2013 PAGE 31

DEVELOPMENTS & RESULTS - SALOME Description of the meshing

Beryllium part Lithium part Building mesh groups for

component one group for each subsolid Constraints for

part coarse on beryllium and lithium parts quadratic

32 DEVELOPMENTS & RESULTS - SALOME Description of the meshing result Meshing half of the box Meshing only : Metallic part Beryllium part Lithium part Building mesh groups for computation : one group for each part one group for each component one group for each subsolid Constraints for meshing : less than D elements fine on metallic part coarse on beryllium and lithium parts quadratic meshing 19 NOVEMBRE ème journée utilisateurs Salomé, 21 November 2013 PAGE 32

33 DEVELOPMENTS & RESULTS - SALOME Algorithm Geom part Opening the step file Distributing subsolids in function of their material Cutting half the shape Automatic reparing the cutted shape Preparing the repared shape for computation Creating geometrical 2D and 3D groups Generating the Helium lines Creating geometrical 1D groups 19 NOVEMBRE ème journée utilisateurs Salomé, 21 November 2013 PAGE 33

34 DEVELOPMENTS & RESULTS - SALOME Algorithm Meshing part Meshing the final shape Converting to quadratic Creating mesh groups Saving in a med file Converting to sauv format 19 NOVEMBRE ème journée utilisateurs Salomé, 21 November 2013 PAGE 34

35 DEVELOPMENTS & RESULTS - SALOME Improvements realized in SALOME for project Retrieving the path of a pipe Defining a point strictly on a face 19 NOVEMBRE ème journée utilisateurs Salomé, 21 November 2013 PAGE 35

36 DEVELOPMENTS & RESULTS CAST3M DEVELOPMENTS AND RESULTS WP12 DTM01 IRM, Garching 30 Oct 2012, SALOME Platform PAGE 36

37 DEVELOPMENTS & RESULTS CAST3M Development of the MODE ADVECTION 'ADVECTION' it represents an exchange with a fluid moving inside a pipe. The temperature of the fluid is constant in a cross section of the pipe. This is a non linear part of the heat transport behavior Supported by SEG2 type elements (only with linear elements). 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 37

38 DEVELOPMENTS & RESULTS CAST3M Development of the MODE ADVECTION 'ADVECTION' it represents an exchange with a fluid moving inside a pipe. The temperature of the fluid is constant in a cross section of the pipe. This is a non linear part of the heat transport behavior Supported by SEG2 type elements (only with linear elements). Development of an operator automatically linking the temperature of the fluid line with the one of the wall option TUYA (e.g. RIG1 = RELA TUYA GEOM1 GEOM2 CRIT); Assuming a pipe (GEOM1 seg2 or seg3) in which the temperature of the fluid can be considered constant this operation will impose that each points of GEO2 has the same temperature as its projected point on the line. The distance between one point and its projected one must be less than CRIT. T1 T1 T2 T1 T2 T1 T1 T2 T2 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 38

DEVELOPMENTS & RESULTS RESULTS HCPB TBM CAD test case calculation using Cast3M Non-Linear calculation Material properties as fonction of temperature Possibility to use transient thermal loading (as

39 DEVELOPMENTS & RESULTS RESULTS HCPB TBM CAD test case calculation using Cast3M Non-Linear calculation Material properties as fonction of temperature Possibility to use transient thermal loading (as fonction of time) tetrahedral elements (linear) Thermal models : Advection, Convection & Conduction Heat transfert coefficient (as fonction of temperature) using Gnielinsky formulae Mass flow distribution -> 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 39

40 DEVELOPMENTS & RESULTS RESULTS HCPB TBM CAD test case calculation using Cast3M Loading conditions : Temperature inlet in the first manifold 300 C set to a virtual node inside the manifold Heat flux on the First Wall 0.5 MW/m2 Power deposition 3D nuclear heating data (ASCI file) - Could be integrated in SHIVA since it comes from Neutronic analysis. 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 40

41 DEVELOPMENTS & RESULTS RESULTS Thermal analysis results : Module. 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 41

42 DEVELOPMENTS & RESULTS RESULTS Thermal analysis results : First Wall. 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 42

43 DEVELOPMENTS & RESULTS RESULTS Thermal analysis results : VSG, CAP, HSG

44 DEVELOPMENTS & RESULTS RESULTS Mechanical analysis results : Module. 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 44

45 CONCLUSION AND PERSPECTIVES WP12 DTM01 IRM, Garching 30 Oct 2012, SALOME Platform PAGE 45

46 CONCLUSION Performed activities Salomé script and development Import of the CAD model Modify of the CAD for thermal model (1D lines, 2D surfaces, ) Meshing CAST3M development New operator developed and improved (ADVE, RELA) Thermal analysis Simplified CFD Material properties as function of temperature Structural Analysis Potentials of SALOME platform (SHIVA) for IVC multiphysics calculation platform Thermal results +/- 5 % compared to ANSYS WB results for the test case. Used in 2013 for the thermal analysis of DEMO HCLL blanket concept Flexibility and adaptability of tools and development teams Possibility of implement ad hoc models (e.g. for EM calculations) Level of user friendliness far from ANSYS platform 7 ème journée utilisateurs Salomé, 21 November 2013 PAGE 46

47 PERSPECTIVE Application based on SALOME : SHIVA Export CAD with.xao files Create dialog boxes for data settings Launch computation with Cast3M Use YACS for chaining Use MEDCoupling / MEDLoader library to : manipulate meshes and fields exchange data possibility to use some tools (Ex : REMAPPER) to send some informations between 2 meshes 19 NOVEMBRE ème journée utilisateurs Salomé, 21 November 2013 PAGE 47

Thank you for your attention PAGE 48 WP12 DTM01 IRM, Garching 30 Oct 2012, SALOME Platform Commissariat à l énergie atomique et aux énergies alternatives Centre de Saclay

48 Thank you for your attention PAGE 48 WP12 DTM01 IRM, Garching 30 Oct 2012, SALOME Platform Commissariat à l énergie atomique et aux énergies alternatives Centre de Saclay Gif-sur-Yvette Cedex T. +33 (0) F. +33 (0) DEN DM2S SEMT Etablissement public à caractère industriel et commercial RCS Paris B

Segmentation of internal components and impact on maintenance for DEMO

Segmentation of internal components and impact on maintenance for DEMO F. Cismondi E. Magnani L.V. Boccaccini This presentation includes the works done under EFDA Contracts by CEA, EFET, FZK, and HAS 1