CD_EBHelmholtzRobinDomainBC.H

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
/*
  @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