NekMeshGenerator.h

Go to the documentation of this file.

/********************************************************************/
/*                  SOFTWARE COPYRIGHT NOTIFICATION                 */
/*                             Cardinal                             */
/*                                                                  */
/*                  (c) 2021 UChicago Argonne, LLC                  */
/*                        ALL RIGHTS RESERVED                       */
/*                                                                  */
/*                 Prepared by UChicago Argonne, LLC                */
/*               Under Contract No. DE-AC02-06CH11357               */
/*                With the U. S. Department of Energy               */
/*                                                                  */
/*             Prepared by Battelle Energy Alliance, LLC            */
/*               Under Contract No. DE-AC07-05ID14517               */
/*                With the U. S. Department of Energy               */
/*                                                                  */
/*                 See LICENSE for full restrictions                */
/********************************************************************/
 
#pragma once
 
#include "MeshGenerator.h"
 
class NekMeshGenerator : public MeshGenerator
{
public:
  static InputParameters validParams();
 
  NekMeshGenerator(const InputParameters & parameters);
 
  std::unique_ptr<MeshBase> generate() override;
 
 
  void initializeElemData(std::unique_ptr<MeshBase> & mesh);
 
  void checkElementType(std::unique_ptr<MeshBase> & mesh);
 
  virtual std::pair<unsigned int, unsigned int> pairedNodesAboutMidPoint(const unsigned int & node_id) const;
 
  unsigned int getFaceNode(const unsigned int & primary_face) const;
 
  virtual bool isCornerNode(const unsigned int & node) const;
 
  void moveNodes(std::unique_ptr<MeshBase> & mesh, std::vector<Real> & polygon_layer_smoothing);
 
  const Elem * getNextLayerElem(const Elem & elem, const unsigned int & touching_face,
    unsigned int & next_touching_face) const;
 
  unsigned int midPointNodeIndex(const unsigned int & face_id, const unsigned int & face_node) const;
 
  unsigned int pairedFaceNode(const unsigned int & node_id, const unsigned int & face_id) const;
 
  std::vector<Elem *> getBoundaryLayerElems(Elem * elem, const unsigned int & n_layers,
    const unsigned int & primary_face) const;
 
  void moveElem(Elem * elem, const unsigned int & boundary_index, const unsigned int & primary_face,
    const std::vector<Real> & polygon_layer_smoothing);
 
  bool isNearCorner(const Point & pt) const;
 
  unsigned int getNodeIndex(const Elem * elem, const Point & pt) const;
 
  Point getClosestOrigin(const unsigned int & index, const Point & pt) const;
 
  void checkPointLength(const std::vector<std::vector<Real>> & points, std::string name) const;
 
  void adjustMidPointNode(const unsigned int & node_id, Elem * elem) const;
 
  BoundaryID getBoundaryID(const BoundaryName & name, const MeshBase & mesh) const;
 
  Point projectPoint(const Point & origin, const Point & pt) const;
 
  Point adjustPointToCircle(const unsigned int & node_id, Elem * elem, const Real & radius, const Point & origin) const;
 
protected:
  std::unique_ptr<MeshBase> & _input;
 
  const MooseEnum & _geometry_type;
 
  const MooseEnum & _axis;
 
  const bool & _curve_corners;
 
  const Real & _rotation_angle;
 
  const bool & _retain_original_elem_type;
 
  const bool _has_moving_boundary;
 
  unsigned int _n_start_nodes;
 
  unsigned int _n_end_nodes;
 
  unsigned int _n_start_nodes_per_side;
 
  unsigned int _n_sides;
 
  unsigned int _n_corner_nodes;
 
  std::vector<BoundaryID> _moving_boundary;
 
  std::vector<Real> _radius;
 
  std::vector<std::vector<Real>> _origin;
 
  std::vector<unsigned int> _layers;
 
  Real _corner_radius;
 
  std::vector<Point> _polygon_corners;
 
  Real _max_corner_distance;
 
  unsigned int _n_noncorner_boundaries;
 
  std::set<BoundaryID> _boundaries_to_rebuild;
 
  std::map<BoundaryID, BoundaryName> _boundary_id_to_name;
 
  std::vector<std::vector<unsigned int>> _side_ids;
 
  std::vector<std::vector<unsigned int>> _corner_nodes;
 
  std::vector<std::vector<std::pair<unsigned int, unsigned int>>> _across_pair;
 
  std::vector<unsigned int> _across_face;
 
  ElemType _etype;
 
  std::vector<std::vector<unsigned int>> _side_nodes_map;
 
  std::vector<std::vector<unsigned int>> _face_nodes_map;
};
class NekMeshGenerator : public MeshGenerator {…};