×

一种使用自动微分工具开发离散伴随解算器的方法。 (英语) Zbl 1184.76832号

摘要:提出了一种快速开发计算流体动力学(CFD)模型伴随求解器的方法。该方法依赖于使用自动微分(AD)工具来几乎完全自动化离散伴随解算器的开发过程。该方法用于为两种不同的3D CFD解算器生成伴随代码:以单元为中心的单块单处理器模式的Euler解算器和以顶点为中心的多块磁流体力学(MHD)解算器。我们没有使用AD区分CFD解算器的整个源代码,而是选择性地将其应用于生成代码,以计算通量雅可比矩阵的转置以及使用伴随方法计算灵敏度所需的其他偏导数。然后使用科学计算可移植可扩展工具包(PETSc)库求解离散伴随方程。选择性应用AD是这种新方法的主要思想,我们称之为AD伴随(adjoint)。ADjoint方法的优点是它适用于任何一组控制方程和目标函数,并且它与通过原始离散解算器的精确数值微分计算的梯度完全一致。此外,该方法不需要手动微分,从而避免了开发偏微分方程的离散伴随求解器通常需要的较长开发时间,以及在复杂守恒定律系统的微分过程中使用的必要近似所导致的错误。这些优点是以增加离散伴随解算器的内存需求为代价的。然而,考虑到并行计算机中通常可用的内存量以及多核处理器数量增加的趋势,与所展示的非常显著的优势相比,这一劣势相当小。使用离散伴随求解器获得的阻力系数和升力系数对不同参数的敏感性与复步法和有限差分法产生的基准结果非常一致。此外,对于所使用的两个CFD解算器,该方法的总体性能优于大多数传统伴随方法。

MSC公司:

76N25号 可压缩流体和气体动力学的流量控制与优化
76周05 磁流体力学和电流体力学
76M99型 流体力学基本方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Balay S.,《科学计算中的现代软件工具学报》第163页–(1997)·doi:10.1007/978-1-4612-1986-68
[2] Balay,S.、Buschelman,K.、Eijkhout,V.、Gropp,W.D.、Kaushik,D.、Knepley,M.G.、McInnes,L.C.、Smith,B.F.和Zhang,H.,PETSc用户手册,技术报告ANL-95/11–2.3.0版,阿贡国家实验室。2004
[3] Carle,A.和Fagan,M.,ADIFOR 3.0概述,技术报告CAAM-TR-00-02,莱斯大学2000
[4] Carle A.,AIAA 1998–4807号文件,1998–(1998)
[5] 内政部:10.1006/jcph.1994.1057·Zbl 0832.65098号 ·doi:10.1006/jcph.1994.1057
[6] 内政部:10.1006/jcph.1998.6114·Zbl 0921.65059号 ·doi:10.1006/jcph.1998.6114
[7] 内政部:10.1002/fld.885·Zbl 1134.76431号 ·doi:10.1002/fld.885
[8] Driver J.,第44届AIAA航空航天科学会议记录和展览(2006年)
[9] Dwight R.P.,第44届AIAA航空航天科学会议记录和展览(2006年)
[10] DOI:10.1016/j.未来.2004.11.005·doi:10.1016/j.future.2004.11.005
[11] Faure,C.和Papegay,Y.,《奥德赛1.6版:语言参考手册》,INRIA Rapport Technique 211
[12] 第四S.A.,《自动区分:从模拟到优化》,第153页–(2002年)
[13] Gaitonde D.V.,第43届AIAA航空航天科学会议记录和展览(2005年)
[14] Giering R.,GAMM 2002年会议记录(2002年)
[15] DOI:10.1016/j.未来.2004.11.003·doi:10.1016/j.未来.2004.11.003
[16] DOI:10.1023/A:1011430410075·Zbl 0996.76023号 ·doi:10.1023/A:1011430410075
[17] Gockenbach M.S.,理解TAMC生成的代码(2000)
[18] Griewank A.,《衍生品评估》(2000年)
[19] Hascoöt,L.和Pascual,V.,TAPENADE 2.1用户指南,技术报告300,INRIA 2004
[20] Hascoöt L.,《自动区分:应用、理论和工具》(2005)
[21] DOI:10.1016/j.未来.2004.11.010·doi:10.1016/j.future.2004.11.010
[22] Horwedel J.E.,《算法的自动区分:理论、实现和应用》,第243页–(1991年)
[23] DOI:10.1007/BF01061285·Zbl 0676.76055号 ·doi:10.1007/BF01061285
[24] Lee K.H.,第44届AIAA航空航天科学会议记录和展览(2006年)
[25] Marta A.C.,第三届AIAA流量控制会议记录(2006年)
[26] Marta A.C.,第14届AIAA/AHI国际空间飞机和高超音速系统与技术会议记录(2006年)
[27] Martins J.R.R.A.,第39届航空航天科学会议论文集(2001年)
[28] 内政部:10.1145/838250.838251·Zbl 1072.65027号 ·doi:10.1145/838250.838251
[29] 内政部:10.2514/1.11478·数字对象标识代码:10.2514/1.11478
[30] DOI:10.1023/B:OPTE.0000048536.47956.62·兹比尔1145.76418 ·doi:10.1023/B:OPTE.000048536.47956.62
[31] Martins J.R.R.A.,第二届AIAA多学科设计优化专家会议记录(2006年)
[32] Martins J.R.R.A.,第十一届AIAA/ISSMO多学科分析与优化会议记录(2006年)
[33] DOI:10.1016/j.jcp.2004.03.001·Zbl 1071.65025号 ·doi:10.1016/j.jcp.2004.03.001
[34] DOI:10.1023/B:JOMP.000027955.75872.3f·Zbl 1085.76050号 ·doi:10.1023/B:JOMP.000027955.75872.3f
[35] Mohammadi B.,《计算微分:技术、应用和工具》,第309页–(1996)
[36] 国际超级计算机应用杂志。高性能计算。第8页159页–(1994)
[37] 内政部:10.1145/1114268.1114270·Zbl 1136.65310号 ·数字对象标识代码:10.1145/1114268.1114270
[38] 内政部:10.2514/1.10415·数字对象标识代码:10.2514/1.10415
[39] DOI:10.1017/S002211207400223·Zbl 0281.76020号 ·doi:10.1017/S002211207400223
[40] 内政部:10.2514/2.2413·数字对象标识代码:10.2514/2.2413
[41] 内政部:10.2514/2.2414·数字对象标识代码:10.2514/2.2414
[42] 内政部:10.1137/0907058·Zbl 0599.65018号 ·doi:10.1137/0907058
[43] 内政部:10.1006/jcph.1996.0252·Zbl 0933.76070号 ·doi:10.1006/jcph.1996.0252
[44] Van der Weide E.,第44届AIAA航空航天科学会议记录和展览(2006年)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。