ESTS0386 MAGNUM2D.
MAGNUM-2D, Heat Transport and Groundwater Flow in Fractured Porous Media
NAME OR DESIGNATION OF PROGRAM, COMPUTER, DESCRIPTION OF PROGRAM OR FUNCTION, METHOD OF SOLUTION, RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM, TYPICAL RUNNING TIME, UNUSUAL FEATURES OF THE PROGRAM, RELATED AND AUXILIARY PROGRAMS, STATUS, REFERENCES, MACHINE REQUIREMENTS, LANGUAGE, OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED, OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS, NAME AND ESTABLISHMENT OF AUTHORS, MATERIAL, CATEGORIES
[ top ]
[ top ]
To submit a request, click below on the link of the version you wish to order.
Only liaison officers are authorised to submit online requests. Rules for requesters are
available here.
| Program name | Package id | Status | Status date |
|---|---|---|---|
| MAGNUM-2D | ESTS0386/01 | Arrived | 18-APR-2001 |
Machines used:
| Package ID | Orig. computer | Test computer |
|---|---|---|
| ESTS0386/01 | CRAY X-MP |
[ top ]
3. DESCRIPTION OF PROGRAM OR FUNCTION
MAGNUM2D was developed to analyze thermally driven fluid motion in the deep basalts below the Paco Basin at the Westinghouse Hanford Site. Has been used in the Basalt Waste Isolation Project to simulate nonisothermal groundwater flow in a heterogeneous anisotropic medium and heat transport in a water-rock system near a high level nuclear waste repository. Allows three representations of the hydrogeologic system: an equivalent porous continuum, a system of discrete, unfilled, and inter- connecting fractures separated by impervious rock mass, and a low permeability porous continuum with several discrete, unfilled fractures traversing the medium. The calculations assume local thermodynamic equilibrium between the rock and groundwater, non- isothermal Darcian flow in the continuum portions of the rock, and nonisothermal Poiseuille flow in discrete unfilled fractures. In addition, the code accounts for thermal loading within the elements, zero normal gradient and fixed boundary conditions for both temperature and hydraulic head, and simulation of the temperature and flow independently. The Q2DGEOM preprocessor was developed to generate, modify, plot and verify quadratic two dimensional finite element geometries. The BCGEN preprocessor generates the boundary conditions for head and temperature and ICGEN generates the initial conditions. The GRIDDER postprocessor interpolates nonregularly spaced nodal flow and temperature data onto a regular rectangular grid. CONTOUR plots and labels contour lines for a function of two variables and PARAM plots cross sections and time histories for a function of time and one or two spatial variables. NPRINT generates data tables that display the data along horizontal or vertical cross sections. VELPLT differentiates the hydraulic head and buoyancy data and plots the velocity vectors. The PATH postprocessor plots flow paths and computes the corresponding travel times.
MAGNUM2D was developed to analyze thermally driven fluid motion in the deep basalts below the Paco Basin at the Westinghouse Hanford Site. Has been used in the Basalt Waste Isolation Project to simulate nonisothermal groundwater flow in a heterogeneous anisotropic medium and heat transport in a water-rock system near a high level nuclear waste repository. Allows three representations of the hydrogeologic system: an equivalent porous continuum, a system of discrete, unfilled, and inter- connecting fractures separated by impervious rock mass, and a low permeability porous continuum with several discrete, unfilled fractures traversing the medium. The calculations assume local thermodynamic equilibrium between the rock and groundwater, non- isothermal Darcian flow in the continuum portions of the rock, and nonisothermal Poiseuille flow in discrete unfilled fractures. In addition, the code accounts for thermal loading within the elements, zero normal gradient and fixed boundary conditions for both temperature and hydraulic head, and simulation of the temperature and flow independently. The Q2DGEOM preprocessor was developed to generate, modify, plot and verify quadratic two dimensional finite element geometries. The BCGEN preprocessor generates the boundary conditions for head and temperature and ICGEN generates the initial conditions. The GRIDDER postprocessor interpolates nonregularly spaced nodal flow and temperature data onto a regular rectangular grid. CONTOUR plots and labels contour lines for a function of two variables and PARAM plots cross sections and time histories for a function of time and one or two spatial variables. NPRINT generates data tables that display the data along horizontal or vertical cross sections. VELPLT differentiates the hydraulic head and buoyancy data and plots the velocity vectors. The PATH postprocessor plots flow paths and computes the corresponding travel times.
[ top ]
4. METHOD OF SOLUTION
MAGNUM2D is a two-dimensional finite element numerical model for transient or steady state analysis of coupled processes of advective and conductive heat transfers in a porous continuum and buoyancy driven groundwater flow in a saturated, fractured, porous medium. The governing equations consist of a set of coupled, quasilinear partial differential equations that are based on the physical laws of fluid continuity and conservation of momentum and energy. The solutions are based on a Galerkin finite element scheme, using a block diagonal frontal solution technique. The computational network consists of both quadrilateral and triangular elements to model the continuum regions of the spatial domain, and line elements to represent the discrete fractures within the domain. A Newton-Raphson algorithm is embedded in the Galerkin functional to formulate the problems in terms of the incremental changes in the dependent variables. This method allows substantial flexibility in discretization of the spatial domain. The flow equations are solved simultaneously, and the flow and heat transfer equations are solved in alternating sequence. The simulation domain may be modeled in two Cartesian dimensions or in a radial coordinate system that is symmetric around the vertical axis.
MAGNUM2D is a two-dimensional finite element numerical model for transient or steady state analysis of coupled processes of advective and conductive heat transfers in a porous continuum and buoyancy driven groundwater flow in a saturated, fractured, porous medium. The governing equations consist of a set of coupled, quasilinear partial differential equations that are based on the physical laws of fluid continuity and conservation of momentum and energy. The solutions are based on a Galerkin finite element scheme, using a block diagonal frontal solution technique. The computational network consists of both quadrilateral and triangular elements to model the continuum regions of the spatial domain, and line elements to represent the discrete fractures within the domain. A Newton-Raphson algorithm is embedded in the Galerkin functional to formulate the problems in terms of the incremental changes in the dependent variables. This method allows substantial flexibility in discretization of the spatial domain. The flow equations are solved simultaneously, and the flow and heat transfer equations are solved in alternating sequence. The simulation domain may be modeled in two Cartesian dimensions or in a radial coordinate system that is symmetric around the vertical axis.
[ top ]
5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM
Maxima of 3000 nodes, 960 elements, 225 heat source elements, 50 thermal load points, 25 material types, 25 time intervals, 5 binary input files. The number of nodes define the element: 8 nodes specify a quadrilateral, 6 nodes specify a triangular shape and 3 nodes specify a linear element. The element aspect ratio affects the accuracy of the results and should be kept at a value below 100. The reference temperature must be between 0.0 and 300.0 degrees C.
MAGNUM2D requires the proprietary CA-DISSPLA graphics library (see http://www.gaeinc.com) to utilize the graphics post processor; this program is not included.
Maxima of 3000 nodes, 960 elements, 225 heat source elements, 50 thermal load points, 25 material types, 25 time intervals, 5 binary input files. The number of nodes define the element: 8 nodes specify a quadrilateral, 6 nodes specify a triangular shape and 3 nodes specify a linear element. The element aspect ratio affects the accuracy of the results and should be kept at a value below 100. The reference temperature must be between 0.0 and 300.0 degrees C.
MAGNUM2D requires the proprietary CA-DISSPLA graphics library (see http://www.gaeinc.com) to utilize the graphics post processor; this program is not included.
[ top ]
[ top ]
[ top ]
[ top ]
10. REFERENCES
- R.G. Baca, B. Sagar, and D.W. Langford,
MAGNUM2D - A Finite Element Model for Coupled Heat Transport and
Ground Water Flow in Fractured Porous Media: Mathematical Theory,
Numerical Techniques, and Computational Test
WHC-EP-0023 (December 1987).
- L.L. Eyler and M.J. Budden,
Verification and Benchmarking of MAGNUM2D - A Finite Element
Computer Code for Flow and Heat Transfer in Fractured Porous Media
PNL-5237 (March 1985).
- D. Langford,
MAGNUM2D Computer Code: Requirements Document
Westinghouse Hanford Report (August 28, 1987).
- R.G. Baca, B. Sagar, and D.W. Langford,
MAGNUM2D - A Finite Element Model for Coupled Heat Transport and
Ground Water Flow in Fractured Porous Media: Mathematical Theory,
Numerical Techniques, and Computational Test
WHC-EP-0023 (December 1987).
- L.L. Eyler and M.J. Budden,
Verification and Benchmarking of MAGNUM2D - A Finite Element
Computer Code for Flow and Heat Transfer in Fractured Porous Media
PNL-5237 (March 1985).
- D. Langford,
MAGNUM2D Computer Code: Requirements Document
Westinghouse Hanford Report (August 28, 1987).
ESTS0386/01, included references:
- R.L. England, K.J. Ekblad and R.G. Baca:MAGNUM-2D Computer Code - Users Guide
RHO-BW-CR--143P (January 1985)
[ top ]
[ top ]
[ top ]
[ top ]
[ top ]
ESTS0386/01
source program mag tapeCOMDECK MAGNUM Update Routine COMDECKS SRCTPsource program mag tapeMAGNUM.F MAGNUM Source File SRCTP
source program mag tapeCOMDECK FDMIO Update Routine COMDECKS SRCTP
source program mag tapeFDMIO.F FDMIO Source File SRCTP
source program mag tapeCOMMON.F PLTLIB Include File SRCTP
source program mag tapePLTLIB.F PLOT Library Source SRCTP
source program mag tapeCOMMON.F O2DIO Include File SRCTP
source program mag tapeO2DIO.F 2D Quadtratic Elements Routines SRCTP
source program mag tapeCOMMON.F O2DUTIL Include File SRCTP
source program mag tapeO2DUTIL.F 2D Geometry Utilities SRCTP
source program mag tapeCOMMON.F USRLIB Include File SRCTP
source program mag tapeUSRLIB.F General Library Source SRCTP
source program mag tapeCOMMON.F BCGEN.F Include File SRCTP
source program mag tapeBCGEN.F Boundary Value Generator Source SRCTP
source program mag tapeCOMMON.F CONTOURE.F Include File SRCTP
source program mag tapeCONTOURE.F CONTOUR Line Drawing Routine SRCTP
source program mag tapeFEMESH.F MESH Generator Source SRCTP
source program mag tapeCOMMON.F GRIDDER.F Include File SRCTP
source program mag tapeGRIDDER.F Elements to Grid Transl. Source SRCTP
source program mag tapeCOMMON.F ICGEN Include File SRCTP
source program mag tapeICGEN.F Initial Condition Generator SRCTP
source program mag tapeCOMMON.F MESHER Include File SRCTP
source program mag tapeMESHER.F MESH Optimizing Routine SRCTP
source program mag tapeCOMMON.F NPRINT Include File SRCTP
source program mag tapeNPRINT.F CHAINT/MAGNUM Print Routines SRCTP
source program mag tapeCOMMON.F PARAM Include File SRCTP
source program mag tapePARAM.F 2D Data Plot Routine SRCTP
source program mag tapeCOMMON.F VELPLT Include File SRCTP
source program mag tapeVELPLT.F VELOVITY Vector Plotting Routine SRCTP
report RHO-BW-CR--143P (January 1985) REPPT
Keywords: finite element method, flow models, fluid flow, ground water, heat transfer, hydraulics, partial differential equations, porous materials, radioactive waste storage, radionuclide migration, temperature distribution.