跳到主要内容
研究论文

基于格子Boltzmann方法的湿土边界层水汽输送高效GPU模拟器:开发和实验验证

作者信息和声明
出版:2023年5月15日出版历史
跳过抽象节

摘要

摘要

我们提出了一种有效的计算方法来模拟多孔介质边界层单相自由流中的组分输运。基于该方法的数值模型通过气候控制风洞与土壤试验台耦合的实验数据进行了验证。所开发的建模方法是基于用于模拟流体流动的格子Boltzmann方法(LBM)和用于求解组分输运的混合网格有限元方法(MHFEM)的组合。这两种方法单独或耦合时都完全在GPU加速器上实现,以利用其计算能力并避免GPU和CPU之间通过PCI-E总线进行缓慢通信所造成的硬件限制。为了利用现代超级计算机上可用的大量计算资源,基于域分解和消息传递接口(MPI)将该实现扩展到分布式多GPU计算。我们描述了计算方法背后的数学细节,主要关注耦合机制。求解器的性能在现代高性能计算系统上进行了演示。本文验证了流动和输运模拟结果,并将其与在暴露于稳定气流的平坦部分饱和土层上进行的实验速度和相对湿度测量值进行了比较。在几个不同的实验场景中,通过向水流中引入立方体钝体,证明了模型的鲁棒性和灵活性。实验测量值可在公共数据集中获得,可作为未来研究的基准。最后,我们讨论了该模型的潜在改进以及未来的实验工作。

参考文献

  1. [1]班杰斯·W·。,哈特曼·R。,Kanschat G.公司。,处理。II–通用面向对象有限元库,ACM事务处理。数学。柔和。 33(4) (2007)24/1第24页,共27页,10.1145/1268776.1268779.谷歌学者谷歌学者数字图书馆数字图书馆
  2. [2]巴斯蒂安·P。,布拉特·M。,工程师C。,德纳A。,克罗夫科恩·R。,Kuttanikkad股份有限公司。,Ohlberger M.公司。,砂光机O。,分布式统一数值环境(DUNE),单位:程序。第19届仿真技术研讨会,2006.谷歌学者谷歌学者
  3. [3]Jasak H。,杰姆科夫A。,图科维奇Z。,等人。,OpenFOAM:一个用于复杂物理模拟的C++库,数值动力学耦合方法国际研讨会,第1000卷,IUC公司,克罗地亚杜布罗夫尼克,2007,第页。120.谷歌学者谷歌学者
  4. [4]普鲁斯·K。,奥尔登堡C.M。,莫里迪斯G。,TOUGH2用户指南版本2,技术代表。 劳伦斯伯克利国家实验室(LBNL),加州伯克利(美国),1999.谷歌学者谷歌学者
  5. [5]Syamlal M。,罗杰斯·W·。,奥布莱恩·T·J。,MFIX文档理论指南,技术代表。 USDOE Morgantown能源技术中心,WV(美国),1993.谷歌学者谷歌学者
  6. [6]ANSYS公司。,ANSYS Fluent 12.0文档,ANSYS公司。,2009,https://www.afs.enea.it/project/neptunius/docs/fluent/.谷歌学者谷歌学者
  7. [7]多物理场耦合分析,COMSOL Multiphysics®简介,多物理场耦合分析,马萨诸塞州伯灵顿,1998.谷歌学者谷歌学者
  8. [8]Krüger T。,Kusumaatmaja H。,库兹明A。,Shardt O。,席尔瓦·G·。,Viggen E.M.公司。,格子Boltzmann方法,施普林格,2017,10.1007/978-3-319-44649-3.谷歌学者谷歌学者
  9. [9]布雷齐·F·。,福汀M。,混合和混合有限元方法,计算数学中的Springer级数,施普林格Verlag,1991,10.1007/978-1-4612-3172-1.谷歌学者谷歌学者
  10. [10]富奇克·R。,克林科夫斯克J。,索洛夫斯克J。,Oberhuber T。,米基什卡·J。,非均质多孔介质组分两相流多维混合有限元方法及其在GPU上的并行实现,计算。物理学。Commun公司。 238(2019)165180,2016年10月10日/j.cpc.2018.12.004.谷歌学者谷歌学者
  11. [11]半径F.A。,Suciu N.公司。,霍夫曼J。,沃格尔A。,Kolditz O。,公园中心。,阿廷格·S。,非均质含水层中污染物运移数值模拟的准确性比较研究,高级水资源。 34(1) (2011)4761,2016年10月10日/j.advwatres.2010.09.012.谷歌学者谷歌学者
  12. [12]于C。,马拉克普尔K。,惠格J.M。,混合杂交有限元法与标准有限元法在超大变形溶胀模拟中的比较,计算。机械。 66(2) (2020)287309,10.1007/s00466-020-01851-z.谷歌学者谷歌学者数字图书馆数字图书馆
  13. [13] 消息传递接口论坛,MPI:消息传递接口标准3.1版(2015年6月) :https://www.mpi-forum.org/docs/mpi-3.1/mpi31-report.pdf.谷歌学者谷歌学者
  14. [14]Brutsaert W。,蒸发到大气中:理论、历史和应用,施普林格,1982.谷歌学者谷歌学者
  15. [15]石原Y。,Shimojuma E.公司。,原田H。,裸露土壤下的水蒸气传输,蒸发受湍流表面风的影响,J.水文学。 131(1992)63104,10.1016/0022-1694(92)90213-F.谷歌学者谷歌学者
  16. [16]特劳茨A.C。,伊兰加塞卡尔T.H。,豪灵顿S。,持续地上气流存在下裸土蒸发动力学研究的实验测试规模考虑因素,水资源。物件。 54(11) (2018)89638982,2018年10月29日WR023102.谷歌学者谷歌学者
  17. [17]特劳茨A.C。,伊兰加塞卡尔T.H。,豪灵顿S。,磁韩A。,连续尺度多孔介质传热传质模型对应用大气强迫数据空间离散长度尺度的敏感性,水资源。物件。 55(4) (2019)35203540,2018年10月29日WR023923.谷歌学者谷歌学者
  18. [18]特劳茨A.C。,伊兰加塞卡尔T.H。,罗德里格斯-伊图尔贝一世。,水分限制环境中竞争和促进在植物间距流体动力学中的作用,程序。国家。阿卡德。科学。 114(35) (2017)93799384,10.1073/pnas.1706046114.谷歌学者谷歌学者
  19. [19]特劳茨A.C。,伊兰加塞卡尔T.H。,罗德里格斯-伊图尔贝一世。,见鬼K。,赫尔米格·R。,使用植物模拟物研究土壤-土壤-大气耦合过程的实验方法的发展,水资源。物件。 53(4) (2017)33193340,10.1002/2016WR019884.谷歌学者谷歌学者
  20. [20]特劳茨A.C。,伊兰加塞卡尔T.H。,两个石灰石块的蒸发蒸腾(工厂模拟)实验,门德利数据 2022,10.17632/6油炸4xzgh.1.谷歌学者谷歌学者
  21. [21]狮子P.-L。,马斯穆迪N。,粘性可压缩流体的不可压缩极限,数学杂志。Pures应用程序。 77(6) (1998)585627,10.1016/S0021-7824(98)80139-6.谷歌学者谷歌学者
  22. [22]Desjardins B。,格雷尼尔E。,全空间粘性可压缩流动的低马赫数极限,R.Soc.伦敦。程序。,序列号。A、 数学。物理学。工程科学。 455(1986) (1999)22712279,10.1098/rspa.1999.0403.谷歌学者谷歌学者
  23. [23]鸟R.B。,斯图尔特·W.E。,Lightfoot E.N.公司。,运输现象,第2版,John Wiley&Sons公司。,2002,10.1115/1.1424298.谷歌学者谷歌学者
  24. [24]工程工具箱:空气-动力和运动粘度.[在线;访问日期:2022年3月27日](2003) https://www.engineeringtoolbox.com/air-absolute-kinematic-visocity-d_601.html.谷歌学者谷歌学者
  25. [25]Massman W.J.公司。,H2O、CO2、CH4、CO、O3、SO2、NH3、N2O、NO和NO2在STP附近空气、O2和N2中的分子扩散率综述,大气。环境。 32(6) (1998)11111127,10.1016/S1352-2310(97)00391-9.谷歌学者谷歌学者
  26. [26]康S.K。,哈桑Y.A。,格子Boltzmann方法中格子模型在壁面湍流模拟中的作用,J.计算。物理学。 232(1) (2013)100117,2016年10月10日/j.jcp.2012.07.023.谷歌学者谷歌学者数字图书馆数字图书馆
  27. [27]盖尔·M·。,Schönherr M。,帕斯夸利A。,克拉夫茨克·M·。,三维累积格子Boltzmann方程:理论与验证,计算。数学。应用。 70(4) (2015)507547,2016年10月10日/j.camwa.2015.05.001.谷歌学者谷歌学者
  28. [28]Kumar P。,库彻K。,莫纳·M·。,Radespiel R。,克拉夫茨克·M·。,盖尔·M·。,用累积格子Boltzmann DNS/LES和实验验证多孔平板湍流的VRANS模型,J.多孔介质 21(5) (2018),10.1615/JPorMedia.v21.i5.60.谷歌学者谷歌学者
  29. [29]彭C。,日内瓦N。,郭Z。,王丽萍。,用格子Boltzmann方法直接数值模拟湍流管道流动,J.计算。物理学。 357(2018)1642,2016年10月10日/j.jcp.2017.11.040.谷歌学者谷歌学者数字图书馆数字图书馆
  30. [30]威特曼·M·。,泽塞T。,海格·G。,Wellein G。,大规模格子Boltzmann模拟的区域分解和局部优化,计算。流体 80(2013)283289,2016年10月10日/j.compfluid.2012.002007.谷歌学者谷歌学者
  31. [31]扎基洛夫A。,佩雷佩尔基纳A。,列夫琴科五世。,基尔科夫S。,流技术:揭示格子Boltzmann方法的自然并发性,J.超级计算机。 77(10) (2021)1191111929,2007年10月17日/11227-021-03762-z.谷歌学者谷歌学者数字图书馆数字图书馆
  32. [32]富奇克·R。,艾希勒·P·。,斯特拉卡R。,鲍斯普。,克林科夫斯克J。,Oberhuber T。,二维和三维浸没边界格子Boltzmann方法的最优节点间距,计算。数学。应用。 77(4) (2019)11441162,2016年10月10日/j.camwa.2018.10.045.谷歌学者谷歌学者数字图书馆数字图书馆
  33. [33]艾希勒·P·。,Fučík R。,斯特拉卡R。,流体-结构相互作用的浸没边界格子Boltzmann方法计算研究,离散连续。动态。系统。,序列号。S公司 14() (2021)819,10.3934/dcdss.2020349.谷歌学者谷歌学者
  34. [34]艾希勒·P·。,福卡五世。,富奇克·R。,粗糙表面湍流的累积格子Boltzmann模拟,计算。数学。应用。 92(2021)3747,2016年10月10日/j.camwa.2021.03.016.谷歌学者谷歌学者
  35. [35]盖尔·M·。,帕斯卡利A。,Schönherr M。,四阶精确扩散累积量格子Boltzmann方法的参数化——第一部分:推导和验证,J.计算。物理学。 348(2017)862888,2016年10月10日/j.jcp.2017.05.040.谷歌学者谷歌学者数字图书馆数字图书馆
  36. [36]盖尔·M·。,帕斯卡利A。,Schönherr M。,四阶精确扩散累积点阵Boltzmann方法的参数化第二部分:阻力危机下绕球体流动的应用,J.计算。物理学。 348(2017)889898,2016年10月10日/j.jcp.2017.07.004.谷歌学者谷歌学者数字图书馆数字图书馆
  37. [37]Latt J。,肖伯德B。,马拉菠菜。,Deville M。,米克勒A。,格子Boltzmann方法中的直线速度边界,物理学。版本E 77(2008),10.1103/物理版次E.77.056703.谷歌学者谷歌学者
  38. [38]穆罕默德A.A。,Succi S。,格子Boltzmann流体动力学模拟中平衡边界条件的注记,欧洲物理学。J.规格顶部。 171(1) (2009)213221,10.1140/epjst/e2009-01031-9.谷歌学者谷歌学者
  39. [39]豪斯曼M。,巴雷托交流。,Kouyi G.L.公司。,Rivière N。,Nirschl H。,克劳斯·M.J。,用格子Boltzmann方法对高雷诺数湍流通道流动的大涡模拟与壁面模型耦合——在科里奥利质量流量计中的应用,计算。数学。应用。 78(10) (2019)32853302,2016年10月10日/j.camwa.2019.04.033.谷歌学者谷歌学者数字图书馆数字图书馆
  40. [40]白色F.M。,粘性流体流动,麦格劳-希尔,纽约,2005.谷歌学者谷歌学者
  41. [41]戴维森·L·。,比尔森·M·。,混合LES-RANS使用合成湍流涨落在界面区域内强制,国际热流学杂志 27(6) (2006)10281042,2016年10月10日/j.ijheatfluidflow.2006.02.025.谷歌学者谷歌学者
  42. [42]戴维森·L·。,使用各向同性合成波动作为非定常模拟的入口边界条件,高级申请。流体力学。 1(2007)135.谷歌学者谷歌学者
  43. [43]戴维森·L·。,混合LES-RANS:包括再循环在内的流动入口边界条件,单位:第五届湍流和剪切流现象国际研讨会,Begel House公司。,2007,第页。689694.谷歌学者谷歌学者
  44. [44]戴维森·L·。,流体力学、湍流和湍流建模,2018.谷歌学者谷歌学者
  45. [45]索洛夫斯克J。,富奇克·R。,普拉宾M.R。,伊兰加塞卡尔T.H。,米基什卡·J。,相间传质对浅层含水层中结构捕集的气态二氧化碳衰减的尺寸效应,J.计算。物理学。 405(2020),2016年10月10日/j.jcp.2019.109178.谷歌学者谷歌学者数字图书馆数字图书馆
  46. [46]阿斯卡尔·A.H。,伊兰加塞卡尔T.H。,特劳茨A。,索洛夫斯克J。,张勇。,富奇克·R。,探索源条件不确定性对CO2地质封存中远场盐水泄漏羽流预测的影响:将中等规模实验室测试与数值模拟相结合,水资源。物件。 57(9) (2021),10.1029/2021WR029679.谷歌学者谷歌学者
  47. [47]布伦纳·F·。,半径F.A。,克纳布纳P。,求解运输问题的迎风混合有限元方法分析,SIAM J.数字。分析。 52(1) (2014)83102,10.1137/130908191.谷歌学者谷歌学者数字图书馆数字图书馆
  48. [48]拉尔森M.G。,尼克拉森·A.J。,基于不连续富集的连续Galerkin方法的守恒通量,卡尔科洛 41(2) (2004)6576,2007年10月10日/BF02637255.谷歌学者谷歌学者
  49. [49]太阳S。,惠勒M.F。,与运输兼容的速度数据投影,计算。方法应用。机械。工程师。 195(7–8) (2006)653673,2016年10月10日/j.cma.2005.02.011.谷歌学者谷歌学者
  50. [50]Kees首席执行官。,Farting M.W.公司。,道森公司。,变饱和流的局部保守稳定有限元方法,计算。方法应用。机械。工程师。 197(51–52) (2008)46104625,2016年10月10日/j.cma.2008.06.005.谷歌学者谷歌学者
  51. [51]Oberhuber T。,克林科夫斯克J。,富奇克·R。,TNL:现代并行体系结构的数字库,理工学院学报。 61() (2021)122134,10.14311/AP.2021.61.0122.谷歌学者谷歌学者
  52. [52]NVIDIA公司,CUDA工具包文档,11.6版,英伟迪亚,2022,https://docs.nvidia.com/cuda/archive/11.6.0/.谷歌学者谷歌学者
  53. [53]韦尔福德B。,关于计算平方和乘积修正和的方法的注记,技术计量学 4() (1962)419420,10.1080/00401706.1962.10490022.谷歌学者谷歌学者
  54. [54]林R.F。,几种样本均值和方差计算算法的比较,美国统计协会。 69(348) (1974)859866,10.1080/01621459.1974.10480219.谷歌学者谷歌学者
  55. [55]Chan T.F.先生。,Golub G.H.公司。,勒韦克·R.J。,计算样本方差的算法:分析和建议,美国统计局。 37() (1983)242247,10.1080/00031305.1983.10483115.谷歌学者谷歌学者
  56. [56]杜奇J.C.I。,寻找水文规律,水资源。物件。 22(1986)46秒58秒,10.1029/WR022i09Sp0046S.谷歌学者谷歌学者
  57. [57]Wolfe股份有限公司。,尼克·W.G。,稀疏植被在风蚀中的保护作用,掠夺。物理学。地理。 17(1) (1993)5068,10.1177/030913339301700104.谷歌学者谷歌学者

建议

评论

登录选项

检查您是否可以通过登录凭据或您的机构访问本文。

登录

完全访问权限

  • 文章指标

    • 下载量(最近12个月)0
    • 下载次数(最近6周)0

    其他指标