boundary_conditions.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/exceptions.h>
 
#include <array>
 
namespace MeltPoolDG::SpeciesTransport
{
 
  enum class BoundaryConditionType
  {
    dirichlet,
    neumann
  };
 
  template <std::size_t n_partial_densities,
            typename VectorizedArrayType,
            typename ReadOnlyDofView,
            typename WritableDofView>
  void
  set_boundary_value_and_gradient(
    const WritableDofView                                      &w_p,
    const ReadOnlyDofView                                      &w_m,
    const BoundaryConditionType                                 boundary_type,
    const std::array<VectorizedArrayType, n_partial_densities> &dirichlet_boundary_values = {})
  {
    const auto set_to_value = [](const ReadOnlyDofView     &src,
                                 const WritableDofView     &dst,
                                 const VectorizedArrayType &scaling_factor = 1.) {
      for (unsigned int species = 0; species < n_partial_densities; ++species)
        {
          dst.partial_density(species) = src.partial_density(species) * scaling_factor;
        }
    };
 
    const auto set_to_gradient = [](const ReadOnlyDofView     &src,
                                    const WritableDofView     &dst,
                                    const VectorizedArrayType &scaling_factor = 1.) {
      for (unsigned int species = 0; species < n_partial_densities; ++species)
        dst.grad_partial_density(species) = src.grad_partial_density(species) * scaling_factor;
    };
 
    if (boundary_type == BoundaryConditionType::dirichlet)
      {
        for (unsigned int species = 0; species < n_partial_densities; ++species)
          w_p.partial_density(species) = dirichlet_boundary_values[species];
 
        set_to_gradient(w_m, w_p);
      }
    else if (boundary_type == BoundaryConditionType::neumann)
      {
        set_to_value(w_m, w_p);
        set_to_gradient(w_m, w_p, -1.);
      }
    else
      {
        AssertThrow(false, dealii::ExcMessage("Unknown boundary type for species transport."));
      }
  }
} // namespace MeltPoolDG::SpeciesTransport