@@ -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 , Area, Density = 0.0 , GradTemperature = 0.0 , WallDistMod, FrictionVel,
23972397 UnitNormal[3 ] = {0.0 }, TauElem[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;
@@ -2577,25 +2577,16 @@ void CFVMFlowSolverBase<V, FlowRegime>::Friction_Forces(const CGeometry* geometr
25772577
25782578 } else {
25792579
2580- unsigned short iSpecies, nSpecies = config->GetnSpecies ();
2581-
25822580 const auto & thermal_conductivity_tr = nodes->GetThermalConductivity (iPoint);
25832581 const auto & thermal_conductivity_ve = nodes->GetThermalConductivity_ve (iPoint);
25842582 const auto & Grad_PrimVar = nodes->GetGradient_Primitive (iPoint);
2585- const auto & Ds = nodes->GetDiffusionCoeff (iPoint);
2586- const auto & hs = nodes->GetEnthalpys (iPoint);
2587-
2588- /* --- Compute enthalpy transport to surface due to mass diffusion ---*/
2589- su2double sumJhs = 0.0 ;
2590- for (iSpecies = 0 ; iSpecies < nSpecies; iSpecies++)
2591- sumJhs += Ds[iSpecies]*hs[iSpecies]*GeometryToolbox::DotProduct (nDim, Grad_PrimVar[RHOS_INDEX+iSpecies], UnitNormal);
25922583
25932584 su2double dTn = GeometryToolbox::DotProduct (nDim, Grad_PrimVar[T_INDEX], UnitNormal);
25942585 su2double dTven = GeometryToolbox::DotProduct (nDim, Grad_PrimVar[TVE_INDEX], UnitNormal);
25952586
25962587 /* --- Surface energy balance: trans-rot heat flux, vib-el heat flux,
25972588 enthalpy transport due to mass diffusion ---*/
2598- HeatFlux[iMarker][iVertex] = thermal_conductivity_tr*dTn + thermal_conductivity_ve*dTven + sumJhs ;
2589+ HeatFlux[iMarker][iVertex] = thermal_conductivity_tr*dTn + thermal_conductivity_ve*dTven;
25992590 }
26002591
26012592 /* --- Note that y+, and heat are computed at the
0 commit comments