RayTracingViewFactor.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
#include "RayTracingViewFactor.h"
 
#include "ViewFactorRayStudy.h"
 
registerMooseObject("HeatTransferApp", RayTracingViewFactor);
 
InputParameters
RayTracingViewFactor::validParams()
{
  InputParameters params = ViewFactorBase::validParams();
  params.addRequiredParam<UserObjectName>("ray_study_name",
                                          "Name of the view factor ray study UO.");
  params.addClassDescription("Computes view factors for arbitrary geometries using raytracing.");
  return params;
}
 
RayTracingViewFactor::RayTracingViewFactor(const InputParameters & parameters)
  : ViewFactorBase(parameters), _ray_study(getUserObject<ViewFactorRayStudy>("ray_study_name"))
{
  if (_mesh.dimension() == 1)
    mooseError("View factor calculations do not support 1D");
}
 
void
RayTracingViewFactor::threadJoinViewFactor(const UserObject & y)
{
  const auto & vf = static_cast<const RayTracingViewFactor &>(y);
  for (unsigned int i = 0; i < _n_sides; ++i)
    for (unsigned int j = 0; j < _n_sides; ++j)
      _view_factors[i][j] += vf._view_factors[i][j];
}
 
void
RayTracingViewFactor::finalizeViewFactor()
{
  for (unsigned int i = 0; i < _n_sides; ++i)
    gatherSum(_view_factors[i]);
 
  // get the _view_factors from ray study
  for (const auto & from_name : boundaryNames())
    for (const auto & to_name : boundaryNames())
    {
      unsigned int from_index = getSideNameIndex(from_name);
      unsigned int to_index = getSideNameIndex(to_name);
      BoundaryID from_id = _mesh.getBoundaryID(from_name);
      BoundaryID to_id = _mesh.getBoundaryID(to_name);
      _view_factors[from_index][to_index] = _ray_study.viewFactorInfo(from_id, to_id);
    }
 
  // divide view_factor Fij by Ai and pi
  const Real divisor = _mesh.dimension() == 2 ? 2 : libMesh::pi;
  for (unsigned int i = 0; i < _n_sides; ++i)
  {
    const Real factor = 1. / (_areas[i] * divisor);
    for (auto & vf : _view_factors[i])
      vf *= factor;
  }
}