We use a preconditioned conjugate-gradient method [7] to solve our discretized diffusion equations. The preconditioner is completely analogous to that used for the 2-D local support-operators scheme. [1] It is obtained simply by setting the off-diagonal elements of all the S -matrices, defined by Eq. (19), to zero. This makes the W -matrices diagonal, which ultimately results in a 7-point pure cell-center diffusion equation. That is to say that the balance equation for each cell in the mesh interior contains only the cell-center intensity in that cell together with the cell-center intensities in the six cells sharing a face with that cell. Each face-center intensity is expressed in terms of the cell-center intensities in the two cells that share that face-center intensity. For instance, if we set the off-diagonal elements of the S -matrices to zero, Eqs. (26) and (27) yield:
and respectively, where= 2s_{L, L}^{LBD}V^{LBD} + s_{L, L}^{LTD}V^{LTD} + s_{L, L}^{LBU}V^{LBU} + s_{L, L}^{LTU}V^{LTU} , | (57) |
= 2s_{R, R}^{RBD}V^{RBD} + s_{R, R}^{RTD}V^{RTD} + s_{R, R}^{RBU}V^{RBU} + s_{R, R}^{RTU}V^{RTU} . | (58) |
To obtain the analog of Eq. (61) for a cell face on the outer mesh boundary, we again consider a cell (1, j, k) , whose left face is on the boundary with its other faces in the mesh interior. Substituting from Eqs. (53) and (55) into Eq. (54), we obtain the equation for :
where Solving Eq. (65) for . we get Substituting from Eq. (67) into Eqs. (53) and (55), we obtain the desired expression for the face-area flux component on a boundary face: where and is given by Eq. (63) evaluated with i = 0 and Eq. (66). This completes the derivation of the approximate cell-center diffusion scheme used to precondition the full cell-center/cell-face scheme.Once the 7-point discretization equations have been solved for the cell-center intensities, the face-center intensities can be directly calculated. In particular, Eq. (60) and its analogs for the Bottom/Top and Down/Up faces are used to calculate the face-center intensities on the mesh interior, while Eq. (65) and its analogs for the Bottom/Top and Down/Up faces are used to calculate the face-center intensities on the outer mesh boundary. The 7-point cell-center discretization yields an SPD coefficient matrix. Most importantly, since the S -matrices are diagonal when the mesh is orthogonal, it follows that the 7-point cell-center discretization is equivalent to the full cell-center/face-center discretization whenever the mesh is orthogonal. Thus our preconditioner can be expected to be very effective if the mesh is not too skewed. Our preconditioning system costs much less to solve than the full system because the coefficient matrix of the preconditioning system has roughly one-fourth as many rows and one-sixth as many elements as the full-system coefficient matrix. Computational results presented in the next section confirm this expectation.
In closing this section we note that when the 7-point system is used for preconditioning purposes, an inhomogeneous source term will generally appear in both the cell-center and face-center intensity equations. We did not include such a source in our derivation of the face-center intensity equations because they do not appear in standard calculations. One must remember to include these sources before the face-center intensities are eliminated to obtain the 7-point cell-center system. This matter is extensively discussed for the 2-D case in [1].