一种适用于浸入有限元方法的网格自适应方法

doi:10.13700/j.bh.1001-5965.2019.0269

北京航空航天大学学报 ›› 2020, Vol. 46 ›› Issue (3): 588-597.doi: 10.13700/j.bh.1001-5965.2019.0269

一种适用于浸入有限元方法的网格自适应方法

张华¹, 白俊强^1,2, 乔磊¹, 刘艳³

1. 西北工业大学航空学院, 西安 710072;
2. 西北工业大学无人系统技术研究院, 西安 710072;
3. 航空工业沈阳飞机设计研究所, 沈阳 110035

收稿日期:2019-05-31 发布日期:2020-03-28
通讯作者: 白俊强 E-mail:junqiang@nwpu.edu.cn
作者简介:张华,女,硕士研究生。主要研究方向:浸入类算法、有限元方法、网格自适应技术;白俊强,男,博士,教授,博士生导师。主要研究方向:飞行器设计、计算流体力学;乔磊,男,博士,助理研究员。主要研究方向:计算流体力学;刘艳,女,博士,高级工程师。主要研究方向:气动弹性算法。
基金资助:
国家自然科学基金（11802245，11702284）

An adaptive mesh refinement method based on immersed finite element method

ZHANG Hua¹, BAI Junqiang^1,2, QIAO Lei¹, LIU Yan³

1. School of Aeronautics, Northwestern Polytechnical University, Xi'an 710072, China;
2. Unmanned System Research Institute, Northwestern Polytechnical University, Xi'an 710072, China;
3. AVIC Shenyang Aircraft Research Institute, Shenyang 110035, China

Received:2019-05-31 Published:2020-03-28
Supported by:
National Natural Science Foundation of China (11802245,11702284)

摘要/Abstract

摘要： 针对动边界流固耦合的数值模拟问题，基于浸入有限元方法提出了一种耦合流场特征和几何特征的笛卡儿网格局部加密自适应方法，克服了单个自适应指示因子无法精确捕捉固体运动的特征的不足。在耦合自适应策略中，分别以流场涡量和固体位置作为流场和几何信息指示因子来驱动网格自适应。通过方腔顶盖驱动圆盘流动算例，以圆盘体积守恒和特征点的运动轨迹验证耦合自适应方法的优势。计算结果表明：仅基于流动特征的自适应不能很好地保证圆盘的体积守恒；仅基于几何特征的自适应无法有效追踪圆盘的轨迹；而耦合自适应策略能同时较好地保证两项指标的计算精度，在保证总体计算自由度不变的情况下，圆盘区域速度散度2-范数降低了一个数量级，圆盘的轨迹误差2-范数降低了2个数量级。

关键词: 网格自适应, 浸入有限元方法, 浸入边界法, 流固耦合, 笛卡儿网格

Abstract: For the numerical simulation of fluid-structure interaction with moving boundary, a local Cartesian mesh adaptation method coupling flow field features and geometric features is developed based on immersed finite element method. This method overcomes the inaccuracy of simulating solid motion with a single adaptive indicator. In the coupling adaptation, the vorticity is used as the adaptive indicator factor for flow field, and the solid position is used as the indicator for the geometric feature to drive mesh adaptation. The advantages of the coupling adaptive strategy are verified by a numerical example, disk entrained in a lid-driven cavity flow, with volume conservation of the disk and some points' motion trajectory on disks. The computational results show that the volume conservation of the disk cannot be well guaranteed only by the adaptation based on flow characteristics; the trajectory tracking of the disk cannot be effectively achieved only by the geometry-based adaptation; but the coupling adaptation strategy in this paper can ensure the accuracy of the two indexes at the same time. When the overall computational degrees of freedom remain constant, the 2-norm of divergence of velocity can be reduced by one order of magnitude and the trajectory error 2-norm of the disk is reduced by two orders of magnitude.

Key words: mesh adaptation, immersed finite element method, immersed boundary method, fluid-structure interaction, Cartesian mesh

中图分类号:

张华, 白俊强, 乔磊, 刘艳. 一种适用于浸入有限元方法的网格自适应方法[J]. 北京航空航天大学学报, 2020, 46(3): 588-597.

ZHANG Hua, BAI Junqiang, QIAO Lei, LIU Yan. An adaptive mesh refinement method based on immersed finite element method[J]. Journal of Beijing University of Aeronautics and Astronautics, 2020, 46(3): 588-597.

参考文献

[1] PESKIN C.Flow patterns around heart valves:A numerical method[J].Journal of Computational Physics,1972,10(2):252-271.
[2] BORAZJANI I,GE L,SOTIROPOULOS F.High-resolution fluid-structure interaction simulations of flow through a bi-leaflet mechanical heart valve in an anatomic aorta[J].Annual Review of Biomedical Engineering,2010,38:326-344.
[3] SOTIROPOULOS F,BORAZJANI I.A review of state-of-the-art numerical methods for simulating flow through mechanical heart valves[J].Medical & Biological Engineering & Computing,2009,47(3):245-256.
[4] ZHU L,PESKIN C S.Simulation of a flapping flexible filament in a flowing soap film by the immersed boundary method[J].Journal of Computational Physics,2002,179(2):452-468.
[5] ZHANG J,CHILDRESS S,LIBCHABER A,et al.Flexible filaments in a flowing soap film as a model for one-dimensional flags in a two-dimensional wind[J].Nature,2000,408(6814):835-839.
[6] MITTAL R,SESHADRI V,UDAYKUMAR H S.Flutter,tumble and vortex induced autorotation[J].Theoretical and Computational Fluid Dynamics,2004,17(3):165-170.
[7] GOLDSTEIN D,HANDLER R,SIROVICH L.Modeling a no-slip flow boundary with an external force field[J].Journal of Computational Physics,1993,105(2):354-366.
[8] GHIAS R,MITTAL R,DONG H.A sharp interface immersed boundary method for compressible viscous flows[J].Journal of Computational Physics,2007,225(1):528-533.
[9] MITTAL R,DONG H.A versatile sharp interface immersed boundary method for incompressible flows with complex boundaries[J].Journal of Computational Physics,2008,227(10):4825-4852.
[10] 秦望龙,吕宏强,伍贻兆.基于混合网格的高阶间断有限元黏流数值解法[J].力学学报,2013,45(6):987-991.QIN W L,LV H Q,WU Y Z.High-order discontinuous Galerkin solution of N-S equations on hybrid mesh[J].Chinese Journal of Theoretical and Applied Mechanics,2013,45(6):987-991(in Chinese).
[11] 于剑,阎超.Navier-Stokes方程间断Galerkin有限元方法研究[J].力学学报,2010,42(5):962-970.YU J,YAN C.Study on discontinuous Galerkin method for Navier-Stokes equations[J].Chinese Journal of Theoretical and Applied Mechanics,2010,42(5):962-970(in Chinese).
[12] BOFFI D,GASTALDI L.A finite element approach for the immersed boundary method[J].Computers and Structures,2003,81(8-11):491-501.
[13] WANG X,LIU W K.Extended immersed boundary method using FEM and RKPM[J].Computer Methods in Applied Mechanics and Engineering,2004,193(12-14):1305-1321.
[14] ZHANG L,GERSTENBERGER A,WANG X,et al.Immersed finite element method[J].Computer Methods in Applied Mechanics and Engineering,2004,193(21-22):2051-2067.
[15] HELTAI L,FRANCESCO C.Variational implementation of immersed finite element methods[J].Computer Methods in Applied Mechanics and Engineering,2012,229-232:110-127.
[16] ROY S,HELTAI L,COSTANZO F.Benchmarking the immersed finite element method for fluid-structure interaction problems[J].Computers and Mathematics with Applications,2015,69(10):1167-1188.
[17] 吴泽艳,王立峰,武哲.基于简单WENO-间断Galerkin的Euler方程自适应计算[J].北京航空航天大学学报,2016, 42(4):806-814.WU Z Y,WANG L F,WU Z.Adaptive simple WENO limiter-discontinuous Galerkin method for Euler equations[J].Journal of Beijing University of Aeronautics and Astronautics, 2016,42(4):806-814(in Chinese).
[18] BERGER M J,OLIGER J.Adaptive mesh refinement for hyperbolic partial differential equations[J].Journal of Computational Physics,1984,53(3):484-512.
[19] ABGRALL R,BEAUGENDRE H,DOBRZYNSKI C.An immersed boundary method using unstructured anisotropic mesh adaptation combined with level-sets and penalization techniques[J].Journal of Computational Physics,2014,257:83-101.
[20] PENG B,ZHOU C H.An approach of dynamic mesh adaptation for simulating 3-dimensional unsteady moving-immersed-boundary flows[J].International Journal for Numerical Methods in Fluids,2018,87(4):180-201.
[21] XU L C,TIAN F B,YOUNG J,et al.A novel geometry-adaptive Cartesian grid based immersed boundary-lattice Boltzmann method for fluid-structure interactions at moderate and high Reynolds numbers[J].Journal of Computational Physics,2018,375:22-56.

一种适用于浸入有限元方法的网格自适应方法

An adaptive mesh refinement method based on immersed finite element method

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	宣传伟, 韩景龙. 重叠网格中隐式装配策略的改进[J]. 北京航空航天大学学报, 2020, 46(2): 350-358.
[2]	郝帅, 马铁林, 甘文彪, 李根. 氢动力无人机大展弦比机翼静气弹特性分析[J]. 北京航空航天大学学报, 2017, 43(8): 1670-1676.
[3]	倪奇峰, 侯安平, 刘若阳, 周拜豪, 张明明. 大攻角状态压气机分离流及叶片动力响应特性[J]. 北京航空航天大学学报, 2017, 43(7): 1410-1418.
[4]	孟繁超, 董素君, 江泓升, 王浚. 长时间流固耦合传热过程的快速算法[J]. 北京航空航天大学学报, 2017, 43(6): 1224-1230.
[5]	刘森, 张怀亮, 彭欢. 基础振动诱发的流体-管道轴向耦合振动特性[J]. 北京航空航天大学学报, 2016, 42(3): 610-618.
[6]	杨慧, 李振鹏. 转静干涉对转子叶片颤振特性的影响[J]. 北京航空航天大学学报, 2016, 42(2): 258-264.
[7]	于志强, 佟浩, 李勇, 柴永生, 王志强, 李朝将. 电火花穿孔加工中管状电极振动幅值特性[J]. 北京航空航天大学学报, 2016, 42(11): 2516-2523.
[8]	杨永强, 马云鹏, 武哲. 基于流固耦合的充气翼内压对翼面变形影响分析[J]. 北京航空航天大学学报, 2014, 40(2): 188-192.
[9]	张春, 曹义华. 基于弱耦合的翼伞气动变形数值模拟[J]. 北京航空航天大学学报, 2013, 39(5): 605-609.
[10]	郑赟, 杨慧. 跨音速风扇全环叶片颤振特性的流固耦合分析[J]. 北京航空航天大学学报, 2013, 39(5): 626-630.
[11]	张小伟, 王延荣, 谢胜百, 单鹏. 弹性整流罩分离的流固耦合仿真方法[J]. 北京航空航天大学学报, 2009, 35(8): 976-979.
[12]	宋乾福, 曹义华, 蒋崇文. 半球形伞流固耦合工程算法[J]. 北京航空航天大学学报, 2009, 35(1): 96-99.
[13]	潘星,曹义华. 降落伞开伞过程的多结点模型仿真[J]. 北京航空航天大学学报, 2008, 34(02): 188-192.
[14]	陈卫卫,孙凌玉,张广越,李飞. 流体驱动高速机构中缓冲装置的优化设计[J]. 北京航空航天大学学报, 2008, 34(02): 206-209.
[15]	王侃, 曹义华, 于子文, 潘星. 降落伞流固耦合问题的数值模拟和流场分析[J]. 北京航空航天大学学报, 2007, 33(09): 1029-1032.