Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "QuadInterWrapper1PhaseProblem.h"
      11             : 
      12             : registerMooseObject("SubChannelApp", QuadInterWrapper1PhaseProblem);
      13             : 
      14             : InputParameters
      15          66 : QuadInterWrapper1PhaseProblem::validParams()
      16             : {
      17          66 :   InputParameters params = InterWrapper1PhaseProblem::validParams();
      18          66 :   params.addClassDescription(
      19             :       "Solver class for interwrapper of assemblies in a square-lattice arrangement");
      20          66 :   return params;
      21           0 : }
      22             : 
      23          33 : QuadInterWrapper1PhaseProblem::QuadInterWrapper1PhaseProblem(const InputParameters & params)
      24          33 :   : InterWrapper1PhaseProblem(params)
      25             : {
      26          33 : }
      27             : 
      28             : double
      29     1405512 : QuadInterWrapper1PhaseProblem::computeFrictionFactor(Real Re)
      30             : {
      31             :   Real a, b;
      32     1405512 :   if (Re < 1)
      33             :   {
      34             :     return 64.0;
      35             :   }
      36     1405512 :   else if (Re >= 1 and Re < 5000)
      37             :   {
      38             :     a = 64.0;
      39             :     b = -1.0;
      40             :   }
      41      639000 :   else if (Re >= 5000 and Re < 30000)
      42             :   {
      43             :     a = 0.316;
      44             :     b = -0.25;
      45             :   }
      46             :   else
      47             :   {
      48             :     a = 0.184;
      49             :     b = -0.20;
      50             :   }
      51     1405512 :   return a * std::pow(Re, b);
      52             : }