CD_CdrCTU.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_CDRCTU_H
#define CD_CDRCTU_H
 
// Our includes
#include <CD_CdrMultigrid.H>
#include <CD_CdrGodunov.H>
#include <CD_NamespaceHeader.H>
 
class CdrCTU : public CdrMultigrid
{
public:
  CdrCTU();
 
  virtual ~CdrCTU();
 
  virtual void
  parseOptions() override;
 
  virtual void
  parseRuntimeOptions() override;
 
  virtual void
  advectToFaces(EBAMRFluxData& a_facePhi, const EBAMRCellData& a_cellPhi, const Real a_dt) override;
 
  virtual Real
  computeAdvectionDt() override;
 
  virtual void
  allocate() override;
 
  virtual void
  regrid(const int a_lbase, const int a_oldFinestLevel, const int a_newFinestLevel) override;
 
  virtual void
  preRegrid(const int a_lbase, const int a_oldFinestLevel) override;
 
protected:
  Vector<RefCountedPtr<LayoutData<VoFIterator>>> m_grownVofIter;
 
  void
  defineIterators();
 
  enum class Limiter
  {
    None,
    MinMod,
    Superbee,
    MonotonizedCentral,
  };
 
  Limiter m_limiter;
 
  bool m_useCTU;
 
  virtual void
  computeNormalSlopes(EBCellFAB&           a_normalSlopes,
                      const EBCellFAB&     a_cellPhi,
                      const Box&           a_cellBox,
                      const ProblemDomain& a_domain,
                      const int            a_level,
                      const DataIndex&     a_dit);
 
  virtual void
  upwind(EBFluxFAB&           a_facePhi,
         const EBCellFAB&     a_normalSlopes,
         const EBCellFAB&     a_cellPhi,
         const EBCellFAB&     a_cellVel,
         const EBFluxFAB&     a_faceVel,
         const ProblemDomain& a_domain,
         const Box&           a_cellBox,
         const int&           a_level,
         const DataIndex&     a_dit,
         const Real&          a_dt);
 
  virtual void
  parseSlopeLimiter();
 
  static Real
  minmod(const Real& dwl, const Real& dwr) noexcept;
 
  static Real
  superbee(const Real& dwl, const Real& dwr) noexcept;
 
  static Real
  monotonizedCentral(const Real& dwl, const Real& dwr) noexcept;
};
 
#include <CD_NamespaceFooter.H>
 
#endif