diff --git a/test/porousmediumflow/1p/implicit/1ptestproblem.hh b/test/porousmediumflow/1p/implicit/1ptestproblem.hh
index fdad5160ec2cb9894b09bd0f682e3db11a791c13..dc893c2572a3be07652fc97701f59deb9ee67e85 100644
--- a/test/porousmediumflow/1p/implicit/1ptestproblem.hh
+++ b/test/porousmediumflow/1p/implicit/1ptestproblem.hh
@@ -37,6 +37,10 @@
#include "1ptestspatialparams.hh"
+#if FORCHHEIMER
+#include <dumux/discretization/forchheimerslaw.hh>
+#endif
+
namespace Dumux
{
template <class TypeTag>
@@ -71,6 +75,10 @@ SET_PROP(OnePTestTypeTag, SpatialParams)
using Scalar = typename GET_PROP_TYPE(TypeTag, Scalar);
using type = OnePTestSpatialParams<FVGridGeometry, Scalar>;
};
+
+#ifdef FORCHHEIMER
+SET_TYPE_PROP(OnePTestTypeTag, AdvectionType, ForchheimersLaw<TypeTag>);
+#endif
}
/*!
diff --git a/test/porousmediumflow/1p/implicit/CMakeLists.txt b/test/porousmediumflow/1p/implicit/CMakeLists.txt
index 268cfce8876ce6ec5725970d02e2cb398d5b67a4..fcf04f8f5583f31905c9a1de7a5df05b754cc966 100644
--- a/test/porousmediumflow/1p/implicit/CMakeLists.txt
+++ b/test/porousmediumflow/1p/implicit/CMakeLists.txt
@@ -37,6 +37,16 @@ dune_add_test(NAME test_1pbox
${CMAKE_CURRENT_BINARY_DIR}/1ptestbox-00001.vtu
--command "${CMAKE_CURRENT_BINARY_DIR}/test_1pbox test_1pfv.input -Problem.Name 1ptestbox")
+dune_add_test(NAME test_1pforchheimercctpfa
+ SOURCES test_1pfv.cc
+ COMPILE_DEFINITIONS TYPETAG=OnePTestCCTpfaTypeTag FORCHHEIMER=1
+ COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+ CMD_ARGS --script fuzzy
+ --files ${CMAKE_SOURCE_DIR}/test/references/1pforchheimercctpfa.vtu
+ ${CMAKE_CURRENT_BINARY_DIR}/1ptestforchheimertpfa-00001.vtu
+ --command "${CMAKE_CURRENT_BINARY_DIR}/test_1pforchheimercctpfa"
+ --zeroThreshold {"velocity_liq \(m/s\)":1e-12})
+
# a gstat test (becaue it's a random permeability field we can't test against a reference solution)
dune_add_test(NAME test_1pwithgstat
SOURCES test_1pfv.cc
diff --git a/test/porousmediumflow/1p/implicit/test_1pforchheimercctpfa.input b/test/porousmediumflow/1p/implicit/test_1pforchheimercctpfa.input
new file mode 100644
index 0000000000000000000000000000000000000000..4a13c84901c1e1ce6d3deae077c80be973d598a8
--- /dev/null
+++ b/test/porousmediumflow/1p/implicit/test_1pforchheimercctpfa.input
@@ -0,0 +1,23 @@
+[TimeLoop]
+DtInitial = 1 # [s]
+TEnd = 1 # [s]
+
+[Grid]
+LowerLeft = 0 0
+UpperRight = 1 1
+Cells = 2 10
+
+[Problem]
+Name = 1ptestforchheimertpfa # name passed to the output routines
+EnableGravity = false
+
+[SpatialParams]
+LensLowerLeft = 0.2 0.2
+LensUpperRight = 0.8 0.8
+
+Permeability = 1e-7 # [m^2]
+PermeabilityLens = 1e-7 # [m^2]
+ForchCoeff = 0.55
+
+[Vtk]
+AddVelocity = true
diff --git a/test/references/1pforchheimercctpfa.vtu b/test/references/1pforchheimercctpfa.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..2137998a675f8842813bc64a88a1983ef05742ca
--- /dev/null
+++ b/test/references/1pforchheimercctpfa.vtu
@@ -0,0 +1,56 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+ <UnstructuredGrid>
+ <Piece NumberOfCells="20" NumberOfPoints="33">
+ <CellData Scalars="p" Vectors="velocity_liq (m/s)">
+ <DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii">
+ 195000 195000 185000 185000 175000 175000 165000 165000 155000 155000 145000 145000
+ 135000 135000 125000 125000 115000 115000 105000 105000
+ </DataArray>
+ <DataArray type="Float32" Name="velocity_liq (m/s)" NumberOfComponents="3" format="ascii">
+ 0 0.236926 0 -0 0.236926 0 0 0.236926 0 -0 0.236926 0
+ 0 0.236926 0 -0 0.236926 0 0 0.236926 0 -0 0.236926 0
+ 5.82077e-15 0.236926 0 5.82077e-15 0.236926 0 0 0.236926 0 -0 0.236926 0
+ 0 0.236926 0 -0 0.236926 0 0 0.236926 0 -0 0.236926 0
+ 0 0.236926 0 -0 0.236926 0 0 0.236926 0 -0 0.236926 0
+ </DataArray>
+ <DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0
+ </DataArray>
+ </CellData>
+ <Points>
+ <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
+ 0 0 0 0.5 0 0 0 0.1 0 0.5 0.1 0
+ 1 0 0 1 0.1 0 0 0.2 0 0.5 0.2 0
+ 1 0.2 0 0 0.3 0 0.5 0.3 0 1 0.3 0
+ 0 0.4 0 0.5 0.4 0 1 0.4 0 0 0.5 0
+ 0.5 0.5 0 1 0.5 0 0 0.6 0 0.5 0.6 0
+ 1 0.6 0 0 0.7 0 0.5 0.7 0 1 0.7 0
+ 0 0.8 0 0.5 0.8 0 1 0.8 0 0 0.9 0
+ 0.5 0.9 0 1 0.9 0 0 1 0 0.5 1 0
+ 1 1 0
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+ 0 1 3 2 1 4 5 3 2 3 7 6
+ 3 5 8 7 6 7 10 9 7 8 11 10
+ 9 10 13 12 10 11 14 13 12 13 16 15
+ 13 14 17 16 15 16 19 18 16 17 20 19
+ 18 19 22 21 19 20 23 22 21 22 25 24
+ 22 23 26 25 24 25 28 27 25 26 29 28
+ 27 28 31 30 28 29 32 31
+ </DataArray>
+ <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii">
+ 4 8 12 16 20 24 28 32 36 40 44 48
+ 52 56 60 64 68 72 76 80
+ </DataArray>
+ <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii">
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>