研究论文

技术：：块预处理能力及在Navier-Stokes和MHD中的具体应用示例

作者:埃里克·C。赛尔,约翰·N。阴影、和雷蒙德·S。图米纳罗作者信息和声明

SIAM科学计算杂志,体积38,问题5

页S307系列-S331型

https://doi.org/10.1137/15M1017946

出版:2016年10月27日出版历史

摘要

本文描述了Teko的设计，这是一个用于实现高级块预处理程序的面向对象C++库。解释并显示了阐明块预处理库和技术需求的数学设计标准，以激励Teko的结构。例如，Teko的主要设计选择是强烈反映预处理程序的数学陈述，以减少开发负担，并允许专注于数值。解释了其他机制，为开发具有优化生产能力的Teko区块预处理能力提供了途径。最后，对Teko在流体流动和磁流体力学方面的应用进行了演示。除了突出Teko库的特点外，这些新结果还说明了最近预处理技术的发展对高级离散化方法的有效性。

工具书类

[1]

R.Bartlett，Thyra线性算子和向量抽象数值算法开发和互操作的接口和支持软件概述，技术报告SAND2007-5984，桑迪亚国家实验室，2007年。

[2]

M.Benzi、M.Ng、Q.Niu和Z.Wang，不可压缩Navier-Stokes方程的松弛维分解预条件，J.计算。物理。，230（2011年），第6185-6202页。

[3]

M.Benzi和M.A.Olshanskii，基于增广拉格朗日方法的Oseen问题，SIAM J.科学。计算。，28（2006），第2095-2113页。

[4]

M.Benzi和Z.Wang，不可压Navier-Stokes方程修正增广拉格朗日预条件的并行实现，数字。《算法》，64（2013），第73-84页。

[5]

J.Brown、M.G.Knepley、D.A.May、L.C.McInnes和B.Smith，多物理的可组合线性解算器，《2012年第11届并行和分布式计算国际研讨会论文集》，IEEE，华盛顿特区，2012年，第55-62页。

[6]

拉查科，三维粘阻磁流体力学的最优并行全隐式Newton-Krylov解算器，物理。Plasmas，15（2008），056103。

[7]

L.Chacón等人，三维磁流体力学可伸缩并行隐式求解器《物理学杂志》。Conf.序列号。，125 (2008), 012041.

[8]

L.Chacoín和D.A.Knoll，二维高-$\beta$Hall MHD隐式非线性求解器，J.计算。物理。，188（2003），第573--592页。

[9]

L.Chacoín、D.A.Knoll和J.M.Finn，隐式非线性降阻MHD求解器，J.计算。物理。，178（2002），第15-36页。

[10]

E.Chow和M.A.Heroux，面向对象的块预处理框架，ACM变速器。数学。《软件》，24（1998），第159-183页。

[11]

E.C.Cyr、J.N.Shadid和R.S.Tuminaro，Navier-Stokes方程的稳定性和可伸缩块预处理，J.计算。物理。，231（2012），第345-363页。

[12]

E.C.Cyr、J.N.Shadid、R.S.Tuminaro、R.P.Pawlowski和L.ChacoóN，二维不可压缩（约化）电阻磁流体力学的一种新的近似块分解预处理器，SIAM J.科学。计算。，35（2013年），第B701-B730页。

[13]

H.Elman、V.E.Howle、J.Shadid、R.Shuttleworth和R.Tuminaro，基于近似交换子的块预条件器，SIAM J.科学。计算。，27（2006），第1651--1668页。

[14]

H.Elman、V.E.Howle、J.Shadid、R.Shuttleworth和R.Tuminaro，不可压缩Navier-Stokes方程并行块多级预处理的分类与比较，J.计算。物理。，227（2008），第1790--1808页。

[15]

H.Elman、D.Silvester和A.Wathen，有限元和快速迭代解及其在不可压缩流体动力学中的应用，牛津大学出版社，牛津，2005年。

[16]

H.Elman和R.S.Tuminaro，Navier-Stokes方程近似交换子预条件中的边界条件，电子。事务处理。数字。分析。，35（2009年），第257--280页。

[17]

H.Elman、V.E.Howle、J.N.Shadid和R.S.Tuminaro，三维不可压Navier-Stokes方程的并行块多级预处理，J.计算。物理。，187（2003），第504-523页。

[18]

E.Gamma、R.Helm、R.Johnson和J.Vlisseds，设计模式：可重用面向对象软件的元素，Addison Wesley Professional，马萨诸塞州雷丁市，1995年。

[19]

M.W.Gee、C.M.Siefert、J.J.Hu、R.S.Tuminaro和M.G.Sala，ML 5.0平滑聚合用户指南，技术报告SAND2006-2649，桑迪亚国家实验室，2006年。

[20]

M.D.Gunzburger，粘性不可压流动的有限元方法《学术出版社》，波士顿，1989年。

[21]

M.A.Heroux、R.A.Bartlett、V.E.Howle、R.J.Hoekstra、J.J.Hu、T.G.Kolda、R.B.Lehoucq、K.R.Long、R.P.Pawlowski、E.T.Phipps、A.G.Salinger、H.K.Thornquist、R.S.Tuminaro、J.M.Willenbring、A.Williams和K.S.Stanley，Trilinos项目概述，ACM变速器。数学。《软件》，31（2005），第397--423页。

[22]

V.E.Howle、R.C.Kirby、K.Long、B.Brennan和K.Kennedy，Playa：高性能可编程线性代数《科学规划》，20（2012），第257至273页。

[23]

P.Lin、M.Bettencourt、S.Domino、T.Fisher、M.Hoemmen、J.Hu、E.Phipps、A.Prokopenko、S.Rajamanickam、C.Siefert和S.Kennon，第二代Trilinos对低马赫数流体的极值模拟，并行处理Lett。，24 (2014), 1442005.

[24]

P.T.Lin、M.Sala、J.N.Shadid和R.S.Tuminaro，不可压缩流动和输运的全耦合代数多级区域分解预条件的性能，《国际期刊数字》。方法。《工程师》，67（2006），第208-225页。

[25]

R.C.Martin，开闭原理《更多C++宝石》，剑桥大学出版社，纽约，1996年，第97-112页。

[26]

D.A.May和L.Moresi，计算地球动力学中Stokes流问题的预处理迭代方法，物理。地球行星。In.，171（2008），第33-47页。

[27]

B.迈耶，面向对象的软件构造《Prentice-Hall计算机科学国际丛书》，新泽西州上鞍河Prentice-Hall出版社，1988年。

[28]

M.F.Murphy、G.H.Golub和A.J.Wathen，关于不定线性系统预处理的注记，SIAM J.科学。计算。，21（2000），第1969--1972页。

[29]

S.V.Patankar和D.B.Spalding，三维抛物流中热量、质量和动量传递的计算方法《国际传热传质杂志》，第15卷（1972年），第1787-1806页。

[30]

E.G.Phillips、E.C.Cyr和J.N.Shadid，非定常Navier-Stokes的块预条件研究，载于CSRI《2010年夏季论文集》，桑迪亚国家实验室计算机科学研究所，2010年，第47-58页。

[31]

E.G.Phillips、H.C.Elman、E.C.Cyr、J.N.Shadid和R.P.Pawlowski，稳态磁流体力学精确罚公式的块预处理，SIAM J.科学。计算。，36（2014年），第B930--B951页。

[32]

M.Rajan、D.W.Doerfler、P.T.Lin、S.D.Hammond、R.F.Barrett和C.T.Vaughan，新的NNSA Tri-Lab Linux容量群集2具有前所未有的可扩展性和性能，《2012 SC Companion:High Performance Computing，Networking，Storage and Analysis》（SCC），2012年，第417-425页。

[33]

J.N.Shadid、R.P.Pawlowski、J.W.Banks、L.ChacoíN、P.T.Lin和R.S.Tuminaro，基于稳定有限元方法的可扩展全隐式全耦合电阻磁流体力学公式，J.计算。物理。，229（2010），第7649-7671页。

[34]

J.N.Shadid、R.P.Pawlowski、E.C.Cyr、R.S.Tuminaro、L.Chacon和P.D.Weber，具有稳定有限元和完全耦合Newton-Krylov-AMG的可伸缩隐式不可压缩电阻磁流体动力学,计算。应用程序中的方法。机械。《工程师》，304（2016），第1-25页。

建议

二维不可压缩（约化）电阻磁流体力学的一种新的近似块分解预条件

单流体粘电阻MHD模型描述了带电流体在电磁场影响下的动力学。该模型具有强耦合、高度非线性的特点，其物理机制跨越范围广。。。
阅读更多信息
固定MHD精确惩罚公式的块预条件

磁流体动力学（MHD）方程用于模拟液态金属等应用中的导电流体和等离子体。这个非自伴非线性偏微分方程组耦合了流体的Navier-Stokes方程还有。。。
阅读更多信息
关于非对称问题的固定点、Krylov和$2乘2$块预条件

具有$2乘以2$块结构的矩阵的解产生于计算数学的许多领域，如PDE基于混合有限元方法的离散化，约束优化问题，或隐式或稳态处理。。。
阅读更多信息

评论

信息和贡献者

问询处

发布于

封面图片SIAM科学计算杂志

SIAM科学计算杂志第38卷第5期

^†两个主题专题：CSE软件和CSE大数据

2016

1789页

国际标准编号：1064-8275

内政部：10.1137/sjoce3.38.5

©2016，工业和应用数学学会。

出版商

工业和应用数学学会

美国

出版历史

出版：2016年10月27日

作者标记

作者标记

限定符

研究文章

贡献者

其他指标

查看文章指标

文献计量学和引文

文献计量学

文章指标

0
引文总数
0
总下载次数

下载次数（过去12个月）0
下载次数（最近6周）0

其他指标

查看作者指标

引文

视图选项

查看选项

获取访问权限

登录选项

检查您是否可以通过登录凭据或您的机构访问本文。

完全访问权限

获取此出版物

媒体

数字

其他

桌子

查看问题目录