SCMTriDuctMeshGenerator.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 "MeshGenerator.h"

class SCMTriDuctMeshGenerator : public MeshGenerator
{
public:
  SCMTriDuctMeshGenerator(const InputParameters & parameters);
  static InputParameters validParams();
  std::unique_ptr<MeshBase> generate() override;

protected:
  size_t
  ductPointIndex(unsigned int points_per_layer, unsigned int layer, unsigned int point) const;

  void ductCorners(std::vector<Point> & corners, Real flat_to_flat, const Point & center) const;

  void ductCrossSec(std::vector<Point> & cross_sec,
                    const std::vector<Point> & corners,
                    unsigned int nrings,
                    Real pitch) const;

  void ductPoints(std::vector<Point> & points,
                  const std::vector<Point> & cross_sec,
                  const std::vector<Real> & z_layers) const;

  void ductElems(std::vector<std::vector<size_t>> & elem_point_indices,
                 unsigned int n_layers,
                 unsigned int points_per_layer) const;

  void buildDuct(std::unique_ptr<MeshBase> & mesh,
                 std::vector<Node *> & duct_nodes,
                 const std::vector<Point> & points,
                 const std::vector<std::vector<size_t>> & elem_point_indices,
                 SubdomainID block) const;

  std::unique_ptr<MeshBase> & _input;
  const unsigned int _n_cells;
  std::vector<Real> _z_grid;
  const Real _unheated_length_entry;
  const Real _heated_length;
  const Real _unheated_length_exit;
  const unsigned int _block_id;
  const Real _pitch;
  const unsigned int _n_rings;
  const Real _flat_to_flat;
};