diff --git a/web/content/docs/benchmarks/thermo-hydro-mechanics/heating-homogeneous-domain/Szz.png b/web/content/docs/benchmarks/thermo-hydro-mechanics/heating-homogeneous-domain/Szz.png
new file mode 100644
index 0000000000000000000000000000000000000000..aa88967773b47f44dc6b3f7b625b7fb8e8c5d89d
Binary files /dev/null and b/web/content/docs/benchmarks/thermo-hydro-mechanics/heating-homogeneous-domain/Szz.png differ
diff --git a/web/content/docs/benchmarks/thermo-hydro-mechanics/heating-homogeneous-domain/domain.png b/web/content/docs/benchmarks/thermo-hydro-mechanics/heating-homogeneous-domain/domain.png
new file mode 100644
index 0000000000000000000000000000000000000000..6d547363790e5c47527b0b19a669c559272baeb4
Binary files /dev/null and b/web/content/docs/benchmarks/thermo-hydro-mechanics/heating-homogeneous-domain/domain.png differ
diff --git a/web/content/docs/benchmarks/thermo-hydro-mechanics/heating-homogeneous-domain/index.md b/web/content/docs/benchmarks/thermo-hydro-mechanics/heating-homogeneous-domain/index.md
new file mode 100644
index 0000000000000000000000000000000000000000..9dd64420245a5417e63d6822b0330b0becc14a4d
--- /dev/null
+++ b/web/content/docs/benchmarks/thermo-hydro-mechanics/heating-homogeneous-domain/index.md
@@ -0,0 +1,58 @@
++++
+project = ["ThermoHydroMechanics/HeatingHomogeneousDomain/hex_THM.prj"]
+author = "Wenqing Wang"
+title = "Heating of a homogeneous volume"
+date = 2023-02-10T16:12:48+01:00
+weight = 100
+image = "domain.png"
++++
+
+{{< data-link >}}
+This example is one of the mini-benchmarks of  FE-Modelling Task Force
+(by Andrés Alcolea et. al.). The example is aimed to simulate the coupled THM
+ processes in a fully water saturated unit cubic porous medium ($[0, 1]^3\,\text{m}^3$)
+ with a linear homogeneous temperature increment from
+ 20$^{\circ}$C to 30$^{\circ}$C in 100 days.
+
+The gravity is not considered in all balance equations.
+ Since the temperature is homogeneous, the specific heat capacity is set to zero
+ and thermal conductivity can be any non-zero number. The liquid density is
+ given as
+ $$\rho_L = \rho_0 \exp (\beta_L (p_L-p_0)+\alpha_L^T(T-T_{ref})) $$
+ with
+
+- $\rho_0=1002.6$ kg/m$^3$ the initial liquid density,
+- $\beta_L=4.5\cdot10^{-10}$ Pa the liquid compressibility,
+- $\alpha_L^T=-2.0\cdot 10^{-4} \text{ K}^{-1}$ the liquid thermal expansivity,
+- $T_{ref} = 273.15 \text{ K}$ the reference temperature.
+
+ While the liquid viscosity is defined as
+  $$\mu_L = \text{A}\exp(\text{B}/T)$$
+with  A=$2.1\cdot 10^{-6} \text{ Pa}\cdot\text{s}$, and B=1808.5 K.
+
+ The other material parameters are given below:
+
+| Property               | Value | Unit        |
+|------------------------|-------|-------------|
+| Young's modulus        | 1     | GPa         |
+| Poisson ratio          | 0.35  | -           |
+| Solid thermal expansion | $3 \cdot 10^{-6}$  |$\text{K}^{-1}$     |
+| Biot's coefficient   | 0.96111  | -     |
+| Porosity   | 0.1  | -     |
+| Intrinsic permeability |$3.0 \cdot 10^{-20}$ | m$^2$|
+
+Initially, the temperature is  $20 ^{\circ}$C, the pore pressure is
+ $2\cdot 10^6$ Pa, and all effective stress components are zero.
+
+At the boundary surfaces, there is no heat or flow flux, and the normal
+ displacement is fixed to zero.
+
+As a CTest, only 5 time steps with a fixed time step size $1.728\cdot 10^4$ s
+ are computed. If the example is run up to 50 time steps, corresponding to 100
+ days of simulation time, the variations of pressure and effective
+ stress can be obtained as that are shown in the two figures below:
+
+<figure>
+<img  src="p.png" alt="drawing" width="400"/>
+<img  src="Szz.png" alt="drawing" width="415"/>
+</figure>
diff --git a/web/content/docs/benchmarks/thermo-hydro-mechanics/heating-homogeneous-domain/p.png b/web/content/docs/benchmarks/thermo-hydro-mechanics/heating-homogeneous-domain/p.png
new file mode 100644
index 0000000000000000000000000000000000000000..e7ea7be191def3a0713e39f3dda2c50dbe539475
Binary files /dev/null and b/web/content/docs/benchmarks/thermo-hydro-mechanics/heating-homogeneous-domain/p.png differ