NuclearMaterial.h

Go to the documentation of this file.

#pragma once
 
#include "GenericConstantMaterial.h"
#include "SplineInterpolation.h"
#include "BicubicSplineInterpolation.h"
#include "MonotoneCubicInterpolation.h"
#include "LinearInterpolation.h"
#include "json.h"
 
class NuclearMaterial : public GenericConstantMaterial
{
public:
  NuclearMaterial(const InputParameters & parameters);
 
  static InputParameters validParams();
 
  // correspond to descriptions above
  enum INTERPOLATOR
  {
    BICUBIC,
    SPLINE,
    MONOTONE_CUBIC,
    LINEAR,
    NONE,
    LSQ
  };
 
  // returns a MooseEnum corresponding to the above enum
  static MooseEnum interpTypes()
  {
    return MooseEnum("bicubic=0 spline=1"
                     " monotone_cubic=2 linear=3 none=4 least_squares=5");
  }
 
protected:
  virtual void dummyComputeQpProperties();
  virtual void splineComputeQpProperties();
  virtual void monotoneCubicComputeQpProperties();
  virtual void linearComputeQpProperties();
  virtual void preComputeQpProperties();
 
  const VariableValue & _temperature;
 
  // Number of neutron groups
  unsigned int _num_groups;
 
  // Number of precursor groups
  unsigned int _num_precursor_groups;
 
  // Group constant MaterialProperty(s)
  MaterialProperty<std::vector<Real>> & _remxs;
  MaterialProperty<std::vector<Real>> & _fissxs;
  MaterialProperty<std::vector<Real>> & _nsf;
  MaterialProperty<std::vector<Real>> & _fisse;
  MaterialProperty<std::vector<Real>> & _diffcoef;
  MaterialProperty<std::vector<Real>> & _recipvel;
  MaterialProperty<std::vector<Real>> & _chi_t;
  MaterialProperty<std::vector<Real>> & _chi_p;
  MaterialProperty<std::vector<Real>> & _chi_d;
  MaterialProperty<std::vector<Real>> & _gtransfxs;
  MaterialProperty<std::vector<Real>> & _beta_eff;
  MaterialProperty<Real> & _beta;
  MaterialProperty<std::vector<Real>> & _decay_constant;
  MaterialProperty<std::vector<Real>> & _d_remxs_d_temp;
  MaterialProperty<std::vector<Real>> & _d_fissxs_d_temp;
  MaterialProperty<std::vector<Real>> & _d_nsf_d_temp;
  MaterialProperty<std::vector<Real>> & _d_fisse_d_temp;
  MaterialProperty<std::vector<Real>> & _d_diffcoef_d_temp;
  MaterialProperty<std::vector<Real>> & _d_recipvel_d_temp;
  MaterialProperty<std::vector<Real>> & _d_chi_t_d_temp;
  MaterialProperty<std::vector<Real>> & _d_chi_p_d_temp;
  MaterialProperty<std::vector<Real>> & _d_chi_d_d_temp;
  MaterialProperty<std::vector<Real>> & _d_gtransfxs_d_temp;
  MaterialProperty<std::vector<Real>> & _d_beta_eff_d_temp;
  MaterialProperty<Real> & _d_beta_d_temp;
  MaterialProperty<std::vector<Real>> & _d_decay_constant_d_temp;
 
  // Group constant interpolation type
  MooseEnum _interp_type;
 
  // Vector of group constant names
  std::vector<std::string> _xsec_names{"REMXS",
                                       "FISSXS",
                                       "NSF",
                                       "FISSE",
                                       "DIFFCOEF",
                                       "RECIPVEL",
                                       "CHI_T",
                                       "CHI_P",
                                       "CHI_D",
                                       "GTRANSFXS",
                                       "BETA_EFF",
                                       "DECAY_CONSTANT"};
 
  // Map of group constant names to group constant values
  std::map<std::string, std::vector<std::vector<Real>>> _xsec_map;
 
  // Group constant interpolators
  std::map<std::string, std::vector<SplineInterpolation>> _xsec_spline_interpolators;
  std::map<std::string, std::vector<MonotoneCubicInterpolation>> _xsec_monotone_cubic_interpolators;
  std::map<std::string, std::vector<LinearInterpolation>> _xsec_linear_interpolators;
  std::map<std::string, std::vector<BicubicSplineInterpolation>> _xsec_bicubic_spline_interpolators;
 
  // Map of group constant names to number of neutron/precursor groups
  std::map<std::string, int> _vec_lengths;
 
  // Vector of temperature values
  std::vector<double> _XsTemperature;
 
  std::vector<std::vector<Real>> _remxs_consts = std::vector<std::vector<Real>>(2);
  std::vector<std::vector<Real>> _fissxs_consts = std::vector<std::vector<Real>>(2);
  std::vector<std::vector<Real>> _nubar_consts = std::vector<std::vector<Real>>(2);
  std::vector<std::vector<Real>> _nsf_consts = std::vector<std::vector<Real>>(2);
  std::vector<std::vector<Real>> _fisse_consts = std::vector<std::vector<Real>>(2);
  std::vector<std::vector<Real>> _diffcoeff_consts = std::vector<std::vector<Real>>(2);
  std::vector<std::vector<Real>> _recipvel_consts = std::vector<std::vector<Real>>(2);
  std::vector<std::vector<Real>> _chi_t_consts = std::vector<std::vector<Real>>(2);
  std::vector<std::vector<Real>> _chi_p_consts = std::vector<std::vector<Real>>(2);
  std::vector<std::vector<Real>> _chi_d_consts = std::vector<std::vector<Real>>(2);
  std::vector<std::vector<Real>> _gtransfxs_consts = std::vector<std::vector<Real>>(2);
  std::vector<std::vector<Real>> _beta_eff_consts = std::vector<std::vector<Real>>(2);
  std::vector<std::vector<Real>> _decay_constants_consts = std::vector<std::vector<Real>>(2);
};