PorousFlowPorosityConst.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 "PorousFlowPorosityConst.h"
 
registerMooseObject("PorousFlowApp", PorousFlowPorosityConst);
registerMooseObject("PorousFlowApp", ADPorousFlowPorosityConst);
 
template <bool is_ad>
InputParameters
PorousFlowPorosityConstTempl<is_ad>::validParams()
{
  InputParameters params = PorousFlowPorosityBaseTempl<is_ad>::validParams();
  params.addRequiredCoupledVar(
      "porosity",
      "The porosity (assumed indepenent of porepressure, temperature, "
      "strain, etc, for this material).  This should be a real number, or "
      "a constant monomial variable (not a linear lagrange or other kind of variable).  If it is "
      "less than porosity_min, porosity is set to porosity_min instead.");
  // The porosity derivatives are always zero in this Material, so zero_modifier has no effect
  params.suppressParameter<Real>("zero_modifier");
  params.addClassDescription("This Material calculates the porosity assuming it is constant");
  return params;
}
 
template <bool is_ad>
PorousFlowPorosityConstTempl<is_ad>::PorousFlowPorosityConstTempl(
    const InputParameters & parameters)
  : PorousFlowPorosityBaseTempl<is_ad>(parameters), _input_porosity(coupledValue("porosity"))
{
}
 
template <bool is_ad>
void
PorousFlowPorosityConstTempl<is_ad>::initQpStatefulProperties()
{
  // note the [0] below: _phi0 is a constant monomial and we use [0] regardless of _nodal_material
  _porosity[_qp] = _input_porosity[0];
 
  if (!is_ad)
  {
    // The derivatives are zero for all time
    (*_dporosity_dvar)[_qp].assign(_num_var, 0.0);
    (*_dporosity_dgradvar)[_qp].assign(_num_var, RealGradient());
  }
 
  this->applyPorosityMin();
}
 
template <bool is_ad>
void
PorousFlowPorosityConstTempl<is_ad>::computeQpProperties()
{
  initQpStatefulProperties();
}
 
template class PorousFlowPorosityConstTempl<false>;
template class PorousFlowPorosityConstTempl<true>;