BreakMeshByBlockGenerator.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"
#include <unordered_set>

/*
 * A mesh generator to split a mesh by a set of blocks
 */
class BreakMeshByBlockGenerator : public MeshGenerator
{
public:
  static InputParameters validParams();

  BreakMeshByBlockGenerator(const InputParameters & parameters);

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

private:
  void checkInputParameter();

  void findBoundaryName(const MeshBase & mesh,
                        subdomain_id_type primaryBlockID,
                        subdomain_id_type secondaryBlockID,
                        BoundaryName & boundary_name,
                        boundary_id_type boundaryID,
                        BoundaryInfo & boundary_info);

  std::set<std::pair<std::string, BoundaryID>> _bName_bID_set;

  BoundaryName generateBoundaryName(const MeshBase & mesh,
                                    subdomain_id_type primaryBlockID,
                                    subdomain_id_type secondaryBlockID);

  void mapBoundaryIdAndBoundaryName(boundary_id_type boundaryID, const std::string & boundaryName);

  void addInterface(MeshBase & mesh);

  std::unordered_map<dof_id_type, std::set<subdomain_id_type>> syncConnectedBlocks(
      const std::unordered_map<dof_id_type, std::vector<dof_id_type>> & node_to_elem_map,
      const MeshBase & mesh);

  subdomain_id_type blockRestrictedElementSubdomainID(const Elem * elem);

  bool findBlockPairs(subdomain_id_type block_one, subdomain_id_type block_two);

  std::unique_ptr<MeshBase> & _input;
  std::string _interface_name;
  bool _split_interface;
  std::unordered_set<std::pair<SubdomainID, SubdomainID>> _block_pairs;
  std::unordered_set<SubdomainID> _block_set;
  const bool _block_pairs_restricted;
  const bool _surrounding_blocks_restricted;
  const bool _add_transition_interface;
  const bool _split_transition_interface;
  const BoundaryName _interface_transition_name;
  const bool _add_interface_on_two_sides;

  using NodeConnectedBlocksPair = std::pair<dof_id_type, std::vector<subdomain_id_type>>;
  using SubdomainPair = std::pair<subdomain_id_type, subdomain_id_type>;
  using ElemSide = std::tuple<const Elem *, unsigned int>;

  std::set<SubdomainPair> _neighboring_block_list;

  std::unordered_map<SubdomainPair, std::set<ElemSide>> _subid_pairs_to_sides;

  std::unordered_map<std::pair<subdomain_id_type, subdomain_id_type>, boundary_id_type>
      _subid_pairs_to_boundary_id;
};