×

一个简单有效的不可压缩Navier-Stokes解算器,用于截断区域上的非定常复杂几何流。 (英语) Zbl 1390.76562号

摘要:基于投影方法的不可压缩Navier-Stokes解算器通常需要对压力类变量的泊松方程进行昂贵的数值解。这通常涉及基于迭代和多重网格方法的线性系统求解器,这可能会限制扩展到大量处理器的能力。人工压缩法(ACM)[A.J.乔林,J.计算。物理学。135,第2期,118-125页(1997年;Zbl 0899.76283号)]将压力的时间导数引入连续性方程的不可压缩形式,创建一个耦合的闭合双曲系统,该系统不需要泊松方程解,并允许显式时间推进和局部模板数值方法。理论上,这种方案应该能够在大量CPU、GPU或CPU-GPU混合体系结构上很好地扩展。最初的ACM仅适用于稳定流动,双时间步长通常用于时间精确的模拟。最近,在[“用于显式模拟不可压缩流的人工压缩性的熵阻尼形式”中,Phys.Rev.(3)E 87,No.1,Article ID 013309,12 p.(2013;doi:10.1103/physreve.87.013309)],J.R.公司。克劳森提出了熵阻尼人工可压缩性(EDAC)方法,该方法适用于定常和非定常流动,无需双重时间步长。用有限差分MacCormack方法对二维盖驱动腔和周期性双剪切层问题进行了测试,并用有限元方法对方柱绕流进行了测试。在本研究中,我们使用新的优化高阶中心有限差分格式和显式四阶Runge-Kutta方法离散EDAC公式。这与浸没边界方法相结合,可以有效地处理复杂的几何形状,并与新的鲁棒流出边界条件相结合,可以在截断域上进行更高的雷诺数模拟。对二维和三维泰勒-格林涡问题和盖驱动空腔问题进行了验证研究。使用涡流粘性亚脊尺度模型对三维情况进行大涡模拟。最后,给出了一个在球体上流动的应用程序,以突出边界条件,并对三维盖驱动腔与传统的不可压缩Navier-Stokes解算器进行了性能比较。总的来说,EDAC公式和离散化的组合既有效又经济。

MSC公司:

76M20码 有限差分方法在流体力学问题中的应用
76M10个 有限元方法在流体力学问题中的应用
6500万06 含偏微分方程初值和初边值问题的有限差分方法
65M60毫米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
76-04 流体力学相关问题的软件、源代码等
76D05型 不可压缩粘性流体的Navier-Stokes方程

软件:

炒作
PDF格式BibTeX公司 XML格式引用
全文: DOI程序

参考文献:

[1] Anupindi,K。;德洛姆,Y。;Frankel,S。;Shetty,D.,一种用于复杂动脉血流动力学大涡模拟的新型多块浸没边界法,《计算物理杂志》,254200-218,(2013)·Zbl 1349.76112号
[2] 转向架,C。;Bailly,C.,流量和噪声计算的一系列低色散和低耗散显式格式,《计算物理杂志》,194194-214,(2004)·兹比尔1042.76044
[3] Brachet,M.,泰勒-格林涡旋三维湍流的直接数值模拟,流体动力学研究,8,1-9,(1991)
[4] Chen,K。;Pletcher,R.,使用强隐式计算程序模拟三维液体晃动流动,AIAA J,31,901-910,(1993)·Zbl 0774.76057号
[5] Choi,D。;Merkle,C.,时间迭代格式在不可压缩流动中的应用,AIAA J,231518-1524,(1985)·Zbl 0571.76016号
[6] Chorin,A.,求解不可压缩粘性流问题的数值方法,计算机物理杂志,135118-125,(1967)·Zbl 0899.76283号
[7] Clausen,J.,显式模拟不可压缩流的人工压缩性熵阻尼形式,Phys Rev E,87,13309-13321,(2013)
[8] 德洛姆,Y。;Anupindi,K。;科尔洛,A。;Shetty,D。;罗德菲尔德,M。;Chen,J.,动力方丹血流动力学的大涡模拟,生物医学杂志,46408-422,(2013)
[9] 德洛姆,Y。;Anupindi,K。;Frankel,S.,FDA理想医疗设备的大涡模拟,心血管工程技术,4,4,(2013)
[10] Dong,S。;卡尼亚达基斯,G。;Chryssostomidis,C.,多个截断无界区域上不可压缩流模拟的稳健且准确的流出边界条件,计算物理杂志,261,83-105,(2014)·Zbl 1349.76569号
[11] Fadlun,E。;Verzicco,R。;奥兰迪,P。;Mohd-Yusof,J.,三维复杂流动模拟的组合浸没边界有限差分方法,计算物理杂志,161,35-60,(2000)·Zbl 0972.76073号
[12] Falgout R,Jones J.Multigrid关于大规模并行架构的研究·Zbl 0972.65110号
[13] Falgout R,Yang。Hypre:一个高性能预处理程序库。2002. ·Zbl 1056.65046号
[14] http://www.acts.nersc.gov/hypre/。
[15] Ghia,美国。;Ghia,K。;Shin,C.,使用Navier-Stokes方程和多重网格方法求解不可压缩流的高分辨率解,《计算物理杂志》,48,387-411,(1982)·Zbl 0511.76031号
[16] Hafez M,South J,Murman E.跨声速全势方程数值解的人工压缩方法。美国汽车协会期刊17:838-844·Zbl 0409.76013号
[17] Hartwich,P。;Hsu,C.,三维不可压缩Navier-Stokes方程的高分辨率迎风格式,AIAA J,26,1321-1328,(1988)·Zbl 0664.76034号
[18] 哈特维奇,P。;徐,C。;Liu,C.,通量差分裂三维Navier-Stokes方程的矢量化隐式算法,流体工程杂志,110,297-305,(1988)
[19] 桥本,T。;塔诺,I。;Yasuda,T。;田中,Y。;Morinishi,K。;Satofuka,N.,在GPU上使用动力学简化的局部Navier-Stokes方程对非定常粘性不可压缩流动进行高阶数值模拟,计算流体,110,108-113,(2015)·兹比尔1390.76718
[20] 夸克,D。;Chang,J。;柄,P。;Chakravarthy,S.,使用基本变量的三维不可压缩Navier-Stokes流解算器,AIAA J,24,390-396,(1986)·Zbl 0587.76044号
[21] Leriche,E。;Gavrilakis,S.,《盖驱动立体腔中流动的直接数值模拟》,《物理流体》,第12期,第1363-1376页,(2000年)·Zbl 1149.76452号
[22] 林德斯,V。;Nordstrom,J.,《空间中心差分算子的一致最佳波数近似》,《复合物理杂志》,300,695-709,(2015)·Zbl 1349.65566号
[23] Madsen,P。;Shaffer,H.,《人工压缩性的讨论》,《海岸工程》,53,93-98,(2006)
[24] 小仆,M。;Brown,D.,欠分辨率二维不可压缩流动模拟的性能,II,J Comput Phys,138734-765,(1997)·Zbl 0914.76063号
[25] 马克·A。;van Wachem,B.,新型隐式二阶精确浸没边界法的推导与验证,计算物理杂志,22760-6680,(2008)·Zbl 1338.76098号
[26] Pan D,Chakravarthy S.不可压缩流的统一公式。1989年AIAA。
[27] Pope,S.,《湍流》(2000),剑桥大学出版社·Zbl 0966.76002号
[28] Rogers S,Kwak D,Kiris C.稳态和时间相关问题的不可压缩navier-stokes方程的数值解。1989年AIAA。
[29] Tsynkov,S.,无界域上问题的数值解:综述,应用数值数学,27465-532,(1998)·Zbl 0939.76077号
[30] Shetty,D。;费希尔,T。;Chunekar,A。;Frankel,S.,《完全非均匀湍流的高阶不可压缩大涡模拟》,《计算物理杂志》,2298802-8822,(2010)·Zbl 1282.76115号
[31] 斯特格,J。;Kulter,P.,《涡尾迹计算的隐式有限差分程序》,AIAA J,15,581-590,(1977)·Zbl 0381.76030号
[32] Turkel,E.,《流体动力学预处理方法综述》,应用数值数学,12,257-284,(1993)·Zbl 0770.76048号
[33] Vrahliotis,S。;Pappou,T。;Tsangaris,S.,混合网格非定常不可压缩流的人工压缩三维Navier-Stokes解算器,Eng Appl Comput Fluid Mech,6248-270,(2012)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。