Massively parallel computers will play an increasingly dominant role in hydrological computing. One such computer is the Connection Machine model CM-2, a single-instruction stream, multiple-data stream computer with up to 65,536 processors, as much as 8 gigabytes (Gbyte) of random access memory distributed among the processors, and a FORTRAN compiler based on the proposed FORTRAN-90 standard. One-, two-, and three-dimensional examples from hydrology are used in this paper to present a tutorial on programming for the CM-2. The problem of saturated, steady flow in a randomly heterogeneous three-dimensional porous medium is explored here in some detail. A diagonally preconditioned conjugate gradient (DPCG( iterative solver is applied to this problem for up to 1283 nodes. Less than l min of CM-2 time is required to reduce the error by a factor of 10−6 for a 128 × 128 × 128 grid with heterogeneous hydraulic conductivity. Measured CPU times for the DPCG method are significantly smaller than those reported in the literature for a polynomial PCG solver applied to the same domain with different boundary conditions and executed on a Cray X-MP/48 and an Alliant FX/8. The measured performance is also much greater than that reported in the literature for a strongly implicit procedure solver applied to a similar problem on a Cray 2. The need for continued development of massively parallel algorithms, including effective iterative solution of linear systems of equations and problems with irregular domains, is indicated.
[1]
John P. Hayes,et al.
Computer architecture and organization; (2nd ed.)
,
1988
.
[2]
R W Hockney,et al.
Computer Simulation Using Particles
,
1966
.
[3]
Gehan A. J. Amaratunga,et al.
Parallel computational techniques for simulating dopant diffusion in silicon
,
1989
.
[4]
Ted Belytschko,et al.
Finite element analysis on the connection machine
,
1990
.
[5]
S. Ashby,et al.
A numerical investigation of the conjugate gradient method as applied to three‐dimensional groundwater flow problems in randomly heterogeneous porous media
,
1989
.
[6]
Giuseppe Gambolati,et al.
Is a simple diagonal scaling the best preconditioner for conjugate gradients on supercomputers
,
1990
.
[7]
D. E. Dougherty,et al.
Particle-grid methods for reacting flows in porous media with application to Fisher's equation
,
1992
.
[8]
John von Neumann,et al.
Theory Of Self Reproducing Automata
,
1967
.
[9]
David M. Young,et al.
A historical overview of iterative methods
,
1989
.
[10]
G. Birkhoff,et al.
Numerical Solution of Elliptic Problems
,
1984
.
[11]
S. Lennart Johnsson,et al.
Data structures and algorithms for the finite element method on a data parallel supercomputer
,
1990
.
[12]
Bruce M. Boghosian,et al.
Computational physics on the Connection Machine
,
1989
.