CD_MFHelmholtzRobinEBBC.H

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
/*
  @file   CD_MFHelmholtzRobinEBBC.H
  @brief  Declaration of a boundary condition class for MFHelmholtzOp
  @author Robert Marskar
*/
 
#ifndef CD_MFHELMHOLTZROBINEBBC_H
#define CD_MFHELMHOLTZROBINEBBC_H
 
// Std includes
#include <functional>
 
// Our includes
#include <CD_VofUtils.H>
#include <CD_MFHelmholtzEBBC.H>
#include <CD_NamespaceHeader.H>
 
class MFHelmholtzRobinEBBC : public MFHelmholtzEBBC
{
public:
  MFHelmholtzRobinEBBC(const int a_phase, const RefCountedPtr<MFHelmholtzJumpBC>& a_jumpBC);
 
  virtual ~MFHelmholtzRobinEBBC();
 
  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);
 
protected:
  /*
    @brief Use constant for BC
  */
  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
  defineSinglePhase() override;
 
  virtual void
  applyEBFluxSinglePhase(VoFIterator&           a_singlePhaseVofs,
                         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;
 
  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