彼得·卢卡斯；范祖伊伦（Alexander H。；比尔，赫斯特 使用无雅可比牛顿-克利洛夫算法进行快速非定常流动计算。 （英语） 兹比尔1427.76214 J.计算。物理学。 229，第24号，9201-9215（2010）. 摘要：尽管过去几十年来计算机能力和数值算法取得了进步，但非定常流动问题的解决方案仍然需要大量的计算时间。特别是对于大雷诺数流动，常用于求解非线性方程组的非线性多重网格收敛速度较慢。这种求解方法不能很好地处理由大纵横比单元和湍流引起的刚度。在之前的工作中，我们表明，对于二维大雷诺数，与标准非线性多重网格相比，无雅可比矩阵的Newton-Krylov（jfnk）算法，通过近似匹配目标残差算子的雅可比函数的近似因式分解预处理，可以使速度提高10倍，非恒定流计算。本文的目的是证明jfnk算法也适用于处理由最大纵横比、网格密度、物理时间步长和雷诺数引起的刚度。与标准的非线性多重网格相比，可以实现高达25倍的加速。 引用于4文件 MSC公司： 76米99 流体力学基本方法 关键词：无Jacobian牛顿-克利洛夫；非恒定流；不同流量和网格参数引起的刚度；高阶隐式时间积分 软件：ITSOL公司；FLUENT公司；伊卢包；斯帕拉尔-奥尔马拉斯；FUN3D公司 PDF格式BibTeX公司 XML格式引用 \textit{P.Lucas}等人，J.Compute。物理学。229，第24号，9201--9215（2010；Zbl 1427.76214） 全文： 内政部 参考文献： [1] R.P.L.Nijssen，风力涡轮机转子叶片复合材料的疲劳寿命预测和强度退化，博士论文，荷兰代尔夫特理工大学，2006年11月，<darenet.nl；R.P.L.Nijssen，风力涡轮机转子叶片复合材料的疲劳寿命预测和强度退化，博士论文，荷兰代尔夫特理工大学，2006年11月，<darenet.nl [2] H.F.Veldkamp，《风能的机会：风力涡轮机疲劳设计的概率方法》，荷兰代尔夫特理工大学博士论文，2006年10月，darenet.nl；H.F.Veldkamp，《风能的机会：风力涡轮机疲劳设计的概率方法》，荷兰代尔夫特理工大学博士论文，2006年10月，darenet.nl [3] Fluent Inc.，Fluent 6.1用户指南，技术代表，黎巴嫩，NH 037662003年2月。；Fluent Inc.，Fluent 6.1用户指南，技术代表，黎巴嫩，NH 037662003年2月。 [4] 乐趣乐趣<http://fun3d.larc.nasa.gov>;乐趣乐趣<http://fun3d.larc.nasa.gov> [5] Numeca International，用户手册FINE™/Hexa v2.5（包括Hexstream）-Flow Integrated Environment，技术代表，布鲁塞尔，2007年12月<http://www.numeca.be网站>; Numeca International，用户手册FINE™/Hexa v2.5（包括Hexstream）-Flow Integrated Environment，技术代表，布鲁塞尔，2007年12月<http://www.numeca.be网站> [6] 卢卡斯，P。；Bijl，H。；van Zuijlen，A.H.，非结构化网格上的高效非定常高雷诺数流动计算，计算机与流体，39，271-282（2010）·兹比尔1242.76175 [7] J.V.Lassaline，D.W.Zingg，《应用于凝聚多重网格的方向粗化和线隐式平滑的研究》，载于：第16届AIAA计算流体动力学会议，美国佛罗里达州奥兰多，2003年，AIAA论文2003-3435。；J.V.Lassalize，D.W.Zingg，《定向粗化和线隐平滑应用于凝聚多重网格的研究》，载于：第16届AIAA计算流体动力学会议，美国佛罗里达州奥兰多，2003年，AIAA论文2003-3435。 [8] D.J.Mavrilis，《高雷诺数粘性流的定向凝聚多重网格技术》，技术代表98-7，ICASE，弗吉尼亚州汉普顿，NASA/CR-1998-206911，1998年1月。；D.J.Mavrilis，《高雷诺数粘性流定向凝聚多重网格技术》，技术代表98-7，ICASE，弗吉尼亚州汉普顿，NASA/CR-1998-206911，1998年1月。 [9] D.J.Mavrilis，《各向异性非结构网格上粘性流体求解器的多重网格策略》，技术代表98-6，ICASE，弗吉尼亚州汉普顿，NASA/CR-1998-2069101998年1月。；D.J.Mavrilis，《各向异性非结构网格上粘性流求解器的多重网格策略》，技术代表98-6，ICASE，弗吉尼亚州汉普顿，NASA/CR-1998-2069101998年1月·兹伯利0926.76066 [10] Bijl，H。；Carpenter，M.H.，使用高阶时间积分方案进行非定常流动计算的迭代求解技术，《国际流体数值方法杂志》，47，857-862（2005）·Zbl 1134.76418号 [11] S.Isono，D.W.Zingg，《应用于非定常流的四阶隐式时间推进的Runge Kutta Newton Krylov算法》，载于：第42届美国航空航天协会航空航天与科学会议和展览，美国内华达州雷诺市，2004年，美国航空航天协会论文2004-433。；S.Isono，D.W.Zingg，应用于非定常流的四阶隐式时间推进的Runge-Kutta-Newton-Krylov算法，载于：第42届美国航空航天与科学会议和展览，美国内华达州里诺，2004年，美国航空航天局论文2004-433。 [12] 乔蒂普拉萨德，G。；Mavrilis，D.J。；Caughey，D.A.，非结构网格上非定常Navier-Stokes方程的高阶时间积分格式，计算物理杂志，191，542-566（2003）·Zbl 1134.76428号 [13] 秦，N。；Ludlow，D.K。；Shaw，S.T.，无矩阵预条件牛顿/gmres公司非定常Navier-Stokes解的方法，国际流体数值方法杂志，33，223-248（2000）·Zbl 0976.76049号 [14] Rumpfkeil，M.P。；Zingg，D.W.，离散伴随方法对非定常流的最优控制，优化与工程，11，1，5-22（2010）·Zbl 1273.76097号 [15] M.Tabesh，D.W.Zingg，使用Newton-Krylov算法的高效隐式时间推进方法，载于：第47届AIAA航空航天科学会议，包括新视野论坛和航空航天博览会，美国佛罗里达州奥兰多，2009年，AIAA论文2009-164。；M.Tabesh，D.W.Zingg，使用Newton-Krylov算法的高效隐式时间推进方法，载于：第47届AIAA航空航天科学会议，包括新视野论坛和航空航天博览会，美国佛罗里达州奥兰多，2009年，AIAA论文2009-164。 [16] 卢卡斯，P。；A.H.van，Zuijlen。；Bijl，H.，A预处理jfnk公司应用于非定常不可压缩流动和流固耦合问题的算法，CMES：工程与科学中的计算机建模，59，1，79-106（2010）·Zbl 1231.74111号 [17] Chisholm，T.T。；Zingg，D.W.，《可压缩湍流的无雅可比牛顿-克利洛夫算法》，计算物理杂志，2283490-3507（2009）·Zbl 1319.76033号 [18] P.R.Spalart，S.R.Allmaras，《气动流动的单方程湍流模型》，载于：美国航空航天局第30届航空航天科学会议与展览，美国内华达州雷诺，1992年，美国航空航天协会论文92-0439。；P.R.Spalart，S.R.Allmaras，《气动流动的单方程湍流模型》，载于：AIAA第30届航空航天科学会议与展览，美国内华达州雷诺，1992年，AIAA论文92-0439。 [19] A.Jameson，W.Schmidt，E.Turkel，用有限体积法和Runge-Kutta时间步进格式对Euler方程进行数值求解，载于：第14届AIAA流体和等离子体动力学会议，美国加利福尼亚州帕洛阿尔托，1981年，第23-25页，AIAA论文81-1259。；A.Jameson，W.Schmidt，E.Turkel，用有限体积法和Runge-Kutta时间步进格式对Euler方程进行数值求解，载于：第14届AIAA流体和等离子体动力学会议，美国加利福尼亚州帕洛阿尔托，1981年，第23-25页，AIAA论文81-1259。 [20] Bijl，H。；Carpenter，M.H。；Vatsa，V.N。；Kennedy，C.A.，非定常可压缩Navier-Stokes方程的隐式时间积分格式：层流，计算物理杂志，179，313-329（2002）·Zbl 1060.76079号 [21] Carpenter，M.H。；肯尼迪，C.A。；Bijl，H。；Viken，S.A。；Vatsa，V.N.，流体力学应用的四阶Runge-Kutta格式，科学计算杂志，25157-194（2005）·Zbl 1203.76112号 [22] 肯尼迪，C。；Carpenter，M.H.，对流-扩散-反应方程的加法Runge-Kutta格式，应用数值数学，44139-181（2003）·Zbl 1013.65103号 [23] A.Patel，Dédevelopment D'un solveur公司rans公司; A.Patel，Dédevelopment D'un solveur公司rans公司 [24] 德怀特，R。；卢卡斯，P。；Bijl，H.，（Blockley，R.；Shyy，W.，《CFD中的求解方法和加速技术》，受邀为《航空航天工程百科全书》（2010）撰稿，John Wiley and Sons Ltd.） [25] 周，E。；Saad，Y.，不定矩阵ILU预条件的实验研究，计算与应用数学杂志，86，387-414（1997）·Zbl 0891.65028号 [26] Knoll，D.A。；Keyes，D.E.，无Jacobian牛顿-克利洛夫方法：方法和应用调查，计算物理杂志，193，357-397（2004）·Zbl 1036.65045号 [27] 萨阿德，Y。；M.H.舒尔茨。，gmres公司：求解非对称线性系统的广义最小残差算法，SIAM科学与统计计算杂志，7，3，856-869（1986）·Zbl 0599.65018号 [28] Pueyo，A。；Zingg，D.W.，空气动力学计算的高效Newton-Krylov求解器，美国航空航天学会期刊，36,11191991-1997（1998） [29] P.Lucas，《计算非结构网格上非定常亚音速流的有效数值方法》，荷兰代尔夫特理工大学博士论文，2010年6月<http://repository.tudelft.nl>; P.Lucas，《计算非结构网格上非定常亚音速流的有效数值方法》，荷兰代尔夫特理工大学博士论文，2010年6月<http://repository.tudelft.nl> [30] Y.Saad等人，Itsol-一般稀疏线性方程组的迭代求解器库，2006年<http://www-users.cs.umn.edu/∼;saad/software/ITSOL/>；Y.Saad等人，Itsol-一般稀疏线性方程组的迭代求解器库，2006<http://www-users.cs.umn.edu/∼;saad/software/ITSOL/> [31] M.Bollhöfer，Y.Saad，O.Schenk，Ilupack-预处理软件包，2006年<http://www.math.tu-berlin.de/ilupack网站>; M.Bollhöfer，Y.Saad，O.Schenk，Ilupack-预处理软件包，2006年<http://www.math.tu-berlin.de/ilupack网站> [32] Wong，P。；Zingg，D.W.，使用Newton-Krylov方法在非结构化网格上进行三维空气动力学计算，计算机与流体，37，107-120（2008）·Zbl 1194.76167号 [33] W.A.Timmer，风力涡轮机应用翼型上高迎角气流的一些方面，载于：EWEC 2001年欧洲风能会议，2001年。；W.A.Timmer，风力涡轮机应用翼型上高迎角气流的一些方面，载于：EWEC 2001年欧洲风能会议，2001年。 [34] 卢卡斯，P。；A.H.van，Zuijlen。；Bijl，H.，一种基于解决方案的网格自适应自动方法，用于提高给定数量网格单元的数值精度——应用于六面体网格上的稳定流动，CMES：工程与科学中的计算机建模，41，2，147-175（2009） 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。