One-way coupling of UR factor between flow and turb

bmunguia · bmunguia · commit 10e9e0ecc844 · 2020-08-03T11:12:07.000-07:00
diff --git a/SU2_CFD/src/solvers/CTurbSolver.cpp b/SU2_CFD/src/solvers/CTurbSolver.cpp
@@ -682,6 +682,10 @@ void CTurbSolver::ComputeUnderRelaxationFactor(CSolver **solver_container, const
       }
     }
 
+    /* Choose the minimum factor between mean flow and turbulence. */
+
+    localUnderRelaxation = min(localUnderRelaxation, solver_container[FLOW_SOL]->GetNodes()->GetUnderRelaxation(iPoint));
+
     /* Threshold the relaxation factor in the event that there is
      a very small value. This helps avoid catastrophic crashes due
      to non-realizable states by canceling the update. */

Original file line number	Diff line number	Diff line change
`@@ -682,6 +682,10 @@ void CTurbSolver::ComputeUnderRelaxationFactor(CSolver **solver_container, const`
`682`	`682`	`}`
`683`	`683`	`}`
`684`	`684`
	`685`	`+ /* Choose the minimum factor between mean flow and turbulence. */`
	`686`	`+`
	`687`	`+ localUnderRelaxation = min(localUnderRelaxation, solver_container[FLOW_SOL]->GetNodes()->GetUnderRelaxation(iPoint));`
	`688`	`+`
`685`	`689`	`/* Threshold the relaxation factor in the event that there is`
`686`	`690`	`a very small value. This helps avoid catastrophic crashes due`
`687`	`691`	`to non-realizable states by canceling the update. */`