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();
    using std::abs;
 
    const auto delta_face = (*grad_phi_upwind) * (face_centroid - cell_centroid);
    const auto delta_max = abs(max_value - phi_upwind) + 1e-10;
    const auto delta_min = abs(min_value - phi_upwind) + 1e-10;
 
    const auto rf = delta_face >= 0 ? abs(delta_face) / delta_max : 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;
};
}
}