Update kEpsilonViscosityAux.C idaholab#27800 idaholab#27887 idaholab#…

…27888

modules/navier_stokes/src/auxkernels/kEpsilonViscosityAux.C

@@ -117,15 +117,14 @@ kEpsilonViscosityAux::computeValue()
 
     // Switch for determining the near wall quantities
     // wall_treatment can be: "eq_newton eq_incremental eq_linearized neq"
-    ADReal u_tau;
     ADReal y_plus;
     ADReal mut_log; // turbulent log-layer viscosity
     ADReal mu_wall; // total wall viscosity to obtain the shear stress at the wall
 
     if (_wall_treatment == "eq_newton")
     {
       // Full Newton-Raphson solve to find the wall quantities from the law of the wall
-      u_tau = NS::findUStar(mu, rho, parallel_speed, min_wall_dist);
+      const auto u_tau = NS::findUStar(mu, rho, parallel_speed, min_wall_dist);
       y_plus = min_wall_dist * u_tau * rho / mu;
       mu_wall = rho * Utility::pow<2>(u_tau) * min_wall_dist / parallel_speed;
       mut_log = mu_wall - mu;
@@ -146,7 +145,7 @@ kEpsilonViscosityAux::computeValue()
                          (std::log(NS::E_turb_constant * std::max(min_wall_dist, 1.0) / mu) + 1.0);
       const ADReal c_c = parallel_speed;
 
-      u_tau = (-b_c + std::sqrt(std::pow(b_c, 2) + 4.0 * a_c * c_c)) / (2.0 * a_c);
+      const auto u_tau = (-b_c + std::sqrt(std::pow(b_c, 2) + 4.0 * a_c * c_c)) / (2.0 * a_c);
       y_plus = min_wall_dist * u_tau * rho / mu;
       mu_wall = rho * Utility::pow<2>(u_tau) * min_wall_dist / parallel_speed;
       mut_log = mu_wall - mu;