Merge pull request #575 from pecos-hybrid/bugfix_viscosity_in_FD

Change dynamic viscosity in FD hybrid central / upwind blending to kinematic viscosity
su2code · Sep 6, 2018 · d00c3b3 · d00c3b3
2 parents a0ecddf + f3c985d
commit d00c3b3
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/SU2_CFD/src/variable_direct_mean.cpp b/SU2_CFD/src/variable_direct_mean.cpp
@@ -677,7 +677,7 @@ void CNSVariable::SetRoe_Dissipation_FD(su2double val_wall_dist){
 
   static const su2double k2 = pow(0.41,2.0);
 
-  su2double uijuij = 0, r_d;
+  su2double uijuij = 0;
   unsigned short iDim, jDim;
 
   AD::StartPreacc();
@@ -697,7 +697,9 @@ void CNSVariable::SetRoe_Dissipation_FD(su2double val_wall_dist){
   uijuij=sqrt(fabs(uijuij));
   uijuij=max(uijuij,1e-10);
 
-  r_d = (GetEddyViscosity()+GetLaminarViscosity())/(uijuij*k2*pow(val_wall_dist, 2.0));
+  const su2double nu = GetLaminarViscosity()/GetDensity();
+  const su2double nu_t = GetEddyViscosity()/GetDensity();
+  const su2double r_d = (nu + nu_t)/(uijuij*k2*pow(val_wall_dist, 2.0));
 
   Roe_Dissipation = 1.0-tanh(pow(8.0*r_d,3.0));