CoupledBEKinetic.C

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
 
#include "CoupledBEKinetic.h"
 
registerMooseObject("ChemicalReactionsApp", CoupledBEKinetic);
 
template <>
InputParameters
validParams<CoupledBEKinetic>()
{
  InputParameters params = validParams<TimeDerivative>();
  params.addRequiredParam<std::vector<Real>>("weight",
                                             "The weight of kinetic species concentration");
  params.addCoupledVar("v", "List of kinetic species being coupled by concentration");
  params.addClassDescription("Derivative of kinetic species concentration wrt time");
  return params;
}
 
CoupledBEKinetic::CoupledBEKinetic(const InputParameters & parameters)
  : TimeDerivative(parameters),
    _porosity(getMaterialProperty<Real>("porosity")),
    _weight(getParam<std::vector<Real>>("weight"))
{
  const unsigned int n = coupledComponents("v");
  _vals.resize(n);
  _vals_old.resize(n);
 
  for (unsigned int i = 0; i < n; ++i)
  {
    _vals[i] = &coupledValue("v", i);
    _vals_old[i] = &coupledValueOld("v", i);
  }
}
 
Real
CoupledBEKinetic::computeQpResidual()
{
  Real assemble_conc = 0.0;
 
  for (MooseIndex(_vals) i = 0; i < _vals.size(); ++i)
    assemble_conc += _weight[i] * ((*_vals[i])[_qp] - (*_vals_old[i])[_qp]) / _dt;
 
  return _porosity[_qp] * _test[_i][_qp] * assemble_conc;
}