KokkosLinearFVFunctorDirichletBC.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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 "KokkosLinearFVBoundaryCondition.h"
#include "KokkosParsedFunction.h"

class KokkosLinearFVFunctorDirichletBC : public Moose::Kokkos::LinearFVBoundaryCondition
{
public:
  static InputParameters validParams();

  KokkosLinearFVFunctorDirichletBC(const InputParameters & parameters);

  template <typename Derived>
  KOKKOS_FUNCTION BoundaryRelation computeBoundaryValue(const FVDatum & datum) const;
  template <typename Derived>
  KOKKOS_FUNCTION BoundaryRelation computeBoundaryNormalGradient(const FVDatum & datum) const;

private:
  const Moose::Kokkos::ReferenceWrapper<const KokkosParsedFunction> _functor;
};

template <typename Derived>
KOKKOS_FUNCTION KokkosLinearFVFunctorDirichletBC::BoundaryRelation
KokkosLinearFVFunctorDirichletBC::computeBoundaryValue(const FVDatum & datum) const
{
  return {0, _functor->value(_t, datum.faceCentroid())};
}

template <typename Derived>
KOKKOS_FUNCTION KokkosLinearFVFunctorDirichletBC::BoundaryRelation
KokkosLinearFVFunctorDirichletBC::computeBoundaryNormalGradient(const FVDatum & datum) const
{
  const auto distance = datum.faceDCFMag();
  return {-1 / distance, _functor->value(_t, datum.faceCentroid()) / distance};
}