blob: f4ba048a83d432fb33d7c315d47f9459c74cc61e (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
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: ;
|
