×
Supriyo Ghosh;克里斯托弗·纽曼。;玛丽安·弗朗索瓦（Marianne M.Francois）。

图萨（Tusas）：耦合相场方程的完全隐式并行方法。 （英语） Zbl 07516818号

J.计算。物理学。 448，文章ID 110734，19 p.（2022）.
总结：我们开发了一种全耦合、全隐式方法，用于金属和合金凝固的相场建模。纯金属和金属合金凝固的预测模拟仍然是材料科学领域的一个重大挑战，因为凝固过程中微观结构的形成对固体材料的性能和性能起着至关重要的作用。我们的模拟方法包括在微尺度上对全耦合非线性偏微分方程组进行有限元空间离散，并使用预处理的无雅可比牛顿-克利洛夫方法对其进行隐式时间处理。由于基于代数多重网格和块分解的有效预处理策略，该方法在算法上是可扩展的，并且是有效的。我们在开源中实现了这种方法图萨（Tusas）框架，这是一个用C++开发的通用、灵活的工具，用于求解非线性偏微分方程耦合系统。我们从算法可扩展性和效率方面分析了该方法的性能，而图萨（Tusas）在新兴异构体系结构上的并行可伸缩性和效率方面进行了介绍。我们证明，现代算法、离散化、计算科学和异构硬件为添加剂制造过程中微观结构演变的预测相场模拟提供了一条稳健的途径。

MSC公司：

65传真 数值线性代数
65华夏 非线性代数或超越方程
7.6亿 流体力学基本方法

关键词：

隐式方法;JFNK方法;预处理;相场模拟;凝固

软件：

磁粉探伤;净CDF;阿美索斯2;高密度纤维;增强C++库;阿兹特科;炒作;驼鹿;CUDA公司;日晷;菲比;埃克索斯二世;速度3D;最小均方误差;迈克尔;相位传感器;图萨（Tusas）;特里利诺斯;促进;C制动器;毫升;梅鲁;科科斯;硝基苯砜;拜洛什
PDF格式BibTeX公司 XML格式引用
\textit{S.Ghosh}等人，J.Compute。物理学。448，文章ID 110734，19 p.（2022；Zbl 07516818）
全文： 内政部 arXiv公司

参考文献：

[1] 拉帕兹，M。；Dantzig，J.A.，《凝固，工程科学》（2009），EPFL出版社：洛桑EPFL·Zbl 1194.76001号
[2] Kurz，W。；Fisher，D.J.，《固化基础》（1998），Trans-Tech出版物：苏黎世Trans-Tec出版物
[3] 弗朗索瓦，M.M。；Sun，A。；King，W.E。；新泽西州亨森。；托雷特，D。；Bronkhorst，C.A。；卡尔森，N.N。；纽曼，C.K。；豪特，T.S。；Bakosi，J.，《金属添加剂制造工艺的建模：挑战与机遇》，Curr。操作。固态材料。科学。，21 (2017)
[4] Ghosh，S.，《镍高温合金激光添加剂制造过程中凝固的预测建模：最新发展，未来方向》，Mater。Res.Express，5，第012001条pp.（2018）
[5] Keller，T。；林德沃尔，G。；Ghosh，S。；马，L。；莱恩，B。；张，F。；卡特纳，U.R。；Lass，E.A。；海格尔，J.C。；Idell，Y。；威廉姆斯，M.E。；艾伦·A·J。；盖耶，J.E。；Levine，L.E.，有限元、相场和基于CALPHAD的方法在镍基高温合金增材制造中的应用，材料学报。，139, 244-253 (2017)
[6] 纪毅。；Chen，L。；Chen，L.-Q.，使用相场建模理解金属合金添加剂制造过程中的微观结构演变，（添加剂制造的热-机械建模（2018），Elsevier:Elsevier NY），93-116
[7] 阿查里亚，R。；沙龙，J.A。；Staroselsky，A.，《激光粉末床熔化过程中微观结构的预测》，《材料学报》。，124, 360-371 (2017)
[8] 昆丁，J。；Mushongera，L。；Emmerich，H.，《Inconel 718高温合金快速凝固过程中微观组织形成的相场模拟》，《材料学报》。，95, 343-356 (2015)
[9] Ghosh，S。；马，L。；Ofori-Opoku，N。；Guyer，J.E.，关于激光沉积Ni-Nb合金中胞状枝晶的一次间距和微观偏析，模型。模拟。马特。科学。工程，25，第065002条pp.（2017）
[10] Ghosh，S。；马，L。；莱文，L.E。；瑞克·R·E。；斯托特，M.R。；海格尔，J.C。；Guyer，J.E.，Inconel 625中的单道熔池测量和微观结构，JOM，1-6（2018）
[11] Ghosh，S。；McReynolds，K。；盖耶，J.E。；Banerjee，D.，Ti-6Al-4V激光沉积过程中温度、应力和微观结构场的模拟，模型。模拟。马特。科学。Eng.，26，文章075005第（2018）页
[12] 卡拉亚吉兹，K。；约翰逊，L。；种子，R。；阿塔里五世。；张，B。；黄，X。；Ghosh，S。；Duong，T。；卡拉曼，I。；Elwany，A。；Arroyave，R.，附加制造条件下Ni-Nb的有限界面耗散相场建模，材料学报。，185, 320-339 (2020)
[13] Chen，L.Q.，微观结构演化的相场模型，年。修订版材料。研究，32，113-140（2002）
[14] 博廷格，W.J。；沃伦，J.A。；贝克曼，C。；Karma，A.，凝固相场模拟，阿努。修订版材料。决议，32，163-194（2002）
[15] 北卡罗来纳州莫兰斯。；布兰帕因，B。；Wollants，P.，《微观结构演化的相场建模简介》，卡尔帕德，32，268-294（2008）
[16] 斯坦巴赫，I.，材料科学中的相场模型，模型。模拟。马特。科学。工程，17，第073001条pp.（2009）
[17] Lambert，J.，《常微分系统的数值方法：初值问题》（1991），John Wiley&Sons:John Willey&Sons Hoboken，NJ·Zbl 0745.65049号
[18] 布伦纳，S。；Scott，R.，《有限元方法的数学理论》，第15卷（2007），Springer科学与商业媒体：Springer Science&Business Media New York
[19] Braess，D.，《有限元：固体力学中的理论、快速求解和应用》（2001），剑桥大学出版社：剑桥大学出版社·Zbl 0976.65099号
[20] Hübner，B。；Walhorn，E。；Dinkler，D.，使用时空有限元进行流体-结构相互作用的整体方法，计算。方法应用。机械。工程，1932087-2104（2004）·兹比尔1067.74575
[21] Michler，C。；赫尔肖夫，S。；Van Brummelen，E。；De Borst，R.，流体-结构相互作用的整体方法，计算。流体，33839-848（2004）·Zbl 1053.76042号
[22] 马蒂斯·H·G。；Steindorf，J.，流体-结构相互作用的分区强耦合算法，计算。结构。，81005-812（2003年）
[23] Erbts，P。；哈特曼，S。；Düster，A.，《电热-机械问题的分区求解方法》，Arch。申请。机械。，85, 1075-1101 (2015) ·Zbl 1347.74025号
[24] 乔希，V。；Jaiman，R.K.，保守和正保持Allen-Cahn相场模型的自适应变分程序，J.Compute。物理。，366, 478-504 (2018) ·Zbl 1406.76049号
[25] 乔希，V。；Jaiman，R.K.，两相流相场建模的保正保守变分格式，J.Compute。物理。，360, 137-166 (2018) ·Zbl 1391.76651号
[26] Knoll，D.A。；Keyes，D.E.，《无雅可比牛顿-克利洛夫方法：方法和应用调查》，J.Compute。物理。，193, 357-397 (2004) ·Zbl 1036.65045号
[27] Kelley，C.T.，《线性和非线性方程的迭代方法》（1995），SIAM:SIAM Philadephia，PA·Zbl 0832.65046号
[28] Knoll，D.A。；Rider，W.J.，《多重网格预处理Newton-Krylov方法》，SIAM J.Sci。计算。，21691-710（2000年）·Zbl 0952.65102号
[29] Saad，Y.，《稀疏线性系统的迭代方法》，《计算机科学中的PWS系列》（1995），PWS出版公司：PWS出版社，马萨诸塞州波士顿·Zbl 1002.65042号
[30] 乔希，V。；Jaiman，R.K.，两相流-结构相互作用耦合分析的混合变分Allen-Cahn/ALE格式，国际期刊数值。方法工程，117，405-429（2019）
[31] Knoll，D.A。；Mousseau，V.A。；Chacón，L。；Reisner，J.M.，《刚性波系统精确时间积分的无雅可比牛顿-克利洛夫方法》，J.Sci。计算。，25, 213-230 (2005) ·Zbl 1203.65071号
[32] 帕克，H。；努尔加列夫，R.R。；马提诺，R.C。；Knoll，D.A.，《关于基于物理的Navier-Stokes方程预处理》，J.Compute。物理。，228, 9131-9146 (2010) ·Zbl 1395.65028号
[33] Mousseau，V.A。；Knoll，D.A。；Reisner，J.M.，带科里奥利力的二维浅水方程的隐式非线性一致方法，Mon。Weather Rev.，1302611-2625（2002年）
[34] Knoll，D.A。；Chac'on，L。；Margolin，L.G。；Mousseau，V.A.，《关于多时间尺度系统时间积分的平衡近似》，J.Compute。物理。，185, 583-611 (2003) ·Zbl 1047.76074号
[35] 格罗普，W。；Lusk，大肠杆菌。；Skjellum，A.，《使用MPI：具有消息传递接口的可移植并行编程》，第1卷（1999年），麻省理工学院出版社
[36] 达贡，L。；Menon，R.，OpenMP：用于共享内存编程的行业标准API，IEEE Compute。科学。工程师，5，46-55（1998）
[37] NVIDIA Corporation，NVIDIA-CUDA C编程指南，2010年，3.2版。
[38] 雅各布森，D。；Thibault，J。；Senocak，I.，多GPU集群上大规模并行不可压缩流计算的MPI-CUDA实现，（第48届AIAA航空航天科学会议，包括新视野论坛和航空航天博览会（2010）），522
[39] Glaser，J。；Nguyen，T.D。；安德森，J.A。；Lui，P。；斯皮加，F。；Millan，J.A。；莫尔斯特区。；Glotzer，S.C.，GPU上通用分子动力学模拟的强缩放，计算。物理学。社区。，192, 97-107 (2015)
[40] Komatitsch，D。；Erlebacher，G。；Göddeke，D。；Michéa，D.，大型GPU集群上MPI高阶有限元地震波传播建模，J.Compute。物理。，229, 7692-7714 (2010) ·兹比尔1194.86019
[41] 杨，C。；徐，Q。；Liu，B.，GPU加速镍基高温合金中枝晶生长的三维相场模拟，Comput。马特。科学。，136, 133-143 (2017)
[42] Yamanaka，A。；青木，T。；小川，S。；Takaki，T.，二元合金枝晶凝固的GPU加速相场模拟，J.Cryst。增长，31840-45（2011）
[43] Sakane，S。；Takaki，T。；Ohno，M。；Shimokawabe，T。；Aoki，T.，二元合金定向凝固过程中枝晶竞争生长的GPU加速三维相场模拟，IOP会议系列：材料科学与工程，第84卷，012063（2015），IOP出版社
[44] D.加斯顿。；纽曼，C。；Hansen，G。；Lebrun-Grandie，D.，MOOSE:非线性方程耦合系统的并行计算框架，Nucl。工程设计。，239, 1768-1778 (2009)
[45] （2019），PACE3D（3D晶体演化并行算法）
[46] MMSP，中尺度微观结构模拟项目（2014）
[47] MICRESS，微观结构演化模拟软件（2019）
[48] OpenPhase（2007）
[49] （2019），PhasePot微观结构模拟工具
[50] 盖耶，J.E。；惠勒博士。；Warren，J.A.，FiPy：Python偏微分方程，计算。科学。工程，11，6-15（2009）
[51] 德维特，S。；Rudraraju，S。；蒙蒂尔，D。；安德鲁斯·W·B。；Thornton，K.，PRISMS-PF：无矩阵有限元法相场建模的一般框架，npj Compute。材料。，6（2020）
[52] 多尔，M.R。；法特伯特，J.-L。；Wickett，M。；Belak，J。；Turchi，P.，多晶材料相场模型求解的数值算法，J.Compute。物理。，229, 626-641 (2010) ·Zbl 1183.82084号
[53] 拉德哈克里希南，B。；戈蒂，S.B。；特纳，J.A。；阿查里亚，R。；沙龙，J.A。；Staroselsky，A。；El-Wardany，T.，在激光粉末床熔化期间使用替代Ni-Fe-Nb合金对IN718微观结构演变的相场模拟，金属，9，14（2019）
[54] 朱，C。；刘杰。；朱，M。；Feng，L.，二元合金中的Multi-GPU混合编程加速三维相场模型，AIP Adv.，8，文章035312 pp.（2018）
[55] Shimokawabe，T。；青木，T。；Takaki，T。；Endo，T。；Yamanaka，A。；Maruyama，N。；努卡达，A。；Matsuoka，S.，TSUBAME 2.0超级计算机上枝晶凝固的Peta-scale相场模拟，（2011年高性能计算、网络、存储和分析国际会议论文集（2011），ACM），3
[56] Caginalp，G.，Stefan和Hele-Shaw型模型作为相场方程的渐近极限，Phys。版本A，39，5887（1989）·Zbl 1027.80505号
[57] Sekerka，R.F.，《形态学：从尖锐界面到相场模型》，J.Cryst。增长，264530-540（2004）
[58] Arfken，G.B。；韦伯，H.J。；Harris，F.E.，《物理学家的数学方法》（2013年），爱思唯尔出版社：纽约·Zbl 1239.00005
[59] 尤因，G.M.，《应用变分法》（1969年），W.W.诺顿公司：W.W.Norton公司，纽约·兹比尔0198.44501
[60] 艾伦，S.M。；Cahn，J.W.，《反相边界运动微观理论及其在反相畴粗化中的应用》，《金属学报》。，27, 1085-1095 (1979)
[61] 卡恩，J.W。；Hilliard，J.E.，非均匀系统的自由能。I.界面自由能，J.Chem。物理。，28, 258-267 (1958) ·Zbl 1431.35066号
[62] 丹尼斯，J.E。；Schnabel，R.B.，《无约束优化和非线性方程的数值方法》（1996），SIAM:SIAM Philadelphia，PA·Zbl 0847.65038号
[63] Pawlowski，R.P。；沙迪德，J.N。；Simonis，J.P。；Walker，H.F.，Newton-Krylov方法的全球化技术及其在Navier-Stokes方程完全耦合解中的应用，SIAM Rev.，48，700-721（2006）·Zbl 1110.65039号
[64] 佩妮斯，M。；Walker，H.F.，NITSOL:非线性系统的牛顿迭代求解器，SIAM J.Sci。计算。，19, 302-318 (1998) ·Zbl 0916.65049号
[65] Brown，P.N。；Saad，Y.，非线性方程组的混合Krylov方法，SIAM J.Sci。统计计算。，11, 450-481 (1990) ·Zbl 0708.65049号
[66] Dembo，R。；艾森斯塔特，S.C。；Steihaug，T.，不精确牛顿方法，SIAM J.Numer。分析。，19, 400-408 (1982) ·Zbl 0478.65030号
[67] 艾森斯塔特，S.C。；Walker，H.F.，在不精确牛顿方法中选择强迫项，SIAM J.Sci。计算。，17, 16-32 (1996) ·Zbl 0845.65021号
[68] Chacón，L。；Knoll，D.A。；Finn，J.M.，《隐式、非线性降低电阻MHD求解器》，J.Compute。物理。，178, 15-36 (2002) ·Zbl 1139.76328号
[69] 北普罗瓦塔斯。；Elder，K.，《材料科学与工程中的相场方法》（2011），John Wiley&Sons：新泽西州John Willey&Sons Hoboken
[70] Briggs，W.L。；McCormick，S.F.，《多重网格教程》，第72卷（2000），SIAM：宾夕法尼亚州费城SIAM·Zbl 0958.65128号
[71] Hackbusch，W.，《多网格方法与应用》，第4卷（2013年），Springer科学与商业媒体：Springer科技与商业媒体纽约
[72] 美国特罗滕贝格。；Oosterlee，C.W。；Schüller Multigrid，A.（2000），学术出版社：纽约学术出版社
[73] 吉，M。；西弗特，C。；胡，J。；杜米纳罗，R。；Sala，M.，ML 5.0平滑聚合用户指南（2006），桑迪亚国家实验室，技术报告SAND2006-2649
[74] Berger-Vergiat，L。；Glusa，C.A。；胡建杰。；Mayr，M。；普罗科彭科，A。；西弗特，C.M。；杜米纳罗，R.S。；Wiesner，T.A.，MueLu用户指南（2019），桑迪亚国家实验室，技术报告SAND2019-0537
[75] Berger-Vergiat，L。；Glusa，C.A。；胡建杰。；Mayr，M。；普罗科彭科，A。；西弗特，C.M。；杜米纳罗，R.S。；Wiesner，T.A.，MueLu多重网格框架（2019年）
[76] Heroux，M.A.，Trilinos项目概述，ACM Trans。数学。软质。，31, 397-423 (2005) ·Zbl 1136.65354号
[77] 万亿（2019）
[78] 图萨（2020）
[79] 出埃及记II（2007）
[80] 斯科夫，L。；Yarberry，V.，EXODUS II：有限元数据模型（1996），Sandia国家实验室，技术报告SAND92-2137
[81] NOX和LOCA，面向对象的非线性解算器和延续程序包（2010年）
[82] Bavier，E。；霍姆曼，M。；拉贾马尼卡姆，S。；Thornquist，H.、Amesos2和Belos:大型稀疏线性系统的直接和迭代求解器，科学。程序。，20, 241-255 (2012)
[83] M.A.Heroux，AztecOO用户指南。
[84] Rythmos（2019年）
[85] D.加斯顿。；纽曼，C。；Hansen，G.，MOOSE:非线性方程耦合系统的并行计算框架，（美国核学会2009年数学、计算方法和反应堆物理进展国际会议。美国核学会09年数学、计算机方法和反应堆物理学进展国际会议，纽约萨拉托加斯普林斯（2009））
[86] CMake（2019年）
[87] Boost C++库（2019年）
[88] 分层数据格式，第5版（1997年）
[89] 雷沃，R。；Davis，G.，《数据管理：NetCDF:科学数据访问接口》，IEEE Comput。图表。申请。，10, 76-82 (1990)
[90] （2020），Kokkos C++性能可移植编程生态系统：编程模型-并行执行和内存抽象
[91] 爱德华兹，H.C。；Trott，C.R。；Sunderland，D.，Kokkos：通过多态内存访问模式实现多核性能可移植性，J.Parallel Distrib.Comput。，74, 3202-3216 (2014)
[92] R·D·法尔古特。；Yang，U.M.，HYPRE:高性能预处理程序库，（国际计算科学会议（3）（2002）），632-641·Zbl 1056.65046号
[93] 辛德马什，A.C。；Brown，P.N。；格兰特，K.E。；Lee，S.L。；塞尔维亚人，R。；Shumaker，D.E。；Woodward，C.S.，SUNDIALS：非线性和微分/代数方程求解器套件，ACM Trans。数学。软质。，31, 363-396 (2005) ·Zbl 1136.65329号
[94] Jokisaari，A.M。；沃里斯，P.W。；盖耶，J.E。；沃伦，J.A。；Heinonen，O.G.，枝晶生长和线性弹性的相场基准问题，计算。马特。科学。，149, 336-347 (2018)
[95] 康明斯，S.J。；Quirk，J.J。；Kothe，D.B.，《微观结构凝固建模相场技术的探索》，（2002年美国机械工程师协会凝固微观结构发展会议记录（2002年），美国机械工程师学会：美国机械工程师委员会新奥尔良）
[96] Kessler博士。；Levine，H.，《树突生长中的速度选择》，Phys。B版，33，7867（1986）
[97] 埃切巴里亚，B。；Folch，R。；Karma，A。；Plapp，M.，合金凝固的定量相场模型，物理。E版，70，第061604条，pp.（2004）
[98] Farzadi，A。；Do-Quang，M。；Serajzadeh，S。；科卡比，A。；Amberg，G.，Al-Cu合金焊缝凝固微观结构的相场模拟，模型。模拟。马特。科学。工程，16，1-18（2008）
[99] Karma，A.，合金凝固定量建模的相场公式，Phys。修订稿。，87，第115701条pp.（2001）
[100] Ghosh，S。；奥福里·奥波库，北。；Guyer，J.E.，镍基高温合金激光粉末沉积过程中γ-Ni细胞生长的模拟和分析，计算。马特。科学。，144, 256-264 (2018)
[101] Ghosh，S。；斯托特，M.R。；莱文，L.E。；Guyer，J.E.，《激光沉积Ni-矩阵微结构中富铌液滴的形成》，Scr。材料。，146, 36-40 (2018)
[102] 峰会（2018）
[103] Sierra（2018年）
[104] 刘，H。；杨，B。；Chen，Z.，在NVIDIA GPU上加速代数多重网格求解器，计算。数学。申请。，70, 1162-1181 (2015) ·Zbl 1443.65046号
[105] Haase，G。；利伯曼，M。；C.C.道格拉斯。；Plank，G.，图形处理单元上的并行代数多重网格求解器，（高性能计算与应用（2010），施普林格：施普林格纽约），38-47
[106] 里希特，C。；Schöps，S。；Clemens，M.，离散椭圆场问题代数多重网格预条件的GPU加速，IEEE Trans。马格纳。，50, 461-464 (2014)
[107] 甘达姆，R。；Esler，K。；Zhang，Y.，GPU加速聚合代数多重网格方法，计算。数学。申请。，68, 1151-1160 (2014) ·Zbl 1367.65049号
[108] 杜米纳罗，R.S。；沙迪德，J.N。；Cyr，E.C。；林，P。；Pawlowski，R.P。；菲利普斯，E.G。；Wiesner，T。；阿德勒，J。；Benson，T。；Chacon，L。；法雷尔，P。；Rappaport，A。；麦克拉克伦，S。；胡建杰。；Berger-Vergiat，L。；Glusa，C.A。；Siefert，C.，可伸缩预处理程序，用于改进磁流体力学应用的解决时间（2019年）
[109] 科特，D。；Lee，S。；Qualters，I.，美国Exascale计算，计算机。科学。工程师，21，17-29（2018）
[110] Alexander，F.，《Exascale应用：游戏中的皮肤》，Philos。事务处理。R.Soc.A，378（2020年）·Zbl 1462.65230号
此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配，并且可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。