electrostatic_configuration.h

Go to the documentation of this file.

//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
// Copyright (C) 2026 by the ryujin authors
//
 
#pragma once
 
#include <ryujin/base/compile_time_options.h>
#include <ryujin/base/convenience_macros.h>
#include <ryujin/discretization/discretization.h>
 
#include <deal.II/base/exceptions.h>
#include <deal.II/base/parameter_acceptor.h>
#include <deal.II/base/tensor.h>
 
#include <string>
 
namespace ryujin
{
  template <int dim, typename Number = double>
  class ElectrostaticConfiguration : public dealii::ParameterAcceptor
  {
  public:
    using curl_type = dealii::Tensor<1, dim == 2 ? 1 : dim, Number>;
 
    ElectrostaticConfiguration(const std::string &name,
                               const std::string &subsection)
        : ParameterAcceptor(subsection + "/" + name)
        , name_(name)
    {
      dirichlet_boundaries_.insert({Boundary::do_nothing,
                                    Boundary::slip,
                                    Boundary::no_slip,
                                    Boundary::dirichlet,
                                    Boundary::dynamic,
                                    Boundary::dirichlet_momentum,
                                    Boundary::dirichlet_velocity});
      this->add_parameter(
          "dirichlet boundaries",
          dirichlet_boundaries_,
          "A list of hyperbolic boundary types where homogeneous boundary "
          "conditions will be enforced on the potential.");
 
      is_time_dependent_ = false;
    }
 
    virtual double background_density(const dealii::Point<dim> &point,
                                      Number t) const = 0;
 
    virtual curl_type magnetic_field(const dealii::Point<dim> &point,
                                     Number t) const = 0;
 
    ACCESSOR_READ_ONLY(name)
 
    
    ACCESSOR_READ_ONLY(dirichlet_boundaries)
 
    
    ACCESSOR_READ_ONLY(is_time_dependent)
 
  protected:
    bool is_time_dependent_;
 
  private:
    const std::string name_;
    std::set<dealii::types::boundary_id> dirichlet_boundaries_;
  };
} // namespace ryujin