An accelerated Poisson solver based on multidomain spectral discretization

dc.citation.journalTitleBIT Numerical Mathematics
dc.contributor.authorBabb, Tracy
dc.contributor.authorGillman, Adrianna
dc.contributor.authorHao, Sijia
dc.contributor.authorMartinsson, Per-Gunnar
dc.date.accessioned2018-09-27T17:33:25Z
dc.date.available2018-09-27T17:33:25Z
dc.date.issued2018
dc.description.abstractThis paper presents a numerical method for variable coefficient elliptic PDEs with mostly smooth solutions on two dimensional domains. The method works best for domains that can readily be mapped onto a rectangle, or a collection of nonoverlapping rectangles. The PDE is discretized via a multi-domain spectral collocation method of high local order (order 30 and higher have been tested and work well). Local mesh refinement results in highly accurate solutions even in the presence of local irregular behavior due to corner singularities, localized loads, etc. The system of linear equations attained upon discretization is solved using a direct (as opposed to iterative) solver with O(N1.5)O(N1.5) complexity for the factorization stage and O(NlogN)O(Nlogā”N) complexity for the solve. The scheme is ideally suited for executing the elliptic solve required when parabolic problems are discretized via time-implicit techniques. In situations where the geometry remains unchanged between time-steps, very fast execution speeds are obtained since the solution operator for each implicit solve can be pre-computed.
dc.identifier.citationBabb, Tracy, Gillman, Adrianna, Hao, Sijia, et al.. "An accelerated Poisson solver based on multidomain spectral discretization." <i>BIT Numerical Mathematics,</i> (2018) Springer: https://doi.org/10.1007/s10543-018-0714-0.
dc.identifier.digitalBabb2018
dc.identifier.doihttps://doi.org/10.1007/s10543-018-0714-0
dc.identifier.urihttps://hdl.handle.net/1911/102722
dc.language.isoeng
dc.publisherSpringer
dc.rightsThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
dc.rights.urihttp://creativecommons.org/licenses/by/4.0/
dc.titleAn accelerated Poisson solver based on multidomain spectral discretization
dc.typeJournal article
dc.type.dcmiText
dc.type.publicationpublisher version
Files
Original bundle
Now showing 1 - 1 of 1
Loading...
Thumbnail Image
Name:
Babb2018.pdf
Size:
1.44 MB
Format:
Adobe Portable Document Format