Phys. Rev. E 63, 066701 (2001) (8 pages)
Full text: PDF (852 kB)
A highly efficient algorithm for the reconstruction of microstructures of heterogeneous media from spatial correlation functions is presented. Since many experimental techniques yield two-point correlation functions, the restoration of heterogeneous structures, such as composites, porous materials, microemulsions, ceramics, or polymer blends, is an inverse problem of fundamental importance. Similar to previously proposed algorithms, the new method relies on Monte Carlo optimization, representing the microstructure on a discrete grid. An efficient way to update the correlation functions after local changes to the structure is introduced. In addition, the rate of convergence is substantially enhanced by selective Monte Carlo moves at interfaces. Speedups over prior methods of more than two orders of magnitude are thus achieved. Moreover, an improved minimization protocol leads to additional gains. The algorithm is ideally suited for implementation on parallel computers. The increase in efficiency brings new classes of problems within the realm of the tractable, notably those involving several different structural length scales and/or components.
©2001 The American Physical Society
02.70.-c | Computational techniques; | |
61.43.Bn | Structural modeling | |
05.10.Ln | Computational methods in statistical physics and nonlinear dynamics | |
81.07.-b | Nanoscale materials and structures |
Keywords: heterogeneous medium; two-point correlation functions; reconstruction; stereology; mathematical morphology;