数学物理学报(英文版) ›› 2023, Vol. 43 ›› Issue (5): 2026-2042.doi: 10.1007/s10473-023-0506-8

• • 上一篇    下一篇

A MIXED FINITE ELEMENT AND CHARACTERISTIC MIXED FINITE ELEMENT FOR INCOMPRESSIBLE MISCIBLE DARCY-FORCHHEIMER DISPLACEMENT AND NUMERICAL ANALYSIS*

Yirang Yuan1, Changfeng Li2,†, Tongjun Sun1, Qing Yang3   

  1. 1. Institute of Mathematics, Shandong University, Jinan 250100, China;
    2. School of Economics, Shandong University, Jinan 250100, China;
    3. School of Mathematical Sciences, Shandong Normal University, Jinan 250014, China
  • 收稿日期:2022-02-14 修回日期:2023-05-05 发布日期:2023-10-25

A MIXED FINITE ELEMENT AND CHARACTERISTIC MIXED FINITE ELEMENT FOR INCOMPRESSIBLE MISCIBLE DARCY-FORCHHEIMER DISPLACEMENT AND NUMERICAL ANALYSIS*

Yirang Yuan1, Changfeng Li2,†, Tongjun Sun1, Qing Yang3   

  1. 1. Institute of Mathematics, Shandong University, Jinan 250100, China;
    2. School of Economics, Shandong University, Jinan 250100, China;
    3. School of Mathematical Sciences, Shandong Normal University, Jinan 250014, China
  • Received:2022-02-14 Revised:2023-05-05 Published:2023-10-25
  • Contact: †Changfeng Li, E-mail:cfli@sdu.edu.cn
  • About author:Yirang Yuan, E-mail: yryuan@sdu.edu.cn; Tongjun Sun, E-mail:tjsun@sdu.edu.cn; Qing Yang, E-mail: sd_yangq@163.com
  • Supported by:
    Natural Science Foundation of Shandong Province (ZR2021MA019).

摘要: In this paper a mixed finite element-characteristic mixed finite element method is discussed to simulate an incompressible miscible Darcy-Forchheimer problem. The flow equation is solved by a mixed finite element and the approximation accuracy of Darch-Forchheimer velocity is improved one order. The concentration equation is solved by the method of mixed finite element, where the convection is discretized along the characteristic direction and the diffusion is discretized by the zero-order mixed finite element method. The characteristics can confirm strong stability at sharp fronts and avoids numerical dispersion and nonphysical oscillation. In actual computations the characteristics adopts a large time step without any loss of accuracy. The scalar unknowns and its adjoint vector function are obtained simultaneously and the law of mass conservation holds in every element by the zero-order mixed finite element discretization of diffusion flux. In order to derive the optimal $3/2$-order error estimate in $L^2$ norm, a post-processing technique is included in the approximation to the scalar unknowns. Numerical experiments are illustrated finally to validate theoretical analysis and efficiency. This method can be used to solve such an important problem.

关键词: Darcy-Forchheimer miscible displacement, mixed element-characteristic mixed element-postprocessing scheme, local conservation of mass, $3/2$-order error estimates in $L^2$-norm, numerical computation

Abstract: In this paper a mixed finite element-characteristic mixed finite element method is discussed to simulate an incompressible miscible Darcy-Forchheimer problem. The flow equation is solved by a mixed finite element and the approximation accuracy of Darch-Forchheimer velocity is improved one order. The concentration equation is solved by the method of mixed finite element, where the convection is discretized along the characteristic direction and the diffusion is discretized by the zero-order mixed finite element method. The characteristics can confirm strong stability at sharp fronts and avoids numerical dispersion and nonphysical oscillation. In actual computations the characteristics adopts a large time step without any loss of accuracy. The scalar unknowns and its adjoint vector function are obtained simultaneously and the law of mass conservation holds in every element by the zero-order mixed finite element discretization of diffusion flux. In order to derive the optimal $3/2$-order error estimate in $L^2$ norm, a post-processing technique is included in the approximation to the scalar unknowns. Numerical experiments are illustrated finally to validate theoretical analysis and efficiency. This method can be used to solve such an important problem.

Key words: Darcy-Forchheimer miscible displacement, mixed element-characteristic mixed element-postprocessing scheme, local conservation of mass, $3/2$-order error estimates in $L^2$-norm, numerical computation

中图分类号: 

  • 65M12