用于模拟大型熔丝网络断裂的高效块循环预处理器
材料科学
2009-11-11 v1 统计力学
摘要
与临界点附近迭代求解器相关的{\it 临界慢化}常常阻碍使用离散晶格网络进行大规模断裂数值模拟。本文提出了一种用于迭代求解器的块循环预处理器,用于模拟无序准脆性材料中使用大型离散晶格网络的渐进断裂。本算法每次迭代的平均计算成本为 ,其中刚度矩阵 被划分为 -by- 个块,每个块是 -by- 矩阵, 表示求解具有 -by- 稠密矩阵块的分块对角矩阵的运算次数。使用块循环预处理器的该算法比傅里叶加速预条件共轭梯度 (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