A three-dimensional, three-phase numerical model is presented for simulating the movement of immiscible fluids, including
nonaqueous-phase liquids (NAPLs), through porous media. The model is designed to simulate soil flume experiments and for practical
application to a wide variety of contamination scenarios involving light or dense NAPLs in heterogeneous subsurface systems.
The model is derived for the three-phase flow of water, NAPL, and air in porous media. The basic governing equations are based
upon the mass conservation of the constitutents within the phases. The descretization chosen to transform the governing equations
into the approximating equations, although logically regular, is very general. The approximating equations are a set of simultaneous
coupled nonlinear equations which are solved by the Newton-Raphson method. The linear system solutions needed for the Newton-Raphson
method are obtained using a matrix of preconditioner/accelerator iterative methods.
Because of the special way the governing equations are implemented, the model is capable of simulating many of the phenomena
considered necessary for the sucessful simulation of field problems including entry pressure phenomena, entrapment, and preferential
flow paths. The model is verified by comparing it with several exact analytic test solutions and three soil flume experiments
involving the introduction and movement of light nonaqueous-phase liquid (LNAPL) or dense nonaqueous-phase liquid (DNAPL)
in heterogeneous sand containing a watertable.
This revised version was published online in August 2006 with corrections to the Cover Date. 相似文献
The solubility of Ti- and P-rich accessory minerals has been examined as a function of pressure and K2O/Na2O ratio in two series of highly evolved silicate systems. These systems correspond to (a) alkaline, varying from alkaline to peralkaline with increasing K2O/Na2O ratio; and (b) strongly metaluminous (essentially trondhjemitic at the lowest K2O/Na2O ratio) and remaining metaluminous with increasing K2O/Na2O ratio (to 3). The experiments were conducted at a fixed temperature of 1000 °C, with water contents varying from 5 wt.% at low pressure (0.5 GPa), increasing through 5–10 wt.% at 1.5–2.5 GPa to 10 wt.% at 3.5 GPa. Pressure was extended outside the normal crustal range, so that the results may also be applied to derivation of hydrous silicic melts from subducted oceanic crust.
For the alkaline composition series, the TiO2 content of the melt at Ti-rich mineral saturation decreases with increasing pressure but is unchanged with increasing K content (at fixed pressure). The P2O5 content of the alkaline melts at apatite saturation increases with increased pressure at 3.5 GPa only, but decreases with increasing K content (and peralkalinity). For the metaluminous composition series (termed as “trondhjemite-based series” (T series)), the TiO2 content of the melt at Ti-rich mineral saturation decreases with increasing pressure and with increasing K content (at fixed pressure). The P2O5 content of the T series melts at apatite saturation is unchanged with increasing pressure, but decreases with increasing K content. The contrasting results for P and Ti saturation levels, as a function of pressure in both compositions, point to contrasting behaviour of Ti and P in the structure of evolved silicate melts. Ti content at Ti-rich mineral saturation is lower in the alkaline compared with the T series at 0.5 GPa, but is similar at higher pressures, whereas P content at apatite saturation is lower in the T series at all pressures studied. The results have application to A-type granite suites that are alkaline to peralkaline, and to I-type metaluminous suites that frequently exhibit differing K2O/Na2O ratios from one suite to another. 相似文献