PorousFlowPermeabilityConst.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 "PorousFlowPermeabilityConst.h"
 
registerMooseObject("PorousFlowApp", PorousFlowPermeabilityConst);
 
template <>
InputParameters
validParams<PorousFlowPermeabilityConst>()
{
  InputParameters params = validParams<PorousFlowPermeabilityBase>();
  params.addRequiredParam<RealTensorValue>(
      "permeability",
      "The permeability tensor (usually in m^2), which is assumed constant for this material");
  params.addClassDescription(
      "This Material calculates the permeability tensor assuming it is constant");
  return params;
}
 
PorousFlowPermeabilityConst::PorousFlowPermeabilityConst(const InputParameters & parameters)
  : PorousFlowPermeabilityBase(parameters),
    _input_permeability(getParam<RealTensorValue>("permeability"))
{
}
 
void
PorousFlowPermeabilityConst::computeQpProperties()
{
  _permeability_qp[_qp] = _input_permeability;
  _dpermeability_qp_dvar[_qp].assign(_num_var, RealTensorValue());
  _dpermeability_qp_dgradvar[_qp].resize(LIBMESH_DIM);
  for (unsigned i = 0; i < LIBMESH_DIM; ++i)
    _dpermeability_qp_dgradvar[_qp][i].assign(_num_var, RealTensorValue());
}