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 "RadiativeP1DiffusionCoefficientMaterial.h"
      11             : 
      12             : registerMooseObject("HeatTransferApp", RadiativeP1DiffusionCoefficientMaterial);
      13             : registerMooseObject("HeatTransferApp", ADRadiativeP1DiffusionCoefficientMaterial);
      14             : 
      15             : template <bool is_ad>
      16             : InputParameters
      17          38 : RadiativeP1DiffusionCoefficientMaterialTempl<is_ad>::validParams()
      18             : {
      19          38 :   InputParameters params = FunctorMaterial::validParams();
      20          38 :   params.addClassDescription("Computes the P1 diffusion coefficient from the opacity and effective "
      21             :                              "scattering cross section.");
      22          76 :   params.addRequiredParam<MooseFunctorName>("opacity", "Opacity.");
      23          76 :   params.addParam<MooseFunctorName>("sigma_scat_eff", 0.0, "Effective P1 scatterig cross section.");
      24          76 :   params.addRequiredParam<std::string>("P1_diff_coef_name",
      25             :                                        "Name given to the P1 diffusion coefficient.");
      26          38 :   return params;
      27           0 : }
      28             : 
      29             : template <bool is_ad>
      30          20 : RadiativeP1DiffusionCoefficientMaterialTempl<is_ad>::RadiativeP1DiffusionCoefficientMaterialTempl(
      31             :     const InputParameters & parameters)
      32             :   : FunctorMaterial(parameters),
      33          20 :     _opacity(getFunctor<GenericReal<is_ad>>("opacity")),
      34          60 :     _sigma_scat_eff(getFunctor<GenericReal<is_ad>>("sigma_scat_eff"))
      35             : {
      36          80 :   addFunctorProperty<GenericReal<is_ad>>(
      37          40 :       getParam<std::string>("P1_diff_coef_name"),
      38          24 :       [this](const auto & r, const auto & t) -> GenericReal<is_ad>
      39             :       {
      40          24 :         const auto kappa = _opacity(r, t);
      41          24 :         const auto sigma_scat = _sigma_scat_eff(r, t);
      42          24 :         return 1.0 / (3.0 * kappa + sigma_scat);
      43             :       });
      44          40 : }
      45             : 
      46             : template class RadiativeP1DiffusionCoefficientMaterialTempl<false>;
      47             : template class RadiativeP1DiffusionCoefficientMaterialTempl<true>;