中文

用于模拟大型熔丝网络断裂的高效块循环预处理器

材料科学 2009-11-11 v1 统计力学

摘要

与临界点附近迭代求解器相关的{\it 临界慢化}常常阻碍使用离散晶格网络进行大规模断裂数值模拟。本文提出了一种用于迭代求解器的块循环预处理器,用于模拟无序准脆性材料中使用大型离散晶格网络的渐进断裂。本算法每次迭代的平均计算成本为 O(rslogs)+delopsO(rs log s) + delops,其中刚度矩阵 A{\bf A} 被划分为 rr-by-rr 个块,每个块是 ss-by-ss 矩阵,delopsdelops 表示求解具有 rr-by-rr 稠密矩阵块的分块对角矩阵的运算次数。使用块循环预处理器的该算法比傅里叶加速预条件共轭梯度 (PCG) 算法更快,并缓解了在临界点附近尤为严重的{\it 临界慢化}。使用随机电阻网络的数值结果证实了本算法的效率。

关键词

引用

@article{arxiv.cond-mat/0503722,
  title  = {An Efficient Block Circulant Preconditioner For Simulating Fracture Using Large Fuse Networks},
  author = {Phani Kumar V. V. Nukala and Srdjan Simunovic},
  journal= {arXiv preprint arXiv:cond-mat/0503722},
  year   = {2009}
}

备注

16 pages including 2 figures