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 "PorousFlowDesorpedMassTimeDerivative.h"
      11             : 
      12             : #include "MooseVariable.h"
      13             : 
      14             : #include "libmesh/quadrature.h"
      15             : 
      16             : registerMooseObject("PorousFlowApp", PorousFlowDesorpedMassTimeDerivative);
      17             : 
      18             : InputParameters
      19         183 : PorousFlowDesorpedMassTimeDerivative::validParams()
      20             : {
      21         183 :   InputParameters params = TimeKernel::validParams();
      22         366 :   params.addRequiredParam<UserObjectName>(
      23             :       "PorousFlowDictator", "The UserObject that holds the list of PorousFlow variable names.");
      24         366 :   params.addRequiredCoupledVar(
      25             :       "conc_var", "The variable that represents the concentration of desorped species");
      26         183 :   params.addClassDescription("Desorped component mass derivative wrt time.");
      27         183 :   return params;
      28           0 : }
      29             : 
      30          99 : PorousFlowDesorpedMassTimeDerivative::PorousFlowDesorpedMassTimeDerivative(
      31          99 :     const InputParameters & parameters)
      32             :   : TimeKernel(parameters),
      33          99 :     _dictator(getUserObject<PorousFlowDictator>("PorousFlowDictator")),
      34          99 :     _conc_var_number(coupled("conc_var")),
      35          99 :     _conc(coupledValue("conc_var")),
      36          99 :     _conc_old(coupledValueOld("conc_var")),
      37         198 :     _porosity(getMaterialProperty<Real>("PorousFlow_porosity_qp")),
      38         198 :     _porosity_old(getMaterialPropertyOld<Real>("PorousFlow_porosity_qp")),
      39         198 :     _dporosity_dvar(getMaterialProperty<std::vector<Real>>("dPorousFlow_porosity_qp_dvar")),
      40          99 :     _dporosity_dgradvar(
      41         198 :         getMaterialProperty<std::vector<RealGradient>>("dPorousFlow_porosity_qp_dgradvar"))
      42             : {
      43          99 : }
      44             : 
      45             : Real
      46        8680 : PorousFlowDesorpedMassTimeDerivative::computeQpResidual()
      47             : {
      48        8680 :   Real c = (1.0 - _porosity[_qp]) * _conc[_qp];
      49        8680 :   Real c_old = (1.0 - _porosity_old[_qp]) * _conc_old[_qp];
      50        8680 :   return _test[_i][_qp] * (c - c_old) / _dt;
      51             : }
      52             : 
      53             : Real
      54       41640 : PorousFlowDesorpedMassTimeDerivative::computeQpJacobian()
      55             : {
      56       41640 :   return computeQpJac(_var.number());
      57             : }
      58             : 
      59             : Real
      60       11520 : PorousFlowDesorpedMassTimeDerivative::computeQpOffDiagJacobian(unsigned int jvar)
      61             : {
      62       11520 :   return computeQpJac(jvar);
      63             : }
      64             : 
      65             : Real
      66       53160 : PorousFlowDesorpedMassTimeDerivative::computeQpJac(unsigned int jvar) const
      67             : {
      68             :   Real deriv = 0.0;
      69             : 
      70       53160 :   if (jvar == _conc_var_number)
      71        5800 :     deriv = (1.0 - _porosity[_qp]) * _phi[_j][_qp];
      72             : 
      73       53160 :   if (_dictator.notPorousFlowVariable(jvar))
      74           0 :     return _test[_i][_qp] * deriv / _dt;
      75       53160 :   const unsigned int pvar = _dictator.porousFlowVariableNum(jvar);
      76             : 
      77       53160 :   deriv -= _dporosity_dgradvar[_qp][pvar] * _grad_phi[_j][_qp] * _conc[_qp];
      78       53160 :   deriv -= _dporosity_dvar[_qp][pvar] * _phi[_j][_qp] * _conc[_qp];
      79             : 
      80       53160 :   return _test[_i][_qp] * deriv / _dt;
      81             : }