DesorptionFromMatrix.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 "DesorptionFromMatrix.h"
 
#include <iostream>
 
registerMooseObject("ChemicalReactionsApp", DesorptionFromMatrix);
 
template <>
InputParameters
validParams<DesorptionFromMatrix>()
{
  InputParameters params = validParams<Kernel>();
  params.addRequiredCoupledVar(
      "pressure_var",
      "Variable representing the porepressure of the fluid adsorbed into the matrix");
  params.addClassDescription("Mass flow rate from the matrix to the porespace.  Add this to "
                             "TimeDerivative kernel to get complete DE for the fluid adsorbed "
                             "in the matrix");
  return params;
}
 
DesorptionFromMatrix::DesorptionFromMatrix(const InputParameters & parameters)
  : Kernel(parameters),
    _pressure_var(coupled("pressure_var")),
    _mass_rate_from_matrix(getMaterialProperty<Real>("mass_rate_from_matrix")),
    _dmass_rate_from_matrix_dC(getMaterialProperty<Real>("dmass_rate_from_matrix_dC")),
    _dmass_rate_from_matrix_dp(getMaterialProperty<Real>("dmass_rate_from_matrix_dp"))
{
}
 
Real
DesorptionFromMatrix::computeQpResidual()
{
  return _test[_i][_qp] * _mass_rate_from_matrix[_qp];
}
 
Real
DesorptionFromMatrix::computeQpJacobian()
{
  return _test[_i][_qp] * _dmass_rate_from_matrix_dC[_qp] * _phi[_j][_qp];
}
 
Real
DesorptionFromMatrix::computeQpOffDiagJacobian(unsigned int jvar)
{
  if (jvar != _pressure_var)
    return 0.0;
  return _test[_i][_qp] * _dmass_rate_from_matrix_dp[_qp] * _phi[_j][_qp];
}