diff --git a/ProcessLib/ComponentTransport/Tests.cmake b/ProcessLib/ComponentTransport/Tests.cmake
index c0075687e6764c9ae877d6894649ad245d4c6f4e..564882f4a5231767a11c56982ceaf37a105dee8e 100644
--- a/ProcessLib/ComponentTransport/Tests.cmake
+++ b/ProcessLib/ComponentTransport/Tests.cmake
@@ -193,6 +193,18 @@ AddTest(
VIS DiffusionAndStorageAndAdvectionAndDispersionHalf_pcs_0_ts_672_t_900.000000.vtu
)
+AddTest(
+ NAME 3D_ComponentTransport_surfaceflux
+ PATH Parabolic/ComponentTransport/SimpleSynthetics
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS surfaceflux_component-transport_cube_1e3.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT OGS_USE_MPI
+ DIFF_DATA
+ cube_1x1x1_hex_1e3_complete_surface.vtu flux_1e3_t_1.000000.vtu flux flux 1e-10 1e-16
+)
+
AddTest(
NAME LARGE_2D_ComponentTransport_Goswami
PATH Parabolic/ComponentTransport/goswami
diff --git a/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/cube_1x1x1_hex_1e3.vtu b/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/cube_1x1x1_hex_1e3.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..5f491d77c01fe3246feba3443d5548791e934e82
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/cube_1x1x1_hex_1e3.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:61a70995845833f16ba3b2edf0312eca29868be20f2cc7015df06d1f35eed669
+size 22546
diff --git a/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/cube_1x1x1_hex_1e3_complete_surface.vtu b/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/cube_1x1x1_hex_1e3_complete_surface.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..8a2407e5f3183a2091ee355e342e9c8ab505e457
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/cube_1x1x1_hex_1e3_complete_surface.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:bcb45ed892b69e6ddf23038d5d6dbe292a2f95448b97e6401336bbdb4f8982e5
+size 79181
diff --git a/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/cube_1x1x1_hex_1e3_left.vtu b/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/cube_1x1x1_hex_1e3_left.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..b9be1a60800686b8a6c073613bd5169f6f69ccad
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/cube_1x1x1_hex_1e3_left.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:074db4f230b18b0fc6300e7b49669b55cd2e06b948b1344302305a59261de79d
+size 12897
diff --git a/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/cube_1x1x1_hex_1e3_right.vtu b/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/cube_1x1x1_hex_1e3_right.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..8c9a97647805e21f6efe641cec909daae8b5f394
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/cube_1x1x1_hex_1e3_right.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:8427200cf2b82af3762994c1e259839aa5b502a8d5b8a7dce60919f70fa04e41
+size 12898
diff --git a/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/surfaceflux_component-transport_cube_1e3.prj b/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/surfaceflux_component-transport_cube_1e3.prj
new file mode 100644
index 0000000000000000000000000000000000000000..e0a53624697ec24f631eacd9d1d32af6b456ab7d
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/surfaceflux_component-transport_cube_1e3.prj
@@ -0,0 +1,240 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <meshes>
+ <mesh>cube_1x1x1_hex_1e3.vtu</mesh>
+ <mesh>cube_1x1x1_hex_1e3_left.vtu</mesh>
+ <mesh>cube_1x1x1_hex_1e3_right.vtu</mesh>
+ <mesh>cube_1x1x1_hex_1e3_complete_surface.vtu</mesh>
+ </meshes>
+ <processes>
+ <process>
+ <name>hc</name>
+ <type>ComponentTransport</type>
+ <integration_order>2</integration_order>
+ <process_variables>
+ <concentration>concentration</concentration>
+ <pressure>pressure</pressure>
+ </process_variables>
+ <fluid>
+ <density>
+ <type>ConcentrationDependent</type>
+ <reference_density>1</reference_density>
+ <reference_concentration>0</reference_concentration>
+ <fluid_density_difference_ratio>0.0</fluid_density_difference_ratio>
+ </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>constant_porosity_parameter</porosity_parameter>
+ </porosity>
+ <storage>
+ <type>Constant</type>
+ <value>0</value>
+ </storage>
+ </porous_medium>
+ </porous_medium>
+ <fluid_reference_density>rho_fluid</fluid_reference_density>
+ <solute_dispersivity_longitudinal>beta_l</solute_dispersivity_longitudinal>
+ <solute_dispersivity_transverse>beta_t</solute_dispersivity_transverse>
+ <molecular_diffusion_coefficient>Dm</molecular_diffusion_coefficient>
+ <retardation_factor>retardation</retardation_factor>
+ <decay_rate>decay</decay_rate>
+ <specific_body_force>0 0 0</specific_body_force>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="darcy_velocity" output_name="darcy_velocity"/>
+ </secondary_variables>
+ <calculatesurfaceflux>
+ <mesh>cube_1x1x1_hex_1e3_complete_surface</mesh>
+ <property_name>flux</property_name>
+ <output_mesh>flux_1e3.vtu</output_mesh>
+ </calculatesurfaceflux>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="hc">
+ <nonlinear_solver>basic_picard</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1e-1</reltol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <output>
+ <variables>
+ <variable>concentration</variable>
+ <variable>pressure</variable>
+ <variable>darcy_velocity</variable>
+ </variables>
+ </output>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0.0</t_initial>
+ <t_end>1</t_end>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>DiffusionOnly_surfaceflux</prefix>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <each_steps>1</each_steps>
+ </pair>
+ </timesteps>
+ </output>
+ </time_loop>
+ <parameters>
+ <parameter>
+ <name>kappa1</name>
+ <type>Constant</type>
+ <values>1.239e-7 0 0 0 1.239e-7 0 0 0 1.239e-7</values>
+ </parameter>
+ <parameter>
+ <name>constant_porosity_parameter</name>
+ <type>Constant</type>
+ <value>0.2</value>
+ </parameter>
+ <parameter>
+ <name>rho_fluid</name>
+ <type>Constant</type>
+ <value>1000</value>
+ </parameter>
+ <parameter>
+ <name>Dm</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>retardation</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>decay</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>beta_l</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>beta_t</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>c0</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>p0</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>p_Dirichlet_left</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>p_Dirichlet_right</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>c_Dirichlet_left</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>c_Dirichlet_right</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>concentration</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>c0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>cube_1x1x1_hex_1e3_left</mesh>
+ <type>Dirichlet</type>
+ <parameter>c_Dirichlet_left</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>cube_1x1x1_hex_1e3_right</mesh>
+ <type>Dirichlet</type>
+ <parameter>c_Dirichlet_right</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>p0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>cube_1x1x1_hex_1e3_left</mesh>
+ <type>Dirichlet</type>
+ <parameter>p_Dirichlet_left</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>cube_1x1x1_hex_1e3_right</mesh>
+ <type>Dirichlet</type>
+ <parameter>p_Dirichlet_right</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_picard</name>
+ <type>Picard</type>
+ <max_iter>2</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-16 -maxiter 20000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>ILUT</precon_type>
+ <max_iteration_step>20000</max_iteration_step>
+ <error_tolerance>1e-8</error_tolerance>
+ </eigen>
+ <petsc>
+ <prefix>hc</prefix>
+ <parameters>-hc_ksp_type bcgs -hc_pc_type bjacobi -hc_ksp_rtol 1e-8 -hc_ksp_max_it 20000</parameters>
+ </petsc>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>