geometry_disk.h

Go to the documentation of this file.

//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
// [LANL Copyright Statement]
// Copyright (C) 2024 by the ryujin authors
// Copyright (C) 2024 by Triad National Security, LLC
//
 
#pragma once
 
#include "geometry_common_includes.h"
 
namespace ryujin
{
  namespace Geometries
  {
    template <int dim>
    class Disk : public Geometry<dim>
    {
    public:
      Disk(const std::string subsection)
          : Geometry<dim>("disk", subsection)
      {
        balanced_ = true;
        this->add_parameter("balanced",
                            balanced_,
                            "Use GridGenerator::hyper_ball_balanced() instead "
                            "of the older GridGenerator::hyper_ball()");
 
        radius_ = 1.2;
        this->add_parameter("radius", radius_, "radius of disk");
 
        boundary_ = Boundary::dirichlet;
        this->add_parameter("boundary condition",
                            boundary_,
                            "Type of boundary condition enforced on the "
                            "boundary of the disk/ball");
      }
 
      void create_triangulation(
          typename Geometry<dim>::Triangulation &triangulation) final
      {
        if (balanced_) {
          GridGenerator::hyper_ball_balanced(
              triangulation, dealii::Point<dim>(), radius_);
        } else {
          GridGenerator::hyper_ball(
              triangulation, dealii::Point<dim>(), radius_);
        }
 
        /*
         * Set boundary ids:
         */
        for (auto cell : triangulation.active_cell_iterators()) {
          for (auto f : cell->face_indices()) {
            const auto face = cell->face(f);
 
            if (!face->at_boundary())
              continue;
 
            face->set_boundary_id(boundary_);
          }
        }
      }
 
    private:
      bool balanced_;
      double radius_;
      Boundary boundary_;
    };
  } /* namespace Geometries */
} /* namespace ryujin */