Initial commit from Polytechnique Montreal

1 files changed, 238 insertions, 0 deletions

diff --git a/Donjon/data/thm_salt_lf_trans.x2m b/Donjon/data/thm_salt_lf_trans.x2m
new file mode 100644
index 0000000..e65b777
--- /dev/null
+++ b/Donjon/data/thm_salt_lf_trans.x2m
@@ -0,0 +1,238 @@
+*=====================================================================*
+* FILE:           thm_salt_lf_trans.x2m
+* TYPE:           DONJON input file
+* AUTHOR:         Cristian Garrido Tamm (cristian.garrido@idom.com)
+* PURPOSE:        Single pincell surrounded by FLiBe as coolant to
+*                 test the Molten Salt thermalhydraulics in THM
+* CREATION DATE:  Wednesday Aug 16, 2023 21:13:01 CEST
+* LAST MODIFIED:  Tuesday Jul 30, 2024 19:12:05 CEST
+*=====================================================================*
+*---
+* DEFINE MODULES
+*---
+MODULE      GEO: RESINI: TRIVAT: TRIVAA: NCR: FLUD: FLPOW: MACINI:
+            USPLIT: TINST: UTL: DELETE: ABORT: THM: END: ;
+*---
+* DEFINE PROCEDURES
+*---
+PROCEDURE   assertS ;
+*---
+* DEFINE DATA STRUCTURES
+*---
+LINKED_LIST Geom Track Flux Cpo Matex Lib MicroF MacroF Macro2 Fmap
+            System Power Burnup Thm DONNEES OLCM ;
+*---
+* DEFINE VARIABLES
+*---
+* Constants
+STRING DOUBLE_LINE := 
+"====================================================================" ;
+STRING SINGLE_LINE := 
+"--------------------------------------------------------------------" ;
+LOGICAL True False := $True_L $False_L ;
+STRING FILENAME := 'thm_salt_of' ;
+INTEGER MaxR := 10000 ;
+REAL Pi := 3.14159 ;
+* Input data
+REAL maxh := 100.0 ; ! cm
+REAL Cote := 1.26  ; ! cm
+REAL Rfuel := 4.55063E-1 ; ! cm
+REAL Rclad := 5.00184E-1 ; ! cm
+REAL z1 := maxh 20.0 / ;
+REAL z2 := maxh 20.0 / z1 + ;
+REAL z3 := maxh 20.0 / z2 + ;
+REAL z4 := maxh 20.0 / z3 + ;
+REAL z5 := maxh 20.0 / z4 + ;
+REAL z6 := maxh 20.0 / z5 + ;
+REAL z7 := maxh 20.0 / z6 + ;
+REAL z8 := maxh 20.0 / z7 + ;
+REAL z9 := maxh 20.0 / z8 + ;
+REAL z10 := maxh 20.0 / z9 + ;
+REAL z11 := maxh 20.0 / z10 + ;
+REAL z12 := maxh 20.0 / z11 + ;
+REAL z13 := maxh 20.0 / z12 + ;
+REAL z14 := maxh 20.0 / z13 + ;
+REAL z15 := maxh 20.0 / z14 + ;
+REAL z16 := maxh 20.0 / z15 + ;
+REAL z17 := maxh 20.0 / z16 + ;
+REAL z18 := maxh 20.0 / z17 + ;
+REAL z19 := maxh 20.0 / z18 + ;
+REAL z20 := maxh 20.0 / z19 + ;
+REAL Fuelpwr := 100.0 ; ! kW/l or W/cm3
+REAL Ucool := 3.0 ; ! m/s
+*REAL Mass := 4.655149E-3 ; ! kg
+REAL Tfuel := 800.15 ; (*fuel temperature (K)*)
+REAL Tcool :=  800.15 ; (*coolant temperature (K)*)
+* Options
+* Variables
+INTEGER istep := 0 ; ! step counter
+REAL t tend dt := 0.0 1.0 0.1 ; ! time, end time, time step (s)
+REAL Toutlet ;
+
+*---
+* BEGIN CALCULATION
+*---
+ECHO DOUBLE_LINE ;
+ECHO 'BEGIN INPUT ' FILENAME ;
+ECHO DOUBLE_LINE ;
+
+REAL Volume := Pi Rfuel * Rfuel * maxh * ; ! cm3
+ECHO "Volume =" Volume "cm3" ;
+REAL powi := Fuelpwr Volume * 1.E6 / ;
+* EVALUATE powi := powi 1E3 * 20.0 / ;
+ECHO "total reactor power=" powi "MW" ;
+
+Geom := GEO: :: CAR3D 1 1 20
+   X- REFL X+ REFL    Y- REFL Y+ REFL    Z- REFL Z+ REFL
+   MESHX 0.0 <<Cote>>
+   MESHY 0.0 <<Cote>>
+   MESHZ  0.0 5.000 10.000 15.000 20.000 25.000 30.000 35.000 40.000
+ 45.000 50.000 55.000 60.000 65.000 70.000 75.000 80.000 85.000 90.000
+ 95.000 100.000
+   MIX
+   PLANE 1
+      1
+   PLANE 2 SAME 1
+   PLANE 3 SAME 1
+   PLANE 4 SAME 1
+   PLANE 5 SAME 1
+   PLANE 6 SAME 1
+   PLANE 7 SAME 1
+   PLANE 8 SAME 1
+   PLANE 9 SAME 1
+   PLANE 10 SAME 1
+   PLANE 11 SAME 1
+   PLANE 12 SAME 1
+   PLANE 13 SAME 1
+   PLANE 14 SAME 1
+   PLANE 15 SAME 1
+   PLANE 16 SAME 1
+   PLANE 17 SAME 1
+   PLANE 18 SAME 1
+   PLANE 19 SAME 1
+   PLANE 20 SAME 1
+
+;
+
+Geom Matex := USPLIT: Geom :: NGRP 2 MAXR <<MaxR>>
+               NFUEL 1  FMIX  1
+            ;
+
+*--
+* Fuel map definition
+*--
+Fmap Matex := RESINI: Matex ::
+      ::: GEO: CAR3D 1 1 20
+                EDIT  0
+                X- REFL X+ REFL    Y- REFL Y+ REFL    Z- REFL Z+ REFL
+   MESHX 0.0 <<Cote>>
+   MESHY 0.0 <<Cote>>
+   MESHZ  0.0 5.000 10.000 15.000 20.000 25.000 30.000 35.000 40.000
+ 45.000 50.000 55.000 60.000 65.000 70.000 75.000 80.000 85.000 90.000
+ 95.000 100.000
+   MIX
+   PLANE 1
+      1
+   PLANE 2 SAME 1
+   PLANE 3 SAME 1
+   PLANE 4 SAME 1
+   PLANE 5 SAME 1
+   PLANE 6 SAME 1
+   PLANE 7 SAME 1
+   PLANE 8 SAME 1
+   PLANE 9 SAME 1
+   PLANE 10 SAME 1
+   PLANE 11 SAME 1
+   PLANE 12 SAME 1
+   PLANE 13 SAME 1
+   PLANE 14 SAME 1
+   PLANE 15 SAME 1
+   PLANE 16 SAME 1
+   PLANE 17 SAME 1
+   PLANE 18 SAME 1
+   PLANE 19 SAME 1
+   PLANE 20 SAME 1
+
+;
+!
+NXNAME '01' NYNAME  'A'
+NCOMB 1
+B-ZONE 1
+
+ADD-PARAM PNAME 'T-FUEL' PARKEY 'TFA' GLOBAL
+ADD-PARAM PNAME 'T-COOL' PARKEY 'TCA' GLOBAL
+BTYPE INST-BURN
+INST-BVAL CHAN 0.0
+REACTOR-POW <<powi>> 
+AXIAL-PFORM  1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0
+ 1.0 1.0 1.0 1.0 1.0
+*BUNDLE-POW SAME <<powi>>
+SET-PARAM 'T-FUEL' <<Tfuel>>
+SET-PARAM 'T-COOL' <<Tcool>>
+*FUEL WEIGHT <<Mass>>
+;
+
+* UTL: Fmap :: STEP UP PARAM STEP AT 4 DIR IMPR P-NAME * ;
+*--
+* THM single-stage calculation
+*--
+Thm Fmap := THM: Fmap ::
+    EDIT 5
+    FLUID SALT LiF-BeF2 0.66-0.34 
+    FUEL SALT LiF-NaF-KF 0.465-0.115-0.42
+    FPUISS 1.0
+    INLET 101325.0 800.15
+    CRITFL 9E12
+    SPEED <<Ucool>>
+    ASSMB 1 0
+    RADIUS 4.55063E-3 4.55063E-3 5.00184E-3 5.00184E-3 (* m *)
+    RODMESH 15 20
+    MONO
+    POWER-LAW 1E3 12
+    0.0 0.5
+    0.1 1.0
+    0.2 2.0
+    0.3 3.0
+    0.4 4.0
+    0.5 5.0
+    0.6 5.1
+    0.7 5.2
+    0.8 5.3
+    0.9 5.4
+    1.0 5.5
+    2.0 5.5
+;
+
+WHILE t tend <= DO
+  EVALUATE istep := istep 1 + ;
+  ECHO "istep=" istep "t=" t "dt=" dt ;
+
+  Thm Fmap := THM: Thm Fmap ::
+    EDIT 1
+    FLUID SALT LiF-BeF2 0.66-0.34 
+    FUEL SALT LiF-NaF-KF 0.465-0.115-0.42
+    TIME <<t>> <<dt>>
+  ;
+  EVALUATE t := t dt + ;
+ENDWHILE ;
+
+*--
+* Dump THM object
+*--
+UTL: Thm :: DIR DUMP ;
+UTL: Fmap :: DIR DUMP ;
+REAL tc ;
+OLCM := Fmap :: STEP UP PARAM STEP AT 2 ;
+assertS OLCM :: "P-VALUE" 5 8.071872E2 ;
+assertS OLCM :: "P-VALUE" 10 8.144556E2 ;
+
+*---
+* END CALCULATION
+*---
+ECHO DOUBLE_LINE ;
+ECHO 'END INPUT ' FILENAME ;
+ECHO DOUBLE_LINE ;
+ECHO "test thm_salt_lf_trans completed" ;
+
+END: ;
+QUIT .