×

用于引力动力学多尺度模拟的非侵入式分层耦合策略。 (英语) Zbl 1454.85005号

摘要:分层代码耦合策略使得将单个数值解算器的结果合并为一个自洽辛解成为可能。我们探讨了允许这种耦合策略非侵入的可能性。在这种情况下,底层的数值实现不受耦合本身的影响,但其功能在接口中进行。这种方法对于求解大范围自引力系统的运动方程是有效的。我们采用一个专门的集成器来解决问题的每个特定部分,并将结果合并为一个自我协调的解决方案。特别是,我们探索了将嵌入宏观系统中的一个或多个微观系统的演化组合在一起的可能性。这里给出的概括包括高阶耦合策略(从经典的2阶到10阶),但我们还演示了如何嵌套多个桥,以及如何在桥的时间步长引入额外的过程来丰富物理,例如通过合并耗散过程或。这种增强允许在经典牛顿体积分器中包含额外的过程,而不会更改底层代码。这些附加过程包括例如雅科夫斯基效应、动力摩擦或相对论动力学。其中一些过程适用于所有粒子,而其他过程仅适用于子集。
所提出的方法是非侵入性的,即底层方法在不更改代码的情况下保持可操作性(除了添加get和set-functions以启用桥操作符)。因此,基本积分器继续使用其内部时间步长进行操作,并保留其局部优化和并行性。多个桥可以分层嵌套和耦合,从而可以构建多个嵌套增强桥的复杂环境。虽然耦合拓扑可能变得相当复杂,但我们引入了层次耦合语言(HCL),这是一种可以描述复杂桥拓扑的元语言。元语言旨在激发对更复杂层次结构的讨论,其中桥接操作符作为模式引入。我们给出了其中几种情况的示例应用程序,并讨论了这些积分器可以应用的条件。当我们将该方法应用于模拟银河系(100 kpc-空间尺度和10 Gyr-时间尺度)中的恒星群中的行星系统(非盟空间尺度和年时尺度)时,典型应用在时间和空间尺度上的范围超过10个数量级。

MSC公司:

85甲15 星系和恒星结构
85-10 天文学和天体物理学相关问题的数学建模或模拟
2015年1月70日 天体力学
85-08 天文学和天体物理学相关问题的计算方法
65页第10页 含辛积分器哈密顿系统的数值方法
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 图灵,A.M.,《关于可计算数字及其在entscheidungsproblem中的应用》,《伦敦数学学会学报》,第2期,第42期,第230-265页(1936年)·JFM 62.1059.03标准
[2] Oberkampf,W.L。;Roy,C.J.,《科学计算中的验证与确认》,第1版(2010年),剑桥大学出版社:美国纽约州剑桥大学出版社·Zbl 1211.68499号
[3] Moore,G.E.,《将更多元件注入集成电路》,《电子学》,第38期,第8-12页(1965年)
[4] Felleisen,M.,《论编程语言的表达能力》,《科学计算程序》,第17、1、35-75页(1991年)·Zbl 0745.68033号
[5] Hoekstra,A.G。;Zwart,S.P。;Coveney,P.V.,《多尺度建模、仿真和计算:从桌面到exascale》,Philos Trans R Soc A:数学物理工程科学,377,2142,20180355(2019)
[6] Manser,C.J.,在白矮星周围的碎片盘中运行的星子,科学,3646435,66-69(2019)
[7] Fowler,M.,《重构:改进现有代码的设计》(1999),Addison-Wesley:Addison-Whesley Boston,MA,USA·Zbl 1020.68632号
[8] 565页。(ISBN:9781584889090,ISBN:10584889098),第147-164页
[9] 格伦·D·。;波特吉斯·兹瓦特,S。;石山,T。;Makino,J.,在行星范围的分布式超级计算机上进行高性能引力n体模拟,计算科学与Discov,4,1,015001(2011)
[10] 牛顿,I.,《自然哲学与数学原理》,1(1687)·Zbl 0050.00201号
[11] 瓦尔顿,M。;Karttunen,H.,《三体问题》(2006),剑桥大学出版社·邮编1093.70001
[12] Portegies Zwart,S.F。;Boekholt,T.C.,牛顿运动方程微观时间可逆性的数值验证:对抗指数发散,Commun非线性科学数值模拟,61160-166(2018)·Zbl 1473.70057号
[13] 巴恩斯,J。;Hut,P.,分层o(nlogn)力计算算法,国家自然科学院,324446-449(1986)
[14] Fujii,M。;岩川,M。;Y.Funato。;Makino,J.,BRIDGE:一种直接树混合n体算法,用于星团及其母星系的完全自洽模拟,Publ Astr Soc Jpn,59,1095(2007)
[15] Wisdom,J。;霍尔曼,M.,n体问题的辛映射,AJ,102,1528-1538(1991)
[16] Makino,J。;Aarseth,S.J.,《关于引力多体问题中采用Ahmad-Cohen方案的hermite积分器》,Publ Astr Soc Jpn,44,141-151(1992)
[17] 米科拉,S。;Aarseth,S.J.,《n体模拟中二进制的有效积分方法》,《新天文》,3309-320(1998)
[18] 岩川,M。;波特吉斯·兹瓦特,S。;Makino,J.,用于模拟稠密恒星星团系统的支持GPU的粒子-粒子无粒子方案,《计算天体物理宇宙》,2,6(2015)
[19] 贝多夫,J。;Portegies Zwart,S.,《许多微小合并对致密静止星系尺寸增长的影响》,MNRAS,431767-780(2013)
[20] 第97-104页
[21] Murray C.D.,Dermott S.F.,《太阳系动力学》。1999年,C.D.Murray,S.F.Dermott·兹比尔0957.70002
[22] 第54-65页
[23] 本德,M。;Heenen,P.H。;Reinhard,P.G.,核结构的自洽平均场模型,《物理评论》,75,1,121-180(2003)
[24] 梅隆,Y。;Li,B。;霍利·伯克曼,K。;Spurzem,R.,《无碰撞恒星动力学模拟的膨胀技术》,ApJ,792,98(2014)
[25] Dehnen,W.,《恒星动力学的快速多极方法》,《计算天体物理学-宇宙》,第1期,第1页(2014年)
[26] Milgrom,M.,牛顿动力学的一种修改,作为隐藏质量假说的可能替代方案,ApJ,270365-370(1983)
[27] Will,C.M.,《将后牛顿效应纳入n-体动力学》,《物理学评论D》,89,4,044043(2014)
[28] 医学博士Duez。;Zlochower,Y.,《21世纪紧凑双星的数值相对论》,Rep Prog Phys,82,1,016902(2019)
[29] 斯皮策,L.J。;Hart,M.H.,《随机引力遭遇和球面系统的演化》。ii、。模型,ApJ,166,483(1971)
[30] 苏里亚诺,S.S。;李振英。;克拉斯诺波尔斯基,R。;铃木,T.K。;Shang,H.,风生非理想磁流体盘环和间隙的形成:三维模拟,MNRAS,484,107-124(2019)
[31] 第31卷·Zbl 1094.65125号
[32] Yoshida,H.,高阶辛积分器的构造,Phys-Lett A,150262-268(1990)
[33] Sofroniou,M。;Spaletta,G.,对称积分器对称组成常数的推导,Optim Methods Softw,20,4-5597-613(2005)·Zbl 1082.65070号
[34] Dehnen,W.,一种分层的力计算算法,《计算物理杂志》,179,27-42(2002)·Zbl 1130.85303号
[35] Minchev,I。;Famaey,B.,《星系盘径向迁移的新机制:螺旋-基底共振重叠》,ApJ,722112-121(2010)
[36] Martínez-Barbosa,C.,《太阳诞生星团的演化以及与gaia寻找太阳兄弟姐妹》,MNRAS,4571062-1075(2016)
[37] Pfenniger,D。;Friedli,D.,与三维n体模拟相关的计算问题,aap,270,561-572(1993)
[38] Rubinam,D.P.,小行星的辐射自转和自转,伊卡洛斯,148,2-11(2000)
[39] 沃克鲁利克·D.、波特克·W.F.、切斯利·S.R.、谢瑞斯·D.J.、斯塔特勒·T.S.《雅科夫斯基和约克效应》。2015年,Michel P,DeMeo F.E.,Bottke W.F.,509-531。
[40] 波特吉斯·兹瓦特,S。;McMillan,S.,《天体物理食谱》;AMUSE艺术(2018)
[41] 赫尔斯特罗姆,C。;Mikkola,S.,速度相关扰动的少体问题中的显式算法正则化,天体力学,106143-156(2010)·邮编:1223.70056
[42] 岩川,M。;Oshino,S。;Fujii,M.S。;Hori,Y.,PENTACLE:行星形成的平行粒子-粒子-无粒子代码,Publ Astr Soc Jpn,69,5,81(2017)
[43] Hozumi,S。;岩川,M。;Nitadori,K.,使用基于粒子的方法模拟无碰撞恒星系统合并的平均场方法,ApJ,875,1,20(2019)
[44] 波特吉斯·兹瓦特,S。;麦克米兰,S.L.W。;van Elteren,E。;Pelupessy,I。;de Vries,N.,《使用分层可互换软件接口的多物理模拟》,Comput Phys Commun,183456-468(2013)
[45] Pelupessy,F.I.,《天体物理多用途软件环境》,A&A,557,A84(2013)
[46] Pelupessy,I.,《海洋学多用途软件环境》,Geosci Model Dev Discussion,2016,1-36(2016)
[47] 波特吉斯·兹瓦特,S。;Bédorf,J.,《使用GPU实现天体物理学中计算引力动力学的模拟》,《计算机》,48,11,50-58(2015)
[48] Portegies Zwart,S.,《未来的计算天体物理学》,《科学》,361,6406,979-980(2018)·Zbl 1416.85001号
[49] 普卢默,H.C.,《球状星团中的分布问题》,MNRAS,71,460-470(1911)
[50] 夏尔马,S。;布兰德·霍桑,J。;约翰斯顿,K.V。;Binney,J.,《银河系:生成银河系综合调查的代码》(2011),天体物理学源代码库
[51] Pelupessy,F.I。;Portegies Zwart,S.,嵌入式星团的演化,MNRAS,4201503-1517(2012)
[52] 第267卷,第233页
[53] Boekholt,T。;Portegies Zwart,S.,《关于n体模拟的可靠性》,《计算天体物理学宇宙》,第2期,第2页(2015年)
[54] Verlet,L.,经典流体的计算机“实验”。i.lennard-jones分子的热力学性质,Phys Rev,159,1,98-103(1967)
[55] King,I.R.,《星团的结构》。iii.一些简单的数据模型,AJ,71,64-75(1966)
[56] Barnes,J.E.,《作为平滑操作的重力软化》,MNRAS,425,2,1104-1120(2012)
[57] Portegies Zwart,S.F。;Boekholt,T.C.N.,牛顿运动方程微观时间可逆性的数值验证:对抗指数发散,Commun非线性科学数值模拟,61160-166(2018)·Zbl 1473.70057号
[58] 1987年,191页。
[59] 麦克米兰,S.L.W。;Portegies Zwart,S.F.,银河系中心附近星团的命运。i.分析考虑,ApJ,596,314-322(2003)
[60] 斯宾纳托,P.F。;Fellhauer,M。;Portegies Zwart,S.F.,黑洞螺旋进入星系中心的效率,MNRAS,344,22-32(2003)
[61] 艾伦,C。;Santillan,A.,用于轨道计算的星系质量分布改进模型,Rev Mex Astron Astrofis,22255-263(1991)
[62] 费雷尔,N.M.,《纯粹应用数学》,14,1(1877)
[63] Antoja,T.,用RAVE追踪穿过星系的大力神流对星系棒的约束,A&A,563,A60(2014)
[64] 马丁内斯·巴博萨,哥伦比亚特区。;Brown,A.G.A。;Portegies Zwart,S.,《银河系中太阳的径向迁移:统计研究》,MNRAS,446823-841(2015)
[65] 波特吉斯·兹瓦特,S。;Boekholt,T.,《关于用直接n体方法模拟自重系统所需的最小精度》,ApJL,785,L3(2014)
[66] van Elteren,A。;波特吉斯·兹瓦特,S。;Pelupessy,I。;蔡明霞。;McMillan,S.L.W.,年轻和稠密星团中行星系统的生存能力,A&A,624,A120(2019)
[67] 哈默斯,A.S。;Portegies Zwart,S.F。;Merritt,D.,超大质量黑洞附近恒星的相对论动力学,MNRAS,443,355-387(2014)
[68] Kokubo,E。;Ida,S.,《原行星的寡头增长》,伊卡洛斯,131,171-178(1998)
[69] Pelupessy,F.I。;Jänes,J。;Portegies Zwart,S.,《N-Body积分器与分层分割的单个时间步长》,New Astron,17,711-719(2012)
[70] Portegies Zwart,S.,AMUSE:天体物理多用途软件环境(2011),天体物理源代码库
[71] Portegies Zwart S.等人。AMUSE:天体物理多用途软件环境。2018年b。
[72] Bovy,J.,galpy:银河动力学包(2014),天体物理学源代码库
[73] 第45卷,第303.01页
[74] Hunter,J.D.,《Matplotlib:二维图形环境》,《计算科学与工程》,第9期,第90-95页(2007年)
[75] 第1卷
[76] van Rossum,G.,《扩展和嵌入python解释器》,报告CS-R9527(1995)
此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配,并且可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。