diff --git a/Tests/Data/Parabolic/T/3D_2U_BHE/3D_2U_BHE.prj b/Tests/Data/Parabolic/T/3D_2U_BHE/3D_2U_BHE.prj
index 8d4d98cff9705c2e1cb954651fcbdb963949881b..38fb238de539765531efa362fb9ac7850031bfcb 100644
--- a/Tests/Data/Parabolic/T/3D_2U_BHE/3D_2U_BHE.prj
+++ b/Tests/Data/Parabolic/T/3D_2U_BHE/3D_2U_BHE.prj
@@ -15,8 +15,8 @@
                 <borehole_heat_exchanger>
                     <type>2U</type>
                     <flow_and_temperature_control>
-                        <type>TemperatureCurveConstantFlow</type>
-                        <flow_rate>2.0e-4</flow_rate>
+                        <type>TemperatureCurveFlowCurve</type>
+                        <flow_rate_curve>flow_rate</flow_rate_curve><!-- flow rate within each U-pipe -->
                         <temperature_curve>inflow_temperature</temperature_curve>
                     </flow_and_temperature_control>
                     <borehole>
@@ -236,6 +236,11 @@
         </linear_solver>
     </linear_solvers>
     <curves>
+        <curve>
+            <name>flow_rate</name>
+            <coords>0  3600</coords>
+            <values>2.0e-4 2.0e-4</values>
+        </curve>
         <curve>
             <name>inflow_temperature</name>
             <coords>0 60 120 180 240 300 360 420 480 540 600 1200 1800 2400 3000 3600</coords>
diff --git a/web/content/docs/benchmarks/heat-transport-bhe/3D_2U_BHE.md b/web/content/docs/benchmarks/heat-transport-bhe/3D_2U_BHE.md
index 7b7a2483be0bce9779f0ba12f3c9398eeb437517..c76803e626bddd9ce1a930c9b780396971b67983 100644
--- a/web/content/docs/benchmarks/heat-transport-bhe/3D_2U_BHE.md
+++ b/web/content/docs/benchmarks/heat-transport-bhe/3D_2U_BHE.md
@@ -32,7 +32,7 @@ The numerical model was established based on dual continuum method developed by
 
 ## OGS Input Files
 
-The detailed input parameters can be seen from the 3D_2U_BHE.prj file. The inflow temperature of the BHE, which was imposed as boundary condition of the BHE is shown in Figure 1. All the initial temperatures are set as 22 $^{\circ}$C.
+The detailed input parameters can be seen from the 3D_2U_BHE.prj file. The inflow temperature of the BHE, which was imposed as boundary condition of the BHE is shown in Figure 1. All the initial temperatures are set as 22 $^{\circ}$C. The flow rate within each U-pipe is set to $2.0\times10^{-4}$ $\mathrm{m^{3} s^{-1}}$ during the whole simulation time.
 
 {{< img src="../3D_2U_BHE_figures/In_out_temperature_comparison.png" width="200">}}