廉、严建;Bui，Ha H。;阮，Giang D。;Tran、Hieu T。;哈克，阿萨杜尔 考虑各向异性扩散的非饱和多孔介质瞬态渗流的通用SPH框架。 （英语） Zbl 1507.76164号 计算。方法应用。机械。工程师。 387，文章ID 114169，43 p.（2021）. 摘要：光滑粒子流体力学（SPH）方法已成为解决许多具有挑战性的工程应用的流行数值工具。然而，它在求解各向异性瞬态非饱和饱和渗流问题中的应用在文献中尚属空白。在解决这一重要问题所涉及的几个挑战中，缺乏一个稳健的通用SPH公式来获得无序粒子系统二阶导数的精确近似值是阻碍这一情况发生的关键障碍。在本文中，我们首先提出了一个通用的二阶导数SPH公式，该公式能够考虑各向异性扩散，并证明了新公式在二阶导数方面优于现有的SPH公式并且在高度无序的粒子系统上实现了高精度。随后，我们利用这一新开发的SPH公式首次处理了多孔介质中的瞬态渗流，考虑了各向异性流动和从非饱和状态到饱和状态的完全时间依赖转换以及副变量。这自然导致了一个通用的SPH框架，用于使用一组拉格朗日粒子（或拉格朗夫离散化）解决非饱和/饱和多孔介质中的瞬态渗流问题。还提出了用于解决SPH中一般渗流问题的稳健边界条件及其实现，使所提出的SPH模型能够自动捕获渗流面（包括无侧限渗流的自由面），而无需任何困难或特殊处理。此外，首次提出了一种利用Voronoi细分技术生成适用于大规模问题SPH模拟的随机粒子系统的简单方法，用于研究瞬态渗流问题。通过几个验证实例，证明了所提出的SPH框架在解决复杂渗流问题（包括多分辨率粒子配置问题）方面的能力。结果表明，新提出的SPH模型可以应用于求解复杂的瞬态渗流问题。这标志着SPH未来发展的一个重要里程碑，以解决多孔介质大变形和破坏的完全耦合流动变形问题，这对于最先进的基于网格的数值方法来说是一个挑战。 引用于6文件 MSC公司： 76米28 粒子法和晶格气体法 76兰特 扩散 76S05号 多孔介质中的流动；过滤；渗流 关键词：SPH公司;瞬态渗流;各向异性扩散;非饱和介质;Voronoi细分 软件：PolyMesher公司 PDF格式BibTeX公司 XML格式引用 \textit{Y.Lian}等人，计算。方法应用。机械。工程387，文章ID 114169，43 p.（2021；Zbl 1507.76164） 全文： 内政部 参考文献： [1] Pedroso，D.M.，非饱和多孔介质中瞬态渗流的解决方案，计算。方法应用。机械。工程，285791-816（2015）·Zbl 1423.76437号 [2] 郑浩。；Liu，D.F。；Lee，C.F。；Tham，L.G.，自由表面渗流问题的Signorini型新公式，国际。J.数字。方法工程，64，1-16（2005）·Zbl 1073.76069号 [3] 郑浩。；Dai，H.C。；Liu，D.F.，多孔介质无侧限渗流问题的变分不等式公式，应用。数学。型号。，33, 437-450 (2009) ·Zbl 1167.76370号 [4] 莱西，S.J。；Prevost，J.H.，《穿过多孔介质的流动：定位自由表面的程序》，国际期刊编号。分析。方法地质力学。，11, 585-601 (1987) ·兹比尔062373115 [5] 博尔贾，R.I。；Kishnani，S.S.，关于用牛顿法求解椭圆自由边界问题，计算。方法应用。机械。工程，88，341-361（1991）·Zbl 0745.76066号 [6] Oden，J.T。；Kikuchi，N.，变分不等式理论及其在多孔介质流动问题中的应用，国际。工程科学杂志。，18, 1173-1284 (1980) ·Zbl 0444.76069号 [7] Lee，K.K。；Leap，D.I.，使用无边界坐标系方法模拟自由表面和渗流面，J.Hydrol。，196, 297-309 (1997) [8] Bardet，J.P。；Tobita，T.，解决自由表面渗流问题的实用方法，计算。岩土工程。，29, 451-475 (2002) [9] Brezis，H.，《关于生态保护问题的新提法》，C.R.Acad。科学。巴黎，287，711-714（1978），https://ci.nii.ac.jp/naid/80013858036（2021年5月17日访问）·Zbl 0391.76072号 [10] 马里兰州埃雷罗斯。；Mabssout，M。；Pastor，M.，水平集方法在多孔介质流动中移动界面和自由表面问题的应用，计算。方法应用。机械。工程，195，1-25（2006）·Zbl 1091.76054号 [11] Bui，H.H。；Nguyen，G.D.，通过可变形多孔介质模拟流动的流固耦合SPH方法，国际固体结构杂志。，125244-264（2017） [12] L.B.Lucy，裂变假设检验的数值方法，82（1977）1013-1024。 [13] A.R.Gingold，J.J.Monaghan，《平滑粒子流体动力学：非球形恒星的理论和应用》，181（1977）375-389·Zbl 0421.76032号 [14] 哈迪，S。；施莱纳，M。；Janoske，U.，《使用有限粒子方法增强小尺度浅水方程的平滑粒子流体动力学》，计算。方法应用。机械。工程，344，360-375（2019）·Zbl 1440.76109号 [15] 刘凯。；王，Y。；黄，M。；袁伟华，用随机广义插值材料点法进行边坡破坏后分析，国际地质力学杂志。，21，第04021015条pp.（2021） [16] 胡，W。；潘·W。；拉赫什，M。；田强。；胡，H。；Negrut，D.，一致多分辨率平滑粒子流体动力学方法，计算机。方法应用。机械。工程，324，278-299（2017）·Zbl 1439.76128号 [17] Tafuni，A。；Domínguez，J.M。；瓦康迪奥，R。；Crespo，A.J.C.，平滑粒子流体动力学GPU模型中开放边界条件处理的通用算法，计算。方法应用。机械。工程，342，604-624（2018）·Zbl 1440.76125号 [18] Khayyer，A。；H.Gotoh。；Shimizu，Y.，两种粒子正则化方案的精度和守恒性的比较研究，以及ISPH背景下优化粒子移位方案的建议，J.Compute。物理。，332, 236-256 (2017) ·Zbl 1378.76094号 [19] 森川，D。；Senadheera，H。；Asai，M.，用于大规模模拟的多GPU环境中的显式不可压缩平滑粒子流体动力学，计算。部分。机械。，8, 493-510 (2021) [20] 季，Z。；Fu，L。；胡晓云。；Adams，N.A.，SPH的一种新的多分辨率并行框架，计算。方法应用。机械。工程师，346，1156-1178（2019）·Zbl 1440.76112号 [21] 胡晓云。；Adams，N.A.，《不可压缩多相SPH方法》，J.Compute。物理。，227, 264-278 (2007) ·Zbl 1126.76045号 [22] 森川，D.S。；Asai，M.，耦合整体拉格朗日SPH-EISPH，用于与大型变形超弹性固体的流体-结构相互作用，计算。方法应用。机械。工程，381，第113832条pp.（2021）·Zbl 1506.76129号 [23] Monteleone，A。；De Marchis，M。；米利奇，B。；Napoli，E.，高度复杂流动的平滑粒子流体动力学模拟的多域方法，计算。方法应用。机械。工程师，340956-977（2018）·Zbl 1440.76118号 [24] 拉赫什，M。；Pazouki，A。；塞尔维亚人，R。；Negrut，D.，使用半隐式积分方案求解基于SPH的流体-固体相互作用问题，计算。方法应用。机械。工程，345100-122（2019）·Zbl 1440.74128号 [25] Antoci，C。；加拉蒂，M。；Sibilla，S.，SPH流体-结构相互作用的数值模拟，计算。结构。，85, 879-890 (2007) [26] Khayyer，A。；Y.清水。；H.Gotoh。；Nagashima，K.，用于复合结构水弹性FSI的耦合不可压缩SPH-Hamilton SPH解算器，应用。数学。型号。，94242-271（2021）·Zbl 1481.74155号 [27] Khayyer，A。；Y.清水。；H.Gotoh。；Hattori，S.，《用于准确高效模拟海洋工程中水弹性流体-结构相互作用的多分辨率ISPH-SPH》，海洋工程，226，第108652页，（2021） [28] 张，Z.L。；Long，T。；Chang，J.Z。；Liu，M.B.，用于模拟具有大流体变形的流体-结构相互作用问题的平滑粒子元方法（SPEM），计算。方法应用。机械。工程，356261-293（2019）·Zbl 1441.76097号 [29] Asai，M。；李毅。；钱德拉，B。；Takase，S.，使用耦合稳定ISPH-DEM与多体接触的能量跟踪脉冲方法相结合的流体-刚体相互作用模拟和验证，计算。方法应用。机械。工程，377，第113681条pp.（2021）·Zbl 1506.74103号 [30] 张，C。；Rezavand，M。；Hu，X.，流体-结构相互作用的多分辨率SPH方法，J.Compute。物理。，429 (2021) ·Zbl 07500760号 [31] 巴，K。；Gakwaya，A.，大变形问题中固体力学的热机械全拉格朗日SPH公式，计算。方法应用。机械。工程，342458-473（2018）·Zbl 1440.74011号 [32] Bui，H.H。；福川，R。；Sako，K。；Ohno，S.，拉格朗日无网格颗粒法（SPH），用于使用弹塑性土本构模型的岩土材料的大变形和破坏流动，国际J·数值。分析。方法地质力学。，32, 1537-1570 (2008) ·Zbl 1273.74563号 [33] Nguyen，C.T。；Nguyen，C.T。；Bui，H.H。；Nguyen，G.D。；Fukagawa，R.，《使用粘性阻尼和应力调节模拟颗粒流的基于SPH的新方法》，《滑坡》，第14期，第69-81页（2017年） [34] Chalk，C.M。；Pastor，M。；Peakall，J。；博曼，D.J。；雪橇，P.A。；墨菲，W。；Fuentes，R.，《应力-颗粒光滑颗粒流体动力学：边坡破坏和破坏后行为的应用》，计算。方法应用。机械。工程，366（2020）·Zbl 1442.74124号 [35] R.I.Borja，采用乘法弹塑性的干颗粒流连续流体动力学，3（2018）。http://dx.doi.org/10.1007/s11440-018-0700-3。 [36] 杨，E。；Bui，H.H。；De Sterck，H。；Nguyen，G.D。；Bouazza，A.，地球物理颗粒流预测的可扩展并行计算SPH框架，计算。岩土工程。，121，第103474条pp.（2020） [37] Fávero Neto，A.H。；Askarinejad，A。；斯普林曼，S.M。；Borja，R.I.，使用更新的拉格朗日连续体粒子法模拟仪器化测试边坡上的泥石流，《岩土学报》。，15, 2757-2777 (2020) [38] 黄，C。；Liu，M.B.，用混合光滑颗粒流体动力学模拟含水沉积物，计算。机械。，66, 877-891 (2020) ·Zbl 1462.76138号 [39] Bui，H.H。；福川，R。；Sako，K。；Wells，J.C.，《利用弹塑性光滑粒子流体动力学（SPH）进行边坡稳定性分析和不连续边坡破坏模拟》，《岩土工程》，61，565-574（2011） [40] 何，X。；Liang，D。；Bolton，M.D.，《切坡滑坡径流：无网格模拟》，《岩土工程》，68，50-63（2018） [41] Bui，H.H。；Fukagawa，R.，《饱和土壤的改进SPH方法及其在研究路堤破坏机理中的应用：静水孔压案例》，国际J·数值。分析。方法地质力学。，37, 31-50 (2013) [42] Pastor，M。；布兰科·T。；哈达德，B。；Petrone，S.公司。；桑切斯·莫尔斯，M。；Drempetic，V。；Issler，D。；克罗斯塔，G.B。；卡西尼，L。；索比诺，G。；Cuomo，S.，《SPH深度积分模型在滑坡径流分析中的应用》，《滑坡》，第11793-812页（2014年） [43] 布兰科·T。；Pastor，M.，《用于动力学大变形问题的稳定化runge-kutta，taylor平滑粒子流体动力学算法》，国际。J.数字。方法工程，91，1427-1458（2012） [44] 赵，S。；Bui，H.H。；莱米亚莱，V。；Nguyen，G.D。；Darve，F.，《SPH中柔性受限边界条件建模的通用方法及其应用》，国际J·数值。分析。方法地质力学。，43, 1005-1031 (2019) [45] Chen，W。；Qiu，T.，使用SPH方法模拟地震引起的边坡变形，Int.J.Numer。分析。方法地质力学。，38, 297-330 (2014) [46] 杨，E。；Bui，H.H。；Nguyen，G.D。；Choi，C.E。；Ng、C.W.W。；De Sterck，H。；Bouazza，A.，《颗粒流对刚性控制结构冲击机理的数值研究》，《岩土学报》。(2021) [47] 詹，L。；彭，C。；张，B。；Wu，W.，用混合接触法研究土壤和刚体系统之间动态相互作用的SPH框架，国际J·数值。分析。方法地质力学。，44, 1446-1471 (2020) [48] 谢赫，B。；邱，T。；Ahmadipur，A.，《模拟摩擦土-结构相互作用中SPH边界方法的比较》，《岩土学报》。，1-20 (2020) [49] Bui，H.H。；Kodikara，J.K。；Bouazza，A。；哈克，A。；Ranjith，P.G.，节段式挡土墙系统大变形和破坏后分析的新计算方法，国际J·数值。分析。方法地质力学。，38, 1321-1340 (2014) [50] Tran，H.T。；王，Y。；Nguyen，G.D。；Kodikara，J。；桑切斯，M。；Bui，H.H.，使用尺寸相关的SPH计算方法模拟粘性土中的3D干燥开裂，计算。岩土工程。，116，第103209条pp.（2019） [51] Tran，H.T。；Nguyen，N.H.T。；Nguyen，G.D。；Bui，H.H.，《粘性土中收缩诱发裂缝的无网格SPH模型》，Lect。注释Civ。工程，54，889-894（2020） [52] 王，Y。；Tran，H.T。；Nguyen，G.D。；Ranjith，P.G。；Bui，H.H.，《使用SPH颗粒和嵌入的断裂过程区模拟混合模式断裂》，国际期刊编号。分析。方法地质力学。，44, 1417-1445 (2020) [53] 王，Y。；Bui，H.H。；Nguyen，G.D。；Ranjith，P.G.，《一种新的基于SPH的连续体框架，带有嵌入的断裂过程区，用于模拟岩石断裂，国际固体结构杂志》。，159, 40-57 (2019) [54] 加雷达什，S。；沈，L。；Gan，Y.，爆破破碎岩石力学和水力特性的数值研究，工程计算。，36, 915-929 (2020) [55] 布兰科·T。；Pastor，M.，地质力学耦合问题的稳定分步Runge-Kutta-Taylor SPH算法，计算。方法应用。机械。工程，221-222，41-53（2012）·Zbl 1253.74099号 [56] 科尔马尼，E。；邱，T.，使用平滑粒子流体动力学方法模拟固定多孔介质的渗流，699-708（2017），美国土木工程师学会（ASCE） [57] Bui，H.H。；Sako，K。；Fukagawa，R.，使用平滑粒子流体动力学（SPH）方法对土壤-水相互作用进行数值模拟，J.Terramech。，44333-346（2007年） [58] Brookshaw，L.，粒子模拟中计算辐射热扩散的方法，Astron。澳大利亚州。，6, 207-210 (1985) [59] 克利里，P.W。；Monaghan，J.J.，《使用平滑粒子流体动力学的传导建模》，J.Compute。物理。，148, 227-264 (1999) ·Zbl 0930.76069号 [60] Monaghan，J.J.，《平滑粒子流体动力学》，众议员程序。物理。，68, 1703-1759 (2005) [61] 埃斯帕诺，P。；Revenga，M.，平滑耗散粒子动力学，物理学。修订版E，67，12（2003） [62] 佩特科娃，M。；Springel，V.，宇宙模拟代码小工具中辐射传输的实现，Mon。不是。R.阿斯顿。Soc.，3961383-1403（2009年） [63] 霍普金斯，P.F.，《无网格数值磁流体动力学中的各向异性扩散》，孟买。不是。R.阿斯顿。Soc.，4663387-3405（2017年） [64] Chen，J.K。；Beraun，J.E。；Carney，T.C.，《热传导边值问题的修正光滑粒子法》，国际。J.数字。方法工程，46，231-252（1999）·Zbl 0941.65104号 [65] Bui，H.H。；Nguyen，G.D.，平滑粒子流体力学（SPH）及其在地质力学中的应用：从固体断裂到颗粒行为和多孔介质中的多相流，计算。岩土工程。，138，第104315条pp.（2021） [66] Fatehi，R。；Manzari，M.T.，平滑粒子流体动力学中的误差估计和二阶导数的新方案，计算。数学。申请。，61, 482-498 (2011) ·Zbl 1211.76089号 [67] 迪尔，S。；洛克菲勒，G。；傅兰雅，C.L。；Riethmiller，D。；Statler，T.S.，为平滑粒子流体动力学模拟生成最佳初始条件，Publ。阿童木。澳大利亚州。，32 (2015) [68] 他，S。；Springel，V.，《粒子流体动力学与细分技术》，周一。不是。R.阿斯顿。《社会学杂志》，4062289-2311（2010） [69] Ghaffari，医学硕士。；Xiao，S.，具有应力点和质心voronoi镶嵌（CVT）拓扑优化的光滑粒子流体动力学，Int.J.Comput。方法，13，1-23（2016）·Zbl 1359.74482号 [70] Shobeyri，G。；Ardakani，R.R.，使用voronoi图提高SPH方法的准确性，伊朗。科学杂志。Technol.-技术事务处理。公民。工程，41，345-350（2017） [71] 卞，X。；李，Z。；Karniadakis，G.E.，通过区域分解通过平滑粒子流体动力学进行的多分辨率流动模拟，J.Compute。物理。，297, 132-155 (2015) ·Zbl 1349.76663号 [72] 奥卡，F。；Shahbodagh，B。；Kimoto，S.，《非饱和弹塑性土壤中动态应变局部化的计算模型》，国际J·数值。分析。方法地质力学。，43, 138-165 (2019) [73] 奥卡，F。；Kodaka，T。；铃木，H。；Kim，Y.S。；北西松。；Kimoto，S.，《三轴压缩下非饱和压实粉土特性的试验研究》，《发现的土壤》。，50, 27-44 (2010) [74] 贝尔·J。；Cheng，A.H.-D.，《地下水流动和污染物运移建模》（2010年）·Zbl 1195.76002号 [75] van Genuchten，M.T.，预测非饱和土壤导水率的封闭式方程，土壤科学。《美国社会杂志》，44，892-898（1980） [76] Monaghan，J.J.，《插值的B样条外推》，J.Compute。物理。，60, 253-262 (1985) ·Zbl 0588.41005号 [77] 奥格，G。；多林，M。；亚历山德里尼，B。；Ferrant，P.，《一种改进的SPH方法：走向高阶收敛》，J.Compute。物理。，225, 1472-1492 (2007) ·Zbl 1118.76050号 [78] Randles，P.W。；Libersky，L.D.，《光滑粒子流体动力学：一些最近的改进和应用》，计算。方法应用。机械。工程，139，375-408（1996）·Zbl 0896.73075号 [79] Chen，J.K。；Beraun，J.E.，非线性动力学问题的广义光滑粒子流体动力学方法，计算。方法应用。机械。工程，190225-239（2000）·Zbl 0967.76077号 [80] 杜琪。；费伯，V。；Gunzburg，M.，《质心Voronoi镶嵌：应用与算法》，SIAM Rev.，41637-676（1999）·Zbl 0983.65021号 [81] Talischi，C。；Paulino，G.H。；佩雷拉，A。；Menezes，I.F.M.，PolyMesher:用matlab编写的多边形元素通用网格生成器，Struct。多磁盘。最佳。，45, 309-328 (2012) ·Zbl 1274.74401号 [82] Lloyd，S.P.，PCM中的最小二乘量化，IEEE Trans。Inf.理论，28，129-137（1982）·Zbl 0504.94015号 [83] 刘，M.B。；Liu，G.R.，《光滑粒子流体动力学（SPH）：概述和最新发展》，Arch。计算。方法工程，17，25-76（2010）·Zbl 1348.76117号 [84] Bui，H.H.，《使用弹塑性本构模型研究土工材料大变形和断裂后的拉格朗日无网格颗粒法》（SPH），立命馆大学。，http://r-cube.ritsumei.ac.jp/repo/repository/rcube/9152/k_466_e.htm [85] Gardner，W.R.，《应用于地下水位蒸发的非饱和水分流动方程的一些稳态解》，土壤科学。，85, 228-232 (1958) [86] Polubarinova-Kochina，P.J。；De Wiest，R.J.M.，《地下水运动理论》（1962年），普林塞特。大学出版社，https://lib.ugent.be/catalog/rug01:000200814（2021年5月16日访问）·Zbl 0114.42601号 此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配，并且可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。