{"id":{"repo_id":"uiuc","oai_identifier":"oai:www.ideals.illinois.edu:2142/106232"},"canonical_url":"https://search.dev.ndltd.org/etd/uiuc/oai:www.ideals.illinois.edu:2142/106232","repository":{"repo_id":"uiuc","name":"University of Illinois - Urbana-Champaign","base_url":"https://www.ideals.illinois.edu/oai-pmh"},"display":{"title":"High-order numerical methods for layer potential evaluation","abstract":"This thesis addresses a number of obstacles in the practical realization of integral equation methods for boundary value problems of elliptic partial differential equations. Layer potentials play a central role in the integral equation formulation of boundary value problems. However, the numerical evaluation of layer potentials presents significant practical challenges associated with issues of singular/near-singular quadrature and efficient scaling. Quadrature by Expansion (QBX) is a quadrature method that addresses many of the issues encountered in the evaluation of layer potentials, by providing a high-order and kernel/dimension independent quadrature scheme that is also acceleration-compatible, i.e., able to be integrated with a fast summation scheme to reduce the cost of $O(N^2)$ pairwise interactions to $O(N)$. The focus of the first part of this thesis is an acceleration scheme for QBX based on the Fast Multipole Method (FMM) that provides both high performance and strong accuracy guarantees for mathematical control over the error introduced by acceleration. The scheme extends to both two and three dimensions and Laplace and Helmholtz kernels. The contribution in this thesis also includes a geometry processing framework to prepare arbitrary smooth geometries for accurate quadrature with QBX, a comprehensive complexity and error analysis of the algorithm, and a cost model and study of key optimizations. This thesis also considers the application of layer potentials to the problem of numerical conformal mapping. It is shown how there is a strong connection between the double-layer potential and the Faber/Faber-Laurent polynomials, which are polynomials closely related to the series expansion of the Riemann map. From this, a scheme for computing the Riemann map of a piecewise smooth Jordan domain is devised. This latter work also provides an insight into the convergence behavior of QBX. Finally, this thesis describes a software abstraction for the design and implementation of time integration algorithms for the solution of initial value problems. This abstraction, Dagrt, has two aims. The first is a epresentation that decouples the mathematical specification of a time integration algorithm from its realization in a particular programming model, and the second is to give the user control over implementation details. We demonstrate the capabilities of this abstraction by presenting Leap, a collection of pre-written time integration algorithm specifications that includes complex multistep and multistage methods.","abstract_html":"This thesis addresses a number of obstacles in the practical realization of integral equation methods for boundary value problems of elliptic partial differential equations. Layer potentials play a central role in the integral equation formulation of boundary value problems. However, the numerical evaluation of layer potentials presents significant practical challenges associated with issues of singular/near-singular quadrature and efficient scaling. Quadrature by Expansion (QBX) is a quadrature method that addresses many of the issues encountered in the evaluation of layer potentials, by providing a high-order and kernel/dimension independent quadrature scheme that is also acceleration-compatible, i.e., able to be integrated with a fast summation scheme to reduce the cost of <span class=\"etd-inline-math\">O(N<sup>2</sup>)</span> pairwise interactions to $O(N)$. The focus of the first part of this thesis is an acceleration scheme for QBX based on the Fast Multipole Method (FMM) that provides both high performance and strong accuracy guarantees for mathematical control over the error introduced by acceleration. The scheme extends to both two and three dimensions and Laplace and Helmholtz kernels. The contribution in this thesis also includes a geometry processing framework to prepare arbitrary smooth geometries for accurate quadrature with QBX, a comprehensive complexity and error analysis of the algorithm, and a cost model and study of key optimizations. This thesis also considers the application of layer potentials to the problem of numerical conformal mapping. It is shown how there is a strong connection between the double-layer potential and the Faber/Faber-Laurent polynomials, which are polynomials closely related to the series expansion of the Riemann map. From this, a scheme for computing the Riemann map of a piecewise smooth Jordan domain is devised. This latter work also provides an insight into the convergence behavior of QBX. Finally, this thesis describes a software abstraction for the design and implementation of time integration algorithms for the solution of initial value problems. This abstraction, Dagrt, has two aims. The first is a epresentation that decouples the mathematical specification of a time integration algorithm from its realization in a particular programming model, and the second is to give the user control over implementation details. We demonstrate the capabilities of this abstraction by presenting Leap, a collection of pre-written time integration algorithm specifications that includes complex multistep and multistage methods.","abstract_has_math":true,"creators":["Wala, Mateusz Michal"],"institution":"University of Illinois at Urbana-Champaign","degree_name":"Ph.D.","degree_level":"Dissertation","degree_discipline":"Computer Science","degree_department":null,"school":null,"contributors":["Klöckner, Andreas","Olson, Luke","Fischer, Paul","Greengard, Leslie"],"advisors":[],"committee_chairs":[],"committee_members":[],"year":2020,"date_issued":"2020-03-02T21:58:20Z","date_published":"2020-03-02T21:58:20Z","updated_at":"2026-07-22T22:24:45Z","subjects":["integral equations","layer potentials","fast multipole method","singular integrals","quadrature by expansion","conformal mapping","time integration"],"languages":["en"],"rights":["Copyright 2019 Mateusz Wala"],"rights_urls":[],"identifier_entries":[]},"links":{"outbound_url":"http://hdl.handle.net/2142/106232","outbound_label":"Handle","outbound_source":"dc:identifier"},"metadata_groups":[{"id":"people","label":"People","entries":[{"key":"dc:contributor","label":"Contributor","values":["Klöckner, Andreas","Olson, Luke","Fischer, Paul","Greengard, Leslie"]},{"key":"dc:creator","label":"Author","values":["Wala, Mateusz Michal"]}]},{"id":"academic_context","label":"Academic Context","entries":[{"key":"dc:date","label":"Dc Date","values":["2020-03-02T21:58:20Z","2019-12-02","2019-12"]},{"key":"dc:type","label":"Dc Type","values":["text"]},{"key":"thesis:degree_discipline","label":"Discipline","values":["Computer Science"]},{"key":"thesis:degree_level","label":"Degree Level","values":["Dissertation"]},{"key":"thesis:degree_name","label":"Degree Name","values":["Ph.D."]},{"key":"thesis:institution_name","label":"Thesis Institution Name","values":["University of Illinois at Urbana-Champaign"]}]},{"id":"subjects_keywords","label":"Subjects and Keywords","entries":[{"key":"dc:subject","label":"Dc Subject","values":["integral equations","layer potentials","fast multipole method","singular integrals","quadrature by expansion","conformal mapping","time integration"]}]},{"id":"language_rights","label":"Language and Rights","entries":[{"key":"dc:language","label":"Dc Language","values":["en"]},{"key":"dc:rights","label":"Dc Rights","values":["Copyright 2019 Mateusz Wala"]}]},{"id":"identifiers","label":"Identifiers","entries":[{"key":"dc:identifier","label":"Identifier","values":["http://hdl.handle.net/2142/106232"]}]},{"id":"additional","label":"Additional Metadata","entries":[{"key":"dc:description","label":"Description","values":["This thesis addresses a number of obstacles in the practical realization of integral equation methods for boundary value problems of elliptic partial differential equations. Layer potentials play a central role in the integral equation formulation of boundary value problems. However, the numerical evaluation of layer potentials presents significant practical challenges associated with issues of singular/near-singular quadrature and efficient scaling. Quadrature by Expansion (QBX) is a quadrature method that addresses many of the issues encountered in the evaluation of layer potentials, by providing a high-order and kernel/dimension independent quadrature scheme that is also acceleration-compatible, i.e., able to be integrated with a fast summation scheme to reduce the cost of $O(N^2)$ pairwise interactions to $O(N)$. The focus of the first part of this thesis is an acceleration scheme for QBX based on the Fast Multipole Method (FMM) that provides both high performance and strong accuracy guarantees for mathematical control over the error introduced by acceleration. The scheme extends to both two and three dimensions and Laplace and Helmholtz kernels. The contribution in this thesis also includes a geometry processing framework to prepare arbitrary smooth geometries for accurate quadrature with QBX, a comprehensive complexity and error analysis of the algorithm, and a cost model and study of key optimizations. This thesis also considers the application of layer potentials to the problem of numerical conformal mapping. It is shown how there is a strong connection between the double-layer potential and the Faber/Faber-Laurent polynomials, which are polynomials closely related to the series expansion of the Riemann map. From this, a scheme for computing the Riemann map of a piecewise smooth Jordan domain is devised. This latter work also provides an insight into the convergence behavior of QBX. Finally, this thesis describes a software abstraction for the design and implementation of time integration algorithms for the solution of initial value problems. This abstraction, Dagrt, has two aims. The first is a epresentation that decouples the mathematical specification of a time integration algorithm from its realization in a particular programming model, and the second is to give the user control over implementation details. We demonstrate the capabilities of this abstraction by presenting Leap, a collection of pre-written time integration algorithm specifications that includes complex multistep and multistage methods.","Submission original under an indefinite embargo labeled 'Open Access'. The submission was exported from vireo on 2020-02-28 without embargo terms","The student, Mateusz Wala, accepted the attached license on 2019-12-02 at 15:01.","The student, Mateusz Wala, submitted this Dissertation for approval on 2019-12-02 at 15:11.","This Dissertation was approved for publication on 2019-12-02 at 16:09.","DSpace SAF Submission Ingestion Package generated from Vireo submission #14649 on 2020-02-28 at 17:14:43","Made available in DSpace on 2020-03-02T21:58:20Z (GMT). No. of bitstreams: 2 WALA-DISSERTATION-2019.pdf: 3603663 bytes, checksum: a9073a4c9bcdcb117e5dcbdb0d09d89b (MD5) LICENSE.txt: 4209 bytes, checksum: 20f8e4977e649e9d60f350ef5e7fa397 (MD5) Previous issue date: 2019-12-02"]},{"key":"dc:format","label":"Dc Format","values":["application/pdf"]},{"key":"dc:title","label":"Title","values":["High-order numerical methods for layer potential evaluation"]}]}],"canonical_facts":{"dc:contributor":["Klöckner, Andreas","Olson, Luke","Fischer, Paul","Greengard, Leslie"],"dc:creator":["Wala, Mateusz Michal"],"dc:date":["2020-03-02T21:58:20Z","2019-12-02","2019-12"],"dc:description":["This thesis addresses a number of obstacles in the practical realization of integral equation methods for boundary value problems of elliptic partial differential equations. Layer potentials play a central role in the integral equation formulation of boundary value problems. However, the numerical evaluation of layer potentials presents significant practical challenges associated with issues of singular/near-singular quadrature and efficient scaling. Quadrature by Expansion (QBX) is a quadrature method that addresses many of the issues encountered in the evaluation of layer potentials, by providing a high-order and kernel/dimension independent quadrature scheme that is also acceleration-compatible, i.e., able to be integrated with a fast summation scheme to reduce the cost of $O(N^2)$ pairwise interactions to $O(N)$. The focus of the first part of this thesis is an acceleration scheme for QBX based on the Fast Multipole Method (FMM) that provides both high performance and strong accuracy guarantees for mathematical control over the error introduced by acceleration. The scheme extends to both two and three dimensions and Laplace and Helmholtz kernels. The contribution in this thesis also includes a geometry processing framework to prepare arbitrary smooth geometries for accurate quadrature with QBX, a comprehensive complexity and error analysis of the algorithm, and a cost model and study of key optimizations. This thesis also considers the application of layer potentials to the problem of numerical conformal mapping. It is shown how there is a strong connection between the double-layer potential and the Faber/Faber-Laurent polynomials, which are polynomials closely related to the series expansion of the Riemann map. From this, a scheme for computing the Riemann map of a piecewise smooth Jordan domain is devised. This latter work also provides an insight into the convergence behavior of QBX. Finally, this thesis describes a software abstraction for the design and implementation of time integration algorithms for the solution of initial value problems. This abstraction, Dagrt, has two aims. The first is a epresentation that decouples the mathematical specification of a time integration algorithm from its realization in a particular programming model, and the second is to give the user control over implementation details. We demonstrate the capabilities of this abstraction by presenting Leap, a collection of pre-written time integration algorithm specifications that includes complex multistep and multistage methods.","Submission original under an indefinite embargo labeled 'Open Access'. The submission was exported from vireo on 2020-02-28 without embargo terms","The student, Mateusz Wala, accepted the attached license on 2019-12-02 at 15:01.","The student, Mateusz Wala, submitted this Dissertation for approval on 2019-12-02 at 15:11.","This Dissertation was approved for publication on 2019-12-02 at 16:09.","DSpace SAF Submission Ingestion Package generated from Vireo submission #14649 on 2020-02-28 at 17:14:43","Made available in DSpace on 2020-03-02T21:58:20Z (GMT). No. of bitstreams: 2 WALA-DISSERTATION-2019.pdf: 3603663 bytes, checksum: a9073a4c9bcdcb117e5dcbdb0d09d89b (MD5) LICENSE.txt: 4209 bytes, checksum: 20f8e4977e649e9d60f350ef5e7fa397 (MD5) Previous issue date: 2019-12-02"],"dc:format":["application/pdf"],"dc:identifier":["http://hdl.handle.net/2142/106232"],"dc:language":["en"],"dc:rights":["Copyright 2019 Mateusz Wala"],"dc:subject":["integral equations","layer potentials","fast multipole method","singular integrals","quadrature by expansion","conformal mapping","time integration"],"dc:title":["High-order numerical methods for layer potential evaluation"],"dc:type":["text"],"thesis:degree_discipline":["Computer Science"],"thesis:degree_level":["Dissertation"],"thesis:degree_name":["Ph.D."],"thesis:institution_name":["University of Illinois at Urbana-Champaign"]},"updated_at":"2026-07-22T22:24:45Z"}