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 "ADHeatStructureHeatSource.h"
      11             : 
      12             : registerMooseObject("ThermalHydraulicsApp", ADHeatStructureHeatSource);
      13             : 
      14             : InputParameters
      15         184 : ADHeatStructureHeatSource::validParams()
      16             : {
      17         184 :   InputParameters params = ADKernel::validParams();
      18         184 :   params.addClassDescription("Adds a heat source term for the energy equation");
      19         368 :   params.addRequiredParam<Real>("power_fraction", "The fraction of power used");
      20         368 :   params.addRequiredCoupledVar("total_power", "Total reactor power");
      21         368 :   params.addRequiredParam<Real>("num_units", "The number of units");
      22         368 :   params.addRequiredParam<FunctionName>("power_shape_function",
      23             :                                         "The name of the function that defines the power shape");
      24         368 :   params.addRequiredParam<PostprocessorName>("power_shape_integral_pp",
      25             :                                              "Power shape integral post-processor name");
      26         368 :   params.addParam<Real>("scale", 1.0, "Scaling factor for residual");
      27         368 :   params.declareControllable("power_fraction");
      28         184 :   return params;
      29           0 : }
      30             : 
      31          98 : ADHeatStructureHeatSource::ADHeatStructureHeatSource(const InputParameters & parameters)
      32             :   : ADKernel(parameters),
      33          98 :     _power_fraction(getParam<Real>("power_fraction")),
      34          98 :     _total_power(coupledScalarValue("total_power")),
      35          98 :     _power_shape_function(getFunction("power_shape_function")),
      36          98 :     _power_shape_integral(getPostprocessorValue("power_shape_integral_pp")),
      37         196 :     _scale(getParam<Real>("scale")),
      38         294 :     _num_units(getParam<Real>("num_units"))
      39             : {
      40          98 : }
      41             : 
      42             : ADReal
      43     2428514 : ADHeatStructureHeatSource::computeQpResidual()
      44             : {
      45     2428514 :   const ADReal power = _power_fraction * _total_power[0];
      46             :   const ADReal power_density =
      47     2428514 :       power / (_num_units * _power_shape_integral) * _power_shape_function.value(_t, _q_point[_qp]);
      48     4857028 :   return -_scale * power_density * _test[_i][_qp];
      49             : }