The growth of streamer trees in insulating fluids (a submicrosecond process that triggers high-voltage breakdown) has been simulated with a combination of parallel-coding tools. Large grids and arrays display well the multifractal, self-avoiding character of the streamer trees. Three physical cases have been approximated by different power-law weightings of the statistical growth filter: dense anode trees, in the uniform field; sparse cathode trees (a rarer experimental case); and ultrasparse anode trees (seen in some fluids of higher viscosity). The model is contained in a software package that is written in Fortran 90 with data parallel extensions for distributed execution. These extensions encapsulate an underlying, invisible message-passing environment, thus enabling the solution of memory-intensive problems on a group of limited-memory processors. Block partitioning creates processes of reasonable size, which operate in parallel like small copies of the original code. The user needs only to express his model in transparent array-directed commands; parallel interfacing between blocks is handled invisibly. Breakdown is performed in parallel in each of the local blocks. Results are presented for experiments run on eight and nine nodes of the IBM SP2, and four and eight nodes of the SGI Onyx and Origin, three examples of multiple-processor machines. Display is carried out in three dimensions. Timing of the growth can be shown by color banding or by frame animation of the results. The adequacy of the growth rules and size scaling are tested by comparing the simulations against snapshots from high-voltage discharge events.

1.
L.
Pietronero
and
H. J.
Wiesmann
,
Z. Phys. B
70
,
87
(
1988
).
2.
W. G.
Chadband
,
IEEE Trans. Electr. Insul.
23
,
697
(
1988
).
3.
W. G.
Chadband
and
T. M.
Sufian
,
IEEE Trans. Electr. Insul.
20
,
239
(
1985
).
4.
S. Satpathy, Fractals in Physics (Elsevier Science, Amsterdam, 1986), p. 173.
5.
T. J. Lewis, in The Liquid State and its Electrical Properties, edited by E. E. Kunhardt, L. G. Christophorou, and L. E. Luessen, Proceedings the NATO Advanced Study Institute, 5–17 July 1989, Sintra, Portugal (Plenum, New York, 1988), pp. 431–453.
6.
Larry Davis and Mitchell Murphy of the University of Maryland Institute for Advanced Computing Systems provided access to the CM-2.
7.
W. Gropp, E. Lusk, and A. Skjellum, Using MPI (MIT Press, Cambridge, 1994), especially pp. 60–80.
8.
J. E. Devaney and J. Hagedorn, National Institute of Standards and Technology internal paper, May 7, 1996.
9.
J. Hagedorn and A. Heckert, DPARLIB User’s Guide, draft of NIST software library documentation, March 18, 1997.
10.
L.
Niemyer
,
L.
Pietronero
, and
H. J.
Wiesmann
,
Phys. Rev. Lett.
52
,
1033
(
1984
).
11.
L.
Pietronero
and
H. J.
Wiesmann
,
J. Stat. Phys.
36
,
909
(
1984
).
12.
H. J. Wiesmann and L. Pietronero, in Ref. 4, p. 151.
13.
L. Pietronero, C. Evertsz, and H. J. Wiesmann, in Ref. 4, p. 159.
14.
M. Murat, in Ref. 4, p. 169.
15.
E.
Arian
,
P.
Alstrom
,
A.
Aharony
, and
H. E.
Stanley
,
Phys. Rev. Lett.
63
,
2005
(
1989
).
16.
P. D.
Beale
and
P. M.
Duxbury
,
Phys. Rev. B
37
,
2785
(
1988
).
17.
J. J. Modi, Parallel Algorithms and Matrix Computation (Clarendon, Oxford, 1988), pp. 135–143.
18.
A.
Sanchez
,
F.
Guinea
,
L. M.
Sander
,
V.
Hakim
, and
E.
Louis
,
Phys. Rev. E
48
,
1296
(
1993
). In these treatments, idealizations such as radial growth to a remote spherical-shell boundary, with normalized probabilities, favor growth having a simple, well-defined fractal dimension.
19.
T. A.
Witten
and
L. M.
Sander
,
Phys. Rev. Lett.
47
,
1400
(
1981
);
T. A.
Witten
and
L. M.
Sander
,
Phys. Rev. B
27
,
5686
(
1983
).
20.
A.
Erzan
,
L.
Pietronero
, and
A.
Vespignani
,
Rev. Mod. Phys.
67
,
545
(
1995
).
21.
M. M.
Clark
,
L. M.
Ruff
, and
H. L.
Scott
,
Comput. Phys.
10
,
584
(
1996
).
22.
K. Binder and D. W. Heermann, Monte Carlo Simulation in Statistical Physics (Springer, Berlin, 1988), especially, pp. 28–29.
23.
D. Burns and R. B. Daoud, Supercomputing Symposium ’94, June 1994, Toronto, Canada; code available at http//www.osc.edu/lam.html.
24.
For a similar treatment of a second-order operator, see, for example,
A.
Villareal
and
J. A.
Scales
,
Comput. Phys.
11
,
388
(
1997
), especially p. 392.
25.
I.
Beichl
and
F. E.
Sullivan
,
Comput. Sci. Eng.
4
,
91
(
1997
);
see also
A. B.
Bortz
,
M. H.
Kalos
, and
J. L.
Lebowitz
,
J. Comput. Phys.
17
,
10
(
1975
) especially Appendix B;
J. G.
Amar
and
F.
Family
,
Phys. Rev. Lett.
74
,
2066
(
1995
).
26.
K. L. Stricklett, E. F. Kelley, H. Yamashita, H. Kawai, and C. Fenimore Conference Record, IEEE International Symposium on Electrical Insulation, 1990, pp. 161–164. These data were taken by means of a high-speed “optical storage line” technique which limited the available spatial resolution. Dr. Stricklett made his original data available to the authors.
27.
P. K. Watson, Conference Record of the 10th International Conference on Conduction and Breakdown in Dielectric Liquids, Grenoble, France, 10–14 September 1990, pp. 22–28.
28.
C. Fenimore, Conference Record, IEEE International Symposium on Electrical Insulation, Cambridge, MA, 5–8 June 1988, pp. 27–30.
29.
K. Miyano, K. Yamazawa, M. Pompili, C. Mazzetti, and H. Yamashita, 1997 Annual Report of the Conference on Electrical Insulation and Dielectric Phenomena, Minneapolis, 19–22 October 1997, Vol. II, pp. 640–643 (available from IEEE Service Center, Piscataway, NJ). We refer especially to Figs. 3 and 4 in this paper. The authors associate the properties of the observed streamers with the viscosity of the perfluoro liquids used.
30.
J. C. Adams, W. S. Brainerd, J. T. Martin, B. T. Smith, and J. L. Wagener, Fortran 90 Handbook (McGraw–Hill, New York, 1992).
This content is only available via PDF.
You do not currently have access to this content.