北京大学学报(自然科学版) ›› 2026, Vol. 62 ›› Issue (3): 459-473.DOI: 10.13209/j.0479-8023.2026.005

上一篇    下一篇

适用于周期性问题的Poisson模块库求解器的开发及验证

王雯1, 赵佳敏1, 李想1, 李青2,3,†   

  1. 1. 旱区水工程生态环境全国重点实验室, 西安理工大学, 西安 710048 2. 天目山实验室, 杭州 311115 3. 北京航空航天大学航空科学与工程学院, 北京 100083
  • 收稿日期:2025-03-28 修回日期:2025-08-21 出版日期:2026-05-20 发布日期:2026-05-20
  • 基金资助:
    国家重点研发计划(2022YFF1300803)和天目山实验室自主科研项目(TK-2024-D-004)资助

Development and Validation of a Poisson Solver Library for Periodic Boundary Problems

WANG Wen1, ZHAO Jiamin1, LI Xiang1, LI Qing2,3,†   

  1. 1. State Key Laboratory of Water Engineering Ecology and Environment in Arid Area, Xi’an University of Technology, Xi’an 710048 2. Tianmushan Laboratory, Hangzhou 311115 3. School of Aeronautical Science and Engineering, Beihang University, Beijing 100083
  • Received:2025-03-28 Revised:2025-08-21 Online:2026-05-20 Published:2026-05-20

摘要:

对于周期边界条件Poisson问题, 采用传统的SOR求解算法使Poisson求解器难以快速收敛。为解决这一问题, 针对SOR算法的缺点, 利用多网格加速技术, 开发一种新的有限差分法多网格迭代泊松求解器Poisson-SOR-MG。此外, 提出一种基于1D快速傅里叶变换(FFT)和3D并行通信算法耦合的3D傅里叶谱变换算法, 摒弃繁复的转置技术。在此基础上, 开发和验证基于标准2π周期的Poisson-Fourier求解器, 将Poisson求解器的精度提高至谱精度。然而, 基于标准2π周期的Poisson-Fourier求解器无法应用到广泛的工程问题, 因此, 基于坐标拉伸变换的理论分析, 进一步开发和验证基于一般周期的Poisson-Fourier求解器。

关键词: 多重网格技术, 3D傅里叶谱变换, Poisson求解器, SOR迭代法, 坐标拉伸变换

Abstract:

For Poisson problems with periodic boundary conditions, the use of the traditional SOR solution algorithm does not permit the rapid convergence of the Poisson solver. To address this limitation, a novel multigrid iterative Poisson solver (Poisson-SOR-MG) based on finite difference methods was developed, leveraging multigrid acceleration techniques to overcome the drawbacks of SOR. Additionally, a 3D Fourier spectral transform algorithm was designed by coupling 1D Fast Fourier Transform (FFT) with 3D parallel communication algorithms, eliminating cumbersome transpose techniques. Building on this, a standard 2π-periodic Poisson-Fourier solver was developed and validated, achieving spectral accuracy. However, the standard 2π-periodic Poisson-Fourier solver cannot be directly applied to general engineering problems. To resolve this, theoretical analysis based on coordinate stretching transformations was conducted, leading to the further development and validation of a generalized arbitrary-periodic Poisson-Fourier solver.

Key words: