Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "PorousFlowFullySaturatedHeatAdvection.h"
      11             : 
      12             : registerMooseObject("PorousFlowApp", PorousFlowFullySaturatedHeatAdvection);
      13             : 
      14             : InputParameters
      15         452 : PorousFlowFullySaturatedHeatAdvection::validParams()
      16             : {
      17         452 :   InputParameters params = PorousFlowFullySaturatedDarcyBase::validParams();
      18         452 :   params.addClassDescription("Heat flux that arises from the advection of a fully-saturated single "
      19             :                              "phase fluid.  No upwinding is used");
      20         452 :   return params;
      21           0 : }
      22             : 
      23         243 : PorousFlowFullySaturatedHeatAdvection::PorousFlowFullySaturatedHeatAdvection(
      24         243 :     const InputParameters & parameters)
      25             :   : PorousFlowFullySaturatedDarcyBase(parameters),
      26         243 :     _enthalpy(getMaterialProperty<std::vector<Real>>("PorousFlow_fluid_phase_enthalpy_qp")),
      27         486 :     _denthalpy_dvar(getMaterialProperty<std::vector<std::vector<Real>>>(
      28         243 :         "dPorousFlow_fluid_phase_enthalpy_qp_dvar"))
      29             : {
      30         243 : }
      31             : 
      32             : Real
      33    89754280 : PorousFlowFullySaturatedHeatAdvection::mobility() const
      34             : {
      35             :   const unsigned ph = 0;
      36    89754280 :   return _enthalpy[_qp][ph] * PorousFlowFullySaturatedDarcyBase::mobility();
      37             : }
      38             : 
      39             : Real
      40    85059008 : PorousFlowFullySaturatedHeatAdvection::dmobility(unsigned pvar) const
      41             : {
      42             :   const unsigned ph = 0;
      43    85059008 :   const Real darcy_mob = PorousFlowFullySaturatedDarcyBase::mobility();
      44    85059008 :   const Real ddarcy_mob = PorousFlowFullySaturatedDarcyBase::dmobility(pvar);
      45    85059008 :   return _denthalpy_dvar[_qp][ph][pvar] * darcy_mob + _enthalpy[_qp][ph] * ddarcy_mob;
      46             : }