Adaptive Fast Multipole Regularized Meshless Method for Large-Scale Three Dimensional Potential Problems
-
摘要: 正则化无网格法(regularized meshless method, RMM)是一种新的边界型无网格数值离散方法.该方法克服了近年来引起广泛关注的基本解方法(method of fundamental solutions, MFS)的虚假边界缺陷,继承了其无网格、无数值积分、易实施等优点.另一方面,RMM方法同MFS方法的插值方程都涉及非对称稠密系数矩阵,运用常规代数方程的迭代法求解时都要求O(N2)量级的乘法计算量和存储量.随着问题自由度的增加,该方法的计算量增加极快,效率较低,一般难以计算大规模问题.为了克服这个缺点,利用对角形式的快速多级算法(fast multipole method, FMM)来加速RMM方法,发展了快速多级正则化无网格法(fast multipole regularized mesheless method, FM-RMM).该方法无需数值积分并且具有O(N)量级的计算量和存储量,可有效地求解大规模工程问题.数值算例表明,FM-RMM算法可成功在内存为4GB的Core(TM)Ⅱ台式机上求解高达百万级自由度的三维位势问题.Abstract: The regularized meshless method (RMM) is a new meshless boundary collocation method. This method overcame the perplexing fictitious boundary associated in the method of fundamental solutions (MFS), while inherited all its merits being truly meshless, integration-free, and easy toprogram. Like the MFS, the RMM also produced dense and nonsymmetric coefficient interpolation matrix, which required O(N2) multiplication operations and memory requirements in an iterative solution procedure. Since the calculation operations would dramatically increase with the number of DOF, the RMM was computationally too expensive to solve largescale problems. In order to overcome this bottleneck, this study combined the RMM with the popular diagonal form fast multipole method (FMM) to develop the fast multipole regularized meshless method (FM-RMM). The proposed scheme was integrationfree and meshfree and significantly reduced
-
[1] Young D L, Chen K H, Lee C W.Novel meshless method for solving the potential problems with arbitrary domain[J]. Journal of Computational Physics, 2005, 209(1): 290-321. [2] 赵光明, 宋顺成.无网格Galerkin法与有限元耦合新算法[J].应用数学和力学, 2005, 26(8): 899-904.(ZHAO Guang-ming, SONG Shun-cheng.New algorithm of coupling element-free Galerkin with finite element method[J]. Applied Mathematics and Mechanics(English Edition) , 2005, 26(8): 982988.) [3] 苏成, 韩大建.正交各向异性弹性力学平面问题的样条虚边界元法[J].应用数学和力学, 2002, 23(4): 400-406.(SU Cheng, HAN Da-jian.Elastic analysis of orthotropic plane problems by the spline fictitious boundary element method[J].Applied Mathematics and Mechanics(English Edition) , 2002, 23(4): 446-453.) [4] Fairweather G, Karageorghis A.The method of fundamental solutions for elliptic boundary value problems[J]. Advances in Computational Mathematics, 1998, 9(1): 69-95. [5] Cheng H, Crutchfield W Y, Gimbutas Z, Greengard L, Ethridge J F, Huang J F.A wideband fast multipole method for the Helmholtz equation in three dimensions[J]. Journal of Computational Physics, 2006, 216(1): 300-325. [6] Jin B T, Zheng Y.A meshless method for some inverse problems associated with the Helmholtz equation[J]. Computer Methods in Applied Mechanics and Engineering, 2006, 195(19/22): 2270-2288. [7] Greengard L, Rokhlin V.A fast algorithm for particle simulations[J]. Journal of Computational Physics, 1987, 73(2): 325-348. [8] Carrier J, Greengard L, Rokhlin V.A fast adaptive multipole algorithm for particle simulations[J]. SIAM Journal on Scientific and Statistical Computing, 1988, 9(4): 669-686. [9] Greengard L, Rokhlin V.A new version of the fast multipole method for the Laplace equation in three dimensions[J]. Acta Numerica, 1997, 6: 229-269. [10] Amini S, Harris P J.A comparison between various boundary integral formulations of the exterior acoustic problem[J]. Computer Methods in Applied Mechanics and Engineering, 1990, 84(1): 59-75. [11] 周焕林, 牛忠荣, 王秀喜.位势问题边界元法中几乎奇异积分的正则化[J].应用数学和力学, 2003, 24(10): 10691074.(ZHOU Huan-lin, NIU Zhong-rong, WANG Xiu-xi.Regularization of nearly singular integrals in the boundary element method of potential problems[J]. Applied Mathematics and Mechanics(English Edition) , 2003, 24(10): 1208-1214.) [12] Wang H T, Yao Z H, Wang P B.On the preconditioners for fast multipole boundary element methods for 2D multi-domain elastostatics[J]. Engineering Analysis With Boundary Elements, 2005, 29(7): 673-688. [13] Shen L, Liu Y J.An adaptive fast multipole boundary element method for three-dimensional potential problems[J]. Computational Mechanics, 2007,39(6): 681-691. [14] Bapat M S, Shen L, Liu Y J.Adaptive fast multipole boundary element method for three-dimensional half-space acoustic wave problems[J]. Engineering Analysis With Boundary Elements, 2009, 33(8/9): 1113-1123. [15] Tausch J.The variable order fast multipole method for boundary integral equations of the second kind[J]. Computing, 2004, 72(3): 267-291. [16] ZHAO Li-bin,YAO Zhen-hao.Fast multipole BEM for 3-D elastostatic problems with applications for thin structures[J]. Tsinghua Science & Technology, 2005, 10(1): 67-75. [17] Nishimura N, Yoshida K, Kobayashi S.A fast multipole boundary integral equation method for crack problems in 3D[J]. Engineering Analysis With Boundary Elements, 1999, 23(1): 97-105. [18] Yoshida K, Nishimura N, Kobayashi S.Application of new fast multipole boundary integral equation method to crack problems in 3D[J]. Engineering Analysis With Boundary Elements, 2001, 25(4/5): 239-247. [19] Nishimura N.Fast multipole accelerated boundary integral equation methods[J]. Applied Mechanics Reviews, 2002, 55(4): 299-324. [20] Gu Y, Chen W, Zhang J.Investigation on near-boundary solutions by singular boundary method[J]. Engineering Analysis With Boundary Elements, 2012, 36(8): 1173-1182. [21] Barra L P S, Coutinho A L G A, Mansur W J, Telles J C F.Iterative solution of BEM equations by GMRES algorithm[J]. Computers & Structures, 1992, 44(6): 1249-1253. [22] Greengard L, Huang J F.A new version of the fast multipole method for screened coulomb interactions in three dimensions[J]. Journal of Computational Physics, 2002, 180(2): 642-658.
点击查看大图
计量
- 文章访问数: 1941
- HTML全文浏览量: 78
- PDF下载量: 1491
- 被引次数: 0