Sensitivity of Calculated Hydrological Flows Through Multilayered Hard Rock to Computational Solution Procedures

摘要

Permeability and moisture content curves for partially saturated, fractured, welded tuffs, such as those found near the site of a prospective nuclear waste repository at Yucca Mountain, Nevada, are highly nonlinear. We illustrate, by means of one-dimensional infiltration in multilayered fractured tuff, the numerical instabilities that can arise when analyzing the flow in such porous materials. We compare the results obtained using two numerical schemes for treating the nonlinearities: the first uses a first-order Piccard iteration scheme to attain convergence at each time step; the second uses a second-order Newton-Raphson iteration procedure to attain convergence. The Piccard scheme converges when flow through the rock matrix dominates but does not converge when flow through fractures dominates. On the other hand, the Newton-Raphson scheme converges in both situations. 6 refs., 10 figs. (ERA citation 10:038890)