diff --git a/ProcessLib/TwoPhaseFlowWithPrho/TwoPhaseFlowWithPrhoMaterialProperties.cpp b/ProcessLib/TwoPhaseFlowWithPrho/TwoPhaseFlowWithPrhoMaterialProperties.cpp
index a0d7db9596dbaf6e06673eaeb7401c0381c7dbfe..2a4ae2113db9c31b75afc50983e9b43e8918f2fa 100644
--- a/ProcessLib/TwoPhaseFlowWithPrho/TwoPhaseFlowWithPrhoMaterialProperties.cpp
+++ b/ProcessLib/TwoPhaseFlowWithPrho/TwoPhaseFlowWithPrhoMaterialProperties.cpp
@@ -98,16 +98,6 @@ double TwoPhaseFlowWithPrhoMaterialProperties::getGasDensity(
     return _gas_density->getValue(vars);
 }
 
-double TwoPhaseFlowWithPrhoMaterialProperties::getDerivGasDensity(
-    const double p, const double T) const
-{
-    ArrayType vars;
-    vars[static_cast<int>(MaterialLib::Fluid::PropertyVariableType::T)] = T;
-    vars[static_cast<int>(MaterialLib::Fluid::PropertyVariableType::p)] = p;
-
-    return _gas_density->getdValue(vars,
-                                   MaterialLib::Fluid::PropertyVariableType::p);
-}
 double TwoPhaseFlowWithPrhoMaterialProperties::getLiquidViscosity(
     const double p, const double T) const
 {
diff --git a/ProcessLib/TwoPhaseFlowWithPrho/TwoPhaseFlowWithPrhoMaterialProperties.h b/ProcessLib/TwoPhaseFlowWithPrho/TwoPhaseFlowWithPrhoMaterialProperties.h
index 13c4687ca07c97f1a1c8f45d036fd3bf0ff6963c..66a7cb9befa59c489a7e7408efca3a281ad601fc 100644
--- a/ProcessLib/TwoPhaseFlowWithPrho/TwoPhaseFlowWithPrhoMaterialProperties.h
+++ b/ProcessLib/TwoPhaseFlowWithPrho/TwoPhaseFlowWithPrhoMaterialProperties.h
@@ -89,7 +89,6 @@ public:
     double getGasDensity(const double p, const double T) const;
     double getGasViscosity(const double p, const double T) const;
     double getLiquidViscosity(const double p, const double T) const;
-    double getDerivGasDensity(double const p, double const T) const;
     bool computeConstitutiveRelation(
         double const t,
         ParameterLib::SpatialPosition const& x_position,