TriInterWrapperMesh.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 "InterWrapperMesh.h"

class TriInterWrapperMesh : public InterWrapperMesh
{
public:
  TriInterWrapperMesh(const InputParameters & parameters);
  TriInterWrapperMesh(const TriInterWrapperMesh & other_mesh);
  virtual std::unique_ptr<MooseMesh> safeClone() const override;
  virtual void buildMesh() override;

  virtual const unsigned int & getNumOfAssemblies() const override { return _n_assemblies; }

  virtual Node * getPinNode(unsigned int, unsigned) const override { return nullptr; }

  virtual bool pinMeshExist() const override { return false; }
  virtual bool ductMeshExist() const override { return false; }

  virtual const Real & getDuctToPinGap() const { return _duct_to_pin_gap; }

  virtual const unsigned int & getPinIndex(const unsigned int channel_idx,
                                           const unsigned int neighbor_idx)
  {
    return _chan_to_pin_map[channel_idx][neighbor_idx];
  }

  virtual Node * getChannelNode(unsigned int i_chan, unsigned iz) const override
  {
    return _nodes[i_chan][iz];
  }
  virtual const unsigned int & getNumOfChannels() const override { return _n_channels; }
  virtual const unsigned int & getNumOfGapsPerLayer() const override { return _n_gaps; }
  virtual const std::pair<unsigned int, unsigned int> &
  getGapChannels(unsigned int i_gap) const override
  {
    return _gap_to_chan_map[i_gap];
  }
  virtual const std::vector<unsigned int> & getChannelGaps(unsigned int i_chan) const override
  {
    return _chan_to_gap_map[i_chan];
  }
  virtual const std::vector<double> & getGapMap() const override { return _gij_map; }
  virtual const Real & getCrossflowSign(unsigned int i_chan, unsigned int i_local) const override
  {
    return _sign_id_crossflow_map[i_chan][i_local];
  }

  virtual unsigned int getSubchannelIndexFromPoint(const Point & p) const override;
  virtual unsigned int channelIndex(const Point & point) const override;

  virtual EChannelType getSubchannelType(unsigned int index) const override
  {
    return _subch_type[index];
  }

  virtual Real getGapWidth(unsigned int gap_index) const override { return _gij_map[gap_index]; }

  virtual const std::pair<unsigned int, unsigned int> & getSweepFlowGaps(unsigned int i_chan) const
  {
    return _gap_pairs_sf[i_chan];
  }

  virtual const std::pair<unsigned int, unsigned int> & getSweepFlowChans(unsigned int i_chan) const
  {
    return _chan_pairs_sf[i_chan];
  }

  virtual const std::vector<unsigned int> & getPinChannels(unsigned int i_pin) const override
  {
    return _pin_to_chan_map[i_pin];
  }

  virtual const std::vector<unsigned int> & getChannelPins(unsigned int i_chan) const override
  {
    return _chan_to_pin_map[i_chan];
  }

  virtual const bool & getIsTightSide() const { return _tight_side_bypass; }

  virtual unsigned int getPinIndexFromPoint(const Point & p) const override;
  virtual unsigned int pinIndex(const Point & p) const override;

  void setChannelToDuctMaps(const std::vector<Node *> & duct_nodes);

protected:
  unsigned int _n_rings;
  unsigned int _n_channels;
  Real _flat_to_flat;
  Real _duct_to_pin_gap;
  std::vector<std::vector<Node *>> _nodes;

  std::vector<Node *> _duct_nodes;
  std::map<Node *, Node *> _chan_to_duct_node_map;
  std::map<Node *, Node *> _duct_node_to_chan_map;

  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<Real> _gij_map;
  std::vector<std::vector<Real>> _subchannel_position;
  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 _n_assemblies;
  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;
  bool _tight_side_bypass;

public:
  static InputParameters validParams();

  static void rodPositions(std::vector<Point> & positions,
                           unsigned int nrings,
                           Real assembly_pitch,
                           Point center);

  friend class SCMTriInterWrapperMeshGenerator;

  static const unsigned int N_CORNERS = 6;
};