LinearFVReaction.C

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

#include "LinearFVReaction.h"
#include "Assembly.h"
#include "SubProblem.h"

registerMooseObject("MooseApp", LinearFVReaction);

InputParameters
LinearFVReaction::validParams()
{
  InputParameters params = LinearFVElementalKernel::validParams();
  params.addClassDescription(
      "Represents the matrix and right hand side contributions of a reaction "
      "term ($c u$) in a partial differential equation.");
  params.addParam<MooseFunctorName>("coeff", 1.0, "The reaction coefficient.");
  return params;
}

LinearFVReaction::LinearFVReaction(const InputParameters & params)
  : LinearFVElementalKernel(params), _coefficient(getFunctor<Real>("coeff"))
{
}

Real
LinearFVReaction::computeMatrixContribution()
{
  // The matrix contribution is c_C*V_C
  return _coefficient(makeElemArg(_current_elem_info->elem()), determineState()) *
         _current_elem_volume;
}

Real
LinearFVReaction::computeRightHandSideContribution()
{
  // We don't have any contributions to the right hand side as the reaction term is
  // treated implicitly. If we treated it explicitly, the contribution to the RHS would be
  // c * u * elem_volume
  return 0.0;
}