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