CD_EBHelmholtzNeumannDomainBC.H

Go to the documentation of this file.

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
#ifndef CD_EBHELMHOLTZNEUMANNDOMAINBC_H
#define CD_EBHELMHOLTZNEUMANNDOMAINBC_H
 
// Std includes
#include <functional>
 
// Our includes
#include <CD_EBHelmholtzDomainBC.H>
#include <CD_NamespaceHeader.H>
 
class EBHelmholtzNeumannDomainBC : public EBHelmholtzDomainBC
{
public:
  EBHelmholtzNeumannDomainBC();
 
  EBHelmholtzNeumannDomainBC(const Real a_DphiDn);
 
  EBHelmholtzNeumannDomainBC(const std::function<Real(const RealVect& a_pos)>& a_DphiDn);
 
  virtual ~EBHelmholtzNeumannDomainBC();
 
  void
  setDphiDn(const Real a_DphiDn);
 
  void
  setDphiDn(const std::function<Real(const RealVect& a_pos)>& a_DphiDn);
 
  void
  setBxDphiDn(const Real a_BxDphiDn);
 
  void
  setBxDphiDn(const std::function<Real(const RealVect& a_pos)>& a_BxDphiDn);
 
  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;
 
  Real
  getDiagWeight(const int a_dir, const Side::LoHiSide a_side) const override;
 
protected:
  bool m_multByBco;
 
  bool m_useConstant;
 
  bool m_useFunction;
 
  Real m_constantDphiDn;
 
  std::function<Real(const RealVect& a_pos)> m_functionDphiDn;
};
 
#include <CD_NamespaceFooter.H>
 
#endif