CD_EBHelmholtzEBBC.H

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
/*
  @file   CD_EBHelmholtzEBBC.H
  @brief  Declaration of a base boundary condition class for EB boundary conditions in EBHelmholtzOp
  @author Robert Marskar
*/
 
#ifndef CD_EBHELMHOLTZEBBC_H
#define CD_EBHELMHOLTZEBBC_H
 
// Chombo includes
#include <EBLevelGrid.H>
#include <RefCountedPtr.H>
#include <LevelData.H>
#include <BaseIVFAB.H>
 
// Our includes
#include <CD_Location.H>
#include <CD_NamespaceHeader.H>
 
class EBHelmholtzEBBC
{
public:
  using AmrMask = RefCountedPtr<LevelData<BaseFab<bool>>>;
 
  EBHelmholtzEBBC();
 
  virtual ~EBHelmholtzEBBC();
 
  EBHelmholtzEBBC(const EBHelmholtzEBBC& a_other) = delete;
 
  EBHelmholtzEBBC&
  operator=(const EBHelmholtzEBBC& a_other) = delete;
 
  void
  define(const Location::Cell a_dataLocation,
         const EBLevelGrid&   a_eblg,
         const AmrMask&       a_validCells,
         const RealVect&      a_probLo,
         const Real&          a_dx,
         const int            a_ghostCF);
 
  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 = 0;
 
  const LayoutData<BaseIVFAB<VoFStencil>>&
  getGradPhiStencils() const;
 
protected:
  constexpr static int m_comp = 0;
 
  constexpr static int m_nComp = 1;
 
  Location::Cell m_dataLocation;
 
  int m_ghostCF;
 
  Real m_dx;
 
  RealVect m_probLo;
 
  EBLevelGrid m_eblg;
 
  RefCountedPtr<LevelData<BaseFab<bool>>> m_validCells;
 
  LayoutData<BaseIVFAB<VoFStencil>> m_gradPhiStencils;
 
  virtual void
  define() = 0;
 
  inline bool
  isStencilValidCF(const VoFStencil& a_stencil, const DataIndex& a_dit) const;
 
  inline RealVect
  getBoundaryPosition(const VolIndex& a_vof, const DataIndex& a_dit) const;
 
  inline Real
  applyStencil(const VoFStencil& a_stencil, const EBCellFAB& a_phi) const;
};
 
#include <CD_NamespaceFooter.H>
 
#include <CD_EBHelmholtzEBBCImplem.H>
 
#endif