Search USGS⌕ Search

Geology topics

L.J. Torak

Publications and source records attributed to L.J. Torak.

4 recordsLinked to original sources
USGS1993

A modular finite-element model (MODFE) for areal and axisymmetric ground-water-flow problems, Part 3: Design philosophy and programming details

A MODular Finite-Element, digital-computer program (MODFE) was developed to simulate steady or unsteady-state, two-dimensional or axisymmetric ground-water-flow. The modular structure of MODFE places the computationally independent tasks that are performed routinely by digital-computer programs simulating ground-water flow into separate subroutines, which are executed from the main program by control statements. Each subroutine consists of complete sets of computations, or modules, which are identified by comment statements, and can be modified by the user without affecting unrelated computations elsewhere in the program. Simulation capabilities can be added or modified by either adding or modifying subroutines that perform specific computational tasks, and the modular-program structure allows the user to create versions of MODFE that contain only the simulation capabilities that pertain to the ground-water problem of interest. MODFE is written in a Fortran programming language that makes it virtually device independent and compatible with desk-top personal computers and large mainframes. MODFE uses computer storage and execution time efficiently by taking advantage of symmetry and sparseness within the coefficient matrices of the finite-element equations. Parts of the matrix coefficients are computed and stored as single-subscripted variables, which are assembled into a complete coefficient just prior to solution. Computer storage is reused during simulation to decrease storage requirements. Descriptions of subroutines that execute the computational steps of the modular-program structure are given in tables that cross reference the subroutines with particular versions of MODFE. Programming details of linear and nonlinear hydrologic terms are provided. Structure diagrams for the main programs show the order in which subroutines are executed for each version and illustrate some of the linear and nonlinear versions of MODFE that are possible. Computational aspects of changing stresses and boundary conditions with time and of mass-balance and error terms are given for each hydrologic feature. Program variables are listed and defined according to their occurrence in the main programs and in subroutines. Listings of the main programs and subroutines are given.

Techniques of Water-Resources Investigations↗
USGS1993

Geohydrology and evaluation of water-resource potential of the upper Floridan Aquifer in the Albany area, southwestern Georgia

In the Albany area of southwestern Georgia, the Upper Floridan aquifer lies entirely within the Dougherty Plain district of the Coastal Plain physiographic province, and consists of the Ocala Limestone of late Eocene age. The aquifer is divided throughout most of the study area into an upper and a lower lithologic unit, which creates an upper and a lower water-bearing zone. The lower waterbearing zone consists of alternating layers of sandy limestone and medium-brown, recrystallized dolomitic limestone, and ranges in thickness from about 50 ft to 100 ft. It is highly fractured and exhibits well-developed permeability by solution features that are responsible for transmitting most of the ground water in the aquifer. Transmissivity of the lower water-bearing zone ranges from about 90,000 to 178,000 ft2/d. The upper water-bearing zone is a finely crystallized-to-oolitic, locally dolomitic limestone having an average thickness of about 60 ft. Transmissivities are considerably less in the upper water-bearing zone than in the lower water-bearing zone. The Upper Floridan aquifer is overlain by about 20-120 ft of undifferentiated overburden consisting of fine-to-coarse quartz sand and noncalcareous clay. A clay zone about 10-30 ft thick may be continuous throughout the southwestern part of the Albany area and, where present, causes confinement of the Upper Floridan aquifer and creates perched ground water after periods of heavy rainfall. The Upper Floridan aquifer is confined below by the Lisbon Formation, a mostly dolomitic limestone that contains trace amounts of glauconite. The Lisbon Formation is at least 50 ft thick in the study area and acts as an impermeable base to the Upper Floridan aquifer. The quality of ground water in the Upper Floridan aquifer is suitable for most uses; wells generally yield water of the hard, calcium-bicarbonate type that meets the U.S. Environmental Protection Agency's Primary or Secondary Drinking-Water Regulations. The water-resource potential of the Upper Floridan aquifer was evaluated by compiling results of drilling and aquifer testing in the study area, and by conducting computer simulations of the ground-water flow system under the seasonally low conditions of November 1985, and under conditions of pumping within a 12-mi 2 area located southwest of Albany. Results of test drilling, aquifer testing, and water-quality analyses indicate that, in the area southwest of Albany, geohydrologic conditions in the Upper Floridan aquifer, undifferentiated overburden, and Lisbon Formation were favorable for the aquifer to provide a large quantity of water without having adverse effects on the groundwater system. The confinement of the Upper Floridan aquifer by the undifferentiated overburden and the rural setting of the area of potential development decrease the likelihood that chemical constituents will enter the aquifer during development of the ground-water resources. Computer simulations of ground-water flow in the Upper Floridan aquifer, incorporating conditions for regional flow across model boundaries, leakage from rivers and other surface-water features, and vertical leakage from the undifferentiated overburden, were conducted by using a finite-element model for ground-water flow in two dimensions. Comparison of computed and measured water levels in the Upper Floridan aquifer for November 1985 at 74 locations indicated that computed water levels generally were within 5 ft of the measured values, which is the accuracy to which measured water levels were known. Water-level altitudes ranged from about 260 ft to 130 ft above sea level in the study area during calibration. Aquifer discharge to the Flint River downstream from the Lake Worth dam was computed by the calibrated model to be about 1 billion gallons per day; about 300 million gallons per day (Mgal/d) greater than was measured for similar lowflow conditions. T

Georgia↗
USGS1982

Modifications and corrections to the finite-difference model for simulation of three-dimensional ground-water flow

This report describes modifications incorporated into the finite-difference model for simulation of groundwater flow in three dimensions. These modifications extend the application of this model to simulations involving head-dependent sources and sinks (i.e., rivers, evapotranspiration, and springs or drains). Other modifications are made that enhance the iterative-solution process of the Strongly Implicit Procedure (SIP). An acceleration (or dampening) factor is introduced to the matrix equation that is to be solved, and optional methods of computing iteration parameters are incorporated into the original model. This report also describes corrections to the model that eliminate errors in the equation formulation and in mass-balance computations of certain types of simulations. Additional data-input instructions, definitions of new variables, and a list of program-statement changes are given in an appendix. (USGS)

Water-Resources Investigations Report↗