自然边界元和有限元耦合法的并行化研究.pdfVIP

自然边界元与有限元耦合法的并行化研究· 赵慧明董正筑刘玉庆 (中国矿业大学应用教学与力学系．徐州，221008) j_要本文简要地介绍了自然边界元与有限元耦合法并行化的实现方法，并对并行算法的 基本理论及并行计算的平台PvM和删作了简要地介绍． 关■词 自然边界元，有限元，耦合法．并行lL／ O 引 雷 自然边界元方法是一种有着独特优点的新型的边界元方法，它与有限元方法基于相同 的变分原理，这为二者的耦合提供了坚实的数学理论基础．自然边界元与有限元的藕合法为 解决裂纹问题、不规则区域问题和无限区域问题的求解提供了新的思路．自然边界元与有限 元的耦合法有多种形式：既可以把应用自然边界元方法的区域当做有限元的一个特殊单元， 联合列方程进行求解：又可以参照区域分解的思想，在两类子区域上分别应用自然边界元 方法和有限元方法，然后应用迭代法求解． 在实际工程中，研究对象(如岩石)中含有许多裂纹和微裂纹，边界也很不规范，如 果采用简单的耦合法，则计算效率很低．这时如果采用并行算法，大量的同类计算可蚍并行 完成，如在同一迭代步内，各个子区域内的计算采用并行，这样就可以大大提高计算效率． 1自然边界元法简介 自然边界元方法是根据Green函数和Green公式，将微分方程边值问题归化为边界上 强奇异积分方程(或超奇异积分方程，本方法中称之为自然积分方程)，然后再化成相应的变 分形式在边界上离散化求解的一类数值计算方法．与传统的边界元方法相比，除具有传统边 界元方法的所有优点以外，自然边界元方法还有许多独特的优点： 【1)推导直观自然．自然边界归化时，从Gr∞n函数和Gre％公式出发，可以轻易地得到 Poisson积分公式，然后以Neum∞微分边值算子作用于积分公式两边，再令参变量从区域 内部向边界取极限，即可得自然积分公式． (2)能量不在自然边界归化下，自然积分算子导出的双线性泛函的能量值与区域内微分 算子导出的原问题的双线性泛函有相同的值．即能量泛函值保持不变．因此，原边值问题的 许多有用性质，如双线性型的对称性、强制性等均被保持．由于这一特点，自然积分方程解 的存在唯一性及稳定性等结果也随之而得． (3)计算量小．用自然边界归化得到的刚度矩阵是对称正定的，尤其是边界为圆周时． 刚度矩阵还具有某种循环特性．因此，我们并不需要计算矩阵的每一个分量，而只要计算大 约半行系数就行了．这与传统的边界元方法相比，计算量确实大为减少． ’本课题得到科学与工程计算国家重点实验室支持 294 (4)方程形式唯一．自然积分方程与其归化途径及所选择的基本解无关，它仅由原边值 问怠唯一确定，它准确地反映了边值问题解与互补的微分边值之间的本质关系．我们可以采 用各种不同的途径，如Gr∞n函数法：硒llri日变换法、F0嘶口级数法、复变函数法等来推 导自然积分方程，但殊途同归，对同一边值问题只鸵得到同一个自然积分方程． 鉴于以上优点，可以说，自然边界归化在各种边界归化中占有特殊的地位． 2弗行算法的基本理论 五l并行算法的的目标 并行是指有一个以上的事件在同一时刻或同一时间段内发生．从计算复杂性的角度来 看．一个算法的复杂性可表示为空问复杂性和时间复杂性两个方面．并行算法的目标是尽可 能减少时间复杂性．通常这是通过增加空间复杂性(如增加空间的维数及增加处理器的台数) 来实现的．从算法树的结构来看，通常的串行算法树是“深而窄”；而并行算法树的结构是“浅 而宽”，即每个时刻可容纳的计算量相应增加，使整个算法的步数尽可能减少威者说．通 过增加每个时刻步的算法复杂性来减少整体的时间复杂性． 2．2阿络并行稿程的基本模式 从不同的任务组织的角度看，分布式网络并行编程主要有以下两种模式： (1)星形计算模式． 这种计算模式又可进一步分为两类： 负责进程的生成，初始化，收集并显示计算结果，其余的从进程(slave)执行实际计算。 其负载或者由主进程分配，或者由进程本身分配． b．纯结点模式．这时候所有进程都在执行单个的程序．只有少数进程(经常初始时由 人工指定)同时负责非计算的功能． (2)树形(n∞)计算模式． 在这种计算模式中，进程通常是在计算过程中以树形方式动态生成．这是求解组合优化 问题常用的一种方法，是一种构造性的探索算法，主要思想是对解集合反复进行分支，对每 一个分支计算最优解的界．如果该解符合要求，则继续分支以探索更好的解，直