sphinx.addnodesdocument)}( rawsourcechildren](docutils.nodestarget)}(h.. _CAAScapability:h]
attributes}(ids]classes]names]dupnames]backrefs]refidcaascapabilityutagnameh
lineKparenthhhsource8/Users/john/Documents/SCALE-test/docs/CAAScapability.rstubh section)}(hhh](h title)}(h"MAVRIC Appendix A: CAAS Capabilityh]h Text"MAVRIC Appendix A: CAAS Capability}(hh,h h*hhh!NhNubah}(h]h]h]h]h]uhh(h h%hhh!h"hKubh$)}(hhh](h))}(hIntroductionh]h/Introduction}(hh?h h=hhh!NhNubah}(h]h]h]h]h]uhh(h h:hhh!h"hKubh paragraph)}(hX_Modeling criticality accident alarm systems (CAAS) presents challenges
since the analysis consists of both a criticality problem and a
deep-penetration shielding problem :cite:`peplow_criticality_2009,peplow_criticality_2010`. Modern codes
are typically optimized to handle one of those types of problems but
usually not both. The two problems also differ in scale—the criticality
problem depends on materials relatively close to the fissionable
materials while the shielding problem can cover a much larger range.
SCALE now contains fully three-dimensional tools to perform both parts
of a CAAS analysis.h](h/Modeling criticality accident alarm systems (CAAS) presents challenges
since the analysis consists of both a criticality problem and a
deep-penetration shielding problem }(hModeling criticality accident alarm systems (CAAS) presents challenges
since the analysis consists of both a criticality problem and a
deep-penetration shielding problem h hMhhh!NhNubhpending_xref)}(hpeplow_criticality_2009h]h inline)}(hhZh]h/[peplow_criticality_2009]}(hhh h^ubah}(h]h]h]h]h]uhh\h hXubah}(h]id1ah]bibtexah]h]h] refdomaincitationreftyperef reftargethZrefwarnsupport_smartquotesuhhVh!h"hK h hMhhubhW)}(hpeplow_criticality_2010h]h])}(hh|h]h/[peplow_criticality_2010]}(hhh h~ubah}(h]h]h]h]h]uhh\h hzubah}(h]id2ah]hoah]h]h] refdomainhtreftypehv reftargeth|refwarnsupport_smartquotesuhhVh!h"hK h hMhhubh/X~. Modern codes
are typically optimized to handle one of those types of problems but
usually not both. The two problems also differ in scale—the criticality
problem depends on materials relatively close to the fissionable
materials while the shielding problem can cover a much larger range.
SCALE now contains fully three-dimensional tools to perform both parts
of a CAAS analysis.}(hX~. Modern codes
are typically optimized to handle one of those types of problems but
usually not both. The two problems also differ in scale—the criticality
problem depends on materials relatively close to the fissionable
materials while the shielding problem can cover a much larger range.
SCALE now contains fully three-dimensional tools to perform both parts
of a CAAS analysis.h hMhhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hK h h:hhubhL)}(hX]CAAS analysis can be performed with SCALE using the KENO-VI criticality
code and the MAVRIC shielding sequence. First, the fission distribution
(in space and energy) is determined by KENO-VI. This distribution is
saved to a file using a user-specified three-dimensional mesh grid and
an energy structure from the cross section library (or a user-defined
energy structure). MAVRIC then uses the fission distribution as the
source for a shielding calculation. MAVRIC is designed to implement
advanced variance reduction methods to calculate dose rates or detector
responses for difficult shielding problems.h]h/X]CAAS analysis can be performed with SCALE using the KENO-VI criticality
code and the MAVRIC shielding sequence. First, the fission distribution
(in space and energy) is determined by KENO-VI. This distribution is
saved to a file using a user-specified three-dimensional mesh grid and
an energy structure from the cross section library (or a user-defined
energy structure). MAVRIC then uses the fission distribution as the
source for a shielding calculation. MAVRIC is designed to implement
advanced variance reduction methods to calculate dose rates or detector
responses for difficult shielding problems.}(hhh hhhh!NhNubah}(h]h]h]h]h]uhhKh!h"hKh h:hhubhL)}(hX2For different types of shielding and different combinations of sources
and detector locations, different strategies can be used within
SCALE :cite:`peplow_methods_2012,wetzel_caas_2012`. Due to the way that cross sections for
neutron reactions that create photons are stored in ENDF, some of the
parameters used in the CAAS capability have changed since SCALE 6.1.
Please be sure to understand and follow the new guidance for the correct
accounting of secondary photons from neutron
reactions :cite:`miller_guidance_2013,miller_guide_2013,miller_corrected_2013`.h](h/For different types of shielding and different combinations of sources
and detector locations, different strategies can be used within
SCALE }(hFor different types of shielding and different combinations of sources
and detector locations, different strategies can be used within
SCALE h hhhh!NhNubhW)}(hpeplow_methods_2012h]h])}(hhh]h/[peplow_methods_2012]}(hhh hubah}(h]h]h]h]h]uhh\h hubah}(h]id3ah]hoah]h]h] refdomainhtreftypehv reftargethrefwarnsupport_smartquotesuhhVh!h"hKh hhhubhW)}(hwetzel_caas_2012h]h])}(hhh]h/[wetzel_caas_2012]}(hhh hubah}(h]h]h]h]h]uhh\h hubah}(h]id4ah]hoah]h]h] refdomainhtreftypehv reftargeth،refwarnsupport_smartquotesuhhVh!h"hKh hhhubh/X4. Due to the way that cross sections for
neutron reactions that create photons are stored in ENDF, some of the
parameters used in the CAAS capability have changed since SCALE 6.1.
Please be sure to understand and follow the new guidance for the correct
accounting of secondary photons from neutron
reactions }(hX4. Due to the way that cross sections for
neutron reactions that create photons are stored in ENDF, some of the
parameters used in the CAAS capability have changed since SCALE 6.1.
Please be sure to understand and follow the new guidance for the correct
accounting of secondary photons from neutron
reactions h hhhh!NhNubhW)}(hmiller_guidance_2013h]h])}(hhh]h/[miller_guidance_2013]}(hhh hubah}(h]h]h]h]h]uhh\h hubah}(h]id5ah]hoah]h]h] refdomainhtreftypehv reftargethrefwarnsupport_smartquotesuhhVh!h"hKh hhhubhW)}(hmiller_guide_2013h]h])}(hjh]h/[miller_guide_2013]}(hhh jubah}(h]h]h]h]h]uhh\h jubah}(h]id6ah]hoah]h]h] refdomainhtreftypehv reftargetjrefwarnsupport_smartquotesuhhVh!h"hKh hhhubhW)}(hmiller_corrected_2013h]h])}(hj4h]h/[miller_corrected_2013]}(hhh j6ubah}(h]h]h]h]h]uhh\h j2ubah}(h]id7ah]hoah]h]h] refdomainhtreftypehv reftargetj4refwarnsupport_smartquotesuhhVh!h"hKh hhhubh/.}(h.h hhhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hKh h:hhubeh}(h]introductionah]h]introductionah]h]uhh#h h%hhh!h"hKubh$)}(hhh](h))}(hMethodsh]h/Methods}(hjgh jehhh!NhNubah}(h]h]h]h]h]uhh(h jbhhh!h"hK'ubhL)}(hXPThe CAAS capability in SCALE is a two-step approach using KENO-VI and
MAVRIC. The first step is the determination of the source distribution,
done with the CSAS6 sequence which uses the KENO-VI functional module.
Along with calculating the system k\ :sub:`eff`, KENO-VI has been modified to
now accumulate the fission distribution over the non-skip generations.
This information is collected on a three-dimensional Cartesian mesh that
overlays the physical geometry model and is saved as a Monaco mesh tally
file. A utility program is used to convert the mesh tally into a Monaco
mesh source.h](h/The CAAS capability in SCALE is a two-step approach using KENO-VI and
MAVRIC. The first step is the determination of the source distribution,
done with the CSAS6 sequence which uses the KENO-VI functional module.
Along with calculating the system k }(hThe CAAS capability in SCALE is a two-step approach using KENO-VI and
MAVRIC. The first step is the determination of the source distribution,
done with the CSAS6 sequence which uses the KENO-VI functional module.
Along with calculating the system k\ h jshhh!NhNubh subscript)}(h
:sub:`eff`h]h/eff}(hhh j~ubah}(h]h]h]h]h]uhj|h jsubh/XL, KENO-VI has been modified to
now accumulate the fission distribution over the non-skip generations.
This information is collected on a three-dimensional Cartesian mesh that
overlays the physical geometry model and is saved as a Monaco mesh tally
file. A utility program is used to convert the mesh tally into a Monaco
mesh source.}(hXL, KENO-VI has been modified to
now accumulate the fission distribution over the non-skip generations.
This information is collected on a three-dimensional Cartesian mesh that
overlays the physical geometry model and is saved as a Monaco mesh tally
file. A utility program is used to convert the mesh tally into a Monaco
mesh source.h jshhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hK)h jbhhubhL)}(hXThe mesh source is then used in the second step as the source term in
MAVRIC. The absolute source strength is set by the user based on the
total number of fissions (based on the total power released) during the
criticality excursion. Further neutron multiplication should be
prevented in the MAVRIC transport calculation. Because the fission
neutrons have already been accounted for in the KENO-VI calculation,
failure to suppress neutron multiplication in the MARVIC sequence would
lead to incorrect flux estimates. In addition, if further fissions were
allowed, Monaco would add neutrons to its particle bank faster than they
could be removed (since the system is at or above critical) and the
simulation may never finish.h]h/XThe mesh source is then used in the second step as the source term in
MAVRIC. The absolute source strength is set by the user based on the
total number of fissions (based on the total power released) during the
criticality excursion. Further neutron multiplication should be
prevented in the MAVRIC transport calculation. Because the fission
neutrons have already been accounted for in the KENO-VI calculation,
failure to suppress neutron multiplication in the MARVIC sequence would
lead to incorrect flux estimates. In addition, if further fissions were
allowed, Monaco would add neutrons to its particle bank faster than they
could be removed (since the system is at or above critical) and the
simulation may never finish.}(hjh jhhh!NhNubah}(h]h]h]h]h]uhhKh!h"hK3h jbhhubhL)}(hXFor the transport part of a CAAS analysis, MAVRIC can be optimized to
calculate one specific detector response at one location using CADIS or
can be optimized to calculate multiple responses/locations with roughly
the same relative uncertainty using FW-CADIS. For calculating mesh
tallies of fluxes or dose rates, MAVRIC also uses FW-CADIS to help
balance the Monaco Monte Carlo calculation such that low flux voxels are
computed with about the same relative uncertainty as high flux voxels.h]h/XFor the transport part of a CAAS analysis, MAVRIC can be optimized to
calculate one specific detector response at one location using CADIS or
can be optimized to calculate multiple responses/locations with roughly
the same relative uncertainty using FW-CADIS. For calculating mesh
tallies of fluxes or dose rates, MAVRIC also uses FW-CADIS to help
balance the Monaco Monte Carlo calculation such that low flux voxels are
computed with about the same relative uncertainty as high flux voxels.}(hjh jhhh!NhNubah}(h]h]h]h]h]uhhKh!h"hK?h jbhhubhL)}(hXWith this two-step approach, users will have a great deal of flexibility
in modeling CAAS problems. The CSAS6 step and the MAVRIC step could both
use the same geometry and materials definitions or could have different
levels of detail included in each. The fission source distribution from
one CSAS6 calculation could be used in a number of different MAVRIC
building/detector models, with each MAVRIC calculation being optimized
for a given type of detector.h]h/XWith this two-step approach, users will have a great deal of flexibility
in modeling CAAS problems. The CSAS6 step and the MAVRIC step could both
use the same geometry and materials definitions or could have different
levels of detail included in each. The fission source distribution from
one CSAS6 calculation could be used in a number of different MAVRIC
building/detector models, with each MAVRIC calculation being optimized
for a given type of detector.}(hjh jhhh!NhNubah}(h]h]h]h]h]uhhKh!h"hKGh jbhhubeh}(h]methodsah]h]methodsah]h]uhh#h h%hhh!h"hK'ubh$)}(hhh](h))}(h
User Inputh]h/
User Input}(hjh jhhh!NhNubah}(h]h]h]h]h]uhh(h jhhh!h"hKPubhL)}(hXThe user can create either one input file containing both the CSAS6 and
MAVRIC calculations or can create two input files—one for each sequence.
The materials and geometry for these two models could be the same but do
not have to be. For example, the CSAS6 sequence might only contain the
materials and geometry important to the criticality calculation. Note,
however, that the critical source geometry and materials should be
modeled identically in both problems. This allows users greater
flexibility in modeling their specific problems.h]h/XThe user can create either one input file containing both the CSAS6 and
MAVRIC calculations or can create two input files—one for each sequence.
The materials and geometry for these two models could be the same but do
not have to be. For example, the CSAS6 sequence might only contain the
materials and geometry important to the criticality calculation. Note,
however, that the critical source geometry and materials should be
modeled identically in both problems. This allows users greater
flexibility in modeling their specific problems.}(hjh jhhh!NhNubah}(h]h]h]h]h]uhhKh!h"hKRh jhhubh$)}(hhh](h))}(h
KENO-VI inputh]h/
KENO-VI input}(hjh jhhh!NhNubah}(h]h]h]h]h]uhh(h jhhh!h"hK\ubhL)}(hXFor the criticality problem, the only extra input a user needs to
supply is the keyword “cds=yes” in the parameter block and a spatial
mesh around all of the fissionable materials of the problem in its own
gridGeometry block. Standard input for KENO-VI is described in the
KENO-VI chapter and CSAS6 chapter. The mesh used for the fission
source distribution is input using the “read gridGeometry *id*\ ”
block, where *id* is an identification number for that grid. Note that
only one grid can be specified, but that may change in the future. The
cells of the mesh are specified in each dimension separately by either
(1) listing all of the planes bounding the cells (keyword “xplanes”
followed by an “end”), (2) using keyword “xLinear *n* *a* *b*\ ” to
specify *n* cells between *a* and *b*, or by (3) specifying the
minimum plane, the maximum plane, and how many cells to make in that
dimension (“xmin=”, “xmax=”, “numXCells=”). The keywords “xplanes” and
“xLinear” can be used together and multiple times. Similar keywords
are used for the *y*- and *z*-dimensions. An example CSAS6 input file
that collects the fission distribution information would be as
follows:h](h/XFor the criticality problem, the only extra input a user needs to
supply is the keyword “cds=yes” in the parameter block and a spatial
mesh around all of the fissionable materials of the problem in its own
gridGeometry block. Standard input for KENO-VI is described in the
KENO-VI chapter and CSAS6 chapter. The mesh used for the fission
source distribution is input using the “read gridGeometry }(hXFor the criticality problem, the only extra input a user needs to
supply is the keyword “cds=yes” in the parameter block and a spatial
mesh around all of the fissionable materials of the problem in its own
gridGeometry block. Standard input for KENO-VI is described in the
KENO-VI chapter and CSAS6 chapter. The mesh used for the fission
source distribution is input using the “read gridGeometry h jhhh!NhNubh emphasis)}(h*id*h]h/id}(hhh jubah}(h]h]h]h]h]uhjh jubh/ ”
block, where }(h\ ”
block, where h jhhh!NhNubj)}(h*id*h]h/id}(hhh jubah}(h]h]h]h]h]uhjh jubh/XD is an identification number for that grid. Note that
only one grid can be specified, but that may change in the future. The
cells of the mesh are specified in each dimension separately by either
(1) listing all of the planes bounding the cells (keyword “xplanes”
followed by an “end”), (2) using keyword “xLinear }(hXD is an identification number for that grid. Note that
only one grid can be specified, but that may change in the future. The
cells of the mesh are specified in each dimension separately by either
(1) listing all of the planes bounding the cells (keyword “xplanes”
followed by an “end”), (2) using keyword “xLinear h jhhh!NhNubj)}(h*n*h]h/n}(hhh j*ubah}(h]h]h]h]h]uhjh jubh/ }(h h jhhh!NhNubj)}(h*a*h]h/a}(hhh j=ubah}(h]h]h]h]h]uhjh jubh/ }(hj<h jubj)}(h*b*h]h/b}(hhh jOubah}(h]h]h]h]h]uhjh jubh/ ” to
specify }(h\ ” to
specify h jhhh!NhNubj)}(h*n*h]h/n}(hhh jbubah}(h]h]h]h]h]uhjh jubh/ cells between }(h cells between h jhhh!NhNubj)}(h*a*h]h/a}(hhh juubah}(h]h]h]h]h]uhjh jubh/ and }(h and h jhhh!NhNubj)}(h*b*h]h/b}(hhh jubah}(h]h]h]h]h]uhjh jubh/X, or by (3) specifying the
minimum plane, the maximum plane, and how many cells to make in that
dimension (“xmin=”, “xmax=”, “numXCells=”). The keywords “xplanes” and
“xLinear” can be used together and multiple times. Similar keywords
are used for the }(hX, or by (3) specifying the
minimum plane, the maximum plane, and how many cells to make in that
dimension (“xmin=”, “xmax=”, “numXCells=”). The keywords “xplanes” and
“xLinear” can be used together and multiple times. Similar keywords
are used for the h jhhh!NhNubj)}(h*y*h]h/y}(hhh jubah}(h]h]h]h]h]uhjh jubh/- and }(h- and h jhhh!NhNubj)}(h*z*h]h/z}(hhh jubah}(h]h]h]h]h]uhjh jubh/p-dimensions. An example CSAS6 input file
that collects the fission distribution information would be as
follows:}(hp-dimensions. An example CSAS6 input file
that collects the fission distribution information would be as
follows:h jhhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hK^h jhhubh
literal_block)}(hXc=csas6
CAAS Example
v7-238n
read composition
…
end composition
read parameters
…
cds=yes
end parameters
read geometry
…
end geometry
read gridGeometry 1
title="Mesh for Collecting Fission Source"
xLinear 13 0.0 78.0
yplanes 0 8 16 24 32 34 36 38 40 48 56 64 72 end
zLinear 10 -2.54 77.46
end gridGeometry
end data
endh]h/Xc=csas6
CAAS Example
v7-238n
read composition
…
end composition
read parameters
…
cds=yes
end parameters
read geometry
…
end geometry
read gridGeometry 1
title="Mesh for Collecting Fission Source"
xLinear 13 0.0 78.0
yplanes 0 8 16 24 32 34 36 38 40 48 56 64 72 end
zLinear 10 -2.54 77.46
end gridGeometry
end data
end}(hhh jubah}(h]h]h]h]h]forcehighlight_args} xml:spacepreservelanguagescaleuhjh!h"hKrh jhhubhL)}(hXCThe fission source distribution collected by KENO-VI is saved to a
Monaco mesh tally file and copied back to the home area with the name
“\ *problemName*.fissionSource.3dmap”. This file can be viewed with the
Mesh File Viewer capability of Fulcrum that comes with SCALE. Note that
the finer the mesh spacing is the more generations/histories will have
to be simulated by the criticality calculation in order to reduce the
stochastic uncertainty in each mesh voxel of the distribution.
Regardless of the mesh size, creation of a fission mesh source file will
take more iterations than the number required to find k\ :sub:`eff`. KENO-VI
also saves the value of the system the average number of neutrons per
fission, in a file called “\ *problemName*.kenoNuBar.txt”. This value is
needed to properly determine the source strength.h](h/The fission source distribution collected by KENO-VI is saved to a
Monaco mesh tally file and copied back to the home area with the name
“ }(hThe fission source distribution collected by KENO-VI is saved to a
Monaco mesh tally file and copied back to the home area with the name
“\ h jhhh!NhNubj)}(h
*problemName*h]h/problemName}(hhh jubah}(h]h]h]h]h]uhjh jubh/X.fissionSource.3dmap”. This file can be viewed with the
Mesh File Viewer capability of Fulcrum that comes with SCALE. Note that
the finer the mesh spacing is the more generations/histories will have
to be simulated by the criticality calculation in order to reduce the
stochastic uncertainty in each mesh voxel of the distribution.
Regardless of the mesh size, creation of a fission mesh source file will
take more iterations than the number required to find k }(hX.fissionSource.3dmap”. This file can be viewed with the
Mesh File Viewer capability of Fulcrum that comes with SCALE. Note that
the finer the mesh spacing is the more generations/histories will have
to be simulated by the criticality calculation in order to reduce the
stochastic uncertainty in each mesh voxel of the distribution.
Regardless of the mesh size, creation of a fission mesh source file will
take more iterations than the number required to find k\ h jhhh!NhNubj})}(h
:sub:`eff`h]h/eff}(hhh jubah}(h]h]h]h]h]uhj|h jubh/o. KENO-VI
also saves the value of the system the average number of neutrons per
fission, in a file called “ }(ho. KENO-VI
also saves the value of the system the average number of neutrons per
fission, in a file called “\ h jhhh!NhNubj)}(h
*problemName*h]h/problemName}(hhh j
ubah}(h]h]h]h]h]uhjh jubh/R.kenoNuBar.txt”. This value is
needed to properly determine the source strength.}(hR.kenoNuBar.txt”. This value is
needed to properly determine the source strength.h jhhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hKh jhhubeh}(h]
keno-vi-inputah]h]
keno-vi inputah]h]uhh#h jhhh!h"hK\ubh$)}(hhh](h))}(h$Mesh tally to mesh source conversionh]h/$Mesh tally to mesh source conversion}(hj3h j1hhh!NhNubah}(h]h]h]h]h]uhh(h j.hhh!h"hKubhL)}(hA utility program is used to convert the Monaco mesh tally file into a
Monaco mesh source file. It can be part of the CSAS6 input file. The
user then needs to copy the resulting \*.msm file back to his home area.h]h/A utility program is used to convert the Monaco mesh tally file into a
Monaco mesh source file. It can be part of the CSAS6 input file. The
user then needs to copy the resulting *.msm file back to his home area.}(hA utility program is used to convert the Monaco mesh tally file into a
Monaco mesh source file. It can be part of the CSAS6 input file. The
user then needs to copy the resulting \*.msm file back to his home area.h j?hhh!NhNubah}(h]h]h]h]h]uhhKh!h"hKh j.hhubj)}(hX=csas6
…
end
=mt2msm
'fissionSource.3dmap' ! existing Keno fission source mesh tally
1 ! which family (for Keno files, there is only 1)
-1 ! use the whole family (keep all energy groups)
1 ! particle type for *.msm file (1-neutron, 2-photon)
'fissionSource.msm' ! name of newly created mesh source map file
end
=shell
copy fissionSource.msm “C:\mydocu~1\caasExample”
endh]h/X=csas6
…
end
=mt2msm
'fissionSource.3dmap' ! existing Keno fission source mesh tally
1 ! which family (for Keno files, there is only 1)
-1 ! use the whole family (keep all energy groups)
1 ! particle type for *.msm file (1-neutron, 2-photon)
'fissionSource.msm' ! name of newly created mesh source map file
end
=shell
copy fissionSource.msm “C:\mydocu~1\caasExample”
end}(hhh jNubah}(h]h]h]h]h]forcehighlight_args}jjjscaleuhjh!h"hKh j.hhubhL)}(haDetails on the conversion utility program are contained in :ref:`appendixb` of
the MAVRIC manual.h](h/;Details on the conversion utility program are contained in }(h;Details on the conversion utility program are contained in h j`hhh!NhNubhW)}(h:ref:`appendixb`h]h])}(hjkh]h/ appendixb}(hhh jmubah}(h]h](xrefstdstd-refeh]h]h]uhh\h jiubah}(h]h]h]h]h]refdocCAAScapability refdomainjxreftyperefrefexplicitrefwarn reftarget appendixbuhhVh!h"hKh j`ubh/ of
the MAVRIC manual.}(h of
the MAVRIC manual.h j`hhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hKh j.hhubhL)}(hXIn SCALE 6.1, the fission source distribution mesh tally produced by
KENO contained data representing the number of fissions in each mesh
cell in each energy group. In SCALE 6.2, the data stored was changed to
be the fissions per unit volume – the fission density. This is more
consistent with other mesh tallies from Monaco which store flux or dose
rates that represent averages over the mesh cells. This change also
allows the Mesh File Viewer to display the KENO fission source
distribution better. The mt2msm utility program also changed from SCALE
6.1 to SCALE 6.2 to account for the change in what is stored in the Keno
mesh tally file. Therefore, **KENO-produced fission source mesh tallies
and the mt2msm utility should not be mixed-and-matched across versions
of SCALE**. Doing so would result in the final Monaco mesh source file
being improperly normalized, which would not represent the KENO fission
source distribution and would give incorrect results in subsequent
MAVRIC calculations. Because there is not a specific ‘version flag’ in a
mesh tally file or mesh source map file, the user must ensure that they
have used the same version of SCALE for both the CSAS6 and MAVRIC
sequences any time the CAAS capability is used.h](h/XIn SCALE 6.1, the fission source distribution mesh tally produced by
KENO contained data representing the number of fissions in each mesh
cell in each energy group. In SCALE 6.2, the data stored was changed to
be the fissions per unit volume – the fission density. This is more
consistent with other mesh tallies from Monaco which store flux or dose
rates that represent averages over the mesh cells. This change also
allows the Mesh File Viewer to display the KENO fission source
distribution better. The mt2msm utility program also changed from SCALE
6.1 to SCALE 6.2 to account for the change in what is stored in the Keno
mesh tally file. Therefore, }(hXIn SCALE 6.1, the fission source distribution mesh tally produced by
KENO contained data representing the number of fissions in each mesh
cell in each energy group. In SCALE 6.2, the data stored was changed to
be the fissions per unit volume – the fission density. This is more
consistent with other mesh tallies from Monaco which store flux or dose
rates that represent averages over the mesh cells. This change also
allows the Mesh File Viewer to display the KENO fission source
distribution better. The mt2msm utility program also changed from SCALE
6.1 to SCALE 6.2 to account for the change in what is stored in the Keno
mesh tally file. Therefore, h jhhh!NhNubh strong)}(h}**KENO-produced fission source mesh tallies
and the mt2msm utility should not be mixed-and-matched across versions
of SCALE**h]h/yKENO-produced fission source mesh tallies
and the mt2msm utility should not be mixed-and-matched across versions
of SCALE}(hhh jubah}(h]h]h]h]h]uhjh jubh/X. Doing so would result in the final Monaco mesh source file
being improperly normalized, which would not represent the KENO fission
source distribution and would give incorrect results in subsequent
MAVRIC calculations. Because there is not a specific ‘version flag’ in a
mesh tally file or mesh source map file, the user must ensure that they
have used the same version of SCALE for both the CSAS6 and MAVRIC
sequences any time the CAAS capability is used.}(hX. Doing so would result in the final Monaco mesh source file
being improperly normalized, which would not represent the KENO fission
source distribution and would give incorrect results in subsequent
MAVRIC calculations. Because there is not a specific ‘version flag’ in a
mesh tally file or mesh source map file, the user must ensure that they
have used the same version of SCALE for both the CSAS6 and MAVRIC
sequences any time the CAAS capability is used.h jhhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hKh j.hhubeh}(h]$mesh-tally-to-mesh-source-conversionah]h]$mesh tally to mesh source conversionah]h]uhh#h jhhh!h"hKubh$)}(hhh](h))}(hMAVRIC inputh]h/MAVRIC input}(hjh jhhh!NhNubah}(h]h]h]h]h]uhh(h jhhh!h"hKubhL)}(hXThe input for the MAVRIC portion of the CAAS problem should include the
materials and geometry of the criticality model, use the fission
distribution as a source, set the source strength, and set any optional
modifiers to the source to change its location or add fission photons.
The cross section library used by the MAVRIC calculation does not need
to have the same group structure as the fission distribution. MAVRIC
will automatically convert the fission source group structure to match
the group structure of its cross section library.h]h/XThe input for the MAVRIC portion of the CAAS problem should include the
materials and geometry of the criticality model, use the fission
distribution as a source, set the source strength, and set any optional
modifiers to the source to change its location or add fission photons.
The cross section library used by the MAVRIC calculation does not need
to have the same group structure as the fission distribution. MAVRIC
will automatically convert the fission source group structure to match
the group structure of its cross section library.}(hjh jhhh!NhNubah}(h]h]h]h]h]uhhKh!h"hKh jhhubhL)}(hXpThe shielding calculation needs to specify that the source is the
fission distribution file, which is typically “fissionSource.msm”. The
total source strength can be specified by either the number of fissions
in the criticality accident (fission rate or total number) or by the
number of released neutrons (the fission rate multiplied by :math:`\over{v}` per
fission). The value of :math:`\over{v}` will be read from the file “kenoNuBar.txt” in the
SCALE temporary directory if it is not given in the source input with
the keyword “nu-bar=”. The mesh source can also be placed at different
coordinates in the geometry using the “origin x=\ *x* y=\ *y* z=\ *z*\ ”
keywords, if a different reference frame was used with the criticality
geometry model that created the mesh source. Rotations of mesh sources
are not available at this time. It is also recommended to use filters in
the source block to define the source, such as the ‘mixture=’ filter to
only allow source sampling from a specific mixture since the mesh source
can be transformed from it’s original origin or meshes can cover
non-fissionable materials.h](h/XVThe shielding calculation needs to specify that the source is the
fission distribution file, which is typically “fissionSource.msm”. The
total source strength can be specified by either the number of fissions
in the criticality accident (fission rate or total number) or by the
number of released neutrons (the fission rate multiplied by }(hXVThe shielding calculation needs to specify that the source is the
fission distribution file, which is typically “fissionSource.msm”. The
total source strength can be specified by either the number of fissions
in the criticality accident (fission rate or total number) or by the
number of released neutrons (the fission rate multiplied by h jhhh!NhNubh math)}(h:math:`\over{v}`h]h/\over{v}}(hhh jubah}(h]h]h]h]h]uhjh jubh/ per
fission). The value of }(h per
fission). The value of h jhhh!NhNubj)}(h:math:`\over{v}`h]h/\over{v}}(hhh jubah}(h]h]h]h]h]uhjh jubh/ will be read from the file “kenoNuBar.txt” in the
SCALE temporary directory if it is not given in the source input with
the keyword “nu-bar=”. The mesh source can also be placed at different
coordinates in the geometry using the “origin x= }(h will be read from the file “kenoNuBar.txt” in the
SCALE temporary directory if it is not given in the source input with
the keyword “nu-bar=”. The mesh source can also be placed at different
coordinates in the geometry using the “origin x=\ h jhhh!NhNubj)}(h*x*h]h/x}(hhh jubah}(h]h]h]h]h]uhjh jubh/ y= }(h y=\ h jhhh!NhNubj)}(h*y*h]h/y}(hhh j&ubah}(h]h]h]h]h]uhjh jubh/ z= }(h z=\ h jhhh!NhNubj)}(h*z*h]h/z}(hhh j9ubah}(h]h]h]h]h]uhjh jubh/X ”
keywords, if a different reference frame was used with the criticality
geometry model that created the mesh source. Rotations of mesh sources
are not available at this time. It is also recommended to use filters in
the source block to define the source, such as the ‘mixture=’ filter to
only allow source sampling from a specific mixture since the mesh source
can be transformed from it’s original origin or meshes can cover
non-fissionable materials.}(hX\ ”
keywords, if a different reference frame was used with the criticality
geometry model that created the mesh source. Rotations of mesh sources
are not available at this time. It is also recommended to use filters in
the source block to define the source, such as the ‘mixture=’ filter to
only allow source sampling from a specific mixture since the mesh source
can be transformed from it’s original origin or meshes can cover
non-fissionable materials.h jhhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hKh jhhubhL)}(hFor example, using a KENO-VI fission distribution, placing it somewhere
in the MAVRIC model and setting the source strength (in neutrons/s) to
correspond to 10\ :sup:`17` fission/s would look likeh](h/For example, using a KENO-VI fission distribution, placing it somewhere
in the MAVRIC model and setting the source strength (in neutrons/s) to
correspond to 10 }(hFor example, using a KENO-VI fission distribution, placing it somewhere
in the MAVRIC model and setting the source strength (in neutrons/s) to
correspond to 10\ h jRhhh!NhNubh superscript)}(h :sup:`17`h]h/17}(hhh j]ubah}(h]h]h]h]h]uhj[h jRubh/ fission/s would look like}(h fission/s would look likeh jRhhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hKh jhhubj)}(hX'=shell
copy “C:\mydocu~1\caasExample\kenoInput.kenoNuBar.txt” kenoNuBar.txt
end
=mavric
…
read sources
src 1
meshSourceFile “C:\mydocu~1\caasExample\fissionSource.msm"
origin x=600 y=650 z=400
fissions=1.0e17
end src
end sources
…
end data
endh]h/X'=shell
copy “C:\mydocu~1\caasExample\kenoInput.kenoNuBar.txt” kenoNuBar.txt
end
=mavric
…
read sources
src 1
meshSourceFile “C:\mydocu~1\caasExample\fissionSource.msm"
origin x=600 y=650 z=400
fissions=1.0e17
end src
end sources
…
end data
end}(hhh jvubah}(h]h]h]h]h]forcehighlight_args}jjjscaleuhjh!h"hKh jhhubhL)}(hX_The source strength in neutrons/s will be calculated by MAVRIC to be the
fission rate multiplied by the value of read from the “kenoNuBar.txt”
file. The neutron strength could have alternatively been specified using
the standard source strength keyword “strength=2.5e17” (for an example
with the average number of neutrons per fission of 2.5).h]h/X_The source strength in neutrons/s will be calculated by MAVRIC to be the
fission rate multiplied by the value of read from the “kenoNuBar.txt”
file. The neutron strength could have alternatively been specified using
the standard source strength keyword “strength=2.5e17” (for an example
with the average number of neutrons per fission of 2.5).}(hjh jhhh!NhNubah}(h]h]h]h]h]uhhKh!h"hKh jhhubhL)}(hX The Monte Carlo functional module used by MAVRIC, Monaco, is a
fixed-source code. Unless told otherwise, neutrons will multiply in
fissionable materials. Since all of the neutrons were part of the
source, neutron multiplication should not be allowed and MAVRIC should
be run with the keyword “fissionMult=0” in the parameters block. For
systems at or near critical without the “fissionMult=0” keyword, Monaco
simulations may not end since neutrons will be added to the particle
bank at the same rate they leave the system or get killed.h]h/X The Monte Carlo functional module used by MAVRIC, Monaco, is a
fixed-source code. Unless told otherwise, neutrons will multiply in
fissionable materials. Since all of the neutrons were part of the
source, neutron multiplication should not be allowed and MAVRIC should
be run with the keyword “fissionMult=0” in the parameters block. For
systems at or near critical without the “fissionMult=0” keyword, Monaco
simulations may not end since neutrons will be added to the particle
bank at the same rate they leave the system or get killed.}(hjh jhhh!NhNubah}(h]h]h]h]h]uhhKh!h"hMh jhhubhL)}(hX
The shielding calculation can be run using standard variance reduction
methods (such as path length stretching, user-defined weight windows
based on geometry regions, and user-defined source biasing) or using the
automated tools which employ approximate discrete-ordinates calculations
to determine the space/energy weight windows as well as a biased source
distribution in space and energy. The automated tools can be used to
optimize the shielding calculation to determine one specific tally using
CADIS or several separate tallies or a mesh tally over a large volume of
the problem space using FW-CADIS. When using these advanced variance
reduction methods, remember to include planes in the discrete‑ordinates
mesh definition that correspond to the planes in the fission
distribution that the source is based on. If they are not included,
MAVRIC will resample the fission source on the discrete-ordinates mesh
it is using for the importance map, possibly smearing or reducing the
original resolution of the fission distribution.h]h/X
The shielding calculation can be run using standard variance reduction
methods (such as path length stretching, user-defined weight windows
based on geometry regions, and user-defined source biasing) or using the
automated tools which employ approximate discrete-ordinates calculations
to determine the space/energy weight windows as well as a biased source
distribution in space and energy. The automated tools can be used to
optimize the shielding calculation to determine one specific tally using
CADIS or several separate tallies or a mesh tally over a large volume of
the problem space using FW-CADIS. When using these advanced variance
reduction methods, remember to include planes in the discrete‑ordinates
mesh definition that correspond to the planes in the fission
distribution that the source is based on. If they are not included,
MAVRIC will resample the fission source on the discrete-ordinates mesh
it is using for the importance map, possibly smearing or reducing the
original resolution of the fission distribution.}(hjh jhhh!NhNubah}(h]h]h]h]h]uhhKh!h"hM
h jhhubeh}(h]mavric-inputah]h]mavric inputah]h]uhh#h jhhh!h"hKubeh}(h]
user-inputah]h]
user inputah]h]uhh#h h%hhh!h"hKPubh$)}(hhh](h))}(hExample problemh]h/Example problem}(hjh jhhh!NhNubah}(h]h]h]h]h]uhh(h jhhh!h"hMubhL)}(hXcConsider the Jezebel critical plutonium sphere experiment, shown in
:numref:`fig4a1`, taking place inside a simple fictitious building, shown
in :numref:`fig4a2`. The building has two rooms: an experiment room and a
control room. In the control room there is a criticality alarm detector,
and it is positioned furthest from the entry to the experiment room. For
this example, assume that a criticality excursion results in a total of
10\ :sup:`18` fissions. This example will calculate the neutron and
photon doses seen by a detector in the control room, as well as
calculate a dose map for the entire building.h](h/DConsider the Jezebel critical plutonium sphere experiment, shown in
}(hDConsider the Jezebel critical plutonium sphere experiment, shown in
h jhhh!NhNubhW)}(h:numref:`fig4a1`h]h literal)}(hjh]h/fig4a1}(hhh jubah}(h]h](jwstd
std-numrefeh]h]h]uhjh jubah}(h]h]h]h]h]refdocj refdomainjreftypenumrefrefexplicitrefwarnjfig4a1uhhVh!h"hMh jubh/=, taking place inside a simple fictitious building, shown
in }(h=, taking place inside a simple fictitious building, shown
in h jhhh!NhNubhW)}(h:numref:`fig4a2`h]j)}(hjh]h/fig4a2}(hhh jubah}(h]h](jwstd
std-numrefeh]h]h]uhjh jubah}(h]h]h]h]h]refdocj refdomainjreftypenumrefrefexplicitrefwarnjfig4a2uhhVh!h"hMh jubh/X. The building has two rooms: an experiment room and a
control room. In the control room there is a criticality alarm detector,
and it is positioned furthest from the entry to the experiment room. For
this example, assume that a criticality excursion results in a total of
10 }(hX. The building has two rooms: an experiment room and a
control room. In the control room there is a criticality alarm detector,
and it is positioned furthest from the entry to the experiment room. For
this example, assume that a criticality excursion results in a total of
10\ h jhhh!NhNubj\)}(h :sup:`18`h]h/18}(hhh j(ubah}(h]h]h]h]h]uhj[h jubh/ fissions. This example will calculate the neutron and
photon doses seen by a detector in the control room, as well as
calculate a dose map for the entire building.}(h fissions. This example will calculate the neutron and
photon doses seen by a detector in the control room, as well as
calculate a dose map for the entire building.h jhhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hMh jhhubh)}(h.. _fig4a1:h]h}(h]h]h]h]h]hfig4a1uhh
hM'h jhhh!h"ubh figure)}(hhh](h image)}(h_.. figure:: figs/MAVRIC/fig4a1.png
:align: center
:width: 60 %
Cutaway view of Jezebel.
h]h}(h]h]h]h]h]width60%urifigs/MAVRIC/fig4a1.png
candidates}*j`suhjQh jNh!h"hM-ubh caption)}(hCutaway view of Jezebel.h]h/Cutaway view of Jezebel.}(hjhh jfubah}(h]h]h]h]h]uhjdh!h"hM-h jNubeh}(h](id8jKeh]h]fig4a1ah]h]aligncenteruhjLhM-h jhhh!h"expect_referenced_by_name}jyjAsexpect_referenced_by_id}jKjAsubh)}(h.. _fig4a2:h]h}(h]h]h]h]h]hfig4a2uhh
hM/h jhhh!h"ubjM)}(hhh](jR)}(h`.. figure:: figs/MAVRIC/fig4a2.png
:align: center
:width: 60 %
Simple two-room building.
h]h}(h]h]h]h]h]width60%urifigs/MAVRIC/fig4a2.pngja}jcjsuhjQh jh!h"hM4ubje)}(hSimple two-room building.h]h/Simple two-room building.}(hjh jubah}(h]h]h]h]h]uhjdh!h"hM4h jubeh}(h](id9jeh]h]fig4a2ah]h]j|centeruhjLhM4h jhhh!h"j~}jjsj}jjsubh$)}(hhh](h))}(h3KENO-VI criticality and fission source distributionh]h/3KENO-VI criticality and fission source distribution}(hjh jhhh!NhNubah}(h]h]h]h]h]uhh(h jhhh!h"hM7ubhL)}(hXGFor the criticality calculation, consider just a bare sphere of
plutonium, with a radius of 6.38493 cm. Atom densities (atoms/b·cm) for
the sphere material are :sup:`239`\ Pu 0.037047; :sup:`240`\ Pu
0.0017512; :sup:`241`\ Pu 0.00011674; and Cu 0.0013752. This can be
easily modeled as a sphere at the origin. For collecting the fission
distribution, a uniform mesh grid can be constructed around the sphere,
extending 7 cm in each direction, with a 1×1×1 cm voxel size. The first
portion of the input files ``mavric.caasA.inp`` and ``mavric.caasB.inp`` looks
like the following:h](h/For the criticality calculation, consider just a bare sphere of
plutonium, with a radius of 6.38493 cm. Atom densities (atoms/b·cm) for
the sphere material are }(hFor the criticality calculation, consider just a bare sphere of
plutonium, with a radius of 6.38493 cm. Atom densities (atoms/b·cm) for
the sphere material are h jhhh!NhNubj\)}(h
:sup:`239`h]h/239}(hhh jubah}(h]h]h]h]h]uhj[h jubh/ Pu 0.037047; }(h\ Pu 0.037047; h jhhh!NhNubj\)}(h
:sup:`240`h]h/240}(hhh jubah}(h]h]h]h]h]uhj[h jubh/ Pu
0.0017512; }(h\ Pu
0.0017512; h jhhh!NhNubj\)}(h
:sup:`241`h]h/241}(hhh jubah}(h]h]h]h]h]uhj[h jubh/X! Pu 0.00011674; and Cu 0.0013752. This can be
easily modeled as a sphere at the origin. For collecting the fission
distribution, a uniform mesh grid can be constructed around the sphere,
extending 7 cm in each direction, with a 1×1×1 cm voxel size. The first
portion of the input files }(hX!\ Pu 0.00011674; and Cu 0.0013752. This can be
easily modeled as a sphere at the origin. For collecting the fission
distribution, a uniform mesh grid can be constructed around the sphere,
extending 7 cm in each direction, with a 1×1×1 cm voxel size. The first
portion of the input files h jhhh!NhNubj)}(h``mavric.caasA.inp``h]h/mavric.caasA.inp}(hhh jubah}(h]h]h]h]h]uhjh jubh/ and }(h and h jhhh!NhNubj)}(h``mavric.caasB.inp``h]h/mavric.caasB.inp}(hhh jubah}(h]h]h]h]h]uhjh jubh/ looks
like the following:}(h looks
like the following:h jhhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hM9h jhhubj)}(hX=csas6
Dose Rates from a Jezebel Accident in a Block Building
v7-238n
'-------------------------------------------------------------------------------
' Composition Block
'-------------------------------------------------------------------------------
read composition
Pu-239 1 0 0.037047 end
Pu-240 1 0 0.0017512 end
Pu-241 1 0 0.00011674 end
Cu 1 0 0.0013752 end
end composition
'-------------------------------------------------------------------------------
' Parameters Block
'-------------------------------------------------------------------------------
read parameters
gen=250 npg=200000 nsk=50 htm=no
cds=yes
end parameters
'-------------------------------------------------------------------------------
' Geometry Block - SCALE standard geometry package (SGGP)
'-------------------------------------------------------------------------------
read geometry
global unit 2
sphere 1 6.38493
media 1 1 1 vol=1090.3277
boundary 1
end geometry
'-------------------------------------------------------------------------------
' Grid Block
'-------------------------------------------------------------------------------
read gridGeometry 1
title="Mesh for Collecting Fission Distribution"
xLinear 14 -7.0 7.0
yLinear 14 -7.0 7.0
zLinear 14 -7.0 7.0
end gridGeometry
end data
end
=mt2msm
'fissionSource.3dmap'
1
-1
1
mavric.caas[A/B].fissionSource.msm'
endh]h/X=csas6
Dose Rates from a Jezebel Accident in a Block Building
v7-238n
'-------------------------------------------------------------------------------
' Composition Block
'-------------------------------------------------------------------------------
read composition
Pu-239 1 0 0.037047 end
Pu-240 1 0 0.0017512 end
Pu-241 1 0 0.00011674 end
Cu 1 0 0.0013752 end
end composition
'-------------------------------------------------------------------------------
' Parameters Block
'-------------------------------------------------------------------------------
read parameters
gen=250 npg=200000 nsk=50 htm=no
cds=yes
end parameters
'-------------------------------------------------------------------------------
' Geometry Block - SCALE standard geometry package (SGGP)
'-------------------------------------------------------------------------------
read geometry
global unit 2
sphere 1 6.38493
media 1 1 1 vol=1090.3277
boundary 1
end geometry
'-------------------------------------------------------------------------------
' Grid Block
'-------------------------------------------------------------------------------
read gridGeometry 1
title="Mesh for Collecting Fission Distribution"
xLinear 14 -7.0 7.0
yLinear 14 -7.0 7.0
zLinear 14 -7.0 7.0
end gridGeometry
end data
end
=mt2msm
'fissionSource.3dmap'
1
-1
1
mavric.caas[A/B].fissionSource.msm'
end}(hhh j7ubah}(h]h]h]h]h]jjforcejscalehighlight_args}uhjh!h"hMCh jhhubhL)}(hThe results of this 26 minute calculation are shown in :numref:`tab4a1`,
and details about the calculated fission distribution are shown in
:numref:`fig4a3` and :numref:`fig4a4`.h](h/7The results of this 26 minute calculation are shown in }(h7The results of this 26 minute calculation are shown in h jIhhh!NhNubhW)}(h:numref:`tab4a1`h]j)}(hjTh]h/tab4a1}(hhh jVubah}(h]h](jwstd
std-numrefeh]h]h]uhjh jRubah}(h]h]h]h]h]refdocj refdomainj`reftypenumrefrefexplicitrefwarnjtab4a1uhhVh!h"hMxh jIubh/E,
and details about the calculated fission distribution are shown in
}(hE,
and details about the calculated fission distribution are shown in
h jIhhh!NhNubhW)}(h:numref:`fig4a3`h]j)}(hjyh]h/fig4a3}(hhh j{ubah}(h]h](jwstd
std-numrefeh]h]h]uhjh jwubah}(h]h]h]h]h]refdocj refdomainjreftypenumrefrefexplicitrefwarnjfig4a3uhhVh!h"hMxh jIubh/ and }(h and h jIhhh!NhNubhW)}(h:numref:`fig4a4`h]j)}(hjh]h/fig4a4}(hhh jubah}(h]h](jwstd
std-numrefeh]h]h]uhjh jubah}(h]h]h]h]h]refdocj refdomainjreftypenumrefrefexplicitrefwarnjfig4a4uhhVh!h"hMxh jIubh/.}(hjSh jIhhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hMxh jhhubh)}(h.. _tab4a1:h]h}(h]h]h]h]h]htab4a1uhh
hM|h jhhh!h"ubh table)}(hhh](h))}(h Results of the CSAS6 calculationh]h/ Results of the CSAS6 calculation}(hjh jubah}(h]h]h]h]h]uhh(h!h"hM~h jubh tgroup)}(hhh](h colspec)}(hhh]h}(h]h]h]h]h]colwidthKuhjh jubj)}(hhh]h}(h]h]h]h]h]colwidthKuhjh jubj)}(hhh]h}(h]h]h]h]h]colwidthKuhjh jubj)}(hhh]h}(h]h]h]h]h]colwidthK
uhjh jubh thead)}(hhh]h row)}(hhh](h entry)}(hhh]hL)}(h**Quantity**h]j)}(hj$h]h/Quantity}(hhh j&ubah}(h]h]h]h]h]uhjh j"ubah}(h]h]h]h]h]uhhKh!h"hMh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h **Value**h]j)}(hjDh]h/Value}(hhh jFubah}(h]h]h]h]h]uhjh jBubah}(h]h]h]h]h]uhhKh!h"hMh j?ubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h**Uncertainty**h]j)}(hjdh]h/Uncertainty}(hhh jfubah}(h]h]h]h]h]uhjh jbubah}(h]h]h]h]h]uhhKh!h"hMh j_ubah}(h]h]h]h]h]uhjh jubj)}(hhh]h}(h]h]h]h]h]uhjh jubeh}(h]h]h]h]h]uhjh jubah}(h]h]h]h]h]uhjh jubh tbody)}(hhh](j)}(hhh](j)}(hhh]hL)}(h
k\ :sub:`eff`h](h/k }(hk\ h jubj})}(h
:sub:`eff`h]h/eff}(hhh jubah}(h]h]h]h]h]uhj|h jubeh}(h]h]h]h]h]uhhKh!h"hMh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(hbest estimate system k-effh]h/best estimate system k-eff}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h1.00024h]h/1.00024}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h0.00014h]h/0.00014}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMh jubah}(h]h]h]h]h]uhjh jubeh}(h]h]h]h]h]uhjh jubj)}(hhh](j)}(hhh]hL)}(h:math:`\over{v}`h]j)}(h:math:`\over{v}`h]h/\over{v}}(hhh jubah}(h]h]h]h]h]uhjh jubah}(h]h]h]h]h]uhhKh!h"hMh jubah}(h]h]h]h]h]uhjh j
ubj)}(hhh]hL)}(h
system nu barh]h/
system nu bar}(hj6h j4ubah}(h]h]h]h]h]uhhKh!h"hMh j1ubah}(h]h]h]h]h]uhjh j
ubj)}(hhh]hL)}(h3.15671h]h/3.15671}(hjMh jKubah}(h]h]h]h]h]uhhKh!h"hMh jHubah}(h]h]h]h]h]uhjh j
ubj)}(hhh]hL)}(h4.77938E-05h]h/4.77938E-05}(hjdh jbubah}(h]h]h]h]h]uhhKh!h"hMh j_ubah}(h]h]h]h]h]uhjh j
ubeh}(h]h]h]h]h]uhjh jubeh}(h]h]h]h]h]uhjh jubeh}(h]h]h]h]h]colsKuhjh jubeh}(h](id10jeh]h]tab4a1ah]h]j|defaultuhjh jhhh!h"hNj~}jjsj}jjsubh)}(h.. _fig4a3:h]h}(h]h]h]h]h]hfig4a3uhh
hMh jhhh!h"ubjM)}(hhh](jR)}(h.. figure:: figs/MAVRIC/fig4a3.png
:align: center
:width: 60 %
Fission source spatial distribution for the center horizontal slice.
h]h}(h]h]h]h]h]width60%urifigs/MAVRIC/fig4a3.pngja}jcjsuhjQh jh!h"hMubje)}(hDFission source spatial distribution for the center horizontal slice.h]h/DFission source spatial distribution for the center horizontal slice.}(hjh jubah}(h]h]h]h]h]uhjdh!h"hMh jubeh}(h](id11jeh]h]fig4a3ah]h]j|centeruhjLhMh jhhh!h"j~}jjsj}jjsubh)}(h.. _fig4a4:h]h}(h]h]h]h]h]hfig4a4uhh
hMh jhhh!h"ubjM)}(hhh](jR)}(h.. figure:: figs/MAVRIC/fig4a4.png
:align: center
:width: 60 %
Fission source energy distribution for the center voxel.
h]h}(h]h]h]h]h]width60%urifigs/MAVRIC/fig4a4.pngja}jcjsuhjQh jh!h"hMubje)}(h8Fission source energy distribution for the center voxel.h]h/8Fission source energy distribution for the center voxel.}(hjh jubah}(h]h]h]h]h]uhjdh!h"hMh jubeh}(h](id12jeh]h]fig4a4ah]h]j|centeruhjLhMh jhhh!h"j~}jjsj}jjsubeh}(h]3keno-vi-criticality-and-fission-source-distributionah]h]3keno-vi criticality and fission source distributionah]h]uhh#h jhhh!h"hM7ubh$)}(hhh](h))}(hMAVRIC transport calculationsh]h/MAVRIC transport calculations}(hj
h j hhh!NhNubah}(h]h]h]h]h]uhh(h j hhh!h"hMubhL)}(hXTwo MAVRIC calculations will be done—one that calculates the doses seen
at the detector and one that computes mesh tallies of doses over the
entire building. They will share the same materials, geometry, and
source but will have different tally and variance reduction options.h]h/XTwo MAVRIC calculations will be done—one that calculates the doses seen
at the detector and one that computes mesh tallies of doses over the
entire building. They will share the same materials, geometry, and
source but will have different tally and variance reduction options.}(hj h j hhh!NhNubah}(h]h]h]h]h]uhhKh!h"hMh j hhubhL)}(hXThe two-room building will be a simple model using concrete-block walls,
a concrete floor, and a steel roof, with dimensions shown in
:numref:`fig4a5`. The building exterior dimensions are 1200 cm long,
600 cm wide, and 300 cm high above the ground. The exterior and interior
walls are all made of a double layer of typical concrete blocks (total
of 40 cm thick.) Concrete blocks are typically 39×19×19 cm and weigh
~13.5 kg, since they have a volume fraction of 33.2%. The floor is made
of poured concrete, extending 60 cm into the ground. The roof and the
exterior door (120 cm wide and 210 cm tall) are made of 1/8 in. (0.3175
cm) thick steel. The experiment room on the left connects to the control
room on the right through a maze that prevents radiation streaming.
Assume that the critical experiment was in the center of the experiment
room, 100 cm above the floor. Assume the detector in the control room is
a 30 cm diameter sphere located at position (1145, 55, 285).h](h/The two-room building will be a simple model using concrete-block walls,
a concrete floor, and a steel roof, with dimensions shown in
}(hThe two-room building will be a simple model using concrete-block walls,
a concrete floor, and a steel roof, with dimensions shown in
h j' hhh!NhNubhW)}(h:numref:`fig4a5`h]j)}(hj2 h]h/fig4a5}(hhh j4 ubah}(h]h](jwstd
std-numrefeh]h]h]uhjh j0 ubah}(h]h]h]h]h]refdocj refdomainj> reftypenumrefrefexplicitrefwarnjfig4a5uhhVh!h"hMh j' ubh/X>. The building exterior dimensions are 1200 cm long,
600 cm wide, and 300 cm high above the ground. The exterior and interior
walls are all made of a double layer of typical concrete blocks (total
of 40 cm thick.) Concrete blocks are typically 39×19×19 cm and weigh
~13.5 kg, since they have a volume fraction of 33.2%. The floor is made
of poured concrete, extending 60 cm into the ground. The roof and the
exterior door (120 cm wide and 210 cm tall) are made of 1/8 in. (0.3175
cm) thick steel. The experiment room on the left connects to the control
room on the right through a maze that prevents radiation streaming.
Assume that the critical experiment was in the center of the experiment
room, 100 cm above the floor. Assume the detector in the control room is
a 30 cm diameter sphere located at position (1145, 55, 285).}(hX>. The building exterior dimensions are 1200 cm long,
600 cm wide, and 300 cm high above the ground. The exterior and interior
walls are all made of a double layer of typical concrete blocks (total
of 40 cm thick.) Concrete blocks are typically 39×19×19 cm and weigh
~13.5 kg, since they have a volume fraction of 33.2%. The floor is made
of poured concrete, extending 60 cm into the ground. The roof and the
exterior door (120 cm wide and 210 cm tall) are made of 1/8 in. (0.3175
cm) thick steel. The experiment room on the left connects to the control
room on the right through a maze that prevents radiation streaming.
Assume that the critical experiment was in the center of the experiment
room, 100 cm above the floor. Assume the detector in the control room is
a 30 cm diameter sphere located at position (1145, 55, 285).h j' hhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hMh j hhubh)}(h.. _fig4a5:h]h}(h]h]h]h]h]hfig4a5uhh
hMh j hhh!h"ubjM)}(hhh](jR)}(h.. figure:: figs/MAVRIC/fig4a5.png
:align: center
Coordinates of the floor, walls, ceiling, and door of the simple block building model (in cm).
h]h}(h]h]h]h]h]urifigs/MAVRIC/fig4a5.pngja}jcjt suhjQh jf h!h"hMubje)}(h^Coordinates of the floor, walls, ceiling, and door of the simple block building model (in cm).h]h/^Coordinates of the floor, walls, ceiling, and door of the simple block building model (in cm).}(hjx h jv ubah}(h]h]h]h]h]uhjdh!h"hMh jf ubeh}(h](id13je eh]h]fig4a5ah]h]j|centeruhjLhMh j hhh!h"j~}j j[ sj}je j[ subhL)}(hIf the MAVRIC transport calculation is not in the same file as the CSAS6 calculation, the MAVRIC input would start by moving the KENO-VI results into the SCALE temporary area:h]h/If the MAVRIC transport calculation is not in the same file as the CSAS6 calculation, the MAVRIC input would start by moving the KENO-VI results into the SCALE temporary area:}(hj h j hhh!NhNubah}(h]h]h]h]h]uhhKh!h"hMh j hhubj)}(h=shell
copy %RTNDIR%\caas.kenovi.fissionSource.msm fissionSource.msm
copy %RTNDIR%\caas.kenovi.kenoNuBar.txt kenoNuBar.txt
endh]h/=shell
copy %RTNDIR%\caas.kenovi.fissionSource.msm fissionSource.msm
copy %RTNDIR%\caas.kenovi.kenoNuBar.txt kenoNuBar.txt
end}(hhh j ubah}(h]h]h]h]h]forcehighlight_args}jjjscaleuhjh!h"hMh j hhubhL)}(hThe materials and geometry blocks of the two MAVRIC input files for each of the two calculations,
``smplprbs/caas.mavricA.inp`` and ``smplprbs/caas.mavricB.inp``, look like the following:h](h/bThe materials and geometry blocks of the two MAVRIC input files for each of the two calculations,
}(hbThe materials and geometry blocks of the two MAVRIC input files for each of the two calculations,
h j hhh!NhNubj)}(h``smplprbs/caas.mavricA.inp``h]h/smplprbs/caas.mavricA.inp}(hhh j ubah}(h]h]h]h]h]uhjh j ubh/ and }(h and h j hhh!NhNubj)}(h``smplprbs/caas.mavricB.inp``h]h/smplprbs/caas.mavricB.inp}(hhh j ubah}(h]h]h]h]h]uhjh j ubh/, look like the following:}(h, look like the following:h j hhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hMh j hhubj)}(hX'-------------------------------------------------------------------------------
' Composition Block
'-------------------------------------------------------------------------------
read composition
pu-239 1 0 0.037047 end
pu-240 1 0 0.0017512 end
pu-241 1 0 0.00011674 end
cu 1 0 0.0013752 end
orconcrete 2 1.0 293.0 end
orconcrete 3 0.33198 293.0 end
ss304 4 1.0 293.0 end
end composition
'-------------------------------------------------------------------------------
' Geometry Block - SCALE standard geometry package (SGGP)
'-------------------------------------------------------------------------------
read geometry
global unit 1
com="jezebel"
sphere 1 6.38493 origin x=280 y=300 z=100
com="exterior of the building, roof, floor"
cuboid 10 1200 0 600 0 300.3175 -60.0
cuboid 11 1200 0 600 0 300.3175 300.0
cuboid 12 1200 0 600 0 0.0 -60
com="air space in building - two rooms and maze"
cuboid 20 1160 40 560 40 300 0
com="interior walls to form maze to prevent streaming"
cuboid 21 560 520 560 120 300 0
cuboid 22 680 640 480 40 300 0
com="exterior door"
cuboid 30 1200 1160 360 240 210 0
cuboid 31 1200 1199.6825 360 240 210 0
com="detector sphere"
sphere 40 15.0 origin x=1145 y=55 z=285
com="jezebel"
media 1 1 1 vol=1090.3277
com="walls, roof, floor"
media 3 1 10 -20 -11 -12 -30
media 4 1 11
media 2 1 12
com="air space (void) and maze walls"
media 0 1 20 -21 -22 -40 -11 -12 -1
media 3 1 21 -11 -12
media 3 1 22 -11 -12
com="exterior door"
media 0 1 30 -31
media 4 1 31
com="detector"
media 0 1 40 vol=14137.167
boundary 10
end geometryh]h/X'-------------------------------------------------------------------------------
' Composition Block
'-------------------------------------------------------------------------------
read composition
pu-239 1 0 0.037047 end
pu-240 1 0 0.0017512 end
pu-241 1 0 0.00011674 end
cu 1 0 0.0013752 end
orconcrete 2 1.0 293.0 end
orconcrete 3 0.33198 293.0 end
ss304 4 1.0 293.0 end
end composition
'-------------------------------------------------------------------------------
' Geometry Block - SCALE standard geometry package (SGGP)
'-------------------------------------------------------------------------------
read geometry
global unit 1
com="jezebel"
sphere 1 6.38493 origin x=280 y=300 z=100
com="exterior of the building, roof, floor"
cuboid 10 1200 0 600 0 300.3175 -60.0
cuboid 11 1200 0 600 0 300.3175 300.0
cuboid 12 1200 0 600 0 0.0 -60
com="air space in building - two rooms and maze"
cuboid 20 1160 40 560 40 300 0
com="interior walls to form maze to prevent streaming"
cuboid 21 560 520 560 120 300 0
cuboid 22 680 640 480 40 300 0
com="exterior door"
cuboid 30 1200 1160 360 240 210 0
cuboid 31 1200 1199.6825 360 240 210 0
com="detector sphere"
sphere 40 15.0 origin x=1145 y=55 z=285
com="jezebel"
media 1 1 1 vol=1090.3277
com="walls, roof, floor"
media 3 1 10 -20 -11 -12 -30
media 4 1 11
media 2 1 12
com="air space (void) and maze walls"
media 0 1 20 -21 -22 -40 -11 -12 -1
media 3 1 21 -11 -12
media 3 1 22 -11 -12
com="exterior door"
media 0 1 30 -31
media 4 1 31
com="detector"
media 0 1 40 vol=14137.167
boundary 10
end geometry}(hhh j ubah}(h]h]h]h]h]forcehighlight_args}jjjscaleuhjh!h"hMh j hhubhL)}(hThe response functions used to compute the doses will be the standard
flux-to-dose rate conversion factors for neutrons and photons. These are
defined in the definitions block. Note that these responses have units
of (rem/hr)/(/cm :sup:`2`/s).h](h/The response functions used to compute the doses will be the standard
flux-to-dose rate conversion factors for neutrons and photons. These are
defined in the definitions block. Note that these responses have units
of (rem/hr)/(/cm }(hThe response functions used to compute the doses will be the standard
flux-to-dose rate conversion factors for neutrons and photons. These are
defined in the definitions block. Note that these responses have units
of (rem/hr)/(/cm h j hhh!NhNubj\)}(h:sup:`2`h]h/2}(hhh j ubah}(h]h]h]h]h]uhj[h j ubh//s).}(h/s).h j hhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hMh j hhubj)}(hX'-------------------------------------------------------------------------------
' Definitions Block
'-------------------------------------------------------------------------------
read definitions
response 5
title="ANSI (1977) neutron flux-to-dose-rate"
specialDose=9029
end response
response 6
title="ANSI (1977) photon flux-to-dose-rate"
specialDose=9504
end response
end definitionsh]h/X'-------------------------------------------------------------------------------
' Definitions Block
'-------------------------------------------------------------------------------
read definitions
response 5
title="ANSI (1977) neutron flux-to-dose-rate"
specialDose=9029
end response
response 6
title="ANSI (1977) photon flux-to-dose-rate"
specialDose=9504
end response
end definitions}(hhh j
ubah}(h]h]h]h]h]forcehighlight_args}jjjscaleuhjh!h"hMh j hhubhL)}(hXThe source used by each MAVRIC simulation will be based on the fission
distribution mesh source determined by KENO-VI. The strength of the
source can be specified by the total number of fissions that occurred in
the criticality event. Fission photons will be added for :sup:`239`\ Pu.
MAVRIC will determine the total source strength, including the fission
photons, from the value of saved by KENO-VI and the multiplicity data
from the fission photon data file.h](h/X
The source used by each MAVRIC simulation will be based on the fission
distribution mesh source determined by KENO-VI. The strength of the
source can be specified by the total number of fissions that occurred in
the criticality event. Fission photons will be added for }(hX
The source used by each MAVRIC simulation will be based on the fission
distribution mesh source determined by KENO-VI. The strength of the
source can be specified by the total number of fissions that occurred in
the criticality event. Fission photons will be added for h j*
hhh!NhNubj\)}(h
:sup:`239`h]h/239}(hhh j3
ubah}(h]h]h]h]h]uhj[h j*
ubh/ Pu.
MAVRIC will determine the total source strength, including the fission
photons, from the value of saved by KENO-VI and the multiplicity data
from the fission photon data file.}(h\ Pu.
MAVRIC will determine the total source strength, including the fission
photons, from the value of saved by KENO-VI and the multiplicity data
from the fission photon data file.h j*
hhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hMh j hhubj)}(hXZ'-------------------------------------------------------------------------------
' Sources Block
'-------------------------------------------------------------------------------
read sources
src 1
meshSourceFile=”fissionSource.msm"
origin x=280 y=300 z=100
fissions=1.0e18
mixture=1
end src
end sourcesh]h/XZ'-------------------------------------------------------------------------------
' Sources Block
'-------------------------------------------------------------------------------
read sources
src 1
meshSourceFile=”fissionSource.msm"
origin x=280 y=300 z=100
fissions=1.0e18
mixture=1
end src
end sources}(hhh jL
ubah}(h]h]h]h]h]forcehighlight_args}jjjscaleuhjh!h"hMh j hhubhL)}(hNote that further multiplication needs to be turned off in MAVRIC using
the “fissionMult=0” keyword in the parameter block as shown below.h]h/Note that further multiplication needs to be turned off in MAVRIC using
the “fissionMult=0” keyword in the parameter block as shown below.}(hj`
h j^
hhh!NhNubah}(h]h]h]h]h]uhhKh!h"hM+h j hhubhL)}(hXJFor the responses from the tallies, MAVRIC usually calculates dose rates
(rem/hr) using a source strength in particles/s. For this example
problem, instead of a source rate, we used a total number of particles
(by specifying the number of fissions). Hence, the computed fluxes will
have units of particles/cm\ :sup:`2` and the computed responses using
the standard dose responses from the cross section libraries will have
units of rem s/hr. To get a dose in rem, the responses need to be
multiplied by (3600 s/hr)\ :sup:`-1`. This can be done using the MAVRIC
tally multiplier keyword.h](h/X6For the responses from the tallies, MAVRIC usually calculates dose rates
(rem/hr) using a source strength in particles/s. For this example
problem, instead of a source rate, we used a total number of particles
(by specifying the number of fissions). Hence, the computed fluxes will
have units of particles/cm }(hX6For the responses from the tallies, MAVRIC usually calculates dose rates
(rem/hr) using a source strength in particles/s. For this example
problem, instead of a source rate, we used a total number of particles
(by specifying the number of fissions). Hence, the computed fluxes will
have units of particles/cm\ h jl
hhh!NhNubj\)}(h:sup:`2`h]h/2}(hhh ju
ubah}(h]h]h]h]h]uhj[h jl
ubh/ and the computed responses using
the standard dose responses from the cross section libraries will have
units of rem s/hr. To get a dose in rem, the responses need to be
multiplied by (3600 s/hr) }(h and the computed responses using
the standard dose responses from the cross section libraries will have
units of rem s/hr. To get a dose in rem, the responses need to be
multiplied by (3600 s/hr)\ h jl
hhh!NhNubj\)}(h :sup:`-1`h]h/-1}(hhh j
ubah}(h]h]h]h]h]uhj[h jl
ubh/=. This can be done using the MAVRIC
tally multiplier keyword.}(h=. This can be done using the MAVRIC
tally multiplier keyword.h jl
hhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hM.h j hhubhL)}(hXEach MAVRIC simulation will need a discrete-ordinates mesh. The planes
in each dimension where there are material changes are listed in
:numref:`tab4a2`. In addition to these planes, the discrete-ordinates mesh
should also subdivide the thick shields in the direction of particle
travel. For example, the walls of the maze should be divided to better
model the radiation attenuation through the walls in the Denovo
calculation. The interior walls of the building will reflect particles,
so the first few centimeters are the most important to capture in the
importance map. Mesh planes should also be added that correspond to the
mesh source after it is placed into the geometry model.h](h/Each MAVRIC simulation will need a discrete-ordinates mesh. The planes
in each dimension where there are material changes are listed in
}(hEach MAVRIC simulation will need a discrete-ordinates mesh. The planes
in each dimension where there are material changes are listed in
h j
hhh!NhNubhW)}(h:numref:`tab4a2`h]j)}(hj
h]h/tab4a2}(hhh j
ubah}(h]h](jwstd
std-numrefeh]h]h]uhjh j
ubah}(h]h]h]h]h]refdocj refdomainj
reftypenumrefrefexplicitrefwarnjtab4a2uhhVh!h"hM8h j
ubh/X. In addition to these planes, the discrete-ordinates mesh
should also subdivide the thick shields in the direction of particle
travel. For example, the walls of the maze should be divided to better
model the radiation attenuation through the walls in the Denovo
calculation. The interior walls of the building will reflect particles,
so the first few centimeters are the most important to capture in the
importance map. Mesh planes should also be added that correspond to the
mesh source after it is placed into the geometry model.}(hX. In addition to these planes, the discrete-ordinates mesh
should also subdivide the thick shields in the direction of particle
travel. For example, the walls of the maze should be divided to better
model the radiation attenuation through the walls in the Denovo
calculation. The interior walls of the building will reflect particles,
so the first few centimeters are the most important to capture in the
importance map. Mesh planes should also be added that correspond to the
mesh source after it is placed into the geometry model.h j
hhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hM8h j hhubh)}(h.. _tab4a2:h]h}(h]h]h]h]h]htab4a2uhh
hMCh j hhh!h"ubj)}(hhh](h))}(hMain planes in the buildingh]h/Main planes in the building}(hj
h j
ubah}(h]h]h]h]h]uhh(h!h"hMDh j
ubj)}(hhh](j)}(hhh]h}(h]h]h]h]h]colwidthK uhjh j
ubj)}(hhh]h}(h]h]h]h]h]colwidthKuhjh j
ubj)}(hhh]h}(h]h]h]h]h]colwidthK uhjh j
ubj)}(hhh]j)}(hhh](j)}(hhh]hL)}(h**x**h]j)}(hjh]h/x}(hhh jubah}(h]h]h]h]h]uhjh jubah}(h]h]h]h]h]uhhKh!h"hMHh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h**y**h]j)}(hj=h]h/y}(hhh j?ubah}(h]h]h]h]h]uhjh j;ubah}(h]h]h]h]h]uhhKh!h"hMHh j8ubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h**z**h]j)}(hj]h]h/z}(hhh j_ubah}(h]h]h]h]h]uhjh j[ubah}(h]h]h]h]h]uhhKh!h"hMHh jXubah}(h]h]h]h]h]uhjh jubeh}(h]h]h]h]h]uhjh jubah}(h]h]h]h]h]uhjh j
ubj)}(hhh](j)}(hhh](j)}(hhh]hL)}(h0h]h/0}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMJh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(hjh]h/0}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMJh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h-60h]h/-60}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMJh jubah}(h]h]h]h]h]uhjh jubeh}(h]h]h]h]h]uhjh jubj)}(hhh](j)}(hhh]hL)}(h40h]h/40}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMLh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h40h]h/40}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMLh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(hjh]h/0}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMLh jubah}(h]h]h]h]h]uhjh jubeh}(h]h]h]h]h]uhjh jubj)}(hhh](j)}(hhh]hL)}(h520h]h/520}(hj)h j'ubah}(h]h]h]h]h]uhhKh!h"hMNh j$ubah}(h]h]h]h]h]uhjh j!ubj)}(hhh]hL)}(h120h]h/120}(hj@h j>ubah}(h]h]h]h]h]uhhKh!h"hMNh j;ubah}(h]h]h]h]h]uhjh j!ubj)}(hhh]hL)}(h210h]h/210}(hjWh jUubah}(h]h]h]h]h]uhhKh!h"hMNh jRubah}(h]h]h]h]h]uhjh j!ubeh}(h]h]h]h]h]uhjh jubj)}(hhh](j)}(hhh]hL)}(h560h]h/560}(hjwh juubah}(h]h]h]h]h]uhhKh!h"hMPh jrubah}(h]h]h]h]h]uhjh joubj)}(hhh]hL)}(h240h]h/240}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMPh jubah}(h]h]h]h]h]uhjh joubj)}(hhh]hL)}(h300h]h/300}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMPh jubah}(h]h]h]h]h]uhjh joubeh}(h]h]h]h]h]uhjh jubj)}(hhh](j)}(hhh]hL)}(h640h]h/640}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMRh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h360h]h/360}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMRh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h300.318h]h/300.318}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMRh jubah}(h]h]h]h]h]uhjh jubeh}(h]h]h]h]h]uhjh jubj)}(hhh](j)}(hhh]hL)}(h680h]h/680}(hj
h j
ubah}(h]h]h]h]h]uhhKh!h"hMTh j
ubah}(h]h]h]h]h]uhjh j
ubj)}(hhh]hL)}(h480h]h/480}(hj*
h j(
ubah}(h]h]h]h]h]uhhKh!h"hMTh j%
ubah}(h]h]h]h]h]uhjh j
ubj)}(hhh]h}(h]h]h]h]h]uhjh j
ubeh}(h]h]h]h]h]uhjh jubj)}(hhh](j)}(hhh]hL)}(h1160h]h/1160}(hjS
h jQ
ubah}(h]h]h]h]h]uhhKh!h"hMVh jN
ubah}(h]h]h]h]h]uhjh jK
ubj)}(hhh]hL)}(h560h]h/560}(hjj
h jh
ubah}(h]h]h]h]h]uhhKh!h"hMVh je
ubah}(h]h]h]h]h]uhjh jK
ubj)}(hhh]h}(h]h]h]h]h]uhjh jK
ubeh}(h]h]h]h]h]uhjh jubj)}(hhh](j)}(hhh]hL)}(h1199.68h]h/1199.68}(hj
h j
ubah}(h]h]h]h]h]uhhKh!h"hMXh j
ubah}(h]h]h]h]h]uhjh j
ubj)}(hhh]hL)}(h600h]h/600}(hj
h j
ubah}(h]h]h]h]h]uhhKh!h"hMXh j
ubah}(h]h]h]h]h]uhjh j
ubj)}(hhh]h}(h]h]h]h]h]uhjh j
ubeh}(h]h]h]h]h]uhjh jubj)}(hhh](j)}(hhh]hL)}(h1200h]h/1200}(hj
h j
ubah}(h]h]h]h]h]uhhKh!h"hMZh j
ubah}(h]h]h]h]h]uhjh j
ubj)}(hhh]h}(h]h]h]h]h]uhjh j
ubj)}(hhh]h}(h]h]h]h]h]uhjh j
ubeh}(h]h]h]h]h]uhjh jubeh}(h]h]h]h]h]uhjh j
ubeh}(h]h]h]h]h]colsKuhjh j
ubeh}(h](id14j
eh]h]tab4a2ah]h]j|centeruhjh j hhh!h"hNj~}jj
sj}j
j
subh$)}(hhh](h))}(hDetector doses using CADISh]h/Detector doses using CADIS}(hjh jhhh!NhNubah}(h]h]h]h]h]uhh(h jhhh!h"hM_ubhL)}(hXLThe grid geometry for this calculation should also include planes that bound the adjoint source,
which is the detector area (these values are shown in brackets [] below). The definitions block in
``smplprbs/caas.mavricA.inp`` also includes the location of the center of the detector, which is used in the adjoint source description.h](h/The grid geometry for this calculation should also include planes that bound the adjoint source,
which is the detector area (these values are shown in brackets [] below). The definitions block in
}(hThe grid geometry for this calculation should also include planes that bound the adjoint source,
which is the detector area (these values are shown in brackets [] below). The definitions block in
h j&hhh!NhNubj)}(h``smplprbs/caas.mavricA.inp``h]h/smplprbs/caas.mavricA.inp}(hhh j/ubah}(h]h]h]h]h]uhjh j&ubh/k also includes the location of the center of the detector, which is used in the adjoint source description.}(hk also includes the location of the center of the detector, which is used in the adjoint source description.h j&hhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hMah jhhubj)}(hX.location 1
position 1145 55 285
end location
gridGeometry 1
title="mesh for discrete ordinates 57 x 47 x 31 = 83049"
xplanes 0 10 20 30 35
40 120 160 240
270 272 274 276 278 280 282 284 286 288 290
360 440
520 525 530 550 555
560 600
640 645 650 670 675
680 760 840 920 1000 1080
[1130 1140 1150]
1160 1165 1170 1180 1190
1199.6825
1200 end
yplanes 0 10 20 30 35
40 [50 60] 70
120 125 130 140 200
240 280 290 292 294 296 298 300 302 304 306 308 310
320
360 440 460 470 475
480
560 565 570 580 590
600 end
zplanes -60 -30 -20 -10 -5
0 45
90 92 94 96 98 100 102 104 106 108 110
140 175
210 255 [280 290]
300
300.3175 end
end gridGeometryh]h/X.location 1
position 1145 55 285
end location
gridGeometry 1
title="mesh for discrete ordinates 57 x 47 x 31 = 83049"
xplanes 0 10 20 30 35
40 120 160 240
270 272 274 276 278 280 282 284 286 288 290
360 440
520 525 530 550 555
560 600
640 645 650 670 675
680 760 840 920 1000 1080
[1130 1140 1150]
1160 1165 1170 1180 1190
1199.6825
1200 end
yplanes 0 10 20 30 35
40 [50 60] 70
120 125 130 140 200
240 280 290 292 294 296 298 300 302 304 306 308 310
320
360 440 460 470 475
480
560 565 570 580 590
600 end
zplanes -60 -30 -20 -10 -5
0 45
90 92 94 96 98 100 102 104 106 108 110
140 175
210 255 [280 290]
300
300.3175 end
end gridGeometry}(hhh jHubah}(h]h]h]h]h]forcehighlight_args}jjjscaleuhjh!h"hMeh jhhubhL)}(hXThe tallies are region tallies over the detector region (the 10th media
card in unit 1) using the appropriate response function for the particle
type of the tally. The volume of the detector sphere needs to be listed
in the geometry block so that the fluxes and tallies will be correctly
computed. The importance map uses standard CADIS to bias the particles
towards the detectors, optimizing the calculation of the total dose (by
listing both response functions together, the total response will be
used in the adjoint source).h]h/XThe tallies are region tallies over the detector region (the 10th media
card in unit 1) using the appropriate response function for the particle
type of the tally. The volume of the detector sphere needs to be listed
in the geometry block so that the fluxes and tallies will be correctly
computed. The importance map uses standard CADIS to bias the particles
towards the detectors, optimizing the calculation of the total dose (by
listing both response functions together, the total response will be
used in the adjoint source).}(hj\h jZhhh!NhNubah}(h]h]h]h]h]uhhKh!h"hMh jhhubj)}(hX'-------------------------------------------------------------------------------
' Tallies Block
'-------------------------------------------------------------------------------
read tallies
regionTally 1
title="Doses seen by the detector"
neutron
unit=1 region=10
responseID=5
multiplier=2.777778e-4
end regionTally
regionTally 2
title="Doses seen by the detector"
photon
unit=1 region=10
responseID=6
multiplier=2.777778e-4
end regionTally
end tallies
'-------------------------------------------------------------------------------
' Parameters Block - 3 min batch
'-------------------------------------------------------------------------------
read parameters
randomSeed=3263827
perBatch=654000 batches=40
fissionMult=0
end parameters
'-------------------------------------------------------------------------------
' Importance Map Block
'-------------------------------------------------------------------------------
read importanceMap
gridGeometryID=1
adjointSource 1
locationID=1
responseIDs 5 6 end
end adjointSource
end importanceMaph]h/X'-------------------------------------------------------------------------------
' Tallies Block
'-------------------------------------------------------------------------------
read tallies
regionTally 1
title="Doses seen by the detector"
neutron
unit=1 region=10
responseID=5
multiplier=2.777778e-4
end regionTally
regionTally 2
title="Doses seen by the detector"
photon
unit=1 region=10
responseID=6
multiplier=2.777778e-4
end regionTally
end tallies
'-------------------------------------------------------------------------------
' Parameters Block - 3 min batch
'-------------------------------------------------------------------------------
read parameters
randomSeed=3263827
perBatch=654000 batches=40
fissionMult=0
end parameters
'-------------------------------------------------------------------------------
' Importance Map Block
'-------------------------------------------------------------------------------
read importanceMap
gridGeometryID=1
adjointSource 1
locationID=1
responseIDs 5 6 end
end adjointSource
end importanceMap}(hhh jhubah}(h]h]h]h]h]forcehighlight_args}jjjscaleuhjh!h"hMh jhhubhL)}(hXPThe results of this example problem are shown in :numref:`tab4a3`.
Calculation times were 12 minutes for Denovo and 135 minutes for Monaco.
Note that the uncertainty for the photon dose is much higher than the
neutron dose uncertainty. This is because the simulation was optimized
for the calculation of total dose, and the photon component of the total
dose is less than 2%. Had a separate calculation been done that used an
adjoint source of just the photon response, the photon dose rate
uncertainty would have been much smaller but at the expense of the
neutron dose rate uncertainty. A single calculation could have also been
performed using two adjoint sources, one using the neutron dose response
and one using the photon dose response, and forward weighting to help
calculate each component of dose with more uniform relative
uncertainties.h](h/1The results of this example problem are shown in }(h1The results of this example problem are shown in h jzhhh!NhNubhW)}(h:numref:`tab4a3`h]j)}(hjh]h/tab4a3}(hhh jubah}(h]h](jwstd
std-numrefeh]h]h]uhjh jubah}(h]h]h]h]h]refdocj refdomainjreftypenumrefrefexplicitrefwarnjtab4a3uhhVh!h"hMh jzubh/X.
Calculation times were 12 minutes for Denovo and 135 minutes for Monaco.
Note that the uncertainty for the photon dose is much higher than the
neutron dose uncertainty. This is because the simulation was optimized
for the calculation of total dose, and the photon component of the total
dose is less than 2%. Had a separate calculation been done that used an
adjoint source of just the photon response, the photon dose rate
uncertainty would have been much smaller but at the expense of the
neutron dose rate uncertainty. A single calculation could have also been
performed using two adjoint sources, one using the neutron dose response
and one using the photon dose response, and forward weighting to help
calculate each component of dose with more uniform relative
uncertainties.}(hX.
Calculation times were 12 minutes for Denovo and 135 minutes for Monaco.
Note that the uncertainty for the photon dose is much higher than the
neutron dose uncertainty. This is because the simulation was optimized
for the calculation of total dose, and the photon component of the total
dose is less than 2%. Had a separate calculation been done that used an
adjoint source of just the photon response, the photon dose rate
uncertainty would have been much smaller but at the expense of the
neutron dose rate uncertainty. A single calculation could have also been
performed using two adjoint sources, one using the neutron dose response
and one using the photon dose response, and forward weighting to help
calculate each component of dose with more uniform relative
uncertainties.h jzhhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hMh jhhubh)}(h.. _tab4a3:h]h}(h]h]h]h]h]htab4a3uhh
hMh jhhh!h"ubj)}(hhh](h))}(hDoses seen by the detectorh]h/Doses seen by the detector}(hjh jubah}(h]h]h]h]h]uhh(h!h"hMh jubj)}(hhh](j)}(hhh]h}(h]h]h]h]h]colwidthK
uhjh jubj)}(hhh]h}(h]h]h]h]h]colwidthKuhjh jubj)}(hhh]h}(h]h]h]h]h]colwidthK
uhjh jubj)}(hhh]j)}(hhh](j)}(hhh]hL)}(h** **h]h/** **}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h **Value**h]j)}(hj
h]h/Value}(hhh jubah}(h]h]h]h]h]uhjh jubah}(h]h]h]h]h]uhhKh!h"hMh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h**Rel.**h]j)}(hj-h]h/Rel.}(hhh j/ubah}(h]h]h]h]h]uhjh j+ubah}(h]h]h]h]h]uhhKh!h"hMh j(ubah}(h]h]h]h]h]uhjh jubeh}(h]h]h]h]h]uhjh jubah}(h]h]h]h]h]uhjh jubj)}(hhh](j)}(hhh](j)}(hhh]hL)}(h**Dose**h]j)}(hj_h]h/Dose}(hhh jaubah}(h]h]h]h]h]uhjh j]ubah}(h]h]h]h]h]uhhKh!h"hMh jZubah}(h]h]h]h]h]uhjh jWubj)}(hhh]hL)}(h **(rem)**h]j)}(hjh]h/(rem)}(hhh jubah}(h]h]h]h]h]uhjh j}ubah}(h]h]h]h]h]uhhKh!h"hMh jzubah}(h]h]h]h]h]uhjh jWubj)}(hhh]hL)}(h**Unc.**h]j)}(hjh]h/Unc.}(hhh jubah}(h]h]h]h]h]uhjh jubah}(h]h]h]h]h]uhhKh!h"hMh jubah}(h]h]h]h]h]uhjh jWubeh}(h]h]h]h]h]uhjh jTubj)}(hhh](j)}(hhh]hL)}(hneutronh]h/neutron}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h1539h]h/1539}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h0.78%h]h/0.78%}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMh jubah}(h]h]h]h]h]uhjh jubeh}(h]h]h]h]h]uhjh jTubj)}(hhh](j)}(hhh]hL)}(hphotonh]h/photon}(hjh jubah}(h]h]h]h]h]uhhKh!h"hMh jubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h30.0h]h/30.0}(hj-h j+ubah}(h]h]h]h]h]uhhKh!h"hMh j(ubah}(h]h]h]h]h]uhjh jubj)}(hhh]hL)}(h8.00%h]h/8.00%}(hjDh jBubah}(h]h]h]h]h]uhhKh!h"hMh j?ubah}(h]h]h]h]h]uhjh jubeh}(h]h]h]h]h]uhjh jTubeh}(h]h]h]h]h]uhjh jubeh}(h]h]h]h]h]colsKuhjh jubeh}(h](id15jeh]h]tab4a3ah]h]j|centeruhjh jhhh!h"hNj~}jnjsj}jjsubeh}(h]detector-doses-using-cadisah]h]detector doses using cadisah]h]uhh#h j hhh!h"hM_ubh$)}(hhh](h))}(hDose map using FW-CADISh]h/Dose map using FW-CADIS}(hjh jhhh!NhNubah}(h]h]h]h]h]uhh(h j|hhh!h"hMubhL)}(hThe grid geometry for this calculation does not need extra planes around the detector. The grid geometry in
``smplprbs/caas.mavricB.inp`` looks like the following:h](h/lThe grid geometry for this calculation does not need extra planes around the detector. The grid geometry in
}(hlThe grid geometry for this calculation does not need extra planes around the detector. The grid geometry in
h jhhh!NhNubj)}(h``smplprbs/caas.mavricB.inp``h]h/smplprbs/caas.mavricB.inp}(hhh jubah}(h]h]h]h]h]uhjh jubh/ looks like the following:}(h looks like the following:h jhhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hMh j|hhubj)}(hXgridGeometry 1
title="mesh for discrete ordinates 46x36x23 = 38088"
xplanes 0 10 20 30 35
40 120 160 240
270 272 274 276 278 280 282 284 286 288 290
360 440
520 525 530 550 555
560 600
640 645 650 670 675
680 760 840 920 1000 1080
1160 1165 1170 1180 1190
1199.6825
1200 end
yplanes 0 10 20 30 35
40
120 125 130 140 200
240 280 290 292 294 296 298 300 302 304 306 308 310
320
360 440 460 470 475
480
560 565 570 580 590
600 end
zplanes -60 -30 -20 -10 -5
0 45
90 92 94 96 98 100 102 104 106 108 110
140 175
210 255
300 300.3175 end
end gridGeometryh]h/XgridGeometry 1
title="mesh for discrete ordinates 46x36x23 = 38088"
xplanes 0 10 20 30 35
40 120 160 240
270 272 274 276 278 280 282 284 286 288 290
360 440
520 525 530 550 555
560 600
640 645 650 670 675
680 760 840 920 1000 1080
1160 1165 1170 1180 1190
1199.6825
1200 end
yplanes 0 10 20 30 35
40
120 125 130 140 200
240 280 290 292 294 296 298 300 302 304 306 308 310
320
360 440 460 470 475
480
560 565 570 580 590
600 end
zplanes -60 -30 -20 -10 -5
0 45
90 92 94 96 98 100 102 104 106 108 110
140 175
210 255
300 300.3175 end
end gridGeometry}(hhh jubah}(h]h]h]h]h]forcehighlight_args}jjjscaleuhjh!h"hMh j|hhubhL)}(haA second grid geometry also needs to be added to the definitions block for the mesh tally to use.h]h/aA second grid geometry also needs to be added to the definitions block for the mesh tally to use.}(hjh jhhh!NhNubah}(h]h]h]h]h]uhhKh!h"hMh j|hhubj)}(hgridGeometry 2
title="mesh for uniform mesh tally - 40x40x30 cm voxels"
xLinear 30 0.0 1200.0
yLinear 15 0.0 600.0
zLinear 10 0.0 300.0
end gridGeometryh]h/gridGeometry 2
title="mesh for uniform mesh tally - 40x40x30 cm voxels"
xLinear 30 0.0 1200.0
yLinear 15 0.0 600.0
zLinear 10 0.0 300.0
end gridGeometry}(hhh jubah}(h]h]h]h]h]forcehighlight_args}jjjscaleuhjh!h"hMh j|hhubhL)}(hXThe mesh tallies for each particle type are listed, along with the appropriate response function. The importance map uses FW-CADIS to better spread the particles out over the entire geometry, optimized for the calculation of total dose in the void regions.h]h/XThe mesh tallies for each particle type are listed, along with the appropriate response function. The importance map uses FW-CADIS to better spread the particles out over the entire geometry, optimized for the calculation of total dose in the void regions.}(hjh jhhh!NhNubah}(h]h]h]h]h]uhhKh!h"hMh j|hhubj)}(hX*'-------------------------------------------------------------------------------
' Tallies Block
'-------------------------------------------------------------------------------
read tallies
meshTally 1
title="Neutron doses mapped over the entire building"
neutron
gridGeometryID=2
responseID=5
noGroupFluxes
multiplier=2.777778e-4
end meshTally
meshTally 2
title="Photon doses mapped over the entire building"
photon
gridGeometryID=2
responseID=6
noGroupFluxes
multiplier=2.777778e-4
end meshTally
end tallies
'-------------------------------------------------------------------------------
' Parameters Block - 3 min batch
'-------------------------------------------------------------------------------
read parameters
randomSeed=3263827
perBatch=669000 batches=40
fissionMult=0
end parameters
'-------------------------------------------------------------------------------
' Importance Map Block
'-------------------------------------------------------------------------------
read importanceMap
gridGeometryID=1
adjointSource 1
boundingBox 1200 0 600 0 300.3175 -60.0
responseIDs 5 6 end
mixture=0
end adjointSource
respWeighting
end importanceMaph]h/X*'-------------------------------------------------------------------------------
' Tallies Block
'-------------------------------------------------------------------------------
read tallies
meshTally 1
title="Neutron doses mapped over the entire building"
neutron
gridGeometryID=2
responseID=5
noGroupFluxes
multiplier=2.777778e-4
end meshTally
meshTally 2
title="Photon doses mapped over the entire building"
photon
gridGeometryID=2
responseID=6
noGroupFluxes
multiplier=2.777778e-4
end meshTally
end tallies
'-------------------------------------------------------------------------------
' Parameters Block - 3 min batch
'-------------------------------------------------------------------------------
read parameters
randomSeed=3263827
perBatch=669000 batches=40
fissionMult=0
end parameters
'-------------------------------------------------------------------------------
' Importance Map Block
'-------------------------------------------------------------------------------
read importanceMap
gridGeometryID=1
adjointSource 1
boundingBox 1200 0 600 0 300.3175 -60.0
responseIDs 5 6 end
mixture=0
end adjointSource
respWeighting
end importanceMap}(hhh jubah}(h]h]h]h]h]forcehighlight_args}jjjscaleuhjh!h"hMh j|hhubhL)}(hXThis calculation used 5 minutes for the forward Denovo S\ :sub:`N`
calculation, 6 minutes for the adjoint Denovo, and 120 minutes for the
Monaco Monte Carlo. The resulting mesh tally is shown in
:numref:`fig4a6`. The uncertainties for each voxel are shown in
:numref:`fig4a7`. Note that the dose in the voxel containing the
detector (not shown in the figures) is 1.552×10\ :sup:`3` rem with a
relative uncertainty of 3.8%, closely matching the value calculated
with the first MAVRIC simulation.h](h/:This calculation used 5 minutes for the forward Denovo S }(h:This calculation used 5 minutes for the forward Denovo S\ h jhhh!NhNubj})}(h:sub:`N`h]h/N}(hhh j
ubah}(h]h]h]h]h]uhj|h jubh/
calculation, 6 minutes for the adjoint Denovo, and 120 minutes for the
Monaco Monte Carlo. The resulting mesh tally is shown in
}(h
calculation, 6 minutes for the adjoint Denovo, and 120 minutes for the
Monaco Monte Carlo. The resulting mesh tally is shown in
h jhhh!NhNubhW)}(h:numref:`fig4a6`h]j)}(hjh]h/fig4a6}(hhh j!ubah}(h]h](jwstd
std-numrefeh]h]h]uhjh jubah}(h]h]h]h]h]refdocj refdomainj+reftypenumrefrefexplicitrefwarnjfig4a6uhhVh!h"hM=h jubh/0. The uncertainties for each voxel are shown in
}(h0. The uncertainties for each voxel are shown in
h jhhh!NhNubhW)}(h:numref:`fig4a7`h]j)}(hjDh]h/fig4a7}(hhh jFubah}(h]h](jwstd
std-numrefeh]h]h]uhjh jBubah}(h]h]h]h]h]refdocj refdomainjPreftypenumrefrefexplicitrefwarnjfig4a7uhhVh!h"hM=h jubh/c. Note that the dose in the voxel containing the
detector (not shown in the figures) is 1.552×10 }(hc. Note that the dose in the voxel containing the
detector (not shown in the figures) is 1.552×10\ h jhhh!NhNubj\)}(h:sup:`3`h]h/3}(hhh jgubah}(h]h]h]h]h]uhj[h jubh/q rem with a
relative uncertainty of 3.8%, closely matching the value calculated
with the first MAVRIC simulation.}(hq rem with a
relative uncertainty of 3.8%, closely matching the value calculated
with the first MAVRIC simulation.h jhhh!NhNubeh}(h]h]h]h]h]uhhKh!h"hM=h j|hhubhL)}(hX:Similar to the detector doses above, a single calculation could have
also been performed for the dose maps using two adjoint sources, one
using the neutron dose response and one using the photon dose response,
and forward weighting to help calculate each component of dose with more
uniform relative uncertainties.h]h/X:Similar to the detector doses above, a single calculation could have
also been performed for the dose maps using two adjoint sources, one
using the neutron dose response and one using the photon dose response,
and forward weighting to help calculate each component of dose with more
uniform relative uncertainties.}(hjh jhhh!NhNubah}(h]h]h]h]h]uhhKh!h"hMFh j|hhubh)}(h.. _fig4a6:h]h}(h]h]h]h]h]hfig4a6uhh
hMLh j|hhh!h"ubjM)}(hhh](jR)}(hw.. figure:: figs/MAVRIC/4a6.png
:align: center
Dose (rem) results for the z=100 cm plane (containing the source).
h]h}(h]h]h]h]h]urifigs/MAVRIC/4a6.pngja}jcjsuhjQh jh!h"hMPubje)}(hBDose (rem) results for the z=100 cm plane (containing the source).h]h/BDose (rem) results for the z=100 cm plane (containing the source).}(hjh jubah}(h]h]h]h]h]uhjdh!h"hMPh jubeh}(h](id16jeh]h]fig4a6ah]h]j|centeruhjLhMPh j|hhh!h"j~}jjsj}jjsubh)}(h.. _fig4a7:h]h}(h]h]h]h]h]hfig4a7uhh
hMRh j|hhh!h"ubjM)}(hhh](jR)}(h.. figure:: figs/MAVRIC/fig4a7.png
:align: center
Relative uncertainties in the dose, most less than 1%, for the z=100 cm plane.
h]h}(h]h]h]h]h]urifigs/MAVRIC/fig4a7.pngja}jcjsuhjQh jh!h"hMVubje)}(hNRelative uncertainties in the dose, most less than 1%, for the z=100 cm plane.h]h/NRelative uncertainties in the dose, most less than 1%, for the z=100 cm plane.}(hjh jubah}(h]h]h]h]h]uhjdh!h"hMVh jubeh}(h](id17jeh]h]fig4a7ah]h]j|centeruhjLhMVh j|hhh!h"j~}jjsj}jjsubeh}(h]dose-map-using-fw-cadisah]h]dose map using fw-cadisah]h]uhh#h j hhh!h"hMubeh}(h]mavric-transport-calculationsah]h]mavric transport calculationsah]h]uhh#h jhhh!h"hMubeh}(h]example-problemah]h]example problemah]h]uhh#h h%hhh!h"hMubh$)}(hhh](h))}(hSummaryh]h/Summary}(hjh jhhh!NhNubah}(h]h]h]h]h]uhh(h jhhh!h"hMYubhL)}(hXSCALE now has the capability to do detailed simulations of criticality accident alarm systems. The advanced variance reduction capabilities of the MAVRIC radiation transport sequence allow for the full three-dimensional analysis of CAAS problems in reasonable amounts of computer time. This enables the use of more realistic source definitions, such as a detailed spatial/energy dependent fission source distribution determined by the KENO-VI criticality code, and the critical assembly itself can be included in the transport model.h]h/XSCALE now has the capability to do detailed simulations of criticality accident alarm systems. The advanced variance reduction capabilities of the MAVRIC radiation transport sequence allow for the full three-dimensional analysis of CAAS problems in reasonable amounts of computer time. This enables the use of more realistic source definitions, such as a detailed spatial/energy dependent fission source distribution determined by the KENO-VI criticality code, and the critical assembly itself can be included in the transport model.}(hj!h jhhh!NhNubah}(h]h]h]h]h]uhhKh!h"hM[h jhhubhL)}(hhh](h ht)}(hhh](h label)}(hhh]h/miller_corrected_2013}(hhh j6ubah}(h]h]h]h]h]support_smartquotesuhj4h j1ubhL)}(hhh](h/T.M. Miller and D.E. Peplow.}(hT.M. Miller and D.E. Peplow.h jDubh/ }(hj<h hL)}(hhh](h/L.L. Wetzel, B.}(hL.L. Wetzel, B.h jPubh/ }(h h jPubh/O’Donnell, and W.D. Newmyer.}(hO'Donnell, and W.D. Newmyer.h jPubh/ }(hj<h jPubh/CAAS}(hCAASh jPubh/ detector placement using }(h detector placement using h jPubh/MAVRIC}(hMAVRICh jPubh/.}(hjSh jPubjNj)}(hhh]h/,Transactions of the American Nuclear Society}(h,Transactions of the American Nuclear Societyh jyubah}(h]h]h]h]h]uhjh jPubh/, 107:644–647, January 2012.}(h, 107:644–647, January 2012.h jPubeh}(h]h]h]h]h]uhhKh j0)}(hhh](j5)}(hhh]h/wetzel_caas_2012}(hhh jubah}(h]h]h]h]h]jCuhj4h jubjPeh}(h]wetzel-caas-2012ah]hoah]wetzel_caas_2012ah]h]hadocnamejuhhth j-resolvedKubububh/eCorrected user guidance to perform three-dimensional criticality accident alarm system modeling with }(heCorrected user guidance to perform three-dimensional criticality accident alarm system modeling with h jDubh/SCALE}(hSCALEh jDubh/.}(hjSh jDubjNj)}(hhh]h/,Transactions of the American Nuclear Society}(h,Transactions of the American Nuclear Societyh jubah}(h]h]h]h]h]uhjh jDubh/, 108:498–501, January 2013.}(h, 108:498–501, January 2013.h jDubeh}(h]h]h]h]h]uhhKh j1ubeh}(h]miller-corrected-2013ah]hoah]miller_corrected_2013ah]h]jEajjuhhth j-jKubj0)}(hhh](j5)}(hhh]h/miller_guidance_2013}(hhh jubah}(h]h]h]h]h]jCuhj4h jubhL)}(hhh](h/Thomas}(hThomash jubh/ }(hj\h jPubh/Martin Miller and Douglas}(hMartin Miller and Douglash jubjh/
E. Peplow.}(h
E. Peplow.h jubh/ }(h h jubh/ Guidance }(h Guidance h jubh/ Detailing}(h Detailingh jubh/ }(h h jubh/Methods}(hMethodsh jubh/ to }(h to h jubh/ Calculate}(h Calculateh jubh/ }(h h jubh/Criticality}(hCriticalityh jubh/ }(h h jubh/Accident}(hAccidenth jubh/ }(h h jubh/Alarm}(hAlarmh jubh/ }(h h jubh/Detector}(hDetectorh jubh/ }(hjh jubh/Response}(hResponseh jubh/ and }(h and h jubh/Coverage}(hCoverageh jubh/.}(hjSh jubjNh/^Technical Report, Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States), January 2013.}(h^Technical Report, Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States), January 2013.h jubjNh/URL: }(hURL: h jubh reference)}(hhh]h/#https://www.osti.gov/biblio/1095654}(h#https://www.osti.gov/biblio/1095654h joubah}(h]h]h]h]h]refurijvuhjmh jubh/.}(hjSh jubeh}(h]h]h]h]h]uhhKh jubeh}(h]miller-guidance-2013ah]hoah]miller_guidance_2013ah]h]jajjuhhth j-jKubj0)}(hhh](j5)}(hhh]h/miller_guide_2013}(hhh jubah}(h]h]h]h]h]jCuhj4h jubhL)}(hhh](h/Thomas}(hThomash jubjh/Martin Miller and Douglas}(hMartin Miller and Douglash jubjh/
E. Peplow.}(h
E. Peplow.h jubh/ }(h h jubh/ Guide to }(h Guide to h jubh/
Performing}(h
Performingh jubh/ }(h h jubh/
Computational}(h
Computationalh jubh/ }(h h jubh/Analysis}(hAnalysish jubh/ of }(h of h jubh/Criticality}(hCriticalityh jubh/ }(h h jubh/Accident}(hAccidenth jubh/ }(h h jubh/Alarm}(hAlarmh jubh/ }(hjh jubh/Systems}(hSystemsh jubh/.}(hjSh jubjNh/qTechnical Report ORNL/TM-2013/211, Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States), September 2013.}(hqTechnical Report ORNL/TM-2013/211, Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States), September 2013.h jubjNh/URL: }(hURL: h jubjn)}(hhh]h/#https://www.osti.gov/biblio/1136788}(h#https://www.osti.gov/biblio/1136788h j
ubah}(h]h]h]h]h]refurijuhjmh jubh/, }(h, h jubjn)}(hhh]h/doi:10.2172/1136788}(hdoi:10.2172/1136788h jubah}(h]h]h]h]h]refurihttps://doi.org/10.2172/1136788uhjmh jubh/.}(hjSh jubeh}(h]h]h]h]h]uhhKh jubeh}(h]miller-guide-2013ah]hoah]miller_guide_2013ah]h]j(ajjuhhth j-jKubj0)}(hhh](j5)}(hhh]h/peplow_methods_2012}(hhh jCubah}(h]h]h]h]h]jCuhj4h j@ubhL)}(hhh](h/ Douglas Peplow and Larry Wetzel.}(h Douglas Peplow and Larry Wetzel.h jPubh/ }(h h jPubh/Methods for }(hMethods for h jPubh/Detector}(hDetectorh jPubh/ }(h h jPubh/ Placement}(h Placementh jPubh/ and }(h and h jPubh/Analysis}(hAnalysish jPubh/ of }(h of h jPubh/Criticality}(hCriticalityh jPubh/ }(h h jPubh/Accident}(hAccidenth jPubh/ }(h h jPubh/Alarm}(hAlarmh jPubh/ }(hj\h jPubh/Systems}(hSystemsh jPubh/.}(hjSh jPubjNj)}(hhh]h/,Transactions of the American Nuclear Society}(h,Transactions of the American Nuclear Societyh jubah}(h]h]h]h]h]uhjh jPubh/, January 2012.}(h, January 2012.h jPubeh}(h]h]h]h]h]uhhKh j@ubeh}(h]peplow-methods-2012ah]hoah]peplow_methods_2012ah]h]hajjuhhth j-jKubj0)}(hhh](j5)}(hhh]h/peplow_criticality_2010}(hhh jubah}(h]h]h]h]h]jCuhj4h jubhL)}(hhh](h/Douglas}(hDouglash jubjh/E. Peplow and Jr}(hE. Peplow and Jrh jubjh/ Lester M.}(h Lester M.h jubjh/Petrie.}(hPetrie.h jubh/ }(h h jubh/Criticality }(hCriticality h jubh/Accident}(hAccidenth jubh/ }(h h jubh/Alarm}(hAlarmh jubh/ }(h h jubh/System}(hSystemh jubh/ }(h h jubh/Modeling}(hModelingh jubh/ }(h h jubh/Made}(hMadeh jubh/ }(hjh jubh/Easy}(hEasyh jubh/ with }(h with h jubh/SCALE}(hSCALEh jubh/ 6.1.}(h 6.1.h jubjNj)}(hhh]h/,Transactions of the American Nuclear Society}(h,Transactions of the American Nuclear Societyh j=ubah}(h]h]h]h]h]uhjh jubh/, 102(1):297–299, June 2010.}(h, 102(1):297–299, June 2010.h jubjNh/$Publisher: American Nuclear Society.}(h$Publisher: American Nuclear Society.h jubjNh/URL: }(hURL: h jubjn)}(hhh]h/http://epubs.ans.org/?a=10467}(hhttp://epubs.ans.org/?a=10467h jZubah}(h]h]h]h]h]refurijauhjmh jubh/.}(hjSh jubeh}(h]h]h]h]h]uhhKh jubeh}(h]peplow-criticality-2010ah]hoah]peplow_criticality_2010ah]h]hajjuhhth j-jKubj0)}(hhh](j5)}(hhh]h/peplow_criticality_2009}(hhh j~ubah}(h]h]h]h]h]jCuhj4h j{ubhL)}(hhh](h/Douglas}(hDouglash jubjh/E. Peplow and Lester}(hE. Peplow and Lesterh jubjh/ M. Petrie}(h M. Petrieh jubjh/Jr.}(hJr.h jubh/ }(h h jubh/Criticality }(hCriticality h jubh/Accident}(hAccidenth jubh/ }(h h jubh/Alarm}(hAlarmh jubh/ }(h h jubh/System}(hSystemh jubh/ }(hjh jubh/Modeling}(hModelingh jubh/ with }(h with h jubh/SCALE}(hSCALEh jubh/.}(hjSh jubjNj)}(hhh]h/wProceedings of the 2009 International Conference on Advances in Mathematics, Computational Methods, and Reactor Physics}(hwProceedings of the 2009 International Conference on Advances in Mathematics, Computational Methods, and Reactor Physicsh jubah}(h]h]h]h]h]uhjh jubh/, May 2009.}(h, May 2009.h jubjNh/URL: }(hURL: h jubjn)}(hhh]h/:http://inis.iaea.org/Search/search.aspx?orig_q=RN:40080894}(h:http://inis.iaea.org/Search/search.aspx?orig_q=RN:40080894h jubah}(h]h]h]h]h]refurijuhjmh jubh/.}(hjSh jubeh}(h]h]h]h]h]uhhKh j{ubeh}(h]peplow-criticality-2009ah]hoah]peplow_criticality_2009ah]h]hmajjuhhth j-jKubjeh}(h]$bibtex-bibliography-CAAScapability-0ah]h]h]h]uhhKh jhhh!NhNubeh}(h]summaryah]h]summaryah]h]uhh#h h%hhh!h"hMYubeh}(h](!mavric-appendix-a-caas-capabilityheh]h]("mavric appendix a: caas capabilitycaascapabilityeh]h]uhh#h hhhh!h"hKj~}j*hsj}hhsubeh}(h]h]h]h]h]sourceh"uhhcurrent_sourceNcurrent_lineNsettingsdocutils.frontendValues)}(h(N generatorN datestampNsource_linkN
source_urlN
toc_backlinksjfootnote_backlinksK
sectnum_xformKstrip_commentsNstrip_elements_with_classesN
strip_classesNreport_levelK
halt_levelKexit_status_levelKdebugNwarning_streamN tracebackinput_encoding utf-8-siginput_encoding_error_handlerstrictoutput_encodingutf-8output_encoding_error_handlerjQerror_encodingUTF-8error_encoding_error_handlerbackslashreplace
language_codeenrecord_dependenciesNconfigN id_prefixhauto_id_prefixid
dump_settingsNdump_internalsNdump_transformsNdump_pseudo_xmlNexpose_internalsNstrict_visitorN_disable_configN_sourceh"_destinationN
_config_files]pep_referencesNpep_base_url https://www.python.org/dev/peps/pep_file_url_templatepep-%04drfc_referencesNrfc_base_urlhttps://tools.ietf.org/html/ tab_widthKtrim_footnote_reference_spacefile_insertion_enabledraw_enabledKsyntax_highlightlongsmart_quotessmartquotes_locales]character_level_inline_markupdoctitle_xform
docinfo_xformKsectsubtitle_xformembed_stylesheetcloak_email_addressesenvNubreporterNindirect_targets]substitution_defs}substitution_names}refnames}(peplow_criticality_2009]h citation_reference)}(h[peplow_criticality_2009]_h]h/peplow_criticality_2009}(hhh jubah}(h]hmah]hoah]h]h]hjuhjh hMjKubapeplow_criticality_2010]j)}(h[peplow_criticality_2010]_h]h/peplow_criticality_2010}(hhh jubah}(h]hah]hoah]h]h]hjuuhjh hMjKubapeplow_methods_2012]j)}(h[peplow_methods_2012]_h]h/peplow_methods_2012}(hhh jubah}(h]hah]hoah]h]h]hjuhjh hjKubawetzel_caas_2012]j)}(h[wetzel_caas_2012]_h]h/wetzel_caas_2012}(hhh jubah}(h]hah]hoah]h]h]hjuhjh hjKubamiller_guidance_2013]j)}(h[miller_guidance_2013]_h]h/miller_guidance_2013}(hhh jubah}(h]jah]hoah]h]h]hjuhjh hjKubamiller_guide_2013]j)}(h[miller_guide_2013]_h]h/miller_guide_2013}(hhh jubah}(h]j(ah]hoah]h]h]hj:uhjh hjKubamiller_corrected_2013]j)}(h[miller_corrected_2013]_h]h/miller_corrected_2013}(hhh jubah}(h]jEah]hoah]h]h]hjuhjh hjKubaurefids}(h]hajK]jAaj]jaj]jaj]jaj]jaje ]j[ aj
]j
aj]jaj]jaj]jaunameids}(j*hj)j&j_j\jjjjj+j(jjjjjjjyjKjjj j jjjjjjjjj je jj
jyjvjnjjjjjjjj!jjjjjj=j:jjjxjujjjju nametypes}(j*j)Nj_NjNjNj+NjNjNjNjyjj NjjjjNj jjyNjnjNjjj!Njjj=jjxjjuh}(hh%j&h%j\h:hmjhjhjhjjjj(jjEjjjbjjj(jjj.jjjjjKjNjjj jjjjjjjjj je jf j
j
jvjjjjj|jjjjjjjj1jjj:jjj@jujjj{jjjvjNjjjjjjjjj jf jj
jkjjjjju
footnote_refs}
citation_refs}(j]jaj]jaj]jaj]jaj]jaj]jaj]jau
autofootnotes]autofootnote_refs]symbol_footnotes]symbol_footnote_refs] footnotes] citations](j1jjj@jj{jeautofootnote_startKsymbol_footnote_startK
id_countercollectionsCounter}j_KsRparse_messages]transform_messages](h system_message)}(hhh]hL)}(hhh]h/4Hyperlink target "caascapability" is not referenced.}(hhh j9ubah}(h]h]h]h]h]uhhKh j6ubah}(h]h]h]h]h]levelKtypeINFOsourceh"lineKuhj4ubj5)}(hhh]hL)}(hhh]h/,Hyperlink target "fig4a1" is not referenced.}(hhh jTubah}(h]h]h]h]h]uhhKh jQubah}(h]h]h]h]h]levelKtypejNsourceh"lineM'uhj4ubj5)}(hhh]hL)}(hhh]h/,Hyperlink target "fig4a2" is not referenced.}(hhh jnubah}(h]h]h]h]h]uhhKh jkubah}(h]h]h]h]h]levelKtypejNsourceh"lineM/uhj4ubj5)}(hhh]hL)}(hhh]h/,Hyperlink target "tab4a1" is not referenced.}(hhh jubah}(h]h]h]h]h]uhhKh jubah}(h]h]h]h]h]levelKtypejNsourceh"lineM|uhj4ubj5)}(hhh]hL)}(hhh]h/,Hyperlink target "fig4a3" is not referenced.}(hhh jubah}(h]h]h]h]h]uhhKh jubah}(h]h]h]h]h]levelKtypejNsourceh"lineMuhj4ubj5)}(hhh]hL)}(hhh]h/,Hyperlink target "fig4a4" is not referenced.}(hhh jubah}(h]h]h]h]h]uhhKh jubah}(h]h]h]h]h]levelKtypejNsourceh"lineMuhj4ubj5)}(hhh]hL)}(hhh]h/,Hyperlink target "fig4a5" is not referenced.}(hhh jubah}(h]h]h]h]h]uhhKh jubah}(h]h]h]h]h]levelKtypejNsourceh"lineMuhj4ubj5)}(hhh]hL)}(hhh]h/,Hyperlink target "tab4a2" is not referenced.}(hhh jubah}(h]h]h]h]h]uhhKh jubah}(h]h]h]h]h]levelKtypejNsourceh"lineMCuhj4ubj5)}(hhh]hL)}(hhh]h/,Hyperlink target "tab4a3" is not referenced.}(hhh j
ubah}(h]h]h]h]h]uhhKh jubah}(h]h]h]h]h]levelKtypejNsourceh"lineMuhj4ubj5)}(hhh]hL)}(hhh]h/,Hyperlink target "fig4a6" is not referenced.}(hhh j$ubah}(h]h]h]h]h]uhhKh j!ubah}(h]h]h]h]h]levelKtypejNsourceh"lineMLuhj4ubj5)}(hhh]hL)}(hhh]h/,Hyperlink target "fig4a7" is not referenced.}(hhh j>ubah}(h]h]h]h]h]uhhKh j;ubah}(h]h]h]h]h]levelKtypejNsourceh"lineMRuhj4ubetransformerN
decorationNhhub.