VenkatakrishnanLimiter.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "Limiter.h"
#include "MathFVUtils.h"

namespace Moose
{
namespace FV
{
template <typename T>
class VenkatakrishnanLimiter : public Limiter<T>
{
public:
  T limit(const T & phi_upwind,
          const T & /* phi_downwind */,
          const VectorValue<T> * grad_phi_upwind,
          const VectorValue<T> * /* grad_phi_downwind */,
          const RealVectorValue & /* dCD */,
          const Real & max_value,
          const Real & min_value,
          const FaceInfo * fi,
          const bool & fi_elem_is_upwind) const override final
  {
    const auto face_centroid = fi->faceCentroid();
    const auto cell_centroid = fi_elem_is_upwind ? fi->elemCentroid() : fi->neighborCentroid();

    const auto delta_face = (*grad_phi_upwind) * (face_centroid - cell_centroid);
    const auto delta_max = std::abs(max_value - phi_upwind) + 1e-10;
    const auto delta_min = std::abs(min_value - phi_upwind) + 1e-10;

    const auto rf =
        delta_face >= 0 ? std::abs(delta_face) / delta_max : std::abs(delta_face) / delta_min;

    return (2 * rf + 1.0) / (rf * (2 * rf + 1.0) + 1.0);
  }

  bool constant() const override final { return false; }

  InterpMethod interpMethod() const override final { return InterpMethod::Venkatakrishnan; }

  VenkatakrishnanLimiter() = default;
};
}
}