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];
}