PINSFVRhieChowInterpolator.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 "INSFVRhieChowInterpolator.h"
#include "CellCenteredMapFunctor.h"
#include <unordered_map>
 
class PINSFVRhieChowInterpolator : public INSFVRhieChowInterpolator
{
public:
  static InputParameters validParams();
  PINSFVRhieChowInterpolator(const InputParameters & params);
 
  virtual void meshChanged() override;
  virtual void initialize() override;
  virtual void execute() override;
  virtual void finalize() override;
 
protected:
  const Moose::FunctorBase<ADReal> & epsilon(THREAD_ID tid) const override;
 
  bool isFaceGeometricallyRelevant(const FaceInfo & fi) const;
 
  const Moose::Functor<ADReal> & _eps;
 
  CellCenteredMapFunctor<ADReal, std::unordered_map<dof_id_type, ADReal>> _smoothed_eps;
 
  std::vector<const Moose::Functor<ADReal> *> _epss;
 
  std::vector<const Moose::Functor<ADReal> *> _smoothed_epss;
 
  const unsigned short _smoothing_layers;
 
  std::vector<const FaceInfo *> _geometric_fi;
 
private:
  void pinsfvSetup();
};
 
inline const Moose::FunctorBase<ADReal> &
PINSFVRhieChowInterpolator::epsilon(const THREAD_ID tid) const
{
  if (!_smoothing_layers)
    return *_epss[tid];
  else
    return *_smoothed_epss[tid];
}