多重网格方法求解两类Helmholtz方程.doc

多重网格方法求解两类Helmholtz 方程 赵 婧，聂玉峰 ZHAO Jing, NIE Yufeng 西北工业大学 应用数学系，西安 710129 Department of Applied Mathematics, Northwestern Polytechnical University, Xi’an 710129, China ZHAO Jing, NIE Yufeng. Multigrid method applied to two kinds of Helmholtz equations. Computer Engi- neering and Applications, 2012, 48（22）：41-44. Abstract：The procedure of multigrid method is described in detail and its properties are discussed by applying mul- tigrid to both positive definite Helmholtz equation and indefinite Helmholtz equation. Comparisons of convergence performance are implemented between V cycle, W cycle and F cycle. It is shown that, by solving positive definite Helmholtz equation, multigrid method has high computing efficiency. However, as the value of wave-number in- creases in indefinite Helmholtz equation, the results obtained by multigrid method diverge. The reasons lie in both fine grid smoothing and coarse grid correction, and the improvement of Multigrid is still needed in future. Key words：multigrid; F cycle; positive definite Helmholtz; indefinite Helmholtz 摘 要：详细给出了多重网格方法的实现过程，借助正定 Helmholtz 方程及不定 Helmholtz 方程的求解来探讨多 重网格方法的特性。对多重网格 V 环、W 环以及 F 环三种不同迭代格式的收敛效果进行了对比。通过正定 Helmholtz 方程的求解，发现多重网格的确有很高的计算效率。对于不定 Helmholtz 方程，随着波数的增加，利 用多重网格方法得到结果不收敛，原因出在细网格光滑和粗网格矫正过程。如何针对此问题对多重网格进行 有效改进还有待进一步研究。 关键词：多重网格；F 环；正定Holmholtz 方程；不定Helmholtz 方程 文章编号：1002-8331（2012）22-0041-04 文献标识码：A 中图分类号：TP301.6 1 引言 多重网格方法的基本思想来源于套迭代（nested iteration）、松弛方法的误差光滑以及总下降（total reduction）。 自 从 上 世 纪 60 年 代 Fedorenko[1-2] 和 Bakhvalov[3] 第一次提出严格意义上的多重网格方法 以来，众多学者对此方法进行了多方面的研究。Brandt 于上世纪 70 年代从理论上证明了多重网格方法的可 靠性，并且介绍了非线性多重网格方法（FAS），自适 应技术（MLAT）[4]以及用于误差分析的局部傅里叶分 析方法[5]。此后，如何将多重网格方法推广到更一般 的计算区域，如何将其与有限单元法、有限差分法以 及有限容积法相结合以及它的自适应性、并行性等 研究课题也取得了众多发展[5]。 本文首先详细给出了多重网格的实现过程，然 后利用多重网格方法分别对正定 Helmholtz 方程及不 定 Helmholtz 方程进行求解，其中对比讨论了 V 环，W 环以及 F 环三种不同迭代格式的收敛效果。最后，对 于多重网格方法在求解不定 Helmholtz 方程时失效的  基金项目：国家自然科学基金项目（NoNo。 作者简介：赵婧（1986—），女，硕士研究生，主要研究领域为多重网格方法及其工程应用；聂玉峰（1968—），男，通讯作者，