×
Jean-François Lemieux;普莱斯（Stephen F.Price）。;凯瑟琳·埃文斯。;丹娜·诺尔;安德鲁·塞林格。;霍兰德（David M.Holland）。;安东尼·佩恩。

实现无雅可比牛顿-克利洛夫方法求解一阶冰盖动量平衡。 （英语） Zbl 1229.86013号

J.计算。物理学。 230，第17号，6531-6545（2011）.
小结：为了提高Glimmer-Community冰盖模型（Glimmer-CISM）的数值性能，该模型是社区-地球系统模型（CESM）的陆冰组成部分，我们采用了无雅可比矩阵的Newton-Krylov（JFNK）方法求解一阶冰盖动量方程。我们的JFNK实现基于对现有代码的大量重用。例如，我们的基于物理的预处理程序使用Glimmer-CISM中的原始Picard线性解算器。对于跨越一系列几何体和边界条件的几个测试用例，我们的JFNK实现比Glimmer-CISM中的标准Picard解算器效率高1.8–3.6倍。重要的是，在细化网格时，JFNK相对于Picard解算器的计算增益会增加。通过重新缩放基本边界条件的方程和使用不精确的牛顿方法，JFN解算器全局收敛性得到了显著改善。虽然一组不同的测试用例表明我们的JFNK实现通常是健壮的，但对于某些问题，它可能无法随着分辨率的增加而收敛（就像Picard求解器一样）。通过参数连续化的全球化并没有解决这个问题，未来提高鲁棒性的工作将探索Picard和JFNK的结合以及同伦方法的使用。

引用于1文件

MSC公司：

86A40型 冰川学
86-08 地球物理问题的计算方法

关键词：

冰盖模型;冰流变学;纽顿-克利洛夫;GMRES公司;ILU（仪表着陆装置）

软件：

CESM公司;SMIP-HOM公司;特里利诺斯;JFNK公司
PDF格式BibTeX公司 XML格式引用
\textit{J.-F.Lemieux}等人，《计算杂志》。物理学。230，第17号，6531--6545（2011；Zbl 1229.86013）
全文： 内政部

参考文献：

[1] 乔恩，I。；阿卜杜拉蒂，W。；Fahnestock，M.，《格陵兰雅各布港-伊斯布雷冰川速度的大波动》，《自然》，432608-610（2004）
[2] Howat，I.M。；乔恩，I。；Tulaczyk，S。；Gogineni，S.，《东格陵兰岛海尔海姆冰川的快速后退和加速》，Geophys。Res.Lett.公司。，32，L22502（2005）
[3] Rignot，E。；Kanagratnam，P.，《格陵兰冰盖速度结构的变化》，《科学》，311986-990（2006）
[4] Howat，I.M。；乔恩，I。；Scambos，T.A.，《格陵兰出口冰川冰流量的快速变化》，《科学》，315581559-1561（2007）
[5] A.勒克曼。；Murray，T。；德兰格，R。；Hanna，E.，《东格陵兰岛的快速同步冰动力变化》，Geophys。Res.Lett.公司。，第33页，第03503页（2006年）
[6] Murray，T。；Scharrer，K。；詹姆斯·T·D。；染料，S.R。；汉娜，E。；A.D.布斯。；塞尔姆斯，N。；A.勒克曼。；Hugues，A.L.C。；库克，S。；Huybrechts，P.，《格陵兰东南部冰川动力学的海洋调节假说及其对冰盖质量变化的影响》，J.Geophys。第115号决议，F03026（2010年）
[7] De Angelis，H。；Skvarca，P.，《冰架崩塌后的冰川涌浪》，《科学》，29956121560-1562（2003）
[8] 斯坎波斯，T.A。；Bohlander，J.A。；舒曼，C.A。；Skvarca，P.，《南极Larsen B海湾冰架崩塌后冰川加速和变薄》，Geophys。Res.Lett.公司。，31，L18402（2004）
[9] Rignot，E。；Bamber，J.L。；Van den Broeke，医学博士。；Davis，C。；李毅。；范德伯格，W.J。；Van Meijgaard，E.，《雷达干涉测量和区域气候建模导致的近期南极冰块损失》，《自然地球科学》，1，2，106-110（2008）
[10] Wingham，D.J。；Wallis，D.W。；Shepherd，A.，松岛冰川变薄的时空演变，1995-2006年，地球物理学。Res.Lett.公司。，36、17、L17501（2009）
[11] Holland，D.M。；托马斯·R·H。；de Young，B。；里贝加德，M.H。；Lyberth，B.，温暖的地下海水引发的雅各布港-伊斯布雷加速度，《自然地质学》。，1, 659-664 (2008)
[12] 托马，M。；Jenkins，A。；Holland，D。；Jacobs，S.，《南极阿蒙森海大陆架环极深水侵入体建模》，Geophys。Res.Lett.公司。，18602年1月35日（2008年）
[13] Shepherd，A。；哈伯德，A。；尼诺，P。；金·M。；麦克米兰，M。；Joughin，I.，格陵兰冰盖运动与夏末的每日融化相结合，地球物理。Res.Lett.公司。，36、1、L01501（2009）
[14] 范登·布鲁克，M。；Bamber，J。；Ettema，J。；Rignot，E。；Schrama，E。；范德伯格，W.J。；van Meijgaard，E。；Velicogna，I。；Wouters，B.，《划分格陵兰最近的质量损失》，《科学》，326595984-986（2009）
[15] Schoof，C.，《冰盖接地线动力学：稳态、稳定性和滞后》，J.Geophys。研究，112，F03S28（2007）
[16] Katz，R.F。；Worster，M.G.，《冰层接地线的稳定性》，Proc。R.Soc.A，4661597-1620（2010）·Zbl 1187.86004号
[17] Goldberg，D。；荷兰医学博士。；Schoof，C.，海洋冰盖中的接地线运动和冰架支撑，J.Geophys。Res.，114，F4，F04026（2009）
[18] 小巷，R.B。；Clark，P.U。；Huybrechts，P。；Joughin，I.，《冰盖和海平面变化》，《科学》，310，5747，456-460（2005）
[19] Bamber，J.L。；小巷，R.B。；Joughin，I.，《现代冰盖对外部强迫的快速响应》，《地球》。科学。莱特。，257, 1-2, 1-13 (2007)
[20] Shepherd，A。；Wingham，D.，《南极和格陵兰冰盖的近期海平面贡献》，《科学》，31558181529-1532（2007）
[21] R.Alley，《2007年气候变化：物理科学基础》。第一工作组对政府间气候变化专门委员会的贡献，第1章：决策者摘要。，剑桥大学出版社，2007年。；R.Alley，《2007年气候变化：物理科学基础》。第一工作组对政府间气候变化专门委员会的贡献，第1章：政策制定者摘要。，剑桥大学出版社，2007年。
[22] 鲁特，I.C。；M.哈格顿。；新泽西州赫尔顿。；Payne，A.J.，Glimmer群落冰盖模型，地球物理学杂志。第114号决议，F02004（2009）
[23] Paterson，W.S.B.，《冰川物理学》（1994），佩加蒙：新泽西州塔里敦佩加蒙牛津。
[24] 杜兰德，G。；加利亚迪尼，O。；de Fleurian，B。；Zwinger，T。；Le Meur，E.，《海洋冰盖动力学：滞后和中性平衡》，J.Geophys。决议，114，F03009（2009）
[25] Pattyn，F.，《一个新的三维高阶热力冰盖模型：基本敏感性、冰流发展和冰下湖泊的冰流》，J.Geophys。决议，108，B8，2382（2003）
[26] 斯科夫，C。；Hindmarsh，R.C.A.，《薄膜流动与壁面滑移：高阶冰川流动模型的渐近分析》，Q.J.Mech。申请。数学。，63, 73-114 (2010) ·Zbl 1395.76011号
[27] 杜科维茨，J.K。；价格，S.F。；Lipscomb，W.H.，《基于最小作用原理的冰层动力学一致近似和边界条件》，J.Glaciol。，56, 197, 480-496 (2010)
[28] T.Payne，S.Price，《三维一阶冰流模型：数值实施和基准测试》，编制中。；T.Payne，S.Price，三维一阶冰流模型：数值实施和基准测试，正在准备中。
[29] Pattyn，F。；Perichon，L。；Aschwanden，A。；布鲁尔，B。；de Smedt，B。；加利亚迪尼，O。；古德蒙德松，G.H。；Hindmarsh，R。；哈伯德，A。；约翰逊，J.V。；Kleiner，T。；科诺瓦洛夫，Y。；马丁·C。；Payne，A.J。；波拉德，D。；价格，S。；吕坎普，M。；斋藤，F。；苏切克，O。；杉山，S。；Zwinger，T.，高阶和全斯托克斯冰盖模型的基准实验（ISMIP-HOM），《冰层讨论》。，2, 111-151 (2008)
[30] Hindmarsh，R。；Payne，A.，冰盖方程稳定解的时间步长极限，Ann.Glaciol。，23, 74-85 (1996)
[31] Lemieux，J.-F。；Tremblay，B.，粘塑性海冰模型的数值收敛，J.Geophys。研究，114，C05009（2009）
[32] 帕尼科尼，C。；Putt，M.，多维变饱和水流问题数值解中Picard和Newton迭代法的比较，水资源。研究，30，12，3357-3374（1994）
[33] Knoll，D.A。；Keyes，D.E.，《无雅可比牛顿-克利洛夫方法：方法和应用调查》，J.Compute。物理。，193, 2, 357-397 (2004) ·Zbl 1036.65045号
[34] Reisner，J.M。；Mousseau，V.A。；Wyszogrodzki，A.A。；Knoll，D.A.，《基于物理预处理的隐式平衡飓风模型》，Mon。维阿。修订版，1331003-1022（2005）
[35] 埃文斯，K.J。；Knoll，D.A。；Pernice，M.A.，高效模拟对流和相变的二维算法开发，J.Compute。物理。，219404-417（2006），10.1016/j.jcp.2006.03.025·Zbl 1102.76037号
[36] Lemieux，J.-F。；特伦布雷，B。；塞德拉切克，J。；塔珀，P。；托马斯·S。；Huard，D。；Auclair，J.-P.，用无雅可比牛顿-克利洛夫方法改进粘塑性海冰模型的数值收敛性，J.Compute。物理。，2292840-2852（2010），10.1016/j.jcp.2009.12.01·Zbl 1184.86004号
[37] M.Seager，《大众SLAP》，技术代表UCRL-100267，加州利弗莫尔劳伦斯·利弗莫雷国家实验室（1988）。；M.Seager，《面向大众的SLAP》，技术代表UCRL-100267，加州利弗莫尔劳伦斯·利弗莫雷国家实验室（1988年）。
[38] 萨阿德，Y。；Schultz，M.H.，GMRES:求解非对称线性系统的广义最小残差算法，SIAM J.Sci。统计计算。，7, 3, 856-869 (1986) ·Zbl 0599.65018号
[39] Axelsson，O.，《迭代求解方法》（1994），剑桥大学出版社：剑桥大学出版社，纽约·Zbl 0795.65014号
[40] Y.Saad，稀疏线性系统的迭代方法，PWS，1996。；Y.Saad，稀疏线性系统的迭代方法，PWS，1996年·Zbl 1031.65047号
[41] Saad，Y.，《一种灵活的内外预处理GMRES算法》，SIAM J.Sci。计算。，14, 2, 461-469 (1993) ·Zbl 0780.65022号
[42] 杜米纳罗，R.S。；Walker，H.F。；Shadid，J.N.，关于Newton-GMRES方法中的回溯失败，以及Navier-Stokes方程的演示，J.Compute。物理。，180, 549-558 (2002) ·Zbl 1143.76489号
[43] 艾森斯塔特，南卡罗来纳州。；Walker，H.F.，《用不精确牛顿法选择强迫项》，SIAM J.Sci。计算。，17, 16-32 (1996) ·Zbl 0845.65021号
[44] Bamber，J.L。；Layberry，R.L。；Gogineni，S.P.，《格陵兰冰盖1的新冰层厚度和冰层数据集：测量、数据简化和误差》，J.Geophys。决议，106，D24，33773-33780（2001）
[45] 杜科维茨，J.K。；Baumgardner，J.R.，作为传输/平流算法的增量重映射，J.Compute。物理。，160、318-335（2000）、10.1006/jcph.2000.6465·Zbl 0972.76079号
[46] Lipscomb，W.H。；Hunke，E.C.，《使用增量重映射模拟海冰运输》，周一。维阿。版本：1321341-1354（2004）
[47] Knoll，D.A。；Rider，W.J.，《多重网格预处理Newton-Krylov方法》，SIAM J.Sci。计算。，21, 691-710 (1999) ·Zbl 0952.65102号
[48] 博塔，E.F.F。；Wubs，F.W.，《矩阵重新编号ILU：稀疏矩阵的有效代数多级ILU预条件器》，SIAM J.Matrix Anal。申请。，20, 4, 1007-1026 (1999) ·Zbl 0937.65057号
[49] Heroux，医学硕士。；Bartlett，R.A。；豪，V.E。；霍克斯特拉，R.J。；胡建杰。；科尔达·T·G。；Lehoucq，R.B。；Long，K.R。；Pawlowski，R.P。；Phipps，E.T。；塞林格，A.G。；香港桑奎斯特。；杜米纳罗，R.S。；Willenbring，J.M。；Williams，A.，Trilinos项目概述，ACM Trans。数学。柔软。，31, 3, 397-423 (2005) ·Zbl 1136.65354号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。