CD_EBHelmholtzRobinEBBC.H

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
/*
  @file   CD_EBHelmholtzRobinEBBC.H
  @brief  Declaration of a boundary condition class for EBHelmholtzOp
  @author Robert Marskar
*/
 
#ifndef CD_EBHELMHOLTZROBINEBBC_H
#define CD_EBHELMHOLTZROBINEBBC_H
 
// Std includes
#include <functional>
 
// Our includes
#include <CD_VofUtils.H>
#include <CD_EBHelmholtzEBBC.H>
#include <CD_NamespaceHeader.H>
 
class EBHelmholtzRobinEBBC : public EBHelmholtzEBBC
{
public:
  EBHelmholtzRobinEBBC();
 
  EBHelmholtzRobinEBBC(const int a_order, const int a_weight, const Real a_A, const Real a_B, const Real a_C);
 
  EBHelmholtzRobinEBBC(const int                                         a_order,
                       const int                                         a_weight,
                       const std::function<Real(const RealVect& a_pos)>& a_A,
                       const std::function<Real(const RealVect& a_pos)>& a_B,
                       const std::function<Real(const RealVect& a_pos)>& a_C);
 
  virtual ~EBHelmholtzRobinEBBC();
 
  void
  setOrder(const int a_order);
 
  void
  setWeight(const int a_weight);
 
  virtual void
  setDomainDropOrder(const int a_domainSize);
 
  void
  setCoefficients(const Real a_A, const Real a_B, const Real a_C);
 
  void
  setCoefficients(const std::function<Real(const RealVect& a_pos)>& a_A,
                  const std::function<Real(const RealVect& a_pos)>& a_B,
                  const std::function<Real(const RealVect& a_pos)>& a_C);
 
  virtual void
  applyEBFlux(VoFIterator&           a_vofit,
              EBCellFAB&             a_Lphi,
              const EBCellFAB&       a_phi,
              const BaseIVFAB<Real>& a_Bcoef,
              const DataIndex&       a_dit,
              const Real&            a_beta,
              const bool&            a_homogeneousPhysBC) const override;
 
protected:
  bool m_useConstant;
 
  bool m_useFunction;
 
  int m_order;
 
  int m_weight;
 
  int m_domainDropOrder;
 
  Real m_constantA;
 
  Real m_constantB;
 
  Real m_constantC;
 
  std::function<Real(const RealVect& a_pos)> m_functionA;
 
  std::function<Real(const RealVect& a_pos)> m_functionB;
 
  std::function<Real(const RealVect& a_pos)> m_functionC;
 
  virtual void
  define() override;
 
  VoFStencil
  getInterpolationStencil(const VolIndex&              a_vof,
                          const DataIndex&             a_dit,
                          const VofUtils::Neighborhood a_neighborhood,
                          const int                    a_order) const;
};
 
#include <CD_NamespaceFooter.H>
 
#endif