欢迎光临《计算力学学报》官方网站!
王计真,陈海波,陈红永.基于分区加速和总体共轭梯度法的耦合界面数据传递问题研究[J].计算力学学报,2015,32(2):280~286
基于分区加速和总体共轭梯度法的耦合界面数据传递问题研究
Research of data transfer on coupling interface based on partition acceleration and global CG algorithm
投稿时间:2014-02-26  修订日期:2014-04-10
DOI:10.7511/jslx201502023
中文关键词:  数据传递  耦合动力学  紧支径向基函数  分区加速  总体共轭梯度迭代方法
英文关键词:data exchanges  coupling dynamics  compactly supported radial basis function  partition acceleration  global conjugate iterative method
基金项目:
作者单位E-mail
王计真 中国科学技术大学 近代力学系 中国科学院材料力学行为与设计重点实验室, 合肥 230027  
陈海波 中国科学技术大学 近代力学系 中国科学院材料力学行为与设计重点实验室, 合肥 230027 hbchen@ustc.edu.cn 
陈红永 中国科学技术大学 近代力学系 中国科学院材料力学行为与设计重点实验室, 合肥 230027  
摘要点击次数: 1327
全文下载次数: 1344
中文摘要:
      对于耦合动力学问题的分析过程,在界面上需频繁进行数据交换。为此,基于紧支径向基函数和多项式基函数推导了界面数据传递的插值算法,给出了传递矩阵的具体形式。通过分析时间复杂度,找出该算法在大节点量时效率不高的原因在于径向基矩阵的构造和传递矩阵的计算。为加快径向基矩阵的构造速度,提出分区加速处理以提高相关节点的搜索效率;为避免传递矩阵求解过程中的求逆运算,将其转化为多右端项的大型稀疏对称线性方程组问题,引入多右端项的总体共轭梯度迭代方法求解,并讨论了初始估计矩阵的选取方法。数值算例结果表明,结合使用分区加速原理和总体共轭梯度迭代方法,可在不损失插值精度的前提下显著提高求解效率。
英文摘要:
      Frequent data exchanges take place on the interface in the solution of coupled problems.For this problem, a kind of interpolation algorithm was derived from the compactly supported radial basis functions and the polynomial basis functions;and the specific form of the transfer matrix was given as well.By analyzing the complexity of the algorithm, it was found that constructing the radial basis matrix and computing the transfer matrix consumed most CPU time, causing computational inefficiency when the number of nodes on the coupling interface was large.The partition acceleration principle was proposed to accelerate the construction of the radial basis matrix by improving the efficiency of searching the related nodes.To avoid matrix inversion, the process of computing the transfer matrix is converted to solve a multi-right-hand sparse symmetric linear system of equations, and the global CG iterative method was introduced to deal with it.Finally how to select the initial guess matrix was discussed.Numerical examples show that a combination of the partition acceleration treatment and the global conjugate gradient iterative method can improve numerial efficiency greatly without accuracy loss.
查看全文  查看/发表评论  下载PDF阅读器
您是第5387618位访问者
版权所有:《计算力学学报》编辑部
本系统由 北京勤云科技发展有限公司设计