Hoaglund John R, Pollard David
The Pennsylvania State University, Earth and Mineral Sciences Environment Institute, 2217 Earth-Engineering Sciences, University Park, PA 16802-6813, USA.
Ground Water. 2003 Nov-Dec;41(6):841-6. doi: 10.1111/j.1745-6584.2003.tb02425.x.
Darcy flow equations relating vertical and bedding-parallel flow to vertical and bedding-parallel gradient components are derived for a skewed Cartesian grid in a vertical plane, correcting for structural dip given the principal hydraulic conductivities in bedding-parallel and bedding-orthogonal directions. Incorrect-minus-correct flow error results are presented for ranges of structural dip (0 < or = theta < or = 90) and gradient directions (0 < or = phi < or = 360). The equations can be coded into ground water models (e.g., MODFLOW) that can use a skewed Cartesian coordinate system to simulate flow in structural terrain with deformed bedding planes. Models modified with these equations will require input arrays of strike and dip, and a solver that can handle off-diagonal hydraulic conductivity terms.