Fast Fourier Transform periodic interpolation method for superposition sums in a periodic unit cell

Fangzhou Ai, Vitaliy Lomakin
{"title":"Fast Fourier Transform periodic interpolation method for superposition sums in a periodic unit cell","authors":"Fangzhou Ai, Vitaliy Lomakin","doi":"arxiv-2312.02376","DOIUrl":null,"url":null,"abstract":"We propose a Fast Fourier Transform based Periodic Interpolation Method\n(FFT-PIM), a flexible and computationally efficient approach for computing the\nscalar potential given by a superposition sum in a unit cell of an infinitely\nperiodic array. Under the same umbrella, FFT-PIM allows computing the potential\nfor 1D, 2D, and 3D periodicities for dynamic and static problems, including\nproblems with and without a periodic phase shift. The computational complexity\nof the FFT-PIM is of $O(N \\log N)$ for $N$ spatially coinciding sources and\nobserver points. The FFT-PIM uses rapidly converging series representations of\nthe Green's function serving as a kernel in the superposition sum. Based on\nthese representations, the FFT-PIM splits the potential into its near-zone\ncomponent, which includes a small number of images surrounding the unit cell of\ninterest, and far-zone component, which includes the rest of an infinite number\nof images. The far-zone component is evaluated by projecting the non-uniform\nsources onto a sparse uniform grid, performing superposition sums on this\nsparse grid, and interpolating the potential from the uniform grid to the\nnon-uniform observation points. The near-zone component is evaluated using an\nFFT-based method, which is adapted to efficiently handle non-uniform\nsource-observer distributions within the periodic unit cell. The FFT-PIM can be\nused for a broad range of applications, such as periodic problems involving\nintegral equations in computational electromagnetic and acoustic, micromagnetic\nsolvers, and density functional theory solvers.","PeriodicalId":501061,"journal":{"name":"arXiv - CS - Numerical Analysis","volume":" 9","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2023-12-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - CS - Numerical Analysis","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2312.02376","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0

Abstract

We propose a Fast Fourier Transform based Periodic Interpolation Method (FFT-PIM), a flexible and computationally efficient approach for computing the scalar potential given by a superposition sum in a unit cell of an infinitely periodic array. Under the same umbrella, FFT-PIM allows computing the potential for 1D, 2D, and 3D periodicities for dynamic and static problems, including problems with and without a periodic phase shift. The computational complexity of the FFT-PIM is of $O(N \log N)$ for $N$ spatially coinciding sources and observer points. The FFT-PIM uses rapidly converging series representations of the Green's function serving as a kernel in the superposition sum. Based on these representations, the FFT-PIM splits the potential into its near-zone component, which includes a small number of images surrounding the unit cell of interest, and far-zone component, which includes the rest of an infinite number of images. The far-zone component is evaluated by projecting the non-uniform sources onto a sparse uniform grid, performing superposition sums on this sparse grid, and interpolating the potential from the uniform grid to the non-uniform observation points. The near-zone component is evaluated using an FFT-based method, which is adapted to efficiently handle non-uniform source-observer distributions within the periodic unit cell. The FFT-PIM can be used for a broad range of applications, such as periodic problems involving integral equations in computational electromagnetic and acoustic, micromagnetic solvers, and density functional theory solvers.
查看原文
分享 分享
微信好友 朋友圈 QQ好友 复制链接
本刊更多论文
周期单元格中叠加和的快速傅里叶变换周期插值方法
我们提出了一种基于快速傅立叶变换的周期插值方法(FFT-PIM),这是一种灵活且计算效率高的方法,用于计算由无限周期阵列的单元格中的叠加和给出的标量势。在相同的框架下,FFT-PIM允许计算动态和静态问题的1D、2D和3D周期性的潜力,包括有或没有周期性相移的问题。对于$N$空间重合的源和观测点,FFT-PIM的计算复杂度为$O(N \log N)$。FFT-PIM使用格林函数的快速收敛级数表示作为叠加和的核。基于这些表示,FFT-PIM将势分解为近区分量,其中包括感兴趣的单元格周围的少量图像,以及远区分量,其中包括无限数量图像的其余部分。通过将非均匀源投影到稀疏的均匀网格上,在稀疏网格上进行叠加和,并将均匀网格的势插值到非均匀观测点,来评估远区分量。使用基于fft的方法评估近区分量,该方法适用于有效处理周期单元胞内的非均匀源-观测器分布。FFT-PIM可用于广泛的应用,例如涉及计算电磁和声学积分方程的周期性问题,微磁解算器和密度泛函理论解算器。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 去求助
来源期刊
自引率
0.00%
发文量
0
期刊最新文献
Analysis of the SQP Method for Hyperbolic PDE-Constrained Optimization in Acoustic Full Waveform Inversion Detection of a piecewise linear crack with one incident wave Randomized quasi-Monte Carlo and Owen's boundary growth condition: A spectral analysis Energy stable gradient flow schemes for shape and topology optimization in Navier-Stokes flows Exponential time propagators for elastodynamics
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
已复制链接
已复制链接
快去分享给好友吧!
我知道了
×
扫码分享
扫码分享
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1