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 "SCMPinPowerPostprocessor.h"
      11             : #include "FEProblemBase.h"
      12             : #include "MooseMesh.h"
      13             : #include "SCM.h"
      14             : #include "SubChannel1PhaseProblem.h"
      15             : 
      16             : registerMooseObject("SubChannelApp", SCMPinPowerPostprocessor);
      17             : 
      18             : InputParameters
      19         144 : SCMPinPowerPostprocessor::validParams()
      20             : {
      21         144 :   InputParameters params = GeneralPostprocessor::validParams();
      22         144 :   params.addClassDescription(
      23             :       "Calculates the total heat rate that goes into the coolant from all the heated fuel-pins "
      24             :       "based on aux variable q_prime");
      25         144 :   return params;
      26           0 : }
      27             : 
      28          72 : SCMPinPowerPostprocessor::SCMPinPowerPostprocessor(const InputParameters & parameters)
      29             :   : GeneralPostprocessor(parameters),
      30          72 :     _mesh(SCM::getConstMesh<SubChannelMesh>(_fe_problem.mesh())),
      31          72 :     _value(0)
      32             : {
      33          72 : }
      34             : 
      35             : void
      36          72 : SCMPinPowerPostprocessor::execute()
      37             : {
      38          72 :   const auto scm_problem = dynamic_cast<SubChannel1PhaseProblem *>(&_fe_problem);
      39          72 :   const auto nz = _mesh.getNumOfAxialCells();
      40          72 :   const auto n_channels = _mesh.getNumOfChannels();
      41          72 :   if (!scm_problem)
      42           0 :     mooseError("SCMPinPowerPostprocessor can only be used within a subchannel 1phase problem.");
      43             :   auto power = 0.0;
      44        1236 :   for (unsigned int iz = 1; iz < nz + 1; iz++)
      45       39756 :     for (unsigned int i_ch = 0; i_ch < n_channels; i_ch++)
      46       38592 :       power += scm_problem->computeAddedHeatPin(i_ch, iz);
      47          72 :   _value = power;
      48          72 : }
      49             : 
      50             : Real
      51          72 : SCMPinPowerPostprocessor::getValue() const
      52             : {
      53          72 :   return _value;
      54             : }