The partial differential equation that describes transient flow through a saturated anisotropic medium (Section 2.11) is

(A9.1)

For a horizontal, confined aquifer of thickness, , the two-dimensional form of Eq. (A9.1) reduces to

(A9.2)

where and are the principal components of the transmissivity tensor defined by , and is the storativity defined by . To find the finite-difference equation for an interior node in the nodal grid used to discretize the region of flow, we must replace the partial derivatives in Eq. (A9.2) by differences. Using the definitions developed in Appendix VI and the notation of Section 8.8 and Figure 8.26(c), we can write the finite-difference expression:

(A9.3)

where the subscript means that the bracketed quantity is evaluated at the midpoint between nodes and , and the superscript means that the bracketed quantity is evaluated at time step . We can further approximate terms on the right-hand side of Eq. (A9.3) by

(A9.4A)

(A9.4B)

If we evaluate and by simple averages of the form

(A9.5)

then these expressions can be substituted in Eqs. (A9.4), and Eqs. (A9.4) can in turn be substituted in Eq. (A9.3) to give

(A9.6)

Similarly,

(A9.7)

Finally, we can approximate the right-hand side of Eq. (A9.2) by

(A9.8)

Substituting Eqs. (A9.6), (A9.7), and (A9.8) for the three terms in Eqs. (A9.6) and gathering terms leads to the general finite-difference equation for an internal node in a heterogeneous, anisotropic aquifer:

(A9.9)

where

If the aquifer is homogeneous and isotropic, then for all and for all . Under these conditions, and for a square nodal grid with , the coefficients of Eq. (A9.9) become

If we divide through by , these coefficients are seen to be the same as those developed in a less rigorous way in Section 8.8 and presented in connection with Eq. (8.57).

Trescott et al. (1976) have suggested that there is some advantage to utilizing the harmonic mean rather than the arithmetic mean in Eq. (A9.5). This approach changes the coefficients in the finite-difference equation but does not alter the concepts underlying the development.

Equation (A9.9) is written in terms of the hydraulic head values at five nodes at time step and one node at time step . It is known as a backward-difference approximation. Remson et al. (1971) note that there are some computational advantages to using a central-difference approximation, known as the *Crank-Nicholson scheme*, that utilizes head values at five nodes at time step and five nodes at time step . The alternating-direction implicit procedure (ADIP) utilized by Pinder and Bredehoeft (1968) involves two finite-difference equations, one in the plane and one in the plane. Each uses head values at three nodes at time step and three nodes at time step .