CD_MFHelmholtzEBBC.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_MFHELMHOLTZEBBC_H
#define CD_MFHELMHOLTZEBBC_H
 
// Our includes
#include <CD_EBHelmholtzEBBC.H>
#include <CD_MFHelmholtzJumpBC.H>
#include <CD_VofUtils.H>
#include <CD_NamespaceHeader.H>
 
class MFHelmholtzEBBC : public EBHelmholtzEBBC
{
public:
  MFHelmholtzEBBC() = delete;
 
  MFHelmholtzEBBC(const MFHelmholtzEBBC& a_other) = delete;
 
  MFHelmholtzEBBC(const int a_phase, const RefCountedPtr<MFHelmholtzJumpBC>& a_jumpBC);
 
  virtual ~MFHelmholtzEBBC();
 
  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 override;
 
protected:
  RefCountedPtr<MFHelmholtzJumpBC> m_jumpBC;
 
  int m_phase;
 
  LayoutData<BaseIVFAB<Real>> m_boundaryWeights;
 
  void
  define() override final;
 
  void
  defineMultiPhase();
 
  virtual void
  defineSinglePhase() = 0;
 
  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 = 0;
 
  void
  applyEBFluxMultiPhase(VoFIterator&           a_multiPhaseVofs,
                        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;
 
  virtual bool
  getLeastSquaresBoundaryGradStencil(std::pair<Real, VoFStencil>& a_stencil,
                                     const VolIndex&              a_vof,
                                     const VofUtils::Neighborhood a_neighborhood,
                                     const DataIndex&             a_dit,
                                     const int                    a_order,
                                     const int                    a_weight) const;
};
 
#include <CD_NamespaceFooter.H>
 
#endif