Nondimensionalize Python custom temperatures correctly

j-signorelli · j-signorelli · commit 84f9227ea92c · 2023-07-11T19:06:58.000-05:00
diff --git a/SU2_CFD/src/solvers/CHeatSolver.cpp b/SU2_CFD/src/solvers/CHeatSolver.cpp
@@ -377,7 +377,7 @@ void CHeatSolver::BC_Isothermal_Wall(CGeometry *geometry, CSolver **solver_conta
   SU2_OMP_FOR_STAT(OMP_MIN_SIZE)
   for (auto iVertex = 0ul; iVertex < geometry->nVertex[val_marker]; iVertex++) {
     if (IsPyCustom) {
-      Twall = geometry->GetCustomBoundaryTemperature(val_marker, iVertex);
+      Twall = geometry->GetCustomBoundaryTemperature(val_marker, iVertex) / config->GetTemperature_Ref();
     }
     IsothermalBoundaryCondition(geometry, solver_container[FLOW_SOL], config, val_marker, iVertex, Twall);
   }
diff --git a/SU2_CFD/src/solvers/CNSSolver.cpp b/SU2_CFD/src/solvers/CNSSolver.cpp
@@ -710,7 +710,7 @@ void CNSSolver::BC_Isothermal_Wall_Generic(CGeometry *geometry, CSolver **solver
                                     thermal_conductivity, There, Temperature_Ref);
     }
     else if (config->GetMarker_All_PyCustom(val_marker)) {
-      Twall = geometry->GetCustomBoundaryTemperature(val_marker, iVertex);
+      Twall = geometry->GetCustomBoundaryTemperature(val_marker, iVertex) / Temperature_Ref;
     }
 
     /*--- Compute the normal gradient in temperature using Twall ---*/

Original file line number	Diff line number	Diff line change
`@@ -377,7 +377,7 @@ void CHeatSolver::BC_Isothermal_Wall(CGeometry geometry, CSolver *solver_conta`
`377`	`377`	`SU2_OMP_FOR_STAT(OMP_MIN_SIZE)`
`378`	`378`	`for (auto iVertex = 0ul; iVertex < geometry->nVertex[val_marker]; iVertex++) {`
`379`	`379`	`if (IsPyCustom) {`
`380`		`- Twall = geometry->GetCustomBoundaryTemperature(val_marker, iVertex);`
	`380`	`+ Twall = geometry->GetCustomBoundaryTemperature(val_marker, iVertex) / config->GetTemperature_Ref();`
`381`	`381`	`}`
`382`	`382`	`IsothermalBoundaryCondition(geometry, solver_container[FLOW_SOL], config, val_marker, iVertex, Twall);`
`383`	`383`	`}`
Original file line number	Diff line number	Diff line change
`@@ -710,7 +710,7 @@ void CNSSolver::BC_Isothermal_Wall_Generic(CGeometry geometry, CSolver *solver`
`710`	`710`	`thermal_conductivity, There, Temperature_Ref);`
`711`	`711`	`}`
`712`	`712`	`else if (config->GetMarker_All_PyCustom(val_marker)) {`
`713`		`- Twall = geometry->GetCustomBoundaryTemperature(val_marker, iVertex);`
	`713`	`+ Twall = geometry->GetCustomBoundaryTemperature(val_marker, iVertex) / Temperature_Ref;`
`714`	`714`	`}`
`715`	`715`
`716`	`716`	`/--- Compute the normal gradient in temperature using Twall ---/`