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 "QuadInterWrapperWettedPerimIC.h"
      11             : #include "QuadInterWrapperMesh.h"
      12             : #include "QuadSubChannelBaseIC.h"
      13             : 
      14             : registerMooseObject("SubChannelApp", QuadInterWrapperWettedPerimIC);
      15             : 
      16             : InputParameters
      17          68 : QuadInterWrapperWettedPerimIC::validParams()
      18             : {
      19          68 :   InputParameters params = QuadSubChannelBaseIC::validParams();
      20          68 :   params.addClassDescription("Computes wetted perimeter of inter-wrapper cells in the square "
      21             :                              "lattice subchannel arrangement");
      22          68 :   return params;
      23           0 : }
      24             : 
      25          35 : QuadInterWrapperWettedPerimIC::QuadInterWrapperWettedPerimIC(const InputParameters & params)
      26          35 :   : QuadInterWrapperBaseIC(params)
      27             : {
      28          35 : }
      29             : 
      30             : Real
      31       17280 : QuadInterWrapperWettedPerimIC::value(const Point & p)
      32             : {
      33       17280 :   auto side_x = _mesh.getSideX();
      34       17280 :   auto side_y = _mesh.getSideY();
      35       17280 :   auto square_perimeter = 2.0 * side_x + 2.0 * side_y;
      36       17280 :   auto i = _mesh.getSubchannelIndexFromPoint(p);
      37       17280 :   auto subch_type = _mesh.getSubchannelType(i);
      38             : 
      39       17280 :   if (subch_type == EChannelType::CORNER)
      40             :     // Need to consider specific side - will add later
      41        1920 :     return 0.25 * square_perimeter;
      42       15360 :   else if (subch_type == EChannelType::EDGE)
      43             :     // Need to consider specific side - will add later
      44        7680 :     return 0.75 * square_perimeter;
      45             :   else
      46             :     return square_perimeter;
      47             : }