求解流固耦合问题的一种四步分裂有限元算法-应用数学和力学.PDFVIP

摇 , 34 7 Applied Mathematics and Mechanics摇 摇 2013年7月15 日出版 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 Vol.34,No.7,Jul.15,2013 文章编号:1000鄄0887(2013)07鄄0704鄄10 訫 应用数学和力学编委会,ISSN 1000鄄0887 求解流固耦合问题的一种 四步分裂有限元算法* 1 1,2 1 1 王华坤 洪国军 杨闻宇 喻国良 ,摇 ,摇 ,摇 (1.上海交通大学 船舶海洋与建筑工程学院,上海200240; 2.中交上海航道勘察设计研究院有限公司,上海200120) 摘要:摇 基于arbitrary Lagrangian Eulerian (ALE) 有限元方法,发展了一种求解流固耦合问题的弱 耦合算法.将半隐式四步分裂有限元格式推广至求解ALE描述下的Navier鄄Stokes(N鄄S)方程,并在 动量方程中引入迎风流线(streamline upwind/ Petrov鄄Galerkin,SUPG)稳定项以消除对流引发的速 度场数值振荡;采用Newmark鄄 法对结构方程进行时间离散;运用经典的Galerkin有限元法求解茁 修正的Laplace方程以实现网格更新,每个计算步施加网格总变形量防止结构长时间、大位移运动 时的网格质量恶化.运用上述算法对弹性支撑刚性圆柱体的流致振动问题进行了数值模拟,计算 结果与已有结果相吻合,初步验证了该算法的正确性和有效性. 关摇 键摇 词:摇 ALE有限元;摇 流固耦合;摇 半隐式四步分裂;摇 修正Laplace方程;摇 网格更新 中图分类号:摇 O357.1摇 摇 摇 文献标志码:摇 A DOI:10.3879/ j.issn.1000鄄0887.2013.07.005 引摇 摇 言 近年来,利用有限元方法求解流体与流固耦合问题是工程及科学计算领域的热点方向之 一.有限元方法不仅可以求解具有复杂几何边界流体的流动问题,而且在网格较粗的情况下, 还能通过采用高阶插值函数来提高计算精度.但是,由于不可压缩流体N鄄S方程的混合特性, 当采用不恰当的插值函数时会导致压力场的伪振荡.为此,国内外学者提出了多种稳定化有限 元格式,例如,为避免不可压缩流动 Babu觢ka鄄Brezzi(B鄄B)条件而提出的Galerkin 最小二乘 [1] [2] [3] 法 ,压力稳定化Petrov鄄Galerkin格式 以及有限增量法 等.分裂步骤法也是其中之一,该 方法通过算子分裂解耦求解流场速度和压力,被认为是一种高效、稳定的流体方程数值解法. [4] [5] 分裂法最早是Chorin 在有限差分法中针对不可压缩流动问题提出,随后由Donea等 引入 到有限元公式中,经过多年的发展,已经在众多流动问题得到了广泛应用. 对于分裂步骤法,若在分裂步中完全不考虑压力项,则该算法具有时间方向一阶精度;而 将压力梯度项保留于分裂步则时间方向精度可达到二阶,但会失去压力场稳定性.二阶精度分 裂步骤法成功应用于有限差分法中,但因面临压力场伪振荡问题,直接推广到有限元方法会遇 [6] 到一些很困难.最近Choi等 提出了一种具有二阶时间精度并能够保证压力稳定的四步分裂 * 收稿日期:摇 2013鄄03鄄27;修订日期:摇 2013鄄05鄄28 基金项目:摇 国家高技术研究发展计划项目(2007AA09Z312) 作者简介:摇 王华坤(1985—),男,河北人,博士生 (E鄄mail:hkwang@);