MoltresSNMaterial.h

Go to the documentation of this file.

#pragma once
 
#include "GenericConstantMaterial.h"
#include "SplineInterpolation.h"
#include "MonotoneCubicInterpolation.h"
#include "LinearInterpolation.h"
#include "nlohmann/json.h"
 
class MoltresSNMaterial : public GenericConstantMaterial
{
public:
  MoltresSNMaterial(const InputParameters & parameters);
 
  static InputParameters validParams();
 
  // correspond to descriptions above
  enum INTERPOLATOR
  {
    SPLINE,
    MONOTONE_CUBIC,
    LINEAR,
    NONE
  };
 
  // returns a MooseEnum corresponding to the above enum
  static MooseEnum interpTypes()
  {
    return MooseEnum("spline=0 monotone_cubic=1 linear=2 none=3");
  }
 
protected:
  virtual void dummyComputeQpProperties();
  virtual void splineComputeQpProperties();
  virtual void monotoneCubicComputeQpProperties();
  virtual void linearComputeQpProperties();
  virtual void preComputeQpProperties();
  void Construct(nlohmann::json xs_root);
  virtual void computeQpProperties() override;
 
  const VariableValue & _temperature;
 
  unsigned int _num_groups;
 
  unsigned int _num_precursor_groups;
 
  MaterialProperty<std::vector<Real>> & _totxs;
  MaterialProperty<std::vector<Real>> & _fissxs;
  MaterialProperty<std::vector<Real>> & _nsf;
  MaterialProperty<std::vector<Real>> & _fisse;
  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>> & _scatter;
  MaterialProperty<std::vector<Real>> & _beta_eff;
  MaterialProperty<Real> & _beta;
  MaterialProperty<std::vector<Real>> & _decay_constant;
  MaterialProperty<std::vector<Real>> & _diffcoef;
  MaterialProperty<std::vector<Real>> & _d_totxs_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_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_scatter_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;
  MaterialProperty<std::vector<Real>> & _d_diffcoef_d_temp;
 
  MooseEnum _interp_type;
 
  std::vector<std::string> _group_consts;
 
  std::string _material_key;
 
  const unsigned int _L;
 
  const Real _sigma;
 
  const Real _c;
 
  const MooseEnum _h_type;
 
  MaterialProperty<std::vector<Real>> & _tau_sn;
 
  std::vector<std::string> _xsec_names{"TOTXS",
                                       "FISSXS",
                                       "NSF",
                                       "FISSE",
                                       "RECIPVEL",
                                       "CHI_T",
                                       "CHI_P",
                                       "CHI_D",
                                       "SPN",
                                       "BETA_EFF",
                                       "DECAY_CONSTANT",
                                       "DIFFCOEF"};
 
  std::map<std::string, std::vector<std::vector<Real>>> _xsec_map;
 
  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, int> _vec_lengths;
 
  std::vector<double> _XsTemperature;
};