Home

CodeSaturne tutorial

1. OE I Figure V 33 Solution control Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 75 120 When using Fortran routines 1t is sometimes useful to allocate pre defined user arrays that are present in every sub routine This allocation can be specified in the User arrays item under the Calculation management heading It is not the case in the present calculation The item Memory management allows to set the memory size for the calculation It is the size of the integer and real arrays that will be used to store most of the variables in the Fortran parts of Code_Saturne It is dependent on the number of cells in the mesh In parallel mode it depends on the number of cells treated by each processor and not the total number of cells For this simple case the default values are appropriate Ex Integer array Calculation features a J Mobile mesh Array size NITUSE 0 Turbulence models N Thermal model 0 x Number of cells with halo Radiative transfers o o gt Pa S Physical properties O e Lario eee L Reference values 0 x Number of boundary faces Fluid properties __ Ee J Gravity hydrostatic pres 3 Volume conditions Real array i J Volume regions definition Ba dla Array size NRTUSE 0 joe Initialization Additional scalars 0 x Number of cells with halo Definition and i
2. Figure V 56 Solution control Probes Switch to the Calculation management heading to prepare the launch script and run the calculation Code_Saturne documentation Page 98 120 EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial 3 SOLUTION FOR CASE 3 Only a few elements are different from case 2 In this case the density becomes variable Go to the item Fluid properties under the heading Physical properties and change the nature of the density from constant to user law Click on the highlighted icon and define the user law in the window that pops up Follow the format used in the tab Examples Sx Density 4 Identity and paths E By Calculation environment EER o fa Thermophysical models a va AL i constant E Physical properties 3 user law kg m Reference values i a Fluid properties __ user subroutine usphyv e h Gravity hydrostatic pressure Viscosity B Volume conditions gt E BY Additional scalars constant y amp BY Boundary conditions sr BF Numerical parameters Reference value pu 18 951e 05 Pa s y Calculation control f Calculation management Specific heat constant Reference value Cp 5483 0 J kg K Thermal conductivity constant e Reference value 2 0 02495 Wim K a L User expression Predefined symbols Examples rho Tempt 4 0666e 3 Tempe 5 0754e 2 1000 9
3. The Courant number CFL and Fourier number will be removed from the post processing results Eventually probes will be defined for chronological records following the data given in figure III 4 Then the total pressure will be deactivated from all probes and the Velocity U will be only activated on probes 1 2 6 7 and 8 In addition the domain boundary will be post processed This allows to check the boundary conditions and especially that of the temperature and passive scalar 3 7 Calculation restart After the first run the calculation will be continued for another 400 time steps The calculation restart is managed through the Graphical Interface 3 8 Results Figure III 9 shows the time evolution of temperature recorded on each monitoring probe 300 Probe 5 Probe 2 5 Probe 3 E 200 Probe 4 a Probe 5 5 Probe 6 Probe 7 Probe 8 Probe 9 100 0 10 20 30 Time s Figure 111 9 Time evolution of temperature at monitoring probes for case 3 Figure 111 10 shows the velocity fields colored by temperature in the first run of calculation Figure 11 11 shows the velocity fields in the second calculation restart of the first one this can be very useful to save some disk space if some variables are of no interest as post processing files can be large EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial Code_Saturne documentation Page 26 120 t 0 1s 1 55 T_inlet 24 8C gt
4. _ a T mm Add from Prepocessor listing Import groups and references from Preprocessor listing Figure V 20 Creation of boundary regions symmetry region Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 63 120 The same treatment must be done for the wall conditions All colors 2 3 4 6 and 7 can be grouped in a single boundary zone After defining all the boundary zones the Interface window will look as in figure V 21 Definition of boundary regions _ _ _ eem bre E D walls E REN 2or3or4or6or7 L il Il EI Add from Prepocessor listing Import groups and references from Preprocessor listing Figure V 21 Creation of boundary regions Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 64 120 Now that the boundary zones are defined the boundary conditions assigned to them will be specified Click on the item Boundary conditions to set the inlet boundary conditions for velocity and turbulence As shown on figure V 22 outlet and wall boundary zones also appear in the window Ex Boundary conditions Identity and paths Calculation environment Meshes selection Mesh quality criteria Thermophysical models Calculation features Mobile mesh Turbulence models Thermal model Radiative transfers Conjugate heat transfer Additional scalars Definition and ini
5. Radiative transfers VitesseZ 0 0 m s Conjugate heat transfer Additional scalars Thermal scalar initialization Definition and initialization Physicals properties TempC 38 5 C Physical properties E 4 Reference values Turbulence initialization Fluid properties Gravity hydrostatic pressure Initialization by reference velocity for all zones hd Volume conditions E Mi J Volume regions definition P initialization Reference velocity 0 03183 m s fb Head losses B Boundary conditions BF Numerical parameters B Calculation control By Calculation management EN lele Figure V 67 Thermophysical models Initialization Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 110 120 In the item Fluid properties under the heading Physical properties enter the following information Variable Type Value Density 998 671 kg m Viscosity 0 445 x 107 kg m s Specific Heat 4182 88 J kg 2C Thermal Conductivity Constant 0 601498 W m K71 For density and viscosity the value given here will serve as a reference value see user manual for details Density Identity and paths E Calculation environment user law By Meshes selection Laia Mesh quality criteria O Thermophysical models i Calculation features i Mobile mesh Viscosity Turbulence models l Thermal model user law M R
6. User Scalar Diffusion Coefficient Associated Type of Type of Scalar coefficient value Dscall scalar_2 constant HEGI 8 95e 05 Identity and paths Calculation environment Meshes selection Mesh quality criteria 3 3 Thermophysical models i Calculation features Mobile mesh Turbulence models h Thermal model J Radiative transfers J Conjugate heat transfer 3 gj Additional scalars Definition and initialization MN Physicals properties By Physical properties E fa Volume conditions B Boundary conditions Numerical parameters 6 Bl Calculation control B Calculation management Figure V 43 Additional scalar User scalar physical properties Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 85 120 Create the boundary zones The procedure is the same as in case 1 but the colors are different Note that colors 5 and 32 have completely disappeared in the joining process they are now internal faces and are not considered as boundaries while some boundary faces of color 24 remain Create the inlet outlet and symmetry boundary zones with the following colors e inlet color 1 e outlet color 34 e symmetry colors 8 9 28 29 38 39 EJES Definition of boundary regions PEACE 1 DG il Inlet _ Identity and paths sa Calculation environment Meshes selection _ Mesh quality criteria Thermophysical
7. OK Cancel Figure V 57 Fluid properties Variable density EDF R D Code_Saturne version 2 0 0 rc1 tutorial Code_Saturne documentation Page 99 120 As the density is variable the influence of gravity has to be considered properties go to Gravity hydrostatic pressure and set the value of each component of the gravity vector h Identity and paths B Calculation environment B Thermophysical models E Physical properties i Reference values Fluid properties MN Gravity hydrostatic pressure By Volume conditions B Additional scalars BY Boundary conditions Numerical parameters B Calculation control A B Calculation management E F E al ago Gravity In the heading Physical Gravity along X 0 0 m s Gravity along Y m s Gravity along Z 0 0 m s Hydrostatic pressure Pressure interpolation in stratified flow standard O improved Figure V 58 Fluid properties Gravity Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 100 120 Add a monitoring point close to the entry boundary condition in the Output control item Points X m Y m Z m 2 Identity and paths Loser E amp Calculation environment ARO RAEE Aaa BR Addiionei acera eg Physical properties Reference values Fluid properties 4 Gravity
8. Mesh quality criteria 3 Thermophysical models i Calculation features Mobile mesh J Turbulence models __ Thermal model J Radiative transfers Conjugate heat transfer Outputs listings Output listing at each time step E Post processing Post processing every n time steps gt 2 Fluid domain post processing x Domain boundary post processing x Additional scalars Definition and initialization TE Type of post processing for mesh fixed M oe Physical properties Post processing format EnSight Gold y Reference values l J Fluid properties Options Ma Gravity hydrostatic pressure E Volume conditions format binary hd 4 Volume regions definition polygons display e n Initialization i Head losses polyhedra display y amp Boundary conditions E _ ary big endian E pe Definition of boundary regi _ Boundary conditions E Numerical parameters e J Time step h Equation parameters iek Global parameters 3 Calculation control J Time averages MN Output control iL Volume solution control h Surface solution control Profiles B Calculation management e Lale Figure V 51 Output control post processing Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1
9. By Additional scalars B Calculation control Initialization by reference velocity for all zones E Physical properties Volume regions definition B Boundary conditions Turbulence initialization E y Calculation management Reference velocity 1 0 m s Velocity initialization al Volume conditions Thermal scalar initialization A Initialization c e Numerical parameters ol lt P Ul ag Figure V 11 Initialization of dynamic variables Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 54 120 The initial value for the thermal scalar also appears in the item Definition and Initialization under the heading Additional scalars where more options concerning the scalars can be specified The value of the initial value can be modified in any of the two pages But in case there are additional scalars i e other than the thermal scalar their initialization is only possible in the Additional scalars page Sx Select volume region for initial value h Identity and paths Calculation environment Region selection all_cells y i i Meshes selection Mesh quality criteria Thermophysical models Calculation features Mobile mesh Turbulence models Thermal model Radiative transfers Physical properties E 3 Volume conditions a Volume regi
10. e smooth wall rough wall _ Sliding wall Prescribed flux i Exchan icient Figure V 25 Scalars boundaries walls Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 68 120 Click on inlet to choose the temperature inlet value Here this value is 300 C The default value is left for the outlet Scalars Exchange Coefficient Figure V 26 Scalars boundaries inlet Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 69 120 The calculation parameters need then to be specified under the header Numerical parameters Go to the item Steady flow management to specify the number of iterations 30 in this case The default value of the relaxation coefficient will be kept and the Zero iteration option will not be activated E x Steady flow algorithm management a Identity and paths Calculation environment Relaxation coefficient 0 9 Meshes selection Mesh quality criteria Number of iterations restart included 30 ea 5 Thermophysical models Zero iterations option O of Calculation features Mobile mesh Turbulence models Thermal model i Radiative transfers gj Physical properties i Reference values Fluid properties i Gravity hydrostatic pressure Volume conditions i Volume regions definition Initialization Additional scalars
11. 3 amp Additional scalars EnerTurb Definition and initialization Physicals properties Physical properties ViscTurb _ Reference values _ Fluid properties J Gravity hydrostatic pressure scalar_2 Volume conditions NbCourant Volume regions definition Initialization NbFourier Head losses 12345678 12345678 12345678 12345678 12345678 12345678 12345678 12345678 12345678 12345678 12345678 12345678 Pression VitesseY total_pressure a Dissip TempC me Ge Bel Ge Re Dosi 5 Boundary conditions i Definition of boundary regi Boundary conditions Numerical parameters __ Time step Equation parameters Global parameters Calculation control Time averages Output control A Volume solution control Surface solution control i Ly Profiles B Calculation management Figure V 55 Solution control Output configuration Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 97 120 Delete all the probe numbers for the total_pressure variable No chronological record will be created for this variable As for the VelocitX variable only select probes 1 2 6 7 and 8 Time evolution on the other probes will not be recorded Solution control Print in Post Mame ia Probes listing processing Pression O Fi 12345678 total pressure
12. C The inlet temperature of water in the cold branch is 300 C Water characteristics are considered constant and their values taken at 300 C and 150 x 10 Pa e density p 725 735 kg m73 Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 11 120 e dynamic viscosity u 0 895 x 1074 kg m s7 e specific heat Cp 5 483 LE e Thermal Conductivity 0 02495 W m7 K7 1 4 Mesh characteristics Figure 11 2 shows a global view of the downcomer mesh This two dimensional mesh is composed of 700 cells which is very small compared to those used in real studies This is a deliberate choice so that tutorial calculations run fast inlet WI Figure 11 2 Geometry of the downcomer Note that here the case is two dimensional but Code_Saturne always operates on three dimensional mesh elements cells The present mesh is composed of a layer of hexahedrons created from the 2D mesh shown on figure 11 2 by extrusion elevation in the Z direction The virtual planes parallel to Oxy will have sliding symmetry conditions to account for the two dimensional character of the configuration Type structured mesh Coordinates system cartesian origin on the edge of the main pipe at the outlet level on the nozzle side figure II 1 Mesh generator used
13. EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 56 120 Under the heading Physical properties in the main list the item Reference values allows to set the reference pressure Use the default value of 101300 Pa a x Reference pressure 2 Identity and paths Calculation environment Reference value for total pressure 101300 0 Pa Meshes selection i Mesh quality criteria 5 Thermophysical models Calculation features Mobile mesh Turbulence models Thermal model i Radiative transfers eg Physical properties MN Reference values iL Fluid properties i L Gravity hydrostatic pressure 5 Volume conditions Volume regions definition Initialization eg Additional scalars J Definition and initialization i Ly Physicals properties fi Boundary conditions BF Numerical parameters E iy Calculation control B Calculation management ago Figure V 14 Physical properties reference pressure Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 57 120 Specify the fluid physical characteristics in the item Fluid properties e Density e Viscosity e Specific Heat e Thermal Conductivity In this case they are all constant Density 725 735 kg m73 Viscosity 0 895 x 1074 kg m71 s7 Specific Heat 5483 J kg 2C7 Thermal Conductivity 0 02495 W
14. Mobile mesh Turbulence models Outlet i Thermal model Radiative transfers _ Conjugate heat transfer aj Additional scalars Definition and initialization Physicals properties 3 Physical properties i Reference values Add from Prepocessor listing Fluid properties Gravity hydrostatic pressure Import groups and references from Preprocessor listing Volume conditions Volume regions definition J Initialization i Head losses gj Boundary conditions MN Definition of boundary regi J Boundary conditions By Numerical parameters a Calculation control B Calculation management AA RAY Inlet Figure V 73 Boundary regions EDF R D Code_Saturne version 2 0 0 rc1 tutorial Code_Saturne documentation Page 115 120 For the dynamic boundary conditions the velocity is 0 03183 m s in the z direction and the hydraulic diameter 0 4 m for both inlets Socso Ex Identity and paths Calculation environment Meshes selection Mesh quality criteria Thermophysical models Calculation features Mobile mesh Turbulence models i Thermal model Radiative transfers Conjugate heat transfer ditional scalars Definition and initialization Physicals properties ysical properties Reference values Fluid properties Gravity hydrostatic pressure lume conditions Volume regions definition Initialization Head losses Boundary co
15. Opening a new file Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 45 120 The interface automatically updates the following information Study name e Case name Directory of the case Associated sub directories of the case File Tools Window Help JeBOB E Ex Study T_JUNCTION Case CASE1 XML file Ces Directory of the case ilas ETUDE TUTORIAL_SAT PICS T_JUNCTION CASE1 Associated sub directories of the case By Calculation environment a Thermophysical models ly Physical properties E iy Volume conditions B Additional scalars E Bl Boundary conditions i BY Numerical parameters fi Calculation control Gf Calculation management Data DATA Results RESU User subroutines SRC Running scripts SCRIPTS Directory of meshes MESH Figure V 3 Identity and paths Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 46 120 Save the case to give a name to the new XML file by opening the File menu and clicking on Save as A new window will appear enter the name of the case in File Name then click on Save Remember to save the case regularly throughout the preparation of the calculation X Save File As 2 mnx Look in jeg nome nicol JCASEL DATA Oo n EJ Mm Compu leg THCH ca nicolas n File name casl Files of type Code Saturne GUI file xml ui Cancel 4 F
16. This mesh is composed of 1650 cells which is very small compared to those used in real studies This is a deliberate choice so that tutorial calculations run fast i o GRATA TTT LI OMITE y a TOTO z TTT Y ATTO Di Ul Figure 111 2 View of the full domain mesh with zoom on the joining regions Type block structured mesh Coordinates system cartesian origin on the edge of the main pipe at the outlet level on the nozzle side figure 111 2 Mesh generator used SIMAIL and mesh joining with the Preprocessor of Code_Saturne in order to deal with hanging nodes Color definition see figure 111 3 1 5 Summary of the different calculations Three cases will be studied with this geometry The following table gives a summary of their different characteristics which makes temperature a passive scalar but it is only for simplification purposes Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 17 120 34 38 39 23 Figure III 3 Colors of the boundary faces CASE CASE 2 unsteady flow additionnal passive scala
17. hydrostatic pressure Volume conditions B Boundary conditions Numerical parameters eg Calculation control Time averages I Output control h Volume solution control h Surface solution control iL Profiles Calculation management Figure V 59 New monitoring probe Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 101 120 After completing the interface before running the calculation some Fortran user routines need to be modified Go to the folder SRC REFERENCE base and copy usclim f90 in the SRC directory usclim f90 In this case usclim f90 is used to specify the time dependent boundary condition for the temperature Refer to the comments in the routine or to the Code_Saturne user manual for more information on this routine In our case you need to identify the boundary faces of color 1 The command CALL GETFBR 1 NLELT LSTELT will return an integer NLELT corresponding to the number of boundary faces of color 1 and an integer array LSTELT containing the list of the NLELT boundary faces of color 1 Note that the string 1 can be more complex and combine different colors group references or geometrical criteria with the same syntax as in the Graphical Interface For each boundary face IFAC in the list the Dirichlet value is given in the multi dimension array RCODCL as follows IF TTCABS LT 3 8DO THEN DO IELT 1 NLELT IFAC
18. 0 5 gt X 3 Definition of boundary regions Nature Selection criteria 1 Inlet Symmetry 8 or 9 or 28 or 29 or 38 or 39 Outlet 34 Wall 24 and 0 1 lt x and 0 5 gt L L Add from Prepocessor listing Import groups and references from Preprocessor listing E Figure V 45 Creation of a wall boundary region 3Note that due to the joining process there are in fact no boundary faces of color 24 with X coordinate outside the 0 1 0 5 intervalle The geometrical criterium is therefore not necessary It is presented here to show the capacity of the face selection module Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 87 120 Define the other wall boundary zones The faces of color 6 have to be divided in two separate zones based on a geometrical criterium on Y Definition of boundary regions 4orFor2lor22or23 6 and Y gt 1 6 and Y 1 n 24 and 0 1 lt X and 0 52 X i Add from Prepocessor listing Import groups and references from Preprocessor listing Figure V 46 Creation of wall boundary regions Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 88 120 The dynamic boundary conditions are the same as in case 1 for the inlet and there are still no sliding walls Ces Boundary conditions Identity and paths DORG A _ _ i i i Meshes selection pay e
19. 1 i if gt ponen inlet inlet 1 Mesh quality criteria i 3 Thermophysical models Calculation features Mobile mesh Turbulence models Thermal model Radiative transfers Conjugate heat transfer ditional scalars Definition and initialization Physicals properties ysical properties Reference values Fluid properties Gravity hydrostatic pressure Direction lume conditions Volume regions definition specified coordinates Initialization Head losses 34 20r3 4 or7or2lo0 6 and Y gt 1 6 and Y lt 1 31 or 33 24 and 0 1 lt E 00 NOU NES Boundary conditions X Definition of boundary regi A Boundary conditions B Numerical parameters a Calculation control B Calculation management Calculation by hydraulic diameter v Hydraulic diameter m Figure V 47 Dynamic variables boundary inlet Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 89 120 To configure the scalar boundary conditions on the walls select individually each wall in the item Boundary conditions On all the walls a default homogeneous prescribed flux is set for temperature and prescribed values are specified for the passive scalar according to the following table aILI Prescribed value 0 wall 3 Prescribed value 0 all 6 Boundary conditions inlet outlet 34 wall Zora wall 4orv7or l or 2 wall 6 and Y gt 1 6 and Y lt 1 al or 33 24 and 0
20. 1 lt X E Opt Ch Ol Ri F Smooth or rough wall e smooth wall rough wall Sliding wall Scalars Scalar Name Figure V 48 Scalars boundaries wall_5 Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 90 120 Click on inlet to set the inlet values for the scalars 300 C for temperature and 200 for the passive scalar Exchange Tempc Prescribed scar Preserbed BR Figure V 49 Scalars boundaries inlet Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 91 120 Some calculation parameters now need to be defined Go to the item Time step under the heading Numerical parameters In our case the time step is Uniform and constant Set the number of iterations to 300 and the reference time step to 0 05 s Ex Unsteady flow algorithm management identity and paths Calculation environment Time step option Uniform and constant M _ Meshes selection J Mesh quality criteria i Reference time ste 0 05 s A eg Thermophysical models Calculation features Number of iterations restart included 300 n Mobile mesh ue Turbulence models Time step limitation with 2 Thermal model the local thermal time step Radiative transfers Conjugate heat transfer amp Additional scalars J Definition and initialization i Physicals properties aj Physical p
21. 20 120 25C 8 000e 01 0 000e 00 Figure III 5 View of the boundary domain colored by the scalar_2 variable Case 2 Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 21 120 t 0 1s t 1 1s ittttttttttity SALIDA ARIADNA VAL ELLE ALLELLA LLADIN mann tini ALULLE LLADIN i ini COLADA TO ini MINTI CO00 000010111 VINI STEEPED yay x Temp C x Temp C 3 0008 02 3 000e 02 x 2 300e 02 x 2 300e 02 1 600e 02 1 600e 02 9 000e 01 9 000e 01 2 000e 01 2 000e 01 i Temp C Temp C 3 0008 02 3 000e 02 x 2 300e 02 2 300e 02 1 600e 02 1 600e 02 9 000e 01 9 000e 01 2 000e 01 2 000e 01 y Temp C Temp C 3 0008 02 3 0008 02 x 2 300e 02 2 300e 02 1 600e 02 1 600e 02 9 000e 01 9 000e 01 2 000e 01 2 000e 01 Figure III 6 Water velocity field colored by temperature at different time steps Case 2 Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 22 120 3 CASE 3 Time dependent boundary conditions and vari able fluid density In this case some boundary conditions will be time dependent and some physical characteristics of the fluid will be dependent on the temperature 3 1 Calculation options The options for this case are the same as in case 2 except for the variable fluid density Flow type unsteady flow Time step uniform and constant Turbulence model k e
22. Definition and initialization Physicals properties 5 Boundary conditions Definition of boundary regi Boundary conditions E Numerical parameters A Steady flow management h Equation parameters h Global parameters Bg Calculation control 6 Calculation management 4 ae Figure V 27 Steady flow management Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 70 120 After selecting the item Equation parameters the tab Scheme allows to change different more advanced numerical parameters In this case none of them should be changed from their default value Solver Scheme Identity and paths B Calculation environment Blending Slope Flux Meshes selection Scheme Factor Test 2constructic Mesh quality criteria S Thermophysical models llei Centered Calculation features VitesseY Centered Mobile mesh J Turbulence models Thermal model EnerTurb Upwind i Radiative transfers J Physical properties Reference values Temp C Centered Fluid properties Gravity hydrostatic pres Volume conditions Volume regions definition Initialization Additional scalars Definition and initialization Physicals properties 3 3 Boundary conditions Definition of boundary re J Boundary conditions gt Numerical parameters Steady flow management Equation parameters _ Global parameters Calculation contro
23. LSTELT IELT RCODCL IFAC ISCA 1 1 20 D0 100 DO TTCABS ENDDO ELSE DO IELT 1 NLELT IFAC LSTELT IELT RCODCL IFAC ISCA 1 1 ENDDO ENDIF 400 DO ISCA 1 refers to the first scalar and TTCABS is the current physical time See the example file in the directory examples for the complete usclim f90 file Note that although the inlet boundary conditions for temperature are specified in the usclim f90 file it is necessary to specify them also in the Graphical Interface The value given in the Interface can be anything it will be overwritten by the Fortran routine After updating the Fortran file run the calculation as explained in case 2 Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 102 120 When a calculation is finished Code_Saturne stores all the necessary elements to continue the compu tation in another execution with total continuity These elements are stored in several files grouped in a RESTART xxxxxxxx directory in the RESU directory In this case after the first calculation is finished a second calculation will be run starting from the results of the first one Go directly on the item Start Restart under the heading Calculation management Activate the Cal culation restart by ticking the on box Then click on the folder icon next to it to specify the restart files to use Start Restart a Identity and paths Calculation environment Ca
24. T_inlet 165C gt UITTITTTTTTTIT ANN TTT TTT TTT Temp C Temp C 400 400 305 305 Y 210 210 115 20 x t 3 0s t 4 0s T_inlet 315C gt T_inlet 400C gt Temp C Temp C 400 400 305 305 i nr 210 210 i 115 115 20 20 t 9 0s T_inlet 400C gt t 15 0s T_inlet 400C gt Temp C Temp C 400 NS Y 305 aie 305 210 SR a 210 115 g 115 x 20 i 20 Figure III 10 Water velocity field colored by temperature and inlet temperature value at different time steps first calculation HEE PERERA DON CUE MOM RUD HT MILAN Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 27 120 t 15 1s Temp C Temp C 460 400 205 305 ani 210 115 ps 20 Temp C gt j f N Temp C 210 Pi aw 115 SEE 115 20 l a Figure 111 11 Water velocity field colored by temperature and inlet temperature value at different time steps second calculation Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 28 120 4 CASE 4 Head loss parallelism and spatial average This case will be run in parallel on two processors Head loss will be used to simulate the presence of an obstacle in the flow and the spatial average of the temperature will be calculated at each time step 4 1 Calculation options The options for this case are the same as in case 3 F
25. endian Figure V 32 Output control The Monitoring Points Coordinates tab allows to define specific points in the domain monitoring probes where the time evolution of the different variables will be stored in historic files In this case no monitoring points are defined Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 74 120 The item Volume solution control allows to specify which variable will appear in the output listing in the post processing files or on the monitoring probes In this case the default value is kept where every variable is activated Ex Solution control i Identity and paths Calculation environment Print in Post z Meshes selection Mesh quality criteria l Thermophysical models x Calculation features x _ Mobile mesh J Turbulence models a Thermal model x J Radiative transfers sn Physical properties x Reference values x J Fluid properties pee r Gravity hydrostatic pressure _ 3 Volume conditions isc x Volume regions definition ra Ee x d Initialization 3 Additional scalars x Definition and initialization x J Physicals properties Boundary conditions Definition of boundary regi Boundary conditions sa Numerical parameters Calculation control J Output control A Volume solution control i i J Profiles 3 8 Calculation management
26. first point e X 0 25 m e Y 2 25m eZ O0m Monitoring points coordinates Figure V 53 Output controls monitoring points 1 point Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 95 120 Repeat the procedure for the other probes Their coordinates are indicated in the following table the Z coordinate is always 0 Figure V 54 Output control monitoring points Remember to save the Xml file regularly Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 96 120 Go to the item Volume solution control to define which variables will appear in the listing the post processing and the chronological records Uncheck the boxes in front of the Pressure Tubulent energy and Dissipation variables in the Print in listing column Information on these three variables will not appear in the output listing anymore Uncheck the boxes in front of the Courant number and Fourier number variables in the Post processing column These variables will be removed from the post processing results es Solution control identity and paths Calculation environment i if Meshes selection Print in Post li listing processing x x x J Mesh quality criteria 3 Thermophysical models Calculation features VitesseX Mobile mesh Turbulence models __ Thermal model VitesseZ Radiative transfers J Conjugate heat transfer
27. m7 k7 EJES Density Identity and paths Calculation environment ne y y i J Meshes selection e ee aes Reference value p 725 735 kg m i Calculation features Mobile mesh Viscosity Turbulence models i Thermal model constant y A Radiative transfers S Physical properties Reference value u 8 951e 05 Pa s J Reference values A Fluid properties r Specific heat h Gravity hydrostatic pressure oe Volume conditions _ gt i 4 Volume regions definition 5 J Initialization o Additional scalars Reference value Cp 5483 kg K Definition and initialization J Physicals properties Thermal conductivity H B Boundary conditions A a Be Numerical parameters constant M Da E fly Calculation control B Calculation management Reference value 0 02495 Wim K a lt p gt Figure V 15 Physical properties fluid properties Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 58 120 Set the three components of gravity in the Gravity hydrostatic pressure item In this case since the gravity doesn t have any influence on the flow gravity can be set to 0 As for the pressure interpolation interpolation method keep the standard default value Ex Gravity Identity and paths E 3 Calculati
28. parameters necessary to this study can be defined through the Graphical Interface except the variable fluid characteristics and the advanced post processing features that have to be specified in user routines Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 38 120 Parameters of calculation control Reference time step 100 Maximal Fourier number 60 705 The time step limitation by gravity effects will also be activated 2 5 Output management The standard options for output management will be used Four monitoring points will be created at the following coordinates DE E 0 01534 0 011765 0 031652 0 031652 725 Two vertical temperature profiles will be extracted at the following locations profill6 x 1 6 y 0 2 lt z lt 0 2 m profil32 x 3 2 y 0 0 2 lt z lt 0 2 m 2 6 User routines The following routines have to be copied from the folder SRC REFERENCE base into the folder SRC usdpst f90 usvpst f90 and usmpst f90 In this test case advanced post processing features will be used A clip plane will be created along the symmetry plane of the domain on which the temperature will be written This plane will be added to the standard writer i e it will be an extra part in the standard CHR ENSIGHT case The periodicity of output on the standard writer will be 10 iterations An additional writer will also be created with a periodicity of 5 it
29. part and each variable to write The arguments for PSTEVA are e IPART part number e NAMEVR character string of the name under which the variable will be written e IDIMT dimension of the variable 1 or 3 for scalars or vectors e IENTLA for vectors indicates if the components are interlaced 1 or not 0 e IVARPR shortcut option for specific situations set to 0 here e NTCABS current time step passed to usupst f90 with the right value e TTCABS current physical time passed to usupst f90 with the right value e TRACEL array for variables on cells e TRAFAC array for variables on internal faces e TRAFBR array for boundary faces Part 1 only contains internal faces so only TRAFAC needs to be filled Execute a loop on all the faces from the LSTFAC list For each of them the temperature will be stored in TRAFAC The temperature at each face will be calculated by interpolation from the value at the centers of the two neighboring cells The numbers of the neighbors of face IFAC are IFACEL IFAC 1 and IFACEL IFAC 2 For a proper linear interpolation see in the TEST_CASES directory for the use of the POND parameter yielding the fractionnal position of the face on the line joining the two cell centers Note that in parallel computing the cells on both side of the face can be managed by different pro cessors In order for the interpolation to be correct a parallel synchronization must be done before the loop A similar problem happens
30. present some advanced post processing techniques 1 2 Description of the configuration The configuration is based on a real mock up designed to characterize thermal stratification phenomena and associated fluctuations The geometry is shown on figure IV 1 R 600 SUPERNIMBUS model dimensions mm _ inner diameter 600 zan umana aaa Lodi R 600 aoe Gravi ty 1200 6000 4000 outlet cold inlett hot inlet Figure IV 1 Geometry of the case There are two inlets a hot one in the main pipe and a cold one in the vertical nozzle The volumic flow rate is identical in both inlets It is chosen small enough so that gravity effects are important with respect to inertia forces Therefore cold water creeps backwards from the nozzle towards the elbow until the flow reaches a stable stratified state 1 3 Characteristics Characteristics of the geometry Diameter of the pipe D 0 40 m Characteristics of flow Cold branch volume flow rate Hot branch volume flow rate Cold branch temperature 1 18 26 C Tha 38570 The initial water temperature in the domain is equal to 38 5 C Water specific heat and thermal conductivity are considered constant and calculated at 18 26 C and 10 Pa e heat capacity Cp 4 182 88 Lig Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 36 120 e thermal conductivity A 0 601498 W
31. tutorial documentation Page 93 120 In this case chronological records on specified monitoring probes are needed To define the probes click on the Monitoring points Coordinates tab es Output Control Monitoring Points Coordinates Identity and paths E l 3 Calculation environment Monitoring points recording CES Mesh qualty criteria Monitoring points files at each time step _ v Thermophysical models e Calculation features Monitoring points coordinates Mobile mesh 5 Turbulence models LS Sa CS A Se Thermal model Radiative transfers Conjugate heat transfer Additional scalars Definition and initialization Physicals properties Physical properties Reference values Fluid properties Gravity hydrostatic pressure Volume conditions Volume regions definition Initialization Head losses Boundary conditions Definition of boundary regi Boundary conditions Numerical parameters Time step Equation parameters J Global parameters 3 Calculation control Time averages Output control Volume solution control Surface solution control i Profiles B Calculation management Figure V 52 Output control monitoring points Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 94 120 Click on Add and enter the coordinates of the monitoring points you want to define For the
32. wee ew eee ew ew ee eee ee aw eee eS Li MESH COs lt lt aces e ewe ewe a we we ER Ee eR 2 CASE 5 Stratified junction lt 6 244 444 nm sas OR asa SH ES Dl CALCULATION OFTIONS ee sewa ados 2 2 INITIAL AND BOUNDARY CONDITIONS LL 2 3 VARIABLE DENSITY AND DYNAMIC VISCOSITY 2 0 e e a ee 2a PARAMETERS ek REE ERA Ree Eee ee eee O s 20 OWVIPUT MANAGEMENT sicarios dass a SERRE wed 20 LE ee eh eee EEE eK EE eH da Pd Re ni Er 0 0 akan ec keke Bees hehe adora A V STEP BY STEP SOLUTION L SOLUTION FOR CASEI cocci n deal 22 22 22 23 23 23 24 20 20 28 28 28 29 29 30 30 ol ol 34 39 39 39 30 36 36 36 37 37 37 38 38 39 42 43 Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial Ee age 2 SOLUTION PUR CASE Z2 Liu e a ee ee EE e E a dee a i 80 3 SOLUTION FOR CASE J 4444 86 268 hahahah 98 SULUTIOR FOR CAGES soso ieee tbe see ideale ec bes ge ws 105 S BSOLUTION FOR CASE Pia asa AAA AAA a 109 Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial i age Part I INTRODUCTION Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation 1 Page 8 120 Introduction Code_Saturne is a system designed to solve the Navier Stokes equations in the cases of 2D 2D ax isymmetric or 3D flows Its main module is designed for the simulation of flows which may be steady or unsteady laminar or turbulent incompressible or potentially dilatable isotherma
33. with periodic boundary conditions Hence the calling of routines PARCOM and PERCOM shown in the example in the TEST_CASES directory As for part 2 it contains only cells so only TRACEL need be filled For each cell in the LSTCEL list just set TRACEL to the value of the temperature at the center of the cell
34. 1 2 defines the correspondance between the colors and the type of boundary condition to use 34 Outlet 234672122233133 Wall 24 for 0 1 lt X lt 0 5 8928293839 Symmetry Table 111 2 Boundary faces colors and associated references 3 3 Variable Density In this case the density is a function of temperature the variation law is defined in the Graphical User Interface although it can also be defined in a Fortran user routine The expression is p T AT B C IIL1 where p is the density T is the temperature A 4 0668 x 1073 B 5 0754 x 107 and C 1000 9 In order for the variable density to have an effect on the flow gravity must be set to a non zero value y 9 8le will be specified in the Graphical Interface 3 4 Parameters All the parameters necessary to this study can be defined through the Graphical Interface except the time dependent boundary conditions that have to be specified in user routines Parameters of calculation control Number of iterations 300 Reference time step Output period for post processing files In order to join the separate meshes into a single domain colors 5 24 and 32 will have to be joined through the Graphical Interface 3 5 User routine The following routine has to be copied from the folder SRC REFERENCE base into the folder SRC usclim f90 usclim f90 This routine allows to define advanced boundary conditions on the boundary faces Even if uscli
35. 1 tutorial documentation Page 112 120 The aim of the calculation is to simulate a stratified flow It is therefore necessary to have gravity Set it to the right value in the item Gravity hydrostatic pressure In order to have a sharper stratification the pressure interpolation method will be set to improved Ex Gravity _y Identity and paths a Calculation environment Gravity along X 0 0 m s _ Meshes selection Mesh quality criteria Gravity along Y 0 0 m s Thermophysical models A _ 2 Pr Calculation features Gravity along Z 9 81 m s Mobile mesh i J Turbulence models Hydrostatic pressure J Thermal model Radiative transfers i i Conjugate heat transfer O standard Additional scalars Definition and initialization Physicals properties Physical properties Reference values Fluid properties Gravity hydrostatic pressure Volume conditions Volume regions definition Initialization Head losses H B Boundary conditions fi Numerical parameters E a Calculation control B Calculation management Pressure interpolation in stratified flow improved i Lel Figure V 71 Fluid properties Gravity Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 113 120 Go to the item Definition and initialization under the heading Additiona
36. 7 63s Figure IV 3 Evolution of temperature Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 41 120 Figure IV 4 Sub domain where the temperature is lower than 21 C upper figure and localization in the full domain lower figure Part V STEP BY STEP SOLUTION Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 43 120 1 SOLUTION FOR CASE 1 The first thing to do before running Code_Saturne is to prepare the computation directories In this first example the study directory T_JUNCTION will be created containing a single calculation directory CAS1 This is done by typing the command cree_sat etude T_JUNCTION CASI The mesh files should be copied in the directory MESH The Code_Saturne Graphical Interface is launched by typing the command SaturneGUI in the DATA subdirectory of the CASI directory The following graphic window opens fig V 1 File Tools Window Help JEBHOBb E Welcome to SaturneGUI Figure V 1 User interface Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 44 120 Go to the File menu and click on New file to open a new calculation data file as shown in the figure Ved File Tools Window Help New file Ctrl N ls Open Ctrl 0 Recent file Ah Save Ctrl S Ak Save as GH Close Ctrl W Ee Quit Ctrl Q Welcome to SaturneGUl Figure V 2
37. Code_Saturne version 2 0 0 rc1 tutorial documentation Page 13 120 2 3 Parameters and User routines All parameters necessary to this study can be defined through the Graphical Interface without using any user Fortran files Calculation control parameters Number of iterations Relaxation coefficient 0 9 Output period for post processing files 2 4 Results Figure II 4 presents the results obtained at different iterations in the calculation They were plotted from the post processing files with EnSight Note since the steady flow option has been chosen the evolution of the flow iteration after iteration has no physical meaning It is merely an indication of the rapidity of convergence towards the physical steady state iteration 1 iteration 10 cu FEH CEI 300 Temperature C saes 270 200 208 Y 150 145 100 x 82 20 x iteration 20 Temperature C Temperature C 300 275 Y 250 225 200 x Figure II 4 Water velocity field colored by temperature at different time steps Part III FULL DOMAIN Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 15 120 1 General description 1 1 Objective This aim of this case is to tackle the merging of initially separate meshes into a single fluid domain The questions of mesh joining and hanging nodes will be addressed The test case will then be used to present more complex calculations with
38. EDF R amp D y Y FLUID DYNAMICS POWER GENERATION AND ENVIRONMENT DEPARTMENT SINGLE PHASE THERMAL HYDRAULICS GROUP 6 QUAI WATIER F 78401 CHATOU CEDEX TEL 33 1 30 87 75 40 FAX 33 1 30 87 79 16 Code_Saturne documentation Code_Saturne version 2 0 0 rc1 tutorial contact saturne support Qedf fr http www code saturne org APRIL 2010 EDF 2010 Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial Es age TABLE OF CONTENTS I INTRODUCTION Y E SOC aLaaa OO REE OEE EHO EERE A II SIMPLE JUNCTION TEST CASE 9 1 General description ics eke He wee Ree CHG EARS EHO GEE BS Hs 10 Lol TORTIE ba eh hee Re hee eee Cae oes oe 10 1 2 DESCRIPTION OF THE CONFIGURATION gt s peaos bsos esad puaa edas baei 10 Lo CHARACTER TIOS serrera ia aara e a ei 10 la MESH CHARACTERISTIOS gt o o iege iape kaaa a a a a a Be we a a 11 2 CASE 1 Basic calculation gt s gt s ss rore aa A 11 Sl LAT ATON OFTON seeen ena maaa a dd a we A 11 2 2 INITIAL AND BOUNDARY CONDITIONS aaa a ee 12 29 PARAMETERS AND USER ROUTING cs seos da E I Rw AR ee 13 ARI onerosa rosa REAL LR 13 III FULL DOMAIN 14 1 General description seso sra sa e a e oe Ew Y 15 Li RE lt lt iu lee Gh ewe es he se esa le lg 15 1 2 DESCRIPTION OF THE CONFIGURATION 15 Lar AUARACTREISTI S sass 2 dica AA A ee a e G 15 1 4 MEH CONARACTEMISTIOS gt Lupi GAARA een 16 1 5 SUMMARY OF THE DIFFERENT CALCULATIONS 2 00 eee eee ee eee 16 2 CASE 2
39. Numerical parameters B Calculation control B Calculation management E 0 0 0 0 E E a Ep Figure V 9 Thermal scalar conservation list of models Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 52 120 In the present case select Temperature Celsius degrees Thermal scalar Temperature Celsius Figure V 10 Thermal scalar conservation choice of a model There are no radiative transfers in our case so the item is ignored Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 53 120 To initialize variables at the instant t 0 s go to the item Initialization Here the velocity the thermal scalar and the turbulence can be initialized In this case the default values can be kept zero velocity an initial temperature of 20 C and a turbulence level based on a reference velocity of 1 m s Specific zones can be defined with different initializations In this case only the default all cells is used EJES Select volume zone for initialization Identity and paths Calculation environment Volume zone all_cells v Meshes selection Mesh quality criteria 3 Thermophysical models Calculation features Vitessex 0 0 m s i Mobile mesh E J Turbulence models VitesseY 0 0 m s J Thermal model i Radiative transfers VitesseZ 0 0 m s Temp C 20 0 pr
40. PO des J B Volume conditions pic des Simail NOPO des J B Boundary conditions B Numerical parameters By Calculation control E y Calculation management lx Join meshes Select boundary faces from single meshes to join optional References Groups Reverse IA 29 m 24 32 0 1 _ Subdivide warped faces Correct cell and face orientations Figure V 40 Join a Mesh Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 82 120 In this case Unsteady flow must be selected in the Analysis features item J Identity and paths Calculation environment Meshes selection Mesh quality criteria 3 Thermophysical models MA Calculation features J Mobile mesh Turbulence models Thermal model Radiative transfers Conjugate heat transfer BY Additional scalars B Physical properties amp Volume conditions Boundary conditions B Numerical parameters Calculation control B Calculation management E Ex Steady Unsteady flow algorithm unsteady flow z Eulerian Lagrangian multi phase treatment single Phase Flow y AAA A Atmospheric flows off y E combustion off m Pulverized coal combustion off Electrical models cn EES Figure V 41 Thermophysical models Analysis features Unsteady flow The rest of the heading Thermophysical models is identical t
41. Passive scalar with various boundary conditions and output management 17 il CAMU ATON OFTON keke eae ERE ee eee ei 17 2 2 INITIAL AND BOUNDARY CONDITIONS gt ieissar dedere redere erra 17 2 3 PARAMETERS AND USER ROUTINES gt sosea cacadna aned a rR 18 diet UOTPUT MANAGEMENT io ne ees Eww REE Ee AAA 19 IE RA eee RA A 19 Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial i age 3 CASE 3 Time dependent boundary conditions and variable fluid density nL TADOOLATI N OPTIONS oa iu ed RR i 5 2 INITIAL AND BOUNDARY CONDITIONS er cse ce se ses e e Ew Ho dl VARAPLE DERSITY ean ek ee bee de e x Gwe PRI e e AAA AR A oe SERIE esr ER RRR ESS 2 0 OUTPUT MANAGEMENT 44 6654 Beebe hee eee Hee eRe RHA wae EDDA dr CALCULATION RESTART 4656444408244 484 amp HAH DSR EDOM Dee SA Oe DI na oe we ARE EERE ss e a a A 4 CASE 4 Head loss parallelism and spatial average o 4A L CALCULATION OPTIONS ila e ERED EE RR ee i 4 2 INITIAL AND BOUNDARY CONDITIONS 4 6 ae de e e ol ko VARIABLE DENSITY ciar EA ih Sk BRL irene Lied ae AAA OREO RO ERE ee ee Ree ee ee i PARAM ee cos be ee hw he ee eee eo ee eo ree eo eee EO O Ivi Be ee a RAE AA A AAA AAA a TTPO MANACGOMENT ecc LEE ww A a e ei sn Feet a sica ras LR ARL sI ee IV STRATIFIED JUNCTION 1 General description escocia or leal SG aos hehehehe ewe daa 1 2 DESCRIPTION OF THE CONFIGURATION 1 0 eee ee et ee 0 Lo CHARACTER TOS ce
42. R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 107 120 The calculation of the spatial average is done in the usproj f90 routine Refer to the example file in the directory examples for the complete usproj f90 file The other two changes are controlled in the item Prepare batch analysis To run the calculation on two processors simply change the number of processors indicator to 2 The launch script will automatically deal with the rest EJES Computer selection a gt Identity and paths L sapos B Calculation environment Workstation o By Thermophysical models i i B Additional scalars a a 4 B Physical properties Select the batch script file E runcase B Volume conditions B Boundary conditions Prepare batch calculation H Re Numerical parameters H B Calculation control Number of processors E E E Calculation management AZZ User arrays serra Li Memory management a Start Restart Advanced options gt MS Prepare batch calculation Code_Saturne batch running 4 ag Figure V 65 Number of processors Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 108 120 As seen in paragraph 4 the file moy dat created by usproj f90 will be written in the temporary execution directory It must be identified in the launch script in order to be automatically co
43. Reference values Fluid properties Gravity hydrostatic pressure Add Volume conditions D Volume regions definition Add from Prepocessor listing i Initialization gt 23 Additional scalars Import groups and references from Preprocessor listing i Definition and initialization i LL Physicals properties 5 Boundary conditions MN Definition of boundary regi Ly Boundary conditions E Y Numerical parameters E Bg Calculation control BF Calculation management Figure V 17 Creation of a boundary region Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 60 120 Each boundary must be defined Click on Add to edit a new boundary The boundary faces will be grouped in user defined zones based on their color or on geometrical conditions For each zone a reference number a label a nature and a localization condition must be assigned The different natures that can be assigned are e wall e inlet e symmetry e outlet The Label can be any character string It is used to identify the zone more easily It usually corresponds to the nature of the zone The Zone number can be any integer It will be used by the code to identify the zone No specific order or continuity in the numbering is needed The Localization is used to define the faces that belong to the zone It can be a color number a group reference geometrical conditions on a combination of them rela
44. SIMAIL Color definition see figure 11 3 To specifiy boundary conditions on the boundary faces of the mesh the latter have to be identified It is commonly done by assigning an integer to each of them characteristic of the boundary region they belong to This integer is refered to as color or reference 2 CASE 1 Basic calculation 2 1 Calculation options Most of the options used in this calculation are default options of Code_Saturne Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 12 120 Di Figure 11 3 Colors of the boundary faces Flow type steady flow Turbulence model k e Scalar s 1 temperature Physical properties uniform and constant 2 2 Initial and boundary conditions Initialization none default values The boundary conditions are defined as follows e Flow inlet Dirichlet condition an inlet velocity of 1 m s and an inlet temperature of 300 C are imposed e Outlet default values e Walls default values Figure II 3 shows the colors used for boundary conditions and table 11 1 defines the correspondance between the colors and the type of boundary condition to use Do not forget to enter the value of the hydraulic diameter adapted to the current inlet used for turbulence entry conditions 5 Outlet _ 5 80 Symmetry Table 11 1 Boundary conditions and associated references Code_Saturne EDF R amp D
45. Scalar s 1 temperature 2 passive scalar Physical properties uniform and constant except density Management of monitoring points 3 2 Initial and boundary conditions Initialization 20 C for temperature default value 10 for the passive scalar The boundary conditions are defined in the user interface and depend on the boundary zone The time dependence of the temperature boundary condition implies the use of a Fortran user routine see below e Flow inlet Dirichlet condition an inlet velocity of 1 m s a time dependent inlet temperature and a value of 200 for the passive scalar are imposed e Outlet default value e Walls velocity pressure and thermal scalar default value passive scalar different conditions depending on the color and geometric parameters The boundary conditions for the passive scalar are identical as in case 2 as specified in the following table Wall Diriehlet _0_ wall_3 Dirichlet _0__ The wall_1 to wall_6 regions are defined as follows through color references and geometric local ization Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 23 120 Label Color and geometric parameters wall1 24 and 0 1 lt X and X lt 0 5 wall_2 2 or 3 wall_3 4 or 7 or 21 or 22 or 23 wall_4 6 and Y gt 1 wall_5 6 and Y lt 1 wall_6 31 or 33 Figure 111 3 shows the colors used for boundary conditions and table 11
46. adiative transfers o i e Conjugate heat transfer Reference value u 0 001445 Pa s E 3 Additional scalars Definition and initialization i Physicals properties 6 Physical properties Reference value p 998 671 kg m Specific heat Reference values constant ly Di MN Fluid properties i J Gravity hydrostatic pressure Reference value Cp 4182 88 J kg K 6 Volume conditions Volume regions definition Thermal conductivity Initialization i Head losses constant M gt H B Boundary conditions i B Numerical parameters Reference value 2 0 601498 WIm K ff Calculation control i iy Calculation management Lale Figure V 68 Physical properties fluid properties Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 111 120 For the density and viscosity enter the expressions of the user laws as showed in figures V 69 and V 70 in the windows poping while clicking on the highlighted boxes User expression Predefined symbols rho Tempt 4 0668e 3 Tempt 5 0754e 2 1000 9 Figure V 69 Variable density User expression Predefined symbols m Tempe Tempt 3 4016e 9 Tempt 6 2332e 7 4 557 78 5 1 591536 9 Figure V 70 Variable viscosity Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc
47. aturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 37 120 6 Hot inlet Table IV 1 Boundary faces colors and associated references 2 2 Initial and boundary conditions Initialization temperature initialization at 38 5 C The boundary conditions are defined as follows e Flow inlet Dirichlet condition velocity of 0 03183 m s for both inlets temperature of 38 5 C for the hot inlet temperature of 18 6 C for the cold inlet e Outlet default value e Walls default value Figure IV 2 shows the colors used for boundary conditions and table IV 1 defines the correspondance between the colors and the type of boundary condition to use 2 3 Variable density and dynamic viscosity In this case the density and the dynamic are functions of the temperature the following variation laws are specified in the Graphical User Interface p T AT B 4 C IV 1 where p is the density T is the temperature A 4 0668 x 1073 B 5 0754 x 107 and C 1000 9 For the dynamic viscosity the variation law is u T T AM T BM CM DM IV 2 where u is the dynamic viscosity T is the temperature AM 3 4016 x 107 BM 6 2332 x 107 CM 4 5577 x 107 and DM 1 6935 x 1073 In order for the variable density to have an effect on the flow gravity must be set to a non zero value g 9 81e will be specified in the Graphical Interface 2 4 Parameters All the
48. e of post processing for mesh fixed X Equation parameters Post processing format EnSight Gold Global parameters Calculation control Options Domain boundary post processing _ Time averages I Output control format binary ly Volume solution control T yoons dsplay gt Profiles polygons isplay p alculation management polyhedra display big endian Figure V 79 Output management Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 119 120 Output Control Monitoring Points Coordinates bel Identity and paths Calculation environment E 3 Thermophysical models amp 3 Additional scalars Physical properties E amp Volume conditions Monitoring points coordinates Volume regions definition i J Initialization n X Y AL LIL 1 0 010025 0 01534 0 011765 oundary conditions Definition of boundary regi 3 1 625 0 01534 0 031652 Boundary conditions merical parameters 3 3 225 0 01534 0 031652 Time step Equation parameters 4 3 8726 Global parameters alculation control Time averages PBB Output control Volume solution control Profiles ly Calculation management Monitoring points recording Monitoring points files at each time step 1 i Su o My PH ZO 0 047481 a 1 Add Delete KIO g KIC Figure V 80 Monitoring points For the advanced post proce
49. ed by a user routine 300 The calculation will be run in parallel on 2 processors In order to join the separate meshes into a single domain colors 5 24 and 32 will have to be joined through the Graphical Interface 4 6 User routines The following routines have to be copied from the folder SRC REFERENCE base into the folder SRC usclim f90 usphyv f90 usproj f90 and uskpdc f90 e usclim f90 This routine allows to define advanced boundary conditions on the boundary faces Even if usclim f90 is used all boundary conditions have to be defined in the Graphical User Interface Only the conditions that differ from this first definition need to appear in usclim f90 The boundary conditions defined in usclim f90 will replace those specified in the Graphical Interface In this case the temperature at entry is supposed variable in time following the law T 20 100t for 0 lt t lt 3 8 T 400 for t gt 3 8 111 4 where T is the temperature in C and t is the time in s Probe 9 on the inlet Temperatue C 0 l 2 3 4 5 6 T 8 9 10 11 12 13 14 15 Time s Figure 111 13 Time evolution of the temperature at inlet e usproj f90 This routine is called at the end of each time step and has access to the whole set of variables of the code It is therefore useful for many user specific post processing including the calculation of a spatial Sonly when they appear in the SRC directory will they be taken into accoun
50. ement i User arrays Memory management Start Restart Calculation restart on off Advanced options Prepare batch calculation 4 Figure V 35 Start Restart Code_Saturne documentation Page 77 120 Code_Saturne version 2 0 0 rc1 tutorial EDF R D The final item Prepare batch calculation is used to prepare the launch script and on certain architec tures launch the calculation Calculations can be launched from the Graphical Interface in interactive mode Workstation or in a PBS batch queue Management of chart PBS In this simple case choose the Workstation _ Calculation features a J Mobile mesh Turbulence models J Thermal model i Radiative transfers amp 3 Physical properties Reference values Fluid properties Gravity hydrostatic pres E amp Volume conditions J Volume regions definition Initialization amp Additional scalars Definition and initialization i gt Physicals properties S amp Boundary conditions i J Definition of boundary re i Boundary conditions 6 Numerical parameters i Steady flow management Equation parameters J Global parameters eg Calculation control J Output control Volume solution control Profiles eg Calculation management User arrays h Memory management Start Restart ry Computer selection Workstatio
51. erations It will only contain one part i e one sub mesh the set cells where the temperature is lower than 21 C The temperature will be written on this part The interest of this part is that it is time dependent as for the cells it contains Three Fortran routines will be used e usdpst f90 This routine is called only once at the beginning of the calculation It allows to define the different writers and parts e usmpst f90 This routine is called at each time step It allows to redefine the content of certain parts using any variable especially the temperature for this case e usvpst f90 This routine is called at each time step It allows to specify which variable will be written on which part l only when they appear in the SRC directory will they be taken into account by the code Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 39 120 2 7 Results Figure IV 3 shows the evolution of the temperature in the domain at different time steps The evolution of the stratification is clearly visible Figure IV 4 shows the cells where the temperature is lower than 21 C It is not an isosurface created from the full domain but a visualization of the full sub domain created through the post processing routines Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 40 120 Temperature 38 50 t 17 53s 33 44 28 38 23 32 18 26 t 9
52. he New icon to enter the list of colors to be joined E m Meshes Periodic Boundaries Identity and paths eg Calculation environment MN Meshes selection List of meshes gt i J Mesh quality criteria Meshes Format A B Thermophysical models ic des Simail NOPO des By Physical properties i E B Volume conditions fdc des Simail NOPO des ES 5 B Additional scalars ll downcomer des Simail NOPO des Fa A B Boundary conditions gt f Numerical parameters Add Delete A B Calculation control 3 B Calculation management X Join meshes Select boundary faces from single meshes to join optional References Groups Reverse Fraction o fe C mr New Delete _ Subdivide warped faces Correct cell and face orientations al ago Figure V 39 Meshes list of meshes Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 81 120 Fill the array References for the color reference to be joined 5 24 and 32 different colors can be entered on a single line separated by blanks Sx Meshes Periodic Boundaries List of meshes Identity and paths Calculation environment MN Meshes selection 7 Mesh qualit criteria Meshes Foma B Additional scalars a J B Physical properties fdc des Simail NO
53. igure V 4 Saving the XML file Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 47 120 The next step is to specify the mesh es to be used for the calculation Click on the item Solution Domain under the heading Analysis environment The list of all meshes available in the folder MESH appears in the window List of meshes Delete the mesh es you will not use In this case only the mesh downcomer des is needed _ Meshes Periodic Boundaries a Identity and paths 3 Calculation environment MN Meshes selection Thermophysical models downcomer des Simail NOPO des List of meshes 5 B Physical properties J B Volume conditions fdc des Simail NOPO des By Additional scalars J A Boundary conditions pic des Simail NOPO des amp iy Numerical parameters i i BF Calculation control BF Calculation management Add Delete Join meshes _ Subdivide warped faces Correct cell and face orientations Figure V 5 Meshes list of meshes On this item Solution Domain there are three other tabs e PERIODIC BOUNDARIES e SYRTHES COUPLING e STAND ALONES RUNNING They are not used in this case Keep the default values this operation only deletes the selected entries from the list it does not delete the mesh file in the MESH directory Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1
54. ill be deactivated from all probes and the Velocity U will be only activated on probes 1 2 6 7 and 8 In addition the domain boundary will be post processed This allows to check the boundary conditions and especially that of the temperature and passive scalar 4 8 Results Figure III 15 shows the evolution of the spatial average of the temperature Figure III 16 shows velocity fields colored by temperature The effect of the head loss modeling the obstacle is clearly visible this can be very useful to save some disk space if some variables are of no interest as post processing files can be large Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 32 120 400 300 Temperature spatial mean C to 100 Time s Figure 111 15 Evolution of spatial average of the temperature as a function of time EDF R D Code_Saturne version 2 0 0 rc1 tutorial Code_Saturne documentation Page 33 120 t 5 0s Figure 111 16 Water velocity field colored by temperature t 3 0s t 8 0s t 15 0s MPILILHI MILITIII Part IV STRATIFIED JUNCTION Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 35 120 1 General description 1 1 Objective The aim of this case is to train the Code_Saturne user on a simplified but real 3D computation It corresponds to a stratified flow in a T junction The test case will be used to
55. ion 2 0 0 rc1 tutorial documentation Page 66 120 As for the wall boundary zone the specifications the user might have to give is when the wall is sliding and if the wall is smooth or rough In this case the walls are fixed so the option is not selected and the wall is considered as smooth Note that if one of the walls had been sliding it would have been necessary to isolate the corresponding boundary faces in a specific boundary region Boundary conditions inlet inlet 1 outlet 5 2 orsord4orbor Smooth or rough wall e smooth wall rough wall Sliding wall Exchange Coefficient Figure V 24 Dynamic variables boundary walls Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 67 120 The boundary conditions on the temperature are only applied on inlets outlets and walls For the walls three conditions are available e Prescribed value e Prescribed flux e Exchange Coefficient For the outlet only Prescribed value and Prescribed flux are available but they are taken into account only when the flow re enters from the outlet Otherwise homogeneous Prescribed flux is considered by Code_Saturne For the inlets only Prescribed value is available In this case all walls are adiabatic So the boundary condition for the temperature will be a Prescribed flux set to 0 Boundary conditions inlet inlet 1 outlet 5 E Smooth or rough wall
56. ion of volume regions Label Zone Selection criteria criteria E cells eta etnia Ma 0 4 gt x al Obstacle 2 ere Mina 0 25 gt y ates perce Add Add from Prepocessor listing Import groups and references from Preprocessor listing Figure V 63 Creation of head losses region Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 106 120 To specify the head losses coefficients go to the item Head losses and select the name of the head losses volume region In this example the coefficient is isotropic so that we use the same value for each Q Please note that a 2 x Ki therefore if Ki 10 a 2 10 Sx Select volume zone for head losses J Identity and paths 3 Meshes selection f Mesh quality criteria Thermophysical models 3 Calculation features Mobile mesh Turbulence models Thermal model Obstacle 2 0 2 lt x and 0 4 gt x and 0 75 Tensor coefficients Head losses coefficients Kj 0 5 aj JU hysical properties axx 20000 0 ayy 20000 0 azz 20000 0 Volume conditions J Volume regions definition Initialization Reference frame transformation matrix Definition of boundary regi E Boundary conditions h d Boundary conditions BY Numerical parameters BF Calculation control 5 8 Calculation management Figure V 64 Head losses coefficients Code_Saturne EDF
57. ious conditions will be imposed for the passive scalar as specified in the following table Wall Dirichlet _0_ wall Dirichlet _5_ wall 3 Dirichlet 0 The wall_1 to wall_6 regions are defined as follows through color references and geometric local ization Label Color and geometric parameters wall 1 24 and 0 1 lt X and X lt 0 5 wall_2 2 or 3 wall_3 4 or 7 or 21 or 22 or 23 wall_4 6andY gt 1 wall_5 6and Y lt 1 wall_6 31 or 33 Figure III 3 shows the colors used for boundary conditions and table 11 1 defines the correspondance between the colors and the type of boundary condition to use 23467 21 22 23 24 31 33 Wall 8 9 28 29 38 39 Table 111 1 Boundary faces colors and associated references 2 3 Parameters and User routines All parameters necessary to this study can be defined through the Graphical Interface without using any user Fortran files Calculation control parameters Number of iterations Reference time step Output period for post processing files In order to join the separate meshes into a single domain colors 5 24 and 32 will have to be joined through the Graphical Interface Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 19 120 2 4 Output management In this case different aspects of output management will be addressed By default in the Graphical Interface all variables are set to appear in the listing the po
58. itions the temperature of the cold inlet is 18 6 C and that of the hot inlet is 38 5 C Boundary conditions wall 3 wall 5 cold_inlet 2 hot_inlet inlet 6 outlet outlet 7 We locity KA Dire ct ion _ _ specified coordinates Turbulence Calculation by hydraulic diameter 7 Hydraulic diameter m Scalars Exchange Figure V 76 Temperature boundary conditions Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 117 120 Boundary conditions a wall wall cold inlet inlet ihot_inlet 2 inlet outlet Welocity OO Direction specified coordinates n Turbulence Calculation by hydraulic diameter ha Hydraulic diameter mi Scalars Scalar Name Tempe Prescribed 138 5 Exchange Coefficient Figure V 77 Temperature boundary conditions Tick the appropriate box for the time step to be variable in time and uniform in space In the boxes below enter the following parameters Parameters of calculation control Number of iterations 10 Reference time step Maximal CFL number Maximal Fourier number 60 Minimal time step Maximal time step Time step maximal variation 0 1 And activate the option Time step limitation with the local thermal time step Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 118 120 Unsteady flow algorithm management J Identity and path
59. l Output control Volume solution control Profiles x x VitesseZ Centered Ei Dissip Upwind x Global parameters h Identity and paths i Calculation environment Gradient calculation method i E Se i LP iv Se L Iterative handling of non orthogonalities v Thermophysical models AA AEREA ART O EE Calculation features Multigrid algorithm for pressure xx Mobile mesh AR PRO AAA mE Turbulence models Pseudo coupled velocity pressure solver Lal Thermal model Radiative transfers l l l 3 23 Physical properties Handling of transposed gradient and divergence Reference values source terms in momentum equation Fluid properties _rT_r__ehre_i Gravity hydrostatic pres Extrapolation of pressure gradient lume conditions on pan auna x Neumann 1st order 7 Volume regions definition Initialization ditional scalars Relaxation of pressure increase Definition and initialization Physicals properties 3 3 Boundary conditions Definition of boundary re Boundary conditions merical parameters Steady flow management Equation parameters ISU q TEL Output control Volume solution control Profiles Figure V 29 Numerical parameters Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 71 120 Under the heading Calculation control click on the item Output co
60. l or not Scalars and turbulent fluctuations of scalars can be taken into account The code includes specific modules referred to as specific physics for the treatment of lagrangian particle tracking semi transparent radiative transfer gas pulverized coal and heavy fuel oil combustion electricity effects Joule effect and electric arcs and compressible flows The code also includes an engineering module Matisse for the simulation of nuclear waste surface storage Code_Saturne relies on a finite volume discretization and allows the use of various mesh types which may be hybrid containing several kinds of elements and may have structural non conformities hanging nodes The present document is a tutorial for Code_Saturne version 2 0 0 rc1 It presents five simple test cases and guides the future Code_Saturne user step by step into the preparation and the computation of the cases The test case directories containing the necessary meshes and data are available in the examples directory This tutorial focuses on the procedure and the preparation of the Code_Saturne computations For more elements on the structure of the code and the definition of the different variables it is higly recommended to refer to the user manual Code_Saturne is free software you can redistribute it and or modify it under the terms of the GNU General Public License as published by the Free Software Foundation either version 2 of the License or at your o
61. l scalars to specify the minimal and maximal values for the temperature 18 26 C and 38 5 C Note that the initial value of 38 5 C set earlier is properly taken into account Sx Select volume region for initial value Identity and paths 3 Calculation environment Region selection all_cells i Meshes selection _ Mesh quality criteria Thermophysical models Calculation features Mobile mesh Turbulence models Thermal model Tempc 138 5 J Radiative transfers l E Conjugate heat transfer amp Additional scalars E Definition and initialization d Physicals properties User scalar definition Initial Variance Minimal Maximal of scalar value value no 18 26 138 5 Physical properties i Reference values Fluid properties _ Gravity hydrostatic pressure Add Delete Volume conditions Volume regions definition _ Initialization _ Head losses B Boundary conditions Numerical parameters B Calculation control Calculation management E p E A Figure V 72 Scalar initialization Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 114 120 Create the boundary regions 6 inlet Identity and paths Calculation environment Meshes selection Mesh quality criteria Wall Thermophysical models j Inlet Calculation features
62. lculation restart on off E By Thermophysical models 7 Additional scalars S Physical properties h Reference values J Fluid properties i Gravity hydrostatic pressure Advanced options L E By Volume conditions E B Boundary conditions n l Ji BF Numerical parameters a Calculation control h d Time averages Output control Volume solution control x Surface solution control h d Profiles Calculation management User arrays Memory management MN Start Restart Prepare batch calculation Calculation on frozen velocity and pressure fields _ Ar A e ao Figure V 60 Start Restart Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 103 120 A window opens with the architecture of the study sub directories Open the RESU folder and click on the folder RESTART xxxxxxxx where xxxxxxxx corresponds to the reference of the first calculation Then click on Validate Look in amp home nicolas ETUDE T ULL DOMAIN CASE3 RESU Q o w og W Compu A nicolas E RESTART 03011736 Directory RESTART 03011736 Files of type Directories y Cancel di Figure V 61 Start Restart Selection of the restart directory Code_Saturne documentation Code_Saturne version 2 0 0 rc1 tutorial EDF R D Page 104 120 Go to the Time step item under the heading Numerical pa
63. lon Linear Production 3 Turbulence models Rij epsilon LLR Thermal model Rij epsilon SSG B Physical properties v2f phi model B Volume conditions k omega SST B Additional scalars LES Smagorinsky B Boundary conditions LES classical dynamic model B Numerical parameters LES WALE B Calculation control B Calculation management Fn en Ep a eon E Figure V 7 Turbulence model list of models Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 50 120 In this case the k model is used Turbulence model k epsilon ha Sua Advanced options Figure V 8 Turbulence model choice of a model Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 51 120 For this study the equation for temperature must be solved Click on the item Thermal model to choose between e No thermal scalar e Temperature Celsius degrees e Temperature Kelvin e Enthalpy Identity and paths ii scalar S Calculation environment No thermal scalar i Meshes selection No thermal scalar _y Mesh quality criteria Temperature Celsius S E Thermophysical models Temperature Kelvin i a Fi Calculation features Enthalpy J kg Mobile mesh i Turbulence models MS Thermal model B Physical properties Volume conditions B Additional scalars Boundary conditions B
64. low type unsteady flow Time step uniform and constant Turbulence model k e Scalar s 1 temperature 2 passive scalar Physical properties uniform and constant except density Management of monitoring points 4 2 Initial and boundary conditions Initialization 20 C for temperature default value 10 for the passive scalar The boundary conditions are defined in the user interface and depend on the boundary zone e Flow inlet Dirichlet condition an inlet velocity of 1 m s and a time dependent inlet tem perature are imposed e Outlet default value e Walls velocity pressure and thermal scalar default value passive scalar different conditions depending on the color and geometric parameters The boundary conditions for the passive scalar are identical as in case 2 as specified in the following table fwall I Dirichlet 0 wall_3 Dirichlet _0__ The wall_1 to wall_6 regions are defined as follows through color references and geometric local ization Label Color and geometric parameters wall1 24 and 0 1 lt X and X lt 0 5 wall_2 2 or 3 wall_3 4 or 7 or 21 or 22 or 23 wall_4 6andY gt 1 wall_5 6 and Y lt 1 wall_6 31 or 33 Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 29 120 Figure III 3 shows the colors used for boundary conditions and table 111 3 defines the correspondance between the colors a
65. m f90 only when it appears in the SRC directory will it be taken into account by the code Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 24 120 is used all boundary conditions have to be defined in the Graphical User Interface Only the conditions that differ from this first definition need to appear in usclim f90 The boundary conditions defined in usclim f90 will replace those specified in the Graphical Interface In this case the temperature at entry is supposed variable in time following the law III 2 T 20 100t for 0 lt t lt 3 8 T 400 for t gt 3 8 where T is the temperature in C and t is the time in s Probe 9 on the inlet 0 l 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Time s Figure 111 7 Time evolution of the temperature at inlet 3 6 Output management The output management is the same as in case 2 except that a nineth monitoring point will be added just at the entry to monitor the temperature evolution at inlet Downcomer Vessel s bottom i ama 0 2 005 2 25 0 a 005 05 0 Core s 0 05 os o s oson o 9 05 225 0 Figure 111 8 Position and coordinates of probes in the full domain In this case the Pressure the Tubulent energy and the Dissipation will be removed from the listing file Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 25 120
66. m71 C7 The water density and dynamic viscosity are variable with the temperature The functions are given below 1 4 Mesh characteristics The mesh used in the actual study had 125 000 elements It has been coarsened for this example in order for calculations to run faster The mesh used here contains 16 320 elements Type unstructured mesh Coordinates system cartesian origin on the middle of the horizontal pipe at the intersection with the nozzle Mesh generator used SIMAIL Color definition see figure IV 2 7 2 6 Figure IV 2 Colors of the boundary faces 2 CASE 5 Stratified junction In this case advanced post processing features will be used A specific pos processing sub mesh will be created containing all the cells with a temperature lower than 21 C so that it can be visualized with EnSight for instance The variable temperature will be post processed on this sub mesh A 2D clip plane will also be extracted along the symmetry plane of the domain and temperature will be written on it 2 1 Calculation options The following options are considered for the case Flow type unsteady flow Time step variable in time and uniform in space Turbulence model k e Scalar s temperature Physical properties uniform and constant for specific heat and thermal conductivity and variable for density and dynamic viscosity Pressure interpolation in stratified flows improved Code_S
67. models Calculation features Mobile mesh Turbulence models __ Thermal model _ Radiative transfers _ Conjugate heat transfer Additional scalars Definition and initialization _ Physicals properties Physical properties Reference values J Fluid properties _ Gravity hydrostatic pressure Volume conditions Volume regions definition Initialization i Head losses amp Boundary conditions MN Definition of boundary regi Boundary conditions BY Numerical parameters 6 8 Calculation control By Calculation management outlet 2 134 symmetry 3 Symmetry 8 or 9 or 28 or 29 or 38 or 39 Add Delete Add from Prepocessor listing Import groups and references from Preprocessor listing Figure V 44 Creation of the boundary zones Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 86 120 In this case different conditions are applied for the walls Separate corresponding wall boundary regions must therefore be created following the data in the following table Label Zone Nature Localization wall2 5 wall 2 or 3 wall_3 6 wall 4or 7 or 21 or 22 or 23 wall_4 7 wall 6 and Y gt 1 wall5 8 wall 6 and Y lt 1 wall_6 9 wall 31 or 33 The wall_1 region combines color and geometrical criteria The associated character string to enter in the Selection criteria box is as follows 24 and 0 1 lt X and
68. n Workstation se Cluster with PBS queue system Cluster with LSF queue system Cluster with SGE queue system Prepare batch calculation Number of processors User files a Advanced options Code_Saturne batch running fo Figure V 36 Prepare batch analysis Computer selection Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 78 120 Click on the icon to Select the batch script file to select the launch script The default launch script is named runcase and is located in the SCRIPTS directory Select it and click on Open Remember to save the Xml file before opening the launch script X Search the batch script 9 lolx 000502 Look in 3 home nicolas ETUDE CTION CASE1 SCRIPTS Mm Compu A nicolas runcase runcase help File name runcase Open Files of type All Files I di Figure V 37 Prepare batch analysis Batch script file selection Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 79 120 When the script is selected new options will appear On this calculation the number of processors used will be left to 1 When launching a calculation a temporary directory is created on the machine where the script copies and creates temporary files and from where the Code_Saturne executable is launched Should some user ro
69. nd the type of boundary condition to use o ro y e 2346 7 21 22 23 3133 Wal 24 for 0 1 lt X lt 0 5 8 9 28 29 38 39 Table 111 3 Boundary faces colors and associated references 4 3 Variable Density In this case the density is a function of temperature the variation law is defined in the Graphical User Interface although it can also be defined in a Fortran user routine The expression is p T AT B C III 3 where p is the density T is the temperature A 4 0668 x 1073 B 5 0754 x 107 and C 1000 9 In order for the variable density to have an effect on the flow gravity must be set to a non zero value y 9 81le will be specified in the Graphical Interface 4 4 Head loss To simulate the presence of an obstacle 0 20 m large and 0 5 m high in the vessel a zone of head loss will be created in the domain fig 111 12 The head loss zone is located between the coordinates X 0 2 m and X 0 4 m and Y 0 75 m and Y 0 25 m The head loss coefficient to apply is 10 and is isotropic Ree Lower core RE plate and Ve x e XXX cor Domain MES Dn Frame Obstacle Figure III 12 Full domain geometry with the obstacle Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 30 120 4 5 Parameters All the parameters necessary to this study can be defined through the Graphical Interface However the calculation of the spatial average is defin
70. nditions Definition of boundary regi En Boundary conditions Boundary conditions cold_inlet 1 inlet 2 hot_inlet inlet 6 outlet outlet 7 Velocity Direction specified coordinates gt x oo fly Numerical parameters E fly Calculation control BF Calculation management o E mi E E E Figure V 74 Identity and paths Calculation environment Meshes selection Mesh quality criteria Thermophysical models Calculation features Mobile mesh Turbulence models Thermal model Radiative transfers Conjugate heat transfer ditional scalars Definition and initialization Physicals properties ysical properties Reference values Fluid properties Gravity hydrostatic pressure lume conditions Volume regions definition Initialization Head losses undary conditions Definition of boundary regi Boundary conditions Numerical parameters Calculation control Calculation management Al EU Figure V 75 Calculation by hydraulic diameter v Hydraulic diameter m Dynamic boundary conditions Boundary conditions A wall wall cold_inlet inlet hot_inlet 2 inlet outlet Velocity Direction specified coordinates D Turbulence Calculation by hydraulic diameter D Hydraulic diameter m Dynamic boundary conditions Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 116 120 For the scalar boundary cond
71. nitialization i Physicals properties 10 x Number of interior faces 3 Boundary conditions SS Definition of boundary re 0 x Number of boundary faces Boundary conditions S Numerical parameters i Steady flow management J Equation parameters Global parameters Calculation control Output control Volume solution control Profiles Calculation management no Memory management fa Start Restart z 4 Prepare batch calculation 4 Lal Ea A Figure V 34 User arrays Code_Saturne documentation EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial Page 76 120 The item Start Restart allows to start a new calculation from the results of a former one It is not the case in the present calculation so nothing has to be modified ex _ Calculation features a J Mobile mesh Turbulence models _ Thermal model i Radiative transfers DI a Physical properties Reference values Fluid properties Gravity hydrostatic pres 5 B Volume conditions J Volume regions definition d Initialization Additional scalars Definition and initialization Physicals properties Boundary conditions Definition of boundary re J Boundary conditions Numerical parameters Steady flow management Equation parameters Global parameters Calculation control _ Output control Volume solution control _ Profiles 3 Calculation manag
72. ntrol to change the frequency for the printing of information in the output listing The options are e No output e Output listing at each time step e Output at each n time step the value of n must then be specified Here and in most cases the second option should be chosen Ex Output Control Monitoring Points Coordinates Identity and paths al Calculation environment Meshes selection A quality criteria Output listing at each time step Thermophysical models LT Calculation features Output listing at each time step nn Mobile mesh Output every n time steps 1 Turbulence models Only at the end of calculation Thermal model Radiative transfers Fluid domain post processing x O awl shii Domain boundary post processing i y Reference values Enr gt Gravity hydrostatic pressure Type of post processing for mesh fixed D S E Volume conditions Post processing format EnSight Gold 7 i y Volume regions definition Initialization Options Additional scalars 7 Definition and initialization format binay 7 kraai E J Physicals properties polygons display D Boundary conditions _ J Definition of boundary regi polyhedra display y Boundary conditions big endian _ amp Numerical parameters i Re NN amp Calculation control Outputs listings E Output control iL Volume
73. o case 1 Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 83 120 To add an additional scalar click on the Definition and Initialization item under the Additional scalars heading The characteristics of the thermal scalar are still the same Its initial value is 20 C and it can vary between 0 C and 400 C To create an additional scalar click on Add then enter e its Name scalar_2 e its Initial value 10 e its Minimal value 0 e its Maximal value 400 Select volume region for initial value _ Identity and paths Calculation environment Region selection J Meshes selection __ Mesh quality criteria hermophysical models J Calculation features Initial Variance Minimal Maximal _ Mobile mesh value of scalar value value 2 Thermal model Tempe 20 no o ao O User scalar definition Thermal model Tempc Radiative transfers Conjugate heat transfer A EIT properties LA Physical properties B Calculation control B Calculation management Figure V 42 Additional scalar User scalar definition Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 84 120 In the item Physical properties still under the heading Additional scalars specify the diffusion coeffi cient of this new scalar Click on the scalar name to highlight it then enter the value in the box In this case the value is 0 895 x 1074 m s 1 es
74. on environment Gravity along X 0 0 m s Meshes selection A Mesh quality criteria Gravity along Y 0 0 m s 3 Thermophysical models gt i Calculation features Gravity along Z 0 0 ae Mobile mesh Turbulence models Hydrostatic pressure Thermal model i Radiative transfers Physical properties standard i Reference values Fluid properties MN Gravity hydrostatic pressure 6 Volume conditions LL Volume regions definition i Initialization 3 Additional scalars Definition and initialization eL Physicals properties E a Boundary conditions BF Numerical parameters E fly Calculation control fi Calculation management Pressure interpolation in stratified flow improved l LT Figure V 16 Physical properties gravity hydrostatic pressure Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 59 120 Boundary conditions now need to be defined Go to the item Define boundary regions under the heading Boundary conditions The following window opens fig V 19 Sx o Definition of boundary regions _ Identity and paths Calculation environment Label Zone Nature J Meshes selection abel Zope sce Mesh quality criteria Thermophysical models J Calculation features Mobile mesh Turbulence models Thermal model Radiative transfers 3 Physical properties
75. ons definition Initialization S Additional scalars A Definition and initialization iL Physicals properties l o l ne Boundary conditions Add Delete B Numerical parameters ee B Calculation control Be Calculation management User scalar definition Initial Variance value of scalar Minimal Maximal value value Temp c 20 no le 12 le 12 E E ne e ago Figure V 12 Initialization of the scalar Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 55 120 Click on the thermal scalar in the list to change e its name e its initial value e its minimal value e its maximal value In this case the temperature can vary between 0 C and 400 C After entering the new values click on Modify in order to validate these changes Select volume region for initial value Region selection all cells n User scalar definition Initial Variance Minimal Maximal value of scalar value value Amp E 20 no O 400 Figure V 13 Initialization of the scalar The item Physicals properties under the heading Additional scalars is used to specify the physical properties of the additional scalars i e those that are not the thermal scalar In this case there are no additional scalars the item is therefore unused 2not to be confused with the heading Physical properties in the main list Code_Saturne
76. pied in the RESU directory More precisely a RES_USERS xxxxxxxx directory will be created in the RESU folder in which the file will be copied Click on the icon New user result file to enter the associated dialog window Enter the file name moy dat in the field Further file names can be added When finished click on OK Data files selection Results tiles MH _ moy dat O Cancel Figure V 66 User results files Code_Saturne documentation Page 109 120 EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial SOLUTION FOR CASE 5 The preparation of the calculation for case 5 is very similar to the other cases e Open the Code_Saturne interface e Open a new case Check the name of the mesh Select a k model Use a thermal scalar in Celsius degrees In the item nitialization under the heading Volume conditions set the initial value of the temperature in the domain to 38 5 C Initialize the turbulence with the reference velocity 0 03183 m s E x SR Select volume zone for initialization Identity and paths gt eg Calculation environment Volume zone all_cells Meshes selection A y Mesh quality criteria Velocity initialization E Thermophysical models Calculation features vitesseX 0 0 m s Mobile mesh Turbulence models VitesseY 0 0 m s Thermal model
77. ption any later version Code_Saturne is distributed in the hope that it will be useful but WITHOUT ANY WARRANTY without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE See the GNU General Public License for more details Part Il SIMPLE JUNCTION TEST CASE Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 10 120 1 1 1 Objective General description The aim of this case is to train the user of Code_Saturne on an oversimplified 2D junction including an inlet an outlet walls and symmetries 1 2 Description of the configuration The configuration is two dimensional It consists of a simple junction as shown on figure 11 1 The flow enters through a hot inlet into a cold environment and exits as indicated on the same figure This geometry can be considered as a very rough approximation of the cold branch and the downcomer of the vessel in a nuclear pressurized water reactor The effect of temperature on the fluid density is not taken into account in this first example 0 1 oon infet y Downcomer Frame outlet Figure II 1 Geometry of the downcomer 1 3 Characteristics Characteristics of the geometry and the flow Height of downcomer H 3 00 m Thickness of downcomer Ea 0 10 m Diameter of the cold branch Dy 0 50 m Inlet velocity of fluid Physical characteristics of fluid The initial water temperature in the domain is equal to 20
78. r output management CASE 3 same as case 2 with time dependent boundary conditions fluid density depending on the temperature and calculation restart CASE 4 same as case 3 with head loss parallelism and spatial average 2 CASE 2 Passive scalar with various boundary conditions and output management 2 1 Calculation options Some options are similar to case 1 Turbulence model k e Scalar s 1 temperature Physical properties uniform and constant The new options are Flow type unsteady flow Time step uniform and constant Scalar s 2 passive scalar Management of monitoring points 2 2 Initial and boundary conditions Initialization 20 C for temperature default value 10 for the passive scalar could correspond to a tracer concentration for instance Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 18 120 The boundary conditions are defined in the user interface and depend on the boundary zone 1 e Flow inlet Dirichlet condition an inlet velocity of 1 m s an inlet value of 200 for the passive scalar are imposed an inlet temperature of 300 C and e Outlet default value e Walls velocity pressure and thermal scalar default value passive scalar different conditions depending on the color and geometric parameters In order to test the ability to specify boundary condition regions in the Graphical Interface var
79. r results file The head loss is defined in the Graphical User Interface Go to Volume regions definition under the heading Volume conditions Click on Add unselect Initialization and select Head losses in the box named Nature In the box named Label name the head loss region Define the limits of the head losses region in Selection criteria The associated character string to enter is as foolows 0 2 lt X and 0 4 gt X and 0 75 lt Y and 0 25 gt Y File Tools Window Help GE Seb E es e Definition of volume regions Identity and paths E E Calculation environment Label Zone Nature Selection criteria Meshes selection Mesh quality criteria all_cells 1 Initialization all a Thermophysical models m z PESA 0 2 lt x and 0 4 gt x and 0 75 lt y i Q Calculation features a aUan ui 0 25 gt y Mobile mesh 4 Turbulence models L Thermal model Radiative transfers CF Conjugate heat transfer G By Additional scalars E B Physical properties Volume conditions Add Delete ma volume regions definition Initialization Add from Prepocessor listing _ Head losses ES Boundary conditions Import groups and references from Preprocessor listing Definition of boundary regi Boundary conditions a Numerical parameters amp B Calculation control BF Calculation management dl a E O Definit
80. rameters and change the number of itera tions It must be the total number of iterations from the beginning of the first calculation The first calculation was done with 300 iterations and another 400 iterations are needed for the present case Therefore the value 700 must be entered a x Identity and paths BF Calculation environment H B Thermophysical models E Additional scalars amp Physical properties Reference values Fluid properties i Gravity hydrostatic pressure E B Volume conditions y Boundary conditions d Numerical parameters MN Time step iL Equation parameters Global parameters Calculation control 3 3 Calculation management h d User arrays n i Memory management Start Restart Prepare batch calculation al f gt Unsteady flow algorithm management Time step option Uniform and constant Reference time step Number of iterations restart included 1700 0 05 Time step limitation with the local thermal time step Option zero time step _ Eventually run the calculation Figure V 62 Time step Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 105 120 4 SOLUTION FOR CASE 4 This case is similar to case 3 with the following differences e parallel computation on 2 processors e head loss e calculation of a spatial average e dealing with a use
81. roperties i Reference values Fluid properties i Gravity hydrostatic pressure Volume conditions Volume regions definition Initialization Head losses 5 Boundary conditions Definition of boundary regi Boundary conditions E Numerical parameters MY Time step a 4 Equation parameters i h J Global parameters By Calculation control d Deo Calculation management Option zero time step g Ti Figure V 50 Time step setting No change is needed in the Equation parameters and Global parameters items Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 92 120 Go to the item Output control under the heading Calculation control to set the output parameters Keep the default value for the output listing frequency For the Post processing select the third option output every n time steps and set the value of n to 2 Activate the post processing on the boundary faces by ticking the Domain boundary post processing box The EnSight format file will contain an additional part composed of the boundary faces on which boundary conditions and some other variables can be visualized This allows to check if the boundary conditions for the passive scalar have been properly set Ex Output Control Monitoring Points Coordinates J Identity and paths 3 3 Calculation environment Meshes selection i
82. s Calculation environment Time step option Variable in time and uniform in space D ea Thermophysical models 3 Additional scalars Physical properties Volume conditions Number of iterations restart included Volume regions definition Initialization Head losses Boundary conditions Definition of boundary regi Boundary conditions Minimal time step factor 0 01 amp j Numerical parameters MN Time step Maximal time step factor 70 0 h Equation parameters LA Time step maximal variation 0 1 i Ly Global parameters E BF Calculation control B Calculation management Time step limitation with Pi the local thermal time step Reference time step Maximal CFL number Maximal Fourier number Option zero time step Figure V 78 Time step Set the frequency of post processing files to 10 Create four monitoring probes at the following coordinates Output Control Monitoring Points Coordinates _y Identity and paths a Calculation environment Outputs listings sa Thermophysical models pbb ah dui Output listing at each time step z sa Physical properties Volume conditions Post processing Volume regions definition Bri Initialization Post processing every n time steps 0 Head losses Boundary conditions Fluid domain post processing x Definition of boundary regi Boundary conditions gt E Numerical parameters a J Time step Typ
83. solution control J Profiles f Calculation management dad CD Figure V 30 Output control output listing Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 72 120 For the post processing by default EnSight format files there are three options e Only at the end of calculation e At each time step e Post processing every n time steps In this case we are interested in the evolution of the variables during the calculation so the second option is chosen Output Control Monitoring Points Coordinates Outputs listings Output listing at each time step D Fost processing Only at the end of calculation Only at the end of calculation At each time step Past processing every n time steps La Type of post processing for mesh fixed we Post processing format Ensight Gold Options polygons display polyhedra display big endian Figure V 31 Output control post processing Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 73 120 The other options are kept to their default value Output Control Monitoring Points Coordinates Outputs listings Output listing at each time step M Post processing Fluid domain post processing EJ Domain boundary post processing Type of post processing for mesh fixed we Post processing format EnSight Gold v Options polygons display RA polyhedra display 7 big
84. ssing features copy the three routines usdpst f90 usmpst f90 and usupst f90 in the SRC directory The general content of these routines is described in the user manual or in the examples available in the directory SRC REFERENCE base The modified routines adapted to this test case are available in the examples directory Only the main elements are mentionned here e usdpst f90 This routine is called only once at the beginning of the calculation It allows to define the different writers and parts The first writer is the standard writer which creates the directory CHR ENSIGHT xxxxxxxx It is created by default and has the number 1 Set the number of additional writers NBCAS to 1 For the first and unique additional writer specify the following elements e NOMCAS chr prefix of the EnSight files e NOMREP Tinf21 ensight name of the directory e NOMFMT EnSight Gold format of the post processing e OPTFMT binary format options here binary files A e INDMOD 2 indicates that the parts in this writer will be time dependent in its content e NTCHRL 5 periodicity of output directory TINF21 ENSIGHT xxxxxxxx will be created with the post processing results associated to this writer Set the number of additional parts NBPART to 2 For each part set the number of cells internal faces and boundary faces respectively NLCEL NLFAC NLFBR and the lists LSTCEL LSTFAC and LSTFBR of the elements in the par
85. st processing and the chronological records This default choice can be modified by the user In this case the Pressure the Tubulent energy and the Dissipation will be removed from the listing file The Courant number CFL and Fourier number will be removed from the post processing results Eventually probes will be defined for chronological records following the data given in figure III 4 Then the total pressure will be deactivated for all probes and the Velocity U will only be activated on probes 1 2 6 7 and 8 l Downcomer Z Vessel s bottom i 2 0 05 2 25 0 3 0 05 2 75 0 s 0 05 025 0 6 075 025 0 ar 075 025 0 8 0r on 0 FAP A A Figure 111 4 Position and coordinates of probes in the full domain In addition the domain boundary will be post processed This allows to check the boundary conditions and especially that of the passive scalar 2 5 Results Figure 111 5 shows the boundary domain colored by the passive scalar boundary conditions The different regions of boundary conditions defined earlier can be checked Figure III 6 presents results obtained at different times of the calculation They were plotted from the post processing files with EnSight 3this can be very useful to save some disk space if some variables are of no interest as post processing files can be large Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page
86. start MN Prepare batch calculation al OC lalo User files gt Figure V 38 Prepare batch analysis Execution Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 80 120 2 SOLUTION FOR CASE 2 This case corresponds to a new study in which there will be three calculation cases cases 2 3 and 4 All of them can be created in a single code_saturne create command or additional cases can be added later To test both possibilities first create the study directory with cases directories CAS2 and CAS4 code_saturne create s FULL DOMAIN CAS2 CAS4 then go in the study directory and add the CAS3 directory cd FULL DOMAIN code_saturne create c CAS3 Go to the DATA directory in CAS2 open a new case and select the meshes to use Click on the heading Calculation environment then on the item Meshes selection In this case the three meshes have to be joined So don t delete any mesh and activate the Join meshes option by clicking in the box Additional information appears on the page If it is left untouched the Code_Saturne Preprocessor will test all the boundary faces for potential joining based on geometrical criteria To make mesh joining more efficient this analysis can be restricted to a sub set of boundary faces This is the case in the present calculation since only faces of colors 5 24 and 34 are liable to be joined Click on t
87. t The first part the clip plane will be created by detecting the internal faces which have a center of gravity CDGFAC between 0 01 and 0 01 The second part the cells where the temperature is lower than 21 C will be specified in usmpst f90 Yet it must be initialized in usdpst f90 The easiest is to set NLCEL NCEL total number of cells and when doing so there is no need to specify the LSTCEL array Eventually the different parts must be associated with the different writers through the PSTASS routine Part 1 is associated to the writer 1 and part 2 to the writer 1 parts can only contain similar elements i e combinations of internal and boundary faces are allowed but combina tions of cells and faces are not Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 120 120 e usmpst f90 This routine is called at each time step It allows to redefine the content of certain parts using any variable especially the temperature for this case Only part 2 is concerned A DO ENDDO loop on all the cells allows to identify those where the temperature is lower than 21 C and hence calculate the number of cells NCELPS in the part and the list of cells LSTCEL e usvpst f90 This routine is called at each time step It allows to specify which variable will be written on which part The writing in the post processing files is triggered by the routine PSTEVA that must be called for each
88. t by the code Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 31 120 average in the present case The spatial average of the temperature will be calculated at each time step and the result wrote in a file named moy dat The values are saved in order to draw the time evolution of the average temperature Beware when calculating the average Since the calculation is running in parallel computing the sum of the temperatures on all the cells will only yield for each processor the sum on the cells managed by this processor In order to obtain the full sum the parallelism routine PARSOM must be used see example Note usproj f90 contains many examples They should be removed before running the case 4 7 Output management The output management is the same as in case 3 5 Downcomer a ZZ ves oo s 0 05 2 75 0 Y N 8 6 o o 0_ 7 osos o 7 8 075 075 A o 05 225 0_ VILILILIDIFIFIFIFIFILIFILIPIZIZI AAA Figure 111 14 Position and coordinates of probes in the full domain In this case the Pressure the Tubulent energy and the Dissipation will be removed from the listing file The Courant number CFL and Fourier number will be removed from the post processing results Eventually probes will be defined for chronological records following the data given in figure III 4 Then the total pressure w
89. ted by or or and keywords Definition of boundary regions A Label inlet a pa from Prepocessor listing Import groups and references from Preprocessor listing Figure V 18 Creation of a boundary region Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 61 120 The specification of the inlet condition is detailled in the following pages The settings will be as follows e Label inlet e Zone 1 e Nature inlet e Localization 1 Type all the information in the fields the result diplays as figure V 19 Definition of boundary regions _ _ _ _ __ Label Zone Nature Selection criteria Add Add from Prepocessor listing Import groups and references from Preprocessor listing Figure V 19 Creation of a boundary region Remember to save the Xml file regularly Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 62 120 Do the same thing for the other boundaries In our case colors 8 and 9 are symmetry boundaries One option can be to define a separate zone for each color as follows Label symmetry_l symmetry_2 Zone 3 4 Nature symmetry symmetry Localization 8 9 But it is usually faster to regroup the different colors in one single zone as shown on figure V 20 In our case the localization for this zone is the string 8 or 9 Definition of boundary regions
90. tialization Physicals properties Physical properties Reference values Fluid properties Gravity hydrostatic pressure Volume conditions Volume regions definition Initialization Head losses Boundary conditions Definition of boundary regi MS Boundary conditions Selection criteria inlet J 5 20r3or4or6o0r7 Velocity Direction normal direction to the inlet Turbulence Calculation by hydraulic diameter 7 B Numerical parameters E fly Calculation control o Bg Calculation management Hydraulic diameter 1 0 m Scalars Exchange Scalar Name Coefficient Tempc Prescribed Vo 0 Figure V 22 Dynamic variables boundary conditions Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 65 120 Click on the label inlet In the section Velocity select norm then in the sub section Direction choose specified ccordinates and enter the normal vector components of the inlet velocity For the turbulence chose the inlet condition based on a hydraulic diameter and specify it eX lm eY_ 0m eZ 0m e D 0 5 m Boundary conditions ET 1l inlet 5 2or3or4dor6 or Direction specified coordinates 7 Turbulence Calculation by hydraulic diameter 7 Hydraulic diameter m Figure V 23 Dynamic variables boundary inlet Code_Saturne EDF R amp D Code_Saturne vers
91. time dependent variables and Fortran user routines 1 2 Description of the configuration The fluid domain is composed of three separate meshes very roughly representing elements of a nuclear pressurized water reactor vessel e the downcomer e the vessel s bottom e the lower core plate and core Figure II 1 represents the complete domain The flow circulates from the top left horizontal junction to the right vertical outlet SEN Pe Lower core ELI plate and core KRRRA A AA ALI a Domain Frame Dimensions Figure III 1 Geometry of the complete domain 1 3 Characteristics The characteristics of the geometry and the flow are Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 16 120 Physical characteristics of fluid The initial water temperature in the domain is equal to 20 C The inlet temperature of water in the cold branch is 300 C Water characteristics are considered constant and their values taken at 300 C and 150 x 10 Pa except density which is considered variable in cases 3 and 4 e density p 725 735 kg m73 e dynamic viscosity u 0 895 x 1074 kg m7 s7 e heat capacity Cp 5 483 kg PO Thermal Conductivity 0 02495 W m7 K7 1 4 Mesh characteristics Figure 111 2 shows a global view of the mesh and some details of the joining zones to show that Code_Saturne can deal with hanging nodes
92. tutorial documentation Page 48 120 The item Analysis features under the heading Thermophysical environment allows to define the type of flow to be simulated In this case a steady flow will be chosen EJES Steady Unsteady flow algorithm Identity and paths Calculation environment steady flow n Meshes selection steady flow Mesh quality criteria unsteady flow H Thermophysical models MN Calculation features single Phase Flow Mobile mesh nr Turbulence models Thermal model BY Physical properties B Volume conditions B Additional scalars B Boundary conditions BB Numerical parameters ff off B Calculation control Calculation management Atmospheric flows Gas combustion fF E F E Pulverized coal combustion Figure V 6 Flow type Code_Saturne EDF R amp D Code_Saturne version 2 0 0 rc1 tutorial documentation Page 49 120 The turbulence model is selected in the following list e laminar flow no model e mixing length o k c e k Linear Production e Rij e LLR e Rij e SSG e v2f p model e kw SST e LES Smagorinsky e LES dynamic model Turbulence Model _ _ gt E gt gt gt gt gt 2 R _ _ Identity and paths e 3 Calculation environment i Meshes selection k epsilon ai Mesh quality criteria No model i e laminar flow Thermophysical models Mixing length ne i Calculation features Mobile mesh k epsi
93. utines read or write case specific files they must be copied in the temporary directory or from the temporary directory into the RESU directory The User files icon allows the user to specify user data files in the DATA directory or user result files that will then be copied automatically to or from the temporary directory In this example no user file is needed Finally the Advanced options icon allows to change some more advanced parameters that will not be needed in this simple case Eventually save the Xml file and execute it by clicking on Code_Saturnebatch running The results will be copied in the RESU directory ch ES Computer selection 2 _ Calculation features al _ Mobile mesh Workstation Turbulence models e Select the batch script file QD Physical properties 4 Reference values Prepare batch calculation Fluid properties _ Gravity hydrostatic pres Number of processors 1 E 3 Volume conditions Volume regions definition _ Initialization 6 Additional scalars Advanced options L Definition and initialization J Physicals properties 6 Boundary conditions Code_Saturne Definition of boundary re batch J Boundary conditions running S E Numerical parameters _ Steady flow management _ Equation parameters _ Global parameters 6 Calculation control J Output control Volume solution control _ Profiles eg Calculation management User arrays Memory management Start Re

CodeSaturne tutorial

Contents

Download Pdf Manuals

Related Search

Related Contents