CD_EBHelmholtzRobinDomainBC.H

/* chombo-discharge
 * Copyright © 2021 SINTEF Energy Research.
 * Please refer to Copyright.txt and LICENSE in the chombo-discharge root directory.
 */
 
/*
  @file   CD_EBHelmholtzRobinDomainBC.H
  @brief  Declaration of a Robin boundary condition class for EBHelmholtzOp
  @author Robert Marskar
*/
 
#ifndef CD_EBHelmholtzRobinDomainBC_H
#define CD_EBHelmholtzRobinDomainBC_H
 
// Std includes
#include <functional>
 
// Our includes
#include <CD_EBHelmholtzDomainBC.H>
#include <CD_NamespaceHeader.H>
 
class EBHelmholtzRobinDomainBC : public EBHelmholtzDomainBC
{
public:
  EBHelmholtzRobinDomainBC();
 
  virtual ~EBHelmholtzRobinDomainBC();
 
  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
  getFaceFlux(BaseFab<Real>&        a_faceFlux,
              const BaseFab<Real>&  a_phi,
              const BaseFab<Real>&  a_Bcoef,
              const int&            a_dir,
              const Side::LoHiSide& a_side,
              const DataIndex&      a_dit,
              const bool            a_useHomogeneous) const override;
 
  virtual Real
  getFaceFlux(const VolIndex&       a_vof,
              const EBCellFAB&      a_phi,
              const EBFaceFAB&      a_Bcoef,
              const int&            a_dir,
              const Side::LoHiSide& a_side,
              const DataIndex&      a_dit,
              const bool            a_useHomogeneous) const override;
 
protected:
  bool m_useConstant;
 
  bool m_useFunction;
 
  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;
};
 
#include <CD_NamespaceFooter.H>
 
#endif