i I .", J- SANDIA REPORT SAND85-0008 . UC-70 Unlimited Release Printed June 1987 Nevada Nuclear Waste Storage Investigations Project Sandia Implementation of the TRACR3D Flow and Transport Code Nancy K. Prindle Prepared by Sandia National Laboratories Albuquerque, New Mexico 87185 and Uvermore, California 94550 for the United States Department of Energy under Contract DE-AC04-76DP00789 A~~ et Tn~ C: 3 A~~~~~~ r4 o 13 SF290IQ18-81) '. SF2900QI8-81 ) -I
85
Embed
SAND85-0008, 'Sandia Implementation of the TRACR3D Flow ...SANDIA REPORT i I .", J- SAND85-0008 . UC-70 Unlimited Release Printed June 1987 Nevada Nuclear Waste Storage Investigations
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
i I .", J-
SANDIA REPORTSAND85-0008 . UC-70Unlimited ReleasePrinted June 1987
Sandia Implementation of theTRACR3D Flow and Transport Code
Nancy K. Prindle
Prepared bySandia National LaboratoriesAlbuquerque, New Mexico 87185 and Uvermore, California 94550for the United States Department of Energyunder Contract DE-AC04-76DP00789
A~~
et Tn~C:
3
A~~~~~~ r4o
13
SF290IQ18-81) '.
SF2900QI8-81 ) -I
h
'Prepared by Nevada Nuclear Waste Storage Investigations (NNWSI) Pro-ect particp ants as part of the Civilian Radioactive Waste ManagementProgram (FRWM). The NNWSI Project is managed by the Waste Manage-
ment Project Office (WMPO) of the U. S. Department of Energy, NevadaOperations Office (DOE/NV). NNWSI Project work is sponsored by theOffice of Geologic Repositories (OGR) of the DOE Office of Civilian Radioac-tive Waste Management (OCRWM)W'
Issued by Sandia National Laboratories, operated for the United StatesDepartment of Energy by Sandia Corporation.NOTICE: This report was prepared as an account of work sponsored by anagency of the United States Government Neither the United States Govern-ment nor any agency thereof, nor any of their employees, nor any of theircontractors, subcontractors, or their employees, makes any warranty, expressor implied, or assumes any legal liability or responsibility for the accuracy,comrne ness, or usefulness of any information, apparatus, product, or pro-cess disclosed, or represents that its use would not infringe privately ownedright. Reference herein to any specific commercial product, process, orservice by trade name, trademark, manufacturer, or otherwise, does notnecessay consttute or imply it endorsement, recommendation, or favoringby the United States Government, any agency thereof or any of theircontractors or subcontractors. The views and opinions expressed herein donot necessarily state or reflect those of the United States Government, anyagency thereof or any of their contractors or subcontractors.
Printed in the United States of AmericaAvailable fromNational Technical Information ServiceU.S. Department of Commerce5286 Port Royal RoadSpringfield, VA 22161
SANDIA IMPLEMENTATION OF THE TRACR3D FLOW AND TRANSPORT CODE
Nancy K. PrindleNUWSI Repository Performance Assessments
Division 6312
and
Jeffery FosterTechnadyne Engineering Consultants, Inc.
ABSTRACT
TRACR3D is a computer code developed at Los Alamos National Laboratory tomodel three-dimensional fluid flow and mass transport in a porousmedium. The code is being considered for use in performance assessmentmodeling of a potential repository for high-level radioactive wastes atYucca Mountain in lye County, Nevada.* This report describes modifi-cations made during the installation and implementation of the code atSandia National Laboratories and the results of some test problems. Theoperating system commands necessary to direct TRACR3D on the CRAY-lcomputer at Sandia are included, and a flow chart outlining the codestructure is presented.
*The Nevada Nuclear Waste Storage Investigations (NNWSI) Project, managedby the Nevada Operations Office of the U.S. Department of Energy, isexamining the feasibility of siting a repository for high-level nuclearwastes at Yucca Mountain on and adjacent to the Nevada Test Site. Thiswork, intended to extend our understanding of potential radionuclidetransport in the unsaturated geohydrologic units below Yucca Mountain,was funded by the NNWSI project. The ultimate use of this informationwill be to build reasonable assurance in the probabilistic predictionsof postclosure waste isolation.
ii/iv
Contents
Page
A. '1. INTRODUCTION 1
Purpose 1
Background 3
2. INSTALLATION OF TRACR3D AT SUL 5
3. SUMMARY OF TEST PROBLEMS 9
Problem Descriptions 9
Comparison of Results 9
Input Descriptions 21
4. STEPS IN RUNNING TRACR3D AT SUL 23
Basic Steps 23
Other Helpful Procedures 26
5. STRUCTURE OF TRACR3D 29
Subroutine Calls 29
STOP Statements 29
Flowchart 29
6. SUMMARY OF TRACR3D Status AT SUL 36
REFERENCES 38
APPENDIX A--TRACR3D Code Listing as Stored at SUL
APPENDIX B--Input Files for Test Problems
APPENDIX C--Selected Output from Test Problems
APPENDIX D--Parameters from the UNWSI Reference Information Base
and the KNWSI Technical Data Base
V
Illustrations
Figure Page
1 Liquid Permeabilities as Functions of Saturation for 16Test Problem LS6
2 Gas Permeabilities as Functions of Saturation for Test 17Problem LS6
3 Capillary Pressures as Functions of Saturation for Test 18Problem LS6
4-7 LS6 Comparison (Water Pressure vs. Time) 19
8-11 LS6 Comparison (Air Saturation Vs. Time) 20
12 Flowchart of TRACR3D 32-34
Tables
Table Page
1 Characteristics of TRACR3D 2
2 Descriptions of Test Problems 10
3 Material Properties for Tuff in Test Problem LS6 13
4 Material Properties for Sand in Test Problem LS6 14
5 Material Properties for Gravel in Test Problem LS6 15
6 Subroutine Calls in TRACR3D 30
7 STOP Statements in TRACR3D 31
vi
Acknowledgment~s
Bryan Travis and Steve Hodson completed a series of comparative runsat Los Alamos National Laboratory and provided the authors invaluablehelp with debugging efforts. Jack Gauthier and Wynona Sexson alsoprovided considerable help at Sandia National Laboratories in installingand debugging the TRACR3D code.
vii/viii
1. Introduction
Purpose
TRACR3D is a large (5500 FORTRAN lines) finite-difference code
developed at Los Alamos National Laboratory (LANL) by Dr. Bryan J. Travis
(Travis, 1984). The code, originally developed for applications in the
oil-shale industry, models isothermal, transient fluid and gas flow and
chemical transport in three dimensions in a deformable, heterogeneous,
reactive porous medium. Since its original development, special-purpose
subroutines have been added to the code to model fluid flow and transport
of radioactive chemical species in porous materials containing discrete
fractures. The general characteristics of the TRACR3D code are listed in
Table 1.
The features that have been added to the original code render the
TRACR3D code useful for modeling simultaneous air and water flow and
accompanying radionuclide transport in three-dimensional fractured,
porous media. In particular, TRACR3D has been used at LANL for analyses
of isothermal radionuclide transport in both aqueous and vapor phases
through variably saturated, fractured tuffs such as those that make up
the prospective site for the disposal of high-level radioactive waste at
Yucca Mountain in Nye County, Nevada. The purpose of these analyses was
to investigate the phenomena that affect isothermal radionuclide
transport on relatively small scales (on the order of 1 to 10 meters) by
doing sensitivity studies,. experimental analysis, data interpretation,
and development of "macroscopic" models of bulk behavior. The capability
to do similar analyses in three dimensions can be useful to the
performance assessment group at Sandia National Laboratories (SKL) for
their studies for the Nevada Nuclear Waste Storage Investigations (NWSI)
Project. Thus, a version of TRACR3D was provided to SNL for installation
on the SNL CRAY-1 computer. This version, which corresponds to the
version described in the TRACR3D user's manual, has been stored as T3D784
at SUL.
-1-
Table 1
Characteristics of TRACR3D
DIMENSION
1-, 2-, or 3-D
NUMERICAL TECHNIQUE
Finite difference combined with method of characteristics,orthogonal elements
GOVERNING EQUATIONS
Hydrologic: Air- and water-mass conservation, Forcheimer's equationfor momenta conservation of two fluids, reducing to Richards'equation at low Reynolds numbers
Implicit, iterative, successive overrelaxation techniques forhydrologic solution, explicit or Runge-Kutta for soluteconcentration solutions
BOUNDARY CONDITIONS
Specified pressure, fluid saturation or fluid flux, free flow,ponding, specified concentration or solute flux, leaching, bandrelease; can include radioactive decay
The new parameter NXX was then inserted into the following dimension
statement in the subroutine BONDRY (modifications to the published
version are underlined):
DIMENSION AVXYZ(NXX,3),IJKS(3),BB(10),ITDBDY(l)
The constraints on the relative permeabilities and the modifications
to the PARAMETER list and DIMENSION statements were installed in the
TRACR3D code at LANL, and a new version was supplied to SNL. The 3STE
problem was then run to a normal completion at SNL, and the output
compared exactly with that produced at LANL. The five other test
problems were then run and the results compared. These comparisons are
discussed in the next section and provide further installation testing of
the code. The test problems GENl, GEN2, and PIG3 provide checks on the
numerical accuracy of the code because analytic solutions are available.
The test problems ESTE, VARP, and LS6 provide sensitivity comparisons for
the code by comparing the solutions produced on two different computer
systems.
-8-
3. Summary of Test Problems
Problem Descriptions
The six problems that were run to test the installation of TRACR3D at
SNL and to do initial benchmarking for numerical accuracy are described
in the user's manual (Travis, 1984) and its supporting references (see
REFERENCES). The first problem attempted at SUL was ESTE, which is the
sample problem presented in Appendix C of the published TRACR3D user's
manual. The five other problems described in the user's manual are
analyses done by Travis for verification ind validatLon of TRACR3D. The
input and output to these problems were generously provided to the
authors by Travis to test further the SUL implementation of TRACR3D and
to do some benchmarking compari4ons Table , gives a brief description
of each of these problems. Collectively, :these test problems call on the
subroutines that might be used for-caJiculating isothermal water flow and,- , Jrw A.* ,
radionuclide transport in unsatrpated or saturated, fractured, porous
media.
Comparison of Results
The results generated with TBACR3D at SUL and LAUL have been compared
in several ways. The code calculates not only saturation, pressure, and
concentration fields, but also controls features such as the time steps,
the number of iterations per time step, and residual errors. Variable
property values are calculated. These calculated variables should
compare exactly within the same time steps and cells if the code is
running exactly the same on different systems. As we have already noted,
however, differences in the compilers can cause minor differences in the
results. These differences can be expected to be insignificant, in most
cases. The experiences with running the sample problem ESTE, however,
suggested that it would be prudent to investigate whether there were any
further significant differences in the way the TRACR3D code was running
on the SNL and LARL computer systems. Thus, the number of time steps,
-9-
Table 2
Descriptions of Test Problems
1. ESTEThree-dimensional, transient water flow and transport of anonreactive solute in a low-permeability, partially saturated,porous medium with a single vertical, planar fracture. Thissample problem is fully described in Appendix C of the TRACR3Duser's manual (Travis, 1984).
2. GENiOne-dimensional transport of nonreactive solute by steady-statewater flow and molecular diffusion in a saturated, homogeneous,porous medium (Van Genuchten, 1976).
3. GEN2One-dimensional transport of a radioactive solute by steady-state water flow and molecular diffusion in a saturated,homogeneous, porous medium; equilibrium sorption assumed. Thisproblem was analytically solved by Genuchten (Van Genuchten,1981).
4. VARPOne-dimensional, transient water flow in a homogeneous,saturated porous medium using pressure-dependent permeability(Horning, 1977).
5. PIG3One-dimensional transport of a three-member chain ofradioactive solutes by steady-state water flow in a saturated,homogeneous, porous medium; equilibrium sorption assumed(Pigford, 1980).
6. LS6Three-dimensional, transient water drainage through layers ofcrushed tuff, sand, and gravel using tables for moisturecontents as functions of pressure and the Brooks-Corey modelfor pressure-dependent permeabilities. This problem is ananalysis of a field experiment at LANL (Perkins and Travis,1985).
-10-
time-step sizes, and number of iterations per time step were compared for
each run.' If these compared favorably, the results for the saturation,
pressure, and concentration solutions were compared, depending on the
problem. Finally, mass balances were compared. Selected output from
these problems is given in the tables in Appendix C; a general discussion
of the results and'conclusions based on these results follows in this
section.
For the first four test problems, the variables that control
time-step sizes and numbers of iterations were identical out to four
decimal places, as were calculated water-saturation values. Tabulated
values in both space and time (where appropriate) for pressures and mass
concentrations were the same out to three significant digits. In
particular, water saturations and tracer concentrations as functions of
time and three-dimensional space were the same to three significant
digits for the ESTE problem.' (Complete tables of these results are
published in the TRACR3D user's manual.) For both the GENI and GEN2
problems, the one-dimensional, steady-state results for the tracer
concentrations as'functions of space were also the same to three
significant digits. The transient, three-dimensional solution for
pressures compared similarly well in the VARP problem.
In the fifth problem, PIG3, all but the initial value of the third
member of the decay chain compared exactly to three significant digits.
The initial value for the concentration of the third member of the decay
chain was computed to be -8.54 x 10 at LANL and 2.2 x 10 at
SNL. This discrepancy in the initial value was traced to the subroutine
BOUNDRY, where the Bateman equations for radioactive decay reduce to
several terms that should theoretically cancel to zero. However, because
of the numerical solution scheme, differences were occurring in the 14th
significant digit and beyond, resulting in the nonzero results. If
necessary, a simple correction for this error could be implemented in
TRACR3D by putting a practical lower limit on the variable SUN, below
which SUM is set to zero. However, because of the otherwise excellent
comparison of the PIG3 results, we did not think it necessary to rerun
the calculation at SNL using such a correction.
-11-
The comparisons for the sixth test problem, LS6, were more difficult
to interpret than the others . Results show an excellent match except
for values at or near material interfaces. In LS6, a field experiment ip
modeled where water flows through a vertical caisson 600 cm deep. The
caisson is filled with crushed tuff except at the bottom of the caisson,
where a 9-in. layer of sand overlies a 9-in. layer of gravel. Between
the three distinct materials--tuff, sand, and gravel--there are two
interfaces where material properties change significantly. These
properties, water and gas permeabilities and capillary pressure, are
functions of saturation, as shown in Tables 3 through 5 and in Figures 1
through 3.
In the LS6 problem, a pulse of water was introduced at the top of the
caisson and monitored as it infiltrated downward. In modeling the
movement and distribution of the water pulse in the caisson, the LANL
CRAY and the SUL CRAY used different numbers of iterations within time
steps to converge to a pressure solution. Since TRACR3D automatically
cuts down the time step after a certain number of iterations have been
performed, different time-step patterns were established, making
comparisons between tabulated outputs difficult. The results for water
pressures and air saturations were therefore plotted as functions of time
at four locations--two points in the tuff, the material interface between
the sand and tuff, and the gravel at the boundary node (which is
typically one of the most sensitive nodes in a calculation). These are
illustrated in Figures 4 through 11.
Values obtained for water pressure and air saturation for two
locations in the tuff were almost identical at times near each other
(Figures 4, 5, 8, 9). For the location in the 9-in. layer of sand (at
the interface with tuff), the results also matched well (Figures 6, 10).
However, at first glance, there appeared to be considerable discrepancy
in the results in the gravel at the boundary node (Figures 7, 11). These
discrepancies have been ascribed to sensitivity in the calculation of
permeabilities at this interface. The permeabilities in this region are
extremely nonlinear, changing by more than 14 orders of magnitude over
the water saturation range from 0.0 to 0.1 (see Figure 1). Considering
The TRACR3D code is stored on the CRAY-1 at SNL under the filename
T3D784. This file has been used to run the six problems described in
this report. The successful comparisons between the runs at SNL and LAUL
extend the verification work done at LAmL to the SNL implementation.
File T3D784 can be used to solve the conservation and momentum equations
for one- two- or three-dimensional, transient or steady-state
distributions of water, gas, and radioactive solutes in a reacting,
porous medium. The reactions of the solute with the medium may be either
equilibrium or nonequilibrium.
In the process of installing and implementing the code at SUL, the
authors examined the structure of TRACR3D and performed sensitivity
studies of some of the control parameters. The currently stored version
seems to be free of bugs and to run effectively on the CRAY-1 at SNL.
The code will not be maintained by the authors, however. It is
conceivable (if not probable) that in future applications of the code,
heretofore undiscovered problems may be found. Therefore, some parts of
the code structure (i.e., subroutine calls, program logic flow, and
program STOP statements) have been reproduced in this report to
-36-
supplement the information in the user's manual. In addition, Appendix B
lists the input, UPDATE, and COS control files as they were used to run
the six-test problems; Appendix C lists selected output from the test.
problems. This material can be used to provide supplemental example
problems for new users of TRACR3D, to reproduce the results described in
this report at other installations, and to verify any new versions of the
code that may be developed at SUL.
-37-
REFERENCES
"Computing Division Programmer's Information Manual -- PIN-7, CTSSOVERVIEW," LA-5525-M, Vol. 7, Los Alamos National Laboratories, LosAlamos, NH, September, 1982.
"CRAY-1 OS Version 1 Reference Manual," SR-0011, Cray Research Inc.,Mendota Heights, MN, 1982.
"CRAY-1 Computer Systems Library Reference Manual," SR-0014, CrayResearch Inc., Mendota Heights, MU, June, 1982.
"CRAY-1 Computer Systems FORTRAN (CFT) Reference Manual," SR-0009, CrayResearch Inc., Mendota Heights, MN, November, 1982.
Hayden, N. K., "Benchmarking NNWSI Flow and Transport Codes: COVE 1Results," SAND84-0996, Sandia National Laboratories, Albuquerque, NH,June, 1985.
Horning, U., " A Numerical Solution for the Simulation of UnsteadyGroundwater Flow in Both Saturated and Unsaturated Soils," Soil Sci 124,pp. 140-144 (1977).
Perkins, B. A., and B. J. Travis, "Validation of the TRACR3D Code forSoil Water Flow Under Saturated/Unsaturated Conditions in ThreeExperiments," LA-10263-MS, Los Alamos National Laboratory, Los Alamos,NM, January 1985.
Pigford, T. H., et al., "Migration of Radionuclides Through Sorbing MediaAnalytical Solutions -- II," LBL-11616, Lawrence Berkeley Laboratories,Berkeley, CA, 1980.
Samuelson, N. H., and C.D. Brown, "VAX/VMS Remote Job Entry (RJE) User'sGuide," SAND81-2476, Sandia National Laboratories, Albuquerque, NM,October, 1982.
"SANDIA CRAY-1 Supplements," Computer Consulting and Training Division2614, Sandia National Laboratories, Albuquerque, NM, April, 1984.
Silling, S. A., "Final Technical Position on Documentation of ComputerCodes for High-Level Waste Management," NUREG-0856, U.S. NuclearRegulatory Commission, Washington, DC, 1983.
Travis, B. J., "TRACR3D: A Model of Flow and Transport in Porous/Fractured Media," LA-9667-MS, Los Alamos National Laboratories, LosAlamos, NM, May, 1984.
"UPDATE Reference Manual" Report 60342500, Control Data Corporation,Sunnyvale, CA, 1975.
U.S. DOE NVO (United States Department of Energy Nevada OperationsOffice), "Nevada Nuclear Waste Storage Investigations Quality AssurancePlan Revision 4," NVO-196-17, Las Vegas, NV, January, 1986.
-38-
Van Genuchten, M. Th., "Analytical Solutions for Chemical Transport withSimultaneous Adsorption, Zero-Order Production and First Order Decay," J.Hydrol. 49, 213-233 (1981).
Van Genuchten, M. Th., "On the Accuracy and Efficiency of SeveralNumerical Schemes for Solving the Convective-Dispersive Equation," paperpresented at the International Conference on Finite Elements in WaterResources, Princeton University, Princeton, NJ (1976).
-39-/-40-
Appendix A
TRACR3D Code Listing As Stored at SNL
The TRACR3D code has been stored as permanent file T3D784on the CRAY-1 under COS. This file is in UPDATE format foruse as described in this report. The same file is stored onthe IFS in VAX native text for backup purposes. These filescorrespond to the listing in this appendix. The second fileis accessible by the MASS utility for use on CTSS where it canbe run by the CTSS UPDATE facility or by the utility programHISTORN, as at LANL. In addition, control procedures can bedeveloped to run TRACR3D interactively from CTSS or to submitTRACR3D to CTSS from the user's VAX system.
Since TRACR3D is itself not interactive and program runsare most often lengthy, interactive usage is not a significantadvantage and may require a user-modifiable FORTRAN (non-UPDATE) file. Job submission via UPDATE allows necessary usermodifications while maintaining stricter control over the codefor quality assurance purposes, including audit trails of codemodifications.
The complete TRACR3D code listing for this appendix is 90pages long and has been reproduced on microfiche included atthe end of this report.
A-01/-A-02
Appendix B
Input Files for Test Problems
The tables in this appendix contain listings of the threefiles that were used to run each test problem, the COS controlfile, the UPDATE directives file, and the input data file.
1. ESTECOS Command FileUPDATE FileInput Data File
2. GEN1COS Command FileUPDATE FileInput Data File
3. GEN2COS Command FileUPDATE FileInput Data File
4. VARPCOS Command FileUPDATE FileInput Data File
5. PIG3COS Command FileUPDATE FileInput Data File
6. LS6COS Command FileUPDATE FileInput Data File
B-O1
Table B.1
Sample Problem ESTE
COS Control File
ESTEST2,T599,STSCZ.USER, your username,{Enter NOS Password).SEND, your name, your box number.CHARGE(your charge number)ACCESS, DN=T3D`784,NA.UPDATE, P=T3Y784, IN, F.CFT,I-$CPL,L=0.LUR.DISPOSE,DN=$OUT,SDN=TEMP,DC=ST,DEFER,TEXT=
GEN11,T599,STSCZ.USER, your username,{Enter NOS Password).SEND, your name, your box number.CHARGE(your charge number)ACCESS, DN-T3D784,NA.UPDATE, P-T3D784,IN,F.CFT,I=$CPL,L=O.LDR.DISPOSE,DN=$OUT,SDN=TEMP,DC-STDEFER,TEXT-
GEN22,T599,STSCZ.USER, your username,{Enter NOS Password}.SEND, your name, your box number.CHARGE(your charge number)ACCESS, DN=T3D784,NA.UPDATE, P=T3D784,IN,F.CFTI-$CPLL=O.LDR.DISPOSE,DN=$OUT,SDN-TEMPDC-ST,DEFER,TEXT-
VARPER,T599,STSCZ.USER, your username,{Enter NOS Password}.SEND, your name, your box number.CHARGE(your charge number)ACCESS, DN-T3D784,NA.UPDATE, P=T3D784,IN,F.CFT,I=$CPL,L=O.LDR.DISPOSE, N=$OUT,SDN-TEMP,DC-ST,DEFER,TEXT-
PIG3TRC,T1200,STSCZ.USER, your username,{Enter NOS Password).SEND, your name, your box number.CHARGE(your charge number)ACCESS, DN-T3D784,NA.UPDATE, P-T3D784,IN,F.CFT,I-$CPL,L=O.LDR.DISPOSE,DN-$OUT,SDN-TEMPDC-ST,DEFER,TEXT-
LS6CAIS,T2400,STSCZ.USER, your username,{Enter NOS Password).SEND, your name, your box number.CHARGE(your charge number)ACCESS, DN-T3D784,NA.UPDATE, P-T3D784, IN, F.CFT,I=$CPL,L=O.LDR.DISPOSE,DN=$OUT,SDN-TEMP,DC=ST,DEFER,TEXT-
The amount of output from any run of TRACR3Dis controlled by the user. The output returnedfrom runs of the six test problems presented inthis report amounted to 1400 pages, with a likeamount produced at LANL for comparison. A smallportion of this output is included here for eachproblem as examples for the TRACR3D user. Completelistings of all of this output have been stored inthe NNWSI Records Center.
1. Sample Time-Step Log from ESTE2. ESTE Results3. GEN1 Results4. GEN2 Results5. VARP Results6. PIG3 Results7. LS6 Results
Note: The log file as output by TRACR3D includes onlythe 20 lines (often much more) for the time stepscalculated during that program run. The headingswere added for purposes of this report.
C-02
Table C.2
ESTE Results
Output from problem ESTE is published in itsentirety in Travis (1984) (IA-9667-MS), Appendix Dand is not repeated here. Appendix C of thatreference includes a description of the TRACR3Doutput tables produced by subroutine APRINT ascontrolled by input variables TABTP, IPRNT, JPPNT,and KPRNT. Tables C.3c, C.4b, C.5b, C.6b, and C.7bin this appendix are examples of this type of output.
NOTE: The originals of these tables were of poor quality.
CNANGKS TO INITIAU ADbUOKBE CONCLNRATIDN tG-MI/GN)
CRANRCS TO INITIAU GAS SATURATION0. 2. 0, 20 0. 51. O.OOOOOE#00 0.
CHANGESTINITIAh VRESSURE C_8ATK) 0STRIBUTION0. 2. O. 2. O. 51. 0.1OooE*0o 0.
TIME DEPENDENT BOUNDARY NUMBER AND LOCATION. 10 -1 3BO
. * - L~7ESb ART SET UP FOLLOWS
* I
I
. I- _sI
1
* I
* I-I
i. -I
Table C.3b
First Page of Output from TRACR3D
This page recapitulates the input data for theproblem. (The two warnings from CFT are theresult of two extraneous characters at the endof the FORTRAN source code and have since beendeleted.)
-l* a
4
GENUCHTEN TESTp 1-0&I^qC 514'UU R N CIEE -- JLo SIE 3.
Z LEVEL 16 Z a 560699-
_ _ _ _ ,
MASS CGM)-GASw LRB4,F~,BK 0.00+OOO 0.00001+UO 0.OOO00E+O 0OOOOEOOO O. 0 E0000+00 0,Ou00i0oMASS t1aiM-uIUMU- IRBIYbK 0.0000C+00 0.OOUO-fOO 9.4656L+02 094656ii02 0,O0~i0E+O 090000u.0+0
14i4J1> --. 1450 C5P EXS,1434,2 5,3487 CSP END or JOB _jo
14i43153 5---sJ40@7 CSP14:43.62 5.5481 CSP
-14 43152 5.34Y9 USER JOB NAME - GEN1K6F14143t52 5.3489 USER USER NUMBER - _ _ b6JFOS
-rim 3157 rj5 i UsER TIgE TrCcUYxZ-§N cPU - O0ObTIAoOTU8F.1 514:43852 563489 UbER TIME WAITING TO EXECUTE 0000:00126.23351414152 5.3459 USER TI WAITING FO 1/ --- 000i:o600658v B a14343S52 5.3490 USER TIME WAITING IN INPUT OUEUE - 00012b9:23.212214143852 503490 USER MEMORY * CPU TIME CMWDS*SEC) - 0.4516114i43tS2 5*3490 USER MEMORY * I/O WAIT TIME (MWDS*SEC) ' 0.54705143152 5.3490 USER MINIMUM MEMORY WORDS USED & 5529614:43152 5.3491 USER MAXIMUM MEMORY WORDS USED - 1208321414385J v.3491 USLR DftiK-3MRfSOR7 1VE -14143:52 563491 USER USER I/O REQUESTS - 209 Table C.3d1414355w 5.3491 USER OPEN CALLS - 3214t43t52 5.3491 USER CLOSE CALLS - 18
-1T" 3I 5.5491 w U bLR W dR WREySIE WAES0--- Dayfile from COS for Problem GEN1 -
14143:52 5*3491 USER TEMPORARY DATASET SECTORS USEU - S14143152 33491-T USEp P MANEN bAET SiCT0RS AC1SSL -6
14143:52 8.3491 USER PERMANENT DATASET SECTORS SAVED - O _
14:41852 5.3491 USER SECtURS RECSIVED FRUM FHRCi--ND- --- 014:43:52 5*3492 USER SECTORS QUEUED ro roONT END - 0
7-rAlf7 N. lia isr5 dncrT h-icr- -n---i a l--
Vai.I. a Uww.&"VwtKV - 'rTI c F~rcol71voo &T__ Ads
PANAME $AZIM -WARNING - STATEM&NT NUMBER ON BLANK CARD IGNORED *s*~ww*~w***Us**s6*;p#***$*P0343w"Ip I* a
° LE 382 11MI 0.406211+0* VTx 0*6b1011403 P-l1ER- 6 P-tGL= 0.00 PCL-TOL' 0.76 PASS BAI.-al*297600 Ct- 0.151CICLE 353 Ii?! u.40703E40b COx 0.81721SE03 PFIlEIR3 b P-1GL- 00 PCLuI0I 0.94 PAbb UAL981.e29414 CIM 0,151CYCLE 384 liPF *..40w0*06 Ct- 0.q9gf 60 6 3 P-1!I43 a P-CL' 0.00 PCL-!ItL 0.94 PASS BAL..19295001 CIO 0.190CYCLE 385 TIPF ! ..40l99E*0Q 6T O l E6 I04 P-IESR 8 P-TOLU 0.00 PCL-TLUR 094 PASS B"L.a1.293305 CIO 0018VCYCLE 38h I1PE 'J.41060940b Ct- tt.14121?#04 F-Jllhw 14 P-1CL- 0.00 PCL-ICLI' 0,91 PASS, MALt*291407 CIN 0.307CYCLE 387 lIP! 0.41219E*Oe LT= 0.16946L4U4 P-12E14 29 1-IOLN 000 PCLUXULN 05 PAS 9=,9294755 On 0,591CICLE 38b lIPE L.41402E+06 C7- 0.172%01404 F-JIltl 21 P-tOL 0.00 PCL-TOLa 1.00 PASS BAL.s1.292252 C7' 09436CICL0 389 lipI 1*41l610E+Ob Cia 0,207480L04 F-I1FI. 24 P-lOLN 0.00 PCL-IOL' 0.96 PASS HAL.s1.289267 C!. 0.496CYCLE. 39u lI!PF t.41843E *06 LI. 0.23367L404 F-I1l1CICLF 391 11PE u.42077FO66 CT- 0.233o6E#04 F-IIEICYCLFS 392ll .42311!.Ob Ct- 0.23367t#04 F-IIEICYCLE 391 TIPE o 42b440+06 t-. 0 233671,04 P-niF Table C.7oCYCLE 394 hIA! 0.47b67F!06 Ut- 0.430041403 F-IllMCICLE 395 JIME u.42639E40C CT- 0.5160R43 F-IIE I Diagnostic Output from TRACR3D for Problem LS6CICLIF 390 11FF bej427031406 CI. 0.61925L.03 F-JI1 7tJ.&CYCLE 3497ITFE U.4272VE4to Ct- 0.2w519E.03 P- IttCYCLE 390 IlrF U.42759F*06 CLt 0.31823U*03 F-ITUl Problem LS6 went through 2804 cycles, producing 2804CICLE jq3 TIPPF ',.4279E406 CTs .7 I 0F3 F-Iif lines of output like these as well as some like theCYCLE 400 IIPF u.42843140Of LI. 0.4582SC403 F-11ilCYCLE 40) I1AE U.42b95bE46 LIS 0.549C01403 F-11i interruption shown here from statement number 162.C _CLE 402 TIPF C.42964E+06 CT- 0.6598Fo0 F-1ll (see Table C.5c).CYCLE 403 hIRE C#430u44 ae Ct. 0.79186t#03 F-11EiCYCLE 404 TIPE 0,43139h#V6 Lt- t.95023,#03 F-JilCYCLE 465 21 AF 0.43253E+06 CID 0Ti1 6'Jt.U4 a-u. 11 P-ICL. O.00 PCL-!0I's 0.68 PASS BAL,.1,506071 C!. 0s247
Appendix D
Parameters from the NNWSI Reference Information Baseand the NNWSI Technical Data Base
The six test problems described in this reportwere used to implement the TRACR3D code at SNL andwere not intended to produce any results based onthe Yucca Mountain material properties or directlyapplicable to the planned repository. Therefore,no parameters from the NNWSI Reference InformationBase or NNWSI Technical Data Base were used.
Data used in this report was identical to thatused for the same test problems as run at LANL toprovide comparable output.
4
D-O1/D-02
DISTRIBUTION LIST
a
B. C. Rusche (RW-1)
DirectorOffice of Civilian Radioactive
Waste ManagementU.S. Department of Energy
Forrestal BuildingWashington, DC 20585
Ralph Stein (RW-23)Office of Civilian Radioactive
Waste ManagementU.S. Department of Energy
Forrestal BuildingWashington, DC 20585
T. P. Longo (RW-25)Program Management Division
Office of Geologic Repositories
U.S. Department of Energy
Forrestal BuildingWashington, DC 20585
Cy Klingsberg (RW-24)
Geosciences and Technology Division
Office of Geologic Repositories
U. S. Department of Energy
Forrestal BuildingWashington, DC 20585
J. J. Fiore, (RW-221)office of Civilian Radioactive
Waste ManagementU.S. Department of Energy
Forrestal BuildingWashington, DC 20585
B. G. Gale (RW-223)Office of Civilian Radioactive
Waste ManagementU.S. Department of Energy
Forrestal Building
Washington, DC 20585
M. W. Frei (RW-231)Office of Civilian Radioactive
Waste ManagementU.S. Department of Energy
Forrestal BuildingWashington, DC 20585
R. J. Blaney (RW-22)Program Management Division
Office of Geologic Repositories
U.S. Department of Energy
Forrestal BuildingI Washington, DC 20585
E. S. Burton (RW-25)Siting DivisionOffice of Geologic Repositories
U.S. Department of Energy
Forrestal BuildingWashington, D.C. 20585
R. W. Gale (RW-40)Office of Civilian Radioactive
Waste ManagementU.S. Department of Energy
Forrestal BuildingWashington, D.C. 20585
C. R. Cooley (RW-24)Geosciences & Technology Division
Office of Geologic Repositories
U.S. Department of Energy
Forrestal BuildingWashington, DC 20585
V. J. Cassella@(RW-222)Office of Civilian Radioactive
Waste ManagementU.S. Department of Energy
Forrestal BuildingWashington, DC 20585
J. E. Shaheen (RW-44)Outreach ProgramsOffice of Policy, Integration and
OutreachU.S. Department of Energy
Forrestal BuildingWashington, DC 20585
J. 0. Neff, ManagerSalt Repository Project Office
U.S. Department of Energy
505 King AvenueColumbus, OH 43201
J
E-1
D. C. Newton (RW-23)Engineering & Licensing DivisionOffice of Geologic RepositoriesU.S. Department of EnergyForrestal BuildingWashington, DC 20585
0. L. Olson, ManagerBasalt Waste Isolation Project OfficeRichland Operations OfficeU.S. Department of EnergyPost Office Box 550Richland, WA 99352
D. L. Vieth, Director (4)Waste Management Project OfficeU.S. Department of EnergyPost Office Box 14100Las Vegas, NV 89114
S. A. Mann, ManagerCrystalline Rock Project OfficeU.S. Department of Energy9800 South Cass AvenueArgonne, IL 60439
K. Street, Jr.Lawrence Livermore NationalLaboratory
Post Office Box 808Hail Stop L-209Livermore, CA 94550
L. D. Ramspott (3)Technical Project Officer for NNWSILawrence Livermore NationalLaboratory
P.O. Box 808Mail Stop L-204Livermore, CA 94550
D. F. Miller, DirectorOffice of Public AffairsU.S. Department of EnergyPost Office Box 14100U.S. Department of EnergyLas Vegas, NV 89114
W. J. Purcell (RW-20)Associate DirectorOffice of Civilian RadioactiveWaste Hanagement
U.S. Department of EnergyForrestal BuildingWashington, DC 20585
P. H. Bodin (12)Office of Public AffairsU.S. Department of EnergyPost Office Box 14100Las Vegas, NV 89114
B. W. Church, DirectorHealth Physics DivisionU.S. Department of EnergyPost Office Box 14100Las Vegas, NV 89114
V. H. GlanzmanU.S. Geological SurveyPost Office Box 25046913 Federal CenterDenver, CO 80225
E-2
P. T. PrestholtNRC Site Representative1050 East Flamingo RoadSuite 319Las Vegas, NV 89109
H. E. SpaethTechnical Project Officer for NNWSIScience Applications
International CorporationSuite 407101 Convention Center DriveLas Vegas, NV 89109
SAIC-T&HSS Library (2)Science Applications
International CorporationSuite 407101 Convention Center DriveLas Vegas, NV 89109
W. S. Twenhofel, ConsultantScience Applications
International Corp.820 Estes StreetLakewood, CO 89215
A. E. GurrolaGeneral ManagerEnergy Support DivisionHolmes & Narver, Inc.Hail Stop 580Post Office Box 14340Las Vegas, NV 89114
J. A. Cross, HanagerLas Vegas BranchFenix & Scisson, Inc.Hail Stop 514Post Office Box 14308Las Vegas, NV 89114
Neal Duncan (RW-44)Office of Policy, Integration, and
OutreachU.S. Department of EnergyForrestal BuildingWashington, DC 20585
J. S. WrightTechnical Project Officer for NNWSIWestinghouse Electric CorporationWaste Technology Services DivisionNevada operationsPost Office Box 708Hail Stop 703Mercury, NV 89023
ONWI LibraryBattelle Columbus LaboratoryOffice of Nuclear Waste Isolation505 King AvenueColumbus, OH 43201
W. H. Hewitt, Program ManagerRoy F. Weston, Inc.955 L'Enfant Plaza, Southwest, Suite 800Washington, DC 20024
H. D. CunninghamGeneral ManagerReynolds Electrical &Engineering Co., Inc.
Post Office Box 14400Hail Stop 555Las Vegas, NV 89114
T. Hay, Executive AssistantOffice of the GovernorState of NevadaCapitol ComplexCarson City, NV 89710
R. R. Loux, Jr., Executive Director (3)Nuclear Waste Project OfficeState of NevadaEvergreen Center, Suite 2521802 North Carson StreetCarson City, NV 89701
C. H. Johnson, TechnicalProgram ManagerNuclear Waste Project OfficeState of NevadaEvergreen Center, Suite 2521802 North Carson StreetCarson City, NV 89701
E 3
John FordhamDesert Research InstituteWater Resources CenterPost Office Box 60220Reno, NV 89506
Department of ComprehensivePlanning
Clark County225 Bridger Avenue, 7th FloorLas Vegas, NEV 89155
Lincoln County CommissionLincoln CountyPost Office Box 90Pioche, NV 89043
Community Planning andDevelopment
City of North Las VegasPost Office Box 4086North Las Vegas, NV 89030
City ManagerCity of HendersonHenderson, NV 89015
N. A. NormanProject ManagerBechtel National Inc.P. 0. Box 3965San Francisco, CA 94119
D. K. GartlingR. R. EatonR. W. LynchT. 0. HunterRMWICF22/12144/04-1872/REP-II/QIII72-12143-0/QIIIA. L. StevensC. MbraF. W. BinghamG. E. BarrJ. Foster (6)J. GauthierM. L. GreenB. S. IangkopfW. B. MillerR. R. PetersA. C. PetersonR. W. PrindleN. K. PrindleM. S. TierneyM. WilsonJ. G. YeagerT. E. BlejwasE. A. KlavetterJ. R. Tillerson
6315 S. Sinnock6315 P. C. Kaplan6316 R. B. Pcpe6330 W. D. Weart6331 A. R. Lappin6332 L. D. Tyler6332 WE Library (20)6334 D. R. Anderson6410 N. R. Ortiz3141 S. A. Landenberger3151 W. L. Garner (3)8024 P. W. Dean3154-3 C. H. Dalin (28)
for DOSE/oTI
(5)
i_
E-6
Item Sandia Implementation of the TRACR3D Flow and TransportDescription: Code
Availability: X Publicly Availableo Non-Publicly Available
Sensitivity: X Non-Sensitive o Non-Sensitive-Copyrighto Sensitive o Sensitive-Copyright
Electronic MicroficheMedia Type:(If applicable)
Contact: US Nuclear Regulatory CommissionOffice of Nuclear Materials Safety and SafeguardsYucca Mountain Project Manager
Storage/File US Nuclear Regulatory CommissionLocation: Office of Nuclear Materials Safety and Safeguards
Two White Flint North Room T7- E3411545 Rockville PikeRockville, Maryland 20852-2738