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 precomputeQpResidual(const unsigned int qp, AssemblyDatum & datum) const;
  template <typename Derived>
  KOKKOS_FUNCTION Real precomputeQpJacobian(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::precomputeQpResidual(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::precomputeQpJacobian(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);
}