Line data    Source code

       1             : /****************************************************************/
       2             : /*               DO NOT MODIFY THIS HEADER                      */
       3             : /*                       BlackBear                              */
       4             : /*                                                              */
       5             : /*           (c) 2017 Battelle Energy Alliance, LLC             */
       6             : /*                   ALL RIGHTS RESERVED                        */
       7             : /*                                                              */
       8             : /*          Prepared by Battelle Energy Alliance, LLC           */
       9             : /*            Under Contract No. DE-AC07-05ID14517              */
      10             : /*            With the U. S. Department of Energy               */
      11             : /*                                                              */
      12             : /*            See COPYRIGHT for full restrictions               */
      13             : /****************************************************************/
      14             : 
      15             : #include "ConcreteThermalConvection.h"
      16             : 
      17             : registerMooseObject("BlackBearApp", ConcreteThermalConvection);
      18             : 
      19             : InputParameters
      20         418 : ConcreteThermalConvection::validParams()
      21             : {
      22         418 :   InputParameters params = Kernel::validParams();
      23         418 :   params.addClassDescription(
      24             :       "Convective transport term for heat transfer due to fluid flow in concrete.");
      25         836 :   params.addRequiredCoupledVar("relative_humidity", "nonlinear variable name for rel. humidity");
      26         418 :   return params;
      27           0 : }
      28             : 
      29         220 : ConcreteThermalConvection::ConcreteThermalConvection(const InputParameters & parameters)
      30             :   : Kernel(parameters),
      31         220 :     _cw(getMaterialProperty<Real>("thermal_capacity_of_water")),
      32         440 :     _Dh(getMaterialProperty<Real>("humidity_diffusivity")),
      33         220 :     _grad_rh(coupledGradient("relative_humidity")),
      34         440 :     _rh_var(coupled("relative_humidity"))
      35             : {
      36         220 : }
      37             : 
      38             : Real
      39     2561328 : ConcreteThermalConvection::computeQpResidual()
      40             : {
      41     2561328 :   return _cw[_qp] * _Dh[_qp] * _grad_rh[_qp] * _grad_u[_qp] * _test[_i][_qp];
      42             : }
      43             : 
      44             : Real
      45      815584 : ConcreteThermalConvection::computeQpJacobian()
      46             : {
      47      815584 :   return _cw[_qp] * _Dh[_qp] * _grad_rh[_qp] * _test[_i][_qp] * _grad_phi[_j][_qp];
      48             : }
      49             : 
      50             : Real
      51      807520 : ConcreteThermalConvection::computeQpOffDiagJacobian(unsigned int jvar)
      52             : {
      53      807520 :   if (jvar == _rh_var)
      54      807520 :     return _cw[_qp] * _Dh[_qp] * _grad_u[_qp] * _test[_i][_qp] * _grad_phi[_j][_qp];
      55             : 
      56             :   return 0.0;
      57             : }