Application of a nonlinear slug test model.

C D McElwee

doi:10.1111/j.1745-6584.2001.tb02364.x

Home/Archive/PubMed

PubMed · 11554252

Application of a nonlinear slug test model.

C D McElwee

View original source ↗Export metadata ↓

Abstract

Knowledge of the hydraulic conductivity distribution is of utmost importance in understanding the dynamics of an aquifer and in planning the consequences of any action taken upon that aquifer. Slug tests have been used extensively to measure hydraulic conductivity in the last 50 years since Hvorslev's (1951) work. A general nonlinear model based on the Navier-Stokes equation, nonlinear frictional loss, non-Darcian flow, acceleration effects, radius changes in the wellbore, and a Hvorslev model for the aquifer has been implemented in this work. The nonlinear model has three parameters: beta, which is related primarily to radius changes in the water column; A, which is related to the nonlinear head losses; and K, the hydraulic conductivity. An additional parameter has been added representing the initial velocity of the water column at slug initiation and is incorporated into an analytical solution to generate the first time step before a sequential numerical solution generates the remainder of the time solution. Corrections are made to the model output for acceleration before it is compared to the experimental data. Sensitivity analysis and least squares fitting are used to estimate the aquifer parameters and produce some diagnostic results, which indicate the accuracy of the fit. Finally, an example of field data has been presented to illustrate the application of the model to data sets that exhibit nonlinear behavior. Multiple slug tests should be taken at a given location to test for nonlinear effects and to determine repeatability.

Explore related subjects

Models, Theoretical Pressure Soil Water Water Movements Journal Article

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

C D McElwee. Application of a nonlinear slug test model.. https://doi.org/10.1111/j.1745-6584.2001.tb02364.x

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

PubMed2007

Sensitivity of the solution of the Elder problem to density, velocity and numerical perturbations.

In this paper the Elder problem is studied with the purpose of evaluating the inherent instabilities associated with the numerical solution of this problem. Our focus is first on the question of the existence of a unique numerical solution for this problem, and second on the grid density and fluid density requirements necessary for a unique numerical solution. In particular we have investigated the instability issues associated with the numerical solution of the Elder problem from the following perspectives: (i) physical instability issues associated with density differences; (ii) sensitivity of the numerical solution to idealization irregularities; and, (iii) the importance of a precise velocity field calculation and the association of this process with the grid density levels that is necessary to solve the Elder problem accurately. In the study discussed here we have used a finite element Galerkin model we have developed for solving density-dependent flow and transport problems, which will be identified as TechFlow. In our study, the numerical results of Frolkovic and de Schepper [Frolkovic, P. and H. de Schepper, 2001. Numerical modeling of convection dominated transport coupled with density-driven flow in porous media, Adv. Water Resour., 24, 63-72.] were replicated using the grid density employed in their work. We were also successful in duplicating the same result with a less dense grid but with more computational effort based on a global velocity estimation process we have adopted. Our results indicate that the global velocity estimation approach recommended by Yeh [Yeh, G.-T., 1981. On the computation of Darcian velocity and mass balance in finite element modelling of groundwater flow, Water Resour. Res., 17(5), 1529-1534.] allows the use of less dense grids while obtaining the same accuracy that can be achieved with denser grids. We have also observed that the regularity of the elements in the discretization of the solution domain does make a difference in obtaining a unique stationary solution for this problem. The results of our study also indicate that the density differences are critical in the solution of the Elder problem and that high density differences lead to the physical instability that is inherent with this problem. Other than the physical instability associated with the level of density differences used in the Elder problem, the following two points should be considered in solving the Elder problem in a consistent manner: (i) strict attention should be paid to the vertical grid Peclet number in developing the criteria for convergent grid selection; and, (ii) with a globally continuous velocity calculation stable solutions can be obtained at lower grid densities.

Models, Theoretical↗

PubMed2007

Modular and dynamic functionalization of polymeric scaffolds.

The design and synthesis of multifunctionalized, architecturally controlled polymers is a prerequisite for a variety of future applications of polymeric materials. On the basis of Nature's use of self-assembly in the creation of biomaterials, this Account describes concepts that were developed over the past 5 years that utilize noncovalent interactions such as hydrogen bonding, ionic interactions, electrostatic interactions, metal coordination, and pi-pi stacking in modification of copolymer side-chains to obtain multifunctional polymeric materials, induce polymer morphology changes, and influence bulk-polymer properties.

Models, Theoretical↗

PubMed2006

Eigenvalue spectra of random matrices for neural networks.

The dynamics of neural networks is influenced strongly by the spectrum of eigenvalues of the matrix describing their synaptic connectivity. In large networks, elements of the synaptic connectivity matrix can be chosen randomly from appropriate distributions, making results from random matrix theory highly relevant. Unfortunately, classic results on the eigenvalue spectra of random matrices do not apply to synaptic connectivity matrices because of the constraint that individual neurons are either excitatory or inhibitory. Therefore, we compute eigenvalue spectra of large random matrices with excitatory and inhibitory columns drawn from distributions with different means and equal or different variances.

Models, Theoretical↗