沃纳·威尔第;皮埃尔·凯斯特纳;阿兰·卡塔拉德 具有相变的两相流格子Boltzmann方法的性能可移植性。 （英语） Zbl 1506.76134号 计算。方法应用。机械。工程师。 370，文章ID 113266，29 p.（2020）. 摘要：使用格子Boltzmann方法（LBM）模拟单相或两相流的数值代码被广泛编译并在图形处理单元上运行。然而，这些计算单元必须使用适合这些体系结构的低级语言（例如。CUDA公司对于GPU NVIDIA\（\circledR\）或开放运算语言). 在本文中，我们将重点放在LBM的性能可移植性上，即使用高级抽象编写LB算法的可能性，同时在多核x86、GPU NVIDIA（\circledR）、ARM等多种体系结构上保持高效。为此，LBM的实施是通过开发一个独特的代码来实现的，LBM_saclay公司用C++语言编写，结合库在高性能计算环境中实现性能可移植性。本文使用LBM对具有相变的两相流问题的相场模型进行了模拟。数学模型由不可压缩的Navier-Stokes方程和保守的Allen-Cahn模型耦合而成。该模型最初是在文献中针对不混溶的二元流体开发的，在这里进行了扩展，以模拟液体和气体界面处发生的相变。为此，在热方程中添加了涉及相场时间导数的源项。在相场方程中，添加了一个源项来近似界面处的批量生产速率。进行了几次验证，以逐步检查完整模型的实现情况。最后，针对膜沸腾的物理问题，比较了CPU和GPU平台上的计算时间。 引用于3文件 MSC公司： 76米28 粒子法和晶格气体法 65亿75 涉及偏微分方程的初值和初边值问题的概率方法、粒子方法等 65日元10 特定类别建筑的数值算法 76Txx型 多相多组分流动 关键词：晶格玻尔兹曼方法;相场模型;具有相变的两相流;性能可移植性;Kokkos图书馆;LBM_saclay公司 软件：科科斯;github;旗鱼;CUDA公司 PDF格式BibTeX公司 XML格式引用 \textit{W.Verdier}等人，计算。方法应用。机械。工程370，文章ID 113266，29 p.（2020；Zbl 1506.76134） 全文： 内政部 arXiv公司 参考文献： [1] Krueger，T。；Kusumaatmaja，H。；库兹明，A。；沙特，O。；席尔瓦，G。；Viggen，E.，《格子Boltzmann方法：原理与实践》，《物理学研究生论文》（2016），斯普林格出版社 [2] 郭，Z。；Shu，C.，格子Boltzmann方法及其在工程中的应用（2013），世界科学·Zbl 1278.76001号 [3] 李伟（Li，W.）。；魏，X。；考夫曼，A.，在图形硬件上实现晶格玻尔兹曼计算，Vis。计算。，19, 7, 444-456 (2003) [4] Tölke，J.，使用nVIDIA开发的计算统一设备架构实现格子Boltzmann内核，计算。视觉。科学。，13, 1, 29 (2008) [5] 库兹尼克，F。；奥布莱希特，C。；Rusaouen，G。；Roux，J.-J.，使用GPU计算处理器的基于LBM的流模拟，计算。数学。申请。，59, 7, 2380-2392 (2010) ·Zbl 1193.76112号 [6] 周，H。；莫，G。；Wu，F。；赵，J。；鲁伊，M。；Cen，K.，曲面边界流动的格子Boltzmann方法的GPU实现，计算。方法应用。机械。工程，225-228，65-73（2012）·Zbl 1253.76004号 [7] 奥布莱希特，C。；库兹尼克，F。；图兰乔，B。；Roux，J.-J.，格子Boltzmann方法的Multi-GPU实现，计算。数学。申请。，65，252-261（2013），工程与科学中的介观方法特刊（ICMMES-2010，加拿大埃德蒙顿）·Zbl 1268.76048号 [8] Januszewski，M。；Kostur，M.，Sailfish：晶格玻尔兹曼方法的灵活多GPU实现，Comput。物理学。Comm.，185，9，2350-2368（2014）·Zbl 1344.76004号 [9] 爱德华兹，H.C。；Trott，C.R。；Sunderland，D.，Kokkos：通过多态内存访问模式实现多核性能可移植性，J.Parallel Distrib.Comput。，74, 12, 3202-3216 (2014) [10] Eichstädt，J。；格林，M。；特纳，M。；佩罗，J。；Moxey，D.，《使用架构相关编程模型加速高阶网格优化》，计算。物理学。Comm.，229，36-53（2018） [11] 黄，H。；Sukop，M。；Lu，X.-Y.，多相格子Boltzmann方法。《理论与应用》（2015），威利父子公司 [12] 李强。；罗，K。；康（Kang，Q.）。；何毅。；陈，Q。；Liu，Q.，多相流动和相变传热的格子Boltzmann方法，Prog。能源燃烧。科学。，52, 62-105 (2016) [13] Gunstensen，A.K。；罗斯曼，D.H。；Zaleski，S。；Zanetti，G.，《不混溶流体的格子Boltzmann模型》，Phys。版本A，43，4320-4327（1991） [14] 珊，X。；Chen，H.，模拟多相多组分流动的格子Boltzmann模型，Phys。E版，47，1815-1819（1993） [15] 珊，X。；Chen，H.，通过晶格Boltzmann方程模拟非理想气体和液气相变，物理学。E版，49，2941-2948（1994） [16] 斯威夫特，M.R。；Orlandini，E。；奥斯本，W.R。；Yeomans，J.M.，《液气和二元流体系统的格子Boltzmann模拟》，Phys。E版，54，5041-5052（1996） [17] Jackmin，D.，《使用相场建模计算两相navier-Stokes流》，J.Compute。物理。，155, 96-127 (1999) ·Zbl 0966.76060号 [18] 袁，P。；Schaefer，L.，格子Boltzmann模型中的状态方程，物理学。流体，18042101（2006）·Zbl 1185.76872号 [19] 李，X。；Cheng，P.，Lattice Boltzmann对具有亲水斑点的疏水表面上从滴状冷凝到膜状冷凝转变的模拟，《国际传热杂志》，110，710-722（2017） [20] 李，M。；胡贝尔，C。；Mu，Y。；Tao，W.，Lattice Boltzmann《不凝气体存在下冷凝的模拟》，《国际热质传递杂志》，1091004-1013（2017） [21] 莱克莱尔，S。；雷吉奥，M。；Trépanier，J.-Y.，不混溶两相流格子Boltzmann模型的两个重氯算子的数值评估，Appl。数学。型号。，36, 5, 2237-2252 (2012) ·兹比尔1243.76079 [22] 莱克莱尔，S。；佩莱林，N。；雷吉奥，M。；Trépanier，J.-Y.，在一类格子Boltzmann模型中模拟具有可变密度比的不混溶多相流的增强平衡分布函数，Int.J.Multiph。流量。，57, 159-168 (2013) [23] 安德森，D。；麦克法登，G。；Wheeler，A.，《流体力学中的扩散界面方法》，年。流体力学版次。，30, 139-165 (1998) ·Zbl 1398.76051号 [24] 卡恩，J。；Hilliard，J.，《非均匀系统的自由能》。i.界面自由能，J.Chem。物理。，28, 2, 258-267 (1958) ·Zbl 1431.35066号 [25] Jasnow，D。；Viñals，J.，《热毛细流动的粗粒度描述》，《物理学》。流体，8，3，660-669（1996）·Zbl 1025.76521号 [26] V.Kendon。；Cates，M。；帕戈纳巴拉加，I。；德斯普拉特，J.-C。；Bladon，P.，对称二元流体混合物三维旋节分解中的惯性效应：晶格Boltzmann研究，J.流体力学。，440, 147-203 (2001) ·Zbl 1049.76052号 [27] 郑浩。；舒，C。；Chew，Y.，大密度比多相流的格子Boltzmann模型，J.Compute。物理。，218, 353-371 (2006) ·Zbl 1158.76419号 [28] Lee，T。；Liu，L.，Lattice Boltzmann模拟干表面上的微尺度液滴冲击，J.Compute。物理。，229, 8045-8063 (2010) ·Zbl 1426.76603号 [29] Zu，Y.Q。；He，S.，基于相场的格子Boltzmann模型，用于具有密度和粘度对比的不可压缩二元流体系统，Phys。E版，87，043301（2013） [30] 孙，Y。；Beckermann，C.，使用相场方程的Sharp界面跟踪，J.Comput。物理。，220, 2, 626-653 (2007) ·Zbl 1228.76110号 [31] Chiu，P.-H。；Lin，Y.-T.，求解不可压缩两相流的保守相场方法，J.Compute。物理。，230, 1, 185-204 (2011) ·Zbl 1427.76201号 [32] 盖尔，M。；Fakhari，A。；Lee，T.，界面跟踪方程的保守相场晶格Boltzmann模型，Phys。版本E，91，063309（2015） [33] Fakhari，A。；靠垫，D。；Luo，L.-S.，多相流的加权多重松弛时间格子Boltzmann方法及其在部分聚结叶栅中的应用，J.Compute。物理。，341, 22-43 (2017) ·Zbl 1376.76047号 [34] 米切尔，T。；莱昂纳迪，C。；Fakhari，A.，用于研究高密度比下不混溶流体的三维相场晶格Boltzmann方法的发展，国际J.Multiph。流量。，107, 1-15 (2018) [35] Lee，D。；Kim，J.，保守的艾伦-卡恩和卡恩-希利亚德方程的比较研究，数学。计算。模拟，119，35-56（2016）·Zbl 07313590号 [36] Wang，H.L。；Chai，Z.H。；史，公元前。；Liang，H.，allen-cahn和cahn-hilliard方程格子Boltzmann模型的比较研究，物理学。E版，94033304（2016） [37] 东，Z。；李伟（Li，W.）。；Song，Y.，Lattice Boltzmann模拟上升的蒸汽泡通过过热液体的生长和变形，Numer。热传输。A： 申请。，55, 4, 381-400 (2009) [38] 萨法里，H。；Rahimian，M.H。；Krafczyk，M.，两相流体流动热相变数值模拟的扩展格子Boltzmann方法，Phys。E版，88，013304（2013） [39] Kharangate，C.R。；Mudawar，I.，《沸腾和冷凝计算研究综述》，《国际传热杂志》，108，1164-1196（2017） [40] Folch，R。；Casademunt，J。；埃尔南德斯·马查多，A。；Ramírez-Piscina，L.，具有任意粘度对比度的hele-shaw流动的相场模型。i.理论方法，Phys。版本E，60，1724-1733（1999） [41] Jamet，D。；Misbah，C.，囊泡相场模型的热力学一致图：消除表面张力，Phys。修订版E，78041903（2008） [42] 小林，R.，枝晶生长的建模和数值模拟，Physica D，63，3，410-423（1993）·兹比尔0797.35175 [43] Amaya-Bower，L。；Lee，T.，使用格子Boltzmann方法研究中等雷诺数的单气泡上升动力学，计算与流体，39，7，1191-1207（2010）·Zbl 1242.76264号 [44] Kim，J.，扩散界面模型的连续表面张力公式，J.Compute。物理。，204, 2, 784-804 (2005) ·Zbl 1329.76103号 [45] Delhaye，J.，两相系统中的跳跃条件和熵源。本地即时配方，国际J.Multiph。流量。，1, 3, 395-409 (1974) ·Zbl 0358.76066号 [46] Juric，D。；Tryggvason，G.，《沸腾流动的计算》，《国际多道物理杂志》。流量。，24, 3, 387-410 (1998) ·Zbl 1121.76455号 [47] 萨法里，H。；Rahimian，M.H。；Krafczyk，M.，基于相场多相晶格Boltzmann方法的液滴蒸发一致模拟，Phys。E版，90，033305（2014） [48] Begmohammadi，A。；法哈扎德，M。；Rahimian，M.H.，使用格子Boltzmann方法模拟高密度比下的池沸腾和周期性气泡释放，国际通讯社。热质传递，61，78-87（2015） [49] Karma，A。；Rappel，W.-J.，二维和三维枝晶生长的定量相场建模，物理。E版，57、4、4323-4349（1998年）·Zbl 1086.82558号 [50] Ginzburg，I.，一般平流和各向异性扩散方程的平衡型和连接型格子Boltzmann模型，高级水资源。，28, 11, 1171-1195 (2005) [51] D’Humières，D.，广义格子Boltzmann方程，（Rarefied Gas Dynamics:理论与模拟（1992）），450-458，Prog。宇航员。飞行员。59 [52] d’Humières，d。；金兹堡，I。；Krafczyk，M。；拉勒曼德，P。；Luo，L.-S.，三维多重弛豫时间晶格Boltzmann模型，Philos。事务处理。R.Soc.伦敦。数学。物理学。工程科学。，360, 437-451 (2002) ·兹比尔1001.76081 [53] 何，X。；陈，S。；Doolen，G.D.，《不可压缩极限下格子Boltzmann方法的新型热模型》，J.Compute。物理。，1461282-300（1998年）·Zbl 0919.76068号 [54] Inamuro，T。；Ogata，T。；田岛，S。；Konishi，N.，《大密度差不可压缩两相流的格子Boltzmann方法》，J.Compute。物理。，198, 2, 628-644 (2004) ·Zbl 1116.76415号 [55] Chorin，A.J.，《解决不可压缩粘性流动问题的数值方法》，J.Compute。物理。，2, 1, 12-26 (1967) ·Zbl 0149.44802号 [56] 何，X。；Luo，L.-S.，不可压缩navier-Stokes方程的格子Boltzmann模型，J.Stat.Phys。，88, 3/4, 927-944 (1997) ·Zbl 0939.82042号 [57] 穆罕默德·沙德，M。；Lee，T.，使用尖锐界面能量解算器对沸腾进行相场晶格Boltzmann建模，Phys。E版，96，013306（2017） [58] Ginzburg，I.，格子Boltzmann模型的通用边界条件及其在平流和各向异性色散方程中的应用，高级水资源。，28, 11, 1196-1216 (2005) [59] Lee，T。；Fischer，P.，《消除非理想气体晶格Boltzmann方程法中的寄生电流》，Phys。E版，74046709（2006） [60] Lee，T.，不可压缩性对消除二元流体格子Boltzmann方程法中寄生电流的影响，计算。数学。申请。，58, 987-994 (2009) ·Zbl 1189.76414号 [61] Lee，T。；Lin，C.-L.，用于模拟高密度比下不可压缩两相流的格子Boltzmann方程的稳定离散化，J.Compute。物理。，206, 1, 16-47 (2005) ·Zbl 1087.76089号 [62] Fakhari，A。；米切尔，T。；莱昂纳迪，C。；Bolster，D.，《改进的相场晶格局域性——高密度比下不混溶流体的Boltzmann模型》，Phys。E版，96，053301（2017） [63] kokkos库的兼容性（2020年），https://github.com/kokkos/kokkos/wiki/编译（2020年5月12日可访问Web链接） [64] Calore，E。；加巴纳，A。；斯基法诺，S.F。；Tripiccione，R.，《在格子Boltzmann应用中使用骑士着陆处理器的早期经验》（Wyrzykowski，R.；Dongarra，J.；Deelman，E.；Karczewski，K.，《并行处理与应用数学》（2018），斯普林格国际出版：斯普林格国际出版商会），519-530 [65] T.Z.alesak，S.，流体的全多维通量校正传输算法，J.Compute。物理。，33133-362（1979年）·Zbl 0416.76002号 [66] 哈恩，D。；Özisik，M.，《导热》，752（2012），威利父子公司 [67] Esmaeeli，A。；Tryggvason，G.，《复杂几何体中沸腾计算的前沿跟踪方法》，国际J·多波。流量。，30、7、1037-1050（2004），G.Yadigarouglu教授65岁生日纪念论文集·Zbl 1201.76275号 [68] 潘迪，V。；比斯瓦斯，G。；Dalal，A.，过热和电场对饱和膜沸腾的影响，物理学。流体，28，5，052102（2016） [69] Begmohammadi，A。；拉希米安，M。；法哈扎德，M。；Hatani，M.A.，使用格子Boltzmann方法对单模和多模膜沸腾进行数值模拟，计算。数学。申请。，71, 9, 1861-1874 (2016) ·Zbl 1443.65259号 [70] 胡，A。；Liu，D.，用改进的混合格子Boltzmann模型对壁面沸腾传热进行二维模拟，应用。热量。工程，159113788（2019） [71] 新泽西州辛格。；Premachandran，B.，水平波浪壁上膜沸腾的数值研究，《国际传热传质杂志》，150，119371（2020） 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。