From 7dfcc480ba1e19bd3232349fc733caef94034292 Mon Sep 17 00:00:00 2001 From: stainer_t Date: Mon, 8 Sep 2025 13:48:49 +0200 Subject: Initial commit from Polytechnique Montreal --- Donjon/data/thm_BWR_pdrop_dfm_proc/assertS_low.c2m | 37 +++ Donjon/data/thm_BWR_pdrop_dfm_proc/dfm10c.c2m | 259 ++++++++++++++++++++ Donjon/data/thm_BWR_pdrop_dfm_proc/dfm10s.c2m | 259 ++++++++++++++++++++ Donjon/data/thm_BWR_pdrop_dfm_proc/dfm20c.c2m | 259 ++++++++++++++++++++ Donjon/data/thm_BWR_pdrop_dfm_proc/dfm20s.c2m | 259 ++++++++++++++++++++ Donjon/data/thm_BWR_pdrop_dfm_proc/dfm40c.c2m | 259 ++++++++++++++++++++ Donjon/data/thm_BWR_pdrop_dfm_proc/dfm40s.c2m | 259 ++++++++++++++++++++ Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp10c.c2m | 259 ++++++++++++++++++++ Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp10s.c2m | 259 ++++++++++++++++++++ Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp20c.c2m | 259 ++++++++++++++++++++ Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp20s.c2m | 262 +++++++++++++++++++++ Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp40c.c2m | 259 ++++++++++++++++++++ Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp40s.c2m | 259 ++++++++++++++++++++ 13 files changed, 3148 insertions(+) create mode 100755 Donjon/data/thm_BWR_pdrop_dfm_proc/assertS_low.c2m create mode 100755 Donjon/data/thm_BWR_pdrop_dfm_proc/dfm10c.c2m create mode 100755 Donjon/data/thm_BWR_pdrop_dfm_proc/dfm10s.c2m create mode 100755 Donjon/data/thm_BWR_pdrop_dfm_proc/dfm20c.c2m create mode 100755 Donjon/data/thm_BWR_pdrop_dfm_proc/dfm20s.c2m create mode 100755 Donjon/data/thm_BWR_pdrop_dfm_proc/dfm40c.c2m create mode 100755 Donjon/data/thm_BWR_pdrop_dfm_proc/dfm40s.c2m create mode 100755 Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp10c.c2m create mode 100755 Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp10s.c2m create mode 100755 Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp20c.c2m create mode 100755 Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp20s.c2m create mode 100755 Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp40c.c2m create mode 100755 Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp40s.c2m (limited to 'Donjon/data/thm_BWR_pdrop_dfm_proc') diff --git a/Donjon/data/thm_BWR_pdrop_dfm_proc/assertS_low.c2m b/Donjon/data/thm_BWR_pdrop_dfm_proc/assertS_low.c2m new file mode 100755 index 0000000..a0db294 --- /dev/null +++ b/Donjon/data/thm_BWR_pdrop_dfm_proc/assertS_low.c2m @@ -0,0 +1,37 @@ +* +* Assert procedure for non-regression testing +* Recover a value from a real array +* Special version with user-defined epsilon +* Author: A. Hebert +* +PARAMETER LCMNAM :: ::: LINKED_LIST LCMNAM ; ; +CHARACTER KEY ; +INTEGER ISET IPOS ; +REAL REFVALUE epsilon ; +:: >>KEY<< >>IPOS<< >>REFVALUE<< >>epsilon<< ; +INTEGER ITYLCM ; +REAL VALUE DELTA ; +DOUBLE PRECISION DVALUE ; +MODULE GREP: ABORT: END: ; +* +GREP: LCMNAM :: TYPE <> >>ITYLCM<< ; +IF ITYLCM 2 = THEN + GREP: LCMNAM :: GETVAL <> <> >>VALUE<< ; +ELSEIF ITYLCM 4 = THEN + GREP: LCMNAM :: GETVAL <> <> >>DVALUE<< ; + EVALUATE VALUE := DVALUE D_TO_R ; +ELSE + PRINT "assertS_low: INVALID TYPE=" ITYLCM ; + ABORT: ; +ENDIF ; +EVALUATE DELTA := VALUE REFVALUE - REFVALUE / ABS ; +IF DELTA epsilon < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFVALUE " CALCULATED=" VALUE ; + ABORT: ; +ENDIF ; +END: ; diff --git a/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm10c.c2m b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm10c.c2m new file mode 100755 index 0000000..9e3d07d --- /dev/null +++ b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm10c.c2m @@ -0,0 +1,259 @@ +************************************************************************************* +* Input file : dfm_10_cos +* Test for THM: DFM, NO PDROP +************************************************************************************* + +LINKED_LIST Geom Matex Fmap Thm ; +MODULE GEO: RESINI: USPLIT: THM: GREP: UTL: DELETE: ABORT: END: ; + +PROCEDURE assertS_low ; +REAL Mass := 6.464E-3 ; ! kg + +*---- +* Calculation options +*---- +REAL Tfuel := 900.0 ; (*fuel temperature (K)*) +REAL Tcool := 543.15 ; (*coolant temperature (K)*) +REAL Dfuel := 9.7413951 ; (*fuel density (g/cc)*) +REAL dens_mod_0 := 0.65 ; +REAL powi := 0.01 ; ! 10.0 kW + +INTEGER MaxR := 10000 ; +REAL maxh := 380.0 ; +REAL dz := maxh 40.0 / ; + +REAL z1 := dz ; +REAL z2 := z1 dz + ; +REAL z3 := z2 dz + ; +REAL z4 := z3 dz + ; +REAL z5 := z4 dz + ; +REAL z6 := z5 dz + ; +REAL z7 := z6 dz + ; +REAL z8 := z7 dz + ; +REAL z9 := z8 dz + ; +REAL z10 := z9 dz + ; +REAL z11 := z10 dz + ; +REAL z12 := z11 dz + ; +REAL z13 := z12 dz + ; +REAL z14 := z13 dz + ; +REAL z15 := z14 dz + ; +REAL z16 := z15 dz + ; +REAL z17 := z16 dz + ; +REAL z18 := z17 dz + ; +REAL z19 := z18 dz + ; +REAL z20 := z19 dz + ; +REAL z21 := z20 dz + ; +REAL z22 := z21 dz + ; +REAL z23 := z22 dz + ; +REAL z24 := z23 dz + ; +REAL z25 := z24 dz + ; +REAL z26 := z25 dz + ; +REAL z27 := z26 dz + ; +REAL z28 := z27 dz + ; +REAL z29 := z28 dz + ; +REAL z30 := z29 dz + ; +REAL z31 := z30 dz + ; +REAL z32 := z31 dz + ; +REAL z33 := z32 dz + ; +REAL z34 := z33 dz + ; +REAL z35 := z34 dz + ; +REAL z36 := z35 dz + ; +REAL z37 := z36 dz + ; +REAL z38 := z37 dz + ; +REAL z39 := z38 dz + ; +REAL z40 := z39 dz + ; + +REAL Pitch := 1.295 ; +Geom := GEO: :: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 ; + +Geom Matex := USPLIT: Geom :: NGRP 2 MAXR <> + NFUEL 1 FMIX 1 +; + +*-- +* Fuel map definition +*-- +Fmap Matex := RESINI: Matex :: + ::: GEO: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 +; +! +NXNAME '01' NYNAME 'A' +NCOMB 1 +B-ZONE 1 + +ADD-PARAM PNAME 'T-FUEL' PARKEY 'TFuel' GLOBAL +ADD-PARAM PNAME 'T-COOL' PARKEY 'TCool' GLOBAL +ADD-PARAM PNAME 'D-COOL' PARKEY 'DCool' GLOBAL +BTYPE INST-BURN +INST-BVAL CHAN 0.0 +REACTOR-POW <> AXIAL-PFORM +1.5704 1.5680 1.5631 1.5559 1.5462 1.5342 1.5198 1.5031 1.4840 1.4627 +1.4391 1.4132 1.3852 1.3551 1.3229 1.2886 1.2523 1.2142 1.1741 1.1322 +1.0886 1.0433 0.9964 0.9480 0.8981 0.8468 0.7942 0.7404 0.6855 0.6295 +0.5725 0.5146 0.4559 0.3966 0.3366 0.2761 0.2152 0.1540 0.0925 0.0308 +SET-PARAM 'T-FUEL' <> +SET-PARAM 'T-COOL' <> +SET-PARAM 'D-COOL' <> +FUEL WEIGHT <> +; + +*-- +* THM single-stage calculation +*-- +Thm Fmap := THM: Fmap :: + EDIT 1 + FLUID H2O + FPUISS 1.0 + CRITFL 5.0E7 + INLET 7.20E+06 (*Pa*) 543.15 (*K*) + INLET-Q 8.470E-05 (*m2*) 8.407E-02 (*inlet mass flow rate kg/s*) + ASSMB 1 0 + RADIUS 4.435E-03 4.520E-03 5.140E-03 0.000E+00 (*m*) + RODMESH 15 20 + HGAP 10000.0 + CONDC 0 21.5 KELVIN + CONDF 0 4.18 KELVIN + SAHA + PDROP 0 + DFM 1 +; + +UTL: Thm :: DIR DUMP ; + +assertS_low Thm :: 'AVG-T-FUEL' 1 6.055386E+02 1.0E-3 ; +assertS_low Thm :: 'AVG-T-COOL' 1 5.518742E+02 1.0E-3 ; +assertS_low Thm :: 'MIN-D-COOL' 1 6.177967E-01 1.0E-3 ; +assertS_low Thm :: 'MAX-D-COOL' 1 7.698399E-01 1.0E-3 ; + +Thm := UTL: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' STEP AT 1 DUMP ; +REAL PINLET EPSOUT ; +REAL DELTA ; +REAL REFPVAL := 7.200000E+06 ; +REAL REFEPSVAL := 1.697831E-01 ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'PINLET' 1 1 1 >>PINLET<< ; + +EVALUATE DELTA := PINLET REFPVAL - REFPVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFPVAL " CALCULATED=" PINLET ; + ABORT: ; +ENDIF ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'EPSOUT' 1 1 1 >>EPSOUT<< ; + +EVALUATE DELTA := EPSOUT REFEPSVAL - REFEPSVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFEPSVAL " CALCULATED=" EPSOUT ; + ABORT: ; +ENDIF ; + +ECHO "dfm10c completed" ; +END: ; \ No newline at end of file diff --git a/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm10s.c2m b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm10s.c2m new file mode 100755 index 0000000..643641e --- /dev/null +++ b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm10s.c2m @@ -0,0 +1,259 @@ +************************************************************************************* +* Input file : dfm_10_sine +* Test for THM: DFM, NO PDROP +************************************************************************************* + +LINKED_LIST Geom Matex Fmap Thm ; +MODULE GEO: RESINI: USPLIT: THM: GREP: UTL: DELETE: ABORT: END: ; + +PROCEDURE assertS_low ; +REAL Mass := 6.464E-3 ; ! kg + +*---- +* Calculation options +*---- +REAL Tfuel := 900.0 ; (*fuel temperature (K)*) +REAL Tcool := 543.15 ; (*coolant temperature (K)*) +REAL Dfuel := 9.7413951 ; (*fuel density (g/cc)*) +REAL dens_mod_0 := 0.65 ; +REAL powi := 0.01 ; ! 10.0 kW + +INTEGER MaxR := 10000 ; +REAL maxh := 380.0 ; +REAL dz := maxh 40.0 / ; + +REAL z1 := dz ; +REAL z2 := z1 dz + ; +REAL z3 := z2 dz + ; +REAL z4 := z3 dz + ; +REAL z5 := z4 dz + ; +REAL z6 := z5 dz + ; +REAL z7 := z6 dz + ; +REAL z8 := z7 dz + ; +REAL z9 := z8 dz + ; +REAL z10 := z9 dz + ; +REAL z11 := z10 dz + ; +REAL z12 := z11 dz + ; +REAL z13 := z12 dz + ; +REAL z14 := z13 dz + ; +REAL z15 := z14 dz + ; +REAL z16 := z15 dz + ; +REAL z17 := z16 dz + ; +REAL z18 := z17 dz + ; +REAL z19 := z18 dz + ; +REAL z20 := z19 dz + ; +REAL z21 := z20 dz + ; +REAL z22 := z21 dz + ; +REAL z23 := z22 dz + ; +REAL z24 := z23 dz + ; +REAL z25 := z24 dz + ; +REAL z26 := z25 dz + ; +REAL z27 := z26 dz + ; +REAL z28 := z27 dz + ; +REAL z29 := z28 dz + ; +REAL z30 := z29 dz + ; +REAL z31 := z30 dz + ; +REAL z32 := z31 dz + ; +REAL z33 := z32 dz + ; +REAL z34 := z33 dz + ; +REAL z35 := z34 dz + ; +REAL z36 := z35 dz + ; +REAL z37 := z36 dz + ; +REAL z38 := z37 dz + ; +REAL z39 := z38 dz + ; +REAL z40 := z39 dz + ; + +REAL Pitch := 1.295 ; +Geom := GEO: :: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 ; + +Geom Matex := USPLIT: Geom :: NGRP 2 MAXR <> + NFUEL 1 FMIX 1 +; + +*-- +* Fuel map definition +*-- +Fmap Matex := RESINI: Matex :: + ::: GEO: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 +; +! +NXNAME '01' NYNAME 'A' +NCOMB 1 +B-ZONE 1 + +ADD-PARAM PNAME 'T-FUEL' PARKEY 'TFuel' GLOBAL +ADD-PARAM PNAME 'T-COOL' PARKEY 'TCool' GLOBAL +ADD-PARAM PNAME 'D-COOL' PARKEY 'DCool' GLOBAL +BTYPE INST-BURN +INST-BVAL CHAN 0.0 +REACTOR-POW <> AXIAL-PFORM +0.0617 0.1846 0.3064 0.4263 0.5435 0.6575 0.7673 0.8725 0.9722 1.0660 +1.1532 1.2333 1.3057 1.3702 1.4261 1.4733 1.5114 1.5402 1.5595 1.5692 +1.5692 1.5595 1.5402 1.5114 1.4733 1.4261 1.3702 1.3057 1.2333 1.1532 +1.0660 0.9722 0.8725 0.7673 0.6575 0.5435 0.4263 0.3064 0.1846 0.0617 +SET-PARAM 'T-FUEL' <> +SET-PARAM 'T-COOL' <> +SET-PARAM 'D-COOL' <> +FUEL WEIGHT <> +; + +*-- +* THM single-stage calculation +*-- +Thm Fmap := THM: Fmap :: + EDIT 1 + FLUID H2O + FPUISS 1.0 + CRITFL 5.0E7 + INLET 7.20E+06 (*Pa*) 543.15 (*K*) + INLET-Q 8.470E-05 (*m2*) 8.407E-02 (*inlet mass flow rate kg/s*) + ASSMB 1 0 + RADIUS 4.435E-03 4.520E-03 5.140E-03 0.000E+00 (*m*) + RODMESH 15 20 + HGAP 10000.0 + CONDC 0 21.5 KELVIN + CONDF 0 4.18 KELVIN + SAHA + PDROP 0 + DFM 1 +; + +UTL: Thm :: DIR DUMP ; + +assertS_low Thm :: 'AVG-T-FUEL' 1 6.071170E+02 1.0E-3 ; +assertS_low Thm :: 'AVG-T-COOL' 1 5.541058E+02 1.0E-3 ; +assertS_low Thm :: 'MIN-D-COOL' 1 5.989761E-01 1.0E-3 ; +assertS_low Thm :: 'MAX-D-COOL' 1 7.699062E-01 1.0E-3 ; + +Thm := UTL: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' STEP AT 1 DUMP ; +REAL PINLET EPSOUT ; +REAL DELTA ; +REAL REFPVAL := 7.200000E+06 ; +REAL REFEPSVAL := 1.969337E-01 ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'PINLET' 1 1 1 >>PINLET<< ; + +EVALUATE DELTA := PINLET REFPVAL - REFPVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFPVAL " CALCULATED=" PINLET ; + ABORT: ; +ENDIF ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'EPSOUT' 1 1 1 >>EPSOUT<< ; + +EVALUATE DELTA := EPSOUT REFEPSVAL - REFEPSVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFEPSVAL " CALCULATED=" EPSOUT ; + ABORT: ; +ENDIF ; + +ECHO "dfm10s completed" ; +END: ; \ No newline at end of file diff --git a/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm20c.c2m b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm20c.c2m new file mode 100755 index 0000000..d8a16e4 --- /dev/null +++ b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm20c.c2m @@ -0,0 +1,259 @@ +************************************************************************************* +* Input file : dfm_20_cos +* Test for THM: DFM, NO PDROP +************************************************************************************* + +LINKED_LIST Geom Matex Fmap Thm ; +MODULE GEO: RESINI: USPLIT: THM: GREP: UTL: DELETE: ABORT: END: ; + +PROCEDURE assertS_low ; +REAL Mass := 6.464E-3 ; ! kg + +*---- +* Calculation options +*---- +REAL Tfuel := 900.0 ; (*fuel temperature (K)*) +REAL Tcool := 543.15 ; (*coolant temperature (K)*) +REAL Dfuel := 9.7413951 ; (*fuel density (g/cc)*) +REAL dens_mod_0 := 0.65 ; +REAL powi := 0.02 ; ! 20.0 kW + +INTEGER MaxR := 10000 ; +REAL maxh := 380.0 ; +REAL dz := maxh 40.0 / ; + +REAL z1 := dz ; +REAL z2 := z1 dz + ; +REAL z3 := z2 dz + ; +REAL z4 := z3 dz + ; +REAL z5 := z4 dz + ; +REAL z6 := z5 dz + ; +REAL z7 := z6 dz + ; +REAL z8 := z7 dz + ; +REAL z9 := z8 dz + ; +REAL z10 := z9 dz + ; +REAL z11 := z10 dz + ; +REAL z12 := z11 dz + ; +REAL z13 := z12 dz + ; +REAL z14 := z13 dz + ; +REAL z15 := z14 dz + ; +REAL z16 := z15 dz + ; +REAL z17 := z16 dz + ; +REAL z18 := z17 dz + ; +REAL z19 := z18 dz + ; +REAL z20 := z19 dz + ; +REAL z21 := z20 dz + ; +REAL z22 := z21 dz + ; +REAL z23 := z22 dz + ; +REAL z24 := z23 dz + ; +REAL z25 := z24 dz + ; +REAL z26 := z25 dz + ; +REAL z27 := z26 dz + ; +REAL z28 := z27 dz + ; +REAL z29 := z28 dz + ; +REAL z30 := z29 dz + ; +REAL z31 := z30 dz + ; +REAL z32 := z31 dz + ; +REAL z33 := z32 dz + ; +REAL z34 := z33 dz + ; +REAL z35 := z34 dz + ; +REAL z36 := z35 dz + ; +REAL z37 := z36 dz + ; +REAL z38 := z37 dz + ; +REAL z39 := z38 dz + ; +REAL z40 := z39 dz + ; + +REAL Pitch := 1.295 ; +Geom := GEO: :: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 ; + +Geom Matex := USPLIT: Geom :: NGRP 2 MAXR <> + NFUEL 1 FMIX 1 +; + +*-- +* Fuel map definition +*-- +Fmap Matex := RESINI: Matex :: + ::: GEO: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 +; +! +NXNAME '01' NYNAME 'A' +NCOMB 1 +B-ZONE 1 + +ADD-PARAM PNAME 'T-FUEL' PARKEY 'TFuel' GLOBAL +ADD-PARAM PNAME 'T-COOL' PARKEY 'TCool' GLOBAL +ADD-PARAM PNAME 'D-COOL' PARKEY 'DCool' GLOBAL +BTYPE INST-BURN +INST-BVAL CHAN 0.0 +REACTOR-POW <> AXIAL-PFORM +1.5704 1.5680 1.5631 1.5559 1.5462 1.5342 1.5198 1.5031 1.4840 1.4627 +1.4391 1.4132 1.3852 1.3551 1.3229 1.2886 1.2523 1.2142 1.1741 1.1322 +1.0886 1.0433 0.9964 0.9480 0.8981 0.8468 0.7942 0.7404 0.6855 0.6295 +0.5725 0.5146 0.4559 0.3966 0.3366 0.2761 0.2152 0.1540 0.0925 0.0308 +SET-PARAM 'T-FUEL' <> +SET-PARAM 'T-COOL' <> +SET-PARAM 'D-COOL' <> +FUEL WEIGHT <> +; + +*-- +* THM single-stage calculation +*-- +Thm Fmap := THM: Fmap :: + EDIT 1 + FLUID H2O + FPUISS 1.0 + CRITFL 5.0E7 + INLET 7.20E+06 (*Pa*) 543.15 (*K*) + INLET-Q 8.470E-05 (*m2*) 8.407E-02 (*inlet mass flow rate kg/s*) + ASSMB 1 0 + RADIUS 4.435E-03 4.520E-03 5.140E-03 0.000E+00 (*m*) + RODMESH 15 20 + HGAP 10000.0 + CONDC 0 21.5 KELVIN + CONDF 0 4.18 KELVIN + SAHA + PDROP 0 + DFM 1 +; + +UTL: Thm :: DIR DUMP ; + +assertS_low Thm :: 'AVG-T-FUEL' 1 6.583503E+02 1.0E-3 ; +assertS_low Thm :: 'AVG-T-COOL' 1 5.554367E+02 1.0E-3 ; +assertS_low Thm :: 'MIN-D-COOL' 1 3.812658E-01 1.0E-3 ; +assertS_low Thm :: 'MAX-D-COOL' 1 7.697743E-01 1.0E-3 ; + +Thm := UTL: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' STEP AT 1 DUMP ; +REAL PINLET EPSOUT ; +REAL DELTA ; +REAL REFPVAL := 7.200000E+06 ; +REAL REFEPSVAL := 5.081365E-01 ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'PINLET' 1 1 1 >>PINLET<< ; + +EVALUATE DELTA := PINLET REFPVAL - REFPVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFPVAL " CALCULATED=" PINLET ; + ABORT: ; +ENDIF ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'EPSOUT' 1 1 1 >>EPSOUT<< ; + +EVALUATE DELTA := EPSOUT REFEPSVAL - REFEPSVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFEPSVAL " CALCULATED=" EPSOUT ; + ABORT: ; +ENDIF ; + +ECHO "dfm20c" ; +END: ; \ No newline at end of file diff --git a/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm20s.c2m b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm20s.c2m new file mode 100755 index 0000000..16ebbf7 --- /dev/null +++ b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm20s.c2m @@ -0,0 +1,259 @@ +************************************************************************************* +* Input file : dfm_20_sine +* Test for THM: DFM, NO PDROP +************************************************************************************* + +LINKED_LIST Geom Matex Fmap Thm ; +MODULE GEO: RESINI: USPLIT: THM: GREP: UTL: DELETE: ABORT: END: ; + +PROCEDURE assertS_low ; +REAL Mass := 6.464E-3 ; ! kg + +*---- +* Calculation options +*---- +REAL Tfuel := 900.0 ; (*fuel temperature (K)*) +REAL Tcool := 543.15 ; (*coolant temperature (K)*) +REAL Dfuel := 9.7413951 ; (*fuel density (g/cc)*) +REAL dens_mod_0 := 0.65 ; +REAL powi := 0.02 ; ! 20.0 kW + +INTEGER MaxR := 10000 ; +REAL maxh := 380.0 ; +REAL dz := maxh 40.0 / ; + +REAL z1 := dz ; +REAL z2 := z1 dz + ; +REAL z3 := z2 dz + ; +REAL z4 := z3 dz + ; +REAL z5 := z4 dz + ; +REAL z6 := z5 dz + ; +REAL z7 := z6 dz + ; +REAL z8 := z7 dz + ; +REAL z9 := z8 dz + ; +REAL z10 := z9 dz + ; +REAL z11 := z10 dz + ; +REAL z12 := z11 dz + ; +REAL z13 := z12 dz + ; +REAL z14 := z13 dz + ; +REAL z15 := z14 dz + ; +REAL z16 := z15 dz + ; +REAL z17 := z16 dz + ; +REAL z18 := z17 dz + ; +REAL z19 := z18 dz + ; +REAL z20 := z19 dz + ; +REAL z21 := z20 dz + ; +REAL z22 := z21 dz + ; +REAL z23 := z22 dz + ; +REAL z24 := z23 dz + ; +REAL z25 := z24 dz + ; +REAL z26 := z25 dz + ; +REAL z27 := z26 dz + ; +REAL z28 := z27 dz + ; +REAL z29 := z28 dz + ; +REAL z30 := z29 dz + ; +REAL z31 := z30 dz + ; +REAL z32 := z31 dz + ; +REAL z33 := z32 dz + ; +REAL z34 := z33 dz + ; +REAL z35 := z34 dz + ; +REAL z36 := z35 dz + ; +REAL z37 := z36 dz + ; +REAL z38 := z37 dz + ; +REAL z39 := z38 dz + ; +REAL z40 := z39 dz + ; + +REAL Pitch := 1.295 ; +Geom := GEO: :: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 ; + +Geom Matex := USPLIT: Geom :: NGRP 2 MAXR <> + NFUEL 1 FMIX 1 +; + +*-- +* Fuel map definition +*-- +Fmap Matex := RESINI: Matex :: + ::: GEO: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 +; +! +NXNAME '01' NYNAME 'A' +NCOMB 1 +B-ZONE 1 + +ADD-PARAM PNAME 'T-FUEL' PARKEY 'TFuel' GLOBAL +ADD-PARAM PNAME 'T-COOL' PARKEY 'TCool' GLOBAL +ADD-PARAM PNAME 'D-COOL' PARKEY 'DCool' GLOBAL +BTYPE INST-BURN +INST-BVAL CHAN 0.0 +REACTOR-POW <> AXIAL-PFORM +0.0617 0.1846 0.3064 0.4263 0.5435 0.6575 0.7673 0.8725 0.9722 1.0660 +1.1532 1.2333 1.3057 1.3702 1.4261 1.4733 1.5114 1.5402 1.5595 1.5692 +1.5692 1.5595 1.5402 1.5114 1.4733 1.4261 1.3702 1.3057 1.2333 1.1532 +1.0660 0.9722 0.8725 0.7673 0.6575 0.5435 0.4263 0.3064 0.1846 0.0617 +SET-PARAM 'T-FUEL' <> +SET-PARAM 'T-COOL' <> +SET-PARAM 'D-COOL' <> +FUEL WEIGHT <> +; + +*-- +* THM single-stage calculation +*-- +Thm Fmap := THM: Fmap :: + EDIT 1 + FLUID H2O + FPUISS 1.0 + CRITFL 5.0E7 + INLET 7.20E+06 (*Pa*) 543.15 (*K*) + INLET-Q 8.470E-05 (*m2*) 8.407E-02 (*inlet mass flow rate kg/s*) + ASSMB 1 0 + RADIUS 4.435E-03 4.520E-03 5.140E-03 0.000E+00 (*m*) + RODMESH 15 20 + HGAP 10000.0 + CONDC 0 21.5 KELVIN + CONDF 0 4.18 KELVIN + SAHA + PDROP 0 + DFM 1 +; + +UTL: Thm :: DIR DUMP ; + +assertS_low Thm :: 'AVG-T-FUEL' 1 6.578530E+02 1.0E-3 ; +assertS_low Thm :: 'AVG-T-COOL' 1 5.578923E+02 1.0E-3 ; +assertS_low Thm :: 'MIN-D-COOL' 1 3.727725E-01 1.0E-3 ; +assertS_low Thm :: 'MAX-D-COOL' 1 7.699069E-01 1.0E-3 ; + +Thm := UTL: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' STEP AT 1 DUMP ; +REAL PINLET EPSOUT ; +REAL DELTA ; +REAL REFPVAL := 7.2000000E+06 ; +REAL REFEPSVAL := 5.202804E-01 ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'PINLET' 1 1 1 >>PINLET<< ; + +EVALUATE DELTA := PINLET REFPVAL - REFPVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFPVAL " CALCULATED=" PINLET ; + ABORT: ; +ENDIF ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'EPSOUT' 1 1 1 >>EPSOUT<< ; + +EVALUATE DELTA := EPSOUT REFEPSVAL - REFEPSVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFEPSVAL " CALCULATED=" EPSOUT ; + ABORT: ; +ENDIF ; + +ECHO "dfm20s completed" ; +END: ; \ No newline at end of file diff --git a/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm40c.c2m b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm40c.c2m new file mode 100755 index 0000000..4cd1a46 --- /dev/null +++ b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm40c.c2m @@ -0,0 +1,259 @@ +************************************************************************************* +* Input file : dfm_40_cos +* Test for THM: DFM, NO PDROP +************************************************************************************* + +LINKED_LIST Geom Matex Fmap Thm ; +MODULE GEO: RESINI: USPLIT: THM: GREP: UTL: DELETE: ABORT: END: ; + +PROCEDURE assertS_low ; +REAL Mass := 6.464E-3 ; ! kg + +*---- +* Calculation options +*---- +REAL Tfuel := 900.0 ; (*fuel temperature (K)*) +REAL Tcool := 543.15 ; (*coolant temperature (K)*) +REAL Dfuel := 9.7413951 ; (*fuel density (g/cc)*) +REAL dens_mod_0 := 0.65 ; +REAL powi := 0.04 ; ! 40.0 kW + +INTEGER MaxR := 10000 ; +REAL maxh := 380.0 ; +REAL dz := maxh 40.0 / ; + +REAL z1 := dz ; +REAL z2 := z1 dz + ; +REAL z3 := z2 dz + ; +REAL z4 := z3 dz + ; +REAL z5 := z4 dz + ; +REAL z6 := z5 dz + ; +REAL z7 := z6 dz + ; +REAL z8 := z7 dz + ; +REAL z9 := z8 dz + ; +REAL z10 := z9 dz + ; +REAL z11 := z10 dz + ; +REAL z12 := z11 dz + ; +REAL z13 := z12 dz + ; +REAL z14 := z13 dz + ; +REAL z15 := z14 dz + ; +REAL z16 := z15 dz + ; +REAL z17 := z16 dz + ; +REAL z18 := z17 dz + ; +REAL z19 := z18 dz + ; +REAL z20 := z19 dz + ; +REAL z21 := z20 dz + ; +REAL z22 := z21 dz + ; +REAL z23 := z22 dz + ; +REAL z24 := z23 dz + ; +REAL z25 := z24 dz + ; +REAL z26 := z25 dz + ; +REAL z27 := z26 dz + ; +REAL z28 := z27 dz + ; +REAL z29 := z28 dz + ; +REAL z30 := z29 dz + ; +REAL z31 := z30 dz + ; +REAL z32 := z31 dz + ; +REAL z33 := z32 dz + ; +REAL z34 := z33 dz + ; +REAL z35 := z34 dz + ; +REAL z36 := z35 dz + ; +REAL z37 := z36 dz + ; +REAL z38 := z37 dz + ; +REAL z39 := z38 dz + ; +REAL z40 := z39 dz + ; + +REAL Pitch := 1.295 ; +Geom := GEO: :: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 ; + +Geom Matex := USPLIT: Geom :: NGRP 2 MAXR <> + NFUEL 1 FMIX 1 +; + +*-- +* Fuel map definition +*-- +Fmap Matex := RESINI: Matex :: + ::: GEO: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 +; +! +NXNAME '01' NYNAME 'A' +NCOMB 1 +B-ZONE 1 + +ADD-PARAM PNAME 'T-FUEL' PARKEY 'TFuel' GLOBAL +ADD-PARAM PNAME 'T-COOL' PARKEY 'TCool' GLOBAL +ADD-PARAM PNAME 'D-COOL' PARKEY 'DCool' GLOBAL +BTYPE INST-BURN +INST-BVAL CHAN 0.0 +REACTOR-POW <> AXIAL-PFORM +1.5704 1.5680 1.5631 1.5559 1.5462 1.5342 1.5198 1.5031 1.4840 1.4627 +1.4391 1.4132 1.3852 1.3551 1.3229 1.2886 1.2523 1.2142 1.1741 1.1322 +1.0886 1.0433 0.9964 0.9480 0.8981 0.8468 0.7942 0.7404 0.6855 0.6295 +0.5725 0.5146 0.4559 0.3966 0.3366 0.2761 0.2152 0.1540 0.0925 0.0308 +SET-PARAM 'T-FUEL' <> +SET-PARAM 'T-COOL' <> +SET-PARAM 'D-COOL' <> +FUEL WEIGHT <> +; + +*-- +* THM single-stage calculation +*-- +Thm Fmap := THM: Fmap :: + EDIT 1 + FLUID H2O + FPUISS 1.0 + CRITFL 5.0E7 + INLET 7.20E+06 (*Pa*) 543.15 (*K*) + INLET-Q 8.470E-05 (*m2*) 8.407E-02 (*inlet mass flow rate kg/s*) + ASSMB 1 0 + RADIUS 4.435E-03 4.520E-03 5.140E-03 0.000E+00 (*m*) + RODMESH 15 20 + HGAP 10000.0 + CONDC 0 21.5 KELVIN + CONDF 0 4.18 KELVIN + SAHA + PDROP 0 + DFM 1 +; + +UTL: Thm :: DIR DUMP ; + +assertS_low Thm :: 'AVG-T-FUEL' 1 7.548829E+02 1.0E-3 ; +assertS_low Thm :: 'AVG-T-COOL' 1 5.572057E+02 1.0E-3 ; +assertS_low Thm :: 'MIN-D-COOL' 1 2.574327E-01 1.0E-3 ; +assertS_low Thm :: 'MAX-D-COOL' 1 7.696396E-01 1.0E-3 ; + +Thm := UTL: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' STEP AT 1 DUMP ; +REAL PINLET EPSOUT ; +REAL DELTA ; +REAL REFPVAL := 7.200000E+06 ; +REAL REFEPSVAL := 6.854069E-01 ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'PINLET' 1 1 1 >>PINLET<< ; + +EVALUATE DELTA := PINLET REFPVAL - REFPVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFPVAL " CALCULATED=" PINLET ; + ABORT: ; +ENDIF ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'EPSOUT' 1 1 1 >>EPSOUT<< ; + +EVALUATE DELTA := EPSOUT REFEPSVAL - REFEPSVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFEPSVAL " CALCULATED=" EPSOUT ; + ABORT: ; +ENDIF ; + +ECHO "dfm40c completed" ; +END: ; \ No newline at end of file diff --git a/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm40s.c2m b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm40s.c2m new file mode 100755 index 0000000..8564549 --- /dev/null +++ b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfm40s.c2m @@ -0,0 +1,259 @@ +************************************************************************************* +* Input file : dfm_40_sine +* Test for THM: DFM, NO PDROP +************************************************************************************* + +LINKED_LIST Geom Matex Fmap Thm ; +MODULE GEO: RESINI: USPLIT: THM: GREP: UTL: DELETE: ABORT: END: ; + +PROCEDURE assertS_low ; +REAL Mass := 6.464E-3 ; ! kg + +*---- +* Calculation options +*---- +REAL Tfuel := 900.0 ; (*fuel temperature (K)*) +REAL Tcool := 543.15 ; (*coolant temperature (K)*) +REAL Dfuel := 9.7413951 ; (*fuel density (g/cc)*) +REAL dens_mod_0 := 0.65 ; +REAL powi := 0.04 ; ! 40.0 kW + +INTEGER MaxR := 10000 ; +REAL maxh := 380.0 ; +REAL dz := maxh 40.0 / ; + +REAL z1 := dz ; +REAL z2 := z1 dz + ; +REAL z3 := z2 dz + ; +REAL z4 := z3 dz + ; +REAL z5 := z4 dz + ; +REAL z6 := z5 dz + ; +REAL z7 := z6 dz + ; +REAL z8 := z7 dz + ; +REAL z9 := z8 dz + ; +REAL z10 := z9 dz + ; +REAL z11 := z10 dz + ; +REAL z12 := z11 dz + ; +REAL z13 := z12 dz + ; +REAL z14 := z13 dz + ; +REAL z15 := z14 dz + ; +REAL z16 := z15 dz + ; +REAL z17 := z16 dz + ; +REAL z18 := z17 dz + ; +REAL z19 := z18 dz + ; +REAL z20 := z19 dz + ; +REAL z21 := z20 dz + ; +REAL z22 := z21 dz + ; +REAL z23 := z22 dz + ; +REAL z24 := z23 dz + ; +REAL z25 := z24 dz + ; +REAL z26 := z25 dz + ; +REAL z27 := z26 dz + ; +REAL z28 := z27 dz + ; +REAL z29 := z28 dz + ; +REAL z30 := z29 dz + ; +REAL z31 := z30 dz + ; +REAL z32 := z31 dz + ; +REAL z33 := z32 dz + ; +REAL z34 := z33 dz + ; +REAL z35 := z34 dz + ; +REAL z36 := z35 dz + ; +REAL z37 := z36 dz + ; +REAL z38 := z37 dz + ; +REAL z39 := z38 dz + ; +REAL z40 := z39 dz + ; + +REAL Pitch := 1.295 ; +Geom := GEO: :: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 ; + +Geom Matex := USPLIT: Geom :: NGRP 2 MAXR <> + NFUEL 1 FMIX 1 +; + +*-- +* Fuel map definition +*-- +Fmap Matex := RESINI: Matex :: + ::: GEO: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 +; +! +NXNAME '01' NYNAME 'A' +NCOMB 1 +B-ZONE 1 + +ADD-PARAM PNAME 'T-FUEL' PARKEY 'TFuel' GLOBAL +ADD-PARAM PNAME 'T-COOL' PARKEY 'TCool' GLOBAL +ADD-PARAM PNAME 'D-COOL' PARKEY 'DCool' GLOBAL +BTYPE INST-BURN +INST-BVAL CHAN 0.0 +REACTOR-POW <> AXIAL-PFORM +0.0617 0.1846 0.3064 0.4263 0.5435 0.6575 0.7673 0.8725 0.9722 1.0660 +1.1532 1.2333 1.3057 1.3702 1.4261 1.4733 1.5114 1.5402 1.5595 1.5692 +1.5692 1.5595 1.5402 1.5114 1.4733 1.4261 1.3702 1.3057 1.2333 1.1532 +1.0660 0.9722 0.8725 0.7673 0.6575 0.5435 0.4263 0.3064 0.1846 0.0617 +SET-PARAM 'T-FUEL' <> +SET-PARAM 'T-COOL' <> +SET-PARAM 'D-COOL' <> +FUEL WEIGHT <> +; + +*-- +* THM single-stage calculation +*-- +Thm Fmap := THM: Fmap :: + EDIT 1 + FLUID H2O + FPUISS 1.0 + CRITFL 5.0E7 + INLET 7.20E+06 (*Pa*) 543.15 (*K*) + INLET-Q 8.470E-05 (*m2*) 8.407E-02 (*inlet mass flow rate kg/s*) + ASSMB 1 0 + RADIUS 4.435E-03 4.520E-03 5.140E-03 0.000E+00 (*m*) + RODMESH 15 20 + HGAP 10000.0 + CONDC 0 21.5 KELVIN + CONDF 0 4.18 KELVIN + SAHA + PDROP 0 + DFM 1 +; + +UTL: Thm :: DIR DUMP ; + +assertS_low Thm :: 'AVG-T-FUEL' 1 7.534893E+02 1.0E-3 ; +assertS_low Thm :: 'AVG-T-COOL' 1 5.596287E+02 1.0E-3 ; +assertS_low Thm :: 'MIN-D-COOL' 1 2.529179E-01 1.0E-3 ; +assertS_low Thm :: 'MAX-D-COOL' 1 7.698997E-01 1.0E-3 ; + +Thm := UTL: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' STEP AT 1 DUMP ; +REAL PINLET EPSOUT ; +REAL DELTA ; +REAL REFPVAL := 7.200000E+06 ; +REAL REFEPSVAL := 6.918681E-01 ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'PINLET' 1 1 1 >>PINLET<< ; + +EVALUATE DELTA := PINLET REFPVAL - REFPVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFPVAL " CALCULATED=" PINLET ; + ABORT: ; +ENDIF ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'EPSOUT' 1 1 1 >>EPSOUT<< ; + +EVALUATE DELTA := EPSOUT REFEPSVAL - REFEPSVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFEPSVAL " CALCULATED=" EPSOUT ; + ABORT: ; +ENDIF ; + +ECHO "dfm40s completed" ; +END: ; \ No newline at end of file diff --git a/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp10c.c2m b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp10c.c2m new file mode 100755 index 0000000..f52ac06 --- /dev/null +++ b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp10c.c2m @@ -0,0 +1,259 @@ +************************************************************************************* +* Input file : dfm_pdrop_10_cos +* Test for THM: DFM+PDROP +************************************************************************************* + +LINKED_LIST Geom Matex Fmap Thm ; +MODULE GEO: RESINI: USPLIT: THM: GREP: UTL: DELETE: ABORT: END: ; + +PROCEDURE assertS_low ; +REAL Mass := 6.464E-3 ; ! kg + +*---- +* Calculation options +*---- +REAL Tfuel := 900.0 ; (*fuel temperature (K)*) +REAL Tcool := 543.15 ; (*coolant temperature (K)*) +REAL Dfuel := 9.7413951 ; (*fuel density (g/cc)*) +REAL dens_mod_0 := 0.65 ; +REAL powi := 0.01 ; ! 10.0 kW + +INTEGER MaxR := 10000 ; +REAL maxh := 380.0 ; +REAL dz := maxh 40.0 / ; + +REAL z1 := dz ; +REAL z2 := z1 dz + ; +REAL z3 := z2 dz + ; +REAL z4 := z3 dz + ; +REAL z5 := z4 dz + ; +REAL z6 := z5 dz + ; +REAL z7 := z6 dz + ; +REAL z8 := z7 dz + ; +REAL z9 := z8 dz + ; +REAL z10 := z9 dz + ; +REAL z11 := z10 dz + ; +REAL z12 := z11 dz + ; +REAL z13 := z12 dz + ; +REAL z14 := z13 dz + ; +REAL z15 := z14 dz + ; +REAL z16 := z15 dz + ; +REAL z17 := z16 dz + ; +REAL z18 := z17 dz + ; +REAL z19 := z18 dz + ; +REAL z20 := z19 dz + ; +REAL z21 := z20 dz + ; +REAL z22 := z21 dz + ; +REAL z23 := z22 dz + ; +REAL z24 := z23 dz + ; +REAL z25 := z24 dz + ; +REAL z26 := z25 dz + ; +REAL z27 := z26 dz + ; +REAL z28 := z27 dz + ; +REAL z29 := z28 dz + ; +REAL z30 := z29 dz + ; +REAL z31 := z30 dz + ; +REAL z32 := z31 dz + ; +REAL z33 := z32 dz + ; +REAL z34 := z33 dz + ; +REAL z35 := z34 dz + ; +REAL z36 := z35 dz + ; +REAL z37 := z36 dz + ; +REAL z38 := z37 dz + ; +REAL z39 := z38 dz + ; +REAL z40 := z39 dz + ; + +REAL Pitch := 1.295 ; +Geom := GEO: :: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 ; + +Geom Matex := USPLIT: Geom :: NGRP 2 MAXR <> + NFUEL 1 FMIX 1 +; + +*-- +* Fuel map definition +*-- +Fmap Matex := RESINI: Matex :: + ::: GEO: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 +; +! +NXNAME '01' NYNAME 'A' +NCOMB 1 +B-ZONE 1 + +ADD-PARAM PNAME 'T-FUEL' PARKEY 'TFuel' GLOBAL +ADD-PARAM PNAME 'T-COOL' PARKEY 'TCool' GLOBAL +ADD-PARAM PNAME 'D-COOL' PARKEY 'DCool' GLOBAL +BTYPE INST-BURN +INST-BVAL CHAN 0.0 +REACTOR-POW <> AXIAL-PFORM +1.5704 1.5680 1.5631 1.5559 1.5462 1.5342 1.5198 1.5031 1.4840 1.4627 +1.4391 1.4132 1.3852 1.3551 1.3229 1.2886 1.2523 1.2142 1.1741 1.1322 +1.0886 1.0433 0.9964 0.9480 0.8981 0.8468 0.7942 0.7404 0.6855 0.6295 +0.5725 0.5146 0.4559 0.3966 0.3366 0.2761 0.2152 0.1540 0.0925 0.0308 +SET-PARAM 'T-FUEL' <> +SET-PARAM 'T-COOL' <> +SET-PARAM 'D-COOL' <> +FUEL WEIGHT <> +; + +*-- +* THM single-stage calculation +*-- +Thm Fmap := THM: Fmap :: + EDIT 1 + FLUID H2O + FPUISS 1.0 + CRITFL 5.0E7 + INLET 7.20E+06 (*Pa*) 543.15 (*K*) + INLET-Q 8.470E-05 (*m2*) 8.407E-02 (*inlet mass flow rate kg/s*) + ASSMB 1 0 + RADIUS 4.435E-03 4.520E-03 5.140E-03 0.000E+00 (*m*) + RODMESH 15 20 + HGAP 10000.0 + CONDC 0 21.5 KELVIN + CONDF 0 4.18 KELVIN + SAHA + PDROP 1 + DFM 1 +; + +UTL: Thm :: DIR DUMP ; + +assertS_low Thm :: 'AVG-T-FUEL' 1 6.055386E+02 1.0E-3 ; +assertS_low Thm :: 'AVG-T-COOL' 1 5.518742E+02 1.0E-3 ; +assertS_low Thm :: 'MIN-D-COOL' 1 6.159302E-01 1.0E-3 ; +assertS_low Thm :: 'MAX-D-COOL' 1 7.698969E-01 1.0E-3 ; + +Thm := UTL: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' STEP AT 1 DUMP ; +REAL PINLET EPSOUT ; +REAL DELTA ; +REAL REFPVAL := 7.237462E+06 ; +REAL REFEPSVAL := 1.725661E-01 ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'PINLET' 1 1 1 >>PINLET<< ; + +EVALUATE DELTA := PINLET REFPVAL - REFPVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFPVAL " CALCULATED=" PINLET ; + ABORT: ; +ENDIF ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'EPSOUT' 1 1 1 >>EPSOUT<< ; + +EVALUATE DELTA := EPSOUT REFEPSVAL - REFEPSVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFEPSVAL " CALCULATED=" EPSOUT ; + ABORT: ; +ENDIF ; + +ECHO "dfmp10c completed" ; +END: ; \ No newline at end of file diff --git a/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp10s.c2m b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp10s.c2m new file mode 100755 index 0000000..92d732d --- /dev/null +++ b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp10s.c2m @@ -0,0 +1,259 @@ +************************************************************************************* +* Input file : dfm_pdrop_10_sine +* Test for THM: DFM+PDROP +************************************************************************************* + +LINKED_LIST Geom Matex Fmap Thm ; +MODULE GEO: RESINI: USPLIT: THM: GREP: UTL: DELETE: ABORT: END: ; + +PROCEDURE assertS_low ; +REAL Mass := 6.464E-3 ; ! kg + +*---- +* Calculation options +*---- +REAL Tfuel := 900.0 ; (*fuel temperature (K)*) +REAL Tcool := 543.15 ; (*coolant temperature (K)*) +REAL Dfuel := 9.7413951 ; (*fuel density (g/cc)*) +REAL dens_mod_0 := 0.65 ; +REAL powi := 0.01 ; ! 10.0 kW + +INTEGER MaxR := 10000 ; +REAL maxh := 380.0 ; +REAL dz := maxh 40.0 / ; + +REAL z1 := dz ; +REAL z2 := z1 dz + ; +REAL z3 := z2 dz + ; +REAL z4 := z3 dz + ; +REAL z5 := z4 dz + ; +REAL z6 := z5 dz + ; +REAL z7 := z6 dz + ; +REAL z8 := z7 dz + ; +REAL z9 := z8 dz + ; +REAL z10 := z9 dz + ; +REAL z11 := z10 dz + ; +REAL z12 := z11 dz + ; +REAL z13 := z12 dz + ; +REAL z14 := z13 dz + ; +REAL z15 := z14 dz + ; +REAL z16 := z15 dz + ; +REAL z17 := z16 dz + ; +REAL z18 := z17 dz + ; +REAL z19 := z18 dz + ; +REAL z20 := z19 dz + ; +REAL z21 := z20 dz + ; +REAL z22 := z21 dz + ; +REAL z23 := z22 dz + ; +REAL z24 := z23 dz + ; +REAL z25 := z24 dz + ; +REAL z26 := z25 dz + ; +REAL z27 := z26 dz + ; +REAL z28 := z27 dz + ; +REAL z29 := z28 dz + ; +REAL z30 := z29 dz + ; +REAL z31 := z30 dz + ; +REAL z32 := z31 dz + ; +REAL z33 := z32 dz + ; +REAL z34 := z33 dz + ; +REAL z35 := z34 dz + ; +REAL z36 := z35 dz + ; +REAL z37 := z36 dz + ; +REAL z38 := z37 dz + ; +REAL z39 := z38 dz + ; +REAL z40 := z39 dz + ; + +REAL Pitch := 1.295 ; +Geom := GEO: :: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 ; + +Geom Matex := USPLIT: Geom :: NGRP 2 MAXR <> + NFUEL 1 FMIX 1 +; + +*-- +* Fuel map definition +*-- +Fmap Matex := RESINI: Matex :: + ::: GEO: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 +; +! +NXNAME '01' NYNAME 'A' +NCOMB 1 +B-ZONE 1 + +ADD-PARAM PNAME 'T-FUEL' PARKEY 'TFuel' GLOBAL +ADD-PARAM PNAME 'T-COOL' PARKEY 'TCool' GLOBAL +ADD-PARAM PNAME 'D-COOL' PARKEY 'DCool' GLOBAL +BTYPE INST-BURN +INST-BVAL CHAN 0.0 +REACTOR-POW <> AXIAL-PFORM +0.0617 0.1846 0.3064 0.4263 0.5435 0.6575 0.7673 0.8725 0.9722 1.0660 +1.1532 1.2333 1.3057 1.3702 1.4261 1.4733 1.5114 1.5402 1.5595 1.5692 +1.5692 1.5595 1.5402 1.5114 1.4733 1.4261 1.3702 1.3057 1.2333 1.1532 +1.0660 0.9722 0.8725 0.7673 0.6575 0.5435 0.4263 0.3064 0.1846 0.0617 +SET-PARAM 'T-FUEL' <> +SET-PARAM 'T-COOL' <> +SET-PARAM 'D-COOL' <> +FUEL WEIGHT <> +; + +*-- +* THM single-stage calculation +*-- +Thm Fmap := THM: Fmap :: + EDIT 1 + FLUID H2O + FPUISS 1.0 + CRITFL 5.0E7 + INLET 7.20E+06 (*Pa*) 543.15 (*K*) + INLET-Q 8.470E-05 (*m2*) 8.407E-02 (*inlet mass flow rate kg/s*) + ASSMB 1 0 + RADIUS 4.435E-03 4.520E-03 5.140E-03 0.000E+00 (*m*) + RODMESH 15 20 + HGAP 10000.0 + CONDC 0 21.5 KELVIN + CONDF 0 4.18 KELVIN + SAHA + PDROP 1 + DFM 1 +; + +UTL: Thm :: DIR DUMP ; + +assertS_low Thm :: 'AVG-T-FUEL' 1 6.071170E+02 1.0E-3 ; +assertS_low Thm :: 'AVG-T-COOL' 1 5.541058E+02 1.0E-3 ; +assertS_low Thm :: 'MIN-D-COOL' 1 5.966221E-01 1.0E-3 ; +assertS_low Thm :: 'MAX-D-COOL' 1 7.699603E-01 1.0E-3 ; + +Thm := UTL: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' STEP AT 1 DUMP ; +REAL PINLET EPSOUT ; +REAL DELTA ; +REAL REFPVAL := 7.237578E+06 ; +REAL REFEPSVAL := 2.004285E-01 ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'PINLET' 1 1 1 >>PINLET<< ; + +EVALUATE DELTA := PINLET REFPVAL - REFPVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFPVAL " CALCULATED=" PINLET ; + ABORT: ; +ENDIF ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'EPSOUT' 1 1 1 >>EPSOUT<< ; + +EVALUATE DELTA := EPSOUT REFEPSVAL - REFEPSVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFEPSVAL " CALCULATED=" EPSOUT ; + ABORT: ; +ENDIF ; + +ECHO "dfmp10s completed" ; +END: ; \ No newline at end of file diff --git a/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp20c.c2m b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp20c.c2m new file mode 100755 index 0000000..8d9b26b --- /dev/null +++ b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp20c.c2m @@ -0,0 +1,259 @@ +************************************************************************************* +* Input file : dfm_pdrop_20_cos +* Test for THM: DFM+PDROP +************************************************************************************* + +LINKED_LIST Geom Matex Fmap Thm ; +MODULE GEO: RESINI: USPLIT: THM: GREP: UTL: DELETE: ABORT: END: ; + +PROCEDURE assertS_low ; +REAL Mass := 6.464E-3 ; ! kg + +*---- +* Calculation options +*---- +REAL Tfuel := 900.0 ; (*fuel temperature (K)*) +REAL Tcool := 543.15 ; (*coolant temperature (K)*) +REAL Dfuel := 9.7413951 ; (*fuel density (g/cc)*) +REAL dens_mod_0 := 0.65 ; +REAL powi := 0.02 ; ! 20.0 kW + +INTEGER MaxR := 10000 ; +REAL maxh := 380.0 ; +REAL dz := maxh 40.0 / ; + +REAL z1 := dz ; +REAL z2 := z1 dz + ; +REAL z3 := z2 dz + ; +REAL z4 := z3 dz + ; +REAL z5 := z4 dz + ; +REAL z6 := z5 dz + ; +REAL z7 := z6 dz + ; +REAL z8 := z7 dz + ; +REAL z9 := z8 dz + ; +REAL z10 := z9 dz + ; +REAL z11 := z10 dz + ; +REAL z12 := z11 dz + ; +REAL z13 := z12 dz + ; +REAL z14 := z13 dz + ; +REAL z15 := z14 dz + ; +REAL z16 := z15 dz + ; +REAL z17 := z16 dz + ; +REAL z18 := z17 dz + ; +REAL z19 := z18 dz + ; +REAL z20 := z19 dz + ; +REAL z21 := z20 dz + ; +REAL z22 := z21 dz + ; +REAL z23 := z22 dz + ; +REAL z24 := z23 dz + ; +REAL z25 := z24 dz + ; +REAL z26 := z25 dz + ; +REAL z27 := z26 dz + ; +REAL z28 := z27 dz + ; +REAL z29 := z28 dz + ; +REAL z30 := z29 dz + ; +REAL z31 := z30 dz + ; +REAL z32 := z31 dz + ; +REAL z33 := z32 dz + ; +REAL z34 := z33 dz + ; +REAL z35 := z34 dz + ; +REAL z36 := z35 dz + ; +REAL z37 := z36 dz + ; +REAL z38 := z37 dz + ; +REAL z39 := z38 dz + ; +REAL z40 := z39 dz + ; + +REAL Pitch := 1.295 ; +Geom := GEO: :: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 ; + +Geom Matex := USPLIT: Geom :: NGRP 2 MAXR <> + NFUEL 1 FMIX 1 +; + +*-- +* Fuel map definition +*-- +Fmap Matex := RESINI: Matex :: + ::: GEO: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 +; +! +NXNAME '01' NYNAME 'A' +NCOMB 1 +B-ZONE 1 + +ADD-PARAM PNAME 'T-FUEL' PARKEY 'TFuel' GLOBAL +ADD-PARAM PNAME 'T-COOL' PARKEY 'TCool' GLOBAL +ADD-PARAM PNAME 'D-COOL' PARKEY 'DCool' GLOBAL +BTYPE INST-BURN +INST-BVAL CHAN 0.0 +REACTOR-POW <> AXIAL-PFORM +1.5704 1.5680 1.5631 1.5559 1.5462 1.5342 1.5198 1.5031 1.4840 1.4627 +1.4391 1.4132 1.3852 1.3551 1.3229 1.2886 1.2523 1.2142 1.1741 1.1322 +1.0886 1.0433 0.9964 0.9480 0.8981 0.8468 0.7942 0.7404 0.6855 0.6295 +0.5725 0.5146 0.4559 0.3966 0.3366 0.2761 0.2152 0.1540 0.0925 0.0308 +SET-PARAM 'T-FUEL' <> +SET-PARAM 'T-COOL' <> +SET-PARAM 'D-COOL' <> +FUEL WEIGHT <> +; + +*-- +* THM single-stage calculation +*-- +Thm Fmap := THM: Fmap :: + EDIT 1 + FLUID H2O + FPUISS 1.0 + CRITFL 5.0E7 + INLET 7.20E+06 (*Pa*) 543.15 (*K*) + INLET-Q 8.470E-05 (*m2*) 8.407E-02 (*inlet mass flow rate kg/s*) + ASSMB 1 0 + RADIUS 4.435E-03 4.520E-03 5.140E-03 0.000E+00 (*m*) + RODMESH 15 20 + HGAP 10000.0 + CONDC 0 21.5 KELVIN + CONDF 0 4.18 KELVIN + SAHA + PDROP 1 + DFM 1 +; + +UTL: Thm :: DIR DUMP ; + +assertS_low Thm :: 'AVG-T-FUEL' 1 6.583503E+02 1.0E-3 ; +assertS_low Thm :: 'AVG-T-COOL' 1 5.555917E+02 1.0E-3 ; +assertS_low Thm :: 'MIN-D-COOL' 1 3.704429E-01 1.0E-3 ; +assertS_low Thm :: 'MAX-D-COOL' 1 7.698497E-01 1.0E-3 ; + +Thm := UTL: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' STEP AT 1 DUMP ; +REAL PINLET EPSOUT ; +REAL DELTA ; +REAL REFPVAL := 7.252412E+06 ; +REAL REFEPSVAL := 5.236769E-01 ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'PINLET' 1 1 1 >>PINLET<< ; + +EVALUATE DELTA := PINLET REFPVAL - REFPVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFPVAL " CALCULATED=" PINLET ; + ABORT: ; +ENDIF ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'EPSOUT' 1 1 1 >>EPSOUT<< ; + +EVALUATE DELTA := EPSOUT REFEPSVAL - REFEPSVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFEPSVAL " CALCULATED=" EPSOUT ; + ABORT: ; +ENDIF ; + +ECHO "dfmp20c" ; +END: ; \ No newline at end of file diff --git a/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp20s.c2m b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp20s.c2m new file mode 100755 index 0000000..64ebba9 --- /dev/null +++ b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp20s.c2m @@ -0,0 +1,262 @@ +************************************************************************************* +* Input file : thm_40_h380_sine_pow_20000.0W_SAHA_PDROP1_DFM1.x2m +* +* Test of a BWR single pincell +* +* The aim is to test the THM module of DONJON +**************************************************************************************************** + +LINKED_LIST Geom Matex Fmap Thm ; +MODULE GEO: RESINI: USPLIT: THM: GREP: UTL: DELETE: ABORT: END: ; + +PROCEDURE assertS_low ; +REAL Mass := 6.464E-3 ; ! kg + +*---- +* Calculation options +*---- +REAL Tfuel := 900.0 ; (*fuel temperature (K)*) +REAL Tcool := 543.15 ; (*coolant temperature (K)*) +REAL Dfuel := 9.7413951 ; (*fuel density (g/cc)*) +REAL dens_mod_0 := 0.65 ; +REAL powi := 0.02 ; ! 20.0 kW + +INTEGER MaxR := 10000 ; +REAL maxh := 380.0 ; +REAL dz := maxh 40.0 / ; + +REAL z1 := dz ; +REAL z2 := z1 dz + ; +REAL z3 := z2 dz + ; +REAL z4 := z3 dz + ; +REAL z5 := z4 dz + ; +REAL z6 := z5 dz + ; +REAL z7 := z6 dz + ; +REAL z8 := z7 dz + ; +REAL z9 := z8 dz + ; +REAL z10 := z9 dz + ; +REAL z11 := z10 dz + ; +REAL z12 := z11 dz + ; +REAL z13 := z12 dz + ; +REAL z14 := z13 dz + ; +REAL z15 := z14 dz + ; +REAL z16 := z15 dz + ; +REAL z17 := z16 dz + ; +REAL z18 := z17 dz + ; +REAL z19 := z18 dz + ; +REAL z20 := z19 dz + ; +REAL z21 := z20 dz + ; +REAL z22 := z21 dz + ; +REAL z23 := z22 dz + ; +REAL z24 := z23 dz + ; +REAL z25 := z24 dz + ; +REAL z26 := z25 dz + ; +REAL z27 := z26 dz + ; +REAL z28 := z27 dz + ; +REAL z29 := z28 dz + ; +REAL z30 := z29 dz + ; +REAL z31 := z30 dz + ; +REAL z32 := z31 dz + ; +REAL z33 := z32 dz + ; +REAL z34 := z33 dz + ; +REAL z35 := z34 dz + ; +REAL z36 := z35 dz + ; +REAL z37 := z36 dz + ; +REAL z38 := z37 dz + ; +REAL z39 := z38 dz + ; +REAL z40 := z39 dz + ; + +REAL Pitch := 1.295 ; +Geom := GEO: :: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 ; + +Geom Matex := USPLIT: Geom :: NGRP 2 MAXR <> + NFUEL 1 FMIX 1 +; + +*-- +* Fuel map definition +*-- +Fmap Matex := RESINI: Matex :: + ::: GEO: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 +; +! +NXNAME '01' NYNAME 'A' +NCOMB 1 +B-ZONE 1 + +ADD-PARAM PNAME 'T-FUEL' PARKEY 'TFuel' GLOBAL +ADD-PARAM PNAME 'T-COOL' PARKEY 'TCool' GLOBAL +ADD-PARAM PNAME 'D-COOL' PARKEY 'DCool' GLOBAL +BTYPE INST-BURN +INST-BVAL CHAN 0.0 +REACTOR-POW <> AXIAL-PFORM +0.0617 0.1846 0.3064 0.4263 0.5435 0.6575 0.7673 0.8725 0.9722 1.0660 +1.1532 1.2333 1.3057 1.3702 1.4261 1.4733 1.5114 1.5402 1.5595 1.5692 +1.5692 1.5595 1.5402 1.5114 1.4733 1.4261 1.3702 1.3057 1.2333 1.1532 +1.0660 0.9722 0.8725 0.7673 0.6575 0.5435 0.4263 0.3064 0.1846 0.0617 +SET-PARAM 'T-FUEL' <> +SET-PARAM 'T-COOL' <> +SET-PARAM 'D-COOL' <> +FUEL WEIGHT <> +; + +*-- +* THM single-stage calculation +*-- +Thm Fmap := THM: Fmap :: + EDIT 1 + FLUID H2O + FPUISS 1.0 + CRITFL 5.0E7 + INLET 7.20E+06 (*Pa*) 543.15 (*K*) + INLET-Q 8.470E-05 (*m2*) 8.407E-02 (*inlet mass flow rate kg/s*) + ASSMB 1 0 + RADIUS 4.435E-03 4.520E-03 5.140E-03 0.000E+00 (*m*) + RODMESH 15 20 + HGAP 10000.0 + CONDC 0 21.5 KELVIN + CONDF 0 4.18 KELVIN + SAHA + PDROP 1 + DFM 1 +; + +UTL: Thm :: DIR DUMP ; + +assertS_low Thm :: 'AVG-T-FUEL' 1 6.580398E+02 1.0E-3 ; +assertS_low Thm :: 'AVG-T-COOL' 1 5.580008E+02 1.0E-3 ; +assertS_low Thm :: 'MIN-D-COOL' 1 3.616079E-01 1.0E-3 ; +assertS_low Thm :: 'MAX-D-COOL' 1 7.699757E-01 1.0E-3 ; + +Thm := UTL: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' STEP AT 1 DUMP ; +REAL PINLET EPSOUT ; +REAL DELTA ; +REAL REFPVAL := 7.249683E+06 ; +REAL REFEPSVAL := 5.363159E-01 ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'PINLET' 1 1 1 >>PINLET<< ; + +EVALUATE DELTA := PINLET REFPVAL - REFPVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFPVAL " CALCULATED=" PINLET ; + ABORT: ; +ENDIF ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'EPSOUT' 1 1 1 >>EPSOUT<< ; + +EVALUATE DELTA := EPSOUT REFEPSVAL - REFEPSVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFEPSVAL " CALCULATED=" EPSOUT ; + ABORT: ; +ENDIF ; + +ECHO "dfmp20s" ; +END: ; \ No newline at end of file diff --git a/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp40c.c2m b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp40c.c2m new file mode 100755 index 0000000..0f0ca1f --- /dev/null +++ b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp40c.c2m @@ -0,0 +1,259 @@ +************************************************************************************* +* Input file : dfm_pdrop_40_cos +* Test for THM: DFM+PDROP +************************************************************************************* + +LINKED_LIST Geom Matex Fmap Thm ; +MODULE GEO: RESINI: USPLIT: THM: GREP: UTL: DELETE: ABORT: END: ; + +PROCEDURE assertS_low ; +REAL Mass := 6.464E-3 ; ! kg + +*---- +* Calculation options +*---- +REAL Tfuel := 900.0 ; (*fuel temperature (K)*) +REAL Tcool := 543.15 ; (*coolant temperature (K)*) +REAL Dfuel := 9.7413951 ; (*fuel density (g/cc)*) +REAL dens_mod_0 := 0.65 ; +REAL powi := 0.04 ; ! 40.0 kW + +INTEGER MaxR := 10000 ; +REAL maxh := 380.0 ; +REAL dz := maxh 40.0 / ; + +REAL z1 := dz ; +REAL z2 := z1 dz + ; +REAL z3 := z2 dz + ; +REAL z4 := z3 dz + ; +REAL z5 := z4 dz + ; +REAL z6 := z5 dz + ; +REAL z7 := z6 dz + ; +REAL z8 := z7 dz + ; +REAL z9 := z8 dz + ; +REAL z10 := z9 dz + ; +REAL z11 := z10 dz + ; +REAL z12 := z11 dz + ; +REAL z13 := z12 dz + ; +REAL z14 := z13 dz + ; +REAL z15 := z14 dz + ; +REAL z16 := z15 dz + ; +REAL z17 := z16 dz + ; +REAL z18 := z17 dz + ; +REAL z19 := z18 dz + ; +REAL z20 := z19 dz + ; +REAL z21 := z20 dz + ; +REAL z22 := z21 dz + ; +REAL z23 := z22 dz + ; +REAL z24 := z23 dz + ; +REAL z25 := z24 dz + ; +REAL z26 := z25 dz + ; +REAL z27 := z26 dz + ; +REAL z28 := z27 dz + ; +REAL z29 := z28 dz + ; +REAL z30 := z29 dz + ; +REAL z31 := z30 dz + ; +REAL z32 := z31 dz + ; +REAL z33 := z32 dz + ; +REAL z34 := z33 dz + ; +REAL z35 := z34 dz + ; +REAL z36 := z35 dz + ; +REAL z37 := z36 dz + ; +REAL z38 := z37 dz + ; +REAL z39 := z38 dz + ; +REAL z40 := z39 dz + ; + +REAL Pitch := 1.295 ; +Geom := GEO: :: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 ; + +Geom Matex := USPLIT: Geom :: NGRP 2 MAXR <> + NFUEL 1 FMIX 1 +; + +*-- +* Fuel map definition +*-- +Fmap Matex := RESINI: Matex :: + ::: GEO: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 +; +! +NXNAME '01' NYNAME 'A' +NCOMB 1 +B-ZONE 1 + +ADD-PARAM PNAME 'T-FUEL' PARKEY 'TFuel' GLOBAL +ADD-PARAM PNAME 'T-COOL' PARKEY 'TCool' GLOBAL +ADD-PARAM PNAME 'D-COOL' PARKEY 'DCool' GLOBAL +BTYPE INST-BURN +INST-BVAL CHAN 0.0 +REACTOR-POW <> AXIAL-PFORM +1.5704 1.5680 1.5631 1.5559 1.5462 1.5342 1.5198 1.5031 1.4840 1.4627 +1.4391 1.4132 1.3852 1.3551 1.3229 1.2886 1.2523 1.2142 1.1741 1.1322 +1.0886 1.0433 0.9964 0.9480 0.8981 0.8468 0.7942 0.7404 0.6855 0.6295 +0.5725 0.5146 0.4559 0.3966 0.3366 0.2761 0.2152 0.1540 0.0925 0.0308 +SET-PARAM 'T-FUEL' <> +SET-PARAM 'T-COOL' <> +SET-PARAM 'D-COOL' <> +FUEL WEIGHT <> +; + +*-- +* THM single-stage calculation +*-- +Thm Fmap := THM: Fmap :: + EDIT 1 + FLUID H2O + FPUISS 1.0 + CRITFL 5.0E7 + INLET 7.20E+06 (*Pa*) 543.15 (*K*) + INLET-Q 8.470E-05 (*m2*) 8.407E-02 (*inlet mass flow rate kg/s*) + ASSMB 1 0 + RADIUS 4.435E-03 4.520E-03 5.140E-03 0.000E+00 (*m*) + RODMESH 15 20 + HGAP 10000.0 + CONDC 0 21.5 KELVIN + CONDF 0 4.18 KELVIN + SAHA + PDROP 1 + DFM 1 +; + +UTL: Thm :: DIR DUMP ; + +assertS_low Thm :: 'AVG-T-FUEL' 1 7.541727E+02 1.0E-3 ; +assertS_low Thm :: 'AVG-T-COOL' 1 5.577180E+02 1.0E-3 ; +assertS_low Thm :: 'MIN-D-COOL' 1 2.291111E-01 1.0E-3 ; +assertS_low Thm :: 'MAX-D-COOL' 1 7.697761E-01 1.0E-3 ; + +Thm := UTL: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' STEP AT 1 DUMP ; +REAL PINLET EPSOUT ; +REAL DELTA ; +REAL REFPVAL := 7.297116E+06 ; +REAL REFEPSVAL := 7.259789E-01 ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'PINLET' 1 1 1 >>PINLET<< ; + +EVALUATE DELTA := PINLET REFPVAL - REFPVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFPVAL " CALCULATED=" PINLET ; + ABORT: ; +ENDIF ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'EPSOUT' 1 1 1 >>EPSOUT<< ; + +EVALUATE DELTA := EPSOUT REFEPSVAL - REFEPSVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFEPSVAL " CALCULATED=" EPSOUT ; + ABORT: ; +ENDIF ; + +ECHO "dfmp40c completed" ; +END: ; \ No newline at end of file diff --git a/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp40s.c2m b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp40s.c2m new file mode 100755 index 0000000..58f3e60 --- /dev/null +++ b/Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp40s.c2m @@ -0,0 +1,259 @@ +************************************************************************************* +* Input file : dfm_pdrop_40_sine +* Test for THM: DFM+PDROP +************************************************************************************* + +LINKED_LIST Geom Matex Fmap Thm ; +MODULE GEO: RESINI: USPLIT: THM: GREP: UTL: DELETE: ABORT: END: ; + +PROCEDURE assertS_low ; +REAL Mass := 6.464E-3 ; ! kg + +*---- +* Calculation options +*---- +REAL Tfuel := 900.0 ; (*fuel temperature (K)*) +REAL Tcool := 543.15 ; (*coolant temperature (K)*) +REAL Dfuel := 9.7413951 ; (*fuel density (g/cc)*) +REAL dens_mod_0 := 0.65 ; +REAL powi := 0.04 ; ! 40.0 kW + +INTEGER MaxR := 10000 ; +REAL maxh := 380.0 ; +REAL dz := maxh 40.0 / ; + +REAL z1 := dz ; +REAL z2 := z1 dz + ; +REAL z3 := z2 dz + ; +REAL z4 := z3 dz + ; +REAL z5 := z4 dz + ; +REAL z6 := z5 dz + ; +REAL z7 := z6 dz + ; +REAL z8 := z7 dz + ; +REAL z9 := z8 dz + ; +REAL z10 := z9 dz + ; +REAL z11 := z10 dz + ; +REAL z12 := z11 dz + ; +REAL z13 := z12 dz + ; +REAL z14 := z13 dz + ; +REAL z15 := z14 dz + ; +REAL z16 := z15 dz + ; +REAL z17 := z16 dz + ; +REAL z18 := z17 dz + ; +REAL z19 := z18 dz + ; +REAL z20 := z19 dz + ; +REAL z21 := z20 dz + ; +REAL z22 := z21 dz + ; +REAL z23 := z22 dz + ; +REAL z24 := z23 dz + ; +REAL z25 := z24 dz + ; +REAL z26 := z25 dz + ; +REAL z27 := z26 dz + ; +REAL z28 := z27 dz + ; +REAL z29 := z28 dz + ; +REAL z30 := z29 dz + ; +REAL z31 := z30 dz + ; +REAL z32 := z31 dz + ; +REAL z33 := z32 dz + ; +REAL z34 := z33 dz + ; +REAL z35 := z34 dz + ; +REAL z36 := z35 dz + ; +REAL z37 := z36 dz + ; +REAL z38 := z37 dz + ; +REAL z39 := z38 dz + ; +REAL z40 := z39 dz + ; + +REAL Pitch := 1.295 ; +Geom := GEO: :: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 ; + +Geom Matex := USPLIT: Geom :: NGRP 2 MAXR <> + NFUEL 1 FMIX 1 +; + +*-- +* Fuel map definition +*-- +Fmap Matex := RESINI: Matex :: + ::: GEO: CAR3D 1 1 40 + X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL + MESHX 0.0 <> + MESHY 0.0 <> + MESHZ 0.0 <> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> +<> <> <> <> <> <> <> <> <> <> + 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 +PLANE 21 SAME 1 +PLANE 22 SAME 1 +PLANE 23 SAME 1 +PLANE 24 SAME 1 +PLANE 25 SAME 1 +PLANE 26 SAME 1 +PLANE 27 SAME 1 +PLANE 28 SAME 1 +PLANE 29 SAME 1 +PLANE 30 SAME 1 +PLANE 31 SAME 1 +PLANE 32 SAME 1 +PLANE 33 SAME 1 +PLANE 34 SAME 1 +PLANE 35 SAME 1 +PLANE 36 SAME 1 +PLANE 37 SAME 1 +PLANE 38 SAME 1 +PLANE 39 SAME 1 +PLANE 40 SAME 1 +; +! +NXNAME '01' NYNAME 'A' +NCOMB 1 +B-ZONE 1 + +ADD-PARAM PNAME 'T-FUEL' PARKEY 'TFuel' GLOBAL +ADD-PARAM PNAME 'T-COOL' PARKEY 'TCool' GLOBAL +ADD-PARAM PNAME 'D-COOL' PARKEY 'DCool' GLOBAL +BTYPE INST-BURN +INST-BVAL CHAN 0.0 +REACTOR-POW <> AXIAL-PFORM +0.0617 0.1846 0.3064 0.4263 0.5435 0.6575 0.7673 0.8725 0.9722 1.0660 +1.1532 1.2333 1.3057 1.3702 1.4261 1.4733 1.5114 1.5402 1.5595 1.5692 +1.5692 1.5595 1.5402 1.5114 1.4733 1.4261 1.3702 1.3057 1.2333 1.1532 +1.0660 0.9722 0.8725 0.7673 0.6575 0.5435 0.4263 0.3064 0.1846 0.0617 +SET-PARAM 'T-FUEL' <> +SET-PARAM 'T-COOL' <> +SET-PARAM 'D-COOL' <> +FUEL WEIGHT <> +; + +*-- +* THM single-stage calculation +*-- +Thm Fmap := THM: Fmap :: + EDIT 1 + FLUID H2O + FPUISS 1.0 + CRITFL 5.0E7 + INLET 7.20E+06 (*Pa*) 543.15 (*K*) + INLET-Q 8.470E-05 (*m2*) 8.407E-02 (*inlet mass flow rate kg/s*) + ASSMB 1 0 + RADIUS 4.435E-03 4.520E-03 5.140E-03 0.000E+00 (*m*) + RODMESH 15 20 + HGAP 10000.0 + CONDC 0 21.5 KELVIN + CONDF 0 4.18 KELVIN + SAHA + PDROP 1 + DFM 1 +; + +UTL: Thm :: DIR DUMP ; + +assertS_low Thm :: 'AVG-T-FUEL' 1 7.518594E+02 1.0E-3 ; +assertS_low Thm :: 'AVG-T-COOL' 1 5.600919E+02 1.0E-3 ; +assertS_low Thm :: 'MIN-D-COOL' 1 2.234586E-01 1.0E-3 ; +assertS_low Thm :: 'MAX-D-COOL' 1 7.700225E-01 1.0E-3 ; + +Thm := UTL: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' STEP AT 1 DUMP ; +REAL PINLET EPSOUT ; +REAL DELTA ; +REAL REFPVAL := 7.286176E+06 ; +REAL REFEPSVAL := 7.340708E-01 ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'PINLET' 1 1 1 >>PINLET<< ; + +EVALUATE DELTA := PINLET REFPVAL - REFPVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFPVAL " CALCULATED=" PINLET ; + ABORT: ; +ENDIF ; + +GREP: Thm :: STEP UP 'HISTORY-DATA' STEP UP 'TIMESTEP0000' STEP UP 'CHANNEL' + STEP AT 1 GETVAL 'EPSOUT' 1 1 1 >>EPSOUT<< ; + +EVALUATE DELTA := EPSOUT REFEPSVAL - REFEPSVAL / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFEPSVAL " CALCULATED=" EPSOUT ; + ABORT: ; +ENDIF ; + +ECHO "dfm_pdrop_40_sine completed" ; +END: ; \ No newline at end of file -- cgit v1.2.3