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/dfmp10s.c2m | 259 +++++++++++++++++++++++++ 1 file changed, 259 insertions(+) create mode 100755 Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp10s.c2m (limited to 'Donjon/data/thm_BWR_pdrop_dfm_proc/dfmp10s.c2m') 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 -- cgit v1.2.3