This section compares accurate with inaccurate solution of subdomain problems.
We consider both the cylinder problem from the previous section and
a Poisseuille flow in a rectangle .
The global grid for the cylinder problem consists of 12240 grid cells.
The single-block cylinder grid in Figure 5 consists of
10800 grid cells.
For the Poisseuille problem a Cartesian grid of
cells
is used.
The subdomain problems are solved using GMRES with RILU()
preconditioning [50,48] and a relative
stopping criterion. For
we get the standard ILU
preconditioner [36]
and for
we get the Modified ILU preconditioner [27].
RILU(
) [1] lies in between these two.
With RILUD(
) we mean RILU(
) restricted to the diagonal.
The momentum equations are solved using
a RILUD(
preconditioner and the pressure equations using a
RILU(
). As a short-hand, we will use RILU(
) to mean
RILUD(
) whenever the momentum, and RILU(
) whenever the
pressure equations are involved.
The subdomain solution accuracy is varied. As a special
case the subdomain solution is approximated by means of the inverse of the
RILU(
) [20,50,48] preconditioner, thereby omitting
GMRES for subdomain solution.
The multi-block problem (the outer loop)
is solved up to a relative accuracy of .
In all experiments a Krylov space of dimension 20 is used for both
GMRES and GCR multi-block acceleration and for GMRES subdomain solution.
GMRES always uses a restart after 20 iterations. With GCR, we
investigate both optimized restart, denoted by GCR (restart), and truncation,
denoted by GCR (trunc).
Iteration counts and computing times are given in the tables in the
form time(iteration count). The iteration counts and times are summed over
all time steps taken (10 time steps are used in all examples).
The experiments are run on a HP9000/735 workstation.
In most of the experiments, the multiplicative algorithm is used.
Only section 6.4 examines the additive algorithm.
Section 6.1 examines the effect of lowering the accuracy of
the subdomain solution on the number of iterations and total computing time.
Section 6.2 compares single-block solution time with
multi-block solution time.
Section 6.3 examines the effect of the parameter
in the subdomain RILU(
) preconditioner on convergence using the
RIBLU(
) postconditioner.