UNVolumetricSwellingEigenstrain

Computes the change in volume due to swelling in UN fuel. This class applies a volumetric strain correction before adding the strain from this class to the diagonal entries of the eigenstrain tensor.

Description

This class, UNVolumetricSwellingEigenstrain, calculates the volumetric swelling of UN fuel. An empirical expression for the swelling of UN was determined using data from S. Ross (1989) as $var element = document.getElementById("moose-equation-ab74019c-55e9-498f-ae31-a5b81e625bec");katex.render("s = 0.001 \\times 4.7 \\times 10^{-11} \\times T^{3.12} \\times (100 \\times Bu)^{0.83} \\sqrt{TD}", element, {displayMode:true,throwOnError:false});$ where $var element = document.getElementById("moose-equation-a98359ff-af5a-4585-880d-35b13b9c8228");katex.render("Bu", element, {displayMode:false,throwOnError:false});$ is the instantaneous burnup, $var element = document.getElementById("moose-equation-8e3fd5e7-7887-4bbe-80e1-9ee9847e3ff5");katex.render("T", element, {displayMode:false,throwOnError:false});$ the temperature in Kelvin, and $var element = document.getElementById("moose-equation-63ee2888-a26d-4b33-a873-6b0463c6fd0f");katex.render("TD", element, {displayMode:false,throwOnError:false});$ the initial density as a percent of theoretical. The volumetric swelling is calculated by integrating the incremental strain over burnup.

Automatic differentiation is available using ADUNVolumetricSwellingEigenstrain.

Example Input Syntax

[Materials<<<{"href": "../../../syntax/Materials/index.html"}>>>]
  [swelling]
    type = UNVolumetricSwellingEigenstrain<<<{"description": "Computes the change in volume due to swelling in UN fuel. This class applies a volumetric strain correction before adding the strain from this class to the diagonal entries of the eigenstrain tensor.", "href": "UNVolumetricSwellingEigenstrain.html"}>>>
    temperature<<<{"description": "Coupled temperature in Kelvin"}>>> = temp
    burnup<<<{"description": "Coupled Burnup variable in FIMA"}>>> = burnup
    density_percent<<<{"description": "The initial theoretical density (percent)."}>>> = 95
    eigenstrain_name<<<{"description": "Material property name for the eigenstrain tensor computed by this model. IMPORTANT: The name of this property must also be provided to the strain calculator."}>>> = swelling
  []
[]

Input Parameters

burnupCoupled Burnup variable in FIMA
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:Coupled Burnup variable in FIMA
density_percentThe initial theoretical density (percent).
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:The initial theoretical density (percent).
eigenstrain_nameMaterial property name for the eigenstrain tensor computed by this model. IMPORTANT: The name of this property must also be provided to the strain calculator.
C++ Type:std::string
Controllable:No
Description:Material property name for the eigenstrain tensor computed by this model. IMPORTANT: The name of this property must also be provided to the strain calculator.
temperatureCoupled temperature in Kelvin
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:Coupled temperature in Kelvin

Required Parameters

base_nameOptional parameter that allows the user to define multiple mechanics material systems on the same block, i.e. for multiple phases
C++ Type:std::string
Controllable:No
Description:Optional parameter that allows the user to define multiple mechanics material systems on the same block, i.e. for multiple phases
blockThe list of blocks (ids or names) that this object will be applied
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of blocks (ids or names) that this object will be applied
boundaryThe list of boundaries (ids or names) from the mesh where this object applies
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:The list of boundaries (ids or names) from the mesh where this object applies
computeTrueWhen false, MOOSE will not call compute methods on this material. The user must call computeProperties() after retrieving the MaterialBase via MaterialBasePropertyInterface::getMaterialBase(). Non-computed MaterialBases are not sorted for dependencies.
Default:True
C++ Type:bool
Controllable:No
Description:When false, MOOSE will not call compute methods on this material. The user must call computeProperties() after retrieving the MaterialBase via MaterialBasePropertyInterface::getMaterialBase(). Non-computed MaterialBases are not sorted for dependencies.
constant_onNONEWhen ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped
Default:NONE
C++ Type:MooseEnum
Options:NONE, ELEMENT, SUBDOMAIN
Controllable:No
Description:When ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped
declare_suffixAn optional suffix parameter that can be appended to any declared properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:An optional suffix parameter that can be appended to any declared properties. The suffix will be prepended with a '_' character.

Optional Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Controllable:Yes
Description:Set the enabled status of the MooseObject.
implicitTrueDetermines whether this object is calculated using an implicit or explicit form
Default:True
C++ Type:bool
Controllable:No
Description:Determines whether this object is calculated using an implicit or explicit form
seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
Controllable:No
Description:The seed for the master random number generator
use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

Advanced Parameters

output_propertiesList of material properties, from this material, to output (outputs must also be defined to an output type)
C++ Type:std::vector<std::string>
Controllable:No
Description:List of material properties, from this material, to output (outputs must also be defined to an output type)
outputsnone Vector of output names where you would like to restrict the output of variables(s) associated with this object
Default:none
C++ Type:std::vector<OutputName>
Controllable:No
Description:Vector of output names where you would like to restrict the output of variables(s) associated with this object

Outputs Parameters

prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
Default:False
C++ Type:bool
Controllable:No
Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

Material Property Retrieval Parameters

total_swelling_scaling_factor1Scaling factor to be applied to the swelling strain. Used for sensitivity and calibration studies
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scaling factor to be applied to the swelling strain. Used for sensitivity and calibration studies

Advanced: Scaling Factors Parameters

Input Files

(test/tests/solid_mechanics/un_eigenstrains/ad_un_swelling.i)
(test/tests/solid_mechanics/un_eigenstrains/un_swelling.i)

References

R. Matthews S. Ross, M.El Genk. Uranium Nitride Fuel Swelling Correlation. Technical Report UraniumNitrideSwellingCorrelation, University of New Mexico, 1989.

@TECHREPORT{UraniumNitrideSwellingCorrelation,
    author = "S. Ross, M.El Genk, R. Matthews",
    title = "{Uranium Nitride Fuel Swelling Correlation}",
    institution = "University of New Mexico",
    year = "1989",
    number = "UraniumNitrideSwellingCorrelation"
}

(test/tests/solid_mechanics/un_eigenstrains/un_swelling.i)

#
# Test of swelling of UN (1000K)
#
# For this problem, temperature is fixed at 1000K and burnup is increased from 0 to 0.5 over a time of t=0 to t=1.  There are 10 steps.
#
# The results are as follows:
#                                          BISON                          Analytical Swelling Strain Calc
#                 -------------------------------------------------
#
#
#
#| time           | ave_burnup     | avgFuelTemp    | volume         |           volumetric_strain
#+----------------+----------------+----------------+----------------+
#|   0.000000e+00 |   0.000000e+00 |   1.000000e+03 |  -1.000000e+00 |             0.000000e+00
#|   1.000000e-01 |   5.000000e-03 |   1.000000e+03 |  -1.005903e+00 |             5.903435e-03
#|   2.000000e-01 |   1.000000e-02 |   1.000000e+03 |  -1.010494e+00 |             1.049445e-02
#|   3.000000e-01 |   1.500000e-02 |   1.000000e+03 |  -1.014693e+00 |             1.469317e-00
#|   4.000000e-01 |   2.000000e-02 |   1.000000e+03 |  -1.018656e+00 |             1.865583e-02
#|   5.000000e-01 |   2.500000e-02 |   1.000000e+03 |  -1.022452e+00 |             2.245173e-02
#|   6.000000e-01 |   3.000000e-02 |   1.000000e+03 |  -1.026120e+00 |             2.611983e-02
#|   7.000000e-01 |   3.500000e-02 |   1.000000e+03 |  -1.029685e+00 |             2.968494e-02
#|   8.000000e-01 |   4.000000e-02 |   1.000000e+03 |  -1.033164e+00 |             3.316419e-02
#|   9.000000e-01 |   4.500000e-02 |   1.000000e+03 |  -1.036570e+00 |             3.657009e-02
#|   1.000000e+00 |   5.000000e-02 |   1.000000e+03 |  -1.039912e+00 |             3.991212e-02
#

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  [mesh]
    type = FileMeshGenerator
    file = patch.e
  []
[]

[Variables]

  [disp_x]
    order = FIRST
    family = LAGRANGE
  []

  [disp_y]
    order = FIRST
    family = LAGRANGE
  []

  [disp_z]
    order = FIRST
    family = LAGRANGE
  []

[]

[AuxVariables]

  [burnup]
    order = FIRST
    family = LAGRANGE
  []

  [temp]
    order = FIRST
    family = LAGRANGE
    initial_condition = 1000.0
  []

[]

[Functions]
  [burnupFunc]
    type = PiecewiseLinear
    x = '0 1'
    y = '0 0.05'
  []
[]

[AuxKernels]
  [burnup]
    type = FunctionAux
    variable = burnup
    function = burnupFunc
    execute_on = timestep_begin
  []
[]

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [all]
        strain = FINITE
        incremental = true
        eigenstrain_names = 'swelling'
      []
    []
  []
[]

[BCs]
  [bottom_x]
    type = DirichletBC
    variable = disp_x
    boundary = 10
    value = 0.0
  []
  [bottom_y]
    type = DirichletBC
    variable = disp_y
    boundary = 9
    value = 0.0
  []
  [bottom_z]
    type = DirichletBC
    variable = disp_z
    boundary = 14
    value = 0.0
  []
[]

[Materials]
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 1e11
    poissons_ratio = 0.17
  []
  [stress]
    type = ComputeFiniteStrainElasticStress
  []
  [swelling]
    type = UNVolumetricSwellingEigenstrain
    temperature = temp
    burnup = burnup
    density_percent = 95
    eigenstrain_name = swelling
  []
[]

[Executioner]
   type = Transient

  solve_type = 'PJFNK'


   petsc_options = '-snes_ksp_ew '
   petsc_options_iname = '-ksp_gmres_restart -pc_type -pc_hypre_type'
   petsc_options_value = '70 hypre boomeramg'
   l_max_its = 60
   nl_rel_tol = 1e-8
   nl_abs_tol = 1e-10
   l_tol = 1e-5

   start_time = 0.0
   dt = 0.1
   num_steps = 10

[]

[Postprocessors]
  [avgFuelTemp]
    type = ElementAverageValue
    variable = temp
    block = 4
    execute_on = 'initial timestep_end'
  []
  [volume]
    type = InternalVolume
    boundary = 15
    scale_factor = -1.0
    execute_on = 'initial timestep_end'
  []
  [ave_burnup]
    type = ElementAverageValue
    variable = burnup
    block = 4
  []
[]

[Outputs]
  exodus = true
  print_linear_residuals = false
  hide = 'disp_x disp_y disp_z temp'
[]

(test/tests/solid_mechanics/un_eigenstrains/ad_un_swelling.i)

#
# Test of swelling of UN (1000K)
#
# For this problem, temperature is fixed at 1000K and burnup is increased from 0 to 0.5 over a time of t=0 to t=1.  There are 10 steps.
#
# The results are as follows:
#                                          BISON                          Analytical Swelling Strain Calc
#                 -------------------------------------------------
#
#
#
#| time           | ave_burnup     | avgFuelTemp    | volume         |           volumetric_strain
#+----------------+----------------+----------------+----------------+
#|   0.000000e+00 |   0.000000e+00 |   1.000000e+03 |  -1.000000e+00 |             0.000000e+00
#|   1.000000e-01 |   5.000000e-03 |   1.000000e+03 |  -1.005903e+00 |             5.903435e-03
#|   2.000000e-01 |   1.000000e-02 |   1.000000e+03 |  -1.010494e+00 |             1.049445e-02
#|   3.000000e-01 |   1.500000e-02 |   1.000000e+03 |  -1.014693e+00 |             1.469317e-00
#|   4.000000e-01 |   2.000000e-02 |   1.000000e+03 |  -1.018656e+00 |             1.865583e-02
#|   5.000000e-01 |   2.500000e-02 |   1.000000e+03 |  -1.022452e+00 |             2.245173e-02
#|   6.000000e-01 |   3.000000e-02 |   1.000000e+03 |  -1.026120e+00 |             2.611983e-02
#|   7.000000e-01 |   3.500000e-02 |   1.000000e+03 |  -1.029685e+00 |             2.968494e-02
#|   8.000000e-01 |   4.000000e-02 |   1.000000e+03 |  -1.033164e+00 |             3.316419e-02
#|   9.000000e-01 |   4.500000e-02 |   1.000000e+03 |  -1.036570e+00 |             3.657009e-02
#|   1.000000e+00 |   5.000000e-02 |   1.000000e+03 |  -1.039912e+00 |             3.991212e-02
#

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  [mesh]
    type = FileMeshGenerator
    file = patch.e
  []
[]

[AuxVariables]

  [burnup]
    order = FIRST
    family = LAGRANGE
  []

  [temp]
    order = FIRST
    family = LAGRANGE
    initial_condition = 1000.0
  []

[]

[Functions]
  [burnupFunc]
    type = PiecewiseLinear
    x = '0 1'
    y = '0 0.05'
  []
[]

[AuxKernels]
  [burnup]
    type =FunctionAux
    variable = burnup
    function = burnupFunc
    execute_on = timestep_begin
  []
[]

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [all]
        strain = FINITE
        add_variables = true
        incremental = true
        eigenstrain_names = 'swelling'
        use_automatic_differentiation = true
      []
    []
  []
[]

[BCs]
  [bottom_x]
    type = ADDirichletBC
    variable = disp_x
    boundary = 10
    value = 0.0
  []
  [bottom_y]
    type = ADDirichletBC
    variable = disp_y
    boundary = 9
    value = 0.0
  []
  [bottom_z]
    type = ADDirichletBC
    variable = disp_z
    boundary = 14
    value = 0.0
  []
[]

[Materials]
  [elasticity_tensor]
    type = ADComputeIsotropicElasticityTensor
    youngs_modulus = 1e11
    poissons_ratio = 0.17
  []
  [stress]
    type = ADComputeFiniteStrainElasticStress
  []
  [swelling]
    type = ADUNVolumetricSwellingEigenstrain
    temperature = temp
    burnup = burnup
    density_percent = 95
    eigenstrain_name = swelling
  []
[]

[Executioner]
   type = Transient
   solve_type = 'PJFNK'
   petsc_options = '-snes_ksp_ew '
   petsc_options_iname = '-ksp_gmres_restart -pc_type -pc_hypre_type'
   petsc_options_value = '70 hypre boomeramg'
   l_max_its = 60
   nl_rel_tol = 1e-8
   nl_abs_tol = 1e-10
   l_tol = 1e-5

   start_time = 0.0
   dt = 0.1
   num_steps = 10

[]

[Postprocessors]
  [avgFuelTemp]
    type = ElementAverageValue
    variable = temp
    block = 4
    execute_on = 'initial timestep_end'
  []
  [volume]
    type = InternalVolume
    boundary = 15
    scale_factor = -1.0
    execute_on = 'initial timestep_end'
  []
  [ave_burnup]
    type = ElementAverageValue
    variable = burnup
    block = 4
  []
[]

[Outputs]
  exodus = true
  print_linear_residuals = false
  hide = 'disp_x disp_y disp_z temp'
[]

(test/tests/solid_mechanics/un_eigenstrains/un_swelling.i)

#
# Test of swelling of UN (1000K)
#
# For this problem, temperature is fixed at 1000K and burnup is increased from 0 to 0.5 over a time of t=0 to t=1.  There are 10 steps.
#
# The results are as follows:
#                                          BISON                          Analytical Swelling Strain Calc
#                 -------------------------------------------------
#
#
#
#| time           | ave_burnup     | avgFuelTemp    | volume         |           volumetric_strain
#+----------------+----------------+----------------+----------------+
#|   0.000000e+00 |   0.000000e+00 |   1.000000e+03 |  -1.000000e+00 |             0.000000e+00
#|   1.000000e-01 |   5.000000e-03 |   1.000000e+03 |  -1.005903e+00 |             5.903435e-03
#|   2.000000e-01 |   1.000000e-02 |   1.000000e+03 |  -1.010494e+00 |             1.049445e-02
#|   3.000000e-01 |   1.500000e-02 |   1.000000e+03 |  -1.014693e+00 |             1.469317e-00
#|   4.000000e-01 |   2.000000e-02 |   1.000000e+03 |  -1.018656e+00 |             1.865583e-02
#|   5.000000e-01 |   2.500000e-02 |   1.000000e+03 |  -1.022452e+00 |             2.245173e-02
#|   6.000000e-01 |   3.000000e-02 |   1.000000e+03 |  -1.026120e+00 |             2.611983e-02
#|   7.000000e-01 |   3.500000e-02 |   1.000000e+03 |  -1.029685e+00 |             2.968494e-02
#|   8.000000e-01 |   4.000000e-02 |   1.000000e+03 |  -1.033164e+00 |             3.316419e-02
#|   9.000000e-01 |   4.500000e-02 |   1.000000e+03 |  -1.036570e+00 |             3.657009e-02
#|   1.000000e+00 |   5.000000e-02 |   1.000000e+03 |  -1.039912e+00 |             3.991212e-02
#

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  [mesh]
    type = FileMeshGenerator
    file = patch.e
  []
[]

[Variables]

  [disp_x]
    order = FIRST
    family = LAGRANGE
  []

  [disp_y]
    order = FIRST
    family = LAGRANGE
  []

  [disp_z]
    order = FIRST
    family = LAGRANGE
  []

[]

[AuxVariables]

  [burnup]
    order = FIRST
    family = LAGRANGE
  []

  [temp]
    order = FIRST
    family = LAGRANGE
    initial_condition = 1000.0
  []

[]

[Functions]
  [burnupFunc]
    type = PiecewiseLinear
    x = '0 1'
    y = '0 0.05'
  []
[]

[AuxKernels]
  [burnup]
    type = FunctionAux
    variable = burnup
    function = burnupFunc
    execute_on = timestep_begin
  []
[]

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [all]
        strain = FINITE
        incremental = true
        eigenstrain_names = 'swelling'
      []
    []
  []
[]

[BCs]
  [bottom_x]
    type = DirichletBC
    variable = disp_x
    boundary = 10
    value = 0.0
  []
  [bottom_y]
    type = DirichletBC
    variable = disp_y
    boundary = 9
    value = 0.0
  []
  [bottom_z]
    type = DirichletBC
    variable = disp_z
    boundary = 14
    value = 0.0
  []
[]

[Materials]
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 1e11
    poissons_ratio = 0.17
  []
  [stress]
    type = ComputeFiniteStrainElasticStress
  []
  [swelling]
    type = UNVolumetricSwellingEigenstrain
    temperature = temp
    burnup = burnup
    density_percent = 95
    eigenstrain_name = swelling
  []
[]

[Executioner]
   type = Transient

  solve_type = 'PJFNK'


   petsc_options = '-snes_ksp_ew '
   petsc_options_iname = '-ksp_gmres_restart -pc_type -pc_hypre_type'
   petsc_options_value = '70 hypre boomeramg'
   l_max_its = 60
   nl_rel_tol = 1e-8
   nl_abs_tol = 1e-10
   l_tol = 1e-5

   start_time = 0.0
   dt = 0.1
   num_steps = 10

[]

[Postprocessors]
  [avgFuelTemp]
    type = ElementAverageValue
    variable = temp
    block = 4
    execute_on = 'initial timestep_end'
  []
  [volume]
    type = InternalVolume
    boundary = 15
    scale_factor = -1.0
    execute_on = 'initial timestep_end'
  []
  [ave_burnup]
    type = ElementAverageValue
    variable = burnup
    block = 4
  []
[]

[Outputs]
  exodus = true
  print_linear_residuals = false
  hide = 'disp_x disp_y disp_z temp'
[]

Description
Example Input Syntax
Input Parameters
Input Files
References