download pdf version of PhD book - Universiteit Utrecht
download pdf version of PhD book - Universiteit Utrecht
download pdf version of PhD book - Universiteit Utrecht
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
6. Dispersivity under Partially-Saturated Conditions<br />
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .<br />
performed to determine pore-level distribution <strong>of</strong> each phase. Then, steadystate<br />
flow is established and equations <strong>of</strong> mass balance are solved to calculate<br />
transport properties <strong>of</strong> such a distribution and to obtain BTCs <strong>of</strong> solute concentration.<br />
The results <strong>of</strong> the modeling are also compared with experimental<br />
observations to show the capability <strong>of</strong> this formulation.<br />
The following morphological and computational features are introduced in the<br />
present study<br />
(1) The topology <strong>of</strong> the pore space is modeled using a MDPN which allows<br />
a distribution <strong>of</strong> coordination numbers ranging between one and 26.<br />
(2) To take into account the angularity <strong>of</strong> pores in natural porous media,<br />
pore throats with various cross sections, with a wide range <strong>of</strong> shape factor<br />
values and pore sizes, are used in the network. This includes rectangular,<br />
circular, and various irregular triangular cross sections.<br />
(3) The pore body size distributions are assumed to follow a truncated lognormal<br />
distribution, without any correlation. The pore-throat size distributions<br />
are related to the pore body size distributions.<br />
(4) Both pore bodies and pore throats are considered to have volume. This<br />
means we solve mass balance equations and calculate solute concentrations<br />
and mass fluxes within both pore bodies and pore throats.<br />
(5) As soon as a pore body is (partially) saturated, it will be discretized<br />
further into smaller regions occupied by water, each with its own flow<br />
rate and concentration, in order to capture the effect <strong>of</strong> limited mixing<br />
due to the partial filling <strong>of</strong> the pore.<br />
(6) Upon invasion <strong>of</strong> a pore throat by the non-wetting phase, each edge <strong>of</strong><br />
the pore throat will be considered as a separate domain with its own flow<br />
rate and concentration.<br />
(7) Employing a fully implicit numerical scheme to calculate the unknown<br />
connections, a substitution method is introduced which considerably reduces<br />
the computational time.<br />
(8) Various parameters and relations, including coefficient <strong>of</strong> variation <strong>of</strong> the<br />
velocities field, relative permeability-saturation (k r − S w ) curves, capillary<br />
pressure-saturation (P c − S w ) curves, and fraction <strong>of</strong> percolating<br />
saturated pores are also computed.<br />
132