PorousFlowSinkPTDefiner.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 "PorousFlowSinkPTDefiner.h"
 
template <>
InputParameters
validParams<PorousFlowSinkPTDefiner>()
{
  InputParameters params = validParams<PorousFlowSink>();
  params.addCoupledVar("PT_shift",
                       0.0,
                       "Whenever the sink is an explicit function of porepressure "
                       "(such as a PiecewiseLinear function) the argument of the "
                       "function is set to P - PT_shift instead of simply P.  "
                       "Similarly for temperature.  PT_shift does not enter into "
                       "any use_* calculations.");
  return params;
}
 
PorousFlowSinkPTDefiner::PorousFlowSinkPTDefiner(const InputParameters & parameters)
  : PorousFlowSink(parameters),
    _pp(_involves_fluid ? &getMaterialProperty<std::vector<Real>>("PorousFlow_porepressure_nodal")
                        : nullptr),
    _dpp_dvar(_involves_fluid ? &getMaterialProperty<std::vector<std::vector<Real>>>(
                                    "dPorousFlow_porepressure_nodal_dvar")
                              : nullptr),
    _temp(!_involves_fluid ? &getMaterialProperty<Real>("PorousFlow_temperature_nodal") : nullptr),
    _dtemp_dvar(!_involves_fluid
                    ? &getMaterialProperty<std::vector<Real>>("dPorousFlow_temperature_nodal_dvar")
                    : nullptr),
    _pt_shift(coupledDofValues("PT_shift"))
{
  if (_involves_fluid && (_pp == nullptr || _dpp_dvar == nullptr))
    mooseError("PorousFlowSink: There is no porepressure Material");
  if (!_involves_fluid && (_temp == nullptr || _dtemp_dvar == nullptr))
    mooseError("PorousFlowSink: There is no temperature Material");
}
 
Real
PorousFlowSinkPTDefiner::ptVar() const
{
  if (_involves_fluid)
    return (*_pp)[_i][_ph] - _pt_shift[_i];
  return (*_temp)[_i] - _pt_shift[_i];
}
 
Real
PorousFlowSinkPTDefiner::dptVar(unsigned pvar) const
{
  if (_involves_fluid)
    return (*_dpp_dvar)[_i][_ph][pvar];
  return (*_dtemp_dvar)[_i][pvar];
}