MetallicFuel Actions

For metallic fuel simulations there are 2 distinct actions that can be used. Information on these actions can be found below:

MetallicFuel Example

This tutorial is currently under development.

A pre-existing example already exists in:

[NuclearMaterials<<<{"href": "../syntax/NuclearMaterials/index.html"}>>>]
  physics<<<{"description": "Type(s) of physics used on this block."}>>> = 'Mechanics Thermal'
  fission_operation<<<{"description": "Type of fission occuring in this simulation."}>>> = Normal
  generate_output<<<{"description": "Add scalar quantity output for stress and/or strain"}>>> = 'stress_xx stress_yy stress_zz vonmises_stress
    hydrostatic_stress creep_strain_xx creep_strain_yy creep_strain_zz
    elastic_strain_xx elastic_strain_yy elastic_strain_zz strain_xx
    strain_yy strain_zz hoop_stress'
  add_variables<<<{"description": "Add the displacement variables"}>>> = true
  extra_vector_tags<<<{"description": "The tag names for extra vectors that residual data should be saved into"}>>> = 'ref'
  strain<<<{"description": "Strain formulation"}>>> = FINITE
  decomposition_method<<<{"description": "Methods to calculate the finite strain and rotation increments"}>>> = TaylorExpansion
  stress_free_temperature<<<{"description": "Reference temperature for thermal eigenstrain calculation"}>>> = 295.0
  initial_temperature<<<{"description": "Initial temperature in Kelvins."}>>> = 298
  incremental<<<{"description": "Use incremental or total strain"}>>> = true
  [UPuZr<<<{"href": "../syntax/NuclearMaterials/UPuZr/index.html"}>>>]
    [fuel]
      block<<<{"description": "The list of ids of the blocks (subdomain) that the stress divergence kernels will be applied to"}>>> = pellet
      additional_generate_output<<<{"description": "Add scalar quantity output for stress and/or strain (will be appended to the list in `generate_output`)"}>>> = 'volumetric_strain'
      upuzr_models<<<{"description": "Type(s) of physics models used on this block.   The choices are: FIPD Burnup Elastic Creep Swelling ThermalExpansion"}>>> = 'Elastic Burnup Creep Swelling ThermalExpansion'
      rod_linear_power<<<{"description": "Name of the function describing rod linear power."}>>> = power_history
      axial_power_profile<<<{"description": "Name of the function describing axial power profile."}>>> = axial_peaking_factors
      pellet_radius<<<{"description": "Fuel slug radius."}>>> = 2.195e-03
      bubble_number_density<<<{"description": "Material property name for the number density of intragranular bubbles, [bubbles/m^3]."}>>> = 5e17
      interconnection_initiating_porosity<<<{"description": "Porosity at which fission gas release starts."}>>> = 0.29
      interconnection_terminating_porosity<<<{"description": "Porosity at which fission gas release finishes."}>>> = 0.31
      initial_X_Pu<<<{"description": "Initial atom fraction of plutonium."}>>> = 0.225
      initial_X_Zr<<<{"description": "Initial atom fraction of zirconium."}>>> = 0.163
      critical_porosity<<<{"description": "Porosity at which fission gas release initiates. This needs to be between the interconnection_initiating_porosity and interconnection_terminating_porosity in UPuZrVolumetricSwellingEigenstrain."}>>> = 0.30
      fractional_fgr_initial<<<{"description": "Fraction of fission gas produced that is released once the critical porosity is reached."}>>> = 0.4
      fractional_fgr_post<<<{"description": "Fraction of fission gas produced that is released after interconnection."}>>> = 0.8
      density<<<{"description": "Initial fuel density [kg/m3]."}>>> = 15800
      output_properties<<<{"description": "List of material properties, from this material, to output (outputs must also be defined to an output type)."}>>> = 'porosity gaseous_porosity'
      max_inelastic_increment<<<{"description": "Max inelastic increment."}>>> = 1e-2
    []
  []
  [HT9<<<{"href": "../syntax/NuclearMaterials/HT9/index.html"}>>>]
    [clad]
      block<<<{"description": "The list of ids of the blocks (subdomain) that the stress divergence kernels will be applied to"}>>> = clad
      additional_generate_output<<<{"description": "Add scalar quantity output for stress and/or strain (will be appended to the list in `generate_output`)"}>>> = 'hoop_creep_strain hoop_elastic_strain'
      ht9_models<<<{"description": "Type(s) of physics models used on this block.   The choices are: FIPD Elastic Creep ThermalExpansion"}>>> = 'Elastic Creep ThermalExpansion'
      fast_flux_factor<<<{"description": "Constant multiplied against the function, rod average linear power, or q_variable."}>>> = 2.47e19
    []
  []
[]
(examples/NuclearMaterialActions/MetallicFuel/x441_group_A_nominal_action.i)