@@ -2392,7 +2392,7 @@ void CFVMFlowSolverBase<V, FlowRegime>::Friction_Forces(const CGeometry* geometr
23922392
23932393 unsigned long iVertex, iPoint, iPointNormal;
23942394 unsigned short iMarker, iMarker_Monitoring, iDim, jDim;
2395- unsigned short T_INDEX = 0 , TVE_INDEX = 0 , VEL_INDEX = 0 , RHOS_INDEX = 0 ;
2395+ unsigned short T_INDEX = 0 , TVE_INDEX = 0 , VEL_INDEX = 0 ;
23962396 su2double Viscosity = 0.0 , WallDist[3 ] = {0.0 }, Area, Density = 0.0 , GradTemperature = 0.0 ,
23972397 UnitNormal[3 ] = {0.0 }, TauElem[3 ] = {0.0 }, TauTangent[3 ] = {0.0 }, Tau[3 ][3 ] = {{0.0 }}, Cp,
23982398 thermal_conductivity, MaxNorm = 8.0 , Grad_Vel[3 ][3 ] = {{0.0 }}, Grad_Temp[3 ] = {0.0 }, AxiFactor;
@@ -2572,25 +2572,16 @@ void CFVMFlowSolverBase<V, FlowRegime>::Friction_Forces(const CGeometry* geometr
25722572
25732573 } else {
25742574
2575- unsigned short iSpecies, nSpecies = config->GetnSpecies ();
2576-
25772575 const auto & thermal_conductivity_tr = nodes->GetThermalConductivity (iPoint);
25782576 const auto & thermal_conductivity_ve = nodes->GetThermalConductivity_ve (iPoint);
25792577 const auto & Grad_PrimVar = nodes->GetGradient_Primitive (iPoint);
2580- const auto & Ds = nodes->GetDiffusionCoeff (iPoint);
2581- const auto & hs = nodes->GetEnthalpys (iPoint);
2582-
2583- /* --- Compute enthalpy transport to surface due to mass diffusion ---*/
2584- su2double sumJhs = 0.0 ;
2585- for (iSpecies = 0 ; iSpecies < nSpecies; iSpecies++)
2586- sumJhs += Ds[iSpecies]*hs[iSpecies]*GeometryToolbox::DotProduct (nDim, Grad_PrimVar[RHOS_INDEX+iSpecies], UnitNormal);
25872578
25882579 su2double dTn = GeometryToolbox::DotProduct (nDim, Grad_PrimVar[T_INDEX], UnitNormal);
25892580 su2double dTven = GeometryToolbox::DotProduct (nDim, Grad_PrimVar[TVE_INDEX], UnitNormal);
25902581
25912582 /* --- Surface energy balance: trans-rot heat flux, vib-el heat flux,
25922583 enthalpy transport due to mass diffusion ---*/
2593- HeatFlux[iMarker][iVertex] = thermal_conductivity_tr*dTn + thermal_conductivity_ve*dTven + sumJhs ;
2584+ HeatFlux[iMarker][iVertex] = thermal_conductivity_tr*dTn + thermal_conductivity_ve*dTven;
25942585 }
25952586
25962587 /* --- Note that y+, and heat are computed at the
0 commit comments