diff --git a/dumux/freeflow/staggered/fluxvariables.hh b/dumux/freeflow/staggered/fluxvariables.hh
index a344047526f322bddf33947bf4e79e5f798aebf2..30d48291e5d510f9898c7ceb8a8d87956a9c119a 100644
--- a/dumux/freeflow/staggered/fluxvariables.hh
+++ b/dumux/freeflow/staggered/fluxvariables.hh
@@ -33,7 +33,6 @@ namespace Properties
 {
 // forward declaration
 NEW_PROP_TAG(EnableComponentTransport);
-NEW_PROP_TAG(EnableEnergyBalanceStokes);
 NEW_PROP_TAG(EnableInertiaTerms);
 }
 
diff --git a/dumux/freeflow/staggered/localresidual.hh b/dumux/freeflow/staggered/localresidual.hh
index 54a5030441cde0c9044591ec7aa4e81a361aab58..25d58b08d94be9f4056f5b47f43f32553282252a 100644
--- a/dumux/freeflow/staggered/localresidual.hh
+++ b/dumux/freeflow/staggered/localresidual.hh
@@ -37,7 +37,6 @@ namespace Properties
 {
 // forward declaration
 NEW_PROP_TAG(EnableComponentTransport);
-NEW_PROP_TAG(EnableEnergyBalanceStokes);
 NEW_PROP_TAG(EnableInertiaTerms);
 NEW_PROP_TAG(ReplaceCompEqIdx);
 }
diff --git a/dumux/freeflow/staggered/propertydefaults.hh b/dumux/freeflow/staggered/propertydefaults.hh
index df8c3069f2c18431e09fd724949215f10f42072d..8259358bf997826a49a42298b5a8e3863a60396e 100644
--- a/dumux/freeflow/staggered/propertydefaults.hh
+++ b/dumux/freeflow/staggered/propertydefaults.hh
@@ -119,9 +119,6 @@ SET_BOOL_PROP(NavierStokes, EnableAdvection, true);
 //! The one-phase model has no molecular diffusion
 SET_BOOL_PROP(NavierStokes, EnableMolecularDiffusion, false);
 
-//! Isothermal model by default
-SET_BOOL_PROP(NavierStokes, EnableEnergyBalanceStokes, false);
-
 //! The indices required by the isothermal single-phase model
 SET_TYPE_PROP(NavierStokes, Indices, NavierStokesCommonIndices<TypeTag>);
 
diff --git a/dumux/implicit/properties.hh b/dumux/implicit/properties.hh
index ea9682eb0b05abe7494d8000424036e80f55e2c3..af294c60315c05106efaf6b4727a62796b5b8455 100644
--- a/dumux/implicit/properties.hh
+++ b/dumux/implicit/properties.hh
@@ -97,7 +97,6 @@ NEW_PROP_TAG(MolecularDiffusionType); //! The type for the calculation of the mo
 NEW_PROP_TAG(SolutionDependentMolecularDiffusion); //!< specifies if the parameters for the diffusive fluxes depend on the solution
 NEW_PROP_TAG(EnableEnergyBalance); //! Specifies if the model solves an energy equation
 NEW_PROP_TAG(HeatConductionType); //! The type for the calculation of the heat conduction fluxes
-NEW_PROP_TAG(EnableEnergyBalanceStokes); //! Specifies if the model solves an energy equation
 NEW_PROP_TAG(SolutionDependentHeatConduction); //!< specifies if the parameters for the heat conduction fluxes depend on the solution
 
 // specify if we evaluate the permeability in the volume (for discontinuous fields)