An Efficient ELLAM Implementation for Modeling Solute Transport in Fractured Porous Media

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

The goal of this study is to introduce an adaptation of the Eulerian-Lagrangian localized adjoint method (ELLAM) for the simulation of mass transport in fractured porous media, and to evaluate the performance of ELLAM in such a case. The fractures are represented explicitly using the discrete fracture model. The velocity field was calculated using the mixed hybrid finite element method. A sound ELLAM implementation is developed to address numerical artifacts (spurious oscillations and numerical dispersion) arising from the presence of fractures. The efficiency of the developed ELLAM implementation was further improved by taking advantage of the parallel computing on shared memory architecture for the tasks related to particles tracking and linear system resolving. The performance of ELLAM was tested by comparing its results with the Eulerian discontinuous Galerkin method based on several benchmark problems dealing with different fracture configurations. The results highlight the robustness and accuracy of ELLAM, as it allows the use of large time steps, and overcomes the Courant-Friedrichs-Lewy (CFL) restriction. The outcome also shows that ELLAM is more efficient when fracture density is increased.

Original languageEnglish
Article number46
JournalWater, Air, and Soil Pollution
Volume229
Issue number2
DOIs
Publication statusPublished - 1 Feb 2018

Fingerprint

adjoint method
Solute transport
solute transport
Porous materials
porous medium
modeling
Memory architecture
Galerkin methods
Parallel processing systems
parallel computing
Galerkin method
Linear systems
Mass transfer
mass transport
Acoustic waves
Finite element method
finite element method
artifact
oscillation

Keywords

  • Discontinuous Galerkin method
  • Discrete fracture model
  • Eulerian-Lagrangian method
  • Fractured porous media
  • Solute transport

ASJC Scopus subject areas

  • Environmental Engineering
  • Environmental Chemistry
  • Ecological Modelling
  • Water Science and Technology
  • Pollution

Cite this

@article{d497066d19d44e79a35bea7eb96041a3,
title = "An Efficient ELLAM Implementation for Modeling Solute Transport in Fractured Porous Media",
abstract = "The goal of this study is to introduce an adaptation of the Eulerian-Lagrangian localized adjoint method (ELLAM) for the simulation of mass transport in fractured porous media, and to evaluate the performance of ELLAM in such a case. The fractures are represented explicitly using the discrete fracture model. The velocity field was calculated using the mixed hybrid finite element method. A sound ELLAM implementation is developed to address numerical artifacts (spurious oscillations and numerical dispersion) arising from the presence of fractures. The efficiency of the developed ELLAM implementation was further improved by taking advantage of the parallel computing on shared memory architecture for the tasks related to particles tracking and linear system resolving. The performance of ELLAM was tested by comparing its results with the Eulerian discontinuous Galerkin method based on several benchmark problems dealing with different fracture configurations. The results highlight the robustness and accuracy of ELLAM, as it allows the use of large time steps, and overcomes the Courant-Friedrichs-Lewy (CFL) restriction. The outcome also shows that ELLAM is more efficient when fracture density is increased.",
keywords = "Discontinuous Galerkin method, Discrete fracture model, Eulerian-Lagrangian method, Fractured porous media, Solute transport",
author = "Fanilo Ramasomanana and Marwan Fahs and Husam Baalousha and Nicolas Barth and Said Ahzi",
year = "2018",
month = "2",
day = "1",
doi = "10.1007/s11270-018-3690-8",
language = "English",
volume = "229",
journal = "Water, Air, and Soil Pollution",
issn = "0049-6979",
publisher = "Springer Netherlands",
number = "2",

}

TY - JOUR

T1 - An Efficient ELLAM Implementation for Modeling Solute Transport in Fractured Porous Media

AU - Ramasomanana, Fanilo

AU - Fahs, Marwan

AU - Baalousha, Husam

AU - Barth, Nicolas

AU - Ahzi, Said

PY - 2018/2/1

Y1 - 2018/2/1

N2 - The goal of this study is to introduce an adaptation of the Eulerian-Lagrangian localized adjoint method (ELLAM) for the simulation of mass transport in fractured porous media, and to evaluate the performance of ELLAM in such a case. The fractures are represented explicitly using the discrete fracture model. The velocity field was calculated using the mixed hybrid finite element method. A sound ELLAM implementation is developed to address numerical artifacts (spurious oscillations and numerical dispersion) arising from the presence of fractures. The efficiency of the developed ELLAM implementation was further improved by taking advantage of the parallel computing on shared memory architecture for the tasks related to particles tracking and linear system resolving. The performance of ELLAM was tested by comparing its results with the Eulerian discontinuous Galerkin method based on several benchmark problems dealing with different fracture configurations. The results highlight the robustness and accuracy of ELLAM, as it allows the use of large time steps, and overcomes the Courant-Friedrichs-Lewy (CFL) restriction. The outcome also shows that ELLAM is more efficient when fracture density is increased.

AB - The goal of this study is to introduce an adaptation of the Eulerian-Lagrangian localized adjoint method (ELLAM) for the simulation of mass transport in fractured porous media, and to evaluate the performance of ELLAM in such a case. The fractures are represented explicitly using the discrete fracture model. The velocity field was calculated using the mixed hybrid finite element method. A sound ELLAM implementation is developed to address numerical artifacts (spurious oscillations and numerical dispersion) arising from the presence of fractures. The efficiency of the developed ELLAM implementation was further improved by taking advantage of the parallel computing on shared memory architecture for the tasks related to particles tracking and linear system resolving. The performance of ELLAM was tested by comparing its results with the Eulerian discontinuous Galerkin method based on several benchmark problems dealing with different fracture configurations. The results highlight the robustness and accuracy of ELLAM, as it allows the use of large time steps, and overcomes the Courant-Friedrichs-Lewy (CFL) restriction. The outcome also shows that ELLAM is more efficient when fracture density is increased.

KW - Discontinuous Galerkin method

KW - Discrete fracture model

KW - Eulerian-Lagrangian method

KW - Fractured porous media

KW - Solute transport

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

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

U2 - 10.1007/s11270-018-3690-8

DO - 10.1007/s11270-018-3690-8

M3 - Article

VL - 229

JO - Water, Air, and Soil Pollution

JF - Water, Air, and Soil Pollution

SN - 0049-6979

IS - 2

M1 - 46

ER -