Application of a lumped-process mathematical model to dissolution of non-uniformly distributed immiscible liquid in heterogeneous porous media

J. C. Marble, E. L. DiFilippo, Z. Zhang, G. R. Tick, Mark L Brusseau

Research output: Contribution to journalArticle

11 Citations (Scopus)

Abstract

The use of a lumped-process mathematical model to simulate the complete dissolution of immiscible liquid non-uniformly distributed in physically heterogeneous porous-media systems was investigated. The study focused specifically on systems wherein immiscible liquid was poorly accessible to flowing water. Two representative, idealized scenarios were examined, one wherein immiscible liquid at residual saturation exists within a lower-permeability unit residing in a higher-permeability matrix, and one wherein immiscible liquid at higher saturation (a pool) exists within a higher-permeability unit adjacent to a lower-permeability unit. As expected, effluent concentrations were significantly less than aqueous solubility due to dilution and by-pass flow effects. The measured data were simulated with two mathematical models, one based on a simple description of the system and one based on a more complex description. The permeability field and the distribution of the immiscible-liquid zones were represented explicitly in the more complex, distributed-process model. The dissolution rate coefficient in this case represents only the impact of local-scale (and smaller) processes on dissolution, and the parameter values were accordingly obtained from the results of experiments conducted with one-dimensional, homogeneously-packed columns. In contrast, the system was conceptualized as a pseudo-homogeneous medium with immiscible liquid uniformly distributed throughout the system for the simpler, lumped-process model. With this approach, all factors that influence immiscible-liquid dissolution are incorporated into the calibrated dissolution rate coefficient, which in such cases serves as a composite or lumped term. The calibrated dissolution rate coefficients obtained from the simulations conducted with the lumped-process model were approximately two to three orders-of-magnitude smaller than the independently-determined values used for the simulations conducted with the distributed-process model. This disparity reflects the difference in implicit versus explicit consideration of the larger-scale factors influencing immiscible-liquid dissolution in the systems.

Original languageEnglish (US)
Pages (from-to)1-10
Number of pages10
JournalJournal of Contaminant Hydrology
Volume100
Issue number1-2
DOIs
StatePublished - Aug 20 2008

Fingerprint

Porous materials
porous medium
Dissolution
dissolution
Mathematical models
liquid
Liquids
permeability
saturation
homogeneous medium
Dilution
simulation
Effluents
solubility
dilution
Solubility
effluent
matrix
Water
Composite materials

Keywords

  • Dissolution
  • DNAPL
  • Modeling

ASJC Scopus subject areas

  • Earth-Surface Processes

Cite this

Application of a lumped-process mathematical model to dissolution of non-uniformly distributed immiscible liquid in heterogeneous porous media. / Marble, J. C.; DiFilippo, E. L.; Zhang, Z.; Tick, G. R.; Brusseau, Mark L.

In: Journal of Contaminant Hydrology, Vol. 100, No. 1-2, 20.08.2008, p. 1-10.

Research output: Contribution to journalArticle

@article{decb0ad22de1414485155fd991fc86a0,
title = "Application of a lumped-process mathematical model to dissolution of non-uniformly distributed immiscible liquid in heterogeneous porous media",
abstract = "The use of a lumped-process mathematical model to simulate the complete dissolution of immiscible liquid non-uniformly distributed in physically heterogeneous porous-media systems was investigated. The study focused specifically on systems wherein immiscible liquid was poorly accessible to flowing water. Two representative, idealized scenarios were examined, one wherein immiscible liquid at residual saturation exists within a lower-permeability unit residing in a higher-permeability matrix, and one wherein immiscible liquid at higher saturation (a pool) exists within a higher-permeability unit adjacent to a lower-permeability unit. As expected, effluent concentrations were significantly less than aqueous solubility due to dilution and by-pass flow effects. The measured data were simulated with two mathematical models, one based on a simple description of the system and one based on a more complex description. The permeability field and the distribution of the immiscible-liquid zones were represented explicitly in the more complex, distributed-process model. The dissolution rate coefficient in this case represents only the impact of local-scale (and smaller) processes on dissolution, and the parameter values were accordingly obtained from the results of experiments conducted with one-dimensional, homogeneously-packed columns. In contrast, the system was conceptualized as a pseudo-homogeneous medium with immiscible liquid uniformly distributed throughout the system for the simpler, lumped-process model. With this approach, all factors that influence immiscible-liquid dissolution are incorporated into the calibrated dissolution rate coefficient, which in such cases serves as a composite or lumped term. The calibrated dissolution rate coefficients obtained from the simulations conducted with the lumped-process model were approximately two to three orders-of-magnitude smaller than the independently-determined values used for the simulations conducted with the distributed-process model. This disparity reflects the difference in implicit versus explicit consideration of the larger-scale factors influencing immiscible-liquid dissolution in the systems.",
keywords = "Dissolution, DNAPL, Modeling",
author = "Marble, {J. C.} and DiFilippo, {E. L.} and Z. Zhang and Tick, {G. R.} and Brusseau, {Mark L}",
year = "2008",
month = "8",
day = "20",
doi = "10.1016/j.jconhyd.2008.04.003",
language = "English (US)",
volume = "100",
pages = "1--10",
journal = "Journal of Contaminant Hydrology",
issn = "0169-7722",
publisher = "Elsevier",
number = "1-2",

}

TY - JOUR

T1 - Application of a lumped-process mathematical model to dissolution of non-uniformly distributed immiscible liquid in heterogeneous porous media

AU - Marble, J. C.

AU - DiFilippo, E. L.

AU - Zhang, Z.

AU - Tick, G. R.

AU - Brusseau, Mark L

PY - 2008/8/20

Y1 - 2008/8/20

N2 - The use of a lumped-process mathematical model to simulate the complete dissolution of immiscible liquid non-uniformly distributed in physically heterogeneous porous-media systems was investigated. The study focused specifically on systems wherein immiscible liquid was poorly accessible to flowing water. Two representative, idealized scenarios were examined, one wherein immiscible liquid at residual saturation exists within a lower-permeability unit residing in a higher-permeability matrix, and one wherein immiscible liquid at higher saturation (a pool) exists within a higher-permeability unit adjacent to a lower-permeability unit. As expected, effluent concentrations were significantly less than aqueous solubility due to dilution and by-pass flow effects. The measured data were simulated with two mathematical models, one based on a simple description of the system and one based on a more complex description. The permeability field and the distribution of the immiscible-liquid zones were represented explicitly in the more complex, distributed-process model. The dissolution rate coefficient in this case represents only the impact of local-scale (and smaller) processes on dissolution, and the parameter values were accordingly obtained from the results of experiments conducted with one-dimensional, homogeneously-packed columns. In contrast, the system was conceptualized as a pseudo-homogeneous medium with immiscible liquid uniformly distributed throughout the system for the simpler, lumped-process model. With this approach, all factors that influence immiscible-liquid dissolution are incorporated into the calibrated dissolution rate coefficient, which in such cases serves as a composite or lumped term. The calibrated dissolution rate coefficients obtained from the simulations conducted with the lumped-process model were approximately two to three orders-of-magnitude smaller than the independently-determined values used for the simulations conducted with the distributed-process model. This disparity reflects the difference in implicit versus explicit consideration of the larger-scale factors influencing immiscible-liquid dissolution in the systems.

AB - The use of a lumped-process mathematical model to simulate the complete dissolution of immiscible liquid non-uniformly distributed in physically heterogeneous porous-media systems was investigated. The study focused specifically on systems wherein immiscible liquid was poorly accessible to flowing water. Two representative, idealized scenarios were examined, one wherein immiscible liquid at residual saturation exists within a lower-permeability unit residing in a higher-permeability matrix, and one wherein immiscible liquid at higher saturation (a pool) exists within a higher-permeability unit adjacent to a lower-permeability unit. As expected, effluent concentrations were significantly less than aqueous solubility due to dilution and by-pass flow effects. The measured data were simulated with two mathematical models, one based on a simple description of the system and one based on a more complex description. The permeability field and the distribution of the immiscible-liquid zones were represented explicitly in the more complex, distributed-process model. The dissolution rate coefficient in this case represents only the impact of local-scale (and smaller) processes on dissolution, and the parameter values were accordingly obtained from the results of experiments conducted with one-dimensional, homogeneously-packed columns. In contrast, the system was conceptualized as a pseudo-homogeneous medium with immiscible liquid uniformly distributed throughout the system for the simpler, lumped-process model. With this approach, all factors that influence immiscible-liquid dissolution are incorporated into the calibrated dissolution rate coefficient, which in such cases serves as a composite or lumped term. The calibrated dissolution rate coefficients obtained from the simulations conducted with the lumped-process model were approximately two to three orders-of-magnitude smaller than the independently-determined values used for the simulations conducted with the distributed-process model. This disparity reflects the difference in implicit versus explicit consideration of the larger-scale factors influencing immiscible-liquid dissolution in the systems.

KW - Dissolution

KW - DNAPL

KW - Modeling

UR - http://www.scopus.com/inward/record.url?scp=48549100660&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=48549100660&partnerID=8YFLogxK

U2 - 10.1016/j.jconhyd.2008.04.003

DO - 10.1016/j.jconhyd.2008.04.003

M3 - Article

C2 - 18555558

AN - SCOPUS:48549100660

VL - 100

SP - 1

EP - 10

JO - Journal of Contaminant Hydrology

JF - Journal of Contaminant Hydrology

SN - 0169-7722

IS - 1-2

ER -