M.拉赫什。;杨，L。;胡，W。;内格鲁特，D。 关于使用多体动力学技术来模拟流体动力学和流固相互作用问题。 （英语） Zbl 1490.76168号 多体系统。动态。 53，编号1，29-57（2021）. 摘要：本文将基于多体动力学的流体动力学问题解决方案与计算流体动力学（CFD）社区使用的两种基于拉格朗日的技术进行了比较。基于多体动力学的解决方案具有两个显著属性：通过双边运动约束强制执行不可压缩条件，通过单边运动约束处理与固相的耦合。基于多体动力学的解决方案，在此称为运动约束光滑粒子流体动力学（KCSPH）方法，是解决CFD问题的拉格朗日方法。它依靠光滑粒子流体动力学（SPH）方法离散Navier-Stokes方程中的空间微分算子，并依靠多体动力学的Newton-Euler方程在时间上向前对流SPH粒子。通过与CFD界最常用的两种SPH算法（弱可压缩SPH（WCSPH）和隐式SPH（ISPH）方法）的性能比较，我们证明了基于多体动力学的方法是高效和准确的。比较是结合四个试验进行的：不可压缩性基准试验、溃坝试验、浮动圆筒试验和晃动水槽试验。我们的结论是，对于具有复杂/移动边界的流体-固体相互作用（FSI）问题，KCSPH是传统CFD方法的可靠替代方案。本文使用的解算器和模型在一个名为Chrono的开源软件中公开可用；这些实现使用GPU（用于WCSPH和ISPH）和多核CPU（用于KCSPH）并行计算。 引用于4文件 MSC公司： 76米28 粒子法和晶格气体法 74层10 流固相互作用（包括气动和水弹性、孔隙度等） 70E55型 多体系统动力学 关键词：多体动力学求解技术;运动约束;并行计算;光滑粒子流体力学;流固相互作用 软件：梅沙赫;github;Chrono公司;双球面物理学;CUDA公司 PDF格式BibTeX公司 XML格式引用 \textit{M.Rakhsha}等人，《多体系统》。动态。53，编号1，29-57（2021；Zbl 1490.76168） 全文： 内政部 参考文献： [1] Acary，V.公司。；Brogliato，B.，《非光滑动力系统的数值方法：在力学和电子学中的应用》（2008），柏林：施普林格出版社，柏林·Zbl 1173.74001号 ·doi:10.1007/978-3-540-75392-6 [2] 阿达米，S。；胡，X。；Adams，N.，光滑粒子流体动力学的广义壁面边界条件，J.Compute。物理。，231, 21, 7057-7075 (2012) ·doi:10.1016/j.jcp.2012.05.005 [3] 阿达米，S。；胡，X。；Adams，N.，平滑粒子流体动力学的输运速度公式，J.Compute。物理。，241, 292-307 (2013) ·Zbl 1349.76659号 ·doi:10.1016/j.jcp.2013.01.043 [4] Anitescu，M.，基于优化的非光滑刚体多体动力学仿真，数学。程序。，105, 1, 113-143 (2006) ·Zbl 1085.70008号 ·doi:10.1007/s10107-005-0590-7 [5] Anitescu，M。；Hart，G.D.，《带关节、接触和摩擦的刚性多体动力学的约束稳定时间步长方法》，国际期刊Numer。方法工程，60，14，2335-2371（2004）·Zbl 1075.70501号 ·doi:10.1002/nme.1047 [6] Antoci，C。；加拉蒂，M。；Sibilla，S.，SPH流体-结构相互作用的数值模拟，计算。结构。，85, 11, 879-890 (2007) ·doi:10.1016/j.compstruc.2007.01.002 [7] 安托诺，M。；Colagrossi，A。；Marrone，S.，弱可压缩SPH格式中的数值扩散项，计算。物理学。社区。，183, 12, 2570-2580 (2012) ·Zbl 1507.76152号 ·doi:10.1016/j.cpc.2012.07.006 [8] Asai，M.，Aly，A.M.，Sonoda，Y.，Sakai，Y.：通过放松密度不变条件的稳定不可压缩SPH方法。J.应用。数学。2012 (2012). doi:10.1155/2012/139583·Zbl 1244.76101号 [9] 贝克尔，M。；Teschner，M.，自由表面流的弱可压缩SPH，2007年ACM SIGGRAPH/欧洲制图计算机动画研讨会论文集，209-217（2007） [10] 贝尔泰尔斯·德库比斯，F。；Cadoux，F。；Daviet，G。；Acary，V.，捕捉纤维组件中精确库仑摩擦的非光滑牛顿解算器，ACM Trans。图表。，30, 1, 6 (2011) ·数字对象标识代码：10.1145/1899404.1899410 [11] 博丹，K。；拉库西耶，C。；Servin，M.，约束流体，IEEE Trans。视觉。计算。图表。，18, 3, 516-526 (2012) ·doi:10.1109/TVCG.2011.29 [12] Bonet，J。；Lok，T.S.，光滑粒子流体动力学公式的变分和动量保持方面，计算。方法应用。机械。工程，180，1-2，97-115（1999）·Zbl 0962.76075号 ·doi:10.1016/S0045-7825（99）00051-1 [13] Chorin，A.J.，Navier-Stokes方程的数值解，数学。计算。，22, 104, 745-762 (1968) ·Zbl 0198.50103号 ·doi:10.1090/S0025-5718-1968-0242392-2 [14] Colagrossi，A。；Landrini，M.，用光滑粒子流体动力学对界面流动进行数值模拟，J.Compute。物理。，191, 2, 448-475 (2003) ·Zbl 1028.76039号 ·doi:10.1016/S0021-9991（03）00324-3 [15] 克雷斯波，A.J。；Domínguez，J.M。；罗杰斯，B.D。；戈梅斯·盖斯特拉，M。；Longshaw，S。；卡内拉斯，R。；瓦康迪奥，R。；巴雷罗，A。；García-Feal，O.，DualSPHysics:基于平滑粒子流体动力学（SPH）的开源并行CFD求解器，计算。物理学。社区。，187, 204-216 (2015) ·Zbl 1348.76005号 ·doi:10.1016/j.cpc.2014.10.004 [16] 康明斯，S.J。；Rudman，M.，《SPH投影法》，J.Compute。物理。，152, 2, 584-607 (1999) ·Zbl 0954.76074号 ·doi:10.1006/jcph.1999.6246 [17] 坎达尔，P。；Strack，O.，颗粒组件的离散元模型，《岩土工程》，29，47-65（1979）·doi:10.1680/geot.1979.29.1.47 [18] Dodge，F.T.，《新的“移动容器中液体的动态行为”》（2000年），圣安东尼奥：西南研究所，圣安东尼奥市 [19] Ellero，M。；塞拉诺，M。；Español，P.，《不可压缩光滑粒子流体动力学》，J.Compute。物理。，226, 2, 1731-1752 (2007) ·Zbl 1121.76050号 ·doi:10.1016/j.jp.2007.06.019 [20] 埃斯帕诺，P。；Revenga，M.，平滑耗散粒子动力学，物理学。E版，67，2（2003）·doi:10.1103/PhysRevE.67.026705 [21] Fatehi，R。；Manzari，M.T.，平滑粒子流体动力学中的误差估计和二阶导数的新方案，计算。数学。申请。，61, 2, 482-498 (2011) ·Zbl 1211.76089号 ·doi:10.1016/j.camwa.2010.11.028 [22] Gingold，R.A。；Monaghan，J.J.，《平滑粒子流体动力学理论及其在非球形恒星中的应用》，蒙大拿州。不是。R.阿斯顿。《社会学杂志》，181，1375-389（1977）·Zbl 0421.76032号 ·doi:10.1093/mnras/181.3.375 [23] Gurtin，M.E。；油炸，E。；Anand，L.，《连续统的力学和热力学》（2010），剑桥：剑桥大学出版社，剑桥·doi:10.1017/CBO9780511762956 [24] 海尔，E。；诺塞特，S。；Wanner，G.，《求解常微分方程I：非刚性问题》（2009），柏林：施普林格出版社，柏林·Zbl 1185.65115号 [25] Haug，E.J.，《机械系统的计算机辅助运动学和动力学》第一卷（1989年），《英格伍德悬崖：普伦蒂斯·霍尔》，《英格伍德悬崖》 [26] Heyn，T。；Anitescu，M。；塔索拉。；Negrut，D.，《使用Krylov子空间和谱方法解决多体接触动力学模拟中的互补问题》，Int.J.Numer。方法工程，95，7，541-561（2013）·Zbl 1352.74206号 ·doi:10.1002/nme.4513 [27] 侯赛尼，S.M。；Feng，J.J.，用SPH计算不可压缩流的压力边界条件，J.Compute。物理。，230, 19, 7473-7487 (2011) ·Zbl 1408.76413号 ·doi:10.1016/j.jcp.2011.06.013 [28] 胡，X.Y。；Adams，N.A.，《不可压缩多相SPH方法》，J.Compute。物理。，227, 1, 264-278 (2007) ·Zbl 1126.76045号 ·doi:10.1016/j.jcp.2007.07.013 [29] Hu，W.，Pan，W..，Rakhsha，M.，Negrut，D.：对SPH技术进行概述，以保持2D和3D流体动力学的二阶收敛。技术代表TR-2016-14，基于仿真的工程实验室，威斯康星大学麦迪逊分校（2016） [30] 胡，W。；潘·W。；拉赫什，M。；田强。；胡，H。；Negrut，D.，一种一致的多分辨率平滑粒子流体动力学方法，计算。方法应用。机械。工程，324，278-299（2017）·Zbl 1439.76128号 ·doi:10.1016/j.cma.2017.06.010 [31] Hu，W.，Rakhsha，M.，Yang，L.，Kamrin，K.，Negrut，D.：使用连续体表示和SPH方法模拟颗粒材料动力学及其与运动固体的双向耦合。应用力学和工程中的计算机方法-正在审查00（2021）·Zbl 1502.76120号 [32] 休斯·J·P。；Graham，D.I.，自由表面水流的不可压缩和弱可压缩SPH模型的比较，J.Hydraul。研究，48，S1，105-117（2010）·doi:10.1080/00221686.2010.9641251 [33] 伊姆森，M。；科内利斯，J。；索伦塔勒，B。；Horvath，C。；Teschner，M.，隐式不可压缩SPH，IEEE Trans。视觉。计算。图表。，20, 3, 426-435 (2014) ·doi:10.1109/TVCG.2013.105 [34] 伊斯兰，M.R.I。；Chakraborty，S。；Shaw，A.，《关于光滑粒子流体动力学中的一致性和能量守恒》，国际分子杂志。方法工程，116，9，601-632（2018）·doi:10.1002/nme.5937 [35] 考夫曼，D.M。；Pai，D.K.，不等式约束非光滑哈密顿系统的几何-数值积分，SIAM J.Sci。计算。，34、5、A2670-A2703（2012）·Zbl 1259.65202号 ·数字对象标识代码：10.1137/100800105 [36] 考夫曼，D.M。；铃田，S。；James，D.L。；Pai，D.K.，多体系统中摩擦接触的交错投影，ACM Trans。图表。，27, 5, 164 (2008) ·doi:10.1145/1409060.1409117 [37] 克里米，A。；Rezoug，M。；Khelladi，S。；Nogueira，X。；优雅，M。；Ramírez，L.，《光滑粒子流体动力学：界面多相流体流动模拟的一致模型》，J.Compute。物理。，358, 53-87 (2018) ·Zbl 1381.76287号 ·doi:10.1016/j.jcp.2017.12.006 [38] Lastiwka，M。；巴萨，M。；新泽西州昆兰，《SPH中的可渗透和非反射边界条件》，国际期刊数字。《液体方法》，61，7，709-724（2009）·Zbl 1421.76182号 ·doi:10.1002/fld.1971年 [39] 利伯斯基，L。；Petschek，A。；Carney，T。；希普，J。；Allahdadi，F.，《高应变拉格朗日流体动力学：动态材料响应的三维SPH代码》，J.Compute。物理。，109, 1, 67-75 (1993) ·Zbl 0791.76065号 ·doi:10.1006/jcph.1993.1199 [40] 刘，G。；Liu，M.B.，《平滑粒子流体动力学：无网格粒子方法》（2003），新加坡：世界科学出版社，新加坡·Zbl 1046.76001号 ·doi:10.1142/5340 [41] Lucy，L.B.，《裂变假说测试的数值方法》，Astron。J.，8211013-1024（1977年）·doi:10.1086/112164 [42] J.C.马丁。；Moyce，W.J.，第四部分：刚性水平面上液柱坍塌的实验研究，Philos。事务处理。R.Soc.A，数学。物理学。工程科学。，244, 882, 312-324 (1952) ·doi:10.1098/rsta.1952.0006 [43] Mazhar，H。；Heyn，T。；塔索拉。；Negrut，D.，使用Nesterov方法加速摩擦和接触的多体动力学，ACM Trans。图表。，34, 3 (2015) ·Zbl 1333.68258号 ·doi:10.145/2735627 [44] Mazhar，H。；Pazouki，A。；拉赫什，M。；贾库马尔，P。；Negrut，D.，《通过平滑粒子流体动力学处理流体-固体相互作用问题的微分变分方法》，J.Compute。物理。，371, 92-119 (2018) ·Zbl 1415.70014号 ·doi:10.1016/j.jp.2018.05.013 [45] Monaghan，J.J.，SPH简介，计算。物理学。社区。，48, 1, 89-96 (1988) ·Zbl 0673.76089号 ·doi:10.1016/0010-4655（88）90026-4 [46] Monaghan，J.J.，《用SPH模拟自由表面流动》，J.Compute。物理。，110, 2, 399-406 (1994) ·Zbl 0794.76073号 ·doi:10.1006/jcph.1994.1034 [47] Monaghan，J.J.，《平滑粒子流体动力学》，众议员程序。物理。，68, 1, 1703-1759 (2005) ·Zbl 1160.76399号 ·doi:10.1088/0034-4885/68/8/R01 [48] 莫纳汉，J.J。；Kocharyan，A.，多相流的SPH模拟，计算。物理学。社区。，87, 1, 225-235 (1995) ·Zbl 0923.76195号 ·doi:10.1016/0010-4655（94）00174-Z [49] 莫里斯，J.P。；福克斯·P·J。；Zhu，Y.，使用SPH对低雷诺数不可压缩流动建模，J.Comput。物理。，136, 1, 214-226 (1997) ·Zbl 0889.76066号 ·doi:10.1006/jcph.1997.5776 [50] 奈尔，P。；Pöschel，T.，使用不可压缩SPH的动态毛细管现象，化学。工程科学。，176, 192-204 (2018) ·doi:10.1016/j.ces.2017.10.042 [51] 内格鲁特，D。；塞尔维亚人，R。；Tasora，A.，将摩擦和接触的多体动力学视为微分互补问题，J.Compute。非线性动力学。，13, 1 (2017) ·数字对象标识代码：10.1115/1.4037415 [52] NVIDIA:CUDA编程指南（2015）。可在线访问http://docs.nvidia.com/cuda/cuda-c-programming-guide/index.html [53] OpenMP：规范标准4.5（2017）。可在线访问网址：http://openmp.org/wp/ [54] Pazouki，A。；Negrut，D.，《颗粒特性对管流中悬浮物径向分布影响的数值研究》，计算。流体，108，1-12（2015）·Zbl 1390.76895号 ·doi:10.1016/j.compfluid.2014.11.027 [55] Pazouki，A。；Kwarta，M。；威廉姆斯，K。；利科斯，W。；塞尔维亚人，R。；贾库马尔，P。；Negrut，D.，《柔顺接触与刚性接触：颗粒动力学背景下的比较》，Phys。E版，96（2017）·doi:10.1103/PhysRevE.96.042905 [56] 项目Chrono:Chrono：一个基于物理的动态系统模拟的开源框架。http://projectchrono.org。访问时间：2016-03-07 [57] 项目Chrono开发团队：Chrono：一个基于物理的动态系统模拟的开源框架。https://github.com/projectchrono/chron。访问时间：2019-12-07 [58] 新泽西州昆兰。；巴萨，M。；Lastiwka，M.，《无网格粒子方法中的截断误差》，国际期刊《数值》。方法工程，66，13，2064-2085（2006）·Zbl 1110.76325号 ·数字对象标识代码：10.1002/nme.1617 [59] 拉赫什，M。；Pazouki，A。；塞尔维亚人，R。；Negrut，D.，使用半隐式积分方案求解基于SPH的流体-固体相互作用问题，计算。方法应用。机械。工程师，345100-122（2019）·Zbl 1440.74128号 ·doi:10.1016/j.cma.2018.09.027 [60] Rakhsha，M.，Yang，Y.，Kees，C.，Negrut，D.：拉格朗日与欧拉：自由表面流动和流固相互作用问题两种求解方法的分析。审查中（2019年） [61] Randles，P.W。；Libersky，L.D.，《平滑粒子流体动力学：一些最近的改进和应用》，计算。方法应用。机械。工程，139，1-4，375-408（1996）·兹伯利0896.73075 ·doi:10.1016/S0045-7825（96）01090-0 [62] Saad，Y.：Krylov子空间方法概述及其在控制问题中的应用。美国宇航局艾姆斯研究中心高级计算机科学研究所（1989年）·Zbl 0722.93028号 [63] 沙德鲁，M.S。；Zainali，A。；Yildiz，M。；Suleman，A.，一种稳健的弱可压缩SPH方法及其与不可压缩SPH-的比较，国际J·数值。方法工程，89，8，939-956（2012）·Zbl 1242.76279号 ·doi:10.1002/nme.3267 [64] 邵，S。；Lo，E.，《模拟自由表面牛顿和非牛顿流动的不可压缩SPH方法》，高级水资源。，26, 7, 787-800 (2003) ·doi:10.1016/S0309-1708（03）00030-7 [65] B.史密斯。；考夫曼，D.M。；Vouga，E。；Tamstorf，R。；Grinspin，E.，《关于同时冲击的思考》，ACM Trans。图表。，31, 4 (2012) ·doi:10.1145/2185520.2185602 [66] Stewart，D.E.，摩擦和冲击的刚体动力学，SIAM Rev.，42，1，3-39（2000）·Zbl 0962.70010号 ·doi:10.1137/S0036144599360110 [67] Stewart，D.E。；Trinkle，J.C.，《具有非弹性碰撞和库仑摩擦的刚体动力学的隐式时间步长方案》，国际期刊编号。方法工程，39，2673-2691（1996）·Zbl 0882.70003号 ·doi:10.1002/（SICI）1097-0207（19960815）39:15<2673:：AID-NME972>3.0.CO；2-I型 [68] 武田，H。；Miyama，S.M。；Sekiya，M.，光滑粒子流体动力学粘性流动的数值模拟，Prog。西奥。物理。，92, 5, 939-960 (1994) ·doi:10.1143/ptp/92.5.939 [69] 塔塔科夫斯基，A.M。；Panchenko，A.，《多相流的对偶力平滑粒子流体动力学模型：表面张力和接触线动力学》，J.Compute。物理。，305, 1119-1146 (2016) ·兹比尔1349.76739 ·doi:10.1016/j.jcp.2015.08.037 [70] 塔索拉。；Anitescu，M.，《模拟大颗粒流的凸互补方法》，J.Compute。非线性动力学。，5, 3, 1-10 (2010) ·Zbl 1200.90160号 ·数字对象标识代码：10.1115/1.4001371 [71] 塔索拉。；塞尔维亚人，R。；Mazhar，H。；Pazouki，A。；梅兰兹，D。；Fleischmann，J。；泰勒，M。；杉山，H。；内格鲁特，D。；Kozubek，T.，Chrono:一个开源多物理动力学引擎，《科学与工程中的高性能计算》，19-49（2016），柏林：施普林格出版社，柏林·Zbl 1382.70001号 ·doi:10.1007/978-3-319-40361-82 [72] Trask，N.，Kim，K.，Tartakovsky，A.，Perego，M.，Parks，M.L.：用于模拟几何复杂有界域中单相和多相流的高标度隐式SPH代码。技术代表，桑迪亚国家实验室（SNL-NM），新墨西哥州阿尔伯克基（美国）（2015年） [73] 特拉斯克，N。；马克西，M。；金姆，K。；佩雷戈，M。；Parks，M.L。；Yang，K。；Xu，J.，非定常低雷诺数流动的可缩放一致二阶SPH解算器，计算。方法应用。机械。工程，289，155-178（2015）·Zbl 1423.76372号 ·doi:10.1016/j.cma.2014.12.027 [74] Turek，S。；Hron，J.，弹性物体和层流不可压缩流之间流体-结构相互作用的数值基准建议，流体-结构交互作用，371-385（2006），柏林：施普林格，柏林·Zbl 1323.76049号 ·doi:10.1007/3-540-34596-5_15 [75] Wang，Z.B。；陈，R。；Wang，H。；廖琦（Liao，Q.）。；朱，X。；Li，S.Z.，用于模拟多相流的光滑粒子流体动力学概述，应用。数学。型号。，40, 23, 9625-9655 (2016) ·Zbl 1443.76084号 ·doi:10.1016/j.apm.2016.06.030 [76] Xu，X。；Deng，X.L.，一种改进的弱可压缩SPH方法，用于模拟粘性和粘弹性流体的自由表面流动，计算。物理学。社区。，201, 43-62 (2016) ·Zbl 1348.76119号 ·doi:10.1016/j.cpc.2015.12.016 [77] Xu，R。；斯坦斯比，P。；Laurence，D.，基于投影方法和新方法J.Compute的不可压缩SPH（ISPH）的准确性和稳定性。物理。，228, 18, 6703-6725 (2009) ·Zbl 1261.76047号 ·doi:10.1016/j.jcp.2009.05.032 [78] 杨琼。；琼斯，V。；McCue，L.，使用SPH-FEM模型的自由表面流与可变形结构物的相互作用，海洋工程，55，136-147（2012）·doi:10.1016/j.oceaneng.2012.06.031 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。