CircularBoundaryCorrectionGenerator.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 "MooseEnum.h"
#include "MooseMeshUtils.h"

class CircularBoundaryCorrectionGenerator : public MeshGenerator
{
public:
  static InputParameters validParams();

  CircularBoundaryCorrectionGenerator(const InputParameters & parameters);

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

protected:
  const MeshGeneratorName _input_name;
  const std::vector<BoundaryName> _input_mesh_circular_boundaries;
  const std::vector<Real> _transition_layer_ratios;
  const std::vector<Real> _custom_circular_tolerance;
  std::vector<boundary_id_type> _input_mesh_circular_bids;
  const bool _move_end_nodes_in_span_direction;
  std::unique_ptr<MeshBase> & _input;

  Point circularCenterCalculator(const std::vector<Point> & pts_list,
                                 Real & radius,
                                 const Real tol = 1e-12) const;

  Real generateRadialCorrectionFactor(const std::vector<std::pair<Point, Point>> & bd_side_list,
                                      const Point & circle_center,
                                      const bool is_closed_loop,
                                      const bool move_end_nodes_in_span_direction,
                                      Real & c_coeff,
                                      Real & end_node_disp) const;

  Real sineSummation(const std::vector<Real> & th_list, const Real coeff) const;

  Real sinePrimeSummation(const std::vector<Real> & th_list, const Real coeff) const;
};