# Liakopoulos TH2M/H2M benchmark: mismatch with Grunwald et al. (2022) Fig. 9 (vertical displacement profile)

**URL:** <https://discourse.opengeosys.org/t/liakopoulos-th2m-h2m-benchmark-mismatch-with-grunwald-et-al-2022-fig-9-vertical-displacement-profile/1807>\
**Category:** Help & Support\
**Tags:** ogs6\
**Created:** [February 3, 2026, 11:00am UTC](https://discourse.opengeosys.org/t/liakopoulos-th2m-h2m-benchmark-mismatch-with-grunwald-et-al-2022-fig-9-vertical-displacement-profile/1807 "2026-02-03T11:00:53Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![dbcavalcanti](https://discourse.opengeosys.org/user_avatar/discourse.opengeosys.org/dbcavalcanti/32/1449_2.png) [@dbcavalcanti](https://discourse.opengeosys.org/u/dbcavalcanti)\
**Post date:** [February 3, 2026, 11:00am UTC](https://discourse.opengeosys.org/t/liakopoulos-th2m-h2m-benchmark-mismatch-with-grunwald-et-al-2022-fig-9-vertical-displacement-profile/1807/1 "2026-02-03T11:00:53Z")

</div>

I’m comparing Fig. 9 in Grunwald et al. (2022) (Liakopoulos test) with the results shown in the OGS benchmark notebook/page for the Liakopoulos case.

In particular, the **vertical displacement profile** differs noticeably: not only the magnitude, but the **curvature appears inverted** relative to the paper figure.

Could someone clarify whether:

1. The benchmark notebook/page is intended to reproduce Fig. 9 **exactly** , or if it uses a **different configuration/variant** (e.g., different initial/boundary conditions, output times, sign convention, or plotted quantity)?
2. There were changes in the setup since the paper (material laws, parameters, output selection, etc.) that explain the discrepancy?

References:

- Grunwald et al. (2022), Fig. 9: Non-isothermal two-phase flow in deformable porous media… Geomech. Geophys. Geo-Energy Geo-Resour. 8(3), 107. DOI: 10.1007/s40948-022-00394-2
- OGS benchmark page: H2M/TH2M Liakopoulos (Jupyter): [https://www.opengeosys.org/stable/docs/benchmarks/th2m/ogs-jupyter-lab-h2m-2d-liakopoulos/](https://www.opengeosys.org/stable/docs/benchmarks/th2m/ogs-jupyter-lab-h2m-2d-liakopoulos/)

---

<div class="post-metadata">

**Author:** ![felikskiszkurno](https://discourse.opengeosys.org/user_avatar/discourse.opengeosys.org/felikskiszkurno/32/968_2.png) [@felikskiszkurno](https://discourse.opengeosys.org/u/felikskiszkurno)\
**Post date:** [February 19, 2026, 8:10am UTC](https://discourse.opengeosys.org/t/liakopoulos-th2m-h2m-benchmark-mismatch-with-grunwald-et-al-2022-fig-9-vertical-displacement-profile/1807/2 "2026-02-19T08:10:00Z")

</div>

Hi,  
thank you very much for reporting this!

While I cannot answer your questions yet, I would like to let you know, that we opened an [an issue](https://gitlab.opengeosys.org/ogs/ogs/-/issues/3601) to find the source of the mismatch and fix it.

I will post an update here when the issue is resolved.

Best,  
Feliks Kiszkurno