shiva, application basÉe sur salome pour le … · 2014. 1. 10. · script built for shiva :...

SHIVA, APPLICATION BASÉE SUR SALOME POUR

LE DIMENSIONNEMENT DES COUVERTURES

TRITIGÈNES DES RÉACTEURS À FUSION :

CONCEPT ET PREMIÈRES APPLICATIONS

21 NOVEMBRE 2013

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)

Cast3Mwww-cast3m.cea.fr

CONTENT

CONTEXT

OBJECTIVE

TEST CASE

PLATFORM WORKFLOW

DEVELOPMENTS & RESULTS

7ème journée utilisateurs Salomé, 21 November 2013 | PAGE 2

SALOME developments

CAST3M developments

ResultsCONCLUSION AND PERSPECTIVES

CONTEXT

| PAGE 3

WP12 DTM01 IRM, Garching 30 Oct 2012, SALOME Platform

CONTEXT

FUSION REACTOR : ITER - DEMOnstration REACTOR


ITER

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.

CONTEXT

FUSION REACTOR : ITER - DEMOnstration REACTOR


Blanket

(first structure surrounding the plasma)

ITER

ITER DEMO

Coil Shielding YES YES

Tritium breeding Tests with TBM YES

Power production

- YES

3 goals

Shielding coils and surrounding structures against irradiation damage

Tritium breeding (self-sufficiency)

Power production: Provide high grade heat production and removal

CONTEXT

Present EU DEMO blanket


Poloidal axis

Radial axis

Toroidal axis

Multi Module Segment chosen for blanket modules

Different sizes

Different shapes

CONTEXT



Poloidal axis

Radial axis

Toroidal 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 )


Different sizes

Different shapes

CONTEXT



Poloidal axis

Radial axis

Toroidal 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 )

Tested in the EU ITER TBM programme*


Different sizes

Different shapes

CONTEXT

Present EU DEMO blanket – HCLL concept


HCLL equatorial outboard module

Main featuresEUROFER 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)

CONTEXT

Present EU DEMO blanket – HCPB concept


HCPB outboard blanket module & BU details (adaptation to PPP&T DEMO)

Main featuresEUROFER as structural material

He at 8MPa, Tinl/Tout = 300/500 °C coolant

Ceramic Breeder (CB): Li4SiO4 in form of pebble bed. Li2TiO3 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)

CONTEXT

Present EU DEMO blanket - LoadsKey components since it has to withstand extremely severe operating conditions


Thermal loads:� Heat flux on the FW [MW.m-2]

Fth

CONTEXT



Thermal loads:� Heat flux on the FW [MW.m-2]� Nuclear heating : radial and

poloidal power density [Wm-3]profiles from nuclear analyses

0,00

2,00

4,00

6,00

8,00

10,00

12,00

14,00

16,00

0 100 200 300 400 500 600 700 800 900

Ra

dia

l p

ow

er

de

nsi

ty (

*1

0^

^3

W.m

m^

-3)

R (mm)

H2O (*10^-3 W.mm^-3)

LiPb (*10^-3 W.mm^-3)

Eurofer (*10^-3 W.mm^-3)

FW

CONTEXT



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

CONTEXT








Electro-Magnetic loadsMagneto-hydrodynamics loads

CONTEXT








Electro-Magnetic loadsMagneto-hydrodynamics loads

Used in SHIVA

OBJECTIVE

| PAGE 16


OBJECTIVE

Assessment of the Open Source SALOME software for the use as design platform for In Vessel Components


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

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 hel ium flow pathon 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 temperatureSteady state analyses at first followed by transient analysis taking into account the characteristics of a typical ITER D-T plasma pulses. .

• Thermo-mechanical calculationsResulting 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.


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 toolsability 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.


PLATFORM WORKFLOW

| PAGE 21


PLATFORM WORKFLOW

22


TEST CASE

| PAGE 23


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).


TEST CASE


Fluid (He)

Steel component VSG HSG

CAP BP + Manifolds

FW

TEST CASE


Steel component FWFluid (He)


SALOME DEVELOPMENTS:GEOMETRY AND MESHING

| PAGE 27


DEVELOPMENTS & RESULTS - SALOME

Description

Input: a STEP file

Recommendations for the STEP file

Ouput: a MED file

Scripting

Algorithm

Improvments realized in SALOME



Where SALOME is used ?Preparing the geometry to be meshedManaging the generation of the meshLaunching the computation with Cast3mUsing YACS to build a computation chainVisualizing the result

Script built for SHIVA : Preparing the geometry and managing the generation of the mesh :

Retrieving the STEP file with some naming rulesCreating geometrical groups for boundaries conditionsMeshingCreating mesh groups from geometrical groupsCreating 1D mesh from 3D fluid meshSaving in Cast3M format file

19 NOVEMBRE 2013



Description of the input STEP file

Complete description of the CAD :- metallic parts- liquid parts- void parts

Specific naming rules to respect

19 NOVEMBRE 2013



Recommendations for the STEP file : Naming rules

General naming rules for the sub-solids:FLUID_MAN : CollectorsFLUID_xxxxx : HeliumBERYLLIUM : BerylliumLi4SiO4 : Lithiumxxxxx : Metallic part

Specific naming rules for the metallic and the helium part :Metallic part :- mainComponent_subAreaHelium 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 pipeo numFace = nbFaces => end of the pipe

19 NOVEMBRE 2013



19 NOVEMBRE 2013

Description of the meshing result

Meshing half of the box

Meshing only :Metallic partBeryllium partLithium part

Building mesh groups for computation :one group for each partone group for each componentone group for each subsolid…

Constraints for meshing :less than 1 200 000 3D elementsfine on metallic partcoarse on beryllium and lithium partsquadratic meshing



Algorithm – Geom part

Opening the step fileDistributing subsolids in function of their material

Cutting half the shape

Automatic reparing the cutted shape

Preparing the repared shape for computationCreating geometrical 2D and 3D groups

Generating the Helium linesCreating geometrical 1D groups

19 NOVEMBRE 2013



Algorithm – Meshing part

Meshing the final shape

Converting to quadratic

Creating mesh groups

Saving in a med file

Converting to sauv format

19 NOVEMBRE 2013



Improvements realized in SALOME for project

Retrieving the path of a pipe

Defining a point strictly on a face

19 NOVEMBRE 2013



CAST3M DEVELOPMENTSAND RESULTS

| PAGE 36


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 behaviorSupported by SEG2 type elements (only with linear elements).


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 behaviorSupported 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.

T2

T1

T1

T1

T1T1

T2T2

T2


DEVELOPMENTS & RESULTS – RESULTS

HCPB TBM CAD test case calculation using Cast3MNon-Linear calculation� Material properties as fonction of temperature� Possibility to use transient thermal loading (as fonction of time)

1 113 405 tetrahedral elements (linear)

Thermal models : Advection, Convection & Conduction

Heat transfert coefficient (as fonction of temperature) using Gnielinsky formulae

Mass flow distribution ->



HCPB TBM CAD test case calculation using Cast3M

Loading conditions :

� Temperature inlet in the first manifold300 °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.



Thermal analysis results : Module.




Thermal analysis results : First Wall.


Thermal analysis results : VSG, CAP, HSG


Mechanical analysis results : Module.


CONCLUSION AND PERSPECTIVES

| PAGE 45


CONCLUSION

Performed activitiesSalomé script and development� Import of the CAD model � Modify of the CAD for thermal model (1D lines, 2D surfaces, …)� MeshingCAST3M 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 multiphysicscalculation platform

Thermal results +/- 5 % compared to ANSYS WB results for the test case.Used in 2013 for the thermal analysis of DEMO HCLL blanket conceptFlexibility 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


PERSPECTIVE

Application based on SALOME : SHIVAExport CAD with .XAO filesCreate dialog boxes for data settingsLaunch computation with Cast3MUse YACS for chaining

Use MEDCoupling / MEDLoader library to :manipulate meshes and fieldsexchange datapossibility to use some tools (Ex : REMAPPER) to send some informations between 2 meshes

19 NOVEMBRE 2013


DENDM2SSEMT

Commissariat à l’énergie atomique et aux énergies alternatives

Centre de Saclay | 91191 Gif-sur-Yvette Cedex

T. +33 (0)1 69 08 96 56 | F. +33 (0)1 69 08 69 05

Etablissement public à caractère industriel et commercial | RCS Paris B 775 685 019

| PAGE 48


Thank you for your attention

shiva, application basÉe sur salome pour le … · 2014. 1. 10. · script built for shiva :...

Documents