希拉·马查特;杰罗姆·卡雷鲁 平衡地球化学计算线性求解器的比较。 （英语） Zbl 1415.86001号 计算。地质科学。 21，第1期，131-150（2017）. 小结：平衡化学计算和反应迁移建模需要在非常特殊的条件下大量使用线性求解器。待求解的系统很小或很小（（4乘以4）到（20乘以20），有时更大），并且条件很差（条件数高达10^{100}）。这些特定条件从未从线性解算器的鲁棒性、准确性和效率方面进行过研究。在这项工作中，我们提出了要求解的线性系统的特殊性。使用一组化学系统比较了几种直接和迭代求解器，包括或不包括矿物物种的形成。我们证明了直接求解器和迭代求解器可以用于这些问题，并提出了改进化学求解器的计算键。 引用于1文件 MSC公司： 86-08 地球物理问题的计算方法 65平方英尺 线性系统和矩阵反演的直接数值方法 65层35 矩阵范数、条件、缩放的数值计算 80A32型 化学反应流 86A05型 水文学、水文学、海洋学 关键词：地球化学模拟;瞬时平衡化学;线性系统反演;线性解算器;小矩阵;ill-条件矩阵;Newton-Raphson算法 软件：微型SQL;斯佩西;LAPACK公司;高铁;MoMaS公司;帕迪索;连接器;PhreeqcRM公司;FITEQL公司;MINEQL公司;WATSPEC公司;水文地球化学;PHREEQC公司;FEFLOW公司 PDF格式BibTeX公司 XML格式引用 \textit｛H.Machat｝和\textit｛J.Carrayrou｝，计算机。地质科学。21，第1号，131--150（2017；Zbl 1415.86001） 全文： 内政部 哈尔 参考文献： [1] Walter，A.L.等人：地下水中多组分反应迁移的建模。2.受酸性矿山尾矿排放影响的含水层中的金属流动性。水资源。第30（11）号决议，3149-3158（1994）·doi:10.1029/94WR00954 [2] Arora，B.等人：湖泊沉积物中重金属循环的反应迁移基准计算地球科学（2014）·Zbl 1425.76235号 [3] 温特，L。；Leclercq，S。；Lee，J.，《评估核玻璃在粘土地下处置中二氧化硅控制过程的耐久性》（2005年），根特；比利时 [4] Hoteit，H.、Ackerer，P.、Mose，R.：核废料处理模拟：耦合测试案例。计算。地质科学。8(2), 99-124 (2004) ·Zbl 1060.86005号 ·doi:10.1023/B:COMG.0000035074.37722.71 [5] 汤普森，A.F.B.等人：关于地下核试验对地下水污染的评估。环境。地质。42(2-3), 235-247 (2002) ·doi:10.1007/s00254-001-0493-8 [6] Andre，L.等人：二氧化碳注入Dogger含水层碳酸盐岩储层期间，超临界CO2-液体界面上流体-化学相互作用的数值模拟（法国巴黎盆地）。能量转换。管理。48(6), 1782-1797 (2007) ·doi:10.1016/j.enconman.2007.01.006 [7] Kang，Q.等人：地质CO2封存中反应性运移的孔隙尺度模型。运输。多孔介质82（1），197-213（2010）·doi:10.1007/s11242-009-9443-9 [8] Navarre-Sitchler，A.K.等人：利用高性能计算阐明浅层非均质含水层对CO2泄漏的地球化学响应：对CO2封存监测的影响。高级水资源。53(0), 45-55 (2013) ·doi:10.1016/j.advwatres.2012.10.05 [9] Pruess，K.等人：代码比对建立了CO2地质处置数值模拟模型的可信度。能源29（9-10），1431-1444（2004）·doi:10.1016/j.energy.2004.03.077 [10] Regnault，O.，et al.：练习曲实验二氧化碳超临界与矿物相的反应。二氧化碳封存的地质意义。Compt.公司。伦德斯地质。337（15），1331-1339（2005）·doi:10.1016/j.crte.2005.07.012 [11] Valocchi，A.J.，Street，R.L.，Roberts，P.V.：地下水中离子交换溶质的迁移：色谱理论和现场模拟。水资源。第17号决议，1517-1527（1981）·doi:10.1029/WR017i005p01517 [12] Lichtner，P.C.：水热系统中同步化学反应和质量传输的连续体模型。地理学。科斯莫辛。法案49（3），779-800（1985）·doi:10.1016/0016-7037（85）90172-3 [13] Appelo，C.A.J.：水文地球化学传输建模。继续。信息-通信。水文学。《TNO研究》43，81-104（1990） [14] Yeh，G.T.，Tripathi，V.S.：对反应性多化学成分水文地球化学迁移模型最新发展的批判性评估。水资源。第25号决议，第93-108号决议（1989年）·doi:10.1029/WR025i001p00093 [15] Carrayrou，J.等人：模拟强非线性和异质反应迁移问题的数值方法比较——MoMaS基准案例。计算地球科学14（3），483-502（2010）·Zbl 1426.76723号 ·doi:10.1007/s10596-010-9178-2 [16] Hammond，G.E.，Valocchi，A.J.，Lichtner，P.C.：使用无Jacobian-Free Newton Krylov和算子分裂预处理在并行计算机上模拟多组分反应迁移。摘自：Hassanizadeh，S.M.（编辑）《水科学的发展，水资源的计算方法》，《第十四届水资源计算方法国际会议论文集》（CMWR XIV），第727-734页。爱思唯尔（2002） [17] Carrayrou，J.，Mosé，R.，Behra，P.：求解热力学化学的新高效算法。AIChE J.48（4），894-904（2002）·数字对象标识代码：10.1002/aic.690480423 [18] Amir，L.，Kern，M.：非均匀多孔介质中耦合传输与化学的全局方法。计算。地质科学。14(3), 465-481 (2010) ·Zbl 1425.76231号 ·数字对象标识代码：10.1007/s10596-009-9162-x [19] Quarteroni，A.，Sacco，R.，Saleri，F.：数值数学。收件人：Marsden，J.E.，Sirovich，L.，Antman。S.S.（编辑）应用数学课文。第2版。斯普林格，海德堡（2007）·Zbl 1136.65001号 [20] Axelsson，O.等人：直接解和不完全因式分解预处理共轭梯度法。STW报告中一组线性弹性基准问题代数求解方法的比较。2000年，荷兰奈梅亨天主教大学数学系。第1-36页 [21] Barrett，R.、Berry，M.、Chan，T.F.、Demmel，J.、Donato，J.和Dongarra，J.，Eijkhout，V.、Pozo，R.，Romine，C.、Van Der Vorst，H.：线性系统解的模板：迭代方法的构建块，第二版。SIAM，费城（1994）·Zbl 0814.65030号 ·doi:10.1137/1.9781611971538 [22] Gould，N.I.M.，Hu，Y.，Scott，J.A.：大型稀疏对称线性方程组解的稀疏直接解算器的数值评估。2005年，研究委员会中央实验室委员会·兹比尔1365.65129 [23] Allaire，G.，Kaber，S.M.In:Marsden，J.E.，Sirovich，L.，Antman，S.S.（编辑）：数值线性代数。应用数学课文。施普林格，纽约（2008）·Zbl 1135.65014号 [24] Baldwin，C.等人：2D辐射流体动力学代码中的迭代线性解算器：方法和性能。J.计算。物理。154(1), 1-40 (1999) ·兹比尔0964.76036 ·doi:10.1006/jcph.1999.6290 [25] Chao B.T.，L.H.L.，Scott，E.J.：《关于通过机器计算求解病态、联立、线性代数方程》，载于《伊利诺伊大学公报》。1961年，伊利诺伊大学·Zbl 1060.86005号 [26] Hadjidimos，A.：连续超松弛（SOR）和相关方法。J.计算。申请。数学。123（1-2）、177-199（2000）·Zbl 0965.65052号 ·doi:10.1016/S0377-0427（00）00403-9 [27] Kalambi，I.B.：线性方程解的三种迭代方法的比较。J.应用。科学。环境。管理。12(4), 53-55 (2008) [28] Klisinski，M.，Runesson，K.：改进的有限元计算对称和非对称解算器。高级工程师软件。18(1), 41-51 (1993) ·doi:10.1016/0965-9978（93）90006-F [29] Schenk，O.，Gartner，K.：用PARDISO求解非对称稀疏线性方程组。未来。通用。计算。系统。20(3), 475-487 (2004) ·doi:10.1016/j.future.2003.07.011 [30] Xue，X.J.等：求解病态线性代数系统的直接算法。JCPDS-国际中心衍射。数据高级X射线分析。42, 629-633 (2000) [31] Pyzara，A.，Bylina，B.，Bylina，J.：矩阵条件数对迭代方法收敛性的影响（2011）·Zbl 1170.65022号 [32] Hoffmann，J.、Krasöutle，S.、Knabner，P.：基于简化方案的多孔介质反应输运问题的并行全局隐式二维求解器及其在MoMaS基准问题中的应用。计算。地质科学。14(3), 421-433 (2010) ·Zbl 1425.76240号 ·doi:10.1007/s10596-009-9173-7 [33] Soleymani，F.：求解病态线性系统的新方法。Opuscula数学。33(2), 337-344 (2013) ·Zbl 1284.65057号 ·doi:10.7494/OpMath.2013.33.237 [34] Morel，F.，Morgan，J.：计算水化学体系平衡的数值方法。环境。科学。Technol公司。6(1), 58-67 (1972) ·doi:10.1021/es60060a006 [35] 莫雷尔，F.M.M.：《水化学原理》。Wiley Interscience，纽约（1983年） [36] De Windt，L.等人：应用于UO2氧化溶解和铀迁移的反应迁移模型的相互比较。J.康塔姆。水文学。61(1-4), 303-312 (2003) ·doi:10.1016/S0169-7722（02）00127-4 [37] Jauzein，M.等人：用于模拟多孔介质中传输的灵活计算机代码：冲击。Geoderma 44（2-3），95-113（1989）·doi:10.1016/0016-7061（89）90021-9 [38] Parkhurst，D.L.，Appelo，C.A.J.：PHREEQC用户指南（第2版）——一种用于物种形成、批处理反应、一维迁移和地球化学反演计算的计算机程序。水资源。投资。，编辑。1999年：丹佛。第312页·兹比尔1002.65039 [39] 范德利，J.：国际象棋另一种物种形成和表面络合计算机代码。E.d.M.d.巴黎，编辑。1993年：枫丹白露。第85页·Zbl 1276.65020号 [40] Westall，J.C.：MICROQL:BASIC中的化学平衡程序。BASIC中吸附平衡的计算。S.F.I.o.T.EAWAG，编辑。1979年：杜班多夫，第42页·Zbl 0965.65051号 [41] Westall，J.C.：FITEQL 2.1版。1982年：科瓦利斯 [42] Westall，J.C.，Zachary，J.L.，Morel，F.M.M.：MINEQL：计算水系统化学平衡组成的计算机程序。R.M.P.实验室，编辑。1976年：剑桥。第91页 [43] Walter，L.J.，Wolery，T.J.：用于计算化学物种平衡分布的单调序列算法和FORTRAN IV程序。计算。地质科学。1, 57-63 (1975) ·doi:10.1016/0098-3004（75）90005-9 [44] Wigley，T.M.L.：WATSPEC：测定水溶液平衡形态的计算机程序。B.G.R.G.技术公告。，编辑。1977年，第49页 [45] Jennings，A.A.、Kirkner，D.J.、Theis，T.L.：地下水质量模型中的多组分平衡化学。水资源。第18号决议，1089-1096（1982）·doi:10.1029/WR018i004p01089 [46] Cederberg，A.，Street，R.L.，Leckie，J.O.：多组分系统的地下水质量传输和平衡化学模型。Water Resour Res.211095-1104（1985）·doi:10.1029/WR021i008p01095 [47] Yeh，G.T.，Tripathi，V.S.：模拟活性多物种组分运输的模型：模型开发和演示。水资源。第27（12）号决议，3075-3094（1991）·doi:10.1029/91WR02028 [48] Carrayrou，J.：使用SPECY为MoMaS的反应迁移基准寻找一些参考解决方案。计算。地质科学。14(3), 393-403 (2010) ·Zbl 1425.76235号 ·doi:10.1007/s10596-009-9161-y [49] Brassard，P.，Bodurtha，P.：化学物种形成问题的可行集。计算。地质科学。26(3), 277-291 (2000) ·doi:10.1016/S0098-3004（99）00137-5 [50] Carrayrou，J.，Kern，M.，Knabner，P.：MoMaS的反应性运输基准。计算。地质科学。14（3），385-392（2010）·兹比尔1425.76236 ·doi:10.1007/s10596-009-9157-7 [51] Fendorf，S.E.，Li，G.：亚铁还原铬酸盐的动力学。环境。科学。Technol公司。30(5), 1614-1617 (1996) ·doi:10.1021/es950618m [52] Chilakapati，A.等人：地下水流动、多组分运移和生物地球化学：耦合过程模型的开发和应用。J.康塔姆。水文学。43(3-4), 303-325 (2000) ·doi:10.1016/S0169-7722（99）00107-2 [53] Knight，P.，Ruiz，D.，Ucar，B.：矩阵缩放的对称保持算法。SIAM J.矩阵分析。申请。35(3), 931-955 (2014) ·Zbl 1329.65089号 ·数字对象标识代码：10.1137/10825753 [54] Golub，H.V.，Van Loan，C.F.：矩阵计算。第三版，约翰·霍普金斯大学出版社，巴尔的摩（1996）·Zbl 0865.65009号 [55] Davis，T.A.，Duff，I.S.：非对称稀疏矩阵的组合单面/多面方法。ACM事务处理。数学。柔和。25(1), 1-20 (1999) ·Zbl 0962.65027号 ·doi:10.1145/305658.287640 [56] Woźnicki，Z.：关于求解非对称线性系统的SOR方法的性能。J.计算。申请。数学。137(1), 145-176 (2001) ·Zbl 1002.65039号 ·doi:10.1016/S0377-0427（00）00705-6 [57] Saad，Y.，Van Der Vorst，H.A.：20世纪线性系统的迭代解。J.计算。申请。数学。123(1-2), 1-33 (2000) ·Zbl 0965.65051号 ·doi:10.1016/S0377-0427（00）00412-X [58] Diersch，H.J.G.:FEFLOW参考手册。DHI-WASY GmbH，柏林（2009） [59] Van der Lee，J.等人：反应传输代码HYTEC的介绍和应用。摘自：Hassanizadeh，S.M.（编辑）《水科学发展，水资源计算方法》，《第十四届水资源计算法国际会议论文集》（CMWR XIV），第599-606页。爱思唯尔（2002） [60] Press，W.H.，S.A.T.，Vettering，W.T.，Flannery，B.P.：《FORTRAN中的数值配方：科学计算的艺术》，第二版。，第123-124页。剑桥大学出版社，纽约（1992）·Zbl 0778.65002号 [61] 线性代数包（LAPACK）可从以下地址免费获得：http://www.netlib.org/lapack网站 [62] Kincaid，D.，Cheney，W.：《数值分析：科学计算的数学》，第三版。美国数学学会（2002）·Zbl 0877.65002号 [63] HSL：用于大规模科学计算的Fortran代码集合。http://www.hsl.rl.ac.uk (2013) [64] 第八章非线性方程组。在：计算数学研究，Claude，B.编辑。1997年，爱思唯尔。第287-336页·Zbl 0889.65023号 [65] Soleymani，F.：计算矩阵反演的快速数值算法。国际数学杂志。数学。科学。2012 (2012) ·Zbl 1253.65051号 [66] Soleymani，F.：关于矩阵近似逆的快速迭代方法。Commun公司。韩国数学。Soc.28（2），407-418（2013）·Zbl 1276.65020号 ·doi:10.4134/CKMS.2013.28.2.407 [67] Morin，K.A.：计算水中化学平衡的计算机化方法的简化解释和示例。计算。地质科学。11, 409-416 (1985) ·doi:10.1016/0098-3004（85）90024-X 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。