CD_EBHelmholtzEddingtonSP1DomainBC.H

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
/*
  @file   CD_EBHelmholtzEddingtonSP1DomainBC.H
  @brief  Declaration of a domain boundary condition class for EBHelmholtzOp.
  @author Robert Marskar
*/
 
#ifndef CD_EBHELMHOLTZEDDINGTONDOMAINBC_H
#define CD_EBHELMHOLTZEDDINGTONDOMAINBC_H
 
// Our includes
#include <CD_RtSpecies.H>
#include <CD_EddingtonSP1DomainBc.H>
#include <CD_EBHelmholtzDomainBC.H>
#include <CD_NamespaceHeader.H>
 
class EBHelmholtzEddingtonSP1DomainBC : public EBHelmholtzDomainBC
{
public:
  EBHelmholtzEddingtonSP1DomainBC() = delete;
 
  EBHelmholtzEddingtonSP1DomainBC(const EddingtonSP1DomainBc&     a_eddingtonBCs,
                                  const RefCountedPtr<RtSpecies>& a_species,
                                  const Real                      a_r1,
                                  const Real                      a_r2);
 
  EBHelmholtzEddingtonSP1DomainBC(const EBHelmholtzEddingtonSP1DomainBC& a_other) = delete;
 
  EBHelmholtzEddingtonSP1DomainBC(const EBHelmholtzEddingtonSP1DomainBC&& a_other) = delete;
 
  virtual ~EBHelmholtzEddingtonSP1DomainBC();
 
  EBHelmholtzEddingtonSP1DomainBC&
  operator=(const EBHelmholtzEddingtonSP1DomainBC& a_other) = delete;
 
  EBHelmholtzEddingtonSP1DomainBC&
  operator=(const EBHelmholtzEddingtonSP1DomainBC&& a_other) = delete;
 
  void
  define(const Location::Cell a_dataLocation,
         const EBLevelGrid&   a_eblg,
         const RealVect&      a_probLo,
         const Real           a_dx) override final;
 
  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 final;
 
  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 final;
 
  Real
  getDiagWeight(const int a_dir, const Side::LoHiSide a_side) const override final;
 
protected:
  EddingtonSP1DomainBc m_eddingtonBCs;
 
  RefCountedPtr<RtSpecies> m_species;
 
  Real m_r1;
 
  Real m_r2;
 
  std::map<EddingtonSP1DomainBc::DomainSide, std::shared_ptr<EBHelmholtzDomainBC>> m_bcObjects;
};
 
#include <CD_NamespaceFooter.H>
 
#endif