×

一种用于飞行昆虫湍流多尺度模拟的小波自适应方法。 (英语) Zbl 1473.65101号

摘要:我们提出了一种基于小波的自适应方法,用于计算复杂、依赖时间的几何体中的三维多尺度流,并在大规模并行计算机上实现。虽然我们的重点是模拟扑翼昆虫,但它可以用于其他流动问题。为了避免求解椭圆问题,我们采用人工压缩方法对不可压缩流体进行建模。使用体积惩罚施加无滑移和流入/流出边界条件。控制方程在具有中心有限差分的局部均匀笛卡尔网格上离散,并用四阶Runge-Kutta格式及时积分。将域划分为具有不同分辨率的三次块,并对每个块使用双正交插值小波作为细化指标和预测算子。通过对小波系数进行阈值化,可以生成动态演化的网格,并且自适应策略可以在空间和尺度上跟踪解。块分布在MPI进程之间,网格拓扑使用树状数据结构进行编码。通过分析不同的物理和数值参数,我们可以平衡它们的误差,从而确保最佳收敛,同时最小化计算工作量。不同的验证测试对我们新的开源代码WABBIT的准确性和性能进行了评分。扑动昆虫的流动模拟证明了它对复杂的、受生物启发的问题的适用性。

理学硕士:

6500万06 含偏微分方程初值和初边值问题的有限差分方法
65M50型 涉及偏微分方程初值和初边值问题数值解的网格生成、细化和自适应方法
65M55型 多网格方法;涉及偏微分方程初值和初边值问题的区域分解
65M85型 求解偏微分方程初值和初边值问题的虚拟域方法
65T60型 小波的数值方法
76M20码 有限差分方法在流体力学问题中的应用
第76天05 不可压缩粘性流体的Navier-Stokes方程
76Z10号 水和空气中的生物推进
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] S.Abhyankar、J.Brown、E.M Constantinescu、D.Ghosh、B.F.Smith和H.Zhang。Petsc/ts:一个现代的可扩展ODE/DAE解算器库。arXiv:1806.014372018年。
[2] P.Angot、C.Bruneau和P.Fabrie。一种考虑不可压缩粘性流中障碍物的惩罚方法。数字。数学。,81:497-520, 1999. ·Zbl 0921.76168号
[3] E.Arquis和J.-P.Caltagirone。表面条件流体动力学与界面环境流体环境孔隙:对流性质的应用。C.R.学院。科学。巴黎,Sér。II、 1984年,第299页。
[4] R.Becker和R.Rannacher。后验误差估计的最优控制方法,第1-102页。剑桥大学出版社,剑桥,2001年·Zbl 1105.65349号
[5] M.Berger和P.Colella。激波流体动力学的局部自适应网格细化。J.计算。物理。,82:64-84, 1988. ·Zbl 0665.76070号
[6] M.Berger和J.Oliger。双曲型偏微分方程的自适应网格加密。J.计算。物理。,53:484-512, 1984. ·兹伯利0536.65071
[7] A.勃兰特。边值问题的多级自适应解。数学。公司。,31:333-390, 1977. ·Zbl 0373.65054号
[8] C.Burstede、L.C.Wilcox和O.Ghattas。Ô×Ø:八叉树森林上并行自适应网格细化的可缩放算法。SIAM J.科学。计算。,33(3):1103-1133, 2011. ·兹比尔1230.65106
[9] A.J.Chorin。解决不可压缩粘性流动问题的数值方法。J.计算。物理。,2:12-26, 1967. ·Zbl 0149.44802号
[10] A.科恩。数值分析中的小波方法,《数值分析手册》第七卷。爱思唯尔,阿姆斯特丹,2000年。
[11] A.Cohen、W.Dahmen和R.DeVore。椭圆算子方程的自适应小波方法-收敛速度。数学。Comp,70:27-751998年·兹伯利0980.65130
[12] A.Cohen、I.Daubechies和J.C.Feauveau。紧支撑小波的双正交基。普通纯应用程序。数学。,45:485-560, 1992. ·Zbl 0776.42020号
[13] A.Cohen、N.Dyn、S.M.Kaber和M.Postel。标量守恒律三角形的多分辨率格式。J.公司。物理。,161:264-286, 2000. ·Zbl 0959.65105号
[14] 一、Daubechies。紧支撑小波的正交基。普通纯应用程序。数学。,41(7):909-996, 1988. ·兹伯利0644.42026
[15] R.迪丁。多维爆轰结构的并行自适应仿真。科特布斯大学Brandenburgische Technische Universityät Cottbus博士论文,2003年。
[16] R.迪丁。分布式存储计算机amr算法的构造和应用。Tomasz Plewa、Timur Linde和V.Gregory Weirs主编,《自适应网格优化-理论和应用》,第361-372页,柏林,海德堡,2005年。施普林格-柏林-海德堡·Zbl 1065.65114号
[17] R.Deiterding、M.Domingues、S.Gomes、O.Roussel和K.Schneider。自适应多分辨率还是自适应网格细化?二维欧拉方程的案例研究。ESAIM:进展,29:28-422009·Zbl 1301.76058号
[18] R.Deiterding、M.Domingues、S.Gomes、O.Roussel和K.Schneider。自适应多分辨率和自适应网格细化应用于可压缩Euler方程模拟的比较。SIAM J.科学。计算。,38(5):S173-S193,2016年·Zbl 1457.65069号
[19] G.Deslauriers和S.Dubuc。国际刑警组织,1987年·Zbl 0645.42010号
[20] G.Deslauriers和S.Dubuc。对称迭代插值过程。施工。约,5(1):491989年·Zbl 0659.65004号
[21] R.A.德沃尔。非线性近似。《数字学报》,7:51-1501998年·Zbl 0931.65007号
[22] E.Dilek、B.Erzincanli和M.Sahin。悬停果蝇近尾流结构的拉格朗日和欧拉相干结构的数值研究。理论和计算流体动力学,2019年4月。
[23] M.Domingues、S.Gomes、O.Roussel和K.Schneider。自适应多分辨率方法。ESAIM:会议记录,34:1-962011·Zbl 1302.65185号
[24] M.Domingues、S.M.Gomes和L.M.A.Díaz。块结构网格上的自适应小波表示和微分。申请。数字数学。,47:421-437, 2003. ·Zbl 1035.65167号
[25] D L多诺霍。插值小波变换。斯坦福大学统计系预印本,2(3):1-541992。
[26] D.L.Donoho和J.M.Johnstone。通过小波收缩实现理想的空间自适应。《生物特征》,81(3):425-4551994年·Zbl 0815.62019号
[27] M.Dreissigacker先生。分形树产生的湍流-PIV测量和与数值数据的比较。硕士论文,柏林大学,2017年。
[28] T.Engels、D.Kolomenskiy、K.Schneider、M.Farge、F.-O.Lehmann和J.Sesterhenn。大黄蜂扇动翅膀产生的螺旋涡。流体动力学。2018年第50:011419号决议。
[29] T.Engels、D.Kolomenskiy、K.Schneider、M.Farge、F.-O.Lehmann和J.Sesterhenn。系留飞行和自由飞行中飞行昆虫的湍流影响:高分辨率数值
[30] T.Engels等人/Commun。计算。物理。,30(2021),第1118-1149页实验。物理学。《流体版本》,4:0131032019年1月。
[31] T.Engels、D.Kolomenskiy、K.Schneider、F.-O.Lehmann和J.Sesterhenn。大黄蜂在剧烈颠簸中飞行。物理学。修订稿。,116:028103, 2016.
[32] T.Engels、D.Kolomenskiy、K.Schneider和J.Sesterhenn。使用具有体积惩罚的伪谱方法对颤振不稳定性进行二维模拟。《计算机与结构》,122:101-1122012。
[33] T.Engels、D.Kolomenskiy、K.Schneider和J.Sesterhenn。用体积惩罚法对流体-结构相互作用进行数值模拟。J.计算。物理。,281:96-115, 2015. ·Zbl 1351.76202号
[34] T.Engels、D.Kolomenskiy、K.Schneider和J.Sesterhenn。FluSI:一种新颖的用于扑翼昆虫飞行的并行模拟工具,使用具有体积惩罚的傅立叶方法。SIAM J.科学。计算。,38(5):S3-S242016·Zbl 1348.76106号
[35] M.Farge、K.Schneider和N.Kevlahan。使用自适应正交小波基对二维湍流进行非高斯和相干涡模拟。物理学。《流体》,11:2187-2201,1999年·Zbl 1147.76386号
[36] I.加甘蒂尼。表示四叉树的有效方法。ACM通讯,25(12):905-9101982·Zbl 0504.68057号
[37] I.加甘蒂尼。用于快速处理三维对象的线性八叉树。计算机图形和图像处理,20(4):365-3741982。
[38] A.Grossmann和J.Morlet。将Hardy函数分解为常形状的平方可积小波。SIAM J.数学。分析。,15:723-736, 1984. ·Zbl 0578.42007号
[39] J.Guermond和P.Minev。不可压缩Navier-Stokes方程的高阶时间步长。SIAM J.科学。计算。,37(6):A2656-A26812015·Zbl 1382.65315号
[40] A.Harten。离散多分辨率分析和广义小波。申请。数字。分析。,12:153-192, 1993. ·Zbl 0777.65004号
[41] A.Harten。双曲守恒律数值解的多分辨率算法。普通纯应用程序。数学。,48:1305-1342, 1995. ·Zbl 0860.65078号
[42] A.Harten。数据的多分辨率表示:一般框架。SIAM J.Numer。分析。,33:385-394, 1996.
[43] D.希尔伯特。在弗拉切斯特克,你会看到阿比尔顿·埃内尔·利尼。《数学年鉴》,38:459-4601891年。
[44] M.HolmströM。使用插值小波求解双曲偏微分方程。SIAM J.科学。计算。,21(2):405-420, 1999. ·Zbl 0959.65109号
[45] N.Kevlahan和M.Farge。二维湍流中的涡丝:产生、稳定性和效果。流体力学杂志。,346:49-76, 1997. ·兹比尔0910.76026
[46] N.Kevlahan和O.Vasilyev。高雷诺数下流体-结构相互作用的自适应小波配置方法。SIAM J.科学。计算。,26(6):1894-1915, 2005. ·Zbl 1149.76632号
[47] D.Kolomenskiy、J.-C.Nave和K.Schneider。具有多分辨率误差估计的自适应梯度增强水平集方法。科学杂志。计算。,66:116-140, 2016. ·Zbl 1338.35290号
[48] D.Kolomenskiy和K.Schneider。考虑移动固体障碍物体积惩罚的Navier-Stokes方程的傅里叶谱方法。J.计算。物理。,228:5687-5709, 2009. ·Zbl 1169.76045号
[49] J.Liandrat和P.Tchamitchian。使用空间小波近似求解一维正则Burgers方程。ICASE技术报告,1990年。
[50] S.Mallat公司。L2(R)的多分辨率逼近和正交小波基。事务处理。美国数学。Soc.,315:69-871989年·Zbl 0686.42018号
[51] S.Mallat公司。稀疏信号处理的小波之旅。学术出版社,第3版,T.Engels等人/Commun。计算。物理。,30(2021),第1118-1149页,第1149页,2009年。
[52] R.Mittal和G.Iacarino。浸没边界法。每年。流体力学版次。,37:239-261, 2005. ·Zbl 1117.76049号
[53] 穆勒。《守恒定律的自适应多尺度方案》,《计算科学与工程讲义》第27卷,施普林格,海德堡,2003年·Zbl 1016.76004号
[54] R.Nguyen van yen、D.Kolomenskiy和K.Schneider。通过体积惩罚对拉普拉斯和斯托克斯算子与Dirichlet边界条件的逼近:一个特殊的观点。数字。数学。,128:301-338, 2014. ·Zbl 1301.65118号
[55] T.Ohwada和P.Asinari。人工压缩性方法重新审视:不可压缩Navier-Stokes方程的渐近数值方法。J.公司。物理。,229:1698-1723, 2010. ·Zbl 1329.76063号
[56] C.S.Peskin公司。心脏血流的数值分析。J.计算。物理。,25:220-252, 1977. ·Zbl 0403.76100号
[57] C.S.Peskin公司。浸没边界法。《数字学报》,11:479-5172002年·Zbl 1123.74309号
[58] D.Rossinelli、C.Conti和P.Koumoutsakos。图形处理单元和多核中央处理器上的网格-粒子插值。《皇家学会哲学学报A:数学、物理和工程科学》,369(1944):2164-21752011年5月·Zbl 1223.68122号
[59] O.Roussel和K.Schneider。使用自适应多分辨率方法对弱可压缩涡旋混合层进行相干涡模拟。J.计算。物理。,229:2267-2286, 2010. ·Zbl 1303.76137号
[60] O.Roussel、K.Schneider、A.Tsigulin和H.Bockhorn。抛物线偏微分方程的保守全自适应多分辨率算法。J.计算。物理。,188:493-523, 2003. ·Zbl 1022.65093号
[61] K.施耐德。复杂几何中受限流体和等离子体湍流数值模拟的浸没边界法:综述。血浆物理学杂志。,81:435810601, 2015.
[62] K.Schneider和O.Vasilyev。计算流体动力学中的小波方法。每年。流体力学版次。,42:473-503, 2010. ·Zbl 1345.76085号
[63] F.Schornbaum和U.Rüde。超大规模块结构自适应网格细化。SIAM J.科学。计算。,40(3):C358-C3872018年·Zbl 06890193号
[64] 铃木(K.Suzuki)、南美(K.Minami)和稻村(T.Inamuro)。蝶形襟翼翼身模型产生的升力和推力:浸入式边界晶格玻尔兹曼模拟。流体力学杂志。,767:659-695, 2015.
[65] W.瑞典。提升方案:第二代小波的构造。SIAM J.数学。分析。,29(2):511-546, 1989. ·Zbl 2016年11月9日
[66] W.Sweldens和P.Schröder。在家里制作自己的小波。计算机图形学中的小波,1997年。
[67] W.Sweldens和P.Schröder。在家中构建自己的小波,第72-107页。施普林格-柏林-海德堡,柏林,海德堡·Zbl 0963.65153号
[68] C.K.W.Tam和J.C.Webb。计算声学中的色散关系保留有限差分格式。J.计算。物理。,107(2):262-281, 1993. ·Zbl 0790.76057号
[69] H.Truong、T.Engels、D.Kolomenskiy和K.Schneider。质量-弹簧流体-结构相互作用求解器:应用于柔性旋转机翼。计算和流体,200:1044262020·Zbl 1519.76116号
[70] G.Zumbusch。并行多级方法:自适应网格细化和负载平衡。数值数学进展。斯普林格,2003年·Zbl 1089.65127号
[71] T.Engels、D.Kolomenskiy、K.Schneider和J.Sesterhenn。FluSI:一种使用傅里叶方法和体积惩罚的新型昆虫扑翼飞行并行模拟工具。SIAM J.科学。计算。,38(5):S3-S242016·Zbl 1348.76106号
[72] A.George、L.C.Huang、W.-P.Tang和Y.D.Wu。曲线半交错网格上非定常可压缩流动(Re≤9500)的数值模拟。SIAM J.科学。计算。,21(6):2331-2351, 2000. ·兹伯利0989.76057
[73] N.Kevlahan和O.Vasilyev。高雷诺数下流体-结构相互作用的自适应小波配置方法。SIAM J.科学。计算。,26(6):1894-1915, 2005. ·Zbl 1149.76632号
[74] P.Koumoutsakos和A.Leonard。使用涡流方法对冲击启动圆柱体周围的流动进行高分辨率模拟。流体力学杂志。,296:1-38, 1995. ·Zbl 0849.76061号
[75] T.Ohwada和P.Asinari。重温人工压缩性方法:不可压缩Navier-Stokes方程的渐近数值方法。J.公司。物理。,229:1698-1723, 2010. ·Zbl 1329.76063号
[76] C.K.W.Tam和J.C.Webb。计算声学中的色散关系保留有限差分格式。J.计算。物理。,107(2):262-281, 1993. ·Zbl 0790.76057号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。