TriSubChannelMesh.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 "SubChannelMesh.h"
 
#include <map>
#include <vector>
#include <memory>  // std::unique_ptr
#include <utility> // std::pair
 
class TriSubChannelMesh : public SubChannelMesh
{
public:
  TriSubChannelMesh(const InputParameters & parameters);
  TriSubChannelMesh(const TriSubChannelMesh & other_mesh);
 
  std::unique_ptr<MooseMesh> safeClone() const override;
  void buildMesh() override;
 
  void computeAssemblyHydraulicParameters();
 
  Real getSubchannelFlowArea(unsigned int i_chan, Real z) const override;
  Real getSubchannelWettedPerimeter(unsigned int i_chan) const override;
 
  unsigned int getNumOfPins() const override { return processor_id() == 0 ? _npins : 0; }
 
  Node * getPinNode(unsigned int i_pin, unsigned int iz) const override
  {
    return _pin_nodes[i_pin][iz];
  }
 
  const Real & getDuctToPinGap() const { return _duct_to_pin_gap; }
 
  const unsigned int & getNumOfRings() const { return _n_rings; }
 
  const unsigned int & getPinIndex(const unsigned int channel_idx, const unsigned int neighbor_idx)
  {
    return _chan_to_pin_map[channel_idx][neighbor_idx];
  }
 
  const Real & getWireDiameter() const { return _dwire; }
 
  const Real & getFlatToFlat() const { return _flat_to_flat; }
 
  const Real & getWireLeadLength() const { return _hwire; }
 
  Node * getChannelNode(unsigned int i_chan, unsigned int iz) const override
  {
    return _nodes[i_chan][iz];
  }
 
  unsigned int getNumOfChannels() const override { return processor_id() == 0 ? _n_channels : 0; }
 
  unsigned int getNumOfGapsPerLayer() const override { return processor_id() == 0 ? _n_gaps : 0; }
 
  const std::pair<unsigned int, unsigned int> & getGapChannels(unsigned int i_gap) const override
  {
    return _gap_to_chan_map[i_gap];
  }
 
  const std::pair<unsigned int, unsigned int> & getGapPins(unsigned int i_gap) const override
  {
    return _gap_to_pin_map[i_gap];
  }
 
  const std::vector<unsigned int> & getChannelGaps(unsigned int i_chan) const override
  {
    return _chan_to_gap_map[i_chan];
  }
 
  const Real & getCrossflowSign(unsigned int i_chan, unsigned int i_local) const override
  {
    return _sign_id_crossflow_map[i_chan][i_local];
  }
 
  unsigned int getSubchannelIndexFromPoint(const Point & p) const override;
  unsigned int channelIndex(const Point & point) const override;
 
  EChannelType getSubchannelType(unsigned int index) const override { return _subch_type[index]; }
 
  Real getGapWidth(unsigned int axial_index, unsigned int gap_index) const override
  {
    return _gij_map[axial_index][gap_index];
  }
 
  void setGapWidth(unsigned int axial_index, unsigned int gap_index, Real gap_width)
  {
    _gij_map[axial_index][gap_index] = gap_width;
  }
 
  const std::pair<unsigned int, unsigned int> & getSweepFlowGaps(unsigned int i_chan) const
  {
    return _gap_pairs_sf[i_chan];
  }
 
  const std::pair<unsigned int, unsigned int> & getSweepFlowChans(unsigned int i_chan) const
  {
    return _chan_pairs_sf[i_chan];
  }
 
  const std::vector<unsigned int> & getPinChannels(unsigned int i_pin) const override
  {
    return _pin_to_chan_map[i_pin];
  }
 
  const std::vector<unsigned int> & getChannelPins(unsigned int i_chan) const override
  {
    return _chan_to_pin_map[i_chan];
  }
 
  unsigned int getPinIndexFromPoint(const Point & p) const override;
  unsigned int pinIndex(const Point & p) const override;
 
protected:
  unsigned int _n_rings;
  unsigned int _n_channels;
  Real _flat_to_flat;
  Real _dwire;
  Real _hwire;
  Real _duct_to_pin_gap;
 
  std::vector<std::vector<Node *>> _nodes;
  std::vector<std::vector<Node *>> _pin_nodes;
 
  std::vector<std::pair<unsigned int, unsigned int>> _gap_to_chan_map;
  std::vector<std::pair<unsigned int, unsigned int>> _gap_to_pin_map;
  std::vector<std::vector<unsigned int>> _chan_to_gap_map;
 
  std::vector<std::vector<Real>> _sign_id_crossflow_map;
 
  std::vector<std::vector<Real>> _gij_map;
 
  std::vector<Point> _pin_position;
 
  std::vector<std::vector<Real>> _pins_in_rings;
 
  std::vector<std::vector<unsigned int>> _chan_to_pin_map;
 
  unsigned int _npins;
  unsigned int _n_gaps;
 
  std::vector<EChannelType> _subch_type;
  std::vector<EChannelType> _gap_type;
 
  std::vector<std::pair<unsigned int, unsigned int>> _gap_pairs_sf;
  std::vector<std::pair<unsigned int, unsigned int>> _chan_pairs_sf;
 
  std::vector<std::vector<unsigned int>> _pin_to_chan_map;
 
public:
  static InputParameters validParams();
 
  static void
  pinPositions(std::vector<Point> & positions, unsigned int nrings, Real pitch, Point center);
 
  friend class SCMTriAssemblyMeshGenerator;
  friend class SCMTriDuctMeshGenerator;
  friend class SCMDetailedTriAssemblyMeshGenerator;
 
  static const unsigned int N_CORNERS = 6;
};