KokkosConvectiveHeatFluxBC.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 "KokkosIntegratedBC.h"

class KokkosConvectiveHeatFluxBC final
  : public Moose::Kokkos::IntegratedBC<KokkosConvectiveHeatFluxBC>
{
public:
  static InputParameters validParams();

  KokkosConvectiveHeatFluxBC(const InputParameters & parameters);

  KOKKOS_FUNCTION Real computeQpResidual(const unsigned int i,
                                         const unsigned int qp,
                                         ResidualDatum & datum) const;
  KOKKOS_FUNCTION Real computeQpJacobian(const unsigned int i,
                                         const unsigned int j,
                                         const unsigned int qp,
                                         ResidualDatum & datum) const;

private:
  Moose::Kokkos::MaterialProperty<Real> _T_infinity;

  Moose::Kokkos::MaterialProperty<Real> _htc;

  Moose::Kokkos::MaterialProperty<Real> _htc_dT;
};

KOKKOS_FUNCTION inline Real
KokkosConvectiveHeatFluxBC::computeQpResidual(const unsigned int i,
                                              const unsigned int qp,
                                              ResidualDatum & datum) const
{
  return -_test(datum, i, qp) * _htc(datum, qp) * (_T_infinity(datum, qp) - _u(datum, qp));
}

KOKKOS_FUNCTION inline Real
KokkosConvectiveHeatFluxBC::computeQpJacobian(const unsigned int i,
                                              const unsigned int j,
                                              const unsigned int qp,
                                              ResidualDatum & datum) const
{
  return -_test(datum, i, qp) * _phi(datum, j, qp) *
         (-_htc(datum, qp) + _htc_dT(datum, qp) * (_T_infinity(datum, qp) - _u(datum, qp)));
}