diff --git a/test/freeflow/navierstokes/channel/2d/problem.hh b/test/freeflow/navierstokes/channel/2d/problem.hh
index 85df157ed0e6974a45afffc80ae9b232c751df37..344a0526a5bda30a932869be4d7fc9844919db77 100644
--- a/test/freeflow/navierstokes/channel/2d/problem.hh
+++ b/test/freeflow/navierstokes/channel/2d/problem.hh
@@ -148,15 +148,6 @@ public:
     Scalar temperature() const
     { return 273.15 + 10; } // 10C
 
-   /*!
-     * \brief Returns the sources within the domain.
-     *
-     * \param globalPos The global position
-     */
-    NumEqVector sourceAtPos(const GlobalPosition &globalPos) const
-    {
-        return NumEqVector(0.0);
-    }
     // \}
    /*!
      * \name Boundary conditions
diff --git a/test/freeflow/navierstokes/closedsystem/problem.hh b/test/freeflow/navierstokes/closedsystem/problem.hh
index 49d1452457359528e998d01afdf3037a2912410b..1d8af19ef7c544485e139582dd1ad30e1db12646 100644
--- a/test/freeflow/navierstokes/closedsystem/problem.hh
+++ b/test/freeflow/navierstokes/closedsystem/problem.hh
@@ -122,15 +122,6 @@ public:
     Scalar temperature() const
     { return 273.15 + 10; } // 10C
 
-   /*!
-     * \brief Returns the sources within the domain.
-     *
-     * \param globalPos The global position
-     */
-    NumEqVector sourceAtPos(const GlobalPosition &globalPos) const
-    {
-        return NumEqVector(0.0);
-    }
     // \}
    /*!
      * \name Boundary conditions
diff --git a/test/freeflow/navierstokes/kovasznay/problem.hh b/test/freeflow/navierstokes/kovasznay/problem.hh
index 077caa6e5ba5f711cbecfc80e1ce1efa68169322..a24810faa77a496f8389e4aa6ea69c4faef1ddbb 100644
--- a/test/freeflow/navierstokes/kovasznay/problem.hh
+++ b/test/freeflow/navierstokes/kovasznay/problem.hh
@@ -158,17 +158,6 @@ public:
     Scalar temperature() const
     { return 298.0; }
 
-
-   /*!
-     * \brief Returns the sources within the domain.
-     *
-     * \param globalPos The global position
-     */
-    NumEqVector sourceAtPos(const GlobalPosition &globalPos) const
-    {
-        return NumEqVector(0.0);
-    }
-
     // \}
    /*!
      * \name Boundary conditions
diff --git a/test/freeflow/navierstokesnc/channel/problem.hh b/test/freeflow/navierstokesnc/channel/problem.hh
index 6c771311413e2a0868f538ae1bd8b34185e6d78c..24d7ad2860f30bc89737de435280640e3f92e813 100644
--- a/test/freeflow/navierstokesnc/channel/problem.hh
+++ b/test/freeflow/navierstokesnc/channel/problem.hh
@@ -153,16 +153,6 @@ public:
     Scalar temperature() const
     { return 273.15 + 10; } // 10C
 
-   /*!
-     * \brief Returns the sources within the domain.
-     *
-     * \param globalPos The global position
-     */
-    NumEqVector sourceAtPos(const GlobalPosition &globalPos) const
-    {
-        return NumEqVector(0.0);
-    }
-
     // \}
    /*!
      * \name Boundary conditions
diff --git a/test/freeflow/navierstokesnc/densitydrivenflow/problem.hh b/test/freeflow/navierstokesnc/densitydrivenflow/problem.hh
index 0880641825581e1d197b1e1700ae40ff41aa2647..13465cc26c298ff9c301ce6af165169cc3273f7c 100644
--- a/test/freeflow/navierstokesnc/densitydrivenflow/problem.hh
+++ b/test/freeflow/navierstokesnc/densitydrivenflow/problem.hh
@@ -141,16 +141,6 @@ public:
     Scalar temperature() const
     { return 273.15 + 10; } // 10C
 
-   /*!
-     * \brief Returns the sources within the domain.
-     *
-     * \param globalPos The global position
-     */
-    NumEqVector sourceAtPos(const GlobalPosition &globalPos) const
-    {
-        return NumEqVector(0.0);
-    }
-
     // \}
    /*!
      * \name Boundary conditions
diff --git a/test/freeflow/navierstokesnc/maxwellstefan/problem.hh b/test/freeflow/navierstokesnc/maxwellstefan/problem.hh
index 70e632be2c8d0526326cbcffc5e9bfa00c308d18..1a3c1e9901829e4071cf57e881a5080fbf4b0a65 100644
--- a/test/freeflow/navierstokesnc/maxwellstefan/problem.hh
+++ b/test/freeflow/navierstokesnc/maxwellstefan/problem.hh
@@ -344,16 +344,6 @@ public:
     Scalar temperature() const
     { return 273.15 + 10; } // 10C
 
-   /*!
-     * \brief Returns the sources within the domain.
-     *
-     * \param globalPos The global position
-     */
-    NumEqVector sourceAtPos(const GlobalPosition &globalPos) const
-    {
-        return NumEqVector(0.0);
-    }
-
     // \}
    /*!
      * \name Boundary conditions
diff --git a/test/freeflow/rans/problem.hh b/test/freeflow/rans/problem.hh
index c15860dbcfc538d260c14c1b902e64009c1d4343..db7365642f6e408584fdf773ccc3aa526cfbe117 100644
--- a/test/freeflow/rans/problem.hh
+++ b/test/freeflow/rans/problem.hh
@@ -196,15 +196,6 @@ public:
     Scalar temperature() const
     { return inletTemperature_; }
 
-   /*!
-     * \brief Returns the sources within the domain.
-     *
-     * \param globalPos The global position
-     */
-    NumEqVector sourceAtPos(const GlobalPosition &globalPos) const
-    {
-        return NumEqVector(0.0);
-    }
     // \}
 
    /*!
diff --git a/test/freeflow/ransnc/problem.hh b/test/freeflow/ransnc/problem.hh
index 2b1cdf90d79c267e0798fcb022556a159580752f..be507c1f8635a128504f70384848790e35301c99 100644
--- a/test/freeflow/ransnc/problem.hh
+++ b/test/freeflow/ransnc/problem.hh
@@ -195,17 +195,6 @@ public:
     Scalar temperature() const
     { return 273.15 + 10; } // 10C
 
-
-   /*!
-     * \brief Returns the sources within the domain.
-     *
-     * \param globalPos The global position
-     */
-    NumEqVector sourceAtPos(const GlobalPosition &globalPos) const
-    {
-        return NumEqVector(0.0);
-    }
-
     // \}
    /*!
      * \name Boundary conditions