-
Design Analysis Major Revision Cover SheetDesign Analysis (Major
Revision)
Analysis No.:’ PM-1083 Revision: 2 0Title: LLRWSF Shielding
Analysis for HIC or Liner StorageECIECR No.: 09-00301 RevisIon:
0
Last Page No. 6 9 / Att. J, Page J4
Station(s): Peach Bottom Component(s): ‘ LLWRSFUnit No.: 2&3
N/ADiscIpline: MECHDescrip. Code/Keyword: ‘5
Safety/OA Class: 1 Non-Safety RelatedSystem Code: 12
N/AStructure: 13 N/A
CONTROLLED DOCUMENT REFERENCES ‘
Document No.: FromiTo Document No.: FromlToPBAPS Caic PM-0750
From
PBAPS Caic PM-0753 From
PBAPS Cab PM-0758 From
PBAPS=Calc PM-0759 From
PBAPS Drawing M-2406 From
is this Design Analysis Safeguards information? ‘ Yes U No If
yes, see SY-AA-101 -106Does this Design Analysis contain Unverified
Assumptions? ‘‘ Yes EJ No If yes, ATI/AR#:This Design Analysis
SUPERCEDES: 18 None in its entirety.Description of Revision (list
affected pages for partials):Preparer:
Paul Reiche /‘J! 0 7Prr,t Name Sign Name OatMethod of Review: 21
Detailed Review Alternate Calcula io s (atta ed) fl Testing
UReviewer: Andy Woodruffe O’/2. 8/2 c;zoC1Print Name S e 0 eReview
Notes: 23 Independent review Peer review
All inputs, assumptions, approaches, numerical analyses, and
results were independently reviewed and ecked.
External Approver: Don Gardner
Exelon Reviewer: 25 j ,4Ti ?- !/P9nfName Sit,e ‘ tYóte
independent 3 Party Review Reqd? Yes U NoExelon Approver:
Pint Name Sign Name Date
-
CC-AA-309Revision 9
Page l7of 17
ATTACHMENT IOwners Acceptance Review Checklist for External
Design Analysis
Page 1 of I
DESIGN ANALYSIS NO. PM-iC 83 REV:
_____
Yes No N/A1. Do assumptions have sufficient rationale? El El2
Are assumptions compatible with the way the plant is operated and
with the
licensing basis?
3. Do the design inputs have sufficient rationale? E1- El ElAre
design inputs correct and reasonable with critical parameters
identified, if
* appropriate?
Are design inputs compatible with the way the plant is operated
and with thelicensing basis?
6. Are Engineering Judgments clearly documented and
justified?
Are Engineering Judgments compatible with the way the plant is
operatedand with the licensing basis?
8 Do the results and conclusions satisfy the purpose and
objective of theDesign Analysis?
Are the results and conclusions compatible with the way the
plant is operatedand with the licensing basis?
10 Does the Design Analysis include the applicable design basis
Eldocumentation?Have any limitations on the use of the results been
identified and transmitted Elto the appropriate organizations?
12. Are there any unverified assumptions? El13 Do all unverified
assumptions have a tracking and closure mechanism in El• place?
14 Have all affected design analyses been documented on the
Affected El El LiYDocuments List (ADL) for the associated
Configuration Change?Do the sources of inputs and analysis
methodology used meet currenttechnical requirements and regulatory
commitments? (If the input sources or
15. analysis methodology are based on an out-of-date methodology
or code, El Eladditional reconciliation may be required if the site
has since committed to amore recent code)
16 Have vendor supporting technical documents and references
(including GE El El* DRFs) been reviewed when necessary?
17 Have margin impacts been identified and documented
appropriately for any El El* negative impacts (Reference
ER-AA-2007)?
EXELON REVIEWER IJ /t1j DATE
______
Pri fSigii / !/ j
-
CALCULATION NO. PM-1083 REV. NO. 0 PAGE NO. 2 of 8
Table of Contents 1.
Purpose.............................................................................................................................................3
2. Inputs
................................................................................................................................................3
3. Assumptions
.....................................................................................................................................4
3.1 Verified Assumptions
................................................................................................................4
3.2 Unverified Assumptions
............................................................................................................4
4.
References........................................................................................................................................4
5. Identification of Computer Programs
................................................................................................5
6. Method of
Analysis............................................................................................................................5
7. Numerical Analysis
...........................................................................................................................5
7.1 Isotopic Analysis
.......................................................................................................................5
7.2 VISED Geometry and MCNP Model Testing
............................................................................6
7.3 HIC Baseline Model
..................................................................................................................6
7.4 Skyshine Analysis
.....................................................................................................................6
7.5 Wall Shielding
Assessment.......................................................................................................7
7.6 Intermittent Dose
Considerations..............................................................................................7
8. Results
..............................................................................................................................................8
8.1 Dose Rates over Filled
Cells.....................................................................................................8
8.2 Skyshine vs. Distance
...............................................................................................................8
8.3 Wall Shielding
Assessment.......................................................................................................8
9. Conclusions
......................................................................................................................................8
Attachments
A. Container MicroShield Modeling, Including Isotopic Tests B.
VISED MCNP Model Illustrations C. Single Container MCNP Modeling
for Shielding Evaluations D. Single Storage Cell MCNP Modeling for
Dose Rate Evaluations over Cell and for Translated
Source Generation E. MCNP Modeling for Skyshine from Central 15
Cells F. MCNP Modeling for Wall Shielding Evaluation G. Spreadsheet
for Dose Evaluations based on MCNP Results H. LLRWSF – Potential
Nearest Resident Distance Measurement I. Reviewer Notes J. Computer
Disclosure Sheets
-
CALCULATION NO. PM-1083 REV. NO. 0 PAGE NO. 3 of 8
1. Purpose The Peach Bottom Atomic Power Station (PBAPS)
Low-Level Radioactive Waste Storage Facility (LLRWSF) is an
existing facility for the storage of radioactive waste that is
either being staged for shipment for disposal, or for waste storage
if a disposal site is not available. The LLRWSF can also be used
for decay in storage before shipment. Shielding analyses were
performed as part of the original design and again in the early
1990’s (See References 1, 2, and 3). In general, those analyses
were for drum storage filling the storage cells, with source terms
considered appropriate at that time. These historical analyses, and
the resulting design are conservative and yielded capabilities and
storage flexibility that can address a wide variety of waste source
configurations and radioactivity levels. However, the LLRWSF
mission, as reanalyzed in this calculation, is specifically for the
storage of Class B & C resin in Energy Solutions 8-120 High
Integrity Containers (HICs). This may be required for an extended
period of time due to disposal site unavailability. These types of
waste, on average, may also contain more radioactivity that assumed
in original design analyses. Shielding analyses in this calculation
are designed to conservatively encompass not only PBAPS waste but
also for the possibility that waste could be transported from the
Limerick Generating Station (LGS) for storage in the LLRWSF.
Analyses are based on container contact dose rates. This will be
the primary operational parameter used to assure that waste will be
stored such that offsite and onsite dose rate limits are met, and
that the ALARA principle is followed.
2. Inputs
1. Based on recent operating experience, a bounding container
from PBAPS is expected to be dewatered resin with a contact dose
rate of 400 R/hr. This provides for some allowance for the use of a
decay in storage period prior to transport to offsite locations.
Bounding containers from LGS are to be treated as dewatered resin
with a contact dose rate of 250 R/hr, which is considered a
reasonable upper bound of immediately shippable waste with a high
Co-60 content. Since PBAPS waste is bounding, all stored waste is
analyzed as having a contact dose rate of 400 R/hr for central
cells. Dose rates for peripheral cells are to be determined based
on wall shielding evaluations.
2. Concrete shielding density is 140 lb/ft3, conservative and
consistent with existing shielding calculations.
3. Dose rates to offsite receptors are controlled by 40CFR190
and limited to 1 mrem/year. This is a small fraction of the
40CFR190 limit for all nuclear fuel cycle activities of 25
mrem/year.
4. Dose rates in areas outside of, but in the immediate vicinity
of the LLRWSF, will be limited to less than 1 mrem/hour, to assure
that occupation doses in the LLRWSF vicinity will be within 10CFR20
limits and ALARA. This is an increase from the 0.25 mrem/hr
criterion used in original design analyses, but consistent with
values used at other Exelon nuclear power stations.
-
CALCULATION NO. PM-1083 REV. NO. 0 PAGE NO. 4 of 8
3. Assumptions
3.1 Verified Assumptions
1. Based on the fact that Co-60 can frequently dominate waste
contact doses, and because it is the most conservative expected
isotope for shielding purposes, a Co-60 only assumption is used.
For shielding purposes, this is more conservative that isotopic
mixes historically used in LLRWSF shielding analyses (See
Attachment A).
2. Storage Cell configuration and dimensions are taken from
Reference 4. 3. Waste container are simplified as simply waste in a
5 foot diameter cylinder filled to a 6 foot
height, with a specific gravity of 0.9, with radioactivity
levels set to yield the identified contact dose limit. This same
configuration and waste strength is then used in shielding
analysis.
4. Waste cells, when filled; contain 4 layers of 4 containers
each. Based on data collected and discussions with operators, the
bottom of waste packages are set at 22.0, 124.5, 227.0, and 329.5
inches above the floor, which is treated as 36 feet below the
bottom of cell shield plugs. Each container is centered in its
quadrant of the cell. No credit is taken for any shielding effect
of the tables that support these packages.
5. Waste cells are modeled based on a nominal 400 R/hr contact
container. The central 15 cells are explicitly modeled in this
calculation. As a reasonable approximation, Skyshine doses are
considered proportional to the average contact dose rate of
upwardly exposed containers, the number of such containers, and,
effectively, the number of cells in use.
6. Based principally on the assumption 5, 15 full (16) container
waste cells on the periphery, and the 5 half cells (south end) are
approximated proportionately in terms Skyshine effects to the
central cells’ results. This assumption is based on experience with
these types of analyses.
7. Assumptions 5 and 6 are considered extremely conservative
with respect to possible filling strategies for the LLRWSF, for the
following reasons:
a. All containers will not be worst case; and b. Only about 4
containers per year would be anticipated to require long term
storage if
Class B and C waste cannot be shipped for disposal. This results
in substantial decay as waste accumulates.
Nonetheless, this assumption is used since dose criteria are met
even in this extreme.
3.2 Unverified Assumptions None.
4. References
1. PBAPS Calculation PM-753, Rev.0, “Bulk Storage Shielding
Analysis – Direct Dose Rate Analysis”
2. PBAPS Calculation PM-758, Rev.0, “Skyshine Dose Rate
Analysis” 3. PBAPS Calculation PM-759, Rev.0, “Dose Rate for
Intermittent Conditions and Summary of
Line Storage Analysis” 4. PBAPS Drawing M-2406, Rev. 9, “Onsite
Radwaste Storage Facility Plan – El. 135’ 6” and
Section B-B Looking West”
-
CALCULATION NO. PM-1083 REV. NO. 0 PAGE NO. 5 of 8
5. Identification of Computer Programs The following computer
programs are used in this calculation, along with spreadsheets. The
programs below are qualified for safety-related use under the
URS-Washington Division QA procedures.
1. MicroShield: A computer program used for evaluation of some
simple geometry problems. This program uses point-kernel
methodology as does the QAD-CGGP computer program used in Reference
1.
2. MCNP: A computer program used for complex geometry
calculations including Skyshine analyses. This program uses Monte
Carlo methodology as does the MORSE-CGA program used in Reference
2
These newer programs use a simpler and more sophisticated
handling of input that, particularly in the MCNP case, facilitates
more detailed input. Methodologies are effectively the same as
previously used in LLRWSF shielding analyses.
6. Method of Analysis The shielding analyses in this calculation
are based on management of LLRWSF loading principally based on
observed container contact dose rates. The steps in these analyses
are as follows:
A. An individual container MicroShield model is developed, along
with selection of isotopic assumptions. Individual container MCNP
models are developed to determine the source strength needed for a
given contact dose rate.
B. An individual cell is modeled, with one, two, three, and four
high stacking of filled layers. This model provides the basis for
calculating dose rates immediately above the cell, and also creates
a “surface source” from a cell that is used in the overall building
Skyshine analysis.
C. A building model is prepared that models the 35 cells, with
the central 15 having the “surface source” developed in step B.
D. A wall model is developed reflecting a filled cell to
determine contact dose rate limitation for containers in the
outside cells.
E. A spreadsheet is developed to turn the result discussed above
into actual dose rate calculations. This is because the MCNP models
are run on a dose rate per photon/second basis.
F. The dose results for shielding from historically analyzed
transient conditions are adjusted based on new source terms and
container geometries.
G. The above analyses are used to establish limitations on
contact dose rates for containers placed in central and peripheral
cells.
H. The above analyses also provide bases for assuring that dose
rates limits outside of the building are met, and possibly a basis
for Skyshine prediction that could supplement and support
predictions based on exterior Thermoluminescent Dosimeters
(TLDs).
7. Numerical Analysis
7.1 Isotopic Analysis Attachment A contains a MicroShield based
justification for the use of an all-Co-60 based source term for
shielding, rather than the historically used mix from Reference 1
and 2. This assumption is conservative, but not inordinately
so.
-
CALCULATION NO. PM-1083 REV. NO. 0 PAGE NO. 6 of 8
7.2 VISED Geometry and MCNP Model Testing VISED is a visual
editor for MCNP that provides visual images of model geometries,
and also allows particle test illustrating source locations and
scatter events. Attachment B contains these images and are
discussed below with respect to individual models.
7.3 HIC Baseline Model Attachment C contains an MCNP run for
modeling contact dose rates for a single HIC with a homogeneous
waste source 6 feet high and 5 feet in diameter. Page 1 of
Attachment B, top figure, shows the geometry. Surface crossing
detectors are used at 1 inch (nominal contact) and 1 meter from the
waste source. The horizontal lines segment the surface axially to
allow the assessment of axial variation and to assure that peaking
is considered. Surface crossing detectors can be thought of as
tallying particles crossing the surface (angle corrected), dividing
by the surface area to determine a fluence, and then applying a
fluence to dose conversion factor as a function of photon energy.
The net result is a rem/hr per photon/sec value. All MCNP analyses
are run on a single photon per second source basis. This run is
used to establish the amount of activity (photons/sec) necessary to
yield a 400 R/hr contact container. This derivation is shown in the
Spreadsheet in Attachment G.
7.4 Skyshine Analysis Basic Cell Geometry The modeling of the
LLRWSF is a classical “deep penetration” problem, both because of
the 2 foot thick concrete shield overlying the storage cells, and
the thickness of air through which the Skyshine analysis is being
performed. Therefore, for variance reduction, this analysis is
performed in two steps. The first step is to model an individual
interior cell. MCNP runs are made for 1, 2, 3, and 4 layer stacking
and included in Attachment D. For operational information, doses
are also calculated at 1 inch above the cell cover. Page 1 of
Attachment B, bottom figure, shows the geometry. Page 2 of
Attachment B, top figure demonstrates that source particles, in
fact, originate in the containers. Page 2 of Attachment B, bottom
figure, demonstrates that behavior of the photon population being
scattered by walls, and being transported through the cell cover.
Additionally, a surface source is written to a file that collects
data of particles exiting through top of the cell in terms of
particle weight, direction, and energy. Lateral surfaces of the
cell are one half of the wall thickness, and mirrored to simulate
the effect of an infinite array for purpose of the cell being
modeled. This surface source is used in the actual Skyshine
analysis. Basic Skyshine Geometry The LLRWSF is simplified as shown
in Page 3 of Attachment B, top figure. The DAW, adjacent
compartment and the truck bay are conservatively ignored to
maximize dose rates in the immediate vicinity of the LLRWSF. Page 3
of Attachment B, bottom figure demonstrates that the surface source
generated in the cell model is properly translated over the 15
interior cells. The walls and roof above the cell are consistent
with the design of these bodies, as previously analyzed. Page 4 of
Attachment B, top figure shows the Skyshine model geometry. In
essence, spherical surfaces are place concentrically around the
waste center (also the modeled building center). Two types of
surface crossing detectors are used. The first is based on photons
crossing the ground surface, which is segmented by the concentric
spheres. The detector distance is taken as nominally at
-
CALCULATION NO. PM-1083 REV. NO. 0 PAGE NO. 7 of 8
the midpoint distance between adjacent spheres. The second
surface crossing detector is the surface of the spheres themselves,
from the ground up to the waste height. The MCNP run for Skyshine
assessment is in Attachment E, addressing the central 15 cells.
This derivation of doses as a function of distance is shown in the
Spreadsheet in Attachment G. Variance Reduction MCNP is a Monte
Carlo analysis tool, and variance reduction techniques must be used
to obtain acceptable statistics for complex geometries such as
modeled in this analysis.
The first technique used is importance biasing through the
shield cap cell atop the storage vaults. This shield has three inch
segments, and importances such that particle splitting is applied
to maintain the photon population through the shield.
Importances are also used in the Skyshine model to maintain
photon population as a function of distance.
A source translation is also performed to write a surface source
above a storage vault, which is then replicated above other vaults,
and then re-run for Skyshine analysis.
Dose Conversion Factor Dose conversion factors are taken from
ANSI/ANS-6.1.1 (1977) (Table H.2 of the MCNP Manual)
7.5 Wall Shielding Assessment Attachment F, shows the wall
shielding model used to determine a limit on container contact such
that a 1 mrem/hr dose rate would not be exceeded. Geometrically,
this model uses the same arrangement of containers as the Skyshine
analyses, but models a periphery cell, where the waste is somewhat
further from the outside surface due to crane travel limitations.
Here, the outside wall is segmented for variance reduction
importance utilization. Attachment F contains the MCNP run for the
wall shielding assessment. The Attachment G spreadsheet show the
derivation of contact dose limits to assure that contact dose rates
outside of the LLRWSF are less than or equal to 1 mrem/hr. The
Skyshine contribution in the immediate vicinity of the IRSF is
calculated to be negligible.
7.6 Intermittent Dose Considerations Reference 3 evaluated and
provided the following results for various intermittent and upset
conditions for the LLRWSF.
-
CALCULATION NO. PM-1083 REV. NO. 0 PAGE NO. 8 of 8
Because of the new all Co-60 assumption and the increase in
contact dose rate these results should be multiplied by 4.5 (an
extremely conservative factor given only an approximate increase in
contact dose rate of a factor of 2) to increase the source term
effects as detailed in Attachment A. These adjusted results remain
acceptable.
8. Results
8.1 Dose Rates over Filled Cells Dose rates directly above cells
are effectively a function of fill height, number of upwardly
exposed containers and stacked height. This analysis has shown that
it is unlikely that contact dose rates on top of the vaults will
exceed 100 mR/hr. Even when filled four levels high with each HIC
at 400 R/hr the average dose rate above the vaults is slightly
under 100 mR/hr, with localized hot spots marginally over 100
mR/hr.
8.2 Skyshine vs. Distance The spreadsheet in Attachment G
evaluates Skyshine based on the assumption that the central 15
cells are filled with 400 R/hr contact dose rate containers, and
the peripheral cells are filled with 125 R/hr contact dose rate
containers.
8.3 Wall Shielding Assessment The Attachment G spreadsheet shows
that the dose rate on the lower radiation side of the shield wall
remains ~ 1mR/hr even with the outer vaults filled with HICs up to
125 R/hr on contact.
9. Conclusions Analyses in this calculation demonstrate the
following:
1. The central cells could be filled with 400 R/hr contact dose
rate containers, and the peripheral cells filled with 125 R/hr
contact dose rate containers, and the dose rates outside of the
LLRWSF, both onsite and offsite, would meet dose rate limits and
goals. In fact, this is a likely upper bound on expected
containers, and substantial decay in storage would reduce the
likelihood of approaching the calculated exterior dose rates.
2. Dose rates immediately over a cell filled with 400 R/hr
contact containers could approach high radiation conditions (>
100 mrem/hr). Cells filled to less than 4 layers high would have
significantly lower dose rates over the cells, and be expected to
have proportionally lower external dose effects.
3. Intermittent doses are higher than those calculated in
PM-0759 but are still acceptable. 4. The calculated dose rate to
the nearest offsite residence location is less than 1 mrem/year,
even
with the bounding source assumptions. 5. Realistically, dose
rates outside of the building from Skyshine are unlikely to exceed
10% of the
calculated values, given more realistic HIC contact dose rates
and decay effects. 6. Dose rates in the immediate vicinity of the
walls will be proportional to container contact dose
rates in the near vicinity, and are considered predictive,
though multiple containers are assumed to contribute.
7. Results from this calculation could be used to predict doses
to 40CFR190 receptors, but in any case these doses are expected to
be small compared to dose goals and very small compared to
regulatory limits.
-
MicroShield v5.05 (5.05-00211)Raytheon Engineers &
Constructors
Page : 1DOS File: PBAPS.MS5Run Date: April 17, 2009Run Time:
11:51:22 AMDuration: 00:00:15
File Ref: __________________Date: __________________
By: __________________Checked: __________________
Case Title: Case 1Description: Case 1
Geometry: 7 - Cylinder Volume - Side Shields
Source DimensionsHeight 182.88 cm 6 ftRadius 76.2 cm 2 ft 6.0
in
Dose PointsX Y Z
# 1 78.74 cm 91.44 cm 0 cm2 ft 7.0 in 3 ft 0.0 in
ShieldsShield Name Dimension Material DensitySource 3.34e+06 cm³
Water 0.9Transition Air 0.00122Air Gap Air 0.00122
Source InputGrouping Method : Linear Energy
Number of Groups : 25Lower Energy Cutoff : 0.015Photons <
0.015 : Excluded
Library : GroveNuclide curies becquerels µCi/cm³ Bq/cm³Ag-110m
4.5900e+000 1.6983e+011 1.3759e+000 5.0908e+004Ba-137m 6.6620e+001
2.4649e+012 1.9970e+001 7.3889e+005Co-58 1.1453e+001 4.2376e+011
3.4332e+000 1.2703e+005Co-60 1.5100e+002 5.5869e+012 4.5263e+001
1.6747e+006Cs-134 6.0870e+001 2.2522e+012 1.8246e+001
6.7511e+005Cs-137 7.0423e+001 2.6056e+012 2.1110e+001
7.8107e+005Mn-54 7.4431e+000 2.7539e+011 2.2311e+000
8.2552e+004Nb-95 2.5643e+000 9.4879e+010 7.6868e-001
2.8441e+004Sr-92 1.0112e+002 3.7415e+012 3.0313e+001
1.1216e+006Zn-65 2.6075e+002 9.6478e+012 7.8162e+001
2.8920e+006
BuildupThe material reference is : Source
Integration ParametersRadial 15Circumferential 36Y Direction
(axial) 20
ResultsEnergy Activity Fluence Rate Fluence Rate Exposure Rate
Exposure RateMeV photons/sec MeV/cm²/sec MeV/cm²/sec mR/hr
mR/hr
No Buildup With Buildup No Buildup With Buildup0.0328 1.980e+11
2.896e+03 7.504e+03 2.207e+01 5.718e+010.0989 1.624e+07 1.319e+00
8.827e+00 2.018e-03 1.350e-020.2415 1.111e+11 2.898e+04 1.220e+05
5.314e+01 2.238e+020.2769 7.973e+08 2.500e+02 9.663e+02 4.689e-01
1.813e+000.3686 2.070e+08 9.611e+01 3.194e+02 1.862e-01
6.189e-01
Calculation No. PM-1083, Rev. 0 Attachment A, (Page A1 of
A2)
-
Page : 2DOS File: PBAPS.MS5Run Date: April 17, 2009Run Time:
11:51:22 AMDuration: 00:00:15
Energy Activity Fluence Rate Fluence Rate Exposure Rate Exposure
RateMeV photons/sec MeV/cm²/sec MeV/cm²/sec mR/hr mR/hr
No Buildup With Buildup No Buildup With Buildup0.4313 1.308e+11
7.560e+04 2.333e+05 1.479e+02 4.564e+020.5079 4.422e+11 3.217e+05
9.211e+05 6.314e+02 1.808e+030.5672 5.368e+11 4.566e+05 1.248e+06
8.938e+02 2.443e+030.6342 4.602e+12 4.591e+06 1.199e+07 8.929e+03
2.332e+040.6948 6.061e+10 6.894e+04 1.736e+05 1.331e+02
3.351e+020.7944 2.260e+12 3.120e+06 7.465e+06 5.940e+03
1.421e+040.8204 7.122e+11 1.030e+06 2.436e+06 1.954e+03
4.620e+030.9016 1.813e+11 3.010e+05 6.890e+05 5.639e+02
1.291e+030.9534 1.349e+11 2.431e+05 5.461e+05 4.516e+02
1.015e+031.0386 2.252e+10 4.603e+04 1.005e+05 8.429e+01
1.841e+021.1155 4.896e+12 1.112e+07 2.373e+07 2.009e+04
4.286e+041.1726 5.735e+12 1.403e+07 2.946e+07 2.507e+04
5.265e+041.3325 5.587e+12 1.653e+07 3.334e+07 2.868e+04
5.784e+041.3835 3.477e+12 1.088e+07 2.169e+07 1.870e+04
3.728e+041.4982 2.898e+10 1.021e+05 1.989e+05 1.719e+02
3.348e+021.5623 2.004e+09 7.518e+03 1.447e+04 1.251e+01
2.408e+011.6747 2.275e+09 9.473e+03 1.789e+04 1.545e+01
2.919e+01
TOTALS: 2.912e+13 6.296e+07 1.344e+08 1.125e+05 2.410e+05
Calculation No. PM-1083, Rev. 0 Attachment A, (Page A2 of
A2)
paul.reichertText BoxOf the above 241 R/hr dose, Co-60
contributes ~45%. Fission product contribute 1/3, and based on
other mixes evaluated, cannot dependably displace Co-60. Therefore,
an all Co-60 assumption will be used. This conclusion also applies
to otheractivation products.
Since the Co-60 fraction from the 1.17 and 1.33 MEV photons is
1.11E4 mR/hr, and the design analysis is for 400 R/hr,
forconservatism, the source is increased by 400/111 or 3.6.
Additionally the new containers have a 5 foot diameter, and
thecontainers analyzed in PM-759 had a 4.5 foot diameter.
Therefore, the total correction is 3.6*(5/4.5)^2 = 4.5.
-
Single HIC MCNP Geometry: Nominal Container with waste treated
as 5 feet in diameter and 6 feet high. Left view is section through
a container with surface at 1 inch from waste and at 1 meter. The
former is used as a contact. To the right is a plan view. The
detecting surface are on the waste surface, 1 inch, and 1 meter
away. Doses are determined based on particles crossing these
surfaces. The horizontal lines are used to axially segment these
surfaces to calculate axial peaking. This MCNP run forms the basis
for determine photons/sec required for a 400 R/hr contact dose
container.
Single PBAPS LLRWSF Waste Cell MCNP Geometry: Model of one of
the identical central 15 cells. Containers are shown in the left
section view, positioned as provided by the support tables, whose
mass is not credited. The containers in the right, plan view are
positioned as centered in each quadrant. This left view is colored
by importances, illustrating how splitting is used to keep the
particle population us through the overlying 2 foot thick concrete
top shield. A top of this shield is used as to create a surface
source file recording the position, direction, energy, and particle
weight for photons leaving the top. The source is used later to
perform skyshine analyses. Separate runs are made and dose rates
above the cell calculated for 1, 2, 3, and 4 layers of 4
containers.
Calculation No. PM-1083, Rev. 0 Attachment B (Page B1 of B4)
-
Single PBAPS LLRWSF Waste Cell Source Particle Display: This
display verifies that the source particles are, in fact, started in
the waste packages.
Single PBAPS LLRWSF Waste Cell Collision Display: This display
shows photon collisions, in the containers, in surrounding
shielding, and to a limited extent, in the air in the cell. The
display depth is selected so that back wall collisions are not
shown.
Calculation No. PM-1083, Rev. 0 Attachment B (Page B2 of B4)
-
PBAPS LLRWSF Simplified Geometry for Skyshine Analysis (Inside
Building): These section and plan views show the LLRWSF, simplified
to only include the storage cell area. This maximizes the dose
rates calculated near the building
Surface Source Effect Illustration While cells and walls are
shown, the effect of the filled cells is based on an array of the
surface sources derived in the cell model. The plan view to the
left is below the two foot thick cell shield caps. The plan view to
the right shows a source particle display that illustrates the
source array effects. The net effect of this array is checked to
confirm that it is consistent with the results for doses over
individual cells.
Calculation No. PM-1083, Rev. 0 Attachment B (Page B3 of B4)
-
PBAPS LLRWSF Simplified Geometry for Skyshine Analysis (Outside
Building): To calculate skyshine dose effects a cylinder and a
series of spherical surfaces surround the building. These surfaces
are used both for particle splitting and as surface crossing
detectors. The horizontal ground surfaces is segmented by the
cylinder and spheres to yield average doses in rings surrounding
the LLRWSF. The spherical surfaces are also used, up to a limited
height, to calculated doses vs. distances .
Calculation No. PM-1083, Rev. 0 Attachment B (Page B4 of B4)
-
HIC.o 5/31/2009
Thread Name & Version = mcnp5_RSICC, 4.23_
._ _ _ ._ ._ |_| | | (_ | | |_) _)
|+--------------------------------------------------------------------+|
This program was prepared by the Regents of the University of
||California at Los Alamos National Laboratory (the University)
under || contract number W-7405-ENG-36 with the U.S. Department of
Energy ||(DoE). The University has certain rights in the program
pursuant to|| the contract and the program should not be copied or
distributed || outside your organization. All rights in the program
are reserved ||by the DoE and the University. Neither the U.S.
Government nor the || University makes any warranty, express or
implied, or assumes any || liability or responsibility for the use
of this software.
|+--------------------------------------------------------------------+
1mcnp version 5 ld=07-10-03 05/31/09
11:36:12*************************************************************************
probid = 05/31/09 11:36:12
inp = HIC outp = HIC.o
1- c Created on: Thursday, May 14, 2009 at 23:162- 6 280 -0.9
-21 11 -12 $ Container3- 7 204 -0.00122 21 -22 11 -12 $ Air Around
Container to 1in4- 8 204 -0.00122 22 -23 11 -12 $ Air Around
Container to 1m5- 31 0 -11 :12 :23 $ Void6-7- 11 pz 0 $ Bottom of
Container8- 12 pz 182.88 $ Top of Container9- 13 pz 20.32 $
Segmenter
10- 14 pz 40.64 $ Segmenter11- 15 pz 60.96 $ Segmenter12- 16 pz
81.28 $ Segmenter13- 17 pz 101.6 $ Segmenter14- 18 pz 121.92 $
Segmenter15- 19 pz 142.24 $ Segmenter16- 20 pz 162.56 $
Segmenter17- 21 cz 76.2 $ Container Radius18- 22 cz 78.74 $ Contact
(1 inch) Receptor19- 23 cz 176.2 $ Receptor at 1 meter20-21- mode
p22- m280 1001.60c -0.111915 $Water (density of 1 assumed)
warning. neutron table inconsistent with mode will be
ignored.23- 8016.60c -0.888085
warning. neutron table inconsistent with mode will be
ignored.24- m228 1001.60c -0.005558 $concrete (ordinary with
ENDF-VI)
warning. material 228 is not used in the problem.25- 8016.60c
-0.498076 11023.60c -0.017101 12000.60c -0.00256526- 13027.60c
-0.045746 14000.60c -0.315092 16000.60c -0.00128327- 19000.60c
-0.019239 20000.60c -0.082941 26054.60c -0.00070728- 26056.60c
-0.01139 26057.60c -0.000265 26058.60c -3.6e-00529- m204 7014.60c
-0.755636 $air (US S. Atm at sea level)
warning. neutron table inconsistent with mode will be
ignored.30- 8016.60c -0.231475 18000.59c -0.012889
warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.
31- imp:p 1 2r 0 $ 6, 31
1Calculation No. PM-1083, Rev. 0 Attachment C (Page C1 of
C9)
-
HIC.o 5/31/2009
32- phys:p 100 1 033- sdef erg=d5 axs=0 0 1 rad=d1 ext=d334- si1
0 76.235- sp1 -21 136- si3 0 182.8837- sp3 0.0 1.038- si5 l 1.1730
1.332539- sp5 1 140- f2:p 21 22 23
warning. without bremsstrahlung, flux estimates will be low.41-
fs2 -13 -14 -15 -16 -17 -18 -19 -2042- f22:p 11 12
warning. without bremsstrahlung, flux estimates will be low.43-
de2 0.01 0.03 0.05 0.07 0.10 0.15 0.20 &44- 0.25 0.30 0.35 0.40
0.45 0.50 0.55 &45- 0.60 0.65 0.70 0.80 1.00 1.40 1.80 &46-
2.20 2.60 2.80 3.25 3.75 4.25 4.75 &47- 5.00 5.25 5.75 6.25
6.75 7.50 9.00 &48- 11.0 13.0 15.049- df2 LIN 3.96E-06 5.82E-07
2.90E-07 2.58E-07 2.83E-07 3.79E-07 5.01E-07 &50- 6.31e-07
7.59E-07 8.78E-07 9.85E-07 1.08E-06 1.17E-06 1.27E-06 &51-
1.36e-06 1.44E-06 1.52E-06 1.68E-06 1.98E-06 2.51E-06 2.99E-06
&52- 3.42e-06 3.82E-06 4.01E-06 4.41E-06 4.83E-06 5.23E-06
5.60E-06 &53- 5.80e-06 6.01E-06 6.37E-06 6.74E-06 7.11E-06
7.66E-06 8.77E-06 &54- 1.03e-05 1.18E-05 1.33E-0555- de22 0.01
0.03 0.05 0.07 0.10 0.15 0.20 &56- 0.25 0.30 0.35 0.40 0.45
0.50 0.55 &57- 0.60 0.65 0.70 0.80 1.00 1.40 1.80 &58- 2.20
2.60 2.80 3.25 3.75 4.25 4.75 &59- 5.00 5.25 5.75 6.25 6.75
7.50 9.00 &60- 11.0 13.0 15.061- df22 LIN 3.96E-06 5.82E-07
2.90E-07 2.58E-07 2.83E-07 3.79E-07 5.01E-07 &62- 6.31e-07
7.59E-07 8.78E-07 9.85E-07 1.08E-06 1.17E-06 1.27E-06 &63-
1.36e-06 1.44E-06 1.52E-06 1.68E-06 1.98E-06 2.51E-06 2.99E-06
&64- 3.42e-06 3.82E-06 4.01E-06 4.41E-06 4.83E-06 5.23E-06
5.60E-06 &65- 5.80e-06 6.01E-06 6.37E-06 6.74E-06 7.11E-06
7.66E-06 8.77E-06 &66- 1.03e-05 1.18E-05 1.33E-0567- nps
25000000
1cells print table 60
atom gram photoncell mat density density volume mass pieces
importance
1 6 280 9.02777E-02 9.00000E-01 3.33600E+06 3.00240E+06 1
1.0000E+002 7 204 5.05147E-05 1.22000E-03 2.26107E+05 2.75850E+02 1
1.0000E+003 8 204 5.05147E-05 1.22000E-03 1.42751E+07 1.74157E+04 1
1.0000E+004 31 0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0
0.0000E+00
total 1.78372E+07 3.02009E+06
minimum source weight = 1.0000E+00 maximum source weight =
1.0000E+00
**************************************************** Random
Number Generator = 1 ** Random Number Seed = 19073486328125 **
Random Number Multiplier = 19073486328125 ** Random Number Adder =
0 ** Random Number Bits Used = 48 *
2Calculation No. PM-1083, Rev. 0 Attachment C (Page C2 of
C9)
-
HIC.o 5/31/2009
* Random Number Stride = 152917
****************************************************
8 warning messages so far.1cross-section tables print table
100
table length
tables from file mcplib04
1000.04p 1898 ENDF/B-VI Release 8 Photoatomic Data for 1-H mat
100 02/07/037000.04p 3194 ENDF/B-VI Release 8 Photoatomic Data for
7-N mat 700 02/07/038000.04p 3272 ENDF/B-VI Release 8 Photoatomic
Data for 8-O mat 800 02/07/03
18000.04p 4696 ENDF/B-VI Release 8 Photoatomic Data for 18-AR
mat1800 02/07/03
total 13060
warning. simple physics turned on for photons > 100 mev.
***********************************************************************************************************************dump
no. 1 on file runtpg nps = 0 coll = 0 ctm = 0.00 nrn = 0
9 warning messages so far.1problem summary
run terminated when 25000000 particle histories were done.+
05/31/09 11:50:54
c Created on: Thursday, May 14, 2009 at 23:16 probid = 05/31/09
11:36:120
photon creation tracks weight energy photon loss tracks weight
energy(per source particle) (per source particle)
source 25000000 1.0000E+00 1.2528E+00 escape 12471779 4.9887E-01
3.2559E-01energy cutoff 0 0. 0.time cutoff 0 0. 0.
weight window 0 0. 0. weight window 0 0. 0.cell importance 0 0.
0. cell importance 0 0. 0.weight cutoff 0 0. 0. weight cutoff 0 0.
0.e or t importance 0 0. 0. e or t importance 0 0. 0.dxtran 0 0. 0.
dxtran 0 0. 0.forced collisions 0 0. 0. forced collisions 0 0.
0.exp. transform 0 0. 0. exp. transform 0 0. 0.from neutrons 0 0.
0. compton scatter 0 0. 9.0258E-01bremsstrahlung 0 0. 0. capture
12536219 5.0145E-01 2.4516E-02p-annihilation 15794 6.3176E-04
3.2283E-04 pair production 7897 3.1588E-04 4.1284E-04photonuclear 0
0. 0. photonuclear abs 0 0. 0.electron x-rays 0 0. 0.1st
fluorescence 101 4.0400E-06 1.1907E-082nd fluorescence 0 0. 0.
total 25015895 1.0006E+00 1.2531E+00 total 25015895 1.0006E+00
1.2531E+00
number of photons banked 7897 average time of (shakes)
cutoffsphoton tracks per source particle 1.0006E+00 escape
3.9186E-01 tco 1.0000E+33photon collisions per source particle
8.8793E+00 capture 3.3692E-01 eco 1.0000E-03total photon collisions
221983479 capture or escape 3.6432E-01 wc1 -5.0000E-01
any termination 3.6421E-01 wc2 -2.5000E-01
3Calculation No. PM-1083, Rev. 0 Attachment C (Page C3 of
C9)
-
HIC.o 5/31/2009
computer time so far in this run 15.07 minutes maximum number
ever in bank 1computer time in mcrun 14.69 minutes bank overflows
to backup file 0source particles per minute 1.7021E+06random
numbers generated 2964429671 most random numbers used was 84853 in
history 3810471
range of sampled source weights = 1.0000E+00 to
1.0000E+001photon activity in each cell print table 126
tracks population collisions collisions number flux average
averagecell entering * weight weighted weighted track weight track
mfp
(per history) energy energy (relative) (cm)
1 6 25017112 25007891 221865993 8.8746E+00 4.7252E-01 4.7252E-01
1.0000E+00 1.0151E+012 7 9029607 9005026 4612 1.8448E-04 6.6245E-01
6.6245E-01 1.0000E+00 9.9804E+033 8 8920886 8913141 112874
4.5150E-03 6.4466E-01 6.4466E-01 1.0000E+00 9.8396E+03
total 42967605 42926058 221983479 8.8793E+00
1tally 2 nps = 25000000tally type 2 particle flux averaged over
a surface.tally for photons
this tally is modified by a dose function.
areassurface: 21 22 23
segment1 9.72878E+03 1.00531E+04 2.24962E+042 9.72878E+03
1.00531E+04 2.24962E+043 9.72878E+03 1.00531E+04 2.24962E+044
9.72878E+03 1.00531E+04 2.24962E+045 9.72878E+03 1.00531E+04
2.24962E+046 9.72878E+03 1.00531E+04 2.24962E+047 9.72878E+03
1.00531E+04 2.24962E+048 9.72878E+03 1.00531E+04 2.24962E+049
9.72878E+03 1.00531E+04 2.24962E+04
surface 21segment: -13
8.15263E-12 0.0018
surface 21segment: 13 -14
1.01840E-11 0.0016
surface 21segment: 13 14 -15
1.06759E-11 0.0016
surface 21segment: 13 14 15 -16
1.08137E-11 0.0016
surface 21segment: 13 14 15 16 -17
1.08174E-11 0.0016
4Calculation No. PM-1083, Rev. 0 Attachment C (Page C4 of
C9)
-
HIC.o 5/31/2009
surface 21segment: 13 14 15 16 17 -18
1.07873E-11 0.0016
surface 21segment: 13 14 15 16 17 18 -19
1.06814E-11 0.0016
surface 21segment: 13 14 15 16 17 18 19 -20
1.01906E-11 0.0016
surface 21segment: 13 14 15 16 17 18 19 20
8.14002E-12 0.0019
surface 22segment: -13
6.64658E-12 0.0016
surface 22segment: 13 -14
8.68549E-12 0.0014
surface 22segment: 13 14 -15
9.23856E-12 0.0014
surface 22segment: 13 14 15 -16
9.42309E-12 0.0014
surface 22segment: 13 14 15 16 -17
9.45131E-12 0.0014
surface 22segment: 13 14 15 16 17 -18
9.42503E-12 0.0014
surface 22segment: 13 14 15 16 17 18 -19
9.26245E-12 0.0014
surface 22segment: 13 14 15 16 17 18 19 -20
8.69113E-12 0.0014
surface 22segment: 13 14 15 16 17 18 19 20
6.63099E-12 0.0016
surface 23segment: -13
1.44668E-12 0.0016
surface 23segment: 13 -14
5Calculation No. PM-1083, Rev. 0 Attachment C (Page C5 of
C9)
-
HIC.o 5/31/2009
1.64225E-12 0.0015
surface 23segment: 13 14 -15
1.80614E-12 0.0014
surface 23segment: 13 14 15 -16
1.89595E-12 0.0014
surface 23segment: 13 14 15 16 -17
1.92995E-12 0.0014
surface 23segment: 13 14 15 16 17 -18
1.89849E-12 0.0014
surface 23segment: 13 14 15 16 17 18 -19
1.80312E-12 0.0014
surface 23segment: 13 14 15 16 17 18 19 -20
1.64810E-12 0.0015
surface 23segment: 13 14 15 16 17 18 19 20
1.45179E-12 0.0016
===================================================================================================================================
results of 10 statistical checks for the estimated answer for
the tally fluctuation chart (tfc) bin of tally 2
tfc bin --mean-- ---------relative error--------- ----variance
of the variance---- --figure of merit-- -pdf-behavior behavior
value decrease decrease rate value decrease decrease rate value
behavior slope
desired random
-
HIC.o 5/31/2009
relative error from zero tallies = 0.0011 relative error from
nonzero scores = 0.0015
number of nonzero history tallies = 759039 efficiency for the
nonzero tallies = 0.0304history number of largest tally = 19679
largest unnormalized history tally = 4.86433E-05(largest
tally)/(average tally) = 6.14240E+02 (largest tally)/(avg nonzero
tally)= 1.86493E+01
(confidence interval shift)/mean = 0.0000 shifted confidence
interval center = 8.14007E-12
if the largest history score sampled so far were to occur on the
next history, the tfc bin quantities would change as follows:
estimated quantities value at nps value at nps+1
value(nps+1)/value(nps)-1.
mean 8.14002E-12 8.14022E-12 0.000025relative error 1.85206E-03
1.85218E-03 0.000063variance of the variance 5.94771E-05
5.94870E-05 0.000166shifted center 8.14007E-12 8.14007E-12
0.000000figure of merit 1.98484E+04 1.98459E+04 -0.000126
the 100 largest history tallies appear to have a maximum value
of about 4.86433E-05the large score tail of the empirical history
score probability density function appears to have no unsampled
regions.
fom = (histories/minute)*(f(x) signal-to-noise ratio)**2 =
(1.702E+06)*( 1.080E-01)**2 = (1.702E+06)*(1.166E-02) =
1.985E+04
1tally 22 nps = 25000000tally type 2 particle flux averaged over
a surface.tally for photons
this tally is modified by a dose function.
areassurface: 11 12
9.75353E+04 9.75353E+04
surface 113.58407E-12 0.0009
surface 123.59290E-12 0.0009
===================================================================================================================================
results of 10 statistical checks for the estimated answer for
the tally fluctuation chart (tfc) bin of tally 22
tfc bin --mean-- ---------relative error--------- ----variance
of the variance---- --figure of merit-- -pdf-behavior behavior
value decrease decrease rate value decrease decrease rate value
behavior slope
desired random
-
HIC.o 5/31/2009
----- estimated confidence intervals: -----
estimated asymmetric confidence interval(1,2,3 sigma):
3.5807E-12 to 3.5874E-12; 3.5774E-12 to 3.5908E-12; 3.5740E-12 to
3.5941E-12estimated symmetric confidence interval(1,2,3 sigma):
3.5807E-12 to 3.5874E-12; 3.5774E-12 to 3.5908E-12; 3.5740E-12 to
3.5941E-12
1analysis of the results in the tally fluctuation chart bin
(tfc) for tally 22 with nps = 25000000 print table 160
normed average tally per history = 3.58407E-12 unnormed average
tally per history = 3.49573E-07estimated tally relative error =
0.0009 estimated variance of the variance = 0.0000relative error
from zero tallies = 0.0005 relative error from nonzero scores =
0.0008
number of nonzero history tallies = 3258038 efficiency for the
nonzero tallies = 0.1303history number of largest tally = 3414
largest unnormalized history tally = 4.86433E-05(largest
tally)/(average tally) = 1.39150E+02 (largest tally)/(avg nonzero
tally)= 1.81343E+01
(confidence interval shift)/mean = 0.0000 shifted confidence
interval center = 3.58408E-12
if the largest history score sampled so far were to occur on the
next history, the tfc bin quantities would change as follows:
estimated quantities value at nps value at nps+1
value(nps+1)/value(nps)-1.
mean 3.58407E-12 3.58409E-12 0.000006relative error 9.33125E-04
9.33136E-04 0.000012variance of the variance 1.40366E-05
1.40369E-05 0.000017shifted center 3.58408E-12 3.58408E-12
0.000000figure of merit 7.81911E+04 7.81892E+04 -0.000024
the 100 largest history tallies appear to have a maximum value
of about 4.86433E-05the large score tail of the empirical history
score probability density function appears to have no unsampled
regions.
fom = (histories/minute)*(f(x) signal-to-noise ratio)**2 =
(1.702E+06)*( 2.143E-01)**2 = (1.702E+06)*(4.594E-02) =
7.819E+04
1status of the statistical checks used to form confidence
intervals for the mean for each tally bin
tally result of statistical checks for the tfc bin (the first
check not passed is listed) and error magnitude check for all
bins
2 passed the 10 statistical checks for the tally fluctuation
chart bin resultpassed all bin error check: 27 tally bins all have
relative errors less than 0.10 with no zero bins
22 passed the 10 statistical checks for the tally fluctuation
chart bin resultpassed all bin error check: 2 tally bins all have
relative errors less than 0.10 with no zero bins
the 10 statistical checks are only for the tally fluctuation
chart bin and do not apply to other tally bins.
1tally fluctuation charts
tally 2 tally 22nps mean error vov slope fom mean error vov
slope fom
2048000 8.1160E-12 0.0064 0.0007 10.0 20287 3.5815E-12 0.0033
0.0002 10.0 783234096000 8.1082E-12 0.0046 0.0004 10.0 19934
3.5845E-12 0.0023 0.0001 10.0 779826144000 8.1161E-12 0.0037 0.0002
10.0 19909 3.5863E-12 0.0019 0.0001 10.0 776258192000 8.1132E-12
0.0032 0.0002 10.0 19958 3.5908E-12 0.0016 0.0000 10.0 78055
10240000 8.1212E-12 0.0029 0.0001 10.0 19977 3.5913E-12 0.0015
0.0000 10.0 7799612288000 8.1331E-12 0.0026 0.0001 10.0 19917
3.5882E-12 0.0013 0.0000 10.0 78360
8Calculation No. PM-1083, Rev. 0 Attachment C (Page C8 of
C9)
-
HIC.o 5/31/2009
14336000 8.1450E-12 0.0024 0.0001 10.0 19855 3.5882E-12 0.0012
0.0000 10.0 7826916384000 8.1327E-12 0.0023 0.0001 10.0 19849
3.5867E-12 0.0012 0.0000 10.0 7819518432000 8.1367E-12 0.0022
0.0001 10.0 19869 3.5871E-12 0.0011 0.0000 10.0 7820220480000
8.1428E-12 0.0020 0.0001 10.0 19890 3.5843E-12 0.0010 0.0000 10.0
7826022528000 8.1370E-12 0.0019 0.0001 10.0 19891 3.5848E-12 0.0010
0.0000 10.0 7813324576000 8.1375E-12 0.0019 0.0001 10.0 19850
3.5845E-12 0.0009 0.0000 10.0 7816525000000 8.1400E-12 0.0019
0.0001 10.0 19848 3.5841E-12 0.0009 0.0000 10.0 78191
***********************************************************************************************************************dump
no. 2 on file runtpg nps = 25000000 coll = 221983479 ctm = 14.69
nrn = 2964429671
9 warning messages so far.
run terminated when 25000000 particle histories were done.
computer time = 15.07 minutes
mcnp version 5 07-10-03 05/31/09 11:50:54 probid = 05/31/09
11:36:12
9Calculation No. PM-1083, Rev. 0 Attachment C (Page C9 of
C9)
-
PBrs1.o 5/30/2009
Thread Name & Version = mcnp5_RSICC, 4.23_
._ _ _ ._ ._ |_| | | (_ | | |_) _)
|+--------------------------------------------------------------------+|
This program was prepared by the Regents of the University of
||California at Los Alamos National Laboratory (the University)
under || contract number W-7405-ENG-36 with the U.S. Department of
Energy ||(DoE). The University has certain rights in the program
pursuant to|| the contract and the program should not be copied or
distributed || outside your organization. All rights in the program
are reserved ||by the DoE and the University. Neither the U.S.
Government nor the || University makes any warranty, express or
implied, or assumes any || liability or responsibility for the use
of this software.
|+--------------------------------------------------------------------+
1mcnp version 5 ld=07-10-03 05/30/09
19:35:46*************************************************************************
probid = 05/30/09 19:35:46
inp = PBrs1 outp = PBrs1.o
1- c Created on: Thursday, May 14, 2009 at 23:162- 1 228 -2.243
1 -4 5 -8 9 -10 $floor3- 2 228 -2.243 1 -2 5 -8 10 -19 $wall 14- 3
228 -2.243 3 -4 5 -8 10 -19 $wall 25- 4 228 -2.243 2 -3 5 -6 10 -19
$wall 36- 5 228 -2.243 2 -3 7 -8 10 -19 $wall 47- 6 280 -0.9 -21 11
-12 $L1-NE8- 7 280 -0.9 -22 11 -12 $L1-NW9- 8 280 -0.9 -23 11 -12
$L1-SE
10- 9 280 -0.9 -24 11 -12 $L1-SW11- 22 204 -0.00122 2 -3 6 -7 10
-19 #6 #7 #8 #912- 23 228 -2.243 1 -4 5 -8 19 -2513- 24 228 -2.243
1 -4 5 -8 25 -2614- 25 228 -2.243 1 -4 5 -8 26 -2715- 26 228 -2.243
1 -4 5 -8 27 -2816- 27 228 -2.243 1 -4 5 -8 28 -2917- 28 228 -2.243
1 -4 5 -8 29 -3018- 29 228 -2.243 1 -4 5 -8 30 -3119- 30 228 -2.243
1 -4 5 -8 31 -2020- 31 0 -1 :4 :-5 :8 :-9 :2021-22- *1 px
-243.8423- 2 px -226.0624- 3 px 226.0625- *4 px 243.8426- *5 py
-25427- 6 py -236.2228- 7 py 236.2229- *8 py 25430- 9 pz -7.6231-
10 pz 032- 11 pz 55.8833- 12 pz 238.7634- 13 pz 316.2335- 14 pz
499.1136- 15 pz 576.5837- 16 pz 759.46
1Calculation No. PM-1083, Rev. 0 Attachment D (Page D1 of
D75)
-
PBrs1.o 5/30/2009
38- 17 pz 836.9339- 18 pz 1019.8140- 19 pz 1097.2841- 20 pz
1158.2442- 21 c/z -113.03 -118.11 76.243- 22 c/z 113.03 -118.11
76.244- 23 c/z -113.03 118.11 76.245- 24 c/z 113.03 118.11 76.246-
25 pz 1104.947- 26 pz 1112.5248- 27 pz 1120.1449- 28 pz 1127.7650-
29 pz 1135.3851- 30 pz 114352- 31 pz 1150.6253-54- mode p55- m280
1001.60c -0.111915 $Water (density of 1 assumed)
warning. neutron table inconsistent with mode will be
ignored.56- 8016.60c -0.888085
warning. neutron table inconsistent with mode will be
ignored.57- m228 1001.60c -0.005558 $concrete (ordinary with
ENDF-VI)
warning. neutron table inconsistent with mode will be
ignored.58- 8016.60c -0.498076 11023.60c -0.017101 12000.60c
-0.002565
warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.
59- 13027.60c -0.045746 14000.60c -0.315092 16000.60c
-0.001283warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.
60- 19000.60c -0.019239 20000.60c -0.082941 26054.60c
-0.000707warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.
61- 26056.60c -0.01139 26057.60c -0.000265 26058.60c
-3.6e-005warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.
62- m204 7014.60c -0.755636 $air (US S. Atm at sea
level)warning. neutron table inconsistent with mode will be
ignored.
63- 8016.60c -0.231475 18000.59c -0.012889warning. neutron table
inconsistent with mode will be ignored.warning. neutron table
inconsistent with mode will be ignored.
64- imp:p 1 3 3r 1 3r $ 1, 965- 3 6 12 24 48 $ 22, 266- 96 192
384 768 0 $ 27, 367- phys:p 100 1 068- sdef erg=d5 x=d1 y=d2 z=d3
cel=d4 eff=0.0000169- si1 -190 19070- sp1 0.0 1.071- si2 -195
19572- sp2 0.0 1.073- si3 55 24074- sp3 0.0 1.075- si4 l 6 7 8 976-
sp4 1 1 1 177- si5 l 1.1730 1.3325
2Calculation No. PM-1083, Rev. 0 Attachment D (Page D2 of
D75)
-
PBrs1.o 5/30/2009
78- sp5 1 179- ssw 2080- f2:p 19 25 26 27 28 29 30 31 20
warning. without bremsstrahlung, flux estimates will be low.81-
fs2 -2 3 -6 7 -21 -22 -23 -2482- sd2 19329 19329 16077 16077 18241
18241 18241 18241 140636 &83- 19329 19329 16077 16077 18241
18241 18241 18241 140636 &84- 19329 19329 16077 16077 18241
18241 18241 18241 140636 &85- 19329 19329 16077 16077 18241
18241 18241 18241 140636 &86- 19329 19329 16077 16077 18241
18241 18241 18241 140636 &87- 19329 19329 16077 16077 18241
18241 18241 18241 140636 &88- 19329 19329 16077 16077 18241
18241 18241 18241 140636 &89- 19329 19329 16077 16077 18241
18241 18241 18241 140636 &90- 19329 19329 16077 16077 18241
18241 18241 18241 14063691- f22:p 19 25 26 27 28 29 30 31 20
warning. without bremsstrahlung, flux estimates will be low.92-
fs22 -2 3 -6 793- sd22 1.9329 1.9329 1.6077 1.6077 213600 &94-
1.9329 1.9329 1.6077 1.6077 213600 &95- 1.9329 1.9329 1.6077
1.6077 213600 &96- 1.9329 1.9329 1.6077 1.6077 213600 &97-
1.9329 1.9329 1.6077 1.6077 213600 &98- 1.9329 1.9329 1.6077
1.6077 213600 &99- 1.9329 1.9329 1.6077 1.6077 213600 &
100- 1.9329 1.9329 1.6077 1.6077 213600 &101- 1.9329 1.9329
1.6077 1.6077 213600102- de2 0.01 0.03 0.05 0.07 0.10 0.15 0.20
&103- 0.25 0.30 0.35 0.40 0.45 0.50 0.55 &104- 0.60 0.65
0.70 0.80 1.00 1.40 1.80 &105- 2.20 2.60 2.80 3.25 3.75 4.25
4.75 &106- 5.00 5.25 5.75 6.25 6.75 7.50 9.00 &107- 11.0
13.0 15.0108- df2 LIN 3.96E-06 5.82E-07 2.90E-07 2.58E-07 2.83E-07
3.79E-07 5.01E-07 &109- 6.31e-07 7.59E-07 8.78E-07 9.85E-07
1.08E-06 1.17E-06 1.27E-06 &110- 1.36e-06 1.44E-06 1.52E-06
1.68E-06 1.98E-06 2.51E-06 2.99E-06 &111- 3.42e-06 3.82E-06
4.01E-06 4.41E-06 4.83E-06 5.23E-06 5.60E-06 &112- 5.80e-06
6.01E-06 6.37E-06 6.74E-06 7.11E-06 7.66E-06 8.77E-06 &113-
1.03e-05 1.18E-05 1.33E-05114- de22 0.01 0.03 0.05 0.07 0.10 0.15
0.20 &115- 0.25 0.30 0.35 0.40 0.45 0.50 0.55 &116- 0.60
0.65 0.70 0.80 1.00 1.40 1.80 &117- 2.20 2.60 2.80 3.25 3.75
4.25 4.75 &118- 5.00 5.25 5.75 6.25 6.75 7.50 9.00 &119-
11.0 13.0 15.0120- df22 LIN 3.96E-06 5.82E-07 2.90E-07 2.58E-07
2.83E-07 3.79E-07 5.01E-07 &121- 6.31e-07 7.59E-07 8.78E-07
9.85E-07 1.08E-06 1.17E-06 1.27E-06 &122- 1.36e-06 1.44E-06
1.52E-06 1.68E-06 1.98E-06 2.51E-06 2.99E-06 &123- 3.42e-06
3.82E-06 4.01E-06 4.41E-06 4.83E-06 5.23E-06 5.60E-06 &124-
5.80e-06 6.01E-06 6.37E-06 6.74E-06 7.11E-06 7.66E-06 8.77E-06
&125- 1.03e-05 1.18E-05 1.33E-05126- nps 25000000
warning. sum of segment sizes differs from total in 18
cases.1cells print table 60
atom gram photoncell mat density density volume mass pieces
importance
3Calculation No. PM-1083, Rev. 0 Attachment D (Page D3 of
D75)
-
PBrs1.o 5/30/2009
1 1 228 7.19163E-02 2.24300E+00 1.88779E+06 4.23431E+06 0
1.0000E+002 2 228 7.19163E-02 2.24300E+00 9.91090E+06 2.22301E+07 0
3.0000E+003 3 228 7.19163E-02 2.24300E+00 9.91090E+06 2.22301E+07 0
3.0000E+004 4 228 7.19163E-02 2.24300E+00 8.82070E+06 1.97848E+07 0
3.0000E+005 5 228 7.19163E-02 2.24300E+00 8.82070E+06 1.97848E+07 0
3.0000E+006 6 280 9.02777E-02 9.00000E-01 3.33600E+06 3.00240E+06 1
1.0000E+007 7 280 9.02777E-02 9.00000E-01 3.33600E+06 3.00240E+06 1
1.0000E+008 8 280 9.02777E-02 9.00000E-01 3.33600E+06 3.00240E+06 1
1.0000E+009 9 280 9.02777E-02 9.00000E-01 3.33600E+06 3.00240E+06 1
1.0000E+00
10 22 204 5.05147E-05 1.22000E-03 0.00000E+00 0.00000E+00 0
3.0000E+0011 23 228 7.19163E-02 2.24300E+00 1.88779E+06 4.23431E+06
0 6.0000E+0012 24 228 7.19163E-02 2.24300E+00 1.88779E+06
4.23431E+06 0 1.2000E+0113 25 228 7.19163E-02 2.24300E+00
1.88779E+06 4.23431E+06 0 2.4000E+0114 26 228 7.19163E-02
2.24300E+00 1.88779E+06 4.23431E+06 0 4.8000E+0115 27 228
7.19163E-02 2.24300E+00 1.88779E+06 4.23431E+06 0 9.6000E+0116 28
228 7.19163E-02 2.24300E+00 1.88779E+06 4.23431E+06 0 1.9200E+0217
29 228 7.19163E-02 2.24300E+00 1.88779E+06 4.23431E+06 0
3.8400E+0218 30 228 7.19163E-02 2.24300E+00 1.88779E+06 4.23431E+06
0 7.6800E+0219 31 0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0 0.0000E+00
total 6.77973E+07 1.34148E+08
warning. surface 13 is not used for anything.
warning. surface 14 is not used for anything.
warning. surface 15 is not used for anything.
warning. surface 16 is not used for anything.
warning. surface 17 is not used for anything.
warning. surface 18 is not used for anything.
minimum source weight = 1.0000E+00 maximum source weight =
1.0000E+00
**************************************************** Random
Number Generator = 1 ** Random Number Seed = 19073486328125 **
Random Number Multiplier = 19073486328125 ** Random Number Adder =
0 ** Random Number Bits Used = 48 ** Random Number Stride = 152917
****************************************************
27 warning messages so far.1cross-section tables print table
100
table length
tables from file mcplib04
1000.04p 1898 ENDF/B-VI Release 8 Photoatomic Data for 1-H mat
100 02/07/037000.04p 3194 ENDF/B-VI Release 8 Photoatomic Data for
7-N mat 700 02/07/038000.04p 3272 ENDF/B-VI Release 8 Photoatomic
Data for 8-O mat 800 02/07/03
11000.04p 3995 ENDF/B-VI Release 8 Photoatomic Data for 11-NA
mat1100 02/07/0312000.04p 3781 ENDF/B-VI Release 8 Photoatomic Data
for 12-MG mat1200 02/07/03
4Calculation No. PM-1083, Rev. 0 Attachment D (Page D4 of
D75)
-
PBrs1.o 5/30/2009
13000.04p 4846 ENDF/B-VI Release 8 Photoatomic Data for 13-AL
mat1300 02/07/0314000.04p 4792 ENDF/B-VI Release 8 Photoatomic Data
for 14-SI mat1400 02/07/0316000.04p 4654 ENDF/B-VI Release 8
Photoatomic Data for 16-S mat1600 02/07/0318000.04p 4696 ENDF/B-VI
Release 8 Photoatomic Data for 18-AR mat1800 02/07/0319000.04p 5047
ENDF/B-VI Release 8 Photoatomic Data for 19-K mat1900
02/07/0320000.04p 5013 ENDF/B-VI Release 8 Photoatomic Data for
20-CA mat2000 02/07/0326000.04p 5718 ENDF/B-VI Release 8
Photoatomic Data for 26-FE mat2600 02/07/03
total 50906
warning. simple physics turned on for photons > 100 mev.
***********************************************************************************************************************dump
no. 1 on file runtpf nps = 0 coll = 0 ctm = 0.00 nrn = 0
28 warning messages so far.1problem summary
run terminated when 25000000 particle histories were done.+
05/30/09 20:33:22
c Created on: Thursday, May 14, 2009 at 23:16 probid = 05/30/09
19:35:460
photon creation tracks weight energy photon loss tracks weight
energy(per source particle) (per source particle)
source 25000000 1.0000E+00 1.2527E+00 escape 1323013 4.3444E-02
2.5498E-02energy cutoff 0 0. 0.time cutoff 0 0. 0.
weight window 0 0. 0. weight window 0 0. 0.cell importance
35195667 1.8729E-01 8.3630E-02 cell importance 15081919 1.8730E-01
8.3566E-02weight cutoff 0 0. 0. weight cutoff 0 0. 0.e or t
importance 0 0. 0. e or t importance 0 0. 0.dxtran 0 0. 0. dxtran 0
0. 0.forced collisions 0 0. 0. forced collisions 0 0. 0.exp.
transform 0 0. 0. exp. transform 0 0. 0.from neutrons 0 0. 0.
compton scatter 0 0. 1.1699E+00bremsstrahlung 0 0. 0. capture
47057545 9.9971E-01 5.7443E-02p-annihilation 29392 8.3076E-04
4.2453E-04 pair production 14696 4.1538E-04 5.4355E-04photonuclear
0 0. 0. photonuclear abs 0 0. 0.electron x-rays 0 0. 0.1st
fluorescence 3252114 4.2746E-02 1.7592E-042nd fluorescence 0 0.
0.
total 63477173 1.2309E+00 1.3370E+00 total 63477173 1.2309E+00
1.3370E+00
number of photons banked 19897799 average time of (shakes)
cutoffsphoton tracks per source particle 2.5391E+00 escape
5.8895E-01 tco 1.0000E+33photon collisions per source particle
1.8901E+01 capture 7.4525E-01 eco 1.0000E-03total photon collisions
472523354 capture or escape 7.3875E-01 wc1 -5.0000E-01
any termination 7.6180E-01 wc2 -2.5000E-01
computer time so far in this run 57.50 minutes maximum number
ever in bank 11computer time in mcrun 57.18 minutes bank overflows
to backup file 0source particles per minute 4.3718E+05random
numbers generated 6518761511 most random numbers used was 2542674
in history 20649758
warning. random number stride 152917 exceeded 1 times.
5Calculation No. PM-1083, Rev. 0 Attachment D (Page D5 of
D75)
-
PBrs1.o 5/30/2009
range of sampled source weights = 1.0000E+00 to 1.0000E+00
source efficiency = 0.1216 in cell 6
source efficiency = 0.1217 in cell 7
source efficiency = 0.1217 in cell 8
source efficiency = 0.1217 in cell 91photon activity in each
cell print table 126
tracks population collisions collisions number flux average
averagecell entering * weight weighted weighted track weight track
mfp
(per history) energy energy (relative) (cm)
1 1 3309370 3276959 10096714 4.0387E-01 4.4288E-01 4.4288E-01
1.0000E+00 4.3112E+002 2 9076429 8789954 47900528 6.3867E-01
3.9482E-01 3.9482E-01 1.0000E+00 4.0902E+003 3 9075758 8788977
47888930 6.3852E-01 3.9477E-01 3.9477E-01 1.0000E+00 4.0900E+004 4
8611031 8323470 44702937 5.9604E-01 3.9462E-01 3.9462E-01
1.0000E+00 4.0885E+005 5 8608148 8321321 44729490 5.9639E-01
3.9457E-01 3.9457E-01 1.0000E+00 4.0884E+006 6 7153895 7089022
61173452 2.4469E+00 4.5520E-01 4.5520E-01 1.0000E+00 9.9689E+007 7
7150336 7084341 61149505 2.4460E+00 4.5484E-01 4.5484E-01
1.0000E+00 9.9651E+008 8 7149055 7083720 61144583 2.4458E+00
4.5502E-01 4.5502E-01 1.0000E+00 9.9669E+009 9 7153061 7087612
61123219 2.4449E+00 4.5524E-01 4.5524E-01 1.0000E+00 9.9695E+00
10 22 56395741 43170666 1594932 2.1266E-02 4.5286E-01 4.5286E-01
1.0000E+00 8.2606E+0311 23 2254062 2128013 6624199 4.4161E-02
3.2722E-01 3.2722E-01 1.0000E+00 3.7313E+0012 24 1727857 1575707
5293735 1.7646E-02 3.5374E-01 3.5374E-01 1.0000E+00 3.9233E+0013 25
1454357 1321613 4477715 7.4629E-03 3.6395E-01 3.6395E-01 1.0000E+00
3.9914E+0014 26 1258644 1142110 3883679 3.2364E-03 3.6512E-01
3.6512E-01 1.0000E+00 4.0053E+0015 27 1091639 990753 3386252
1.4109E-03 3.6272E-01 3.6272E-01 1.0000E+00 3.9988E+0016 28 943792
856543 2933174 6.1108E-04 3.5902E-01 3.5902E-01 1.0000E+00
3.9838E+0017 29 806041 732327 2511922 2.6166E-04 3.5522E-01
3.5522E-01 1.0000E+00 3.9694E+0018 30 612714 587278 1908388
9.9395E-05 3.7056E-01 3.7056E-01 1.0000E+00 4.0633E+00
total 133831930 118350386 472523354 1.2753E+01
1tally 2 nps = 25000000tally type 2 particle flux averaged over
a surface.tally for photons
this tally is modified by a dose function.
areassurface: 19 25 26 27 28 29 30
segment1 1.93290E+04 1.93290E+04 1.93290E+04 1.93290E+04
1.93290E+04 1.93290E+04 1.93290E+042 1.93290E+04 1.93290E+04
1.93290E+04 1.93290E+04 1.93290E+04 1.93290E+04 1.93290E+043
1.60770E+04 1.60770E+04 1.60770E+04 1.60770E+04 1.60770E+04
1.60770E+04 1.60770E+044 1.60770E+04 1.60770E+04 1.60770E+04
1.60770E+04 1.60770E+04 1.60770E+04 1.60770E+045 1.82410E+04
1.82410E+04 1.82410E+04 1.82410E+04 1.82410E+04 1.82410E+04
1.82410E+046 1.82410E+04 1.82410E+04 1.82410E+04 1.82410E+04
1.82410E+04 1.82410E+04 1.82410E+047 1.82410E+04 1.82410E+04
1.82410E+04 1.82410E+04 1.82410E+04 1.82410E+04 1.82410E+048
1.82410E+04 1.82410E+04 1.82410E+04 1.82410E+04 1.82410E+04
1.82410E+04 1.82410E+049 1.40636E+05 1.40636E+05 1.40636E+05
1.40636E+05 1.40636E+05 1.40636E+05 1.40636E+05
surface: 31 20segment
1 1.93290E+04 1.93290E+042 1.93290E+04 1.93290E+04
6Calculation No. PM-1083, Rev. 0 Attachment D (Page D6 of
D75)
-
PBrs1.o 5/30/2009
3 1.60770E+04 1.60770E+044 1.60770E+04 1.60770E+045 1.82410E+04
1.82410E+046 1.82410E+04 1.82410E+047 1.82410E+04 1.82410E+048
1.82410E+04 1.82410E+049 1.40636E+05 1.40636E+05
surface 19segment: -2
4.93966E-15 0.0244
surface 19segment: 2 3
4.92641E-15 0.0240
surface 19segment: 2 -3 -6
5.86835E-15 0.0248
surface 19segment: 2 -3 6 7
5.89119E-15 0.0244
surface 19segment: 2 -3 6 -7 -21
7.66918E-14 0.0066
surface 19segment: 2 -3 6 -7 21 -22
7.74915E-14 0.0065
surface 19segment: 2 -3 6 -7 21 22 -23
7.75238E-14 0.0066
surface 19segment: 2 -3 6 -7 21 22 23 -24
7.73216E-14 0.0065
surface 19segment: 2 -3 6 -7 21 22 23 24
7.70793E-14 0.0024
surface 25segment: -2
3.90868E-15 0.0241
surface 25segment: 2 3
4.11788E-15 0.0237
surface 25segment: 2 -3 -6
4.79158E-15 0.0242
surface 25segment: 2 -3 6 7
7Calculation No. PM-1083, Rev. 0 Attachment D (Page D7 of
D75)
-
PBrs1.o 5/30/2009
4.67300E-15 0.0241
surface 25segment: 2 -3 6 -7 -21
3.78231E-14 0.0090
surface 25segment: 2 -3 6 -7 21 -22
3.84556E-14 0.0090
surface 25segment: 2 -3 6 -7 21 22 -23
3.82674E-14 0.0090
surface 25segment: 2 -3 6 -7 21 22 23 -24
3.84850E-14 0.0090
surface 25segment: 2 -3 6 -7 21 22 23 24
3.78882E-14 0.0033
surface 26segment: -2
2.31688E-15 0.0244
surface 26segment: 2 3
2.47886E-15 0.0242
surface 26segment: 2 -3 -6
2.87603E-15 0.0247
surface 26segment: 2 -3 6 7
2.89673E-15 0.0250
surface 26segment: 2 -3 6 -7 -21
1.66522E-14 0.0109
surface 26segment: 2 -3 6 -7 21 -22
1.69442E-14 0.0108
surface 26segment: 2 -3 6 -7 21 22 -23
1.67776E-14 0.0108
surface 26segment: 2 -3 6 -7 21 22 23 -24
1.69505E-14 0.0109
surface 26segment: 2 -3 6 -7 21 22 23 24
1.67110E-14 0.0040
8Calculation No. PM-1083, Rev. 0 Attachment D (Page D8 of
D75)
-
PBrs1.o 5/30/2009
surface 27segment: -2
1.27489E-15 0.0254
surface 27segment: 2 3
1.35356E-15 0.0253
surface 27segment: 2 -3 -6
1.59001E-15 0.0259
surface 27segment: 2 -3 6 7
1.51149E-15 0.0259
surface 27segment: 2 -3 6 -7 -21
7.28563E-15 0.0126
surface 27segment: 2 -3 6 -7 21 -22
7.43915E-15 0.0125
surface 27segment: 2 -3 6 -7 21 22 -23
7.33610E-15 0.0125
surface 27segment: 2 -3 6 -7 21 22 23 -24
7.36363E-15 0.0126
surface 27segment: 2 -3 6 -7 21 22 23 24
7.26979E-15 0.0046
surface 28segment: -2
6.72562E-16 0.0260
surface 28segment: 2 3
6.93648E-16 0.0256
surface 28segment: 2 -3 -6
8.04316E-16 0.0276
surface 28segment: 2 -3 6 7
7.58138E-16 0.0271
surface 28segment: 2 -3 6 -7 -21
3.13258E-15 0.0143
surface 28segment: 2 -3 6 -7 21 -22
9Calculation No. PM-1083, Rev. 0 Attachment D (Page D9 of
D75)
-
PBrs1.o 5/30/2009
3.23005E-15 0.0141
surface 28segment: 2 -3 6 -7 21 22 -23
3.19447E-15 0.0142
surface 28segment: 2 -3 6 -7 21 22 23 -24
3.16036E-15 0.0141
surface 28segment: 2 -3 6 -7 21 22 23 24
3.13435E-15 0.0052
surface 29segment: -2
3.19758E-16 0.0270
surface 29segment: 2 3
3.35201E-16 0.0274
surface 29segment: 2 -3 -6
3.84221E-16 0.0283
surface 29segment: 2 -3 6 7
3.65135E-16 0.0277
surface 29segment: 2 -3 6 -7 -21
1.32916E-15 0.0158
surface 29segment: 2 -3 6 -7 21 -22
1.38544E-15 0.0156
surface 29segment: 2 -3 6 -7 21 22 -23
1.36407E-15 0.0157
surface 29segment: 2 -3 6 -7 21 22 23 -24
1.33969E-15 0.0159
surface 29segment: 2 -3 6 -7 21 22 23 24
1.33724E-15 0.0058
surface 30segment: -2
1.51925E-16 0.0285
surface 30segment: 2 3
1.61516E-16 0.0285
10Calculation No. PM-1083, Rev. 0 Attachment D (Page D10 of
D75)
-
PBrs1.o 5/30/2009
surface 30segment: 2 -3 -6
1.81383E-16 0.0296
surface 30segment: 2 -3 6 7
1.74717E-16 0.0296
surface 30segment: 2 -3 6 -7 -21
5.65637E-16 0.0176
surface 30segment: 2 -3 6 -7 21 -22
5.80514E-16 0.0173
surface 30segment: 2 -3 6 -7 21 22 -23
5.80783E-16 0.0172
surface 30segment: 2 -3 6 -7 21 22 23 -24
5.62206E-16 0.0176
surface 30segment: 2 -3 6 -7 21 22 23 24
5.65081E-16 0.0064
surface 31segment: -2
6.74962E-17 0.0304
surface 31segment: 2 3
7.17668E-17 0.0303
surface 31segment: 2 -3 -6
8.40617E-17 0.0314
surface 31segment: 2 -3 6 7
7.78816E-17 0.0314
surface 31segment: 2 -3 6 -7 -21
2.35765E-16 0.0195
surface 31segment: 2 -3 6 -7 21 -22
2.44906E-16 0.0190
surface 31segment: 2 -3 6 -7 21 22 -23
2.43233E-16 0.0190
surface 31segment: 2 -3 6 -7 21 22 23 -24
11Calculation No. PM-1083, Rev. 0 Attachment D (Page D11 of
D75)
-
PBrs1.o 5/30/2009
2.32883E-16 0.0196
surface 31segment: 2 -3 6 -7 21 22 23 24
2.36451E-16 0.0071
surface 20segment: -2
2.39165E-17 0.0332
surface 20segment: 2 3
2.57492E-17 0.0332
surface 20segment: 2 -3 -6
3.05480E-17 0.0347
surface 20segment: 2 -3 6 7
2.79797E-17 0.0336
surface 20segment: 2 -3 6 -7 -21
7.99483E-17 0.0217
surface 20segment: 2 -3 6 -7 21 -22
8.32255E-17 0.0211
surface 20segment: 2 -3 6 -7 21 22 -23
8.27956E-17 0.0210
surface 20segment: 2 -3 6 -7 21 22 23 -24
7.94216E-17 0.0218
surface 20segment: 2 -3 6 -7 21 22 23 24
8.03972E-17 0.0079
===================================================================================================================================
results of 10 statistical checks for the estimated answer for
the tally fluctuation chart (tfc) bin of tally 2
tfc bin --mean-- ---------relative error--------- ----variance
of the variance---- --figure of merit-- -pdf-behavior behavior
value decrease decrease rate value decrease decrease rate value
behavior slope
desired random
-
PBrs1.o 5/30/2009
the results in other bins associated with this tally may not
meet these statistical criteria.
----- estimated confidence intervals: -----
estimated asymmetric confidence interval(1,2,3 sigma):
7.6894E-14 to 7.7265E-14; 7.6709E-14 to 7.7450E-14; 7.6524E-14 to
7.7636E-14estimated symmetric confidence interval(1,2,3 sigma):
7.6894E-14 to 7.7265E-14; 7.6709E-14 to 7.7450E-14; 7.6523E-14 to
7.7635E-14
1analysis of the results in the tally fluctuation chart bin
(tfc) for tally 2 with nps = 25000000 print table 160
normed average tally per history = 7.70793E-14 unnormed average
tally per history = 1.08401E-08estimated tally relative error =
0.0024 estimated variance of the variance = 0.0000relative error
from zero tallies = 0.0016 relative error from nonzero scores =
0.0018
number of nonzero history tallies = 408672 efficiency for the
nonzero tallies = 0.0163history number of largest tally = 4536888
largest unnormalized history tally = 1.28939E-05(largest
tally)/(average tally) = 1.18946E+03 (largest tally)/(avg nonzero
tally)= 1.94440E+01
(confidence interval shift)/mean = 0.0000 shifted confidence
interval center = 7.70796E-14
if the largest history score sampled so far were to occur on the
next history, the tfc bin quantities would change as follows:
estimated quantities value at nps value at nps+1
value(nps+1)/value(nps)-1.
mean 7.70793E-14 7.70830E-14 0.000048relative error 2.40399E-03
2.40435E-03 0.000148variance of the variance 1.51432E-05
1.52842E-05 0.009306shifted center 7.70796E-14 7.70796E-14
0.000000figure of merit 3.02589E+03 3.02500E+03 -0.000296
the estimated inverse power slope of the 198 largest tallies
starting at 5.43839E-06 is 9.5557the history score probability
density function appears to have an unsampled region at the largest
history scores:please examine. see print table 161.
fom = (histories/minute)*(f(x) signal-to-noise ratio)**2 =
(4.372E+05)*( 8.319E-02)**2 = (4.372E+05)*(6.921E-03) =
3.026E+03
1tally 22 nps = 25000000tally type 2 particle flux averaged over
a surface.tally for photons
this tally is modified by a dose function.
areassurface: 19 25 26 27 28 29 30
segment1 1.93290E+00 1.93290E+00 1.93290E+00 1.93290E+00
1.93290E+00 1.93290E+00 1.93290E+002 1.93290E+00 1.93290E+00
1.93290E+00 1.93290E+00 1.93290E+00 1.93290E+00 1.93290E+003
1.60770E+00 1.60770E+00 1.60770E+00 1.60770E+00 1.60770E+00
1.60770E+00 1.60770E+004 1.60770E+00 1.60770E+00 1.60770E+00
1.60770E+00 1.60770E+00 1.60770E+00 1.60770E+005 2.13600E+05
2.13600E+05 2.13600E+05 2.13600E+05 2.13600E+05 2.13600E+05
2.13600E+05
surface: 31 20segment
1 1.93290E+00 1.93290E+002 1.93290E+00 1.93290E+003 1.60770E+00
1.60770E+004 1.60770E+00 1.60770E+005 2.13600E+05 2.13600E+05
13Calculation No. PM-1083, Rev. 0 Attachment D (Page D13 of
D75)
-
PBrs1.o 5/30/2009
surface 19segment: -2
4.93966E-11 0.0244
surface 19segment: 2 3
4.92641E-11 0.0240
surface 19segment: 2 -3 -6
5.86835E-11 0.0248
surface 19segment: 2 -3 6 7
5.89119E-11 0.0244
surface 19segment: 2 -3 6 -7
7.71401E-14 0.0020
surface 25segment: -2
3.90868E-11 0.0241
surface 25segment: 2 3
4.11788E-11 0.0237
surface 25segment: 2 -3 -6
4.79158E-11 0.0242
surface 25segment: 2 -3 6 7
4.67300E-11 0.0241
surface 25segment: 2 -3 6 -7
3.80145E-14 0.0027
surface 26segment: -2
2.31688E-11 0.0244
surface 26segment: 2 3
2.47886E-11 0.0242
surface 26segment: 2 -3 -6
2.87603E-11 0.0247
surface 26segment: 2 -3 6 7
2.89673E-11 0.0250
surface 26
14Calculation No. PM-1083, Rev. 0 Attachment D (Page D14 of
D75)
-
PBrs1.o 5/30/2009
segment: 2 -3 6 -71.67520E-14 0.0033
surface 27segment: -2
1.27489E-11 0.0254
surface 27segment: 2 3
1.35356E-11 0.0253
surface 27segment: 2 -3 -6
1.59001E-11 0.0259
surface 27segment: 2 -3 6 7
1.51149E-11 0.0259
surface 27segment: 2 -3 6 -7
7.29928E-15 0.0038
surface 28segment: -2
6.72562E-12 0.0260
surface 28segment: 2 3
6.93648E-12 0.0256
surface 28segment: 2 -3 -6
8.04316E-12 0.0276
surface 28segment: 2 -3 6 7
7.58138E-12 0.0271
surface 28segment: 2 -3 6 -7
3.14973E-15 0.0043
surface 29segment: -2
3.19758E-12 0.0270
surface 29segment: 2 3
3.35201E-12 0.0274
surface 29segment: 2 -3 -6
3.84221E-12 0.0283
surface 29segment: 2 -3 6 7
3.65135E-12 0.0277
15Calculation No. PM-1083, Rev. 0 Attachment D (Page D15 of
D75)
-
PBrs1.o 5/30/2009
surface 29segment: 2 -3 6 -7
1.34316E-15 0.0048
surface 30segment: -2
1.51925E-12 0.0285
surface 30segment: 2 3
1.61516E-12 0.0285
surface 30segment: 2 -3 -6
1.81383E-12 0.0296
surface 30segment: 2 -3 6 7
1.74717E-12 0.0296
surface 30segment: 2 -3 6 -7
5.67542E-16 0.0053
surface 31segment: -2
6.74962E-13 0.0304
surface 31segment: 2 3
7.17668E-13 0.0303
surface 31segment: 2 -3 -6
8.40617E-13 0.0314
surface 31segment: 2 -3 6 7
7.78816E-13 0.0314
surface 31segment: 2 -3 6 -7
2.37389E-16 0.0059
surface 20segment: -2
2.39165E-13 0.0332
surface 20segment: 2 3
2.57492E-13 0.0332
surface 20segment: 2 -3 -6
3.05480E-13 0.0347
surface 20
16Calculation No. PM-1083, Rev. 0 Attachment D (Page D16 of
D75)
-
PBrs1.o 5/30/2009
segment: 2 -3 6 72.79797E-13 0.0336
surface 20segment: 2 -3 6 -7
8.07219E-17 0.0066
===================================================================================================================================
results of 10 statistical checks for the estimated answer for
the tally fluctuation chart (tfc) bin of tally 22
tfc bin --mean-- ---------relative error--------- ----variance
of the variance---- --figure of merit-- -pdf-behavior behavior
value decrease decrease rate value decrease decrease rate value
behavior slope
desired random
-
PBrs1.o 5/30/2009
fom = (histories/minute)*(f(x) signal-to-noise ratio)**2 =
(4.372E+05)*( 1.017E-01)**2 = (4.372E+05)*(1.034E-02) =
4.521E+03
1status of the statistical checks used to form confidence
intervals for the mean for each tally bin
tally result of statistical checks for the tfc bin (the first
check not passed is listed) and error magnitude check for all
bins
2 passed the 10 statistical checks for the tally fluctuation
chart bin resultpassed all bin error check: 81 tally bins all have
relative errors less than 0.10 with no zero bins
22 passed the 10 statistical checks for the tally fluctuation
chart bin resultpassed all bin error check: 45 tally bins all have
relative errors less than 0.10 with no zero bins
the 10 statistical checks are only for the tally fluctuation
chart bin and do not apply to other tally bins.
1tally fluctuation charts
tally 2 tally 22nps mean error vov slope fom mean error vov
slope fom
2048000 7.7203E-14 0.0084 0.0002 5.3 3015 7.7090E-14 0.0069
0.0001 10.0 45094096000 7.7283E-14 0.0059 0.0001 10.0 3030
7.7090E-14 0.0049 0.0001 10.0 45276144000 7.6807E-14 0.0049 0.0001
10.0 3019 7.6839E-14 0.0040 0.0000 10.0 45228192000 7.7031E-14
0.0042 0.0000 10.0 3023 7.7074E-14 0.0034 0.0000 10.0 4523
10240000 7.6945E-14 0.0038 0.0000 10.0 3020 7.7028E-14 0.0031
0.0000 6.1 452412288000 7.7001E-14 0.0034 0.0000 10.0 3018
7.7040E-14 0.0028 0.0000 4.6 451814336000 7.6905E-14 0.0032 0.0000
10.0 3020 7.6967E-14 0.0026 0.0000 5.0 451816384000 7.7057E-14
0.0030 0.0000 10.0 3022 7.7002E-14 0.0024 0.0000 4.8 451618432000
7.7157E-14 0.0028 0.0000 10.0 3025 7.7109E-14 0.0023 0.0000 5.0
451920480000 7.7184E-14 0.0027 0.0000 10.0 3027 7.7162E-14 0.0022
0.0000 4.8 452222528000 7.7160E-14 0.0025 0.0000 8.6 3026
7.7134E-14 0.0021 0.0000 5.2 452024576000 7.7087E-14 0.0024 0.0000
9.4 3026 7.7133E-14 0.0020 0.0000 4.7 452125000000 7.7079E-14
0.0024 0.0000 9.6 3026 7.7140E-14 0.0020 0.0000 4.9 4521
***********************************************************************************************************************dump
no. 2 on file runtpf nps = 25000000 coll = 472523354 ctm = 57.18
nrn = 6518761511
surface-source file wssb with nps = 25000000 and 237227 tracks
was written forsurface(s) 20
29 warning messages so far.
run terminated when 25000000 particle histories were done.
computer time = 57.52 minutes
mcnp version 5 07-10-03 05/30/09 20:33:36 probid = 05/30/09
19:35:46
18Calculation No. PM-1083, Rev. 0 Attachment D (Page D18 of
D75)
-
PBrs2.o 5/30/2009
Thread Name & Version = mcnp5_RSICC, 4.23_
._ _ _ ._ ._ |_| | | (_ | | |_) _)
|+--------------------------------------------------------------------+|
This program was prepared by the Regents of the University of
||California at Los Alamos National Laboratory (the University)
under || contract number W-7405-ENG-36 with the U.S. Department of
Energy ||(DoE). The University has certain rights in the program
pursuant to|| the contract and the program should not be copied or
distributed || outside your organization. All rights in the program
are reserved ||by the DoE and the University. Neither the U.S.
Government nor the || University makes any warranty, express or
implied, or assumes any || liability or responsibility for the use
of this software.
|+--------------------------------------------------------------------+
1mcnp version 5 ld=07-10-03 05/30/09
17:21:01*************************************************************************
probid = 05/30/09 17:21:01
inp = PBrs2 outp = PBrs2.o
1- c Created on: Thursday, May 14, 2009 at 23:162- 1 228 -2.243
1 -4 5 -8 9 -10 $floor3- 2 228 -2.243 1 -2 5 -8 10 -19 $wall 14- 3
228 -2.243 3 -4 5 -8 10 -19 $wall 25- 4 228 -2.243 2 -3 5 -6 10 -19
$wall 36- 5 228 -2.243 2 -3 7 -8 10 -19 $wall 47- 6 280 -0.9 -21 11
-12 $L1-NE8- 7 280 -0.9 -22 11 -12 $L1-NW9- 8 280 -0.9 -23 11 -12
$L1-SE
10- 9 280 -0.9 -24 11 -12 $L1-SW11- 10 280 -0.9 -21 13 -14
$L2-NE12- 11 280 -0.9 -22 13 -14 $L2-NW13- 12 280 -0.9 -23 13 -14
$L2-SE14- 13 280 -0.9 -24 13 -14 $L2-SW15- 22 204 -0.00122 2 -3 6
-7 10 -19 #6 #7 #8 #9 #10 #11 #12 #1316- 23 228 -2.243 1 -4 5 -8 19
-2517- 24 228 -2.243 1 -4 5 -8 25 -2618- 25 228 -2.243 1 -4 5 -8 26
-2719- 26 228 -2.243 1 -4 5 -8 27 -2820- 27 228 -2.243 1 -4 5 -8 28
-2921- 28 228 -2.243 1 -4 5 -8 29 -3022- 29 228 -2.243 1 -4 5 -8 30
-3123- 30 228 -2.243 1 -4 5 -8 31 -2024- 31 0 -1 :4 :-5 :8 :-9
:2025-26- *1 px -243.8427- 2 px -226.0628- 3 px 226.0629- *4 px
243.8430- *5 py -25431- 6 py -236.2232- 7 py 236.2233- *8 py 25434-
9 pz -7.6235- 10 pz 036- 11 pz 55.8837- 12 pz 238.76
1Calculation No. PM-1083, Rev. 0 Attachment D (Page D19 of
D75)
-
PBrs2.o 5/30/2009
38- 13 pz 316.2339- 14 pz 499.1140- 15 pz 576.5841- 16 pz
759.4642- 17 pz 836.9343- 18 pz 1019.8144- 19 pz 1097.2845- 20 pz
1158.2446- 21 c/z -113.03 -118.11 76.247- 22 c/z 113.03 -118.11
76.248- 23 c/z -113.03 118.11 76.249- 24 c/z 113.03 118.11 76.250-
25 pz 1104.951- 26 pz 1112.5252- 27 pz 1120.1453- 28 pz 1127.7654-
29 pz 1135.3855- 30 pz 114356- 31 pz 1150.6257-58- mode p59- m280
1001.60c -0.111915 $Water (density of 1 assumed)
warning. neutron table inconsistent with mode will be
ignored.60- 8016.60c -0.888085
warning. neutron table inconsistent with mode will be
ignored.61- m228 1001.60c -0.005558 $concrete (ordinary with
ENDF-VI)
warning. neutron table inconsistent with mode will be
ignored.62- 8016.60c -0.498076 11023.60c -0.017101 12000.60c
-0.002565
warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.
63- 13027.60c -0.045746 14000.60c -0.315092 16000.60c
-0.001283warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.
64- 19000.60c -0.019239 20000.60c -0.082941 26054.60c
-0.000707warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.
65- 26056.60c -0.01139 26057.60c -0.000265 26058.60c
-3.6e-005warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.warning. neutron table inconsistent with mode will be
ignored.
66- m204 7014.60c -0.755636 $air (US S. Atm at sea
level)warning. neutron table inconsistent with mode will be
ignored.
67- 8016.60c -0.231475 18000.59c -0.012889warning. neutron table
inconsistent with mode will be ignored.warning. neutron table
inconsistent with mode will be ignored.
68- imp:p 1 3 3r 1 7r $ 1, 1369- 3 6 12 24 48 $ 22, 270- 96 192
384 768 0 $ 27, 371- phys:p 100 1 072- sdef erg=d5 x=d1 y=d2 z=d3
cel=d4 eff=0.0000173- si1 -190 19074- sp1 0.0 1.075- si2 -195
19576- sp2 0.0 1.077- si3 55 500
2Calculation No. PM-1083, Rev. 0 Attachment D (Page D20 of
D75)
-
PBrs2.o 5/30/2009
78- sp3 0.0 1.079- si4 l 6 7 8 9 10 11 12 1380- sp4 1 1 1 1 1 1
1 181- si5 l 1.1730 1.332582- sp5 1 183- ssw 2084- f2:p 19 25 26 27
28 29 30 31 20
warning. without bremsstrahlung, flux estimates will be low.85-
fs2 -2 3 -6 7 -21 -22 -23 -2486- sd2 19329 19329 16077 16077 18241
18241 18241 18241 140636 &87- 19329 19329 16077 16077 18241
18241 18241 18241 140636 &88- 19329 19329 16077 16077 18241
18241 18241 18241 140636 &89- 19329 19329 16077 16077 18241
18241 18241 18241 140636 &90- 19329 19329 16077 16077 18241
18241 18241 18241 140636 &91- 19329 19329 16077 16077 18241
18241 18241 18241 140636 &92- 19329 19329 16077 16077 18241
18241 18241 18241 140636 &93- 19329 19329 16077 16077 18241
18241 18241 18241 140636 &94- 19329 19329 16077 16077 18241
18241 18241 18241 14063695- f22:p 19 25 26 27 28 29 30 31 20
warning. without bremsstrahlung, flux estimates will be low.96-
fs22 -2 3 -6 797- sd22 1.9329 1.9329 1.6077 1.6077 213600 &98-
1.9329 1.9329 1.6077 1.6077 213600 &99- 1.9329 1.9329 1.6077
1.6077 213600 &
100- 1.9329 1.9329 1.6077 1.6077 213600 &101- 1.9329 1.9329
1.6077 1.6077 213600 &102- 1.9329 1.9329 1.6077 1.6077 213600
&103- 1.9329 1.9329 1.6077 1.6077 213600 &104- 1.9329
1.9329 1.6077 1.6077 213600 &105- 1.9329 1.9329 1.6077 1.6077
213600106- de2 0.01 0.03 0.05 0.07 0.10 0.15 0.20 &107- 0.25
0.30 0.35 0.40 0.45 0.50 0.55 &108- 0.60 0.65 0.70 0.80 1.00
1.40 1.80 &109- 2.20 2.60 2.80 3.25 3.75 4.25 4.75 &110-
5.00 5.25 5.75 6.25 6.75 7.50 9.00 &111- 11.0 13.0 15.0112- df2
LIN 3.96E-06 5.82E-07 2.90E-07 2.58E-07 2.83E-07 3.79E-07 5.01E-07
&113- 6.31e-07 7.59E-07 8.78E-07 9.85E-07 1.08E-06 1.17E-06
1.27E-06 &114- 1.36e-06 1.44E-06 1.52E-06 1.68E-06 1.98E-06
2.51E-06 2.99E-06 &115- 3.42e-06 3.82E-06 4.01E-06 4.41E-06
4.83E-06 5.23E-06 5.60E-06 &116- 5.80e-06 6.01E-06 6.37E-06
6.74E-06 7.11E-06 7.66E-06 8.77E-06 &117- 1.03e-05 1.18E-05
1.33E-05118- de22 0.01 0.03 0.05 0.07 0.10 0.15 0.20 &119- 0.25
0.30 0.35 0.40 0.45 0.50 0.55 &120- 0.60 0.65 0.70 0.80 1.00
1.40 1.80 &121- 2.20 2.60 2.80 3.25 3.75 4.25 4.75 &122-
5.00 5.25 5.75 6.25 6.75 7.50 9.00 &123- 11.0 13.0 15.0124-
df22 LIN 3.96E-06 5.82E-07 2.90E-07 2.58E-07 2.83E-07 3.79E-07
5.01E-07 &125- 6.31e-07 7.59E-07 8.78E-07 9.85E-07 1.08E-06
1.17E-06 1.27E-06 &126- 1.36e-06 1.44E-06 1.52E-06 1.68E-06
1.98E-06 2.51E-06 2.99E-06 &127- 3.42e-06 3.82E-06 4.01E-06
4.41E-06 4.83E-06 5.23E-06 5.60E-06 &128- 5.80e-06 6.01E-06
6.37E-06 6.74E-06 7.11E-06 7.66E-06 8.77E-06 &129- 1.03e-05
1.18E-05 1.33E-05130- nps 50000000
warning. sum of segment sizes differs from total in 18
cases.1cells print table 60
3Calculation No. PM-1083, Rev. 0 Attachment D (Page D21 of
D75)
-
PBrs2.o 5/30/2009
atom gram photoncell mat density density volume mass pieces
importance