A {Scharfetter} {Gummel} based {Geothermal} model

Mustapha Zakari. ( 2025 )
in: 2025 {RING} meeting, pages 186--196, ASGA

Abstract

The numerical modeling of subsurface geothermal systems requires strong coupling of porous flow hydraulic, thermal, chemical and porous rock mechanical processes. Numerical models can give first insights on the characteristics of these processes in subsurface models hence helping to investigate the feasibility and efficiency of different geothermal strategies for different subsurface configurations. We propose to develop a numerical model of reasonable accuracy and complexity. The numerical model takes into account porous flow hydraulic and thermal processes. In order to avoid strong time steps limitations, stationary and fully implicit models are investigated. The discretization of the advection diffusion part of the temperature equation uses the Schaerfetter Gummel Method (SGM). Initially developed for drift diffusion density models arising in semi-conductor and plasma discharges modeling, SGM gives a robust and low stencil numerical scheme that allows to automatically switch from a centered scheme to an upwind one by considering advection-diffusion ratio. We present our SGM formulation and show first validation results.

Download / Links

BibTeX Reference

@inproceedings{Zakari2025RM,
 abstract = {The numerical modeling of subsurface geothermal systems requires strong coupling of porous flow hydraulic, thermal, chemical and porous rock mechanical processes. Numerical models can give first insights on the characteristics of these processes in subsurface models hence helping to investigate the feasibility and efficiency of different geothermal strategies for different subsurface configurations. We propose to develop a numerical model of reasonable accuracy and complexity. The numerical model takes into account porous flow hydraulic and thermal processes. In order to avoid strong time steps limitations, stationary and fully implicit models are investigated. The discretization of the advection diffusion part of the temperature equation uses the Schaerfetter Gummel Method (SGM). Initially developed for drift diffusion density models arising in semi-conductor and plasma discharges modeling, SGM gives a robust and low stencil numerical scheme that allows to automatically switch from a centered scheme to an upwind one by considering advection-diffusion ratio. We present our SGM formulation and show first validation results.},
 author = {Zakari, Mustapha},
 booktitle = {2025 {RING} meeting},
 language = {en},
 pages = {186--196},
 publisher = {ASGA},
 title = {A {Scharfetter} {Gummel} based {Geothermal} model},
 year = {2025}
}