SubChannelMesh.h

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

#pragma once

#include <vector>
#include "MooseMesh.h"
#include "SubChannelEnums.h"

class SubChannelMesh : public MooseMesh
{
public:
  SubChannelMesh(const InputParameters & parameters);
  SubChannelMesh(const SubChannelMesh & other_mesh);

  virtual const std::vector<Real> & getZGrid() const { return _z_grid; }

  virtual unsigned int getZIndex(const Point & point) const;

  virtual const std::vector<std::vector<Real>> & getKGrid() const { return _k_grid; }

  virtual const std::vector<Real> & getZBlockage() const { return _z_blockage; }

  virtual const std::vector<unsigned int> & getIndexBlockage() const { return _index_blockage; }

  virtual const std::vector<Real> & getReductionBlockage() const { return _reduction_blockage; }

  virtual const Real & getKij() const { return _kij; }

  virtual const unsigned int & getNumOfAxialCells() const { return _n_cells; }

  virtual Node * getChannelNode(unsigned int i_chan, unsigned iz) const = 0;

  virtual Node * getPinNode(unsigned int i_pin, unsigned iz) const = 0;

  virtual Node * getChanNodeFromDuct(Node * duct_node) = 0;

  virtual Node * getChannelNodeFromDuct(Node * channel_node) = 0;

  virtual Node * getDuctNodeFromChannel(Node * channel_node) = 0;

  virtual const unsigned int & getNumOfChannels() const = 0;

  virtual bool pinMeshExist() const = 0;

  virtual bool ductMeshExist() const = 0;

  virtual const unsigned int & getNumOfGapsPerLayer() const = 0;

  virtual const unsigned int & getNumOfPins() const = 0;

  virtual const std::pair<unsigned int, unsigned int> &
  getGapChannels(unsigned int i_gap) const = 0;

  virtual const std::pair<unsigned int, unsigned int> & getGapPins(unsigned int i_gap) const = 0;

  virtual const std::vector<unsigned int> & getChannelGaps(unsigned int i_chan) const = 0;

  virtual const std::vector<unsigned int> & getPinChannels(unsigned int i_pin) const = 0;

  virtual const std::vector<unsigned int> & getChannelPins(unsigned int i_chan) const = 0;

  virtual const Real & getPitch() const { return _pitch; }

  virtual const Real & getPinDiameter() const { return _pin_diameter; }

  virtual const Real & getCrossflowSign(unsigned int i_chan, unsigned int i_local) const = 0;

  virtual const Real & getHeatedLengthEntry() const { return _unheated_length_entry; }

  virtual const Real & getHeatedLength() const { return _heated_length; }

  virtual const Real & getHeatedLengthExit() const { return _unheated_length_exit; }

  virtual unsigned int getSubchannelIndexFromPoint(const Point & p) const = 0;

  virtual unsigned int channelIndex(const Point & point) const = 0;

  virtual unsigned int getPinIndexFromPoint(const Point & p) const = 0;

  virtual unsigned int pinIndex(const Point & p) const = 0;

  virtual EChannelType getSubchannelType(unsigned int index) const = 0;

  virtual Real getGapWidth(unsigned int axial_index, unsigned int gap_index) const = 0;

  virtual const std::vector<Node *> getDuctNodes() const = 0;

protected:
  Real _unheated_length_entry;
  Real _heated_length;
  Real _unheated_length_exit;
  std::vector<Real> _z_grid;
  std::vector<std::vector<Real>> _k_grid;
  std::vector<Real> _spacer_z;
  std::vector<Real> _spacer_k;
  std::vector<Real> _z_blockage;
  std::vector<unsigned int> _index_blockage;
  std::vector<Real> _reduction_blockage;
  Real _kij;
  Real _pitch;
  Real _pin_diameter;
  unsigned int _n_cells;

public:
  std::vector<std::vector<Real>> _subchannel_position;
  static InputParameters validParams();

  static void generateZGrid(Real unheated_length_entry,
                            Real heated_length,
                            Real unheated_length_exit,
                            unsigned int n_cells,
                            std::vector<Real> & z_grid);
};