Anisotropic mesh adaptation for 3D anisotropic diffusion problems with application to fractured reservoir simulation

Xianping Li, Weizhang Huang

Research output: Contribution to journalArticlepeer-review

4 Scopus citations

Abstract

Anisotropic mesh adaptation is studied for linear finite element solution of 3D anisotropic diffusion problems. The M-uniform mesh approach is used, where an anisotropic adaptive mesh is generated as a uniform one in the metric specified by a tensor. In addition to mesh adaptation, preservation of the maximum principle is also studied. Some new sufficient conditions for maximum principle preservation are developed, and a mesh quality measure is defined to server as a good indicator. Four different metric tensors are investigated: one is the identity matrix, one focuses on minimizing an error bound, another one on preservation of the maximum principle, while the fourth combines both. Numerical examples show that these metric tensors serve their purposes. Particularly, the fourth leads to meshes that improve the satisfaction of the maximum principle by the finite element solution while concentrating elements in regions where the error is large. Application of the anisotropic mesh adaptation to fractured reservoir simulation in petroleum engineering is also investigated, where unphysical solutions can occur and mesh adaptation can help improving the satisfaction of the maximum principle.

Original languageEnglish (US)
Pages (from-to)913-940
Number of pages28
JournalNumerical Mathematics
Volume10
Issue number4
DOIs
StatePublished - Nov 1 2017
Externally publishedYes

Keywords

  • Anisotropic diffusion
  • Anisotropic mesh adaptation
  • Discrete maximum principle
  • Finite element method
  • Petroleum engineering
  • Three dimensional

ASJC Scopus subject areas

  • Modeling and Simulation
  • Control and Optimization
  • Computational Mathematics
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Anisotropic mesh adaptation for 3D anisotropic diffusion problems with application to fractured reservoir simulation'. Together they form a unique fingerprint.

Cite this