*=====================================================================*
* FILE:           thm_salt_of.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 18:42:46 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 + ;
*CGT
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

*---
* 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 1
    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
;

*--
* Dump THM object
*--
UTL: Thm :: DIR DUMP ;
UTL: Fmap :: DIR DUMP ;
OLCM := Fmap :: STEP UP PARAM STEP AT 2 ;
assertS OLCM :: "P-VALUE" 5 8.013987E2 ;
assertS OLCM :: "P-VALUE" 10 8.027861E2 ;

*---
* END CALCULATION
*---
ECHO DOUBLE_LINE ;
ECHO 'END INPUT ' FILENAME ;
ECHO DOUBLE_LINE ;
ECHO "test thm_salt_lf completed" ;

END: ;
QUIT .