CutMeshByLevelSetGeneratorBase.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 "FunctionParserUtils.h"

class CutMeshByLevelSetGeneratorBase : public MeshGenerator, public FunctionParserUtils<false>
{
public:
  static InputParameters validParams();

  CutMeshByLevelSetGeneratorBase(const InputParameters & parameters);

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

  enum class PointLevelSetRelationIndex : short int
  {
    level_set_out_side = 1,
    level_set_in_side = -1,
    on_level_set = 0
  };

protected:
  const MeshGeneratorName _input_name;
  boundary_id_type _cut_face_id;
  const BoundaryName _cut_face_name;
  std::unique_ptr<MeshBase> & _input;
  SymFunctionPtr _func_level_set;

  PointLevelSetRelationIndex pointLevelSetRelation(const Point & point);

  Point pointPairLevelSetInterception(const Point & point1, const Point & point2);

  void tet4ElemCutter(ReplicatedMesh & mesh,
                      const std::vector<libMesh::BoundaryInfo::BCTuple> & bdry_side_list,
                      const dof_id_type elem_id,
                      const subdomain_id_type & block_id_to_remove,
                      std::vector<const Node *> & new_on_plane_nodes);

  const Node * nonDuplicateNodeCreator(ReplicatedMesh & mesh,
                                       std::vector<const Node *> & new_on_plane_nodes,
                                       const Point & new_point) const;

  Real levelSetEvaluator(const Point & point);
};