×

采用无发散速度插值和力扩散的浸没边界法。 (英语) Zbl 1380.76078号

摘要:浸没边界(IB)方法是一个数学框架,用于构建稳健的数值方法,以研究浸没在粘性流体中的弹性结构问题中的流体-结构相互作用。IB公式使用流体的欧拉表示和结构的拉格朗日表示。拉格朗日框架和欧拉框架通过带δ函数核的积分变换耦合。离散的IB方程使用带正则化δ函数核的这些变换的近似值来插值结构的流体速度,并将结构力传播到流体。众所周知,传统IB方法的体积守恒性较差,因为插值的拉格朗日速度场通常不是无发散的,因此这可能导致虚假的体积变化。实际上,在薄结构边界上存在较大压力差的情况下,缺乏体积守恒尤其明显。本文的目的是通过引入速度插值和力扩散格式,使结构运动的插值速度场至少为(mathcal{C}^1),并满足连续无发散条件,从而大大减小IB方法的体积误差,力扩散算子是速度插值算子的伴随。我们通过二维和三维空间的数值实验证实,与其他现有IB方法相比,这种新的IB方法能够在体积守恒方面取得实质性改进,但代价是计算成本略有增加。此外,与传统的IB方法相比,新方法提供了更平滑的拉格朗日力(牵引力)。这里提出的方法仅限于周期计算域。它对非周期域的推广是未来的重要工作。

MSC公司:

76M20码 有限差分方法在流体力学问题中的应用
74层10 流固相互作用(包括气动和水弹性、孔隙度等)
76D05型 不可压缩粘性流体的Navier-Stokes方程
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Bao,Y。;Kaiser,A.D。;Kaye,J。;Peskin,C.S.,具有三个连续导数和改进的平移不变性的高斯型浸入边界核(2015)
[2] Bao,Y。;凯伊,J。;Peskin,C.S.,具有三个连续导数和改进的平移不变性的类高斯浸入边界核,J.Compute。物理。,316, 139-144 (2016) ·Zbl 1349.76585号
[3] Bhalla,A.P.S。;右侧捆。;格里菲斯,B.E。;Patankar,N.A.,《流体-结构与刚体、变形体和弹性体相互作用的统一数学框架和自适应数值方法》,J.Compute。物理。,250, 446-476 (2013) ·Zbl 1349.65403号
[4] Boffi,D。;Gastaldi,L。;赫尔泰,L。;Peskin,C.,关于浸入边界法的超弹性公式,计算。方法应用。机械。工程,197,25-28,2210-2231(2008)·Zbl 1158.74523号
[5] 科尔特斯,R。;Minion,M.,浸没边界问题的斑点投影法,J.Compute。物理。,161, 2, 428-453 (2000) ·Zbl 0962.74078号
[6] 科尔特斯,R。;Peskin,C.S。;Stockie,J.M。;Varela,D.,浸没弹性边界中的参数共振,SIAM J.Appl。数学。,65, 2, 494-520 (2004) ·Zbl 1074.74024号
[7] Devendran,D。;Peskin,C.S.,《基于浸没边界能量的不可压缩粘弹性方法》,J.Compute。物理。,231, 14, 4613-4642 (2012) ·Zbl 1245.76099号
[8] Dutt,A。;Rokhlin,V.,非等空间数据的快速傅立叶变换,SIAM J.Sci。计算。,1368-1393年6月14日(1993年)·Zbl 0791.65108号
[9] Fai,T.G。;格里菲斯,B.E。;Mori,Y。;Peskin,C.S.,《使用快速恒效线性求解器求解变粘度和变密度问题的浸没边界法I:数值方法和结果》,SIAM J.Sci。计算。,35、5、B1132-B1161(2013)·Zbl 1282.76088号
[10] 高,H。;Wang,H。;Berry,C。;罗,X。;Griffith,B.E.,基于准静态图像的左心室舒张负荷下浸没边界无限元模型,国际期刊数字。方法生物识别。工程,30,11,1199-1222(2014)
[11] A.戈扎。;利斯卡,S。;莫利,B。;Colonius,T.,《用浸没边界法精确计算表面应力和力》,J.Compute。物理。,321, 860-873 (2016) ·Zbl 1349.76466号
[12] Greengard,L。;Lee,J.,《加速非均匀快速傅里叶变换》,SIAM Rev.,46,3,443-454(2004)·Zbl 1064.65156号
[13] Griffith,B.,用投影法作为预条件,求解不可压缩Navier-Stokes方程的一种精确而有效的方法,J.Compute。物理。,228, 20, 7565-7595 (2009) ·Zbl 1391.76474号
[14] Griffith,B.,《生理驱动和负载条件下主动脉瓣动力学的浸没边界模型》,国际期刊数值。方法生物识别。工程,28,317-345(2012)·Zbl 1243.92017年
[15] 格里菲斯,B。;Hornung,R。;McQueen,D。;Peskin,C.,浸没边界法的一种自适应形式上的二阶精确版本,J.Compute。物理。,223、1、10-49(2007),软件可从·Zbl 1163.76041号
[16] 格里菲斯,B。;罗,X。;McQueen,D。;Peskin,C.,《使用浸没边界法模拟天然和人工心脏瓣膜的流体动力学》,《国际应用杂志》。机械。,1, 01, 137-177 (2009)
[17] Griffith,B.E.,《关于浸没边界法的体积守恒》,Commun。计算。物理。,12, 2, 401-432 (2012) ·Zbl 1373.74098号
[18] 格里菲斯,B.E。;Luo,X.,混合有限差分/有限元浸没边界法,国际J数值。方法生物识别。工程(2017)
[19] 格里菲斯,B.E。;Peskin,C.S.,关于浸入边界法的精度阶数:足够光滑问题的高阶收敛率,J.Comput。物理。,208,175-105(2005年)·Zbl 1115.76386号
[20] Kallemov,B。;Bhalla,A.P.S。;格里菲斯,B.E。;Donev,A.,《刚体浸没边界法》,Commun。申请。数学。计算。科学。,11、1、79-141(2016),软件可在·Zbl 1382.76191号
[21] Kim,Y。;Lai,M.-C。;Peskin,C.S.,用浸没边界法对二维泡沫进行数值模拟,J.Compute。物理。,229, 13, 5194-5207 (2010) ·Zbl 1346.76139号
[22] 金,Y。;Lai,M.-C。;Peskin,C.S。;Seol,Y.,《浸没边界法三维泡沫的数值模拟》,J.Compute。物理。,269, 1-21 (2014) ·Zbl 1349.76619号
[23] Ko,W。;Stockie,J.M.,“浸没弹性边界中的参数共振”修正(2012)
[24] Ko,W。;Stockie,J.M.,球形浸没弹性壳的参数共振,SIAM J.Appl。数学。,76, 1, 58-86 (2016) ·Zbl 1339.74008号
[25] 赖,M.-c。;Li,Z.,关于三维Navier-Stokes方程涉及浸入式移动膜的跳跃条件的注释,Appl。数学。莱特。,14, 2, 149-154 (2001) ·Zbl 1013.76021号
[26] Lai,M.-C。;Peskin,C.S.,《一种形式上具有二阶精度和降低数值粘性的浸没边界法》,J.Compute。物理。,160, 2, 705-719 (2000) ·Zbl 0954.76066号
[27] Lee,L。;LeVeque,R.J.,《不可压缩Navier-Stokes方程的浸没界面法》,SIAM J.Sci。计算。,25, 3, 832-856 (2003) ·Zbl 1163.65322号
[28] 李,Z。;Lai,M.-C.,具有奇异力的Navier-Stokes方程的浸没界面法,J.Compute。物理。,171, 2, 822-842 (2001) ·Zbl 1065.76568号
[29] 里斯卡,S。;Colonius,T.,《使用格点格林函数求解外部不可压缩粘性流的快速浸没边界法》,J.Compute。物理。,331, 257-279 (2017) ·Zbl 1378.76083号
[30] 卢氏,E。;Peskin,C.S.,包含大量相互作用的微型游泳运动员的主动悬架的建模和模拟,计算机。结构。,122239-248(2013年)
[31] McQueen,D。;Peskin,C.,用于计算跳动哺乳动物心脏血流的浸没边界法的共享内存并行向量实现,J.Supercomput。,11, 3, 213-236 (1997)
[32] Peskin,C.S.,《心脏瓣膜周围的流动模式:数值方法》,J.Compute。物理。,10, 2, 252-271 (1972) ·Zbl 0244.9202号
[33] Peskin,C.S.,《心脏血流的数值分析》,J.Compute。物理。,25, 3, 220-252 (1977) ·Zbl 0403.76100号
[34] Peskin,C.S.,浸没边界法,数值学报。,2002年1月11日,479-517(2003)·Zbl 1123.74309号
[35] Peskin,C.S。;Printz,B.F.,《浸没弹性边界流动计算中的改进体积守恒》,J.Compute。物理。,105, 1, 33-46 (1993) ·兹比尔0762.92011
[36] 罗马,A.M。;Peskin,C.S。;Berger,M.J.,浸没边界法的自适应版本,J.Compute。物理。,153, 509-534 (1999) ·Zbl 0953.76069号
[37] Stein,D.B。;盖伊,R.D。;Thomases,B.,浸没边界光滑扩展:用傅里叶谱方法求解任意光滑区域上{PDE}的高阶方法,J.Compute。物理。,304, 252-274 (2016) ·Zbl 1349.65656号
[38] Stein,D.B。;盖伊,R.D。;Thomases,B.,浸没边界光滑扩展(IBSE):求解任意光滑区域中不可压缩流动的高阶方法,J.Compute。物理。,335, 155-178 (2017) ·Zbl 1375.76038号
[39] 斯特里查尔斯基,W。;Guy,R.D.,泡状细胞的细胞内压力动力学,生物物理学。J.,110,5,1168-1179(2016)
[40] Unser,M.,《样条:信号和图像处理的完美结合》,IEEE信号处理。Mag.,16,6,22-38(1999)
[41] Unser,M。;阿尔德鲁比,A。;Eden,M.,关于B样条小波到Gabor函数的渐近收敛性,IEEE Trans。《信息论》,38,2864-872(1992)·Zbl 0757.41022号
[42] Usayaga,F.B。;贝尔,J.B。;Delgado Buscalioni,R。;Donev,A。;Fai,T.G。;格里菲斯,B.E。;Peskin,C.S.,波动流体动力学的交错方案,多尺度模型。模拟。,10, 4, 1369-1408 (2012) ·Zbl 1310.76108号
[43] Balboa Usabiaga,F。;Kallemov,B。;Delmotte,B。;Bhalla,A.P.S。;格里菲斯,B.E。;Donev,A.,《被动和主动刚性颗粒悬浮液的流体动力学:刚性多气泡方法》,Commun。申请。数学。计算。科学。,11,21217-296(2016),软件可在
[44] 威廉姆斯,H.A。;Fauci,L.J。;Gaver,D.P.,使用浸没边界法评估界面流体动态应力,离散Contin。动态。系统。,序列号。B、 11、2、519(2009)·Zbl 1277.76132号
[45] 杨,X。;张,X。;李,Z。;He,G.-W.,《离散δ函数的平滑技术及其在移动边界模拟中浸没边界法的应用》,J.Compute。物理。,228, 20, 7821-7836 (2009) ·Zbl 1391.76590号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。