AzimuthalBlockSplitGenerator.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 "PolygonMeshGeneratorBase.h"
#include "MooseEnum.h"

class AzimuthalBlockSplitGenerator : public PolygonMeshGeneratorBase
{
public:
  static InputParameters validParams();

  AzimuthalBlockSplitGenerator(const InputParameters & parameters);

  std::unique_ptr<MeshBase> generate() override;

protected:
  const MeshGeneratorName _input_name;
  std::vector<subdomain_id_type> _old_block_ids;
  const std::vector<subdomain_id_type> _new_block_ids;
  const std::vector<SubdomainName> _new_block_names;
  const bool _preserve_volumes;
  Real _start_angle;
  const Real _angle_range;
  std::vector<Real> & _azimuthal_angle_meta;
  std::unique_ptr<MeshBase> & _input;
  std::vector<unsigned int> _num_sectors_per_side_meta;
  Real _end_angle;

  void angleIdentifier(const Real & terminal_angle,
                       Real & original_down,
                       std::vector<Real>::iterator & original_down_it,
                       Real & original_up,
                       std::vector<Real>::iterator & original_up_it,
                       std::vector<Real> & azimuthal_angles_vtx);

  void angleModifier(const Real & side_angular_shift,
                     const Real & side_angular_range,
                     const Real & azi_tol,
                     const Real & terminal_angle,
                     const Real & original_down,
                     std::vector<Real>::iterator & original_down_it,
                     const Real & original_up,
                     std::vector<Real>::iterator & original_up_it,
                     Real & azi_to_keep,
                     Real & azi_to_mod);

  void nodeModifier(ReplicatedMesh & mesh,
                    const std::vector<std::pair<Real, dof_id_type>> & node_id_mod,
                    const std::vector<std::pair<Real, dof_id_type>> & node_id_keep,
                    std::vector<Real> & circular_rad_list,
                    std::vector<Real> & non_circular_rad_list,
                    const std::vector<std::tuple<dof_id_type, boundary_id_type>> & node_list,
                    const Real & term_angle,
                    const bool & external_block_change,
                    const Real & rad_tol);

  Point midPointCorrector(const Point vertex_0,
                          const Point vertex_1,
                          const Real max_radius,
                          const Real rad_tol);
};