PrimaryConvection.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 "PrimaryConvection.h"
 
registerMooseObject("ChemicalReactionsApp", PrimaryConvection);
 
InputParameters
PrimaryConvection::validParams()
{
  InputParameters params = Kernel::validParams();
  params.addRequiredCoupledVar("p", "Pressure");
  RealVectorValue g(0, 0, 0);
  params.addParam<RealVectorValue>("gravity", g, "Gravity vector (default is (0, 0, 0))");
  params.addClassDescription("Convection of primary species");
  return params;
}
 
PrimaryConvection::PrimaryConvection(const InputParameters & parameters)
  : DerivativeMaterialInterface<Kernel>(parameters),
    _cond(getMaterialProperty<Real>("conductivity")),
    _gravity(getParam<RealVectorValue>("gravity")),
    _density(getDefaultMaterialProperty<Real>("density")),
    _grad_p(coupledGradient("p")),
    _pvar(coupled("p"))
 
{
}
 
Real
PrimaryConvection::computeQpResidual()
{
  RealVectorValue darcy_vel = -_cond[_qp] * (_grad_p[_qp] - _density[_qp] * _gravity);
 
  return _test[_i][_qp] * (darcy_vel * _grad_u[_qp]);
}
 
Real
PrimaryConvection::computeQpJacobian()
{
  RealVectorValue darcy_vel = -_cond[_qp] * (_grad_p[_qp] - _density[_qp] * _gravity);
 
  return _test[_i][_qp] * (darcy_vel * _grad_phi[_j][_qp]);
}
 
Real
PrimaryConvection::computeQpOffDiagJacobian(unsigned int jvar)
{
  if (jvar == _pvar)
  {
    RealVectorValue ddarcy_vel_dp = -_cond[_qp] * _grad_phi[_j][_qp];
    return _test[_i][_qp] * (ddarcy_vel_dp * _grad_u[_qp]);
  }
  else
    return 0.0;
}