HighBurnup LWRNuclearMaterials

This tutorial is currently under development.

schooltip:Try the Action Wizard

There is a recently developed Action Wizard that will create the NuclearMaterials block based on user input. Read through this tutorial, then use it to experiment with the action.

A pre-existing assessment case already exists in:

[NuclearMaterials<<<{"href": "../../syntax/NuclearMaterials/index.html"}>>>]
  generate_output<<<{"description": "Add scalar quantity output for stress and/or strain"}>>> = 'stress_xx stress_yy stress_zz vonmises_stress'
  fission_operation<<<{"description": "Type of fission occuring in this simulation."}>>> = HighBurnup
  physics<<<{"description": "Type(s) of physics used on this block."}>>> = 'Mechanics Thermal'
  initial_temperature<<<{"description": "Initial temperature in Kelvins."}>>> = 300
  stress_free_temperature<<<{"description": "Reference temperature for thermal eigenstrain calculation"}>>> = 300
  strain<<<{"description": "Strain formulation"}>>> = FINITE
  [UO2<<<{"href": "../../syntax/NuclearMaterials/UO2/index.html"}>>>]
    [fuel]
      block<<<{"description": "The list of ids of the blocks (subdomain) that the stress divergence kernels will be applied to"}>>> = pellet_type_1
      uo2_models<<<{"description": "Type(s) of physics models used on this block.   The choices are: Burnup Elastic Creep Relocation Swelling ThermalExpansion HighBurnupStructureFormation"}>>> = 'Burnup Elastic Relocation Swelling ThermalExpansion'
      isotopes<<<{"description": "Fuel isotopes: Gd155 Gd157 U235 U238 Pu239 Pu240 Pu241 Pu242.  Number of entries must match number of entries in isotope_fractions."}>>> = 'U238 U235'
      isotope_fractions<<<{"description": "The isotope fractions associated with the 'isotopes' input line.  Must sum to 1.0."}>>> = '0.9293 0.0707'
      burnup_relocation_stop<<<{"description": "Burnup at which relocation strain stops (in FIMA)"}>>> = 0.035
      fuel_pin_geometry<<<{"description": "Name of the UserObject that reads the pin geometry from the mesh."}>>> = pin_geometry
      rod_ave_lin_pow<<<{"description": "The rod power history function."}>>> = power_history
      axial_power_profile<<<{"description": "The axial power peaking function."}>>> = axial_peaking_factors
      additional_generate_output<<<{"description": "Add scalar quantity output for stress and/or strain (will be appended to the list in `generate_output`)"}>>> = hydrostatic_stress
      fuel_volume_ratio<<<{"description": "Reduction factor for deviation from right circular cylinder representation of the fuel mesh.  The ratio of actual volume to right circular cylinder volume."}>>> = 1.0
      # The mesh is not a right cyclinder due to
      # chamfering and dishing, following the
      # established HBEP_BK363.i results, this
      # is kept as fuel_volume_ratio=1 to match
      # test results
      extra_vector_tags<<<{"description": "The tag names for extra vectors that residual data should be saved into"}>>> = 'ref'
    []
  []
  [ZirconiumAlloy<<<{"href": "../../syntax/NuclearMaterials/ZirconiumAlloy/index.html"}>>>]
    [clad]
      block<<<{"description": "The list of ids of the blocks (subdomain) that the stress divergence kernels will be applied to"}>>> = clad
      cladding_models<<<{"description": "Type(s) of physics models used on this block.   The choices are: Creep Elastic Plasticity ZryOxidation IrradiationGrowth ThermalExpansion ZrPhase ZryCladdingFailure"}>>> = 'Elastic Creep IrradiationGrowth ThermalExpansion'
      additional_generate_output<<<{"description": "Add scalar quantity output for stress and/or strain (will be appended to the list in `generate_output`)"}>>> = 'creep_strain_xx creep_strain_yy
        creep_strain_xy hoop_creep_strain'
      flux_function<<<{"description": "Function which contains the neutron flux data."}>>> = flux
      extra_vector_tags<<<{"description": "The tag names for extra vectors that residual data should be saved into"}>>> = 'ref'
    []
  []
[]
(assessment/LWR/validation/HBEP/analysis/BK363/HBEP_BK363_action.i)