UMSI 2000 Annual Report: John L. Nieber, Associate Fellow Previous Page  |  Table of Contents  |  Next Page

John L. Nieber, Associate Fellow


Simulation of Flow and Solute Transport and Dispersion in Heterogeneous Porous Media

This work involves simulating the processes of flow and transport in porous media with consideration for geochemical and biochemical reactions. The research is being performed in cooperation with researchers at the Engineering Research and Development Center in Vicksburg, Mississippi.

Research Group

Abdelkarim Abulaban, Research Associate
Ronnie Daanen, Graduate Student Researcher
Anita Gruber, Graduate Student Researcher
Debasmita Misra, Research Associate
Cam Nguyen, Research Associate
Hung Viet Nguyen, Graduate Student Researcher
Paul Oduro, Graduate Student Researcher

Two specific problems are being addressed. The first deals with the characterization of the dispersion process based on the various levels of spatial heterogeneity that exist in natural porous media. Solute transport is simulated using a particle method in which the solute mass is represented by a very large ensemble of discrete solute parcels (or walkers) that move in the flow field. Both a deterministic step and a random step are imposed to define the parcel pathway. Applying this method, these researchers are characterizing the velocity covariance function (and apparent dispersion coefficient) that results from the transport of a conservative (non-reactive) solute within a steady-state heterogeneous flow field. Relating the velocity covariance function (and apparent dispersion coefficient) to measures of heterogeneity will be of value to practitioners. From this work, a low-resolution model is being developed to simulate solute transport at a much smaller computational effort than required for high-resolution simulations. The low-resolution model will be parameterized using the velocity covariance data derived from the high-resolution simulations.


Sample concentration breakthrough curves for different combinations of the adsorption-desorption isotherms
For the second problem, the processes of geochemical reactions and biochemical reactions have been incorporated into the simulation of solute transport. Again, the particle method is used to simulate the solute migration. Fairly simple, but quite significant process of equilibrium sorption is being examined for aqueous solutes to porous media solids. Of particular interest is the character of the equilibrium sorption isotherm. When this isotherm is nonlinear, the transport of a solute is affected significantly. We have also considered the condition where adsorption and desorption are not given by the same isotherm. This sorption hysteresis has been found to dramatically affect the movement and spread of solute plumes in saturated porous media.

In both of these problems, it is necessary to perform simulations on highly resolved computational domains, with upwards of 106 grid cells to define the flow field for two-dimensional domains. Three-dimensional flow fields will have much larger requirements. According to the theory for the particle method, the larger the number of parcels, the more accurate the solution. Exact solutions require an infinite number of particles, but in practice the number of particles can be much less, depending on the particular problem being simulated. These researchers have found that a particle number on the order of 106 is generally sufficient. The requirements to perform computations on these highly resolved grids with large numbers of particles naturally leads to the use of the computational resources available at the Supercomputing Institute.


Previous Page  |  Table of Contents  |  Next Page