sa-neg fix

pcarruscag · pcarruscag · commit bd4021020147 · 2022-02-27T12:29:15.000Z
diff --git a/SU2_CFD/include/numerics/turbulent/turb_sources.hpp b/SU2_CFD/include/numerics/turbulent/turb_sources.hpp
@@ -418,7 +418,8 @@ struct Neg {
 
   static void ComputeDestruction(const su2double& nue, const CSAVariables& var, su2double& destruction,
                                  su2double& jacobian) {
-    const su2double dD_dnu = var.cw1 * nue / var.dist_i_2;
+    /*--- The destruction when nue < 0 is added instead of the usual subtraction, hence the negative sign. ---*/
+    const su2double dD_dnu = -var.cw1 * nue / var.dist_i_2;
     destruction = dD_dnu * nue;
     jacobian -= 2 * dD_dnu;
   }

Original file line number	Diff line number	Diff line change
`@@ -418,7 +418,8 @@ struct Neg {`
`418`	`418`
`419`	`419`	`static void ComputeDestruction(const su2double& nue, const CSAVariables& var, su2double& destruction,`
`420`	`420`	`su2double& jacobian) {`
`421`		`- const su2double dD_dnu = var.cw1 * nue / var.dist_i_2;`
	`421`	`+ /--- The destruction when nue < 0 is added instead of the usual subtraction, hence the negative sign. ---/`
	`422`	`+ const su2double dD_dnu = -var.cw1 * nue / var.dist_i_2;`
`422`	`423`	`destruction = dD_dnu * nue;`
`423`	`424`	`jacobian -= 2 * dD_dnu;`
`424`	`425`	`}`