PorousFlowHeatMassTransfer.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 "PorousFlowHeatMassTransfer.h"

#include "MooseVariable.h"

registerMooseObject("PorousFlowApp", PorousFlowHeatMassTransfer);

InputParameters
PorousFlowHeatMassTransfer::validParams()
{
  InputParameters params = Kernel::validParams();

  params.addClassDescription(
      "Calculate heat or mass transfer from a coupled variable v to the variable u. "
      "No mass lumping is performed here.");
  params.addRequiredCoupledVar(
      "v", "The variable which is tranfsered to u using a transfer coefficient");
  params.addCoupledVar("transfer_coefficient",
                       1.0,
                       "Transfer coefficient for heat or mass transferred between variables");

  return params;
}

PorousFlowHeatMassTransfer::PorousFlowHeatMassTransfer(const InputParameters & parameters)
  : Kernel(parameters),
    _v_var(coupled("v")),
    _v(coupledValue("v")),
    _coef_var(coupledValue("transfer_coefficient"))
{
}

Real
PorousFlowHeatMassTransfer::computeQpResidual()
{
  return _coef_var[_qp] * (_u[_qp] - _v[_qp]) * _test[_i][_qp];
}

Real
PorousFlowHeatMassTransfer::computeQpJacobian()
{
  return jac(_var.number());
}

Real
PorousFlowHeatMassTransfer::computeQpOffDiagJacobian(unsigned int jvar)
{
  return jac(jvar);
}

Real
PorousFlowHeatMassTransfer::jac(unsigned int jvar) const
{
  if (jvar == _var.number())
    return _coef_var[_qp] * _phi[_j][_qp] * _test[_i][_qp];
  else if (jvar == _v_var)
    return -_coef_var[_qp] * _phi[_j][_qp] * _test[_i][_qp];
  return 0.0;
}