diff --git a/ProcessLib/ComponentTransport/Tests.cmake b/ProcessLib/ComponentTransport/Tests.cmake
index b7e977d61dc52a0d9144fa1298a370f03621ab86..91dd68099a69e2338eb2bc18e354c057ae6fe1ec 100644
--- a/ProcessLib/ComponentTransport/Tests.cmake
+++ b/ProcessLib/ComponentTransport/Tests.cmake
@@ -853,3 +853,47 @@ AddTest(
vdbc_pcs_0_ts_21990_t_210000.000000_expected.vtu vdbc_pcs_0_ts_21990_t_210000.000000.vtu concentration Si 1e-5 1e-4
vdbc_pcs_0_ts_25990_t_250000.000000_expected.vtu vdbc_pcs_0_ts_25990_t_250000.000000.vtu concentration Si 1e-5 1e-4
)
+
+AddTest(
+ NAME HC_Theis
+ PATH Parabolic/ComponentTransport/HC_Theis
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS HC_theis.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT OGS_USE_MPI
+ RUNTIME 90
+ DIFF_DATA
+ out_HC_theis_pcs_0_ts_0_t_0.000000.vtu
+ out_HC_theis_pcs_0_ts_0_t_0.000000.vtu pressure pressure 1e-3 1e-6
+ out_HC_theis_pcs_0_ts_0_t_0.000000.vtu
+ out_HC_theis_pcs_0_ts_0_t_0.000000.vtu darcy_velocity darcy_velocity 1e-10 1e-6
+ out_HC_theis_pcs_0_ts_10_t_100.000000.vtu
+ out_HC_theis_pcs_0_ts_10_t_100.000000.vtu pressure pressure 1e-3 1e-6
+ out_HC_theis_pcs_0_ts_10_t_100.000000.vtu
+ out_HC_theis_pcs_0_ts_10_t_100.000000.vtu darcy_velocity darcy_velocity 1e-10 1e-6
+ out_HC_theis_pcs_0_ts_20_t_400.000000.vtu
+ out_HC_theis_pcs_0_ts_20_t_400.000000.vtu pressure pressure 1e-3 1e-6
+ out_HC_theis_pcs_0_ts_20_t_400.000000.vtu
+ out_HC_theis_pcs_0_ts_20_t_400.000000.vtu darcy_velocity darcy_velocity 1e-10 1e-6
+ out_HC_theis_pcs_0_ts_30_t_1000.000000.vtu
+ out_HC_theis_pcs_0_ts_30_t_1000.000000.vtu pressure pressure 1e-3 1e-6
+ out_HC_theis_pcs_0_ts_30_t_1000.000000.vtu
+ out_HC_theis_pcs_0_ts_30_t_1000.000000.vtu darcy_velocity darcy_velocity 1e-10 1e-6
+ out_HC_theis_pcs_0_ts_40_t_2000.000000.vtu
+ out_HC_theis_pcs_0_ts_40_t_2000.000000.vtu pressure pressure 1e-3 1e-6
+ out_HC_theis_pcs_0_ts_40_t_2000.000000.vtu
+ out_HC_theis_pcs_0_ts_40_t_2000.000000.vtu darcy_velocity darcy_velocity 1e-10 1e-6
+ out_HC_theis_pcs_0_ts_60_t_20000.000000.vtu
+ out_HC_theis_pcs_0_ts_60_t_20000.000000.vtu pressure pressure 1e-3 1e-6
+ out_HC_theis_pcs_0_ts_60_t_20000.000000.vtu
+ out_HC_theis_pcs_0_ts_60_t_20000.000000.vtu darcy_velocity darcy_velocity 1e-10 1e-6
+ out_HC_theis_pcs_0_ts_70_t_70000.000000.vtu
+ out_HC_theis_pcs_0_ts_70_t_70000.000000.vtu pressure pressure 1e-3 1e-6
+ out_HC_theis_pcs_0_ts_70_t_70000.000000.vtu
+ out_HC_theis_pcs_0_ts_70_t_70000.000000.vtu darcy_velocity darcy_velocity 1e-10 1e-6
+ out_HC_theis_pcs_0_ts_73_t_100000.000000.vtu
+ out_HC_theis_pcs_0_ts_73_t_100000.000000.vtu pressure pressure 1e-3 1e-6
+ out_HC_theis_pcs_0_ts_73_t_100000.000000.vtu
+ out_HC_theis_pcs_0_ts_73_t_100000.000000.vtu darcy_velocity darcy_velocity 1e-10 1e-6
+)
\ No newline at end of file
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/HC_theis.prj b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/HC_theis.prj
new file mode 100644
index 0000000000000000000000000000000000000000..f5ec6c80e0c15647dcc73712f79ffe2c6095a1a8
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/HC_theis.prj
@@ -0,0 +1,302 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <meshes>
+ <mesh>slice_reorder.vtu</mesh>
+ <mesh>inflow_surf.vtu</mesh>
+ <mesh>abstraction_whole_surf.vtu</mesh>
+ </meshes>
+
+ <processes>
+ <process>
+ <name>ComponentTransport</name>
+ <type>ComponentTransport</type>
+ <integration_order>2</integration_order>
+ <process_variables>
+ <concentration>concentration</concentration>
+ <pressure>pressure</pressure>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="darcy_velocity" output_name="darcy_velocity"/>
+ </secondary_variables>
+ <!-- <calculatesurfaceflux>
+ <mesh>abstraction_whole_surf</mesh>
+ <property_name>flux</property_name>
+ <output_mesh>out_abstraction_whole_surf_flux.vtu</output_mesh>
+ </calculatesurfaceflux> -->
+ <fluid>
+ <density>
+ <type>ConcentrationAndPressureDependent</type>
+ <fluid_density_concentration_difference_ratio>0</fluid_density_concentration_difference_ratio>
+ <fluid_density_pressure_difference_ratio>5.0968E-08</fluid_density_pressure_difference_ratio>
+ <reference_concentration>0</reference_concentration>
+ <reference_pressure>0</reference_pressure>
+ <reference_density>1000</reference_density>
+ </density>
+ <viscosity>
+ <type>Constant</type>
+ <value>1.0e-3</value>
+ </viscosity>
+ </fluid>
+ <porous_medium>
+ <porous_medium id="0">
+ <permeability>
+ <type>Constant</type>
+ <permeability_tensor_entries>kappa1</permeability_tensor_entries>
+ </permeability>
+ <porosity>
+ <type>Constant</type>
+ <porosity_parameter>porosity1</porosity_parameter>
+ </porosity>
+ <!-- The implementation of the ComponentTransport process
+ is without Boussinesq approximations, thus, 'storage' is not
+ used here, but should be set through S = thickness * gravity * porosity * d\rho / dp / rho. -->
+ <storage>
+ <type>Constant</type>
+ <value>0</value>
+ </storage>
+ </porous_medium>
+ </porous_medium>
+ <fluid_reference_density>rho_fluid</fluid_reference_density>
+ <retardation_factor>retardation</retardation_factor>
+ <decay_rate>decay</decay_rate>
+ <specific_body_force>0 0 -9.81</specific_body_force>
+ </process>
+ </processes>
+ <media>
+ <medium id="0">
+ <phases>
+ <phase>
+ <type>AqueousLiquid</type>
+ <components>
+ <component>
+ <name>concentration</name>
+ <properties>
+ <property>
+ <name>molecular_diffusion</name>
+ <type>Constant</type>
+ <value>2e-9</value>
+ </property>
+ </properties>
+ </component>
+ </components>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>longitudinal_dispersivity</name>
+ <type>Constant</type>
+ <value>0</value>
+ </property>
+ <property>
+ <name>transversal_dispersivity</name>
+ <type>Constant</type>
+ <value>0</value>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <time_loop>
+ <processes>
+ <process ref="ComponentTransport">
+ <nonlinear_solver>basic_picard</nonlinear_solver>
+ <convergence_criterion>
+ <type>PerComponentDeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltols>1e-5 5e-6</reltols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial> 0.0 </t_initial>
+ <t_end> 1e5 </t_end>
+ <timesteps>
+ <pair>
+ <repeat>10</repeat>
+ <delta_t>10</delta_t>
+ </pair>
+ <pair>
+ <repeat>5</repeat>
+ <delta_t>20</delta_t>
+ </pair>
+ <pair>
+ <repeat>10</repeat>
+ <delta_t>40</delta_t>
+ </pair>
+ <pair>
+ <repeat>5</repeat>
+ <delta_t>80</delta_t>
+ </pair>
+ <pair>
+ <repeat>10</repeat>
+ <delta_t>100</delta_t>
+ </pair>
+ <pair>
+ <repeat>5</repeat>
+ <delta_t>200</delta_t>
+ </pair>
+ <pair>
+ <repeat>5</repeat>
+ <delta_t>400</delta_t>
+ </pair>
+ <pair>
+ <repeat>5</repeat>
+ <delta_t>1000</delta_t>
+ </pair>
+ <pair>
+ <repeat>5</repeat>
+ <delta_t>2000</delta_t>
+ </pair>
+ <pair>
+ <repeat>10</repeat>
+ <delta_t>5000</delta_t>
+ </pair>
+ <pair>
+ <repeat>5</repeat>
+ <delta_t>10000</delta_t>
+ </pair>
+ <pair>
+ <repeat>5</repeat>
+ <delta_t>20000</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>out_HC_theis</prefix>
+ <timesteps>
+ <pair>
+ <repeat> 1 </repeat>
+ <each_steps> 10 </each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable> pressure </variable>
+ <variable> concentration </variable>
+ <variable> darcy_velocity </variable>
+ </variables>
+ </output>
+ </time_loop>
+
+ <parameters>
+ <parameter>
+ <name>rho_fluid</name>
+ <type>Constant</type>
+ <value>1000</value>
+ </parameter>
+ <parameter>
+ <name>retardation</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>decay</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+
+ <parameter>
+ <name>porosity1</name>
+ <type>Constant</type>
+ <value>0.2</value>
+ </parameter>
+ <parameter>
+ <name>kappa1</name>
+ <type>Constant</type>
+ <values>7.6453E-13 0 0 0 7.6453E-13 0 0 0 7.6453E-13</values>
+ </parameter>
+
+ <parameter>
+ <name>c_ini</name>
+ <type>MeshNode</type>
+ <mesh>slice_reorder</mesh>
+ <field_name>c_ini</field_name>
+ </parameter>
+ <parameter>
+ <name>p_ini</name>
+ <type>MeshNode</type>
+ <mesh>slice_reorder</mesh>
+ <field_name>p_ini</field_name>
+ </parameter>
+ <parameter>
+ <name>c_0</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>p_BC_Dirichlet</name>
+ <type>MeshNode</type>
+ <mesh>inflow_surf</mesh>
+ <field_name>p_BC_Dirichlet</field_name>
+ </parameter>
+ <parameter>
+ <name>p_BC_Neumann_c</name>
+ <type>Constant</type>
+ <value>-2.83542E-03</value>
+ </parameter>
+
+
+ </parameters>
+
+ <process_variables>
+ <process_variable>
+ <name>concentration</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>c_ini</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <type>Dirichlet</type>
+ <mesh>inflow_surf</mesh>
+ <parameter>c_0</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>p_ini</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <type>Dirichlet</type>
+ <mesh>inflow_surf</mesh>
+ <parameter>p_BC_Dirichlet</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <type>Neumann</type>
+ <mesh>abstraction_whole_surf</mesh>
+ <parameter>p_BC_Neumann_c</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_picard</name>
+ <type>Picard</type>
+ <max_iter>30</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <lis>-i cg -p jacobi -tol 1e-10 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>CG</solver_type>
+ <precon_type>DIAGONAL</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-10</error_tolerance>
+ </eigen>
+ <petsc>
+ <prefix>lf</prefix>
+ <parameters>-lf_ksp_type cg -lf_pc_type bjacobi -lf_ksp_rtol 1e-16 -lf_ksp_max_it 10000</parameters>
+ </petsc>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/abstraction_whole_surf.vtu b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/abstraction_whole_surf.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..79c0ff13078fbadc6b0b6f73aba343dc86f8147b
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/abstraction_whole_surf.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:4c73a8b4660971a61bd186fb837e306e6590d635a7f1054cf086a22098c8b70b
+size 9506
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/analytical_theis.xlsx b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/analytical_theis.xlsx
new file mode 100644
index 0000000000000000000000000000000000000000..d8b84cbf0f4b0d6901d39a6eff891ca8457ceab2
Binary files /dev/null and b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/analytical_theis.xlsx differ
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/drawdown_s-30.txt b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/drawdown_s-30.txt
new file mode 100644
index 0000000000000000000000000000000000000000..82a5fc7ab5923fa0382575df63df4ad30f4c77ba
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/drawdown_s-30.txt
@@ -0,0 +1,25 @@
+t,s
+500,7.7188116E-05
+850,1.2380776E-03
+1445,8.1109036E-03
+2456.5,2.8216319E-02
+4176.05,6.6501297E-02
+7099.285,1.2267865E-01
+12068.7845,1.9318485E-01
+20516.93365,2.7381783E-01
+34878.78721,3.6109168E-01
+59293.93825,4.5253086E-01
+100799.695,5.4651475E-01
+171359.4815,6.4202926E-01
+291311.1186,7.3845604E-01
+495228.9016,8.3542360E-01
+841889.1328,9.3271071E-01
+1431211.526,1.0301863E+00
+2433059.594,1.1277729E+00
+4136201.309,1.2254250E+00
+7031542.226,1.3231155E+00
+11953621.78,1.4208286E+00
+20321157.03,1.5185551E+00
+34545966.96,1.6162894E+00
+58728143.83,1.7140283E+00
+99837844.5,1.8117699E+00
\ No newline at end of file
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/drawdown_t-5e3.txt b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/drawdown_t-5e3.txt
new file mode 100644
index 0000000000000000000000000000000000000000..76739684204319cb059205a7237ceacab24c54b7
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/drawdown_t-5e3.txt
@@ -0,0 +1,202 @@
+r,s
+0.5,1.49622E+00
+1,1.24094E+00
+1.5,1.09172E+00
+2,9.85946E-01
+2.5,9.04012E-01
+3,8.37179E-01
+3.5,7.80786E-01
+4,7.32049E-01
+4.5,6.89175E-01
+5,6.50938E-01
+5.5,6.16463E-01
+6,5.85105E-01
+6.5,5.56374E-01
+7,5.29888E-01
+7.5,5.05345E-01
+8,4.82501E-01
+8.5,4.61156E-01
+9,4.41146E-01
+9.5,4.22331E-01
+10,4.04594E-01
+10.5,3.87834E-01
+11,3.71966E-01
+11.5,3.56914E-01
+12,3.42613E-01
+12.5,3.29006E-01
+13,3.16040E-01
+13.5,3.03671E-01
+14,2.91859E-01
+14.5,2.80568E-01
+15,2.69764E-01
+15.5,2.59419E-01
+16,2.49505E-01
+16.5,2.39998E-01
+17,2.30876E-01
+17.5,2.22118E-01
+18,2.13705E-01
+18.5,2.05621E-01
+19,1.97848E-01
+19.5,1.90373E-01
+20,1.83181E-01
+20.5,1.76260E-01
+21,1.69597E-01
+21.5,1.63182E-01
+22,1.57003E-01
+22.5,1.51051E-01
+23,1.45316E-01
+23.5,1.39790E-01
+24,1.34465E-01
+24.5,1.29332E-01
+25,1.24384E-01
+25.5,1.19615E-01
+26,1.15017E-01
+26.5,1.10584E-01
+27,1.06310E-01
+27.5,1.02189E-01
+28,9.82161E-02
+28.5,9.43856E-02
+29,9.06926E-02
+29.5,8.71323E-02
+30,8.36999E-02
+30.5,8.03912E-02
+31,7.72017E-02
+31.5,7.41275E-02
+32,7.11646E-02
+32.5,6.83091E-02
+33,6.55576E-02
+33.5,6.29064E-02
+34,6.03521E-02
+34.5,5.78916E-02
+35,5.55215E-02
+35.5,5.32390E-02
+36,5.10411E-02
+36.5,4.89249E-02
+37,4.68876E-02
+37.5,4.49267E-02
+38,4.30395E-02
+38.5,4.12236E-02
+39,3.94765E-02
+39.5,3.77959E-02
+40,3.61796E-02
+40.5,3.46254E-02
+41,3.31312E-02
+41.5,3.16949E-02
+42,3.03145E-02
+42.5,2.89880E-02
+43,2.77138E-02
+43.5,2.64898E-02
+44,2.53145E-02
+44.5,2.41859E-02
+45,2.31027E-02
+45.5,2.20630E-02
+46,2.10655E-02
+46.5,2.01085E-02
+47,1.91907E-02
+47.5,1.83105E-02
+48,1.74668E-02
+48.5,1.66581E-02
+49,1.58831E-02
+49.5,1.51407E-02
+50,1.44296E-02
+50.5,1.37486E-02
+51,1.30967E-02
+51.5,1.24727E-02
+52,1.18756E-02
+52.5,1.13044E-02
+53,1.07580E-02
+53.5,1.02356E-02
+54,9.73611E-03
+54.5,9.25876E-03
+55,8.80265E-03
+55.5,8.36693E-03
+56,7.95081E-03
+56.5,7.55350E-03
+57,7.17425E-03
+57.5,6.81233E-03
+58,6.46704E-03
+58.5,6.13769E-03
+59,5.82363E-03
+59.5,5.52424E-03
+60,5.23890E-03
+60.5,4.96702E-03
+61,4.70804E-03
+61.5,4.46140E-03
+62,4.22660E-03
+62.5,4.00310E-03
+63,3.79044E-03
+63.5,3.58814E-03
+64,3.39574E-03
+64.5,3.21281E-03
+65,3.03894E-03
+65.5,2.87371E-03
+66,2.71675E-03
+66.5,2.56767E-03
+67,2.42613E-03
+67.5,2.29178E-03
+68,2.16429E-03
+68.5,2.04334E-03
+69,1.92863E-03
+69.5,1.81987E-03
+70,1.71678E-03
+70.5,1.61908E-03
+71,1.52654E-03
+71.5,1.43889E-03
+72,1.35590E-03
+72.5,1.27735E-03
+73,1.20301E-03
+73.5,1.13270E-03
+74,1.06620E-03
+74.5,1.00332E-03
+75,9.43897E-04
+75.5,8.87744E-04
+76,8.34700E-04
+76.5,7.84608E-04
+77,7.37316E-04
+77.5,6.92682E-04
+78,6.50568E-04
+78.5,6.10844E-04
+79,5.73385E-04
+79.5,5.38072E-04
+80,5.04793E-04
+80.5,4.73440E-04
+81,4.43910E-04
+81.5,4.16106E-04
+82,3.89935E-04
+82.5,3.65310E-04
+83,3.42146E-04
+83.5,3.20364E-04
+84,2.99890E-04
+84.5,2.80651E-04
+85,2.62581E-04
+85.5,2.45615E-04
+86,2.29695E-04
+86.5,2.14763E-04
+87,2.00767E-04
+87.5,1.87658E-04
+88,1.75389E-04
+88.5,1.63920E-04
+89,1.53211E-04
+89.5,1.43229E-04
+90,1.33943E-04
+90.5,1.25326E-04
+91,1.17359E-04
+91.5,1.10024E-04
+92,1.03313E-04
+92.5,9.72210E-05
+93,9.17531E-05
+93.5,8.69220E-05
+94,8.27511E-05
+94.5,7.92755E-05
+95,7.65445E-05
+95.5,7.46240E-05
+96,7.35996E-05
+96.5,7.35805E-05
+97,7.47035E-05
+97.5,7.71390E-05
+98,8.10971E-05
+98.5,8.68351E-05
+99,9.46670E-05
+99.5,1.04974E-04
+100,1.18220E-04
+100.5,1.34961E-04
\ No newline at end of file
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/inflow_surf.vtu b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/inflow_surf.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..b2fffd1b6708a755ec981e3283998f94ebbcfd09
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/inflow_surf.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:f2b29416fdf4271d8e9908cba7cb8b443c66d39694702496fbb70d5e441b123e
+size 37891
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_0_t_0.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_0_t_0.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..15662e9cf0b5fa86f86a85d2e1ecefa60aaf54c9
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_0_t_0.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:3fd64ba69a006f30ef0decff82c83a3a194101a94607e082e66aad2dcfa5a4bb
+size 1253719
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_10_t_100.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_10_t_100.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..9916035db0a2093868f235026f840b801b1a1961
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_10_t_100.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:f33f2bb4a3ad67b9ed9a27e9674e96200166df0834c2982042f865edeeed3794
+size 1498401
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_20_t_400.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_20_t_400.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..9ef8a0e4584ce14630cc8c31b63979acedfe9b92
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_20_t_400.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:d0f0a04f7e5c6e0c9c8cf805b2f41890c5e3f4f369a00798ba4a549b26e219c2
+size 1522353
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_30_t_1000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_30_t_1000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..4813d79136a68bc37f721c33995fdbfb62de7f88
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_30_t_1000.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:e7f8d77218cb0429f591bbbd28f678efe1f024c25b6e554cf7bc7f7a40c73365
+size 1539968
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_40_t_2000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_40_t_2000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..a858880d717c58a8cccd231b0aba8ee956e8883c
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_40_t_2000.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:116b5083806edf70fc46efb05ab4a672cf805a29910094e16811195573eabbb4
+size 1544993
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_50_t_5000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_50_t_5000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..ced504d208ed9d9e7bddebe21173ab41ad6df3de
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_50_t_5000.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:b2d8b10700e27e29c1a57f49123e23455b9d23d708e373655e77dee7e67c6f74
+size 1547141
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_60_t_20000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_60_t_20000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..8098dffe75e0641e268ac84504acc389e3f08ca2
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_60_t_20000.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:64a68e74adf791a2243d724331ab930b2ee3fedae3e1b78e109e7f8a3c10d59c
+size 1549264
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_70_t_70000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_70_t_70000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..71b3e8f76c6a932f42af12557bde2d030ba9f0b1
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_70_t_70000.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:83614603d6e995408a4f9a497c71556465e5ece0ed158b90b0852a35e8c6dd8c
+size 1550301
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_73_t_100000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_73_t_100000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..babf54e98db8e3c92382903b225373e8231eebad
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/out_HC_theis_pcs_0_ts_73_t_100000.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:66e9ff221f1fa1913a1770ff10d05c1f4f82f4084ee06d75e05c36ffc9b7a834
+size 1549753
diff --git a/Tests/Data/Parabolic/ComponentTransport/HC_Theis/slice_reorder.vtu b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/slice_reorder.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..3184d9a89caf5b026bafaee280991caf84b75808
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/HC_Theis/slice_reorder.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:aae47686ae64a6d738de18b4d10f6af672f042974026f5d4c55ea4e89ac7adfb
+size 3568220
diff --git a/web/content/docs/benchmarks/hydro-component/theis/BCs.png b/web/content/docs/benchmarks/hydro-component/theis/BCs.png
new file mode 100644
index 0000000000000000000000000000000000000000..cb66321ab7c4ee5a4f7c87b8f99d7cc2e15a2b78
--- /dev/null
+++ b/web/content/docs/benchmarks/hydro-component/theis/BCs.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:cb910364c3ab0879c2daf2d30beffbb5de9ab87138bea6c337e64872a696704f
+size 416847
diff --git a/web/content/docs/benchmarks/hydro-component/theis/HC_Theis.pandoc b/web/content/docs/benchmarks/hydro-component/theis/HC_Theis.pandoc
new file mode 100644
index 0000000000000000000000000000000000000000..b45b6f67e2593ab19f49146ceb6fc9d4c0be7cd1
--- /dev/null
+++ b/web/content/docs/benchmarks/hydro-component/theis/HC_Theis.pandoc
@@ -0,0 +1,65 @@
++++
+author = "Marc Walther"
+weight = 142
+project = "Parabolic/ComponentTransport/Theis/HC_theis.prj"
+date = "2019-05-02T14:41:09+01:00"
+title = "Theis solution for well pumping"
+
+[menu]
+ [menu.benchmarks]
+ parent = "hydro-component"
+
++++
+
+{{< data-link >}}
+
+
+## Overview
+
+Abstraction of water from an aquifer is a common procedure, as groundwater is usually better protected from contamination and steadily available even in dry seasons. For a sustainable water usage, the evaluation of the spatio-temporal behaviour of the water level for a specified abstraction rate is of high importance for water management.
+
+This benchmark verifies the groundwater level drawdown for an aquifer that is subjected to pumping.
+
+
+## Problem description
+
+For pumping from a well, Theis (1935) proposed an analytical solution to calculate the water level drawdown over time and distance from the well. A short summary of the solution can be seen in
+web/content/docs/benchmarks/liquid-flow/liquid-flow-theis-problem.pandoc
+
+Here, we verify pumping abstraction with Theis for the `ComponentTransport` process using a 3D setup. The setup is adapted from Walther (2014), section 3.2.
+
+
+### Model setup
+
+The setup comprises a 1/8th slice of a full circle (see figure 1).
+
+{{< img src="../BCs.png" title="Mesh and boundary conditions (BC); blue = outer Dirichlet pressure and concentration BC, red = inner Neumann abstraction BC.">}}
+
+The outer boundary condition is set as Dirichlet with a hydrostatic pressure along the shell surface of the slice equivalent to a head of $h = 0 m$ (i.e. water level equals top of domain). For mass transport, a Dirichlet boundary conditions with concentration $c = 0$ is set at the outer shell. The inner boundary condition is equivalent to the eighth of a total abstraction rate of $Q_t = 15 m^3/d$ for a full cylinder. *NB: In the `ComponentTransport` process, the Neumann BC is given as mass flux and has to be calculated per area, such that the value for the project file is $Q = Q_t / 8 / A \cdot \rho = 2.83542E-03 m^3/s/m^2 \cdot kg/m^3$ with fluid density $\rho = 1000 kg/m^3$ and abstraction area $A = 7.65 m^2$.*
+
+
+The homogeneous, isotropic domain is defined for the radius $1 < r < 100 m$ and a thickness $b = 10 m$. Saturated intrinsic permeability is $\kappa = 7.6453E-13 m^2$ yielding a transmissivity of $T = 7.5E-05 m^2/s$; porosity is $\phi = 0.2$; specific storage is $S_s = 1.0E-03$ and defined through compressibility $\gamma = 5.0968E-08 s^2/m/kg$ (input tag fluid_density_pressure_difference_ratio is $\frac{\partial \rho}{\partial p} = \frac{S_s}{b \phi g \rho}$ with gravitational acceleration $g = 9.81 m^2/s$).
+
+Mass transport properties are irrelevant as no transport processes are calculated.
+
+Initial conditions are $c = 0$ and hydrostatic pressure conditions.
+
+{{< data-link "Link to the HC_theis project file" "Parabolic/ComponentTransport/HC_Theis/HC_theis.prj" >}}
+
+
+### Results
+
+The figure below compares the analytical Theis solution vs. the simulated values from OGS6.
+
+{{< img src="../comparison.png" title="Comparison between numerical (crosses) and analytical (lines) values.">}}
+
+The top figure shows drawdown (i.e. the difference in water level compared to an initial state) over time at distance $r = 30 m$: for a simulation time $t < 40000 s$, the differences between analytical and numerical solutions are marginal; at later simulation times, the drawdown shows lower values than predicted from the analytical solution, as it is influenced by the outer Dirichlet pressure boundary condition.
+
+The lower figure shows drawdown over radius at time $t = 5000 s$; the top figure also signifies this with a green vertical: over the whole domain, the analytical and numerical solutions are coinciding very well. *NB: With higher values of $S_s$, one has to take care to convert the simulated water level to equivalent freshwater head; in this case, fluid density variation is very low such that a conversion is not required.*
+
+
+## References
+
+Theis, C.V., 1935. The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using groundwater storage. Trans. Amer. Geophys. Union 16, 513–514.
+
+Walther, M., 2014. Variable-Density Flow Processes in Porous Media On Small, Medium and Regional Scales (Dissertation). Technische Universität Dresden.
diff --git a/web/content/docs/benchmarks/hydro-component/theis/comparison.png b/web/content/docs/benchmarks/hydro-component/theis/comparison.png
new file mode 100644
index 0000000000000000000000000000000000000000..19a127723600592bd54eae1a677f950cf8c17138
--- /dev/null
+++ b/web/content/docs/benchmarks/hydro-component/theis/comparison.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:f846c43166cbf6379948690261b67f48d6fc44d8082181b2e08f3ff9a3dbeba9
+size 38126