SCMQuadFlowAreaIC.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "SCMQuadFlowAreaIC.h"

#include "SubChannelMesh.h"
#include "QuadSubChannelMesh.h"
#include "SCM.h"

registerMooseObject("SubChannelApp", SCMQuadFlowAreaIC);
registerMooseObjectRenamed("SubChannelApp", QuadFlowAreaIC, "06/30/2025 24:00", SCMQuadFlowAreaIC);

InputParameters
SCMQuadFlowAreaIC::validParams()
{
  InputParameters params = QuadSubChannelBaseIC::validParams();
  params.addClassDescription(
      "Computes subchannel flow area in the square lattice subchannel arrangement");
  return params;
}

SCMQuadFlowAreaIC::SCMQuadFlowAreaIC(const InputParameters & params)
  : QuadSubChannelBaseIC(params), _subchannel_mesh(SCM::getConstMesh<SubChannelMesh>(_mesh))
{
}

Real
SCMQuadFlowAreaIC::value(const Point & p)
{
  Real standard_area, rod_area, additional_area;
  auto pitch = _mesh.getPitch();
  auto pin_diameter = _mesh.getPinDiameter();
  auto gap = _mesh.getGap();
  auto z_blockage = _mesh.getZBlockage();
  auto index_blockage = _mesh.getIndexBlockage();
  auto reduction_blockage = _mesh.getReductionBlockage();
  auto i = _mesh.getSubchannelIndexFromPoint(p);
  auto subch_type = _mesh.getSubchannelType(i);

  if (subch_type == EChannelType::CORNER)
  {
    standard_area = 0.25 * pitch * pitch;
    rod_area = 0.25 * 0.25 * M_PI * pin_diameter * pin_diameter;
    additional_area = pitch * gap + gap * gap;
  }
  else if (subch_type == EChannelType::EDGE)
  {
    standard_area = 0.5 * pitch * pitch;
    rod_area = 0.5 * 0.25 * M_PI * pin_diameter * pin_diameter;
    additional_area = pitch * gap;
  }
  else
  {
    standard_area = pitch * pitch;
    rod_area = 0.25 * M_PI * pin_diameter * pin_diameter;
    additional_area = 0.0;
  }

  auto subchannel_area = standard_area + additional_area - rod_area;

  auto index = 0;
  for (const auto & i_blockage : index_blockage)
  {
    if (i == i_blockage && (p(2) >= z_blockage.front() && p(2) <= z_blockage.back()))
    {
      return reduction_blockage[index] * subchannel_area;
    }
    index++;
  }

  return subchannel_area;
}