在线咨询
中国工业与应用数学学会会刊
主管:中华人民共和国教育部
主办:西安交通大学
ISSN 1005-3085  CN 61-1269/O1

工程数学学报 ›› 2016, Vol. 33 ›› Issue (1): 52-62.doi: 10.3969/j.issn.1005-3085.2016.01.006

• • 上一篇    下一篇

粘弹流动的间断有限元模拟

郭虹平,   欧阳洁,  杨广辉,  周  文   

  1. 西北工业大学应用数学系,西安 710129
  • 收稿日期:2013-11-25 接受日期:2014-05-19 出版日期:2016-02-15 发布日期:2016-04-15
  • 通讯作者: 欧阳洁 E-mail: jieouyang@nwpu.edu.cn
  • 基金资助:
    国家重点基础研究发展计划 (2012CB025903).

Numerical Simulation of Viscoelastic Flows Using Discontinuous Galerkin Finite Element Method

GUO Hong-ping,  OUYANG Jie,  YANG Guang-hui,  ZHOU Wen   

  1. Department of Applied Mathematics, Northwestern Polytechnical University, Xi'an 710129
  • Received:2013-11-25 Accepted:2014-05-19 Online:2016-02-15 Published:2016-04-15
  • Contact: J.Ouyang. E-mail address: jieouyang@nwpu.edu.cn
  • Supported by:
    The State Key Development Program for Basic Research of China (2012CB025903).

摘要: 针对传统有限元法求解Oldroyd-B本构方程时需加入稳定化方案的缺点,本文基于非结构网格给出了统一间断有限元求解框架.该框架包含采用IPDG(interior penalty discontinuous Galerkin)求解质量方程和动量方程,与采用RKDG(Runge-Kutta discontinuous Galerkin)求解本构方程这两个核心.数值结果表明:该方法在求解Oldroyd-B本构方程时无需加入稳定化方案,实施比有限元法简便,且具有较高的计算精度,可有效地模拟包含应力奇异点的复杂粘弹流动问题,进而揭示非牛顿粘弹流动的基本特征.

关键词: 间断有限元, 非结构, 粘弹流体

Abstract:

The traditional finite element method needs to supplement a stabilization scheme to simulate Oldroyd-B viscoelastic flows. To alleviate this issue, a unified discontinuous Galerkin finite element framework based on unstructured grids is proposed in this paper. The system contains two key points: one is using the IPDG (interior penalty discontinuous Galerkin) method to discretize mass and momen-tum equations, and the other is employing the RKDG (Runge-Kutta DG) method to solve the Oldroyd-B constitutive equation. Simulation results reveal the intrinsic characteristics of non-Newtonian viscoelastic fluids and indicate that the approach can effectively overcome the drawback of the traditional finite element method, which redundantly introduces stabilization process in the method. Moreover, these results substantiate that the proposed method is simple to implement, has high accuracy and can be used to simulate complex viscoelastic flows with stress singularity.

Key words: discontinuous Galerkin, unstructured grids, viscoelastic flow

中图分类号: