A fast direct algorithm for the solution of the laplace equation on regions with fractal boundaries

22Citations
Citations of this article
14Readers
Mendeley users who have this article in their library.

This article is free to access.

Abstract

An algorithm is presented for the rapid direct solution of the Laplace equation on regions with fractal boundaries. In a typical application, the numerical simulation has to be on a very large scale involving at least tens of thousands of equations with as many unknowns, in order to obtain any meaningful results. Attempts to use conventional techniques have encountered insurmountable difficulties, due to excessive CPU time requirements of the computations involved. Indeed, conventional direct algorithms for the solution of linear systems require order O(N3) operations for the solution of an N × N-problem, while classical iterative methods require order O(N2) operations, with the constant strongly dependent on the problem in question. In either case, the computational expense is prohibitive for large-scale problems. The direct algorithm of the present paper requires O(N) operations with a constant dependent only on the geometry of the boundary, making it considerably more practical for large-scale problems encountered in the computation of harmonic measure of fractals, complex iteration theory, potential theory, and growth phenomena such as crystallization, electrodeposition, viscous fingering, and diffusion-limited aggregation. © 1994 Academic Press, Inc.

Cite

CITATION STYLE

APA

Jones, P., Ma, J., & Rokhlin, V. (1994). A fast direct algorithm for the solution of the laplace equation on regions with fractal boundaries. Journal of Computational Physics, 113(1), 35–51. https://doi.org/10.1006/jcph.1994.1116

Register to see more suggestions

Mendeley helps you to discover research relevant for your work.

Already have an account?

Save time finding and organizing research with Mendeley

Sign up for free