OpenMCVolumeCalculation.h

Go to the documentation of this file.

/********************************************************************/
/*                  SOFTWARE COPYRIGHT NOTIFICATION                 */
/*                             Cardinal                             */
/*                                                                  */
/*                  (c) 2021 UChicago Argonne, LLC                  */
/*                        ALL RIGHTS RESERVED                       */
/*                                                                  */
/*                 Prepared by UChicago Argonne, LLC                */
/*               Under Contract No. DE-AC02-06CH11357               */
/*                With the U. S. Department of Energy               */
/*                                                                  */
/*             Prepared by Battelle Energy Alliance, LLC            */
/*               Under Contract No. DE-AC07-05ID14517               */
/*                With the U. S. Department of Energy               */
/*                                                                  */
/*                 See LICENSE for full restrictions                */
/********************************************************************/

#pragma once

#include "GeneralUserObject.h"

#include "OpenMCBase.h"
#include "openmc/volume_calc.h"

class OpenMCVolumeCalculation : public GeneralUserObject, public OpenMCBase
{
public:
  static InputParameters validParams();

  OpenMCVolumeCalculation(const InputParameters & parameters);

  virtual void initialize() {}
  virtual void finalize() {}
  virtual void execute() {}

  virtual void initializeVolumeCalculation();

  virtual void computeVolumes();

  virtual void resetVolumeCalculation();

  openmc::Position position(const Point & pt) const;

  void cellVolume(const unsigned int & index, Real & vol, Real & std_dev) const;

protected:
  const unsigned int & _n_samples;

  const MooseEnum _trigger;

  Real _trigger_threshold;

  Real _scaling;

  Point _lower_left;

  Point _upper_right;

  std::unique_ptr<openmc::VolumeCalculation> _volume_calc;

  std::vector<openmc::VolumeCalculation::Result> _results;

  std::map<int, int> _index_to_calc_index;

  unsigned int _calc_index;
};