×

弱可压缩光滑粒子流体力学的拉格朗日自由流边界条件。 (英语) Zbl 07715233号

小结:本文提出了弱可压缩光滑粒子流体力学(WCSPH)的第一个拉格朗日自由流边界条件,显著降低了模拟静态或运动几何体周围流动的区域尺寸效应,提高了计算效率。边界条件是基于以下几种数值技术实现的。首先,时空识别和确保技术能够准确识别表面粒子,而不会产生误判。其次,密度和速度的远场修正精确地针对表面粒子。第三,改进的入口和出口边界保证了平滑的流入剖面和自由流出颗粒,而不存在截断核的问题。通过平板绕流和圆柱绕流的算例验证了本方法的准确性和计算效率。为了证明其适用性,该方法还用更复杂的例子进行了测试,例如流体-结构相互作用(FSI)问题和球体周围的三维流动。

MSC公司:

7.6亿 流体力学基本方法
6500万 偏微分方程、初值和含时初边值问题的数值方法
76天xx 不可压缩粘性流体

软件:

SPHinXsys公司
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Lucy,L.,裂变假说测试的数值方法,Astron。J.,82,1013-1024(1977)
[2] Gingold,R。;Monaghan,J.,光滑粒子流体动力学:理论及其在非球面恒星中的应用,Mon。不是。R.阿斯顿。《社会学杂志》,181,375-389(1977)·Zbl 0421.76032号
[3] Monaghan,J.,《用SPH模拟自由表面流动》,J.Compute。物理。,110, 399-406 (1994) ·Zbl 0794.76073号
[4] 莫里斯,J。;福克斯,P。;Zhu,Y.,使用SPH建模低雷诺数不可压缩流,J.Compute。物理。,136, 214-226 (1997) ·Zbl 0889.76066号
[5] 奔驰,W。;Asphaug,E.,使用光滑粒子流体动力学模拟脆性固体,计算。物理学。社区。,87, 253-265 (1995) ·Zbl 0918.73335号
[6] 格雷,J。;莫纳汉,J。;Swift,R.,SPH弹性动力学,计算。方法应用。机械。工程,1906641-6662(2001)·Zbl 1021.74050号
[7] Antoci,C。;加拉蒂,M。;Sibilla,S.,SPH流体-结构相互作用的数值模拟,计算。结构。,85, 11-14, 879-890 (2007)
[8] Rafiee,A。;Thiagarajan,K.,用于模拟流体-低塑性结构相互作用的SPH投影方法,计算。方法应用。机械。Eng.,198,文章110028 pp.(2009)·Zbl 1228.76117号
[9] 费德里科,I。;马龙,S。;Colagrossi,A。;亚里士多德·F。;Antuono,M.,《通过SPH模型模拟2D明渠水流》,Eur.J.Mech。B、 流体,34,35-46(2012)·Zbl 1258.76123号
[10] Tan,S。;Cheng,N。;谢毅。;Shao,S.,光滑床上明渠流动的不可压缩SPH模拟,J.Hydro-environ。Res.,9340-353(2015年)
[11] Staroszczyk,R.,用于模拟剧烈自由表面流体流动的不可压缩SPH模型,Arch。水力环境。机械。,61, 61-83 (2014)
[12] 马龙,S。;安托诺,M。;Colagrossi,A。;科利奇奥,G。;勒图泽,D。;Graziani,G.,用于模拟剧烈冲击流的Delta-SPH模型,计算。方法应用。机械。工程,200,1526-1542(2011)·Zbl 1228.76116号
[13] 陈,R。;邵,S。;刘,X。;周,X.,浅水SPH模型在山区河流中的应用,J.Appl。流体力学。,8, 863-870 (2015)
[14] 瓦康迪奥,R。;罗杰斯,B。;斯坦斯比,P。;Mignosa,P.,《用于实际洪水模拟的开放边界浅水流SPH建模》,J.Hydraul。工程,138530-541(2012)
[15] 马蒂厄,D。;大卫,L。;Bertrand,A.,《浅水海岸流的SPH建模》,J.Hydraul。决议,48,118-125(2010)
[16] 维奥罗,D。;罗杰斯,B.,《自由表面流动的平滑粒子流体动力学(SPH):过去、现在和未来》,J.Hydraul。决议,54,1-26(2016)
[17] 张,C。;Zhu,Y。;Wu,D。;N.亚当斯。;胡晓霞,《光滑粒子流体力学:方法学发展和最新成果》,J.Hydrodyn。,34, 767-805 (2022)
[18] Luo,M。;Khayyer,A。;Lin,P.,《海洋与海岸工程中的粒子方法》,应用。Ocean Res.,114,第102734条,pp.(2021)
[19] 沙德鲁,M。;奥格,G。;Touze,D.L.,流体流动的平滑粒子流体动力学方法,面向工业应用:动机、现状和挑战,计算。流体,136,11-34(2016)·Zbl 1390.76764号
[20] Tafuni,A。;Domínguez,J。;瓦康迪奥,R。;Crespo,A.,平滑粒子流体动力学GPU模型中开放边界条件处理的通用算法,计算。方法应用。机械。工程,342,604-624(2018)·Zbl 1440.76125号
[21] 张,S。;胡,X。;Mi,C.,基于SPH方法的高速列车气动制动盘制动力研究,Phys。气体,5(0-0 2020)
[22] 贝利斯,A。;Turkel,E.,可压缩流动的远场边界条件,J.Compute。物理。,48, 2, 182-199 (1982) ·Zbl 0494.76072号
[23] 马丁。;巴萨,M。;Quinlan,N.,《SPH中的渗透和非反射边界条件》,国际J·数值。方法流体,61,709-724(2009)·Zbl 1421.76182号
[24] Molteni,D。;Grammauta,R。;Vitanza,E.,《利用平滑粒子流体动力学进行浅波模拟的简单吸收层条件》,海洋工程,62,78-90(2013)
[25] S.Braun,L.Wieth,R.Koch,H.Bauer,《SPH中渗透边界条件的框架:入口、出口、周期性》,2015年。
[26] 瓦康迪奥,R。;罗杰斯,B。;斯坦斯比,P。;Mignosa,P.,《用于实际洪水模拟的开放边界浅水流SPH建模》,J.Hydraul。工程,138530-541(2012)
[27] 费兰德,M。;Joly,A。;骨膜炎,C。;维奥罗,D。;Leroy,A。;莫雷尔,F。;罗杰斯,B.,《使用半分析条件和二维黎曼解算器求解SPH的非定常开放边界》,计算。物理学。社区。,210, 29-44 (2017) ·Zbl 1376.76048号
[28] 马龙,S。;Colagrossi,A。;安托诺,M。;科利奇奥,G。;Graziani,G.,《低雷诺数和中等雷诺数下物体周围粘性流动的精确SPH建模》,J.Compute。物理。,245, 456-475 (2013) ·Zbl 1349.76715号
[29] 王,P。;张,A。;明,F。;Sun,P。;Cheng,H.,光滑粒子流体力学求解的流体动力学新的无反射边界条件,J.流体力学。,860, 81-114 (2019) ·Zbl 1415.76513号
[30] 呼罗珊扎德,南部。;Sousa,J.,《不可压缩SPH的创新开放边界处理》,国际J·数值。《液体方法》,80,161-180(2016)
[31] 阿尔瓦拉多·罗德里格斯,C。;Klapp,J。;西加洛蒂,L。;多明格斯,J。;Sanchez,E.,使用SPH的不可压缩流的无反射出口边界条件,计算。流体,159177-188(2017)·Zbl 1390.76694号
[32] Negi,P。;Ramachandran,P。;Haftu,A.,弱可压缩SPH的改进非反射出口边界条件,计算。方法应用。机械。工程,367,第113119条pp.(2020)·Zbl 1442.76088号
[33] Lee,E。;穆利内克,C。;Xu,R。;维奥罗,D。;劳伦斯,D。;Stansby,P.,SPH网格自由粒子方法的弱可压缩和真正不可压缩算法的比较,J.Comput。物理。,227, 8417-8436 (2008) ·Zbl 1256.76054号
[34] Dilts,G.,移动最小二乘粒子流体动力学-I。一致性和稳定性,国际数值。方法工程,44,1115-1155(1999)·Zbl 0951.76074号
[35] 哈克,A。;Dilts,G.,《粒子方法的三维边界检测》,J.Compute。物理。,226, 1710-1730 (2007) ·Zbl 1173.76393号
[36] Lin,Y。;刘,G。;Wang,G.,一种基于粒子的自由表面检测方法及其在光滑粒子流体动力学(SPH)表面张力效应模拟中的应用,J.Compute。物理。,383, 196-206 (2019) ·Zbl 1451.76093号
[37] 马龙,S。;Colagrossi,A。;勒图泽,D。;Graziani,G.,SPH解算器中的快速自由表面检测和水平集函数定义,J.Comput。物理。,229, 3652-3663 (2010) ·Zbl 1391.76623号
[38] 张,C。;Rezavand先生。;Zhu,Y。;Yu,Y。;Wu,D。;张伟。;Wang,J。;Hu,X.,SPHinXsys:基于平滑粒子流体力学的开放源码多物理和多分辨率库,计算。物理学。社区。,267,第108066条pp.(2021)·Zbl 07695811号
[39] 胡,X。;Adams,N.,《宏观和介观流动的多相SPH方法》,J.Compute。物理。,213, 844-861 (2006) ·Zbl 1136.76419号
[40] 张,C。;Rezavand先生。;胡,X.,弱可压缩光滑粒子流体动力学的双准则时间步长,J.Compute。物理。,404,第109135条pp.(2020)·Zbl 1453.76172号
[41] 张,C。;胡,X。;Adams,N.,平滑粒子流体动力学的广义输运速度公式,J.Compute。物理。,337, 216-232 (2017) ·Zbl 1415.76514号
[42] 张,C。;胡,X。;Adams,N.,基于低耗散黎曼解算器的弱可压缩SPH方法,J.Compute。物理。,335 (2017) ·Zbl 1375.76155号
[43] 张,C。;Rezavand先生。;Hu,X.,流体-结构相互作用的多分辨率SPH方法,J.Compute。物理。,429,第110028条pp.(2021)·Zbl 07500760号
[44] 林德·S。;Xu,R。;斯坦斯比,P。;Rogers,B.,《自由表面流动的不可压缩平滑粒子流体动力学:脉冲流和传播波稳定性和验证的基于扩散的通用算法》,J.Compute。物理。,231, 4, 1499-1523 (2012) ·Zbl 1286.76118号
[45] 奥格,G。;马龙,S。;勒图泽,D。;de Leffe,M.,通过结合准拉格朗日移动传输速度和一致ale公式提高SPH精度,J.Compute。物理。,313, 76-98 (2016) ·Zbl 1349.65537号
[46] 斯基伦,A。;林德·S。;Stansby,P.K。;罗杰斯,B.D.,《时间噪声降低的不可压缩平滑粒子流体力学(SPH)》,广义菲克平滑应用于体水撞击和有效波-体相互作用,计算。方法应用。机械。工程师,265163-173(2013)·Zbl 1286.76120号
[47] Khayyer,A。;H.Gotoh。;Shimizu,Y.,两种粒子正则化方案的精度和守恒性的比较研究,以及ISPH背景下优化粒子移位方案的建议,J.Compute。物理。,332, 236-256 (2017) ·Zbl 1378.76094号
[48] Sun,P。;Colagrossi,A。;马龙,S。;Zhang,A.,delta-plus-SPH模型:进一步改进SPH方案的简单程序,计算。方法应用。机械。工程,31525-49(2017)·Zbl 1439.76133号
[49] Negi,P。;Ramachandran,P.,《无边界二阶收敛弱可压缩光滑粒子流体力学方案的技术》,Phys。流体,34,8,第087125条pp.(2022)
[50] 阿达米,S。;胡,X。;Adams,N.,平滑粒子流体动力学的输运速度公式,J.Compute。物理。,241, 292-307 (2013) ·兹比尔1349.76659
[51] 利特维诺夫,S。;胡,X。;Adams,N.A.,《保守SPH近似的一致性和收敛性》,J.Compute。物理。,301, 394-401 (2015) ·Zbl 1349.76710号
[52] 胡,X。;Adams,N.,《不可压缩湍流的SPH模型》,Proc。IUTAM,18,66-75(2015)
[53] Rezavand先生。;张,C。;Hu,X.,GPU加速光滑粒子流体力学中的广义有效壁面边界条件处理(2021),arXiv预印本
[54] 张,C。;Zhu,Y。;Lyu,X。;Hu,X.,《SPH的一个有效且通用的固体边界条件:多相流和流体-结构相互作用的应用》,《欧洲力学杂志》。B、 流体,94,276-292(2022)·Zbl 1487.76060号
[55] S.Turek,J.Hron,弹性物体和层流不可压缩流之间流体-结构相互作用的数值基准建议,2007年第53卷,第371-385页·Zbl 1323.76049号
[56] Wendland,H.,分段多项式,正定和紧支集最小次径向函数,高级计算。数学。,4, 389-396 (1995) ·Zbl 0838.41014号
[57] Langlois,W。;Deville,M.,粘性流动方程的精确解,105-143(2014)
[58] 新泽西州昆兰。;Basa先生。;Lastiwka,M.,《无网格粒子方法中的截断误差》,国际期刊《数值》。方法工程,66,13,2064-2085(2006)·Zbl 1110.76325号
[59] Springel,V.,《天体物理学中的平滑粒子流体动力学》,年。阿斯顿牧师。天体物理学。,48, 1, 391-430 (2010)
[60] Ellero,M。;Adams,N.,《限制在通道内的周期性圆柱体阵列周围流动的SPH模拟》,国际期刊Numer。方法工程,861027-140(2011)·Zbl 1235.76126号
[61] 马龙,S。;Colagrossi,A。;安托诺,M。;科利奇奥,G。;Graziani,G.,《低雷诺数和中等雷诺数下物体周围粘性流动的精确SPH建模》,J.Compute。物理。,245, 456-475 (2013) ·Zbl 1349.76715号
[62] I.Federico,S.Marrone,A.Colagrossi,F.Aristodemo,P.Veltri,《通过SPH模拟自由表面渠道流动》,2010年·Zbl 1258.76123号
[63] A.Tafuni,J.Domínguez,R.Vacondio,I.Sahin,A.Crespo,大型SPH模拟的开放边界条件,2016年。
[64] White,F.,《粘性流体流动》(2006),McGraw-Hill:McGraw-Hill纽约
[65] Brehm,C。;哈德尔,C。;Fasel,H.,可压缩Navier-Stokes方程的局部稳定浸没边界法,J.Compute。物理。,295, 475-504 (2015) ·兹比尔1349.76437
[66] Al-Marouf,M。;Samtaney,R.,复杂几何中可压缩流动的通用嵌入边界自适应网格方法,J.Compute。物理。,337, 339-378 (2017) ·Zbl 1415.76495号
[67] 刘,C。;郑,X。;Sung,C.,非定常不可压缩流动的预处理多重网格方法,J.Compute。物理。,139, 35-57 (1998) ·Zbl 0908.76064号
[68] 勒·D。;Khoo,B。;Peraire,J.,《涉及刚性和柔性边界的粘性不可压缩流动的浸没界面法》,J.Compute。物理。,220, 109-138 (2006) ·Zbl 1158.74349号
[69] 罗素·D。;Wang,Z.,《在二维不可压缩粘性流中建模多个运动物体的笛卡尔网格方法》,J.Compute。物理。,191, 177-205 (2003) ·Zbl 1160.76389号
[70] Kunihiko,T。;Tim,C.,《浸没边界法:投影法》,J.Compute。物理。,225, 2118-2137 (2007) ·Zbl 1343.76027号
[71] Jin,G。;Braza,M.,不可压缩非定常Navier-Stokes计算的无反射出口边界条件,计算杂志。物理。,107, 239-253 (1993) ·Zbl 0777.76072号
[72] Tritton,D.J.,《低雷诺数圆柱绕流实验》,J.流体力学。,6, 547-567 (1959) ·Zbl 0092.19502号
[73] Bhardwaj,R。;Mittal,R.,《大型流致变形的耦合浸没边界单元单元解算器基准测试》,AIAA J.,50,1638-1642(2012)
[74] 田,F。;戴,H。;罗,H。;多伊尔,J。;卢梭,B.,《涉及大变形的流体-结构相互作用:三维模拟和生物系统应用》,J.Compute。物理。,258, 451-469 (2014) ·Zbl 1349.76274号
[75] Zhu,Y。;张,C。;Yu,Y。;胡,X.,一种适用于任意复杂几何形状的cad-compatible体填充粒子发生器及其在波-结构相互作用中的应用,J.Hydrodyn。,05 (2021)
[76] 坎普雷格尔,R。;Militzer,J。;曼苏尔,S。;Neto,A.,《使用浸没边界法计算中等雷诺数下静止球体的流动》,J.Braz。Soc.机械。科学。工程,31344-352(2009)
[77] Monaghan,J.J.,《SPH可压缩湍流》,周一。不是。R.阿斯顿。Soc.,335,3843-852(2002)
[78] Magnaudet,J。;里韦罗,M。;Fabre,J.,通过刚性球体或球形气泡的加速流动。第1部分。稳定应变流,J.流体力学。,284, 97-135 (1995) ·Zbl 0848.76063号
[79] Tabata,M。;Itakura,K.,《球体阻力系数的精确计算》,《国际计算杂志》。流体动力学。,9, 303-311 (1998) ·Zbl 0917.76040号
[80] Roos,F.W。;Willmarth,W.W.,关于球体和圆盘阻力的一些实验结果,AIAA J.,9,2,285-291(1971)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。