skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Massively parallel simulation of flow and transport in variably saturated porous and fractured media

Abstract

This paper describes a massively parallel simulation method and its application for modeling multiphase flow and multicomponent transport in porous and fractured reservoirs. The parallel-computing method has been implemented into the TOUGH2 code and its numerical performance is tested on a Cray T3E-900 and IBM SP. The efficiency and robustness of the parallel-computing algorithm are demonstrated by completing two simulations with more than one million gridblocks, using site-specific data obtained from a site-characterization study. The first application involves the development of a three-dimensional numerical model for flow in the unsaturated zone of Yucca Mountain, Nevada. The second application is the study of tracer/radionuclide transport through fracture-matrix rocks for the same site. The parallel-computing technique enhances modeling capabilities by achieving several-orders-of-magnitude speedup for large-scale and high resolution modeling studies. The resulting modeling results provide many new insights into flow and transport processes that could not be obtained from simulations using the single-CPU simulator.

Authors:
; ;
Publication Date:
Research Org.:
Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
Sponsoring Org.:
USDOE Office of Civilian Radioactive Waste Management. Yucca Mountain Site Characterization Office (US)
OSTI Identifier:
791815
Report Number(s):
LBNL-49407
R&D Project: 80NB71; TRN: US0200693
DOE Contract Number:
AC03-76SF00098
Resource Type:
Conference
Resource Relation:
Conference: XIV International Conference on Computational Methods in Water Resources, Delft (NL), 06/23/2002--06/28/2002; Other Information: PBD: 15 Jan 2002
Country of Publication:
United States
Language:
English
Subject:
12 MANAGEMENT OF RADIOACTIVE WASTES, AND NON-RADIOACTIVE WASTES FROM NUCLEAR FACILITIES; 11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS; 99 GENERAL AND MISCELLANEOUS//MATHEMATICS, COMPUTING, AND INFORMATION SCIENCE; ALGORITHMS; FRACTURED RESERVOIRS; MULTIPHASE FLOW; PERFORMANCE; COMPUTERIZED SIMULATION; SITE CHARACTERIZATION; ENVIRONMENTAL TRANSPORT; WATER RESOURCES; YUCCA MOUNTAIN; RADIOACTIVE WASTE DISPOSAL; T CODES; RADIONUCLIDE MIGRATION; PARALLEL COMPUTING RESERVOIR SIMULATION

Citation Formats

Wu, Yu-Shu, Zhang, Keni, and Pruess, Karsten. Massively parallel simulation of flow and transport in variably saturated porous and fractured media. United States: N. p., 2002. Web.
Wu, Yu-Shu, Zhang, Keni, & Pruess, Karsten. Massively parallel simulation of flow and transport in variably saturated porous and fractured media. United States.
Wu, Yu-Shu, Zhang, Keni, and Pruess, Karsten. 2002. "Massively parallel simulation of flow and transport in variably saturated porous and fractured media". United States. doi:. https://www.osti.gov/servlets/purl/791815.
@article{osti_791815,
title = {Massively parallel simulation of flow and transport in variably saturated porous and fractured media},
author = {Wu, Yu-Shu and Zhang, Keni and Pruess, Karsten},
abstractNote = {This paper describes a massively parallel simulation method and its application for modeling multiphase flow and multicomponent transport in porous and fractured reservoirs. The parallel-computing method has been implemented into the TOUGH2 code and its numerical performance is tested on a Cray T3E-900 and IBM SP. The efficiency and robustness of the parallel-computing algorithm are demonstrated by completing two simulations with more than one million gridblocks, using site-specific data obtained from a site-characterization study. The first application involves the development of a three-dimensional numerical model for flow in the unsaturated zone of Yucca Mountain, Nevada. The second application is the study of tracer/radionuclide transport through fracture-matrix rocks for the same site. The parallel-computing technique enhances modeling capabilities by achieving several-orders-of-magnitude speedup for large-scale and high resolution modeling studies. The resulting modeling results provide many new insights into flow and transport processes that could not be obtained from simulations using the single-CPU simulator.},
doi = {},
journal = {},
number = ,
volume = ,
place = {United States},
year = 2002,
month = 1
}

Conference:
Other availability
Please see Document Availability for additional information on obtaining the full-text document. Library patrons may search WorldCat to identify libraries that hold this conference proceeding.

Save / Share:
  • The transport of chemicals through solid waste disposal piles presents a potential source of surface and groundwater contamination. In an effort to determine the scope of the potential hazard, many laboratory tests have been devised and used on the waste materials. A prerequisite for such a test would be that it yields reproducible results. An added advantage, however, would be that the mechanisms of leachate generation in the test have some resemblance to those in the field.
  • Improved models of contaminant migration in heterogeneous, variably saturated porous media are required to better define the long-term stewardship requirements for U.S. Department of Energy (DOE) lands and to assist in the design of effective vadose-zone barriers to contaminant migrations. The development of these improved models requires field and laboratory results to evaluate their efficacy. However, controlled laboratory experiments simulating vadose conditions can require extensive period of time, and often are conducted at condition near saturation rather than the much drier conditions common in many contaminated arid vadose zone sites. Collaborative research undertaken by the Idaho National Laboratory (INL) andmore » the University of Idaho as part of this Environmental Management Science Program project focused on the development and evaluation of geocentrifuge techniques and equipment that allows vadose zone experiments to be conducted for relevant conditions in time frames not possible in conventional bench top experiments. A key and novel aspect of the research was the use of the 2-meter radius geocentrifuge capabilities at the Idaho National Laboratory to conduct unsaturated transport experiments. Specifically, the following activities were conducted ** Reviewing of the theory of unsaturated flow in the geocentrifuge to establish the range of centrifuge accelerations/experimental conditions and the translation of centrifuge results to 1 gravity applications. ** Designing, constructing, and testing of in-flight experimental apparatus allowing the replication of traditional bench top unsaturated transport experiments on the geocentrifuge. ** Performing unsaturated 1-dimenstional column geocentrifuge experiments using conservative tracers to evaluate the effects of increased centrifugal acceleration on derived transport properties and assessing the scaling relationships for these properties. Because the application of geocentrifuge techniques to vadose transport is in its infancy experimental apparatus such as pumps, flow meters, columns, fraction collectors, etc. that would reliably function under the increased self weight experienced on the centrifuge had to be developed and tested as part of this project. Although, we initially planed to conduct experiments using reactive tracer and 2-dimensional heterogeneities, the cost and time associated with designing, building, and testing of experimental apparatus limited our experimental program to conservative tracer experiments using 1-dimensional columns. The results we obtained in this study indicate that the geocentrifuge technique is a viable experimental method for the study of subsurface processes where gravitational acceleration is important. The geocentrifuge allows experiments to be completed more quickly than tests conducted at 1-g, can be used to experimentally address important scaling issues, and permits experiments under a range of conditions that would be difficult or impossible using conventional approaches. The application of the geocentrifuge approaches and associated models developed in this project allows more meaningful investigation of DOE relevant vadose-zone issues under scalable conditions in time frames previously not obtainable.« less
  • FACT is a finite element based code designed to model subsurface flow and contaminant transport. It was designed to perform transient three-dimensional calculations that simulate isothermal groundwater flow, moisture movement, and solute transport in variably saturated and fully saturated subsurface porous media. The code is designed specifically to handle complex multi-layer and/or heterogenous aquifer systems in an efficient manner and accommodates a wide range of boundary conditions. Additionally 1-D and 2-D (in Cartesian coordinates) problems are handled in FACT by simply limiting the number of elements in a particular direction(s) to one. The governing equations in FACT are formulated onlymore » in Cartesian coordinates. FACT writes out both ascii and graphical binary files that are TECPLOT-ready. Special features are also available within FACT for handling the typical groundwater modeling needs for remediation efforts at the Savannah River Site.« less
  • FACT is a finite element based code designed to model subsurface flow and contaminant transport. It was designed to perform transient three-dimensional calculations that simulate isothermal groundwater flow, moisture movement, and solute transport in variably saturated and fully saturated subsurface porous media. The code is designed specifically to handle complex multi-layer and/or heterogenous aquifer systems in an efficient manner and accommodates a wide range of boundary conditions. Additionally 1-D and 2-D (in Cartesian coordinates) problems are handled in FACT by simply limiting the number of elements in a particular direction(s) to one. The governing equations in FACT are formulated onlymore » in Cartesian coordinates. FACT writes out both ascii and graphical binary files that are TECPLOT-ready. Special features are also available within FACT for handling the typical groundwater modeling needs for remediation efforts at the Savannah River Site.« less
  • FACT is a finite element based code designed to model subsurface flow and contaminant transport. It was designed to perform transient three-dimensional calculations that simulate isothermal groundwater flow, moisture movement, and solute transport in variably saturated and fully saturated subsurface porous media. The code is designed specifically to handle complex multi-layer and/or heterogenous aquifer systems in an efficient manner and accommodates a wide range of boundary conditions. Additionally 1-D and 2-D (in Cartesian coordinates) problems are handled in FACT by simply limiting the number of elements in a particular direction(s) to one. The governing equations in FACT are formulated onlymore » in Cartesian coordinates. FACT writes out both ascii and graphical binary files that are TECPLOT-ready. Special features are also available within FACT for handling the typical groundwater modeling needs for remediation efforts at the Savannah River Site.« less