KokkosHeatConductionTimeDerivative.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 "KokkosTimeDerivative.h"
 
class KokkosHeatConductionTimeDerivative : public KokkosTimeDerivative
{
public:
  static InputParameters validParams();
 
  KokkosHeatConductionTimeDerivative(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;
 
protected:
  Moose::Kokkos::MaterialProperty<Real> _specific_heat;
  Moose::Kokkos::MaterialProperty<Real> _d_specific_heat_dT;
  Moose::Kokkos::MaterialProperty<Real> _density;
  Moose::Kokkos::MaterialProperty<Real> _d_density_dT;
};
 
template <typename Derived>
KOKKOS_FUNCTION Real
KokkosHeatConductionTimeDerivative::precomputeQpResidual(const unsigned int qp,
                                                         AssemblyDatum & datum) const
{
  return _specific_heat(datum, qp) * _density(datum, qp) *
         KokkosTimeDerivative::precomputeQpResidual<Derived>(qp, datum);
}
 
template <typename Derived>
KOKKOS_FUNCTION Real
KokkosHeatConductionTimeDerivative::precomputeQpJacobian(const unsigned int j,
                                                         const unsigned int qp,
                                                         AssemblyDatum & datum) const
{
  const Real cp = _specific_heat(datum, qp);
  const Real rho = _density(datum, qp);
 
  auto jac = cp * rho * KokkosTimeDerivative::precomputeQpJacobian<Derived>(j, qp, datum);
 
  if (_d_specific_heat_dT)
    jac += _d_specific_heat_dT(datum, qp) * rho * _phi(datum, j, qp) *
           KokkosTimeDerivative::precomputeQpResidual<Derived>(qp, datum);
  if (_d_density_dT)
    jac += cp * _d_density_dT(datum, qp) * _phi(datum, j, qp) *
           KokkosTimeDerivative::precomputeQpResidual<Derived>(qp, datum);
 
  return jac;
}