Skip to content

Commit 21844bd

Browse files
authored
fix error in derivative term(dyp_dup) of spalding wall function
1 parent a0e588d commit 21844bd

File tree

1 file changed

+1
-1
lines changed

1 file changed

+1
-1
lines changed

SU2_CFD/src/solvers/CIncNSSolver.cpp

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -752,7 +752,7 @@ void CIncNSSolver::SetTau_Wall_WF(CGeometry *geometry, CSolver **solver_containe
752752

753753
/* --- Gradient of function defined above wrt U_Tau --- */
754754

755-
const su2double dyp_dup = 1.0 + exp(-kappa * B) * (kappa * exp(kUp) - kappa - kUp - 0.5 * kUp * kUp);
755+
const su2double dyp_dup = 1.0 + kappa * exp(-kappa * B) * (exp(kUp) - 1.0 - kUp - 0.5 * kUp * kUp);
756756
const su2double dup_dutau = - U_Plus / U_Tau;
757757
const su2double grad_diff = Density_Wall * WallDistMod / Lam_Visc_Wall - dyp_dup * dup_dutau;
758758

0 commit comments

Comments
 (0)