Initial commit from Polytechnique Montreal

2 files changed, 285 insertions, 0 deletions

diff --git a/Donjon/data/pulseTHM_0d_proc/PfmapPULSE.c2m b/Donjon/data/pulseTHM_0d_proc/PfmapPULSE.c2m
new file mode 100644
index 0000000..7acd277
--- /dev/null
+++ b/Donjon/data/pulseTHM_0d_proc/PfmapPULSE.c2m
@@ -0,0 +1,113 @@
+****************************************************************
+*                                                              *
+* Procedure :  PfmapPULSE.c2m                                  *
+* Purpose   :  Fuel-map definition for single pincell in a     *
+*              pulse reactor                                   *
+* Author    :  A. Hebert (2017)                                *
+*                                                              *
+* CALL      :  FMAP MATEX := PfmapPULSE MATEX                  *
+*                     :: <<pow0>>  ;                           *
+*                                                              *
+****************************************************************
+
+PARAMETER  FMAP MATEX :: ::: LINKED_LIST FMAP MATEX ; ;
+MODULE     END: RESINI: ;
+REAL pow0 ;
+  :: >>pow0<< ;
+
+REAL maxh := 90.0 ;
+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 + ;
+
+ECHO 'z20=' z20 ;
+
+  FMAP MATEX := RESINI: MATEX :: EDIT 1
+
+          ::: GEO: CAR3D 1 1 20
+  EDIT 0
+  X- REFL   X+ REFL
+  Y- REFL   Y+ REFL
+  Z- VOID   Z+ VOID
+  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
+  MESHX 0.0 1.2584422 SPLITX 2
+  MESHY 0.0 1.2584422 SPLITY 2
+  MESHZ 0.0 <<z1>> <<z2>> <<z3>> <<z4>> <<z5>> <<z6>> <<z7>> <<z8>>
+        <<z9>> <<z10>> <<z11>> <<z12>> <<z13>> <<z14>> <<z15>> <<z16>>
+        <<z17>> <<z18>> <<z19>> <<z20>>
+  ;
+  NXNAME  '1'
+
+  NYNAME  'A'
+  NCOMB ALL
+  ADD-PARAM PNAME 'T-FUEL' PARKEY 'TF' GLOBAL
+  ADD-PARAM PNAME 'T-COOL' PARKEY 'TC' GLOBAL
+  ADD-PARAM PNAME 'D-FUEL' PARKEY 'DF' GLOBAL
+  ADD-PARAM PNAME 'D-COOL' PARKEY 'DC' GLOBAL
+  ADD-PARAM PNAME 'B-ALPHA' PARKEY 'BAR' GLOBAL
+  ADD-PARAM PNAME 'P-HE3' PARKEY 'HE3' GLOBAL
+  BTYPE INST-BURN
+  INST-BVAL SAME 0.0
+  REACTOR-POW <<pow0>>
+  AXIAL-PFORM
+     0.03584503
+     0.04259100
+     0.05006289
+     0.06053855
+     0.07103299
+     0.07919949
+     0.08223456
+     0.08210314
+     0.08027585
+     0.07814192
+     0.07887409
+     0.07456868
+     0.06580141
+     0.04879880
+     0.02381116
+     0.01429921
+     0.01148318
+     0.00873597
+     0.00684610
+     0.00475597
+  ;
+END: ;
+QUIT .
diff --git a/Donjon/data/pulseTHM_0d_proc/PkinPULSE.c2m b/Donjon/data/pulseTHM_0d_proc/PkinPULSE.c2m
new file mode 100644
index 0000000..29074aa
--- /dev/null
+++ b/Donjon/data/pulseTHM_0d_proc/PkinPULSE.c2m
@@ -0,0 +1,172 @@
+****************************************************************
+*                                                              *
+* Procedure :  PkinPULSE.c2m                                   *
+* Purpose   :  Define point kinetic and feedback data for the  *
+*              a pulse reactor                                 *
+* Author    :  A. Hebert (2017)                                *
+*                                                              *
+* CALL      :  FMAP := PkinPULSE                               *
+*                     :: <<dt>> <<pow0>>  ;                    *
+*                                                              *
+****************************************************************
+
+PARAMETER  FMAP :: ::: LINKED_LIST FMAP ; ;
+MODULE     END: PKINI: ;
+REAL dt pow0 ;
+  :: >>dt<< >>pow0<< ;
+
+FMAP := PKINI: FMAP ::
+   EDIT 2
+   TIME -1.0 <<dt>>
+   POWER <<pow0>>
+   LAMBDA 0.0000277
+   EPSILON 1.0E-2
+   NGROUP 6
+   BETAI
+       0.00025014
+       0.00166002
+       0.00148568
+       0.0029941
+       0.0008717
+       0.00031836
+   LAMBDAI
+       0.0124
+       0.0305
+       0.111
+       0.301
+       1.14
+       3.01
+   ALPHA
+   T-FUEL SQDERIV 2
+       280.  -0.0015
+       2000. -0.0015
+   T-COOL DERIV 2
+       280. 0.
+       500. 0.
+   D-FUEL DERIV 2
+       0.   0.0
+       10.0 0.0
+   D-COOL DERIV 2
+       0.  0.
+       1.0 0.
+   B-ALPHA DIRECT 2 (*BCS bar insertion*)
+        0. 0.
+        1. -0.09
+   P-HE3 DIRECT 102
+0. 0.
+0.0508475 -0.000352113
+0.135593 -0.000985915
+0.20339 -0.00161972
+0.288136 -0.00225352
+0.355932 -0.0028169
+0.440678 -0.0034507
+0.525424 -0.00408451
+0.677966 -0.00528169
+0.762712 -0.00591549
+0.847458 -0.0065493
+0.915254 -0.00711268
+1.05085 -0.00802817
+1.13559 -0.00859155
+1.23729 -0.00922535
+1.40678 -0.0102817
+1.50847 -0.0108451
+1.61017 -0.0114085
+1.69492 -0.0120423
+1.79661 -0.0126056
+1.89831 -0.013169
+2.0339 -0.0138732
+2.13559 -0.0144366
+2.33898 -0.0154225
+2.45763 -0.0159155
+2.55932 -0.0164085
+2.67797 -0.0169014
+2.79661 -0.0173944
+2.91525 -0.0178873
+3.0339 -0.0183803
+3.15254 -0.0188028
+3.28814 -0.0192958
+3.40678 -0.0197183
+3.62712 -0.0204225
+3.76271 -0.0208451
+3.89831 -0.0212676
+4.0339 -0.0216901
+4.16949 -0.0221127
+4.32203 -0.0224648
+4.45763 -0.0228873
+4.61017 -0.0232394
+4.74576 -0.0235915
+4.89831 -0.0239437
+5.0339 -0.0242254
+5.18644 -0.0245775
+5.33898 -0.0248592
+5.62712 -0.025493
+5.77966 -0.0257746
+5.9322 -0.0260563
+6.01695 -0.0261972
+6.16949 -0.0264085
+6.33898 -0.0266901
+6.49153 -0.0269718
+6.64407 -0.0271831
+6.81356 -0.0272535
+6.9322 -0.0276056
+7.01695 -0.0277465
+7.16949 -0.0279577
+7.33898 -0.028169
+7.49153 -0.0283803
+7.66102 -0.0285915
+7.81356 -0.0288028
+8.01695 -0.0290845
+8.18644 -0.0292254
+8.33898 -0.0294366
+8.50847 -0.0296479
+8.66102 -0.0298592
+8.81356 -0.0300704
+9.23729 -0.0305634
+9.40678 -0.0307042
+9.55932 -0.0309155
+9.71186 -0.0310563
+9.88136 -0.0312676
+10. -0.0314085
+10.1525 -0.0316197
+10.322 -0.031831
+10.4746 -0.0319718
+10.6441 -0.0321831
+10.7966 -0.0323239
+10.9831 -0.0325352
+11.1525 -0.0327465
+11.3051 -0.0328873
+11.4746 -0.0330986
+11.6271 -0.0333099
+11.7797 -0.0334507
+11.9831 -0.0337324
+12.1525 -0.0338732
+12.3051 -0.0340845
+12.4746 -0.0342254
+12.6441 -0.0343662
+12.7966 -0.0345775
+12.9322 -0.0347887
+12.9831 -0.0347183
+13.1356 -0.0349296
+13.3051 -0.0350704
+13.9831 -0.0356338
+14.1525 -0.0357746
+14.3051 -0.0358451
+14.4746 -0.0359155
+14.6441 -0.0359859
+14.7966 -0.0361268
+15.0000 -0.03615
+ENDA
+
+PTIME
+P-HE3 (*He3 bar time law with B=6.9855 s^-1 *)
+   T-DELT -2.0 0.5 P-VALV 1.66 11.269 0.62984 0.0 6.9855
+B-ALPHA 7 (*BCS bar insertion*)
+   -2.0 0.
+   0. 0.
+   0.0768 0.
+   0.080 0.
+   0.22 0.
+   0.6 1.
+   30. 1.
+ENDP
+;