×

八叉树网格上的稳定节点投影法。 (英文) Zbl 07819049号

摘要:我们提出了一种新的同位投影方法来求解具有任意边界的不可压Navier-Stokes方程。我们的方法使用非分级八叉树网格,其中所有变量都存储在节点上。为了离散粘性和投影步骤,我们使用带尖锐边界处理的超收敛有限差分近似。我们证明了投影在均匀网格上的稳定性,确定了自适应网格上的一个充分稳定性条件,并对这些结果进行了数值验证。通过对不可压缩流体流动的几个典型二维和三维模拟,我们进一步证明了求解器的准确性和功能。总的来说,我们的方法是二阶精度的,允许对任意几何图形进行动态网格自适应,并通过数据搭配减少了代码开发的开销。

MSC公司:

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

参考文献:

[1] M.Aanjaneya。;Patkar,S。;Fedkiw,R.,不可压缩流中气泡的整体质量跟踪公式。J.计算。物理。,17-61 (2013) ·Zbl 1349.76840号
[2] Ahrens,J.P。;Geveci,B。;Law,C.C.W.,Paraview:大数据可视化的最终用户工具
[3] Almgren,A.S。;贝尔,J.B。;科尔拉,P。;Marthaler,T.,复杂几何中不可压缩欧拉方程的笛卡尔网格投影方法。SIAM J.科学。计算。,5, 1289-1309 (1997) ·Zbl 0910.76040号
[4] Almgren,A.S。;贝尔,J.B。;克拉奇菲尔德,W.Y.,近似投影方法:第一部分,无粘分析。SIAM J.科学。计算。,4, 1139-1159 (2000) ·Zbl 0995.76059号
[5] Almgren,A.S。;贝尔,J.B。;Szymczak,W.G.,基于近似投影的不可压缩Navier-Stokes方程的数值方法。SIAM J.科学。计算。,2, 358-369 (1996) ·Zbl 0845.76055号
[6] R.安藤。;Batty,C.,具有自适应表面分辨率的实用八叉树液体模拟器。ACM事务处理。图表。,4(2020年8月)
[7] Aslam,T.D.,多维外推的偏微分方程方法。J.计算。物理。,1, 349-355 (2004) ·Zbl 1036.65002号
[8] 巴蒂,C。;贝塔尔斯,F。;Bridson,R.,精确固液耦合的快速变分框架。ACM事务处理。图表。,3(2007年7月),100-es·兹比尔1392.74035
[9] 伯杰,M.J。;Oliger,J.,双曲偏微分方程的自适应网格精化。J.计算。物理。,3, 484-512 (1984) ·Zbl 0536.65071号
[10] 布拉萨,M。;Chassaing,P。;Minh,H.,圆柱体近尾流中压力场和速度场的数值研究和物理分析。流体力学杂志。,79-130 (1986) ·Zbl 0596.76047号
[11] Brown,D.L。;科尔特斯,R。;Minion,M.L.,不可压缩Navier-Stokes方程的精确投影方法。J.计算。物理。,2, 464-499 (2001) ·Zbl 1153.76339号
[12] Calhoun,D.,求解不规则区域中二维河流功能-危险性方程的笛卡尔网格方法。J.计算。物理。,2, 231-275 (2002) ·Zbl 1130.76371号
[13] Cho,H。;Kang,M.,两相不可压缩Navier-Stokes方程跳跃条件的完全隐式精确处理。J.计算。物理学。(2021) ·Zbl 07515846号
[14] Cho,H。;Park,Y。;Kang,M.,用整体法求解不规则区域和四叉树上的不可压缩Navier-Stokes方程。J.计算。物理学。(2022) ·Zbl 07536793号
[15] Chorin,A.J.,解决不可压缩粘性流动问题的数值方法。J.计算。物理。,1, 12-26 (1967) ·Zbl 0149.44802号
[16] Cleret de Langavant,C。;吉特·A。;米歇尔·泰拉德。;Temprano-Coleto,F。;Gibou,F.,可溶表面活性剂驱动流动的水平集模拟。J.计算。物理。,271-297 (2017) ·Zbl 1380.76079号
[17] 达贡,L。;Menon,R.,Openmp:共享内存编程的行业标准api。IEEE计算。科学。工程师,146-55(1998年1月)
[18] A.杜比。;Almgren,A。;贝尔·J。;Berzins,M。;Brandt,S。;布莱恩,G。;科尔拉,P。;格雷夫斯,D。;Lijewski,M。;Löffler,F.,块结构自适应网格优化包中高级框架的调查。J.平行分布计算。,12, 3217-3227 (2014)
[19] 杜奇,H。;杜斯特,F。;贝克尔,S。;Lienhart,H.,在低Keulegan-Capenter数下绕振荡圆柱体的低雷诺数流动。流体力学杂志。,249-271 (1998) ·Zbl 0922.76024号
[20] 伊根·R。;吉特·A。;Temprano-Coleto,F。;艾萨克,T。;Peaudecerf,F.J。;兰德尔,J.R。;Luzzatto-Fegiz,P。;Burstede,C。;Gibou,F.,平行八叉树网格上不可压缩流动的直接数值模拟。J.计算。物理学。(2021) ·Zbl 07511438号
[21] Engelman,医学硕士。;Jamnia,M.-A.,《圆柱体瞬态流动:基准解决方案》。国际期刊数字。液体方法,7985-1000(1990)
[22] 恩赖特,D。;Nguyen,D。;Gibou,F。;Fedkiw,R.,使用粒子水平集方法和自由表面流动的二阶精确压力边界条件。专题讨论会,A、B和C部分,337-342(2003年7月)
[23] Erturk,E。;科克,T.C。;Gököl,C.,高雷诺数下二维定常不可压缩驱动空腔流的数值解。国际期刊数字。方法流体,7747-774(2005)·Zbl 1071.76038号
[24] 风扇,D。;Yang,L。;王,Z。;Triantafylou,M.S。;Karniadakis,G.E.,实验和模拟中钝体主动流控制的强化学习。程序。国家。阿卡德。科学。,42, 26091-26098 (2020)
[25] Ghia,美国。;Ghia,K.N。;Shin,C.,使用Navier-Stokes方程和多重网格方法求解不可压缩流的高分辨率解。J.计算。物理。,3, 387-411 (1982) ·Zbl 0511.76031号
[26] Gibou,F。;Min,C.,流体/固体相互作用的高效对称正定二阶精确单片解算器。J.计算。物理。,8, 3246-3263 (2012) ·Zbl 1400.76049号
[27] Goda,K.,计算二维或三维空腔流动的隐式差分多步技术。J.计算。物理。,1 (1979-1)
[28] 吉特·A。;米歇尔·泰拉德。;Gibou,F.,在任意几何体和自适应四/八叉树上求解不可压缩Navier-Stokes方程的稳定投影方法。J.计算。物理。,215-238 (2015) ·Zbl 1349.76336号
[29] Gómez,P。;Zanzi,C。;López,J。;Hernández,J.,基于单元顶点/边缘的交错八叉树网格上高密度比界面流动的模拟,在不规则节点处进行二阶离散化。J.计算。物理。,478-507 (2019) ·Zbl 1416.76037号
[30] 哈洛,F.H。;Welch,J.E.,含自由表面流体随时间变化的粘性不可压缩流动的数值计算。物理学。流体,12182-2189(1965)·Zbl 1180.76043号
[31] Harten,A。;Osher,S。;Engquist,B。;Chakravarthy,S.R.,一致高阶精确本质非振荡格式的一些结果。申请。数字。数学。,3-5, 347-377 (1986) ·Zbl 0627.65101号
[32] 徐,M.-C。;卡门斯基,D。;巴齐列夫斯,Y。;萨克斯,M.S。;Hughes,T.J.,《生物假体心脏瓣膜的流体-结构相互作用分析:动脉壁变形的意义》。计算。机械。,1055-1071 (2014) ·Zbl 1311.74039号
[33] 海德,D.A。;Fedkiw,R.,子网格和更高分辨率固体整体固液耦合的统一方法。J.计算。物理。,490-526 (2019) ·Zbl 1452.76255号
[34] 伊利亚迪斯,G。;Anagostopoulos,P.,低Keulegan-Carpenter数和频率参数下绕圆柱的粘性振荡流。国际期刊数字。液体方法,4403-442(1998)·Zbl 0909.76049号
[35] 约翰逊,T.A。;帕特尔,V.C.,流经一个雷诺数高达300的球体。流体力学杂志。,19-70 (1999)
[36] Karniadakis,G.E。;以色列,M。;Orszag,S.A.,不可压缩Navier-Stokes方程的高阶分裂方法。J.计算。物理。,2, 414-443 (1991) ·Zbl 0738.76050号
[37] Kim,J。;Moin,P.,分数步法在不可压缩Navier-Stokes方程中的应用。J.计算。物理。,2, 308-323 (1985) ·Zbl 0582.76038号
[38] 库切洛娃,A。;斯特兰戈,S。;Sukenik,S。;Theillard,M.,蛋白质构象变化的计算模型——应用于开放的sars-cov-2峰。J.计算。物理学。(2021) ·Zbl 07515474号
[39] Kwok,S.,任意网格系统的改进曲线有限差分(CFD)方法。计算。结构。,4, 719-731 (1984) ·Zbl 0532.65067号
[40] 勒克莱尔,B。;哈米利克,A。;Pruppacher,H.,低雷诺数和中雷诺数下球体阻力的数值研究。J.大气。科学。,2, 308-315 (1970)
[41] 刘晓东。;Osher,S。;Chan,T.,加权基本非振荡格式。J.计算。物理。,1, 200-212 (1994) ·Zbl 0811.65076号
[42] 龙,X。;Chen,C.,对流扩散问题的二阶半拉格朗日方法的一般形式。摘要与应用分析(2013),Hindawi·Zbl 1266.65156号
[43] Losasso,F。;Fedkiw,R。;Osher,S.,水平集方法和不可压缩流的空间自适应技术。计算。流体,10995-1010(2006)·Zbl 1177.76295号
[44] Losasso,F。;Gibou,F。;Fedkiw,R.,用八叉树数据结构模拟水和烟雾,457-462
[45] 曼奇尼,M。;米歇尔·泰拉德。;Kim,C.,波动流体动力学方程的投影法。J.计算。物理学。(2022) ·Zbl 07536778号
[46] Marella,S。;Krishnan,S。;刘,H。;Udaykumar,H.,夏普界面笛卡尔网格方法I:一种易于实现的三维移动边界计算技术。J.计算。物理。,1, 1-31 (2005) ·Zbl 1154.76359号
[47] 最小值C。;Gibou,F.,非梯度自适应网格上不可压缩Navier-Stokes方程的二阶精确投影方法。J.计算。物理。,2, 912-929 (2006) ·Zbl 1330.76096号
[48] 最小值C。;Gibou,F.,不规则域上的几何积分及其在水平集方法中的应用。J.计算。物理。,1432-1443 (2007) ·兹比尔1125.65021
[49] 最小值C。;Gibou,F。;Ceniceros,H.D.,非梯度网格上变系数泊松方程的超收敛有限差分格式。J.计算。物理。,1, 123-140 (2006) ·Zbl 1106.65093号
[50] Minion,M.L.,局部细化网格的投影方法。J.计算。物理。,1, 158-178 (1996) ·Zbl 0859.76047号
[51] 米塔尔·R。;Dong,H。;博兹库塔斯,M。;Najjar,F。;瓦尔加斯,A。;von Loebbecke,A.,一种适用于复杂边界不可压缩流的通用锐界面浸没边界法。J.计算。物理。,10, 4825-4852 (2008) ·Zbl 1388.76263号
[52] Mullowney,P.等人。;李,R。;托马斯·S。;阿南坦,S。;Sharma,A。;Rood,J.S。;威廉姆斯,A.B。;Sprague,M.A.,为exascale级风能模拟编制不可压缩流体动力学代码,1-16
[53] Ng,Y.T。;最小值C。;Gibou,F.,单相流的高效流固耦合算法。J.计算。物理。,23, 8807-8829 (2009) ·Zbl 1245.76019号
[54] OpenMP Architecture Review Board(2018年5月),OpenMP应用程序接口版本5.0
[55] Orszag,S.A.公司。;以色列,M。;Deville,M.O.,《不可压缩流动的边界条件》。科学杂志。计算。,75-111 (1986) ·Zbl 0648.76023号
[56] Pazner,W。;Persson,P.-O.,非常高阶间断Galerkin方法的近似张量积预条件。J.计算。物理。,344-369 (2018) ·Zbl 1380.65067号
[57] Popinet,S.,Gerris:复杂几何中不可压缩Euler方程的基于树的自适应求解器。J.计算。物理。,2, 572-600 (2003) ·Zbl 1076.76002号
[58] Ryzhakov,P。;Cotela,J。;罗西,R。;Oñate,E.,高效模拟不可压缩流动的两步整体方法。国际期刊数字。《液体方法》,12919-934(2014)·Zbl 1455.76093号
[59] Samet,H.,《四叉树、八叉树和相关分层数据结构概述》(1988)
[60] Saye,R.I。;Sethian,J.A。;彼得鲁斯基,B。;扎托斯基,A。;卢,X。;Rock,R.,《工业旋转钟雾化的高保真建模的见解》。程序。国家。阿卡德。科学。,4 (2023)
[61] 肖恩鲍姆,F。;Rüde,U.,极端规模块结构自适应网格细化。SIAM J.科学。计算。,3、C358-C387(2018)·Zbl 06890193号
[62] Seo,J.H。;Mittal,R.,一种改进质量守恒和减少虚假压力振荡的尖锐界面浸没边界法。J.计算。物理。,19, 7347-7363 (2011) ·Zbl 1408.76162号
[63] 高桥,T。;Batty,C.,Monolith:强双向刚性-刚性-刚性-流体耦合的单片压力-粘度接触解算器(2020)
[64] Temam,R.,Naver stokes方程解的Sur’近似值与分数方法相同(i)。架构(architecture)。定额。机械。分析。,135-153 (1969) ·兹比尔0195.46001
[65] Temam,R.,Sur l’a近似解navier-stokes par la mémethode des pas fractionnaires方程(ii)。架构(architecture)。定额。机械。分析。,135-153 (1969) ·Zbl 0195.46001号
[66] Theillard,M.,用于水平集表示的卷表示参考图方法。J.计算。物理学。(2021) ·Zbl 07513799号
[67] 米歇尔·泰拉德。;阿隆索·马蒂拉,R。;Saintillan,D.,主动集体运动的几何控制。软物质,363-375(2017)
[68] 米歇尔·泰拉德。;Djodom,L.F。;维耶,J.-L。;Gibou,F.,一种求解不规则区域和自适应网格上线性弹性方程的二阶尖锐数值方法——应用于形状优化。J.计算。物理。,430-448 (2013) ·Zbl 1286.74109号
[69] 米歇尔·泰拉德。;Gibou,F。;Saintillan,D.,不可压缩两相流的Sharp数值模拟。J.计算。物理。,91-118 (2019) ·Zbl 1452.76135号
[70] 米歇尔·泰拉德。;赖克罗夫特,C.H。;Gibou,F.,非分级自适应八叉树和四叉树笛卡尔网格上的多重网格方法。科学杂志。计算。,1-15 (2013) ·Zbl 1273.65188号
[71] 米歇尔·泰拉德。;Saintillan,D.,受限活动流体的计算平均场建模。J.计算。物理学。(2019) ·Zbl 1453.76143号
[72] Van Kan,J.,粘性不可压缩流的二阶精确压力校正方案。SIAM J.科学。统计计算。,3, 870-891 (1986) ·Zbl 0594.76023号
[73] 秀,D。;Karniadakis,G.E.,Navier-Stokes方程的半拉格朗日高阶方法。J.计算。物理。,2, 658-684 (2001) ·Zbl 1028.76026号
[74] 张伟。;Almgren,A。;贝克纳,V。;贝尔·J。;布拉斯科,J。;Chan,C。;Day,M。;弗里森,B。;哥特,K。;Graves,D.,Amrex:块结构自适应网格细化框架。J.开源。软质。,37, 1370 (2019)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。