On the Numerical Evaluation of Electrostatic Fields in Dense Random Dispersions of Cylinders

We consider the interface boundary value problem which arises in the evaluation of electrostatic fields in composite materials consisting of dense random dispersions of cylinders in a uniform background. This is a well-studied problem from the viewpoint of homogenization and effective medium theory, but one for which accurate numerical simulations have been difficult to obtain. Size effects, in particular, have been neglected due to the expense of solving the field equation in the presence of large numbers of close-to-touching inclusions. Such features require very fine discretizations, even with the use of adaptive gridding, and cause the linear systems which arise to be highly ill-conditioned. In this paper, we present a new integral equation method for the solution of the interface problem which uses a recently developed method of images to resolve the close-to-touching interactions and the fast multipole method to compute far field interactions. Only minutes of workstation time are needed to solve the field equation with thousands of inclusions, allowing us to carry out large-scale statistical studies of the effective conductivity of random two-phase materials at a variety of volume fractions and contrast ratios.

[1]  Graeme W. Milton,et al.  Bounds on the transport and optical properties of a two‐component composite material , 1981 .

[2]  Ross C. McPhedran,et al.  Bounds and exact theories for the transport properties of inhomogeneous media , 1981 .

[3]  G. Batchelor,et al.  Transport Properties of Two-Phase Materials with Random Structure , 1974 .

[4]  George Papanicolaou,et al.  Bounds for effective parameters of multicomponent media by analytic continuation , 1985 .

[5]  J. Willis,et al.  Variational and Related Methods for the Overall Properties of Composites , 1981 .

[6]  L. Greengard.  The Rapid Evaluation of Potential Fields in Particle Systems , 1988 .

[7]  G. Milton.  The coherent potential approximation is a realizable effective medium scheme , 1985 .

[8]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[9]  S. Torquato,et al.  Determination of the effective conductivity of heterogeneous media by Brownian motion simulation , 1990 .

[10]  Leslie Greengard,et al.  A renormalization method for the evaluation of lattice sums , 1994 .

[11]  Leslie Greengard,et al.  A numerical study of the ζ2 parameter for random suspensions of disks , 1995 .

[12]  M. A. Jaswon,et al.  Integral equation methods in potential theory and elastostatics , 1977 .

[13]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[14]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[15]  D. A. Dunnett.  Classical Electrodynamics , 2020, Nature.

[16]  Joseph B. Keller,et al.  Conductivity of a Medium Containing a Dense Array of Perfectly Conducting Spheres or Cylinders or Nonconducting Cylinders , 1963 .

[17]  S. Shtrikman,et al.  A Variational Approach to the Theory of the Effective Magnetic Permeability of Multiphase Materials , 1962 .

[18]  D. Kinderlehrer,et al.  Homogenization and effective moduli of materials and media , 1986 .

[19]  Salvatore Torquato,et al.  Random Heterogeneous Media: Microstructure and Improved Bounds on Effective Properties , 1991 .

[20]  A. Sangani,et al.  Transport Processes in Random Arrays of Cylinders. I. Thermal Conduction , 1988 .

[21]  O. D. Kellogg.  Foundations of potential theory , 1934 .

[22]  E. M. Lifshitz,et al.  Electrodynamics of continuous media , 1961 .

[23]  Graeme W. Milton,et al.  Bounds on the Electromagnetic, Elastic, and Other Properties of Two-Component Composites , 1981 .

[24]  Johan Helsing,et al.  Thin Bridges in Isotropic Electrostatics , 1996 .