采用快速多极方法的边界元大规模计算.ppt

内 容 提 要 引言 快速多极边界元法的基本原理 快速多极边界元法的精度和效率验证 二维和三维复合材料模拟中的应用 快速多极边界元大规模并行计算 弹性裂纹体和裂纹扩展模拟 国际上比较成功的其它方面应用 进一步研究方向 边界元法简介 快速多极算法发展概况 快速多极边界元法研究概况 本研究组相关研究概况 边界元法简介 (Boundary Element Method, BEM) 求解连续介质问题的一种数值方法 只需在边界离散单元，降低一维 易于处理无限和半无限域问题 应用于力学、热学、电磁、声学等领域 传统边界元法的缺点：不适合处理大规模问题 方程组 A: 非对称满阵 直接求解方法 (Gauss Elimination) 存储量: O(N 2) 计算量: O(N 3) N: 未知量的个数 (Number of DOFs) 边界元迭代算法研究 (GMRES, CG, CGS, Bi-CGSTAB) Mullen等 (1987) ：将几种迭代法用于Laplace方程的BEM迭代求解 Barra等 (1992) ：将GMRES用于二维弹性力学的BEM迭代求解 Mansur等(1992) ：将用于FEM的迭代算法的非对称版本用于BEM Prasad等 (1994) ：验证了GMRES用于各种BEM格式的迭代求解的效率 迭代算法求解效率 存储量: O (N 2) 计算量: O (kN 2) k: 迭代步数 常规边界元法对解题规模的制约 例如： N = 1,000 双精度 存储量: 8 MB N = 10,000 双精度 存储量: 800 MB 又例如： N = 40,000 双精度 8台PIII-600微机并行求解，需要 43 h (2000) 求解方法已经成为制约边界元法发展和应用的瓶颈！ 如何提高边界元的解题规模？ 新硬件：代价昂贵，无竞争力 新方法： 快速多极算法 (Fast multipole method, 简称FMM) 存储量: O(N) 计算量: O(N) 快速多极算法发展概况 传统的求解 N particle 问题算法的计算量： O(N 2) 静电势 静电力 Barnes, Hut (1986) ： 树结构算法 (tree code) 计算量：O (NlogN) J. Barnes, P. Hut. A hierarchical O (NlogN) force calculation algorithm. Nature, 324, (1986), 446-449. 快速多极算法发展概况 Rokhlin, Greengard (1987)： 快速多极算法 (Fast multipole method, FMM)： O(N) 采用了： 多极展开 (multipole expansion) 局部展开 (local expansion) Carrier等 (1988) 自适应快速多极算法 (Adaptive FMM)： O(N) 快速多极算法发展概况 Greengard, Rokhlin (1997): 3D新型快速多极算法 (New Version FMM) O(N) 引入 指数展开 (Exponential expansion) Hrycak, Rokhlin (1998)： 2D新型快速多极算法 O(N) 快速多极边界元法研究概况 Yamada, Hayami (1995)： 二维弹性力学 复数域Taylor展开 O (NlogN) Pierce, Napier (1995)： 二维弹性力学 实数域Taylor展开 O (NlogN) Hayami, Sauter (1998)： 三维弹性力学 实数域Taylor展开 O (NlogN) Nishimura等(1999)： 三维含裂纹的Laplace方程 球谐函数展开 O(N) Popov, Power (2001)： 三维弹性力学 实数域Taylor展开 O(