HST3D; a computer code for simulation of heat and solute transport in three-dimensional ground-water flow systems

The Heatand Solute-Transport Program (HST3D) simulates ground-water flow and associated heat and solute transport in three dimensions. The HST3D program may be used for analysis of problems such as those related to subsurface-waste injection, landfill leaching, saltwater intrusion, freshwater recharge and recovery, radioactive-waste disposal, hot-water geothermal systems, and subsurface-energy storage. The three governing equations are coupled through the interstitial pore velocity, the dependence of the fluid density on pressure, temperature, and solute-mass fraction, and the dependence of the fluid viscosity on temperature and solute-mass fraction. The solutetransport equation is for only a single, solute species with possible linearequilibrium sorption and linear decay. Finite-difference techniques are used to discretize the governing equations using a point-distributed grid. The flow-, heatand solute-transport equations are solved, in turn, after a partial Gauss-reduction scheme is used to modify them. The modified equations are more tightly coupled and have better stability for the numerical solutions, The basic source-sink term represents wells. A complex well-flow model may be used to simulate specified flow rate and pressure conditions at the land surface or within the aquifer, with or without pressure and flow-rate constraints. Boundary-condition types offered include specified value, specified flux, leakage, heat conduction, an approximate free surface, and two types of aquifer-influence functions. All boundary conditions can be functions of time. Two techniques are available for solution of the finite-difference matrix equations. One technique is a direct-elimination solver, using equations reordered by alternating diagonal planes. The other technique is an iterative solver, using two-line successive overrelaxation. A restart option is available for storing intermediate results and restarting the simulation at an intermediate time with modified boundary conditions. This feature also can be used as protection against computer-system failure. Data input and output may be in metric (SI) units or inch-pound units. Output may include tables of dependent variables and parameters, zoned-contour maps, and plots of the dependent variables versus time. The HST3D program is a descendant of the Survey Waste Injection Program (SWIP) written for the U.S. Geological Survey under contract.

[1]  K. Aziz,et al.  Petroleum Reservoir Simulation , 1979 .

[2]  G. Taylor,et al.  The penetration of a fluid into a porous medium or Hele-Shaw cell containing a more viscous liquid , 1958, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[3]  R. H. Ritchie,et al.  Asymptotic Expansions of Solutions of the Heat Conduction Equation in Internally Bounded Cylindrical Geometry , 1956 .

[4]  G. D. Bennett,et al.  Representation of Multiaquifer Well Effects in Three-Dimensional Ground-Water Flow Simulation , 1982 .

[5]  K. G. Stollenwerk,et al.  Computer model of one-dimensional equilibrium controlled sorption processes , 1984 .

[6]  J. Stoer,et al.  Numerical treatment of ordinary differential equations by extrapolation methods , 1966 .

[7]  Mary P. Anderson,et al.  Using models to simulate the movement of contaminants through groundwater flow systems , 1979 .

[8]  N S Boulton,et al.  THE DRAWDOWN OF THE WATER-TABLE UNDER NON-STEADY CONDITIONS NEAR A PUMPED WELL IN AN UNCONFINED FORMATION. , 1954 .

[9]  G. W. Thomas,et al.  Principles of Hydrocarbon Reservoir Simulation , 1981 .

[10]  R. Banks,et al.  A solution of the differential equation of longitudinal dispersion in porous media , 1961 .

[11]  Some Practical Considerations in the Numerical Solution of Two-Dimensional Reservoir Problems , 1968 .

[12]  John E. Chappelear,et al.  Representing Wells in Numerical Reservoir Simulation: Part 2 - Implementation , 1981 .

[13]  J. S. Aronofsky,et al.  Mobility Ratio - Its Influence on Flood Patterns During Water Encroachment , 1952 .

[14]  J. C. Jaeger,et al.  Conduction of Heat in Solids , 1952 .

[15]  S. Lohman,et al.  Ground-Water Hydraulics , 1972 .

[16]  Mathematical Model of an Unstable Miscible Displacement , 1963 .

[17]  J. W. Mercer,et al.  Parameters and variables appearing in repository-siting models , 1982 .

[18]  B. Finlayson.  The method of weighted residuals and variational principles : with application in fluid mechanics, heat and mass transfer , 1972 .

[19]  R. B. Lantz.  Quantitative Evaluation of Numerical Diffusion (Truncation Error) , 1971 .

[20]  W. F. Stevens.  Applied numerical methods, b. comahan, h.a. luther, and j. o. wilkes, john wiley and sons, new york (1969) , 1970 .

[21]  F. B. Hildebrand,et al.  Introduction To Numerical Analysis , 1957 .

[22]  Thomas A. Prickett,et al.  Selected digital computer techniques for groundwater resource evaluation , 1971 .

[23]  J. Sproston.  Matrix computation for engineers and scientists , 1980 .

[24]  G. Pinder,et al.  Computational Methods in Subsurface Flow , 1983 .

[25]  George M. Hornberger,et al.  Numerical Methods in Subsurface Hydrology , 1972, Soil Science Society of America Journal.

[26]  H. J. Ramey.  Wellbore Heat Transmission , 1962 .

[27]  A. Scheidegger,et al.  THE STATISTICAL BEHAVIOR OF INSTABILITIES IN DISPLACEMENT PROCESSES IN POROUS MEDIA , 1961 .

[28]  John E. Chappelear,et al.  Representing Wells in Numerical Reservoir Simulation: Part 1 - Theory , 1981 .

[29]  E. A. Breitenbach,et al.  Treatment of Individual Wells and Grids in Reservoir Modeling , 1968 .

[30]  Leonard F. Konikow,et al.  Derivation of equations describing solute transport in ground water , 1977 .

[31]  J. Bear.  On the tensor form of dispersion in porous media , 1961 .

[32]  S. W. Lohman,et al.  Nonsteady flow to a well of constant drawdown in an extensive aquifer , 1952 .

[33]  T. Narasimhan,et al.  AN INTEGRATED FINITE DIFFERENCE METHOD FOR ANALYZING FLUID FLOW IN POROUS MEDIA , 1976 .

[34]  J. R. Fanchi.  Analytical representation of the van Everdingen-Hurst aquifer influence functions for reservoir simulation , 1985 .

[35]  Richard S. Varga,et al.  Application of Oscillation Matrices to Diffusion-Convection Equations , 1966 .

[36]  Khalid Aziz,et al.  Use of Irregular Grid in Reservoir Simulation , 1972 .

[37]  Irving H. Shames.  Mechanics of Fluids , 1962 .

[38]  W. Hurst,et al.  The Simplification of the Material Balance Formulas by the Laplace Transformation , 1958 .

[39]  A. Scheidegger.  Growth of Instabilities on Displacement Fronts in Porous Media , 1960 .

[40]  K. H. Coats,et al.  Direct methods in reservoir simulation , 1974 .

[41]  Olaf A. Hougen,et al.  Chemical process principles , 1947 .

[42]  J. W. Mercer,et al.  Theoretical analysis of fluid flow and energy transport in hydrothermal systems , 1977 .

[43]  J. M. Watt.  Numerical Initial Value Problems in Ordinary Differential Equations , 1972 .

[44]  R. Wooding,et al.  The stability of a viscous liquid in a vertical tube containing porous material , 1959, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.