An efficient lumped mixed hybrid finite element formulation for variably saturated groundwater flow

Benjamin Belfort, Fanilo Ramasomanana, Anis Younes, François Lehmann

Research output: Contribution to journalArticle

12 Citations (Scopus)

Abstract

Accurate numerical simulation of infi ltration in the vadose zone remains a challenge, especially when very sharp fronts are modeled. In this work, we used the mixed hybrid fi nite element (MHFE) method, which allows a simultaneous approximation of both pressure head and velocity and can handle general irregular grids with highly heterogeneous permeability. For many problems dealing with unsaturated water fl ow, however, the MHFE solutions exhibit signifi cant unphysical oscillations. To avoid this phenomenon, we developed an effi cient mass-lumping scheme with the MHFE method for solving the mixed form of the Richards equation. In this work, the standard and lumped MHFE formulations were compared and the ability of the lumped formulation to reduce unphysical oscillations was determined for oneand two-dimensional infi ltration problems. Theoretical analysis based on the M-matrix property, which guarantees the discrete maximum principle, and practical test cases performed in this study underline the advantage of using an acute triangulation to completely remove the unphysical oscillations. Indeed, contrary to the standard approach, the lumped formulation satisfi es the M-matrix property without any constraint on the time step size to be used.

Original languageEnglish
Pages (from-to)352-362
Number of pages11
JournalVadose Zone Journal
Volume8
Issue number2
DOIs
Publication statusPublished - May 2009
Externally publishedYes

Fingerprint

groundwater flow
oscillation
Richards' equation
Richards equation
matrix
vadose zone
triangulation
permeability
case studies
methodology
simulation
water
testing
method

ASJC Scopus subject areas

  • Soil Science

Cite this

An efficient lumped mixed hybrid finite element formulation for variably saturated groundwater flow. / Belfort, Benjamin; Ramasomanana, Fanilo; Younes, Anis; Lehmann, François.

In: Vadose Zone Journal, Vol. 8, No. 2, 05.2009, p. 352-362.

Research output: Contribution to journalArticle

Belfort, Benjamin ; Ramasomanana, Fanilo ; Younes, Anis ; Lehmann, François. / An efficient lumped mixed hybrid finite element formulation for variably saturated groundwater flow. In: Vadose Zone Journal. 2009 ; Vol. 8, No. 2. pp. 352-362.
@article{b88267cecbb442ee9fd1dd3ac22eeb56,
title = "An efficient lumped mixed hybrid finite element formulation for variably saturated groundwater flow",
abstract = "Accurate numerical simulation of infi ltration in the vadose zone remains a challenge, especially when very sharp fronts are modeled. In this work, we used the mixed hybrid fi nite element (MHFE) method, which allows a simultaneous approximation of both pressure head and velocity and can handle general irregular grids with highly heterogeneous permeability. For many problems dealing with unsaturated water fl ow, however, the MHFE solutions exhibit signifi cant unphysical oscillations. To avoid this phenomenon, we developed an effi cient mass-lumping scheme with the MHFE method for solving the mixed form of the Richards equation. In this work, the standard and lumped MHFE formulations were compared and the ability of the lumped formulation to reduce unphysical oscillations was determined for oneand two-dimensional infi ltration problems. Theoretical analysis based on the M-matrix property, which guarantees the discrete maximum principle, and practical test cases performed in this study underline the advantage of using an acute triangulation to completely remove the unphysical oscillations. Indeed, contrary to the standard approach, the lumped formulation satisfi es the M-matrix property without any constraint on the time step size to be used.",
author = "Benjamin Belfort and Fanilo Ramasomanana and Anis Younes and Fran{\cc}ois Lehmann",
year = "2009",
month = "5",
doi = "10.2136/vzj2008.0108",
language = "English",
volume = "8",
pages = "352--362",
journal = "Vadose Zone Journal",
issn = "1539-1663",
publisher = "Soil Science Society of America",
number = "2",

}

TY - JOUR

T1 - An efficient lumped mixed hybrid finite element formulation for variably saturated groundwater flow

AU - Belfort, Benjamin

AU - Ramasomanana, Fanilo

AU - Younes, Anis

AU - Lehmann, François

PY - 2009/5

Y1 - 2009/5

N2 - Accurate numerical simulation of infi ltration in the vadose zone remains a challenge, especially when very sharp fronts are modeled. In this work, we used the mixed hybrid fi nite element (MHFE) method, which allows a simultaneous approximation of both pressure head and velocity and can handle general irregular grids with highly heterogeneous permeability. For many problems dealing with unsaturated water fl ow, however, the MHFE solutions exhibit signifi cant unphysical oscillations. To avoid this phenomenon, we developed an effi cient mass-lumping scheme with the MHFE method for solving the mixed form of the Richards equation. In this work, the standard and lumped MHFE formulations were compared and the ability of the lumped formulation to reduce unphysical oscillations was determined for oneand two-dimensional infi ltration problems. Theoretical analysis based on the M-matrix property, which guarantees the discrete maximum principle, and practical test cases performed in this study underline the advantage of using an acute triangulation to completely remove the unphysical oscillations. Indeed, contrary to the standard approach, the lumped formulation satisfi es the M-matrix property without any constraint on the time step size to be used.

AB - Accurate numerical simulation of infi ltration in the vadose zone remains a challenge, especially when very sharp fronts are modeled. In this work, we used the mixed hybrid fi nite element (MHFE) method, which allows a simultaneous approximation of both pressure head and velocity and can handle general irregular grids with highly heterogeneous permeability. For many problems dealing with unsaturated water fl ow, however, the MHFE solutions exhibit signifi cant unphysical oscillations. To avoid this phenomenon, we developed an effi cient mass-lumping scheme with the MHFE method for solving the mixed form of the Richards equation. In this work, the standard and lumped MHFE formulations were compared and the ability of the lumped formulation to reduce unphysical oscillations was determined for oneand two-dimensional infi ltration problems. Theoretical analysis based on the M-matrix property, which guarantees the discrete maximum principle, and practical test cases performed in this study underline the advantage of using an acute triangulation to completely remove the unphysical oscillations. Indeed, contrary to the standard approach, the lumped formulation satisfi es the M-matrix property without any constraint on the time step size to be used.

UR - http://www.scopus.com/inward/record.url?scp=67149108357&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=67149108357&partnerID=8YFLogxK

U2 - 10.2136/vzj2008.0108

DO - 10.2136/vzj2008.0108

M3 - Article

VL - 8

SP - 352

EP - 362

JO - Vadose Zone Journal

JF - Vadose Zone Journal

SN - 1539-1663

IS - 2

ER -