KokkosCoupledConvectiveHeatFluxBC.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 "KokkosIntegratedBCValue.h"

class KokkosCoupledConvectiveHeatFluxBC : public Moose::Kokkos::IntegratedBCValue
{
public:
  static InputParameters validParams();

  KokkosCoupledConvectiveHeatFluxBC(const InputParameters & parameters);

  template <typename Derived>
  KOKKOS_FUNCTION Real computeQpResidual(const unsigned int qp, AssemblyDatum & datum) const;
  template <typename Derived>
  KOKKOS_FUNCTION Real computeQpJacobian(const unsigned int j,
                                         const unsigned int qp,
                                         AssemblyDatum & datum) const;

private:
  const unsigned int _n_components;
  const Moose::Kokkos::VariableValue _T_infinity;
  const Moose::Kokkos::VariableValue _htc;
  const Moose::Kokkos::VariableValue _alpha;
  const Moose::Kokkos::VariableValue _scale_factor;
};

template <typename Derived>
KOKKOS_FUNCTION Real
KokkosCoupledConvectiveHeatFluxBC::computeQpResidual(const unsigned int qp,
                                                     AssemblyDatum & datum) const
{
  Real q = 0;
  Real u = _u(datum, qp);
  for (unsigned int c = 0; c < _n_components; c++)
    q += _alpha(datum, qp, c) * _htc(datum, qp, c) * (u - _T_infinity(datum, qp, c));
  return q * _scale_factor(datum, qp);
}

template <typename Derived>
KOKKOS_FUNCTION Real
KokkosCoupledConvectiveHeatFluxBC::computeQpJacobian(const unsigned int j,
                                                     const unsigned int qp,
                                                     AssemblyDatum & datum) const
{
  Real dq = 0;
  Real phi = _phi(datum, j, qp);
  for (unsigned int c = 0; c < _n_components; c++)
    dq += _alpha(datum, qp, c) * _htc(datum, qp, c) * phi;
  return dq * _scale_factor(datum, qp);
}