diff --git a/CHANGELOG.md b/CHANGELOG.md
index a57eb26ded2578de224f91e0dd136395c401b231..dc2833e9a296485d6ef6dbc63753c9a56c17ac66 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -24,6 +24,7 @@ v5.8.0
 - Make the CWIPI third-party patch work with IntelLLVM (!2053)
 - Use C++17 value template for traits (!2081)
 - Add H1-norm to Error filter (!1834)
+- Delete redundant lines in the Navier-Stokes diffusion functions (!2058)
 
 **CI**
 - Fix CubeAllElements performance test tolerance (!1943)
diff --git a/solvers/CompressibleFlowSolver/EquationSystems/NavierStokesCFE.cpp b/solvers/CompressibleFlowSolver/EquationSystems/NavierStokesCFE.cpp
index 0e390f7cb70e266f3cfc87b396b57c8ac9cbcd8b..9cf1ef68fdeafe41af2420437cf5873abd7eacad 100644
--- a/solvers/CompressibleFlowSolver/EquationSystems/NavierStokesCFE.cpp
+++ b/solvers/CompressibleFlowSolver/EquationSystems/NavierStokesCFE.cpp
@@ -220,10 +220,6 @@ void NavierStokesCFE::v_DoDiffusion(
             inBwd[i]       = Array<OneD, NekDouble>{nTracePts};
         }
 
-        // Extract pressure
-        // (use inarrayDiff[0] as a temporary storage for the pressure)
-        m_varConv->GetPressure(inarray, inarrayDiff[0]);
-
         // Extract temperature
         m_varConv->GetTemperature(inarray, inarrayDiff[nvariables - 2]);
 
@@ -239,9 +235,6 @@ void NavierStokesCFE::v_DoDiffusion(
         }
         else
         {
-            m_varConv->GetPressure(pFwd, inFwd[0]);
-            m_varConv->GetPressure(pBwd, inBwd[0]);
-
             m_varConv->GetTemperature(pFwd, inFwd[nvariables - 2]);
             m_varConv->GetTemperature(pBwd, inBwd[nvariables - 2]);
 
diff --git a/solvers/CompressibleFlowSolver/EquationSystems/NavierStokesImplicitCFE.cpp b/solvers/CompressibleFlowSolver/EquationSystems/NavierStokesImplicitCFE.cpp
index e8c9dc42e5fefc66ae7aee946cc05d6bee3f3df1..fcfbab59d4330b74076c12785e4425c5e282707b 100644
--- a/solvers/CompressibleFlowSolver/EquationSystems/NavierStokesImplicitCFE.cpp
+++ b/solvers/CompressibleFlowSolver/EquationSystems/NavierStokesImplicitCFE.cpp
@@ -164,10 +164,6 @@ void NavierStokesImplicitCFE::v_DoDiffusionCoeff(
             inBwd[i]       = Array<OneD, NekDouble>{nTracePts};
         }
 
-        // Extract pressure
-        //    (use inarrayDiff[0] as a temporary storage for the pressure)
-        m_varConv->GetPressure(inarray, inarrayDiff[0]);
-
         // Extract temperature
         m_varConv->GetTemperature(inarray, inarrayDiff[nvariables - 2]);
 
@@ -183,9 +179,6 @@ void NavierStokesImplicitCFE::v_DoDiffusionCoeff(
         }
         else
         {
-            m_varConv->GetPressure(pFwd, inFwd[0]);
-            m_varConv->GetPressure(pBwd, inBwd[0]);
-
             m_varConv->GetTemperature(pFwd, inFwd[nvariables - 2]);
             m_varConv->GetTemperature(pBwd, inBwd[nvariables - 2]);