MarvelTriFlowAreaIC.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 "MarvelTriFlowAreaIC.h"
#include "TriSubChannelMesh.h"
#include "SCM.h"

registerMooseObject("SubChannelApp", MarvelTriFlowAreaIC);

InputParameters
MarvelTriFlowAreaIC::validParams()
{
  InputParameters params = TriSubChannelBaseIC::validParams();
  params.addClassDescription("Computes flow area of subchannels in a triangular lattice "
                             "arrangement in a MARVEL-type micro-reactor");
  return params;
}

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

Real
MarvelTriFlowAreaIC::value(const Point & p)
{
  Real standard_area, rod_area, wire_area, additional_area;
  auto pitch = _mesh.getPitch();
  auto rod_diameter = _mesh.getPinDiameter();
  auto wire_diameter = _mesh.getWireDiameter();
  auto wire_lead_length = _mesh.getWireLeadLength();
  auto gap = _mesh.getDuctToPinGap();
  auto r_ref = rod_diameter / 2.0 + gap;
  auto z_blockage = _mesh.getZBlockage();
  auto index_blockage = _mesh.getIndexBlockage();
  auto reduction_blockage = _mesh.getReductionBlockage();
  auto theta = std::acos(wire_lead_length /
                         std::sqrt(std::pow(wire_lead_length, 2) +
                                   std::pow(libMesh::pi * (rod_diameter + wire_diameter), 2)));
  auto i = _mesh.getSubchannelIndexFromPoint(p);
  // given the channel number, i, it computes the flow area of the subchannel
  // based on the subchannel type: CENTER, EDGE or CORNER.
  auto subch_type = _mesh.getSubchannelType(i);

  if (subch_type == EChannelType::CENTER)
  {
    standard_area = std::pow(pitch, 2) * std::sqrt(3.0) / 4.0;
    rod_area = libMesh::pi * std::pow(rod_diameter, 2.0) / 8.0;
    additional_area = 0.0;
    wire_area = libMesh::pi * std::pow(wire_diameter, 2.0) / 8.0 / std::cos(theta);
  }
  else if (subch_type == EChannelType::EDGE)
  {
    auto gamma = std::acos(1 - 0.5 * std::pow(pitch / r_ref, 2.0));
    auto alpha = (libMesh::pi - gamma) / 2.0;
    auto sector_angle = libMesh::pi / 2.0 - alpha;
    auto triangle_height = std::sqrt(std::pow(r_ref, 2.0) - std::pow(pitch / 2.0, 2.0));
    auto sector = 0.5 * sector_angle * std::pow(r_ref, 2.0);
    auto triangle = 0.5 * pitch * triangle_height;
    standard_area = 2.0 * sector + triangle;
    rod_area = libMesh::pi * std::pow(rod_diameter, 2.0) / 8.0;
    additional_area = 0.0;
    wire_area = libMesh::pi * std::pow(wire_diameter, 2.0) / 8.0 / std::cos(theta);
  }
  else
  {
    standard_area = (libMesh::pi / 6.0) * std::pow(r_ref, 2.0);
    rod_area = libMesh::pi * std::pow(rod_diameter, 2.0) / 24.0;
    additional_area = 0.0;
    wire_area = libMesh::pi * std::pow(wire_diameter, 2.0) / 24.0 / std::cos(theta);
  }

  auto subchannel_area = standard_area + additional_area - rod_area - wire_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;
}