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 --- Dragon/data/tmatxs2_proc/TMATXS24.c2m | 119 ++++++++++++++++++++++++++++++++++ 1 file changed, 119 insertions(+) create mode 100644 Dragon/data/tmatxs2_proc/TMATXS24.c2m (limited to 'Dragon/data/tmatxs2_proc/TMATXS24.c2m') diff --git a/Dragon/data/tmatxs2_proc/TMATXS24.c2m b/Dragon/data/tmatxs2_proc/TMATXS24.c2m new file mode 100644 index 0000000..f4ba048 --- /dev/null +++ b/Dragon/data/tmatxs2_proc/TMATXS24.c2m @@ -0,0 +1,119 @@ +*---- +* TEST CASE TMATXS24 +* Electron Beam on an intraoperative breast cancer lumpectomy. +* NJOY-[CEPXSR]: electron 3MeV - 100 Groups - P2 non-coupled library. +* This is a highly simplified Mobetron Benchmark. Please, note that +* deterministic parameters are not physical here. +* Author: Ahmed Naceur +*---- +* Define STRUCTURES and MODULES used +*---- +MODULE GEO: LIB: MAC: SNT: ASM: FLU: HEAT: OUT: DELETE: + END: EDI: PSOUR: ; +LINKED_LIST MACRO DEPOS MICRO OLDMACRO GEOM TRACK SYSTEM FLUX + ESOUR EEDIT EDITION EEDITc GWT0 ; +SEQ_ASCII edep_mob.txt :: FILE 'edep_mob.txt' ; +PROCEDURE assertS ; + +* --- Geometry + +GEOM := GEO: :: CAR1D 4 + X- VOID X+ VOID + MIX 1 2 3 1 + MESHX 0.0 .6000 1.2000 1.4250 1.5000 + SPLITX 200 150 150 50 + ; + +* --- Tracking + +TRACK := SNT: GEOM :: + TITLE 'Intraoperative Mobetron Benchmark' + EDIT 0 MAXR 1000000 SCHM 1 DIAM 1 NLIVO SN 8 SCAT 8 EPSI 1.E-5 + MAXI 300 BFPG ; + +* --- Library access + +MICRO := LIB: :: + EDIT 1 + NMIX 3 + ANIS 3 + STERN 1 + MIXS LIB: MATXS2 FIL: _mobetron + * -- water -- * + MIX 1 300.0 1.0 NOGAS + H = h 0.111894 + O = o 0.888106 + * -- aluminium -- * + MIX 2 300.0 2.69890 NOGAS + Al = al 1.00 + * -- steel -- * + MIX 3 300.0 8.060 NOGAS + Fe = fe 0.641922 + Cr = cr 0.1837 + Ni = ni 0.1235 + Mo = mo 0.02386 + Mn = mn 0.01619 + Si = si 0.005093 + Cu = cu 0.001734 + W = w 0.0011 + Co = co 0.001052 + V = v 0.000508 + C = c 0.000445 + N = n 0.0004 + P = p 0.000200 + S = s 0.000175 + As = as 0.0001207 ; + +* --- Library structuring + +MACRO := MICRO :: STEP UP MACROLIB ; + +OLDMACRO := MACRO ; +MACRO := MAC: MACRO OLDMACRO :: + MIX 1 1 OLDL + MIX 2 2 OLDL + MIX 3 3 OLDL ; + +OLDMACRO := DELETE: OLDMACRO ; + +* ---- Electron beam definition + +ESOUR := PSOUR: MACRO TRACK GEOM :: EDIT 1 MONO 1 X- INTG 1 1.0 + DIR 1.0 0.0 0.0 ; + +* ---- Electron flux + +SYSTEM := ASM: MACRO TRACK :: EDIT 1 ARM ; +FLUX := FLU: MACRO TRACK SYSTEM ESOUR :: EDIT 1 TYPE S EXTE 100 ; + +* --- Recover flux and reaction rates + +EEDIT := OUT: FLUX TRACK MACRO GEOM :: EDIT 1 COND INTG NONE ; + +* --- Energy deposition from the e- source + +EEDITc := EEDIT ; +DEPOS := HEAT: EEDIT EEDITc :: EDIT 1 NORM 1.00 2.69890 8.060 ; + +* --- Save data structure for post-treatement + +edep_mob.txt := DEPOS ; + +* --- AssertS + +GWT0 := EEDIT :: STEP UP GROUP STEP AT 1 ; +assertS GWT0 :: 'FLUX-INTG' 1 3.311044E-03 ; +assertS GWT0 :: 'FLUX-INTG' 30 3.885913E-03 ; +assertS GWT0 :: 'FLUX-INTG' 60 4.248134E-03 ; +assertS GWT0 :: 'H-FACTOR' 1 1.504317E+00 ; +assertS GWT0 :: 'H-FACTOR' 30 1.609681E+00 ; +assertS GWT0 :: 'H-FACTOR' 60 1.639779E+00 ; + +* --- Delete data structure + +DEPOS EEDIT FLUX SYSTEM MACRO MICRO TRACK GEOM GWT0 := DELETE: + DEPOS EEDIT FLUX SYSTEM MACRO MICRO TRACK GEOM GWT0 ; + +* +ECHO "test TMATXS24 completed" ; +END: ; -- cgit v1.2.3