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;
};