- fuel_pin_geometryName of fuel pin geometry UserObject which is required for Layered1D or Layered2D calculations.
C++ Type:UserObjectName
Controllable:No
Description:Name of fuel pin geometry UserObject which is required for Layered1D or Layered2D calculations.
- mat_propfis_gas_relThe name of the material property
Default:fis_gas_rel
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:The name of the material property
LayeredElementIntegralFisGasReleasePostprocessor
Reports the fission gas that is released to the plenum in moles for Layered1D or Layered2D geometries.
Description
LayeredElementIntegralFisGasReleasePostprocessor is used to report the fission gas released for Layered1D and Layered2D geometries. This postprocessor simply takes the value of fission gas released calculated in the material model (moles/volume), integrates over the volume of the domain, and reports the total value in moles.
An associated postprocessor is:
Example Input Syntax
[Postprocessors<<<{"href": "../../syntax/Postprocessors/index.html"}>>>]
[gas_released]
type = LayeredElementIntegralFisGasReleasePostprocessor<<<{"description": "Reports the fission gas that is released to the plenum in moles for Layered1D or Layered2D geometries.", "href": "LayeredElementIntegralFisGasReleasePostprocessor.html"}>>>
fuel_pin_geometry<<<{"description": "Name of fuel pin geometry UserObject which is required for Layered1D or Layered2D calculations."}>>> = pin_geometry
[]
[]Input Parameters
- 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
Optional Parameters
- allow_duplicate_execution_on_initialFalseIn the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
Default:False
C++ Type:bool
Controllable:No
Description:In the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
- execute_onTIMESTEP_ENDThe list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.
Default:TIMESTEP_END
C++ Type:ExecFlagEnum
Options:XFEM_MARK, NONE, INITIAL, LINEAR, NONLINEAR_CONVERGENCE, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM, TRANSFER
Controllable:No
Description:The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.
- execution_order_group0Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
Default:0
C++ Type:int
Controllable:No
Description:Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
- force_postauxFalseForces the UserObject to be executed in POSTAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in POSTAUX
- force_preauxFalseForces the UserObject to be executed in PREAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREAUX
- force_preicFalseForces the UserObject to be executed in PREIC during initial setup
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREIC during initial setup
Execution Scheduling 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
- outputsVector of output names where you would like to restrict the output of variables(s) associated with this object
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
- 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
- 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
Input Files
- (assessment/metallic_fuel/EBRII/X441/analysis/group_A/x441_1_5D_A.i)
- (assessment/metallic_fuel/EBRII/X441/analysis/group_E/x441_1_5D_E.i)
- (assessment/metallic_fuel/EBRII/X441/analysis/group_H/x441_1_5D_H.i)
- (test/tests/upuzr_fission_gas_release/nonad/exact_1D.i)
- (assessment/metallic_fuel/EBRII/X441/analysis/group_G/x441_1_5D_G.i)
- (assessment/metallic_fuel/EBRII/X441/analysis/group_D/x441_1_5D_D.i)
- (assessment/metallic_fuel/EBRII/X441/analysis/x441_base_1_5D.i)
- (test/tests/fission_gas_1d/fgr_1d.i)
- (assessment/metallic_fuel/EBRII/X441/analysis/group_B/x441_1_5D_B.i)
- (assessment/metallic_fuel/EBRII/X441/analysis/group_F/x441_1_5D_F.i)
- (test/tests/upuzr_fission_gas_release/ad/exact_1D.i)
- (assessment/metallic_fuel/EBRII/X441/analysis/group_C/x441_1_5D_C.i)
(test/tests/fission_gas_1d/fgr_1d.i)
#This tests the fission gas produced and released postprocessors
#and for 1d and 2d models with the same volume.
#The reported values should be the same for the
#1d and 2d cases. See corresponding excel calculation
#in this directory.
[Problem]
type = ReferenceResidualProblem
reference_vector = 'ref'
extra_tag_vectors = 'ref'
[]
[UserObjects]
[pin_geometry]
type = Layered1DFuelPinGeometry
include_clad = false
mesh_generator = layered1D_mesh
[]
[]
[Mesh]
coord_type = RZ
[layered1D_mesh]
type = Layered1DMeshGenerator
include_clad = false
include_plenum = false
slices_per_block = 1
fuel_height = 3e-3
pellet_outer_radius = 3e-3
pellet_bottom_coor = 0.0
pellet_mesh_density = customize
nx_p = 4
elem_type = EDGE2
[]
[]
[Functions]
[Temp_func]
type = ParsedFunction
expression = '1700'
[]
[Fiss_func]
type = ParsedFunction
expression = '2.5e19'
[]
[Porosity_func]
type = PiecewiseLinear
x = '0 1e7'
y = '0.15 0.241'
[]
[]
[Variables]
[temp]
order = FIRST
family = LAGRANGE
initial_condition = 1700
[]
[]
[AuxVariables]
[porosity_material_aux]
order = CONSTANT
family = MONOMIAL
[]
[]
[Kernels]
[heat]
type = HeatConduction
variable = temp
[]
[]
[AuxKernels]
[porosity_material_aux]
type = MaterialRealAux
variable = porosity_material_aux
property = porosity
execute_on = 'initial timestep_begin'
[]
[]
[BCs]
[temp_bc]
type = FunctionDirichletBC
variable = temp
function = Temp_func
boundary = 10
[]
[]
[Materials]
[fission_rate]
type = GenericFunctionMaterial
prop_names = 'fission_rate thermal_conductivity porosity'
prop_values = 'Fiss_func 5 Porosity_func'
outputs = all
[]
[fission_gas_behavior]
type = UPuZrFissionGasRelease
fission_rate = fission_rate
fractional_yield = 0.3017
[]
[]
[Executioner]
type = Transient
solve_type = 'PJFNK'
l_tol = 1e-4
nl_abs_tol = 1e-5
nl_rel_tol = 1e-5
start_time = 0.0
dt = 5e6
end_time = 1e7
[]
[Postprocessors]
[gas_produced]
type = LayeredElementIntegralFisGasProducePostprocessor
fuel_pin_geometry = pin_geometry
[]
[gas_released]
type = LayeredElementIntegralFisGasReleasePostprocessor
fuel_pin_geometry = pin_geometry
[]
[porosity]
type = ElementAverageValue
variable = porosity_material_aux
[]
[]
[Outputs]
[csv]
type = CSV
execute_on = final
file_base = fgr_1d_out
[]
[]
(assessment/metallic_fuel/EBRII/X441/analysis/group_A/x441_1_5D_A.i)
initial_fuel_density = 15800.0
[GlobalParams]
order = FIRST
family = LAGRANGE
energy_per_fission = 3.2e-11 # J/fission
displacements = 'disp_x'
[]
[Mesh]
coord_type = RZ
# Nominal Design Geometric Parameters (X441)
[layered1D_mesh]
type = Layered1DMeshGenerator
clad_thickness = ${clad_thickness}
pellet_outer_radius = ${pellet_outer_radius}
fuel_height = ${pellet_height}
plenum_height = ${clad_top_gap_height}
clad_gap_width = ${clad_gap_width}
# meshing parameters
clad_mesh_density = customize
pellet_mesh_density = customize
nx_p = 6
nx_c = 4
slices_per_block = 10
elem_type = EDGE2
[]
# mesh options
patch_update_strategy = auto
partitioner = centroid
centroid_partitioner_direction = y
[]
[UserObjects]
[pin_geometry]
type = Layered1DFuelPinGeometry
mesh_generator = layered1D_mesh
[]
[]
[Variables]
[temp]
initial_condition = 298
[]
[]
[AuxVariables]
[disp_y]
[]
[disp_z]
[]
# Aux variables for output
[gap_cond]
order = CONSTANT
family = MONOMIAL
[]
[coolant_htc]
order = CONSTANT
family = MONOMIAL
[]
[cumulative_damage_index]
order = CONSTANT
family = MONOMIAL
[]
[solid_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[gas_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[volumetric_strain]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[hoop_stress]
order = CONSTANT
family = MONOMIAL
[]
[hoop_creep_strain]
order = CONSTANT
family = MONOMIAL
[]
[hoop_elastic_strain]
order = CONSTANT
family = MONOMIAL
[]
[total_hoop_strain]
order = CONSTANT
family = MONOMIAL
[]
[]
[Functions]
[power_history]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '0 44722 44722 0'
[]
[coolant_press_ramp]
type = PiecewiseLinear
x = '0 42000400'
y = '1.0 1.0'
[]
[coolant_temp_ramp]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '298.0 648.0 648.0 350.0'
[]
[axial_peaking_factors]
type = PowerPeakingFunction
fit = EBRII_ROW_4
pellet_length = ${pellet_height}
pellet_y_start = ${pellet_y_start}
[]
[clad_axial_pressure]
type = CladdingAxialPressureFunction
plenum_pressure = plenum_pressure
coolant_pressure = coolant_press_ramp
coolant_pressure_scaling_factor = 0.151e6
fuel_pin_geometry = pin_geometry
[]
[fuel_axial_pressure]
type = ParsedFunction
expression = plenum_pressure
symbol_names = plenum_pressure
symbol_values = plenum_pressure
[]
[]
[Physics/SolidMechanics/Layered1D]
[fuel]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = fuel_axial_pressure
generate_output = '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'
block = fuel
eigenstrain_names = 'fuel_thermal_strain gas_swelling_eigenstrain solid_swelling_eigenstrain'
mesh_generator = layered1D_mesh
[]
[clad]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = clad_axial_pressure
generate_output = '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'
block = clad
eigenstrain_names = 'clad_thermal_eigenstrain'
mesh_generator = layered1D_mesh
[]
[]
[Kernels]
# Define kernels for the various terms in the PDE system
[heat]
type = HeatConduction
variable = temp
[]
[heat_ie]
type = HeatConductionTimeDerivative
variable = temp
[]
[heat_source]
type = NeutronHeatSource
variable = temp
block = fuel
fission_rate = fission_rate
[]
[]
[AuxKernels]
[conductance]
type = MaterialRealAux
property = gap_conductance
variable = gap_cond
boundary = 10
[]
[gas_swell]
type = MaterialRealAux
variable = gas_swell
property = gas_swelling
execute_on = timestep_end
[]
[solid_swell]
type = MaterialRealAux
variable = solid_swell
property = solid_swelling
execute_on = timestep_end
[]
[volumetric_strain]
type = RankTwoScalarAux
rank_two_tensor = total_strain
variable = volumetric_strain
scalar_type = VolumetricStrain
execute_on = timestep_end
block = fuel
[]
[hoop_stress]
type = RankTwoAux
rank_two_tensor = stress
variable = hoop_stress
index_j = 2
index_i = 2
execute_on = timestep_end
[]
[hoop_creep_strain]
type = RankTwoAux
rank_two_tensor = creep_strain
variable = hoop_creep_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[hoop_elastic_strain]
type = RankTwoAux
rank_two_tensor = elastic_strain
variable = hoop_elastic_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[total_hoop_strain]
type = RankTwoAux
rank_two_tensor = total_strain
variable = total_hoop_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[]
[Contact]
[pellet_clad_mechanical]
primary = 5
secondary = 10
penalty = 1e12
model = frictionless
formulation = kinematic
[]
[]
[ThermalContact]
[thermal_contact]
type = GapHeatTransfer
variable = temp
primary = 5
secondary = 10
quadrature = true
gap_conductivity = 61.0
min_gap = ${clad_gap_width}
[]
[]
[BCs]
[no_x_all]
type = DirichletBC
variable = disp_x
boundary = 12
value = 0.0
[]
[Pressure]
[coolantPressure]
boundary = '2'
function = coolant_press_ramp
[]
[]
[PlenumPressure]
[plenumPressure]
boundary = 9
initial_pressure = 0.084e6 # Pa
startup_time = 0
R = 8.3143
temperature = ave_temp_interior
volume = gas_volume
output = plenum_pressure
material_input = fis_gas_released
[]
[]
[]
[CoolantChannel]
[convective_clad_surface]
boundary = '2'
variable = temp
inlet_temperature = coolant_temp_ramp
inlet_pressure = coolant_press_ramp
inlet_massflux = 5261.5 # kg/m^2-sec
coolant_material = sodium
rod_diameter = 5.84e-3 # m
rod_pitch = 7.48e-3 # m (Pitch-to-diameter Ratio = 1.28)
linear_heat_rate = power_history
axial_power_profile = axial_peaking_factors
subchannel_geometry = triangular
[]
[]
[Materials]
[fission_rate]
type = UPuZrFissionRate
rod_linear_power = power_history
axial_power_profile = axial_peaking_factors
pellet_radius = ${pellet_outer_radius}
X_Zr = ${X_Zr}
X_Pu_function = ${X_Pu}
block = fuel
outputs = all
[]
[burnup]
type = UPuZrBurnup
initial_X_Zr = ${X_Zr}
initial_X_Pu = ${X_Pu}
density = ${initial_fuel_density}
block = fuel
outputs = all
[]
[fuel_elasticity_tensor]
type = UPuZrElasticityTensor
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
block = fuel
temperature = temp
[]
[fuel_inelastic_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'fuel_upuzrcreep'
block = fuel
[]
[fuel_upuzrcreep]
type = UPuZrCreepUpdate
block = fuel
temperature = temp
porosity = porosity
max_inelastic_increment = 1e-3
[]
[fuel_thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = fuel
thermal_expansion_coeff = 1.18e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = fuel_thermal_strain
[]
[gas_swelling]
type = UPuZrGaseousEigenstrain
temperature = temp
eigenstrain_name = gas_swelling_eigenstrain
initial_porosity = 0.0
bubble_number_density = 1e20
outputs = all
output_properties = 'porosity gaseous_porosity'
block = fuel
[]
[solid_swelling]
type = BurnupDependentEigenstrain
eigenstrain_name = solid_swelling_eigenstrain
block = fuel
swelling_name = 'solid_swelling'
[]
[metal_fuel_thermal]
type = UPuZrThermal
block = fuel
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
spheat_model = savage
thcond_model = lanl
porosity = porosity
temperature = temp
[]
[fuel_density]
type = StrainAdjustedDensity
block = fuel
strain_free_density = ${initial_fuel_density}
[]
[fission_gas_behavior]
type = UPuZrFissionGasRelease
block = fuel
fission_rate = fission_rate
[]
[clad_elasticity_tensor]
type = ComputeIsotropicElasticityTensor
youngs_modulus = 1.88e11
poissons_ratio = 0.236
block = clad
[]
[clad_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'clad_ht9creep'
block = clad
[]
[fast_flux]
type = FastNeutronFlux
block = clad
factor = 2.47e19
[]
[clad_ht9creep]
type = HT9CreepUpdate
block = clad
temperature = temp
[]
[thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = clad
thermal_expansion_coeff = 1.2e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = clad_thermal_eigenstrain
[]
[clad_thermal]
type = HT9Thermal
block = clad
temperature = temp
[]
[clad_density]
type = StrainAdjustedDensity
block = clad
strain_free_density = 7874.0
[]
[]
[Preconditioning]
[SMP]
type = SMP
full = true
[]
[]
[Executioner]
type = Transient
solve_type = 'PJFNK'
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -ksp_gmres_restart'
petsc_options_value = 'lu superlu_dist 51'
line_search = 'none'
l_max_its = 60
l_tol = 8e-3
nl_max_its = 40
nl_rel_tol = 5e-4
nl_abs_tol = 1e-7
end_time = 42000400
dtmin = 100
dtmax = 5e5
[Quadrature]
order = fifth
side_order = seventh
[]
[TimeStepper]
type = IterationAdaptiveDT
timestep_limiting_postprocessor = creep_timestep
dt = 1e2
time_t = '0 1e5 41990400 42000400'
time_dt = '1e2 1e2 1e2 1e2'
iteration_window = 4
optimal_iterations = 10
[]
[]
[Postprocessors]
[ave_temp_interior] # average temperature of the cladding interior and all pellet exteriors
type = LayeredSideAverageValuePostprocessor
boundary = 9
variable = temp
execute_on = 'initial linear'
fuel_pin_geometry = pin_geometry
[]
[clad_inner_vol] # volume inside of cladding
type = LayeredInternalVolumePostprocessor
boundary = 7
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[pellet_volume] # fuel pellet total volume
type = LayeredInternalVolumePostprocessor
boundary = 8
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[avg_clad_temp] # average temperature of cladding interior
type = LayeredSideAverageValuePostprocessor
boundary = 7
variable = temp
fuel_pin_geometry = pin_geometry
execute_on = 'initial linear'
[]
[gas_volume]
type = LayeredInternalVolumePostprocessor
boundary = 9
execute_on = 'initial linear'
component = 0
out_of_plane_strain = strain_yy
fuel_pin_geometry = pin_geometry
addition = ${gas_addition}
[]
[flux_from_clad] # area integrated heat flux from the cladding
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 5
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[flux_from_fuel] # area integrated heat flux from the fuel
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 10
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[rod_total_power]
type = LayeredElementIntegralPowerPostprocessor
variable = temp
fission_rate = fission_rate
block = fuel
fuel_pin_geometry = pin_geometry
[]
[approx_FCT]
type = AverageNodalVariableValue
boundary = 12
variable = temp
[]
[max_approx_FCT]
type = TimeExtremeValue
value_type = max
postprocessor = approx_FCT
[]
[ave_FST]
type = SideAverageValue
boundary = 10
variable = temp
[]
[max_ave_FST]
type = TimeExtremeValue
value_type = max
postprocessor = ave_FST
[]
[ave_CIT]
type = SideAverageValue
boundary = 5
variable = temp
[]
[max_ave_CIT]
type = TimeExtremeValue
value_type = max
postprocessor = ave_CIT
[]
[peak_clad_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = clad
[]
[peak_fuel_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = fuel
[]
[max_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = max
block = fuel
[]
[min_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = min
block = fuel
[]
[peak_porosity]
type = ElementExtremeValue
variable = porosity
value_type = max
block = fuel
[]
[clad_fuel_gap]
type = NodalExtremeValue
variable = penetration
boundary = 10
[]
[max_cont_press]
type = NodalExtremeValue
variable = contact_pressure
boundary = 10
[]
[LHGR_W_per_cm]
type = FunctionValuePostprocessor
function = power_history
scale_factor = 0.01
[]
[average_burnup]
type = ElementAverageValue
block = fuel
variable = burnup
[]
[fis_gas_produced]
type = LayeredElementIntegralFisGasProducePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
[]
[fis_gas_released]
type = LayeredElementIntegralFisGasReleasePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
execute_on = 'initial timestep_end'
[]
[creep_timestep]
type = MaterialTimeStepPostprocessor
block = fuel
[]
[hydrostatic_stress]
type = ElementAverageValue
variable = hydrostatic_stress
execute_on = 'initial timestep_end'
block = fuel
[]
[solid_swelling]
type = ElementAverageValue
variable = solid_swell
block = fuel
[]
[gas_swelling]
type = ElementAverageValue
variable = gas_swell
block = fuel
[]
[volumetric_strain]
type = ElementAverageValue
variable = volumetric_strain
block = fuel
[]
[fission_rate]
type = ElementAverageValue
variable = fission_rate
block = fuel
[]
[porosity]
type = ElementAverageValue
variable = porosity
block = fuel
[]
[gaseous_porosity]
type = ElementAverageValue
variable = gaseous_porosity
block = fuel
[]
[fis_gas_percent]
type = FGRPercent
fission_gas_released = fis_gas_released
fission_gas_generated = fis_gas_produced
[]
[max_clad_hoop_creep]
type = ElementExtremeValue
value_type = max
block = clad
variable = hoop_creep_strain
[]
[max_total_hoop_strain]
type = ElementExtremeValue
value_type = max
block = clad
variable = total_hoop_strain
[]
[]
[VectorPostprocessors]
[clad_x_disp]
type = NodalValueSampler
variable = disp_x
boundary = 2
sort_by = y
outputs = 'vec1'
[]
[fuel_cl_temp]
type = NodalValueSampler
variable = temp
boundary = 12
sort_by = y
outputs = 'vec2'
[]
[fuel_surf_temp]
type = NodalValueSampler
variable = temp
boundary = 10
sort_by = y
outputs = 'vec3'
[]
[clad_inn_temp]
type = NodalValueSampler
variable = temp
boundary = 5
sort_by = y
outputs = 'vec4'
[]
[clad_out_temp]
type = NodalValueSampler
variable = temp
boundary = 2
sort_by = y
outputs = 'vec5'
[]
[fuel_surf_disp_x]
type = NodalValueSampler
variable = disp_x
boundary = 10
sort_by = y
outputs = 'vec7'
[]
[]
[PerformanceMetricOutputs]
outputs = performance_metrics_file
[]
[Outputs]
time_step_interval = 10
color = true
exodus = true
perf_graph = true
csv = true
sync_times = '1e3 5e3 1e4 5e4 1e5 5e6 1e6 5e6 1e7 2e7 3e7 4e7 41990400 41991000'
file_base = x441_${group_name}_1_5D
[out2]
type = CSV
file_base = x441_${group_name}_1_5D_out2
time_step_interval = 1
[]
[console]
type = Console
max_rows = 25
time_step_interval = 1
output_linear = true
[]
[chkfile]
type = CSV
file_base = x441_${group_name}_1_5D_chkfile
show = 'max_approx_FCT max_ave_FST max_ave_CIT average_burnup fis_gas_percent max_clad_hoop_creep max_total_hoop_strain'
execute_on = 'FINAL'
[]
[performance_metrics_file]
type = CSV
file_base = x441_${group_name}_1_5D_performance_metrics
show = 'simulation_alive_time number_linear_iterations number_nonlinear_iterations time_step_size total_linear_iterations total_nonlinear_iterations physical_memory_use number_dofs number_nonlinear_variables residual_compute_time jacobian_compute_time'
[]
[vec1]
type = CSV
file_base = x441_${group_name}_1_5D_vec1
execute_on = 'FINAL'
[]
[vec2]
type = CSV
file_base = x441_${group_name}_1_5D_vec2
execute_on = 'FINAL'
[]
[vec3]
type = CSV
file_base = x441_${group_name}_1_5D_vec3
execute_on = 'FINAL'
[]
[vec4]
type = CSV
file_base = x441_${group_name}_1_5D_vec4
execute_on = 'FINAL'
[]
[vec5]
type = CSV
file_base = x441_${group_name}_1_5D_vec5
execute_on = 'FINAL'
[]
[vec7]
type = CSV
file_base = x441_${group_name}_1_5D_vec7
execute_on = 'FINAL'
[]
[]
[Debug]
show_var_residual = 'disp_x temp'
show_var_residual_norms = true
[]
(assessment/metallic_fuel/EBRII/X441/analysis/group_E/x441_1_5D_E.i)
initial_fuel_density = 15800.0
[GlobalParams]
order = FIRST
family = LAGRANGE
energy_per_fission = 3.2e-11 # J/fission
displacements = 'disp_x'
[]
[Mesh]
coord_type = RZ
# Nominal Design Geometric Parameters (X441)
[layered1D_mesh]
type = Layered1DMeshGenerator
clad_thickness = ${clad_thickness}
pellet_outer_radius = ${pellet_outer_radius}
fuel_height = ${pellet_height}
plenum_height = ${clad_top_gap_height}
clad_gap_width = ${clad_gap_width}
# meshing parameters
clad_mesh_density = customize
pellet_mesh_density = customize
nx_p = 6
nx_c = 4
slices_per_block = 10
elem_type = EDGE2
[]
# mesh options
patch_update_strategy = auto
partitioner = centroid
centroid_partitioner_direction = y
[]
[UserObjects]
[pin_geometry]
type = Layered1DFuelPinGeometry
mesh_generator = layered1D_mesh
[]
[]
[Variables]
[temp]
initial_condition = 298
[]
[]
[AuxVariables]
[disp_y]
[]
[disp_z]
[]
# Aux variables for output
[gap_cond]
order = CONSTANT
family = MONOMIAL
[]
[coolant_htc]
order = CONSTANT
family = MONOMIAL
[]
[cumulative_damage_index]
order = CONSTANT
family = MONOMIAL
[]
[solid_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[gas_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[volumetric_strain]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[hoop_stress]
order = CONSTANT
family = MONOMIAL
[]
[hoop_creep_strain]
order = CONSTANT
family = MONOMIAL
[]
[hoop_elastic_strain]
order = CONSTANT
family = MONOMIAL
[]
[total_hoop_strain]
order = CONSTANT
family = MONOMIAL
[]
[]
[Functions]
[power_history]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '0 44722 44722 0'
[]
[coolant_press_ramp]
type = PiecewiseLinear
x = '0 42000400'
y = '1.0 1.0'
[]
[coolant_temp_ramp]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '298.0 648.0 648.0 350.0'
[]
[axial_peaking_factors]
type = PowerPeakingFunction
fit = EBRII_ROW_4
pellet_length = ${pellet_height}
pellet_y_start = ${pellet_y_start}
[]
[clad_axial_pressure]
type = CladdingAxialPressureFunction
plenum_pressure = plenum_pressure
coolant_pressure = coolant_press_ramp
coolant_pressure_scaling_factor = 0.151e6
fuel_pin_geometry = pin_geometry
[]
[fuel_axial_pressure]
type = ParsedFunction
expression = plenum_pressure
symbol_names = plenum_pressure
symbol_values = plenum_pressure
[]
[]
[Physics/SolidMechanics/Layered1D]
[fuel]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = fuel_axial_pressure
generate_output = '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'
block = fuel
eigenstrain_names = 'fuel_thermal_strain gas_swelling_eigenstrain solid_swelling_eigenstrain'
mesh_generator = layered1D_mesh
[]
[clad]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = clad_axial_pressure
generate_output = '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'
block = clad
eigenstrain_names = 'clad_thermal_eigenstrain'
mesh_generator = layered1D_mesh
[]
[]
[Kernels]
# Define kernels for the various terms in the PDE system
[heat]
type = HeatConduction
variable = temp
[]
[heat_ie]
type = HeatConductionTimeDerivative
variable = temp
[]
[heat_source]
type = NeutronHeatSource
variable = temp
block = fuel
fission_rate = fission_rate
[]
[]
[AuxKernels]
[conductance]
type = MaterialRealAux
property = gap_conductance
variable = gap_cond
boundary = 10
[]
[gas_swell]
type = MaterialRealAux
variable = gas_swell
property = gas_swelling
execute_on = timestep_end
[]
[solid_swell]
type = MaterialRealAux
variable = solid_swell
property = solid_swelling
execute_on = timestep_end
[]
[volumetric_strain]
type = RankTwoScalarAux
rank_two_tensor = total_strain
variable = volumetric_strain
scalar_type = VolumetricStrain
execute_on = timestep_end
block = fuel
[]
[hoop_stress]
type = RankTwoAux
rank_two_tensor = stress
variable = hoop_stress
index_j = 2
index_i = 2
execute_on = timestep_end
[]
[hoop_creep_strain]
type = RankTwoAux
rank_two_tensor = creep_strain
variable = hoop_creep_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[hoop_elastic_strain]
type = RankTwoAux
rank_two_tensor = elastic_strain
variable = hoop_elastic_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[total_hoop_strain]
type = RankTwoAux
rank_two_tensor = total_strain
variable = total_hoop_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[]
[Contact]
[pellet_clad_mechanical]
primary = 5
secondary = 10
penalty = 1e12
model = frictionless
formulation = kinematic
[]
[]
[ThermalContact]
[thermal_contact]
type = GapHeatTransfer
variable = temp
primary = 5
secondary = 10
quadrature = true
gap_conductivity = 61.0
min_gap = ${clad_gap_width}
[]
[]
[BCs]
[no_x_all]
type = DirichletBC
variable = disp_x
boundary = 12
value = 0.0
[]
[Pressure]
[coolantPressure]
boundary = '2'
function = coolant_press_ramp
[]
[]
[PlenumPressure]
[plenumPressure]
boundary = 9
initial_pressure = 0.084e6 # Pa
startup_time = 0
R = 8.3143
temperature = ave_temp_interior
volume = gas_volume
output = plenum_pressure
material_input = fis_gas_released
[]
[]
[]
[CoolantChannel]
[convective_clad_surface]
boundary = '2'
variable = temp
inlet_temperature = coolant_temp_ramp
inlet_pressure = coolant_press_ramp
inlet_massflux = 5261.5 # kg/m^2-sec
coolant_material = sodium
rod_diameter = 5.84e-3 # m
rod_pitch = 7.48e-3 # m (Pitch-to-diameter Ratio = 1.28)
linear_heat_rate = power_history
axial_power_profile = axial_peaking_factors
subchannel_geometry = triangular
[]
[]
[Materials]
[fission_rate]
type = UPuZrFissionRate
rod_linear_power = power_history
axial_power_profile = axial_peaking_factors
pellet_radius = ${pellet_outer_radius}
X_Zr = ${X_Zr}
X_Pu_function = ${X_Pu}
block = fuel
outputs = all
[]
[burnup]
type = UPuZrBurnup
initial_X_Zr = ${X_Zr}
initial_X_Pu = ${X_Pu}
density = ${initial_fuel_density}
block = fuel
outputs = all
[]
[fuel_elasticity_tensor]
type = UPuZrElasticityTensor
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
block = fuel
temperature = temp
[]
[fuel_inelastic_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'fuel_upuzrcreep'
block = fuel
[]
[fuel_upuzrcreep]
type = UPuZrCreepUpdate
block = fuel
temperature = temp
porosity = porosity
max_inelastic_increment = 1e-3
[]
[fuel_thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = fuel
thermal_expansion_coeff = 1.18e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = fuel_thermal_strain
[]
[gas_swelling]
type = UPuZrGaseousEigenstrain
temperature = temp
eigenstrain_name = gas_swelling_eigenstrain
initial_porosity = 0.0
bubble_number_density = 1e20
outputs = all
output_properties = 'porosity gaseous_porosity'
block = fuel
[]
[solid_swelling]
type = BurnupDependentEigenstrain
eigenstrain_name = solid_swelling_eigenstrain
block = fuel
swelling_name = 'solid_swelling'
[]
[metal_fuel_thermal]
type = UPuZrThermal
block = fuel
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
spheat_model = savage
thcond_model = lanl
porosity = porosity
temperature = temp
[]
[fuel_density]
type = StrainAdjustedDensity
block = fuel
strain_free_density = 15800
[]
[fission_gas_behavior]
type = UPuZrFissionGasRelease
block = fuel
fission_rate = fission_rate
[]
[clad_elasticity_tensor]
type = ComputeIsotropicElasticityTensor
youngs_modulus = 1.88e11
poissons_ratio = 0.236
block = clad
[]
[clad_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'clad_ht9creep'
block = clad
[]
[fast_flux]
type = FastNeutronFlux
block = clad
factor = 2.47e19
[]
[clad_ht9creep]
type = HT9CreepUpdate
block = clad
temperature = temp
[]
[thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = clad
thermal_expansion_coeff = 1.2e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = clad_thermal_eigenstrain
[]
[clad_thermal]
type = HT9Thermal
block = clad
temperature = temp
[]
[clad_density]
type = StrainAdjustedDensity
block = clad
strain_free_density = 7874.0
[]
[]
[Preconditioning]
[SMP]
type = SMP
full = true
[]
[]
[Executioner]
type = Transient
solve_type = 'PJFNK'
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -ksp_gmres_restart'
petsc_options_value = 'lu superlu_dist 51'
line_search = 'none'
l_max_its = 60
l_tol = 8e-3
nl_max_its = 40
nl_rel_tol = 5e-4
nl_abs_tol = 1e-7
end_time = 42000400
dtmin = 100
dtmax = 5e5
[Quadrature]
order = fifth
side_order = seventh
[]
[TimeStepper]
type = IterationAdaptiveDT
timestep_limiting_postprocessor = creep_timestep
dt = 1e2
time_t = '0 1e5 41990400 42000400'
time_dt = '1e2 1e2 1e2 1e2'
iteration_window = 4
optimal_iterations = 10
[]
[]
[Postprocessors]
[ave_temp_interior] # average temperature of the cladding interior and all pellet exteriors
type = LayeredSideAverageValuePostprocessor
boundary = 9
variable = temp
execute_on = 'initial linear'
fuel_pin_geometry = pin_geometry
[]
[clad_inner_vol] # volume inside of cladding
type = LayeredInternalVolumePostprocessor
boundary = 7
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[pellet_volume] # fuel pellet total volume
type = LayeredInternalVolumePostprocessor
boundary = 8
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[avg_clad_temp] # average temperature of cladding interior
type = LayeredSideAverageValuePostprocessor
boundary = 7
variable = temp
fuel_pin_geometry = pin_geometry
execute_on = 'initial linear'
[]
[gas_volume]
type = LayeredInternalVolumePostprocessor
boundary = 9
execute_on = 'initial linear'
component = 0
out_of_plane_strain = strain_yy
fuel_pin_geometry = pin_geometry
addition = ${gas_addition}
[]
[flux_from_clad] # area integrated heat flux from the cladding
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 5
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[flux_from_fuel] # area integrated heat flux from the fuel
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 10
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[rod_total_power]
type = LayeredElementIntegralPowerPostprocessor
variable = temp
fission_rate = fission_rate
block = fuel
fuel_pin_geometry = pin_geometry
[]
[approx_FCT]
type = AverageNodalVariableValue
boundary = 12
variable = temp
[]
[max_approx_FCT]
type = TimeExtremeValue
value_type = max
postprocessor = approx_FCT
[]
[ave_FST]
type = SideAverageValue
boundary = 10
variable = temp
[]
[max_ave_FST]
type = TimeExtremeValue
value_type = max
postprocessor = ave_FST
[]
[ave_CIT]
type = SideAverageValue
boundary = 5
variable = temp
[]
[max_ave_CIT]
type = TimeExtremeValue
value_type = max
postprocessor = ave_CIT
[]
[peak_clad_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = clad
[]
[peak_fuel_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = fuel
[]
[max_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = max
block = fuel
[]
[min_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = min
block = fuel
[]
[peak_porosity]
type = ElementExtremeValue
variable = porosity
value_type = max
block = fuel
[]
[clad_fuel_gap]
type = NodalExtremeValue
variable = penetration
boundary = 10
[]
[max_cont_press]
type = NodalExtremeValue
variable = contact_pressure
boundary = 10
[]
[LHGR_W_per_cm]
type = FunctionValuePostprocessor
function = power_history
scale_factor = 0.01
[]
[average_burnup]
type = ElementAverageValue
block = fuel
variable = burnup
[]
[fis_gas_produced]
type = LayeredElementIntegralFisGasProducePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
[]
[fis_gas_released]
type = LayeredElementIntegralFisGasReleasePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
execute_on = 'initial timestep_end'
[]
[creep_timestep]
type = MaterialTimeStepPostprocessor
block = fuel
[]
[hydrostatic_stress]
type = ElementAverageValue
variable = hydrostatic_stress
execute_on = 'initial timestep_end'
block = fuel
[]
[solid_swelling]
type = ElementAverageValue
variable = solid_swell
block = fuel
[]
[gas_swelling]
type = ElementAverageValue
variable = gas_swell
block = fuel
[]
[volumetric_strain]
type = ElementAverageValue
variable = volumetric_strain
block = fuel
[]
[fission_rate]
type = ElementAverageValue
variable = fission_rate
block = fuel
[]
[porosity]
type = ElementAverageValue
variable = porosity
block = fuel
[]
[gaseous_porosity]
type = ElementAverageValue
variable = gaseous_porosity
block = fuel
[]
[fis_gas_percent]
type = FGRPercent
fission_gas_released = fis_gas_released
fission_gas_generated = fis_gas_produced
[]
[max_clad_hoop_creep]
type = ElementExtremeValue
value_type = max
block = clad
variable = hoop_creep_strain
[]
[max_total_hoop_strain]
type = ElementExtremeValue
value_type = max
block = clad
variable = total_hoop_strain
[]
[]
[VectorPostprocessors]
[clad_x_disp]
type = NodalValueSampler
variable = disp_x
boundary = 2
sort_by = y
outputs = 'vec1'
[]
[fuel_cl_temp]
type = NodalValueSampler
variable = temp
boundary = 12
sort_by = y
outputs = 'vec2'
[]
[fuel_surf_temp]
type = NodalValueSampler
variable = temp
boundary = 10
sort_by = y
outputs = 'vec3'
[]
[clad_inn_temp]
type = NodalValueSampler
variable = temp
boundary = 5
sort_by = y
outputs = 'vec4'
[]
[clad_out_temp]
type = NodalValueSampler
variable = temp
boundary = 2
sort_by = y
outputs = 'vec5'
[]
[fuel_surf_disp_x]
type = NodalValueSampler
variable = disp_x
boundary = 10
sort_by = y
outputs = 'vec7'
[]
[]
[PerformanceMetricOutputs]
outputs = performance_metrics_file
[]
[Outputs]
time_step_interval = 10
color = true
exodus = true
perf_graph = true
csv = true
sync_times = '1e3 5e3 1e4 5e4 1e5 5e6 1e6 5e6 1e7 2e7 3e7 4e7 41990400 41991000'
file_base = x441_${group_name}_1_5D
[out2]
type = CSV
file_base = x441_${group_name}_1_5D_out2
time_step_interval = 1
[]
[console]
type = Console
max_rows = 25
time_step_interval = 1
output_linear = true
[]
[chkfile]
type = CSV
file_base = x441_${group_name}_1_5D_chkfile
show = 'max_approx_FCT max_ave_FST max_ave_CIT average_burnup fis_gas_percent max_clad_hoop_creep max_total_hoop_strain'
execute_on = 'FINAL'
[]
[performance_metrics_file]
type = CSV
file_base = x441_${group_name}_1_5D_performance_metrics
show = 'simulation_alive_time number_linear_iterations number_nonlinear_iterations time_step_size total_linear_iterations total_nonlinear_iterations physical_memory_use number_dofs number_nonlinear_variables residual_compute_time jacobian_compute_time'
[]
[vec1]
type = CSV
file_base = x441_${group_name}_1_5D_vec1
execute_on = 'FINAL'
[]
[vec2]
type = CSV
file_base = x441_${group_name}_1_5D_vec2
execute_on = 'FINAL'
[]
[vec3]
type = CSV
file_base = x441_${group_name}_1_5D_vec3
execute_on = 'FINAL'
[]
[vec4]
type = CSV
file_base = x441_${group_name}_1_5D_vec4
execute_on = 'FINAL'
[]
[vec5]
type = CSV
file_base = x441_${group_name}_1_5D_vec5
execute_on = 'FINAL'
[]
[vec7]
type = CSV
file_base = x441_${group_name}_1_5D_vec7
execute_on = 'FINAL'
[]
[]
[Debug]
show_var_residual = 'disp_x temp'
show_var_residual_norms = true
[]
(assessment/metallic_fuel/EBRII/X441/analysis/group_H/x441_1_5D_H.i)
initial_fuel_density = 15800.0
[GlobalParams]
order = FIRST
family = LAGRANGE
energy_per_fission = 3.2e-11 # J/fission
displacements = 'disp_x'
[]
[Mesh]
coord_type = RZ
# Nominal Design Geometric Parameters (X441)
[layered1D_mesh]
type = Layered1DMeshGenerator
clad_thickness = ${clad_thickness}
pellet_outer_radius = ${pellet_outer_radius}
fuel_height = ${pellet_height}
plenum_height = ${clad_top_gap_height}
clad_gap_width = ${clad_gap_width}
# meshing parameters
clad_mesh_density = customize
pellet_mesh_density = customize
nx_p = 6
nx_c = 4
slices_per_block = 10
elem_type = EDGE2
[]
# mesh options
patch_update_strategy = auto
partitioner = centroid
centroid_partitioner_direction = y
[]
[UserObjects]
[pin_geometry]
type = Layered1DFuelPinGeometry
mesh_generator = layered1D_mesh
[]
[]
[Variables]
[temp]
initial_condition = 298
[]
[]
[AuxVariables]
[disp_y]
[]
[disp_z]
[]
# Aux variables for output
[gap_cond]
order = CONSTANT
family = MONOMIAL
[]
[coolant_htc]
order = CONSTANT
family = MONOMIAL
[]
[cumulative_damage_index]
order = CONSTANT
family = MONOMIAL
[]
[solid_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[gas_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[volumetric_strain]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[hoop_stress]
order = CONSTANT
family = MONOMIAL
[]
[hoop_creep_strain]
order = CONSTANT
family = MONOMIAL
[]
[hoop_elastic_strain]
order = CONSTANT
family = MONOMIAL
[]
[total_hoop_strain]
order = CONSTANT
family = MONOMIAL
[]
[]
[Functions]
[power_history]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '0 44722 44722 0'
[]
[coolant_press_ramp]
type = PiecewiseLinear
x = '0 42000400'
y = '1.0 1.0'
[]
[coolant_temp_ramp]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '298.0 648.0 648.0 350.0'
[]
[axial_peaking_factors]
type = PowerPeakingFunction
fit = EBRII_ROW_4
pellet_length = ${pellet_height}
pellet_y_start = ${pellet_y_start}
[]
[clad_axial_pressure]
type = CladdingAxialPressureFunction
plenum_pressure = plenum_pressure
coolant_pressure = coolant_press_ramp
coolant_pressure_scaling_factor = 0.151e6
fuel_pin_geometry = pin_geometry
[]
[fuel_axial_pressure]
type = ParsedFunction
expression = plenum_pressure
symbol_names = plenum_pressure
symbol_values = plenum_pressure
[]
[]
[Physics/SolidMechanics/Layered1D]
[fuel]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = fuel_axial_pressure
generate_output = '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'
block = fuel
eigenstrain_names = 'fuel_thermal_strain gas_swelling_eigenstrain solid_swelling_eigenstrain'
mesh_generator = layered1D_mesh
[]
[clad]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = clad_axial_pressure
generate_output = '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'
block = clad
eigenstrain_names = 'clad_thermal_eigenstrain'
mesh_generator = layered1D_mesh
[]
[]
[Kernels]
# Define kernels for the various terms in the PDE system
[heat]
type = HeatConduction
variable = temp
[]
[heat_ie]
type = HeatConductionTimeDerivative
variable = temp
[]
[heat_source]
type = NeutronHeatSource
variable = temp
block = fuel
fission_rate = fission_rate
[]
[]
[AuxKernels]
[conductance]
type = MaterialRealAux
property = gap_conductance
variable = gap_cond
boundary = 10
[]
[gas_swell]
type = MaterialRealAux
variable = gas_swell
property = gas_swelling
execute_on = timestep_end
[]
[solid_swell]
type = MaterialRealAux
variable = solid_swell
property = solid_swelling
execute_on = timestep_end
[]
[volumetric_strain]
type = RankTwoScalarAux
rank_two_tensor = total_strain
variable = volumetric_strain
scalar_type = VolumetricStrain
execute_on = timestep_end
block = fuel
[]
[hoop_stress]
type = RankTwoAux
rank_two_tensor = stress
variable = hoop_stress
index_j = 2
index_i = 2
execute_on = timestep_end
[]
[hoop_creep_strain]
type = RankTwoAux
rank_two_tensor = creep_strain
variable = hoop_creep_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[hoop_elastic_strain]
type = RankTwoAux
rank_two_tensor = elastic_strain
variable = hoop_elastic_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[total_hoop_strain]
type = RankTwoAux
rank_two_tensor = total_strain
variable = total_hoop_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[]
[Contact]
[pellet_clad_mechanical]
primary = 5
secondary = 10
penalty = 1e12
model = frictionless
formulation = kinematic
[]
[]
[ThermalContact]
[thermal_contact]
type = GapHeatTransfer
variable = temp
primary = 5
secondary = 10
quadrature = true
gap_conductivity = 61.0
min_gap = ${clad_gap_width}
[]
[]
[BCs]
[no_x_all]
type = DirichletBC
variable = disp_x
boundary = 12
value = 0.0
[]
[Pressure]
[coolantPressure]
boundary = '2'
function = coolant_press_ramp
[]
[]
[PlenumPressure]
[plenumPressure]
boundary = 9
initial_pressure = 0.084e6 # Pa
startup_time = 0
R = 8.3143
temperature = ave_temp_interior
volume = gas_volume
output = plenum_pressure
material_input = fis_gas_released
[]
[]
[]
[CoolantChannel]
[convective_clad_surface]
boundary = '2'
variable = temp
inlet_temperature = coolant_temp_ramp
inlet_pressure = coolant_press_ramp
inlet_massflux = 5261.5 # kg/m^2-sec
coolant_material = sodium
rod_diameter = 5.84e-3 # m
rod_pitch = 7.48e-3 # m (Pitch-to-diameter Ratio = 1.28)
linear_heat_rate = power_history
axial_power_profile = axial_peaking_factors
subchannel_geometry = triangular
[]
[]
[Materials]
[fission_rate]
type = UPuZrFissionRate
rod_linear_power = power_history
axial_power_profile = axial_peaking_factors
pellet_radius = ${pellet_outer_radius}
X_Zr = ${X_Zr}
X_Pu_function = ${X_Pu}
block = fuel
outputs = all
[]
[burnup]
type = UPuZrBurnup
initial_X_Zr = ${X_Zr}
initial_X_Pu = ${X_Pu}
density = ${initial_fuel_density}
block = fuel
outputs = all
[]
[fuel_elasticity_tensor]
type = UPuZrElasticityTensor
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
block = fuel
temperature = temp
[]
[fuel_inelastic_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'fuel_upuzrcreep'
block = fuel
[]
[fuel_upuzrcreep]
type = UPuZrCreepUpdate
block = fuel
temperature = temp
porosity = porosity
max_inelastic_increment = 1e-3
[]
[fuel_thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = fuel
thermal_expansion_coeff = 1.18e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = fuel_thermal_strain
[]
[gas_swelling]
type = UPuZrGaseousEigenstrain
temperature = temp
eigenstrain_name = gas_swelling_eigenstrain
initial_porosity = 0.0
bubble_number_density = 1e20
outputs = all
output_properties = 'porosity gaseous_porosity'
block = fuel
[]
[solid_swelling]
type = BurnupDependentEigenstrain
eigenstrain_name = solid_swelling_eigenstrain
block = fuel
swelling_name = 'solid_swelling'
[]
[metal_fuel_thermal]
type = UPuZrThermal
block = fuel
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
spheat_model = savage
thcond_model = lanl
porosity = porosity
temperature = temp
[]
[fuel_density]
type = StrainAdjustedDensity
block = fuel
strain_free_density = 15800
[]
[fission_gas_behavior]
type = UPuZrFissionGasRelease
block = fuel
fission_rate = fission_rate
[]
[clad_elasticity_tensor]
type = ComputeIsotropicElasticityTensor
youngs_modulus = 1.88e11
poissons_ratio = 0.236
block = clad
[]
[clad_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'clad_ht9creep'
block = clad
[]
[fast_flux]
type = FastNeutronFlux
block = clad
factor = 2.47e19
[]
[clad_ht9creep]
type = HT9CreepUpdate
block = clad
temperature = temp
[]
[thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = clad
thermal_expansion_coeff = 1.2e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = clad_thermal_eigenstrain
[]
[clad_thermal]
type = HT9Thermal
block = clad
temperature = temp
[]
[clad_density]
type = StrainAdjustedDensity
block = clad
strain_free_density = 7874.0
[]
[]
[Preconditioning]
[SMP]
type = SMP
full = true
[]
[]
[Executioner]
type = Transient
solve_type = 'PJFNK'
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -ksp_gmres_restart'
petsc_options_value = 'lu superlu_dist 51'
line_search = 'none'
l_max_its = 60
l_tol = 8e-3
nl_max_its = 40
nl_rel_tol = 5e-4
nl_abs_tol = 1e-7
end_time = 42000400
dtmin = 100
dtmax = 5e5
[Quadrature]
order = fifth
side_order = seventh
[]
[TimeStepper]
type = IterationAdaptiveDT
timestep_limiting_postprocessor = creep_timestep
dt = 1e2
time_t = '0 1e5 41990400 42000400'
time_dt = '1e2 1e2 1e2 1e2'
iteration_window = 4
optimal_iterations = 10
[]
[]
[Postprocessors]
[ave_temp_interior] # average temperature of the cladding interior and all pellet exteriors
type = LayeredSideAverageValuePostprocessor
boundary = 9
variable = temp
execute_on = 'initial linear'
fuel_pin_geometry = pin_geometry
[]
[clad_inner_vol] # volume inside of cladding
type = LayeredInternalVolumePostprocessor
boundary = 7
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[pellet_volume] # fuel pellet total volume
type = LayeredInternalVolumePostprocessor
boundary = 8
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[avg_clad_temp] # average temperature of cladding interior
type = LayeredSideAverageValuePostprocessor
boundary = 7
variable = temp
fuel_pin_geometry = pin_geometry
execute_on = 'initial linear'
[]
[gas_volume]
type = LayeredInternalVolumePostprocessor
boundary = 9
execute_on = 'initial linear'
component = 0
out_of_plane_strain = strain_yy
fuel_pin_geometry = pin_geometry
addition = ${gas_addition}
[]
[flux_from_clad] # area integrated heat flux from the cladding
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 5
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[flux_from_fuel] # area integrated heat flux from the fuel
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 10
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[rod_total_power]
type = LayeredElementIntegralPowerPostprocessor
variable = temp
fission_rate = fission_rate
block = fuel
fuel_pin_geometry = pin_geometry
[]
[approx_FCT]
type = AverageNodalVariableValue
boundary = 12
variable = temp
[]
[max_approx_FCT]
type = TimeExtremeValue
value_type = max
postprocessor = approx_FCT
[]
[ave_FST]
type = SideAverageValue
boundary = 10
variable = temp
[]
[max_ave_FST]
type = TimeExtremeValue
value_type = max
postprocessor = ave_FST
[]
[ave_CIT]
type = SideAverageValue
boundary = 5
variable = temp
[]
[max_ave_CIT]
type = TimeExtremeValue
value_type = max
postprocessor = ave_CIT
[]
[peak_clad_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = clad
[]
[peak_fuel_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = fuel
[]
[max_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = max
block = fuel
[]
[min_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = min
block = fuel
[]
[peak_porosity]
type = ElementExtremeValue
variable = porosity
value_type = max
block = fuel
[]
[clad_fuel_gap]
type = NodalExtremeValue
variable = penetration
boundary = 10
[]
[max_cont_press]
type = NodalExtremeValue
variable = contact_pressure
boundary = 10
[]
[LHGR_W_per_cm]
type = FunctionValuePostprocessor
function = power_history
scale_factor = 0.01
[]
[average_burnup]
type = ElementAverageValue
block = fuel
variable = burnup
[]
[fis_gas_produced]
type = LayeredElementIntegralFisGasProducePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
[]
[fis_gas_released]
type = LayeredElementIntegralFisGasReleasePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
execute_on = 'initial timestep_end'
[]
[creep_timestep]
type = MaterialTimeStepPostprocessor
block = fuel
[]
[hydrostatic_stress]
type = ElementAverageValue
variable = hydrostatic_stress
execute_on = 'initial timestep_end'
block = fuel
[]
[solid_swelling]
type = ElementAverageValue
variable = solid_swell
block = fuel
[]
[gas_swelling]
type = ElementAverageValue
variable = gas_swell
block = fuel
[]
[volumetric_strain]
type = ElementAverageValue
variable = volumetric_strain
block = fuel
[]
[fission_rate]
type = ElementAverageValue
variable = fission_rate
block = fuel
[]
[porosity]
type = ElementAverageValue
variable = porosity
block = fuel
[]
[gaseous_porosity]
type = ElementAverageValue
variable = gaseous_porosity
block = fuel
[]
[fis_gas_percent]
type = FGRPercent
fission_gas_released = fis_gas_released
fission_gas_generated = fis_gas_produced
[]
[max_clad_hoop_creep]
type = ElementExtremeValue
value_type = max
block = clad
variable = hoop_creep_strain
[]
[max_total_hoop_strain]
type = ElementExtremeValue
value_type = max
block = clad
variable = total_hoop_strain
[]
[]
[VectorPostprocessors]
[clad_x_disp]
type = NodalValueSampler
variable = disp_x
boundary = 2
sort_by = y
outputs = 'vec1'
[]
[fuel_cl_temp]
type = NodalValueSampler
variable = temp
boundary = 12
sort_by = y
outputs = 'vec2'
[]
[fuel_surf_temp]
type = NodalValueSampler
variable = temp
boundary = 10
sort_by = y
outputs = 'vec3'
[]
[clad_inn_temp]
type = NodalValueSampler
variable = temp
boundary = 5
sort_by = y
outputs = 'vec4'
[]
[clad_out_temp]
type = NodalValueSampler
variable = temp
boundary = 2
sort_by = y
outputs = 'vec5'
[]
[fuel_surf_disp_x]
type = NodalValueSampler
variable = disp_x
boundary = 10
sort_by = y
outputs = 'vec7'
[]
[]
[PerformanceMetricOutputs]
outputs = performance_metrics_file
[]
[Outputs]
time_step_interval = 10
color = true
exodus = true
perf_graph = true
csv = true
sync_times = '1e3 5e3 1e4 5e4 1e5 5e6 1e6 5e6 1e7 2e7 3e7 4e7 41990400 41991000'
file_base = x441_${group_name}_1_5D
[out2]
type = CSV
file_base = x441_${group_name}_1_5D_out2
time_step_interval = 1
[]
[console]
type = Console
max_rows = 25
time_step_interval = 1
output_linear = true
[]
[chkfile]
type = CSV
file_base = x441_${group_name}_1_5D_chkfile
show = 'max_approx_FCT max_ave_FST max_ave_CIT average_burnup fis_gas_percent max_clad_hoop_creep max_total_hoop_strain'
execute_on = 'FINAL'
[]
[performance_metrics_file]
type = CSV
file_base = x441_${group_name}_1_5D_performance_metrics
show = 'simulation_alive_time number_linear_iterations number_nonlinear_iterations time_step_size total_linear_iterations total_nonlinear_iterations physical_memory_use number_dofs number_nonlinear_variables residual_compute_time jacobian_compute_time'
[]
[vec1]
type = CSV
file_base = x441_${group_name}_1_5D_vec1
execute_on = 'FINAL'
[]
[vec2]
type = CSV
file_base = x441_${group_name}_1_5D_vec2
execute_on = 'FINAL'
[]
[vec3]
type = CSV
file_base = x441_${group_name}_1_5D_vec3
execute_on = 'FINAL'
[]
[vec4]
type = CSV
file_base = x441_${group_name}_1_5D_vec4
execute_on = 'FINAL'
[]
[vec5]
type = CSV
file_base = x441_${group_name}_1_5D_vec5
execute_on = 'FINAL'
[]
[vec7]
type = CSV
file_base = x441_${group_name}_1_5D_vec7
execute_on = 'FINAL'
[]
[]
[Debug]
show_var_residual = 'disp_x temp'
show_var_residual_norms = true
[]
(test/tests/upuzr_fission_gas_release/nonad/exact_1D.i)
# This tests the fission gas produced and released postprocessors
# and for 1d and 2d models with the same volume.
# The reported values should be the same for the
# 1d and 2d cases. See corresponding excel calculation
# in this directory.
[Problem]
solve = false
[]
[UserObjects]
[pin_geometry]
type = Layered1DFuelPinGeometry
include_clad = false
mesh_generator = layered1D_mesh
[]
[]
[Mesh]
coord_type = RZ
[layered1D_mesh]
type = Layered1DMeshGenerator
include_clad = false
include_plenum = false
slices_per_block = 1
fuel_height = 3e-3
pellet_outer_radius = 3e-3
pellet_bottom_coor = 0.0
pellet_mesh_density = customize
nx_p = 4
elem_type = EDGE2
[]
[]
[Functions]
[fission_rate_func]
type = PiecewiseLinear
x = '0 1e7'
y = '1e20 1e18'
[]
[porosity_func]
type = PiecewiseLinear
x = '0 1e7'
y = '0.15 0.25'
[]
[]
[Variables]
[temp]
initial_condition = 1700
[]
[]
[Kernels]
[heat]
type = TimeDerivative
variable = temp
[]
[]
[Materials]
[properties]
type = GenericFunctionMaterial
prop_names = 'fission_rate porosity'
prop_values = 'fission_rate_func porosity_func'
outputs = all
[]
[fission_gas_release]
type = UPuZrFissionGasRelease
fractional_yield = 0.3017
critical_porosity = 0.2
use_old_porosity = false
time_average_fission_rate = false
[]
[]
[Executioner]
type = Transient
dt = 1e6
end_time = 1e7
[]
[Postprocessors]
[fission_rate]
type = ElementAverageValue
variable = fission_rate
[]
[gas_produced]
type = LayeredElementIntegralFisGasProducePostprocessor
fuel_pin_geometry = pin_geometry
[]
[gas_released]
type = LayeredElementIntegralFisGasReleasePostprocessor
fuel_pin_geometry = pin_geometry
[]
[porosity]
type = ElementAverageValue
variable = porosity
[]
[]
[Outputs]
csv = true
file_base = exact_out
[]
(assessment/metallic_fuel/EBRII/X441/analysis/group_G/x441_1_5D_G.i)
initial_fuel_density = 15800.0
[GlobalParams]
order = FIRST
family = LAGRANGE
energy_per_fission = 3.2e-11 # J/fission
displacements = 'disp_x'
[]
[Mesh]
coord_type = RZ
# Nominal Design Geometric Parameters (X441)
[layered1D_mesh]
type = Layered1DMeshGenerator
clad_thickness = ${clad_thickness}
pellet_outer_radius = ${pellet_outer_radius}
fuel_height = ${pellet_height}
plenum_height = ${clad_top_gap_height}
clad_gap_width = ${clad_gap_width}
# meshing parameters
clad_mesh_density = customize
pellet_mesh_density = customize
nx_p = 6
nx_c = 4
slices_per_block = 10
elem_type = EDGE2
[]
# mesh options
patch_update_strategy = auto
partitioner = centroid
centroid_partitioner_direction = y
[]
[UserObjects]
[pin_geometry]
type = Layered1DFuelPinGeometry
mesh_generator = layered1D_mesh
[]
[]
[Variables]
[temp]
initial_condition = 298
[]
[]
[AuxVariables]
[disp_y]
[]
[disp_z]
[]
# Aux variables for output
[gap_cond]
order = CONSTANT
family = MONOMIAL
[]
[coolant_htc]
order = CONSTANT
family = MONOMIAL
[]
[cumulative_damage_index]
order = CONSTANT
family = MONOMIAL
[]
[solid_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[gas_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[volumetric_strain]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[hoop_stress]
order = CONSTANT
family = MONOMIAL
[]
[hoop_creep_strain]
order = CONSTANT
family = MONOMIAL
[]
[hoop_elastic_strain]
order = CONSTANT
family = MONOMIAL
[]
[total_hoop_strain]
order = CONSTANT
family = MONOMIAL
[]
[]
[Functions]
[power_history]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '0 44722 44722 0'
[]
[coolant_press_ramp]
type = PiecewiseLinear
x = '0 42000400'
y = '1.0 1.0'
[]
[coolant_temp_ramp]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '298.0 648.0 648.0 350.0'
[]
[axial_peaking_factors]
type = PowerPeakingFunction
fit = EBRII_ROW_4
pellet_length = ${pellet_height}
pellet_y_start = ${pellet_y_start}
[]
[clad_axial_pressure]
type = CladdingAxialPressureFunction
plenum_pressure = plenum_pressure
coolant_pressure = coolant_press_ramp
coolant_pressure_scaling_factor = 0.151e6
fuel_pin_geometry = pin_geometry
[]
[fuel_axial_pressure]
type = ParsedFunction
expression = plenum_pressure
symbol_names = plenum_pressure
symbol_values = plenum_pressure
[]
[]
[Physics/SolidMechanics/Layered1D]
[fuel]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = fuel_axial_pressure
generate_output = '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'
block = fuel
eigenstrain_names = 'fuel_thermal_strain gas_swelling_eigenstrain solid_swelling_eigenstrain'
mesh_generator = layered1D_mesh
[]
[clad]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = clad_axial_pressure
generate_output = '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'
block = clad
eigenstrain_names = 'clad_thermal_eigenstrain'
mesh_generator = layered1D_mesh
[]
[]
[Kernels]
# Define kernels for the various terms in the PDE system
[heat]
type = HeatConduction
variable = temp
[]
[heat_ie]
type = HeatConductionTimeDerivative
variable = temp
[]
[heat_source]
type = NeutronHeatSource
variable = temp
block = fuel
fission_rate = fission_rate
[]
[]
[AuxKernels]
[conductance]
type = MaterialRealAux
property = gap_conductance
variable = gap_cond
boundary = 10
[]
[gas_swell]
type = MaterialRealAux
variable = gas_swell
property = gas_swelling
execute_on = timestep_end
[]
[solid_swell]
type = MaterialRealAux
variable = solid_swell
property = solid_swelling
execute_on = timestep_end
[]
[volumetric_strain]
type = RankTwoScalarAux
rank_two_tensor = total_strain
variable = volumetric_strain
scalar_type = VolumetricStrain
execute_on = timestep_end
block = fuel
[]
[hoop_stress]
type = RankTwoAux
rank_two_tensor = stress
variable = hoop_stress
index_j = 2
index_i = 2
execute_on = timestep_end
[]
[hoop_creep_strain]
type = RankTwoAux
rank_two_tensor = creep_strain
variable = hoop_creep_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[hoop_elastic_strain]
type = RankTwoAux
rank_two_tensor = elastic_strain
variable = hoop_elastic_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[total_hoop_strain]
type = RankTwoAux
rank_two_tensor = total_strain
variable = total_hoop_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[]
[Contact]
[pellet_clad_mechanical]
primary = 5
secondary = 10
penalty = 1e12
model = frictionless
formulation = kinematic
[]
[]
[ThermalContact]
[thermal_contact]
type = GapHeatTransfer
variable = temp
primary = 5
secondary = 10
quadrature = true
gap_conductivity = 61.0
min_gap = ${clad_gap_width}
[]
[]
[BCs]
[no_x_all]
type = DirichletBC
variable = disp_x
boundary = 12
value = 0.0
[]
[Pressure]
[coolantPressure]
boundary = '2'
function = coolant_press_ramp
[]
[]
[PlenumPressure]
[plenumPressure]
boundary = 9
initial_pressure = 0.084e6 # Pa
startup_time = 0
R = 8.3143
temperature = ave_temp_interior
volume = gas_volume
output = plenum_pressure
material_input = fis_gas_released
[]
[]
[]
[CoolantChannel]
[convective_clad_surface]
boundary = '2'
variable = temp
inlet_temperature = coolant_temp_ramp
inlet_pressure = coolant_press_ramp
inlet_massflux = 5261.5 # kg/m^2-sec
coolant_material = sodium
rod_diameter = 5.84e-3 # m
rod_pitch = 7.48e-3 # m (Pitch-to-diameter Ratio = 1.28)
linear_heat_rate = power_history
axial_power_profile = axial_peaking_factors
subchannel_geometry = triangular
[]
[]
[Materials]
[fission_rate]
type = UPuZrFissionRate
rod_linear_power = power_history
axial_power_profile = axial_peaking_factors
pellet_radius = ${pellet_outer_radius}
X_Zr = ${X_Zr}
X_Pu_function = ${X_Pu}
block = fuel
outputs = all
[]
[burnup]
type = UPuZrBurnup
initial_X_Zr = ${X_Zr}
initial_X_Pu = ${X_Pu}
density = ${initial_fuel_density}
block = fuel
outputs = all
[]
[fuel_elasticity_tensor]
type = UPuZrElasticityTensor
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
block = fuel
temperature = temp
[]
[fuel_inelastic_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'fuel_upuzrcreep'
block = fuel
[]
[fuel_upuzrcreep]
type = UPuZrCreepUpdate
block = fuel
temperature = temp
porosity = porosity
max_inelastic_increment = 1e-3
[]
[fuel_thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = fuel
thermal_expansion_coeff = 1.18e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = fuel_thermal_strain
[]
[gas_swelling]
type = UPuZrGaseousEigenstrain
temperature = temp
eigenstrain_name = gas_swelling_eigenstrain
initial_porosity = 0.0
bubble_number_density = 1e20
outputs = all
output_properties = 'porosity gaseous_porosity'
block = fuel
[]
[solid_swelling]
type = BurnupDependentEigenstrain
eigenstrain_name = solid_swelling_eigenstrain
block = fuel
swelling_name = 'solid_swelling'
[]
[metal_fuel_thermal]
type = UPuZrThermal
block = fuel
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
spheat_model = savage
thcond_model = lanl
porosity = porosity
temperature = temp
[]
[fuel_density]
type = StrainAdjustedDensity
block = fuel
strain_free_density = 15800
[]
[fission_gas_behavior]
type = UPuZrFissionGasRelease
block = fuel
fission_rate = fission_rate
[]
[clad_elasticity_tensor]
type = ComputeIsotropicElasticityTensor
youngs_modulus = 1.88e11
poissons_ratio = 0.236
block = clad
[]
[clad_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'clad_ht9creep'
block = clad
[]
[fast_flux]
type = FastNeutronFlux
block = clad
factor = 2.47e19
[]
[clad_ht9creep]
type = HT9CreepUpdate
block = clad
temperature = temp
[]
[thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = clad
thermal_expansion_coeff = 1.2e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = clad_thermal_eigenstrain
[]
[clad_thermal]
type = HT9Thermal
block = clad
temperature = temp
[]
[clad_density]
type = StrainAdjustedDensity
block = clad
strain_free_density = 7874.0
[]
[]
[Preconditioning]
[SMP]
type = SMP
full = true
[]
[]
[Executioner]
type = Transient
solve_type = 'PJFNK'
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -ksp_gmres_restart'
petsc_options_value = 'lu superlu_dist 51'
line_search = 'none'
l_max_its = 60
l_tol = 8e-3
nl_max_its = 40
nl_rel_tol = 5e-4
nl_abs_tol = 1e-7
end_time = 42000400
dtmin = 100
dtmax = 5e5
[Quadrature]
order = fifth
side_order = seventh
[]
[TimeStepper]
type = IterationAdaptiveDT
timestep_limiting_postprocessor = creep_timestep
dt = 1e2
time_t = '0 1e5 41990400 42000400'
time_dt = '1e2 1e2 1e2 1e2'
iteration_window = 4
optimal_iterations = 10
[]
[]
[Postprocessors]
[ave_temp_interior] # average temperature of the cladding interior and all pellet exteriors
type = LayeredSideAverageValuePostprocessor
boundary = 9
variable = temp
execute_on = 'initial linear'
fuel_pin_geometry = pin_geometry
[]
[clad_inner_vol] # volume inside of cladding
type = LayeredInternalVolumePostprocessor
boundary = 7
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[pellet_volume] # fuel pellet total volume
type = LayeredInternalVolumePostprocessor
boundary = 8
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[avg_clad_temp] # average temperature of cladding interior
type = LayeredSideAverageValuePostprocessor
boundary = 7
variable = temp
fuel_pin_geometry = pin_geometry
execute_on = 'initial linear'
[]
[gas_volume]
type = LayeredInternalVolumePostprocessor
boundary = 9
execute_on = 'initial linear'
component = 0
out_of_plane_strain = strain_yy
fuel_pin_geometry = pin_geometry
addition = ${gas_addition}
[]
[flux_from_clad] # area integrated heat flux from the cladding
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 5
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[flux_from_fuel] # area integrated heat flux from the fuel
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 10
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[rod_total_power]
type = LayeredElementIntegralPowerPostprocessor
variable = temp
fission_rate = fission_rate
block = fuel
fuel_pin_geometry = pin_geometry
[]
[approx_FCT]
type = AverageNodalVariableValue
boundary = 12
variable = temp
[]
[max_approx_FCT]
type = TimeExtremeValue
value_type = max
postprocessor = approx_FCT
[]
[ave_FST]
type = SideAverageValue
boundary = 10
variable = temp
[]
[max_ave_FST]
type = TimeExtremeValue
value_type = max
postprocessor = ave_FST
[]
[ave_CIT]
type = SideAverageValue
boundary = 5
variable = temp
[]
[max_ave_CIT]
type = TimeExtremeValue
value_type = max
postprocessor = ave_CIT
[]
[peak_clad_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = clad
[]
[peak_fuel_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = fuel
[]
[max_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = max
block = fuel
[]
[min_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = min
block = fuel
[]
[peak_porosity]
type = ElementExtremeValue
variable = porosity
value_type = max
block = fuel
[]
[clad_fuel_gap]
type = NodalExtremeValue
variable = penetration
boundary = 10
[]
[max_cont_press]
type = NodalExtremeValue
variable = contact_pressure
boundary = 10
[]
[LHGR_W_per_cm]
type = FunctionValuePostprocessor
function = power_history
scale_factor = 0.01
[]
[average_burnup]
type = ElementAverageValue
block = fuel
variable = burnup
[]
[fis_gas_produced]
type = LayeredElementIntegralFisGasProducePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
[]
[fis_gas_released]
type = LayeredElementIntegralFisGasReleasePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
execute_on = 'initial timestep_end'
[]
[creep_timestep]
type = MaterialTimeStepPostprocessor
block = fuel
[]
[hydrostatic_stress]
type = ElementAverageValue
variable = hydrostatic_stress
execute_on = 'initial timestep_end'
block = fuel
[]
[solid_swelling]
type = ElementAverageValue
variable = solid_swell
block = fuel
[]
[gas_swelling]
type = ElementAverageValue
variable = gas_swell
block = fuel
[]
[volumetric_strain]
type = ElementAverageValue
variable = volumetric_strain
block = fuel
[]
[fission_rate]
type = ElementAverageValue
variable = fission_rate
block = fuel
[]
[porosity]
type = ElementAverageValue
variable = porosity
block = fuel
[]
[gaseous_porosity]
type = ElementAverageValue
variable = gaseous_porosity
block = fuel
[]
[fis_gas_percent]
type = FGRPercent
fission_gas_released = fis_gas_released
fission_gas_generated = fis_gas_produced
[]
[max_clad_hoop_creep]
type = ElementExtremeValue
value_type = max
block = clad
variable = hoop_creep_strain
[]
[max_total_hoop_strain]
type = ElementExtremeValue
value_type = max
block = clad
variable = total_hoop_strain
[]
[]
[VectorPostprocessors]
[clad_x_disp]
type = NodalValueSampler
variable = disp_x
boundary = 2
sort_by = y
outputs = 'vec1'
[]
[fuel_cl_temp]
type = NodalValueSampler
variable = temp
boundary = 12
sort_by = y
outputs = 'vec2'
[]
[fuel_surf_temp]
type = NodalValueSampler
variable = temp
boundary = 10
sort_by = y
outputs = 'vec3'
[]
[clad_inn_temp]
type = NodalValueSampler
variable = temp
boundary = 5
sort_by = y
outputs = 'vec4'
[]
[clad_out_temp]
type = NodalValueSampler
variable = temp
boundary = 2
sort_by = y
outputs = 'vec5'
[]
[fuel_surf_disp_x]
type = NodalValueSampler
variable = disp_x
boundary = 10
sort_by = y
outputs = 'vec7'
[]
[]
[PerformanceMetricOutputs]
outputs = performance_metrics_file
[]
[Outputs]
time_step_interval = 10
color = true
exodus = true
perf_graph = true
csv = true
sync_times = '1e3 5e3 1e4 5e4 1e5 5e6 1e6 5e6 1e7 2e7 3e7 4e7 41990400 41991000'
file_base = x441_${group_name}_1_5D
[out2]
type = CSV
file_base = x441_${group_name}_1_5D_out2
time_step_interval = 1
[]
[console]
type = Console
max_rows = 25
time_step_interval = 1
output_linear = true
[]
[chkfile]
type = CSV
file_base = x441_${group_name}_1_5D_chkfile
show = 'max_approx_FCT max_ave_FST max_ave_CIT average_burnup fis_gas_percent max_clad_hoop_creep max_total_hoop_strain'
execute_on = 'FINAL'
[]
[performance_metrics_file]
type = CSV
file_base = x441_${group_name}_1_5D_performance_metrics
show = 'simulation_alive_time number_linear_iterations number_nonlinear_iterations time_step_size total_linear_iterations total_nonlinear_iterations physical_memory_use number_dofs number_nonlinear_variables residual_compute_time jacobian_compute_time'
[]
[vec1]
type = CSV
file_base = x441_${group_name}_1_5D_vec1
execute_on = 'FINAL'
[]
[vec2]
type = CSV
file_base = x441_${group_name}_1_5D_vec2
execute_on = 'FINAL'
[]
[vec3]
type = CSV
file_base = x441_${group_name}_1_5D_vec3
execute_on = 'FINAL'
[]
[vec4]
type = CSV
file_base = x441_${group_name}_1_5D_vec4
execute_on = 'FINAL'
[]
[vec5]
type = CSV
file_base = x441_${group_name}_1_5D_vec5
execute_on = 'FINAL'
[]
[vec7]
type = CSV
file_base = x441_${group_name}_1_5D_vec7
execute_on = 'FINAL'
[]
[]
[Debug]
show_var_residual = 'disp_x temp'
show_var_residual_norms = true
[]
(assessment/metallic_fuel/EBRII/X441/analysis/group_D/x441_1_5D_D.i)
initial_fuel_density = 15800.0
[GlobalParams]
order = FIRST
family = LAGRANGE
energy_per_fission = 3.2e-11 # J/fission
displacements = 'disp_x'
[]
[Mesh]
coord_type = RZ
# Nominal Design Geometric Parameters (X441)
[layered1D_mesh]
type = Layered1DMeshGenerator
clad_thickness = ${clad_thickness}
pellet_outer_radius = ${pellet_outer_radius}
fuel_height = ${pellet_height}
plenum_height = ${clad_top_gap_height}
clad_gap_width = ${clad_gap_width}
# meshing parameters
clad_mesh_density = customize
pellet_mesh_density = customize
nx_p = 6
nx_c = 4
slices_per_block = 10
elem_type = EDGE2
[]
# mesh options
patch_update_strategy = auto
partitioner = centroid
centroid_partitioner_direction = y
[]
[UserObjects]
[pin_geometry]
type = Layered1DFuelPinGeometry
mesh_generator = layered1D_mesh
[]
[]
[Variables]
[temp]
initial_condition = 298
[]
[]
[AuxVariables]
[disp_y]
[]
[disp_z]
[]
# Aux variables for output
[gap_cond]
order = CONSTANT
family = MONOMIAL
[]
[coolant_htc]
order = CONSTANT
family = MONOMIAL
[]
[cumulative_damage_index]
order = CONSTANT
family = MONOMIAL
[]
[solid_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[gas_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[volumetric_strain]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[hoop_stress]
order = CONSTANT
family = MONOMIAL
[]
[hoop_creep_strain]
order = CONSTANT
family = MONOMIAL
[]
[hoop_elastic_strain]
order = CONSTANT
family = MONOMIAL
[]
[total_hoop_strain]
order = CONSTANT
family = MONOMIAL
[]
[]
[Functions]
[power_history]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '0 44722 44722 0'
[]
[coolant_press_ramp]
type = PiecewiseLinear
x = '0 42000400'
y = '1.0 1.0'
[]
[coolant_temp_ramp]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '298.0 648.0 648.0 350.0'
[]
[axial_peaking_factors]
type = PowerPeakingFunction
fit = EBRII_ROW_4
pellet_length = ${pellet_height}
pellet_y_start = ${pellet_y_start}
[]
[clad_axial_pressure]
type = CladdingAxialPressureFunction
plenum_pressure = plenum_pressure
coolant_pressure = coolant_press_ramp
coolant_pressure_scaling_factor = 0.151e6
fuel_pin_geometry = pin_geometry
[]
[fuel_axial_pressure]
type = ParsedFunction
expression = plenum_pressure
symbol_names = plenum_pressure
symbol_values = plenum_pressure
[]
[]
[Physics/SolidMechanics/Layered1D]
[fuel]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = fuel_axial_pressure
generate_output = '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'
block = fuel
eigenstrain_names = 'fuel_thermal_strain gas_swelling_eigenstrain solid_swelling_eigenstrain'
mesh_generator = layered1D_mesh
[]
[clad]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = clad_axial_pressure
generate_output = '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'
block = clad
eigenstrain_names = 'clad_thermal_eigenstrain'
mesh_generator = layered1D_mesh
[]
[]
[Kernels]
# Define kernels for the various terms in the PDE system
[heat]
type = HeatConduction
variable = temp
[]
[heat_ie]
type = HeatConductionTimeDerivative
variable = temp
[]
[heat_source]
type = NeutronHeatSource
variable = temp
block = fuel
fission_rate = fission_rate
[]
[]
[AuxKernels]
[conductance]
type = MaterialRealAux
property = gap_conductance
variable = gap_cond
boundary = 10
[]
[gas_swell]
type = MaterialRealAux
variable = gas_swell
property = gas_swelling
execute_on = timestep_end
[]
[solid_swell]
type = MaterialRealAux
variable = solid_swell
property = solid_swelling
execute_on = timestep_end
[]
[volumetric_strain]
type = RankTwoScalarAux
rank_two_tensor = total_strain
variable = volumetric_strain
scalar_type = VolumetricStrain
execute_on = timestep_end
block = fuel
[]
[hoop_stress]
type = RankTwoAux
rank_two_tensor = stress
variable = hoop_stress
index_j = 2
index_i = 2
execute_on = timestep_end
[]
[hoop_creep_strain]
type = RankTwoAux
rank_two_tensor = creep_strain
variable = hoop_creep_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[hoop_elastic_strain]
type = RankTwoAux
rank_two_tensor = elastic_strain
variable = hoop_elastic_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[total_hoop_strain]
type = RankTwoAux
rank_two_tensor = total_strain
variable = total_hoop_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[]
[Contact]
[pellet_clad_mechanical]
primary = 5
secondary = 10
penalty = 1e12
model = frictionless
formulation = kinematic
[]
[]
[ThermalContact]
[thermal_contact]
type = GapHeatTransfer
variable = temp
primary = 5
secondary = 10
quadrature = true
gap_conductivity = 61.0
min_gap = ${clad_gap_width}
[]
[]
[BCs]
[no_x_all]
type = DirichletBC
variable = disp_x
boundary = 12
value = 0.0
[]
[Pressure]
[coolantPressure]
boundary = '2'
function = coolant_press_ramp
[]
[]
[PlenumPressure]
[plenumPressure]
boundary = 9
initial_pressure = 0.084e6 # Pa
startup_time = 0
R = 8.3143
temperature = ave_temp_interior
volume = gas_volume
output = plenum_pressure
material_input = fis_gas_released
[]
[]
[]
[CoolantChannel]
[convective_clad_surface]
boundary = '2'
variable = temp
inlet_temperature = coolant_temp_ramp
inlet_pressure = coolant_press_ramp
inlet_massflux = 5261.5 # kg/m^2-sec
coolant_material = sodium
rod_diameter = 5.84e-3 # m
rod_pitch = 7.48e-3 # m (Pitch-to-diameter Ratio = 1.28)
linear_heat_rate = power_history
axial_power_profile = axial_peaking_factors
subchannel_geometry = triangular
[]
[]
[Materials]
[fission_rate]
type = UPuZrFissionRate
rod_linear_power = power_history
axial_power_profile = axial_peaking_factors
pellet_radius = ${pellet_outer_radius}
X_Zr = ${X_Zr}
X_Pu_function = ${X_Pu}
block = fuel
outputs = all
[]
[burnup]
type = UPuZrBurnup
initial_X_Zr = ${X_Zr}
initial_X_Pu = ${X_Pu}
density = ${initial_fuel_density}
block = fuel
outputs = all
[]
[fuel_elasticity_tensor]
type = UPuZrElasticityTensor
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
block = fuel
temperature = temp
[]
[fuel_inelastic_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'fuel_upuzrcreep'
block = fuel
[]
[fuel_upuzrcreep]
type = UPuZrCreepUpdate
block = fuel
temperature = temp
porosity = porosity
max_inelastic_increment = 1e-3
[]
[fuel_thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = fuel
thermal_expansion_coeff = 1.18e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = fuel_thermal_strain
[]
[gas_swelling]
type = UPuZrGaseousEigenstrain
temperature = temp
eigenstrain_name = gas_swelling_eigenstrain
initial_porosity = 0.0
bubble_number_density = 1e20
outputs = all
output_properties = 'porosity gaseous_porosity'
block = fuel
[]
[solid_swelling]
type = BurnupDependentEigenstrain
eigenstrain_name = solid_swelling_eigenstrain
block = fuel
swelling_name = 'solid_swelling'
[]
[metal_fuel_thermal]
type = UPuZrThermal
block = fuel
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
spheat_model = savage
thcond_model = lanl
porosity = porosity
temperature = temp
[]
[fuel_density]
type = StrainAdjustedDensity
block = fuel
strain_free_density = 15800
[]
[fission_gas_behavior]
type = UPuZrFissionGasRelease
block = fuel
fission_rate = fission_rate
[]
[clad_elasticity_tensor]
type = ComputeIsotropicElasticityTensor
youngs_modulus = 1.88e11
poissons_ratio = 0.236
block = clad
[]
[clad_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'clad_ht9creep'
block = clad
[]
[fast_flux]
type = FastNeutronFlux
block = clad
factor = 2.47e19
[]
[clad_ht9creep]
type = HT9CreepUpdate
block = clad
temperature = temp
[]
[thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = clad
thermal_expansion_coeff = 1.2e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = clad_thermal_eigenstrain
[]
[clad_thermal]
type = HT9Thermal
block = clad
temperature = temp
[]
[clad_density]
type = StrainAdjustedDensity
block = clad
strain_free_density = 7874.0
[]
[]
[Preconditioning]
[SMP]
type = SMP
full = true
[]
[]
[Executioner]
type = Transient
solve_type = 'PJFNK'
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -ksp_gmres_restart'
petsc_options_value = 'lu superlu_dist 51'
line_search = 'none'
l_max_its = 60
l_tol = 8e-3
nl_max_its = 40
nl_rel_tol = 5e-4
nl_abs_tol = 1e-7
end_time = 42000400
dtmin = 100
dtmax = 5e5
[Quadrature]
order = fifth
side_order = seventh
[]
[TimeStepper]
type = IterationAdaptiveDT
timestep_limiting_postprocessor = creep_timestep
dt = 1e2
time_t = '0 1e5 41990400 42000400'
time_dt = '1e2 1e2 1e2 1e2'
iteration_window = 4
optimal_iterations = 10
[]
[]
[Postprocessors]
[ave_temp_interior] # average temperature of the cladding interior and all pellet exteriors
type = LayeredSideAverageValuePostprocessor
boundary = 9
variable = temp
execute_on = 'initial linear'
fuel_pin_geometry = pin_geometry
[]
[clad_inner_vol] # volume inside of cladding
type = LayeredInternalVolumePostprocessor
boundary = 7
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[pellet_volume] # fuel pellet total volume
type = LayeredInternalVolumePostprocessor
boundary = 8
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[avg_clad_temp] # average temperature of cladding interior
type = LayeredSideAverageValuePostprocessor
boundary = 7
variable = temp
fuel_pin_geometry = pin_geometry
execute_on = 'initial linear'
[]
[gas_volume]
type = LayeredInternalVolumePostprocessor
boundary = 9
execute_on = 'initial linear'
component = 0
out_of_plane_strain = strain_yy
fuel_pin_geometry = pin_geometry
addition = ${gas_addition}
[]
[flux_from_clad] # area integrated heat flux from the cladding
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 5
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[flux_from_fuel] # area integrated heat flux from the fuel
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 10
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[rod_total_power]
type = LayeredElementIntegralPowerPostprocessor
variable = temp
fission_rate = fission_rate
block = fuel
fuel_pin_geometry = pin_geometry
[]
[approx_FCT]
type = AverageNodalVariableValue
boundary = 12
variable = temp
[]
[max_approx_FCT]
type = TimeExtremeValue
value_type = max
postprocessor = approx_FCT
[]
[ave_FST]
type = SideAverageValue
boundary = 10
variable = temp
[]
[max_ave_FST]
type = TimeExtremeValue
value_type = max
postprocessor = ave_FST
[]
[ave_CIT]
type = SideAverageValue
boundary = 5
variable = temp
[]
[max_ave_CIT]
type = TimeExtremeValue
value_type = max
postprocessor = ave_CIT
[]
[peak_clad_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = clad
[]
[peak_fuel_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = fuel
[]
[max_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = max
block = fuel
[]
[min_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = min
block = fuel
[]
[peak_porosity]
type = ElementExtremeValue
variable = porosity
value_type = max
block = fuel
[]
[clad_fuel_gap]
type = NodalExtremeValue
variable = penetration
boundary = 10
[]
[max_cont_press]
type = NodalExtremeValue
variable = contact_pressure
boundary = 10
[]
[LHGR_W_per_cm]
type = FunctionValuePostprocessor
function = power_history
scale_factor = 0.01
[]
[average_burnup]
type = ElementAverageValue
block = fuel
variable = burnup
[]
[fis_gas_produced]
type = LayeredElementIntegralFisGasProducePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
[]
[fis_gas_released]
type = LayeredElementIntegralFisGasReleasePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
execute_on = 'initial timestep_end'
[]
[creep_timestep]
type = MaterialTimeStepPostprocessor
block = fuel
[]
[hydrostatic_stress]
type = ElementAverageValue
variable = hydrostatic_stress
execute_on = 'initial timestep_end'
block = fuel
[]
[solid_swelling]
type = ElementAverageValue
variable = solid_swell
block = fuel
[]
[gas_swelling]
type = ElementAverageValue
variable = gas_swell
block = fuel
[]
[volumetric_strain]
type = ElementAverageValue
variable = volumetric_strain
block = fuel
[]
[fission_rate]
type = ElementAverageValue
variable = fission_rate
block = fuel
[]
[porosity]
type = ElementAverageValue
variable = porosity
block = fuel
[]
[gaseous_porosity]
type = ElementAverageValue
variable = gaseous_porosity
block = fuel
[]
[fis_gas_percent]
type = FGRPercent
fission_gas_released = fis_gas_released
fission_gas_generated = fis_gas_produced
[]
[max_clad_hoop_creep]
type = ElementExtremeValue
value_type = max
block = clad
variable = hoop_creep_strain
[]
[max_total_hoop_strain]
type = ElementExtremeValue
value_type = max
block = clad
variable = total_hoop_strain
[]
[]
[VectorPostprocessors]
[clad_x_disp]
type = NodalValueSampler
variable = disp_x
boundary = 2
sort_by = y
outputs = 'vec1'
[]
[fuel_cl_temp]
type = NodalValueSampler
variable = temp
boundary = 12
sort_by = y
outputs = 'vec2'
[]
[fuel_surf_temp]
type = NodalValueSampler
variable = temp
boundary = 10
sort_by = y
outputs = 'vec3'
[]
[clad_inn_temp]
type = NodalValueSampler
variable = temp
boundary = 5
sort_by = y
outputs = 'vec4'
[]
[clad_out_temp]
type = NodalValueSampler
variable = temp
boundary = 2
sort_by = y
outputs = 'vec5'
[]
[fuel_surf_disp_x]
type = NodalValueSampler
variable = disp_x
boundary = 10
sort_by = y
outputs = 'vec7'
[]
[]
[PerformanceMetricOutputs]
outputs = performance_metrics_file
[]
[Outputs]
time_step_interval = 10
color = true
exodus = true
perf_graph = true
csv = true
sync_times = '1e3 5e3 1e4 5e4 1e5 5e6 1e6 5e6 1e7 2e7 3e7 4e7 41990400 41991000'
file_base = x441_${group_name}_1_5D
[out2]
type = CSV
file_base = x441_${group_name}_1_5D_out2
time_step_interval = 1
[]
[console]
type = Console
max_rows = 25
time_step_interval = 1
output_linear = true
[]
[chkfile]
type = CSV
file_base = x441_${group_name}_1_5D_chkfile
show = 'max_approx_FCT max_ave_FST max_ave_CIT average_burnup fis_gas_percent max_clad_hoop_creep max_total_hoop_strain'
execute_on = 'FINAL'
[]
[performance_metrics_file]
type = CSV
file_base = x441_${group_name}_1_5D_performance_metrics
show = 'simulation_alive_time number_linear_iterations number_nonlinear_iterations time_step_size total_linear_iterations total_nonlinear_iterations physical_memory_use number_dofs number_nonlinear_variables residual_compute_time jacobian_compute_time'
[]
[vec1]
type = CSV
file_base = x441_${group_name}_1_5D_vec1
execute_on = 'FINAL'
[]
[vec2]
type = CSV
file_base = x441_${group_name}_1_5D_vec2
execute_on = 'FINAL'
[]
[vec3]
type = CSV
file_base = x441_${group_name}_1_5D_vec3
execute_on = 'FINAL'
[]
[vec4]
type = CSV
file_base = x441_${group_name}_1_5D_vec4
execute_on = 'FINAL'
[]
[vec5]
type = CSV
file_base = x441_${group_name}_1_5D_vec5
execute_on = 'FINAL'
[]
[vec7]
type = CSV
file_base = x441_${group_name}_1_5D_vec7
execute_on = 'FINAL'
[]
[]
[Debug]
show_var_residual = 'disp_x temp'
show_var_residual_norms = true
[]
(assessment/metallic_fuel/EBRII/X441/analysis/x441_base_1_5D.i)
initial_fuel_density = 15800.0
[GlobalParams]
order = FIRST
family = LAGRANGE
energy_per_fission = 3.2e-11 # J/fission
displacements = 'disp_x'
[]
[Mesh]
coord_type = RZ
# Nominal Design Geometric Parameters (X441)
[layered1D_mesh]
type = Layered1DMeshGenerator
clad_thickness = ${clad_thickness}
pellet_outer_radius = ${pellet_outer_radius}
fuel_height = ${pellet_height}
plenum_height = ${clad_top_gap_height}
clad_gap_width = ${clad_gap_width}
# meshing parameters
clad_mesh_density = customize
pellet_mesh_density = customize
nx_p = 6
nx_c = 4
slices_per_block = 10
elem_type = EDGE2
[]
# mesh options
patch_update_strategy = auto
partitioner = centroid
centroid_partitioner_direction = y
[]
[UserObjects]
[pin_geometry]
type = Layered1DFuelPinGeometry
mesh_generator = layered1D_mesh
[]
[]
[Variables]
[temp]
initial_condition = 298
[]
[]
[AuxVariables]
[disp_y]
[]
[disp_z]
[]
# Aux variables for output
[gap_cond]
order = CONSTANT
family = MONOMIAL
[]
[coolant_htc]
order = CONSTANT
family = MONOMIAL
[]
[cumulative_damage_index]
order = CONSTANT
family = MONOMIAL
[]
[solid_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[gas_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[volumetric_strain]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[hoop_stress]
order = CONSTANT
family = MONOMIAL
[]
[hoop_creep_strain]
order = CONSTANT
family = MONOMIAL
[]
[hoop_elastic_strain]
order = CONSTANT
family = MONOMIAL
[]
[total_hoop_strain]
order = CONSTANT
family = MONOMIAL
[]
[]
[Functions]
[power_history]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '0 44722 44722 0'
[]
[coolant_press_ramp]
type = PiecewiseLinear
x = '0 42000400'
y = '1.0 1.0'
[]
[coolant_temp_ramp]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '298.0 648.0 648.0 350.0'
[]
[axial_peaking_factors]
type = PowerPeakingFunction
fit = EBRII_ROW_4
pellet_length = ${pellet_height}
pellet_y_start = ${pellet_y_start}
[]
[clad_axial_pressure]
type = CladdingAxialPressureFunction
plenum_pressure = plenum_pressure
coolant_pressure = coolant_press_ramp
coolant_pressure_scaling_factor = 0.151e6
fuel_pin_geometry = pin_geometry
[]
[fuel_axial_pressure]
type = ParsedFunction
expression = plenum_pressure
symbol_names = plenum_pressure
symbol_values = plenum_pressure
[]
[]
[Physics/SolidMechanics/Layered1D]
[fuel]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = fuel_axial_pressure
generate_output = 'effective_creep_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'
block = fuel
eigenstrain_names = 'fuel_thermal_strain gas_swelling_eigenstrain solid_swelling_eigenstrain'
mesh_generator = layered1D_mesh
[]
[clad]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = clad_axial_pressure
generate_output = 'effective_creep_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'
block = clad
eigenstrain_names = 'clad_thermal_eigenstrain'
mesh_generator = layered1D_mesh
[]
[]
[Kernels]
# Define kernels for the various terms in the PDE system
[heat]
type = HeatConduction
variable = temp
[]
[heat_ie]
type = HeatConductionTimeDerivative
variable = temp
[]
[heat_source]
type = NeutronHeatSource
variable = temp
block = fuel
fission_rate = fission_rate
[]
[]
[AuxKernels]
[conductance]
type = MaterialRealAux
property = gap_conductance
variable = gap_cond
boundary = 10
[]
[gas_swell]
type = MaterialRealAux
variable = gas_swell
property = gas_swelling
execute_on = timestep_end
[]
[solid_swell]
type = MaterialRealAux
variable = solid_swell
property = solid_swelling
execute_on = timestep_end
[]
[volumetric_strain]
type = RankTwoScalarAux
rank_two_tensor = total_strain
variable = volumetric_strain
scalar_type = VolumetricStrain
execute_on = timestep_end
block = fuel
[]
[hoop_stress]
type = RankTwoAux
rank_two_tensor = stress
variable = hoop_stress
index_j = 2
index_i = 2
execute_on = timestep_end
[]
[hoop_creep_strain]
type = RankTwoAux
rank_two_tensor = creep_strain
variable = hoop_creep_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[hoop_elastic_strain]
type = RankTwoAux
rank_two_tensor = elastic_strain
variable = hoop_elastic_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[total_hoop_strain]
type = RankTwoAux
rank_two_tensor = total_strain
variable = total_hoop_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[]
[Contact]
[pellet_clad_mechanical]
primary = 5
secondary = 10
penalty = 1e12
model = frictionless
formulation = kinematic
[]
[]
[ThermalContact]
[thermal_contact]
type = GapHeatTransfer
variable = temp
primary = 5
secondary = 10
quadrature = true
gap_conductivity = 61.0
min_gap = ${clad_gap_width}
[]
[]
[BCs]
[no_x_all]
type = DirichletBC
variable = disp_x
boundary = 12
value = 0.0
[]
[Pressure]
[coolantPressure]
boundary = '2'
function = coolant_press_ramp
[]
[]
[PlenumPressure]
[plenumPressure]
boundary = 9
initial_pressure = 0.084e6 # Pa
startup_time = 0
R = 8.3143
temperature = ave_temp_interior
volume = gas_volume
output = plenum_pressure
material_input = fis_gas_released
[]
[]
[]
[CoolantChannel]
[convective_clad_surface]
boundary = '2'
variable = temp
inlet_temperature = coolant_temp_ramp
inlet_pressure = coolant_press_ramp
inlet_massflux = 5261.5 # kg/m^2-sec
coolant_material = sodium
rod_diameter = 5.84e-3 # m
rod_pitch = 7.48e-3 # m (Pitch-to-diameter Ratio = 1.28)
linear_heat_rate = power_history
axial_power_profile = axial_peaking_factors
subchannel_geometry = triangular
[]
[]
[Materials]
[fission_rate]
type = UPuZrFissionRate
rod_linear_power = power_history
axial_power_profile = axial_peaking_factors
pellet_radius = ${pellet_outer_radius}
X_Zr = ${X_Zr}
X_Pu_function = ${X_Pu}
block = fuel
outputs = all
[]
[burnup]
type = UPuZrBurnup
initial_X_Zr = ${X_Zr}
initial_X_Pu = ${X_Pu}
density = ${initial_fuel_density}
block = fuel
outputs = all
[]
[fuel_elasticity_tensor]
type = UPuZrElasticityTensor
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
block = fuel
temperature = temp
[]
[fuel_inelastic_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'fuel_upuzrcreep'
block = fuel
[]
[fuel_upuzrcreep]
type = UPuZrCreepUpdate
block = fuel
temperature = temp
porosity = porosity
max_inelastic_increment = 1e-3
[]
[fuel_thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = fuel
thermal_expansion_coeff = 1.18e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = fuel_thermal_strain
[]
[gas_swelling]
type = UPuZrGaseousEigenstrain
temperature = temp
eigenstrain_name = gas_swelling_eigenstrain
initial_porosity = 0.0
bubble_number_density = 1e20
outputs = all
output_properties = 'porosity gaseous_porosity'
block = fuel
[]
[solid_swelling]
type = BurnupDependentEigenstrain
eigenstrain_name = solid_swelling_eigenstrain
block = fuel
swelling_name = 'solid_swelling'
[]
[metal_fuel_thermal]
type = UPuZrThermal
block = fuel
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
spheat_model = savage
thcond_model = lanl
porosity = porosity
temperature = temp
[]
[fuel_density]
type = StrainAdjustedDensity
block = fuel
strain_free_density = ${initial_fuel_density}
[]
[fission_gas_behavior]
type = UPuZrFissionGasRelease
block = fuel
fission_rate = fission_rate
[]
[clad_elasticity_tensor]
type = ComputeIsotropicElasticityTensor
youngs_modulus = 1.88e11
poissons_ratio = 0.236
block = clad
[]
[clad_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'clad_ht9creep'
block = clad
[]
[fast_flux]
type = FastNeutronFlux
block = clad
factor = 2.47e19
[]
[clad_ht9creep]
type = HT9CreepUpdate
block = clad
temperature = temp
[]
[thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = clad
thermal_expansion_coeff = 1.2e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = clad_thermal_eigenstrain
[]
[clad_thermal]
type = HT9Thermal
block = clad
temperature = temp
[]
[clad_density]
type = StrainAdjustedDensity
block = clad
strain_free_density = 7874.0
[]
[]
[Preconditioning]
[SMP]
type = SMP
full = true
[]
[]
[Executioner]
type = Transient
solve_type = 'PJFNK'
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -ksp_gmres_restart'
petsc_options_value = 'lu superlu_dist 51'
line_search = 'none'
l_max_its = 60
l_tol = 8e-3
nl_max_its = 40
nl_rel_tol = 5e-4
nl_abs_tol = 1e-7
end_time = 42000400
dtmin = 100
dtmax = 5e5
[Quadrature]
order = fifth
side_order = seventh
[]
[TimeStepper]
type = IterationAdaptiveDT
timestep_limiting_postprocessor = creep_timestep
dt = 1e2
time_t = '0 1e5 41990400 42000400'
time_dt = '1e2 1e2 1e2 1e2'
iteration_window = 4
optimal_iterations = 10
[]
[]
[Postprocessors]
[ave_temp_interior] # average temperature of the cladding interior and all pellet exteriors
type = LayeredSideAverageValuePostprocessor
boundary = 9
variable = temp
execute_on = 'initial linear'
fuel_pin_geometry = pin_geometry
[]
[clad_inner_vol] # volume inside of cladding
type = LayeredInternalVolumePostprocessor
boundary = 7
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[pellet_volume] # fuel pellet total volume
type = LayeredInternalVolumePostprocessor
boundary = 8
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[avg_clad_temp] # average temperature of cladding interior
type = LayeredSideAverageValuePostprocessor
boundary = 7
variable = temp
fuel_pin_geometry = pin_geometry
execute_on = 'initial linear'
[]
[gas_volume]
type = LayeredInternalVolumePostprocessor
boundary = 9
execute_on = 'initial linear'
component = 0
out_of_plane_strain = strain_yy
fuel_pin_geometry = pin_geometry
addition = ${gas_addition}
[]
[flux_from_clad] # area integrated heat flux from the cladding
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 5
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[flux_from_fuel] # area integrated heat flux from the fuel
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 10
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[rod_total_power]
type = LayeredElementIntegralPowerPostprocessor
variable = temp
fission_rate = fission_rate
block = fuel
fuel_pin_geometry = pin_geometry
[]
[approx_FCT]
type = AverageNodalVariableValue
boundary = 12
variable = temp
[]
[max_approx_FCT]
type = TimeExtremeValue
value_type = max
postprocessor = approx_FCT
[]
[ave_FST]
type = SideAverageValue
boundary = 10
variable = temp
[]
[max_ave_FST]
type = TimeExtremeValue
value_type = max
postprocessor = ave_FST
[]
[ave_CIT]
type = SideAverageValue
boundary = 5
variable = temp
[]
[max_ave_CIT]
type = TimeExtremeValue
value_type = max
postprocessor = ave_CIT
[]
[peak_clad_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = clad
[]
[peak_fuel_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = fuel
[]
[max_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = max
block = fuel
[]
[min_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = min
block = fuel
[]
[peak_porosity]
type = ElementExtremeValue
variable = porosity
value_type = max
block = fuel
[]
[clad_fuel_gap]
type = NodalExtremeValue
variable = penetration
boundary = 10
[]
[max_cont_press]
type = NodalExtremeValue
variable = contact_pressure
boundary = 10
[]
[LHGR_W_per_cm]
type = FunctionValuePostprocessor
function = power_history
scale_factor = 0.01
[]
[average_burnup]
type = ElementAverageValue
block = fuel
variable = burnup
[]
[fis_gas_produced]
type = LayeredElementIntegralFisGasProducePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
[]
[fis_gas_released]
type = LayeredElementIntegralFisGasReleasePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
execute_on = 'initial timestep_end'
[]
[creep_timestep]
type = MaterialTimeStepPostprocessor
block = fuel
[]
[hydrostatic_stress]
type = ElementAverageValue
variable = hydrostatic_stress
execute_on = 'initial timestep_end'
block = fuel
[]
[solid_swelling]
type = ElementAverageValue
variable = solid_swell
block = fuel
[]
[gas_swelling]
type = ElementAverageValue
variable = gas_swell
block = fuel
[]
[volumetric_strain]
type = ElementAverageValue
variable = volumetric_strain
block = fuel
[]
[fission_rate]
type = ElementAverageValue
variable = fission_rate
block = fuel
[]
[porosity]
type = ElementAverageValue
variable = porosity
block = fuel
[]
[gaseous_porosity]
type = ElementAverageValue
variable = gaseous_porosity
block = fuel
[]
[fis_gas_percent]
type = FGRPercent
fission_gas_released = fis_gas_released
fission_gas_generated = fis_gas_produced
[]
[max_clad_hoop_creep]
type = ElementExtremeValue
value_type = max
block = clad
variable = hoop_creep_strain
[]
[max_total_hoop_strain]
type = ElementExtremeValue
value_type = max
block = clad
variable = total_hoop_strain
[]
[]
[VectorPostprocessors]
[clad_x_disp]
type = NodalValueSampler
variable = disp_x
boundary = 2
sort_by = y
outputs = 'vec1'
[]
[fuel_cl_temp]
type = NodalValueSampler
variable = temp
boundary = 12
sort_by = y
outputs = 'vec2'
[]
[fuel_surf_temp]
type = NodalValueSampler
variable = temp
boundary = 10
sort_by = y
outputs = 'vec3'
[]
[clad_inn_temp]
type = NodalValueSampler
variable = temp
boundary = 5
sort_by = y
outputs = 'vec4'
[]
[clad_out_temp]
type = NodalValueSampler
variable = temp
boundary = 2
sort_by = y
outputs = 'vec5'
[]
[fuel_surf_disp_x]
type = NodalValueSampler
variable = disp_x
boundary = 10
sort_by = y
outputs = 'vec7'
[]
[]
[PerformanceMetricOutputs]
outputs = performance_metrics_file
[]
[Outputs]
time_step_interval = 10
color = true
exodus = true
perf_graph = true
csv = true
sync_times = '1e3 5e3 1e4 5e4 1e5 5e6 1e6 5e6 1e7 2e7 3e7 4e7 41990400 41991000'
file_base = x441_${group_name}_1_5D
[out2]
type = CSV
file_base = x441_${group_name}_1_5D_out2
time_step_interval = 1
[]
[console]
type = Console
max_rows = 25
time_step_interval = 1
output_linear = true
[]
[chkfile]
type = CSV
file_base = x441_${group_name}_1_5D_chkfile
show = 'max_approx_FCT max_ave_FST max_ave_CIT average_burnup fis_gas_percent max_clad_hoop_creep max_total_hoop_strain'
execute_on = 'FINAL'
[]
[performance_metrics_file]
type = CSV
file_base = x441_${group_name}_1_5D_performance_metrics
show = 'simulation_alive_time number_linear_iterations number_nonlinear_iterations time_step_size total_linear_iterations total_nonlinear_iterations physical_memory_use number_dofs number_nonlinear_variables residual_compute_time jacobian_compute_time'
[]
[vec1]
type = CSV
file_base = x441_${group_name}_1_5D_vec1
execute_on = 'FINAL'
[]
[vec2]
type = CSV
file_base = x441_${group_name}_1_5D_vec2
execute_on = 'FINAL'
[]
[vec3]
type = CSV
file_base = x441_${group_name}_1_5D_vec3
execute_on = 'FINAL'
[]
[vec4]
type = CSV
file_base = x441_${group_name}_1_5D_vec4
execute_on = 'FINAL'
[]
[vec5]
type = CSV
file_base = x441_${group_name}_1_5D_vec5
execute_on = 'FINAL'
[]
[vec7]
type = CSV
file_base = x441_${group_name}_1_5D_vec7
execute_on = 'FINAL'
[]
[]
[Debug]
show_var_residual = 'disp_x temp'
show_var_residual_norms = true
[]
(test/tests/fission_gas_1d/fgr_1d.i)
#This tests the fission gas produced and released postprocessors
#and for 1d and 2d models with the same volume.
#The reported values should be the same for the
#1d and 2d cases. See corresponding excel calculation
#in this directory.
[Problem]
type = ReferenceResidualProblem
reference_vector = 'ref'
extra_tag_vectors = 'ref'
[]
[UserObjects]
[pin_geometry]
type = Layered1DFuelPinGeometry
include_clad = false
mesh_generator = layered1D_mesh
[]
[]
[Mesh]
coord_type = RZ
[layered1D_mesh]
type = Layered1DMeshGenerator
include_clad = false
include_plenum = false
slices_per_block = 1
fuel_height = 3e-3
pellet_outer_radius = 3e-3
pellet_bottom_coor = 0.0
pellet_mesh_density = customize
nx_p = 4
elem_type = EDGE2
[]
[]
[Functions]
[Temp_func]
type = ParsedFunction
expression = '1700'
[]
[Fiss_func]
type = ParsedFunction
expression = '2.5e19'
[]
[Porosity_func]
type = PiecewiseLinear
x = '0 1e7'
y = '0.15 0.241'
[]
[]
[Variables]
[temp]
order = FIRST
family = LAGRANGE
initial_condition = 1700
[]
[]
[AuxVariables]
[porosity_material_aux]
order = CONSTANT
family = MONOMIAL
[]
[]
[Kernels]
[heat]
type = HeatConduction
variable = temp
[]
[]
[AuxKernels]
[porosity_material_aux]
type = MaterialRealAux
variable = porosity_material_aux
property = porosity
execute_on = 'initial timestep_begin'
[]
[]
[BCs]
[temp_bc]
type = FunctionDirichletBC
variable = temp
function = Temp_func
boundary = 10
[]
[]
[Materials]
[fission_rate]
type = GenericFunctionMaterial
prop_names = 'fission_rate thermal_conductivity porosity'
prop_values = 'Fiss_func 5 Porosity_func'
outputs = all
[]
[fission_gas_behavior]
type = UPuZrFissionGasRelease
fission_rate = fission_rate
fractional_yield = 0.3017
[]
[]
[Executioner]
type = Transient
solve_type = 'PJFNK'
l_tol = 1e-4
nl_abs_tol = 1e-5
nl_rel_tol = 1e-5
start_time = 0.0
dt = 5e6
end_time = 1e7
[]
[Postprocessors]
[gas_produced]
type = LayeredElementIntegralFisGasProducePostprocessor
fuel_pin_geometry = pin_geometry
[]
[gas_released]
type = LayeredElementIntegralFisGasReleasePostprocessor
fuel_pin_geometry = pin_geometry
[]
[porosity]
type = ElementAverageValue
variable = porosity_material_aux
[]
[]
[Outputs]
[csv]
type = CSV
execute_on = final
file_base = fgr_1d_out
[]
[]
(assessment/metallic_fuel/EBRII/X441/analysis/group_B/x441_1_5D_B.i)
initial_fuel_density = 15800.0
[GlobalParams]
order = FIRST
family = LAGRANGE
energy_per_fission = 3.2e-11 # J/fission
displacements = 'disp_x'
[]
[Mesh]
coord_type = RZ
# Nominal Design Geometric Parameters (X441)
[layered1D_mesh]
type = Layered1DMeshGenerator
clad_thickness = ${clad_thickness}
pellet_outer_radius = ${pellet_outer_radius}
fuel_height = ${pellet_height}
plenum_height = ${clad_top_gap_height}
clad_gap_width = ${clad_gap_width}
# meshing parameters
clad_mesh_density = customize
pellet_mesh_density = customize
nx_p = 6
nx_c = 4
slices_per_block = 10
elem_type = EDGE2
[]
# mesh options
patch_update_strategy = auto
partitioner = centroid
centroid_partitioner_direction = y
[]
[UserObjects]
[pin_geometry]
type = Layered1DFuelPinGeometry
mesh_generator = layered1D_mesh
[]
[]
[Variables]
[temp]
initial_condition = 298
[]
[]
[AuxVariables]
[disp_y]
[]
[disp_z]
[]
# Aux variables for output
[gap_cond]
order = CONSTANT
family = MONOMIAL
[]
[coolant_htc]
order = CONSTANT
family = MONOMIAL
[]
[cumulative_damage_index]
order = CONSTANT
family = MONOMIAL
[]
[solid_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[gas_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[volumetric_strain]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[hoop_stress]
order = CONSTANT
family = MONOMIAL
[]
[hoop_creep_strain]
order = CONSTANT
family = MONOMIAL
[]
[hoop_elastic_strain]
order = CONSTANT
family = MONOMIAL
[]
[total_hoop_strain]
order = CONSTANT
family = MONOMIAL
[]
[]
[Functions]
[power_history]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '0 44722 44722 0'
[]
[coolant_press_ramp]
type = PiecewiseLinear
x = '0 42000400'
y = '1.0 1.0'
[]
[coolant_temp_ramp]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '298.0 648.0 648.0 350.0'
[]
[axial_peaking_factors]
type = PowerPeakingFunction
fit = EBRII_ROW_4
pellet_length = ${pellet_height}
pellet_y_start = ${pellet_y_start}
[]
[clad_axial_pressure]
type = CladdingAxialPressureFunction
plenum_pressure = plenum_pressure
coolant_pressure = coolant_press_ramp
coolant_pressure_scaling_factor = 0.151e6
fuel_pin_geometry = pin_geometry
[]
[fuel_axial_pressure]
type = ParsedFunction
expression = plenum_pressure
symbol_names = plenum_pressure
symbol_values = plenum_pressure
[]
[]
[Physics/SolidMechanics/Layered1D]
[fuel]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = fuel_axial_pressure
generate_output = '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'
block = fuel
eigenstrain_names = 'fuel_thermal_strain gas_swelling_eigenstrain solid_swelling_eigenstrain'
mesh_generator = layered1D_mesh
[]
[clad]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = clad_axial_pressure
generate_output = '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'
block = clad
eigenstrain_names = 'clad_thermal_eigenstrain'
mesh_generator = layered1D_mesh
[]
[]
[Kernels]
# Define kernels for the various terms in the PDE system
[heat]
type = HeatConduction
variable = temp
[]
[heat_ie]
type = HeatConductionTimeDerivative
variable = temp
[]
[heat_source]
type = NeutronHeatSource
variable = temp
block = fuel
fission_rate = fission_rate
[]
[]
[AuxKernels]
[conductance]
type = MaterialRealAux
property = gap_conductance
variable = gap_cond
boundary = 10
[]
[gas_swell]
type = MaterialRealAux
variable = gas_swell
property = gas_swelling
execute_on = timestep_end
[]
[solid_swell]
type = MaterialRealAux
variable = solid_swell
property = solid_swelling
execute_on = timestep_end
[]
[volumetric_strain]
type = RankTwoScalarAux
rank_two_tensor = total_strain
variable = volumetric_strain
scalar_type = VolumetricStrain
execute_on = timestep_end
block = fuel
[]
[hoop_stress]
type = RankTwoAux
rank_two_tensor = stress
variable = hoop_stress
index_j = 2
index_i = 2
execute_on = timestep_end
[]
[hoop_creep_strain]
type = RankTwoAux
rank_two_tensor = creep_strain
variable = hoop_creep_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[hoop_elastic_strain]
type = RankTwoAux
rank_two_tensor = elastic_strain
variable = hoop_elastic_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[total_hoop_strain]
type = RankTwoAux
rank_two_tensor = total_strain
variable = total_hoop_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[]
[Contact]
[pellet_clad_mechanical]
primary = 5
secondary = 10
penalty = 1e12
model = frictionless
formulation = kinematic
[]
[]
[ThermalContact]
[thermal_contact]
type = GapHeatTransfer
variable = temp
primary = 5
secondary = 10
quadrature = true
gap_conductivity = 61.0
min_gap = ${clad_gap_width}
[]
[]
[BCs]
[no_x_all]
type = DirichletBC
variable = disp_x
boundary = 12
value = 0.0
[]
[Pressure]
[coolantPressure]
boundary = '2'
function = coolant_press_ramp
[]
[]
[PlenumPressure]
[plenumPressure]
boundary = 9
initial_pressure = 0.084e6 # Pa
startup_time = 0
R = 8.3143
temperature = ave_temp_interior
volume = gas_volume
output = plenum_pressure
material_input = fis_gas_released
[]
[]
[]
[CoolantChannel]
[convective_clad_surface]
boundary = '2'
variable = temp
inlet_temperature = coolant_temp_ramp
inlet_pressure = coolant_press_ramp
inlet_massflux = 5261.5 # kg/m^2-sec
coolant_material = sodium
rod_diameter = 5.84e-3 # m
rod_pitch = 7.48e-3 # m (Pitch-to-diameter Ratio = 1.28)
linear_heat_rate = power_history
axial_power_profile = axial_peaking_factors
subchannel_geometry = triangular
[]
[]
[Materials]
[fission_rate]
type = UPuZrFissionRate
rod_linear_power = power_history
axial_power_profile = axial_peaking_factors
pellet_radius = ${pellet_outer_radius}
X_Zr = ${X_Zr}
X_Pu_function = ${X_Pu}
block = fuel
outputs = all
[]
[burnup]
type = UPuZrBurnup
initial_X_Zr = ${X_Zr}
initial_X_Pu = ${X_Pu}
density = ${initial_fuel_density}
block = fuel
outputs = all
[]
[fuel_elasticity_tensor]
type = UPuZrElasticityTensor
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
block = fuel
temperature = temp
[]
[fuel_inelastic_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'fuel_upuzrcreep'
block = fuel
[]
[fuel_upuzrcreep]
type = UPuZrCreepUpdate
block = fuel
temperature = temp
porosity = porosity
max_inelastic_increment = 1e-3
[]
[fuel_thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = fuel
thermal_expansion_coeff = 1.18e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = fuel_thermal_strain
[]
[gas_swelling]
type = UPuZrGaseousEigenstrain
temperature = temp
eigenstrain_name = gas_swelling_eigenstrain
initial_porosity = 0.0
bubble_number_density = 1e20
outputs = all
output_properties = 'porosity gaseous_porosity'
block = fuel
[]
[solid_swelling]
type = BurnupDependentEigenstrain
eigenstrain_name = solid_swelling_eigenstrain
block = fuel
swelling_name = 'solid_swelling'
[]
[metal_fuel_thermal]
type = UPuZrThermal
block = fuel
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
spheat_model = savage
thcond_model = lanl
porosity = porosity
temperature = temp
[]
[fuel_density]
type = StrainAdjustedDensity
block = fuel
strain_free_density = 15800
[]
[fission_gas_behavior]
type = UPuZrFissionGasRelease
block = fuel
fission_rate = fission_rate
[]
[clad_elasticity_tensor]
type = ComputeIsotropicElasticityTensor
youngs_modulus = 1.88e11
poissons_ratio = 0.236
block = clad
[]
[clad_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'clad_ht9creep'
block = clad
[]
[fast_flux]
type = FastNeutronFlux
block = clad
factor = 2.47e19
[]
[clad_ht9creep]
type = HT9CreepUpdate
block = clad
temperature = temp
[]
[thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = clad
thermal_expansion_coeff = 1.2e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = clad_thermal_eigenstrain
[]
[clad_thermal]
type = HT9Thermal
block = clad
temperature = temp
[]
[clad_density]
type = StrainAdjustedDensity
block = clad
strain_free_density = 7874.0
[]
[]
[Preconditioning]
[SMP]
type = SMP
full = true
[]
[]
[Executioner]
type = Transient
solve_type = 'PJFNK'
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -ksp_gmres_restart'
petsc_options_value = 'lu superlu_dist 51'
line_search = 'none'
l_max_its = 60
l_tol = 8e-3
nl_max_its = 40
nl_rel_tol = 5e-4
nl_abs_tol = 1e-7
end_time = 42000400
dtmin = 100
dtmax = 5e5
[Quadrature]
order = fifth
side_order = seventh
[]
[TimeStepper]
type = IterationAdaptiveDT
timestep_limiting_postprocessor = creep_timestep
dt = 1e2
time_t = '0 1e5 41990400 42000400'
time_dt = '1e2 1e2 1e2 1e2'
iteration_window = 4
optimal_iterations = 10
[]
[]
[Postprocessors]
[ave_temp_interior] # average temperature of the cladding interior and all pellet exteriors
type = LayeredSideAverageValuePostprocessor
boundary = 9
variable = temp
execute_on = 'initial linear'
fuel_pin_geometry = pin_geometry
[]
[clad_inner_vol] # volume inside of cladding
type = LayeredInternalVolumePostprocessor
boundary = 7
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[pellet_volume] # fuel pellet total volume
type = LayeredInternalVolumePostprocessor
boundary = 8
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[avg_clad_temp] # average temperature of cladding interior
type = LayeredSideAverageValuePostprocessor
boundary = 7
variable = temp
fuel_pin_geometry = pin_geometry
execute_on = 'initial linear'
[]
[gas_volume]
type = LayeredInternalVolumePostprocessor
boundary = 9
execute_on = 'initial linear'
component = 0
out_of_plane_strain = strain_yy
fuel_pin_geometry = pin_geometry
addition = ${gas_addition}
[]
[flux_from_clad] # area integrated heat flux from the cladding
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 5
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[flux_from_fuel] # area integrated heat flux from the fuel
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 10
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[rod_total_power]
type = LayeredElementIntegralPowerPostprocessor
variable = temp
fission_rate = fission_rate
block = fuel
fuel_pin_geometry = pin_geometry
[]
[approx_FCT]
type = AverageNodalVariableValue
boundary = 12
variable = temp
[]
[max_approx_FCT]
type = TimeExtremeValue
value_type = max
postprocessor = approx_FCT
[]
[ave_FST]
type = SideAverageValue
boundary = 10
variable = temp
[]
[max_ave_FST]
type = TimeExtremeValue
value_type = max
postprocessor = ave_FST
[]
[ave_CIT]
type = SideAverageValue
boundary = 5
variable = temp
[]
[max_ave_CIT]
type = TimeExtremeValue
value_type = max
postprocessor = ave_CIT
[]
[peak_clad_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = clad
[]
[peak_fuel_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = fuel
[]
[max_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = max
block = fuel
[]
[min_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = min
block = fuel
[]
[peak_porosity]
type = ElementExtremeValue
variable = porosity
value_type = max
block = fuel
[]
[clad_fuel_gap]
type = NodalExtremeValue
variable = penetration
boundary = 10
[]
[max_cont_press]
type = NodalExtremeValue
variable = contact_pressure
boundary = 10
[]
[LHGR_W_per_cm]
type = FunctionValuePostprocessor
function = power_history
scale_factor = 0.01
[]
[average_burnup]
type = ElementAverageValue
block = fuel
variable = burnup
[]
[fis_gas_produced]
type = LayeredElementIntegralFisGasProducePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
[]
[fis_gas_released]
type = LayeredElementIntegralFisGasReleasePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
execute_on = 'initial timestep_end'
[]
[creep_timestep]
type = MaterialTimeStepPostprocessor
block = fuel
[]
[hydrostatic_stress]
type = ElementAverageValue
variable = hydrostatic_stress
execute_on = 'initial timestep_end'
block = fuel
[]
[solid_swelling]
type = ElementAverageValue
variable = solid_swell
block = fuel
[]
[gas_swelling]
type = ElementAverageValue
variable = gas_swell
block = fuel
[]
[volumetric_strain]
type = ElementAverageValue
variable = volumetric_strain
block = fuel
[]
[fission_rate]
type = ElementAverageValue
variable = fission_rate
block = fuel
[]
[porosity]
type = ElementAverageValue
variable = porosity
block = fuel
[]
[gaseous_porosity]
type = ElementAverageValue
variable = gaseous_porosity
block = fuel
[]
[fis_gas_percent]
type = FGRPercent
fission_gas_released = fis_gas_released
fission_gas_generated = fis_gas_produced
[]
[max_clad_hoop_creep]
type = ElementExtremeValue
value_type = max
block = clad
variable = hoop_creep_strain
[]
[max_total_hoop_strain]
type = ElementExtremeValue
value_type = max
block = clad
variable = total_hoop_strain
[]
[]
[VectorPostprocessors]
[clad_x_disp]
type = NodalValueSampler
variable = disp_x
boundary = 2
sort_by = y
outputs = 'vec1'
[]
[fuel_cl_temp]
type = NodalValueSampler
variable = temp
boundary = 12
sort_by = y
outputs = 'vec2'
[]
[fuel_surf_temp]
type = NodalValueSampler
variable = temp
boundary = 10
sort_by = y
outputs = 'vec3'
[]
[clad_inn_temp]
type = NodalValueSampler
variable = temp
boundary = 5
sort_by = y
outputs = 'vec4'
[]
[clad_out_temp]
type = NodalValueSampler
variable = temp
boundary = 2
sort_by = y
outputs = 'vec5'
[]
[fuel_surf_disp_x]
type = NodalValueSampler
variable = disp_x
boundary = 10
sort_by = y
outputs = 'vec7'
[]
[]
[PerformanceMetricOutputs]
outputs = performance_metrics_file
[]
[Outputs]
time_step_interval = 10
color = true
exodus = true
perf_graph = true
csv = true
sync_times = '1e3 5e3 1e4 5e4 1e5 5e6 1e6 5e6 1e7 2e7 3e7 4e7 41990400 41991000'
file_base = x441_${group_name}_1_5D
[out2]
type = CSV
file_base = x441_${group_name}_1_5D_out2
time_step_interval = 1
[]
[console]
type = Console
max_rows = 25
time_step_interval = 1
output_linear = true
[]
[chkfile]
type = CSV
file_base = x441_${group_name}_1_5D_chkfile
show = 'max_approx_FCT max_ave_FST max_ave_CIT average_burnup fis_gas_percent max_clad_hoop_creep max_total_hoop_strain'
execute_on = 'FINAL'
[]
[performance_metrics_file]
type = CSV
file_base = x441_${group_name}_1_5D_performance_metrics
show = 'simulation_alive_time number_linear_iterations number_nonlinear_iterations time_step_size total_linear_iterations total_nonlinear_iterations physical_memory_use number_dofs number_nonlinear_variables residual_compute_time jacobian_compute_time'
[]
[vec1]
type = CSV
file_base = x441_${group_name}_1_5D_vec1
execute_on = 'FINAL'
[]
[vec2]
type = CSV
file_base = x441_${group_name}_1_5D_vec2
execute_on = 'FINAL'
[]
[vec3]
type = CSV
file_base = x441_${group_name}_1_5D_vec3
execute_on = 'FINAL'
[]
[vec4]
type = CSV
file_base = x441_${group_name}_1_5D_vec4
execute_on = 'FINAL'
[]
[vec5]
type = CSV
file_base = x441_${group_name}_1_5D_vec5
execute_on = 'FINAL'
[]
[vec7]
type = CSV
file_base = x441_${group_name}_1_5D_vec7
execute_on = 'FINAL'
[]
[]
[Debug]
show_var_residual = 'disp_x temp'
show_var_residual_norms = true
[]
(assessment/metallic_fuel/EBRII/X441/analysis/group_F/x441_1_5D_F.i)
initial_fuel_density = 15800.0
[GlobalParams]
order = FIRST
family = LAGRANGE
energy_per_fission = 3.2e-11 # J/fission
displacements = 'disp_x'
[]
[Mesh]
coord_type = RZ
# Nominal Design Geometric Parameters (X441)
[layered1D_mesh]
type = Layered1DMeshGenerator
clad_thickness = ${clad_thickness}
pellet_outer_radius = ${pellet_outer_radius}
fuel_height = ${pellet_height}
plenum_height = ${clad_top_gap_height}
clad_gap_width = ${clad_gap_width}
# meshing parameters
clad_mesh_density = customize
pellet_mesh_density = customize
nx_p = 6
nx_c = 4
slices_per_block = 10
elem_type = EDGE2
[]
# mesh options
patch_update_strategy = auto
partitioner = centroid
centroid_partitioner_direction = y
[]
[UserObjects]
[pin_geometry]
type = Layered1DFuelPinGeometry
mesh_generator = layered1D_mesh
[]
[]
[Variables]
[temp]
initial_condition = 298
[]
[]
[AuxVariables]
[disp_y]
[]
[disp_z]
[]
# Aux variables for output
[gap_cond]
order = CONSTANT
family = MONOMIAL
[]
[coolant_htc]
order = CONSTANT
family = MONOMIAL
[]
[cumulative_damage_index]
order = CONSTANT
family = MONOMIAL
[]
[solid_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[gas_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[volumetric_strain]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[hoop_stress]
order = CONSTANT
family = MONOMIAL
[]
[hoop_creep_strain]
order = CONSTANT
family = MONOMIAL
[]
[hoop_elastic_strain]
order = CONSTANT
family = MONOMIAL
[]
[total_hoop_strain]
order = CONSTANT
family = MONOMIAL
[]
[]
[Functions]
[power_history]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '0 44722 44722 0'
[]
[coolant_press_ramp]
type = PiecewiseLinear
x = '0 42000400'
y = '1.0 1.0'
[]
[coolant_temp_ramp]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '298.0 648.0 648.0 350.0'
[]
[axial_peaking_factors]
type = PowerPeakingFunction
fit = EBRII_ROW_4
pellet_length = ${pellet_height}
pellet_y_start = ${pellet_y_start}
[]
[clad_axial_pressure]
type = CladdingAxialPressureFunction
plenum_pressure = plenum_pressure
coolant_pressure = coolant_press_ramp
coolant_pressure_scaling_factor = 0.151e6
fuel_pin_geometry = pin_geometry
[]
[fuel_axial_pressure]
type = ParsedFunction
expression = plenum_pressure
symbol_names = plenum_pressure
symbol_values = plenum_pressure
[]
[]
[Physics/SolidMechanics/Layered1D]
[fuel]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = fuel_axial_pressure
generate_output = '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'
block = fuel
eigenstrain_names = 'fuel_thermal_strain gas_swelling_eigenstrain solid_swelling_eigenstrain'
mesh_generator = layered1D_mesh
[]
[clad]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = clad_axial_pressure
generate_output = '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'
block = clad
eigenstrain_names = 'clad_thermal_eigenstrain'
mesh_generator = layered1D_mesh
[]
[]
[Kernels]
# Define kernels for the various terms in the PDE system
[heat]
type = HeatConduction
variable = temp
[]
[heat_ie]
type = HeatConductionTimeDerivative
variable = temp
[]
[heat_source]
type = NeutronHeatSource
variable = temp
block = fuel
fission_rate = fission_rate
[]
[]
[AuxKernels]
[conductance]
type = MaterialRealAux
property = gap_conductance
variable = gap_cond
boundary = 10
[]
[gas_swell]
type = MaterialRealAux
variable = gas_swell
property = gas_swelling
execute_on = timestep_end
[]
[solid_swell]
type = MaterialRealAux
variable = solid_swell
property = solid_swelling
execute_on = timestep_end
[]
[volumetric_strain]
type = RankTwoScalarAux
rank_two_tensor = total_strain
variable = volumetric_strain
scalar_type = VolumetricStrain
execute_on = timestep_end
block = fuel
[]
[hoop_stress]
type = RankTwoAux
rank_two_tensor = stress
variable = hoop_stress
index_j = 2
index_i = 2
execute_on = timestep_end
[]
[hoop_creep_strain]
type = RankTwoAux
rank_two_tensor = creep_strain
variable = hoop_creep_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[hoop_elastic_strain]
type = RankTwoAux
rank_two_tensor = elastic_strain
variable = hoop_elastic_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[total_hoop_strain]
type = RankTwoAux
rank_two_tensor = total_strain
variable = total_hoop_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[]
[Contact]
[pellet_clad_mechanical]
primary = 5
secondary = 10
penalty = 1e12
model = frictionless
formulation = kinematic
[]
[]
[ThermalContact]
[thermal_contact]
type = GapHeatTransfer
variable = temp
primary = 5
secondary = 10
quadrature = true
gap_conductivity = 61.0
min_gap = ${clad_gap_width}
[]
[]
[BCs]
[no_x_all]
type = DirichletBC
variable = disp_x
boundary = 12
value = 0.0
[]
[Pressure]
[coolantPressure]
boundary = '2'
function = coolant_press_ramp
[]
[]
[PlenumPressure]
[plenumPressure]
boundary = 9
initial_pressure = 0.084e6 # Pa
startup_time = 0
R = 8.3143
temperature = ave_temp_interior
volume = gas_volume
output = plenum_pressure
material_input = fis_gas_released
[]
[]
[]
[CoolantChannel]
[convective_clad_surface]
boundary = '2'
variable = temp
inlet_temperature = coolant_temp_ramp
inlet_pressure = coolant_press_ramp
inlet_massflux = 5261.5 # kg/m^2-sec
coolant_material = sodium
rod_diameter = 5.84e-3 # m
rod_pitch = 7.48e-3 # m (Pitch-to-diameter Ratio = 1.28)
linear_heat_rate = power_history
axial_power_profile = axial_peaking_factors
subchannel_geometry = triangular
[]
[]
[Materials]
[fission_rate]
type = UPuZrFissionRate
rod_linear_power = power_history
axial_power_profile = axial_peaking_factors
pellet_radius = ${pellet_outer_radius}
X_Zr = ${X_Zr}
X_Pu_function = ${X_Pu}
block = fuel
outputs = all
[]
[burnup]
type = UPuZrBurnup
initial_X_Zr = ${X_Zr}
initial_X_Pu = ${X_Pu}
density = ${initial_fuel_density}
block = fuel
outputs = all
[]
[fuel_elasticity_tensor]
type = UPuZrElasticityTensor
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
block = fuel
temperature = temp
[]
[fuel_inelastic_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'fuel_upuzrcreep'
block = fuel
[]
[fuel_upuzrcreep]
type = UPuZrCreepUpdate
block = fuel
temperature = temp
porosity = porosity
max_inelastic_increment = 1e-3
[]
[fuel_thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = fuel
thermal_expansion_coeff = 1.18e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = fuel_thermal_strain
[]
[gas_swelling]
type = UPuZrGaseousEigenstrain
temperature = temp
eigenstrain_name = gas_swelling_eigenstrain
initial_porosity = 0.0
bubble_number_density = 1e20
outputs = all
output_properties = 'porosity gaseous_porosity'
block = fuel
[]
[solid_swelling]
type = BurnupDependentEigenstrain
eigenstrain_name = solid_swelling_eigenstrain
block = fuel
swelling_name = 'solid_swelling'
[]
[metal_fuel_thermal]
type = UPuZrThermal
block = fuel
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
spheat_model = savage
thcond_model = lanl
porosity = porosity
temperature = temp
[]
[fuel_density]
type = StrainAdjustedDensity
block = fuel
strain_free_density = 15800
[]
[fission_gas_behavior]
type = UPuZrFissionGasRelease
block = fuel
fission_rate = fission_rate
[]
[clad_elasticity_tensor]
type = ComputeIsotropicElasticityTensor
youngs_modulus = 1.88e11
poissons_ratio = 0.236
block = clad
[]
[clad_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'clad_ht9creep'
block = clad
[]
[fast_flux]
type = FastNeutronFlux
block = clad
factor = 2.47e19
[]
[clad_ht9creep]
type = HT9CreepUpdate
block = clad
temperature = temp
[]
[thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = clad
thermal_expansion_coeff = 1.2e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = clad_thermal_eigenstrain
[]
[clad_thermal]
type = HT9Thermal
block = clad
temperature = temp
[]
[clad_density]
type = StrainAdjustedDensity
block = clad
strain_free_density = 7874.0
[]
[]
[Preconditioning]
[SMP]
type = SMP
full = true
[]
[]
[Executioner]
type = Transient
solve_type = 'PJFNK'
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -ksp_gmres_restart'
petsc_options_value = 'lu superlu_dist 51'
line_search = 'none'
l_max_its = 60
l_tol = 8e-3
nl_max_its = 40
nl_rel_tol = 5e-4
nl_abs_tol = 1e-7
end_time = 42000400
dtmin = 100
dtmax = 5e5
[Quadrature]
order = fifth
side_order = seventh
[]
[TimeStepper]
type = IterationAdaptiveDT
timestep_limiting_postprocessor = creep_timestep
dt = 1e2
time_t = '0 1e5 41990400 42000400'
time_dt = '1e2 1e2 1e2 1e2'
iteration_window = 4
optimal_iterations = 10
[]
[]
[Postprocessors]
[ave_temp_interior] # average temperature of the cladding interior and all pellet exteriors
type = LayeredSideAverageValuePostprocessor
boundary = 9
variable = temp
execute_on = 'initial linear'
fuel_pin_geometry = pin_geometry
[]
[clad_inner_vol] # volume inside of cladding
type = LayeredInternalVolumePostprocessor
boundary = 7
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[pellet_volume] # fuel pellet total volume
type = LayeredInternalVolumePostprocessor
boundary = 8
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[avg_clad_temp] # average temperature of cladding interior
type = LayeredSideAverageValuePostprocessor
boundary = 7
variable = temp
fuel_pin_geometry = pin_geometry
execute_on = 'initial linear'
[]
[gas_volume]
type = LayeredInternalVolumePostprocessor
boundary = 9
execute_on = 'initial linear'
component = 0
out_of_plane_strain = strain_yy
fuel_pin_geometry = pin_geometry
addition = ${gas_addition}
[]
[flux_from_clad] # area integrated heat flux from the cladding
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 5
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[flux_from_fuel] # area integrated heat flux from the fuel
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 10
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[rod_total_power]
type = LayeredElementIntegralPowerPostprocessor
variable = temp
fission_rate = fission_rate
block = fuel
fuel_pin_geometry = pin_geometry
[]
[approx_FCT]
type = AverageNodalVariableValue
boundary = 12
variable = temp
[]
[max_approx_FCT]
type = TimeExtremeValue
value_type = max
postprocessor = approx_FCT
[]
[ave_FST]
type = SideAverageValue
boundary = 10
variable = temp
[]
[max_ave_FST]
type = TimeExtremeValue
value_type = max
postprocessor = ave_FST
[]
[ave_CIT]
type = SideAverageValue
boundary = 5
variable = temp
[]
[max_ave_CIT]
type = TimeExtremeValue
value_type = max
postprocessor = ave_CIT
[]
[peak_clad_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = clad
[]
[peak_fuel_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = fuel
[]
[max_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = max
block = fuel
[]
[min_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = min
block = fuel
[]
[peak_porosity]
type = ElementExtremeValue
variable = porosity
value_type = max
block = fuel
[]
[clad_fuel_gap]
type = NodalExtremeValue
variable = penetration
boundary = 10
[]
[max_cont_press]
type = NodalExtremeValue
variable = contact_pressure
boundary = 10
[]
[LHGR_W_per_cm]
type = FunctionValuePostprocessor
function = power_history
scale_factor = 0.01
[]
[average_burnup]
type = ElementAverageValue
block = fuel
variable = burnup
[]
[fis_gas_produced]
type = LayeredElementIntegralFisGasProducePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
[]
[fis_gas_released]
type = LayeredElementIntegralFisGasReleasePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
execute_on = 'initial timestep_end'
[]
[creep_timestep]
type = MaterialTimeStepPostprocessor
block = fuel
[]
[hydrostatic_stress]
type = ElementAverageValue
variable = hydrostatic_stress
execute_on = 'initial timestep_end'
block = fuel
[]
[solid_swelling]
type = ElementAverageValue
variable = solid_swell
block = fuel
[]
[gas_swelling]
type = ElementAverageValue
variable = gas_swell
block = fuel
[]
[volumetric_strain]
type = ElementAverageValue
variable = volumetric_strain
block = fuel
[]
[fission_rate]
type = ElementAverageValue
variable = fission_rate
block = fuel
[]
[porosity]
type = ElementAverageValue
variable = porosity
block = fuel
[]
[gaseous_porosity]
type = ElementAverageValue
variable = gaseous_porosity
block = fuel
[]
[fis_gas_percent]
type = FGRPercent
fission_gas_released = fis_gas_released
fission_gas_generated = fis_gas_produced
[]
[max_clad_hoop_creep]
type = ElementExtremeValue
value_type = max
block = clad
variable = hoop_creep_strain
[]
[max_total_hoop_strain]
type = ElementExtremeValue
value_type = max
block = clad
variable = total_hoop_strain
[]
[]
[VectorPostprocessors]
[clad_x_disp]
type = NodalValueSampler
variable = disp_x
boundary = 2
sort_by = y
outputs = 'vec1'
[]
[fuel_cl_temp]
type = NodalValueSampler
variable = temp
boundary = 12
sort_by = y
outputs = 'vec2'
[]
[fuel_surf_temp]
type = NodalValueSampler
variable = temp
boundary = 10
sort_by = y
outputs = 'vec3'
[]
[clad_inn_temp]
type = NodalValueSampler
variable = temp
boundary = 5
sort_by = y
outputs = 'vec4'
[]
[clad_out_temp]
type = NodalValueSampler
variable = temp
boundary = 2
sort_by = y
outputs = 'vec5'
[]
[fuel_surf_disp_x]
type = NodalValueSampler
variable = disp_x
boundary = 10
sort_by = y
outputs = 'vec7'
[]
[]
[PerformanceMetricOutputs]
outputs = performance_metrics_file
[]
[Outputs]
time_step_interval = 10
color = true
exodus = true
perf_graph = true
csv = true
sync_times = '1e3 5e3 1e4 5e4 1e5 5e6 1e6 5e6 1e7 2e7 3e7 4e7 41990400 41991000'
file_base = x441_${group_name}_1_5D
[out2]
type = CSV
file_base = x441_${group_name}_1_5D_out2
time_step_interval = 1
[]
[console]
type = Console
max_rows = 25
time_step_interval = 1
output_linear = true
[]
[chkfile]
type = CSV
file_base = x441_${group_name}_1_5D_chkfile
show = 'max_approx_FCT max_ave_FST max_ave_CIT average_burnup fis_gas_percent max_clad_hoop_creep max_total_hoop_strain'
execute_on = 'FINAL'
[]
[performance_metrics_file]
type = CSV
file_base = x441_${group_name}_1_5D_performance_metrics
show = 'simulation_alive_time number_linear_iterations number_nonlinear_iterations time_step_size total_linear_iterations total_nonlinear_iterations physical_memory_use number_dofs number_nonlinear_variables residual_compute_time jacobian_compute_time'
[]
[vec1]
type = CSV
file_base = x441_${group_name}_1_5D_vec1
execute_on = 'FINAL'
[]
[vec2]
type = CSV
file_base = x441_${group_name}_1_5D_vec2
execute_on = 'FINAL'
[]
[vec3]
type = CSV
file_base = x441_${group_name}_1_5D_vec3
execute_on = 'FINAL'
[]
[vec4]
type = CSV
file_base = x441_${group_name}_1_5D_vec4
execute_on = 'FINAL'
[]
[vec5]
type = CSV
file_base = x441_${group_name}_1_5D_vec5
execute_on = 'FINAL'
[]
[vec7]
type = CSV
file_base = x441_${group_name}_1_5D_vec7
execute_on = 'FINAL'
[]
[]
[Debug]
show_var_residual = 'disp_x temp'
show_var_residual_norms = true
[]
(test/tests/upuzr_fission_gas_release/ad/exact_1D.i)
# This tests the fission gas produced and released postprocessors
# and for 1d and 2d models with the same volume.
# The reported values should be the same for the
# 1d and 2d cases. See corresponding excel calculation
# in this directory.
[Problem]
solve = false
[]
[UserObjects]
[pin_geometry]
type = Layered1DFuelPinGeometry
include_clad = false
mesh_generator = layered1D_mesh
[]
[]
[Mesh]
coord_type = RZ
[layered1D_mesh]
type = Layered1DMeshGenerator
include_clad = false
include_plenum = false
slices_per_block = 1
fuel_height = 3e-3
pellet_outer_radius = 3e-3
pellet_bottom_coor = 0.0
pellet_mesh_density = customize
nx_p = 4
elem_type = EDGE2
[]
[]
[Functions]
[fission_rate_func]
type = PiecewiseLinear
x = '0 1e7'
y = '1e20 1e18'
[]
[porosity_func]
type = PiecewiseLinear
x = '0 1e7'
y = '0.15 0.25'
[]
[]
[Variables]
[temp]
initial_condition = 1700
[]
[]
[Kernels]
[heat]
type = ADTimeDerivative
variable = temp
[]
[]
[Materials]
[ad_properties]
type = ADGenericFunctionMaterial
prop_names = 'fission_rate'
prop_values = 'fission_rate_func'
outputs = all
[]
[properties]
type = ADGenericFunctionMaterial
prop_names = 'porosity'
prop_values = 'porosity_func'
outputs = all
[]
[fission_gas_release]
type = ADUPuZrFissionGasRelease
fission_rate = fission_rate
fractional_yield = 0.3017
critical_porosity = 0.2
use_old_porosity = false
time_average_fission_rate = false
[]
[]
[Executioner]
type = Transient
dt = 1e6
end_time = 1e7
[]
[Postprocessors]
[fission_rate]
type = ElementAverageValue
variable = fission_rate
[]
[gas_produced]
type = ADLayeredElementIntegralFisGasProducePostprocessor
fuel_pin_geometry = pin_geometry
[]
[gas_released]
type = ADLayeredElementIntegralFisGasReleasePostprocessor
fuel_pin_geometry = pin_geometry
[]
[porosity]
type = ElementAverageValue
variable = porosity
[]
[]
[Outputs]
csv = true
file_base = exact_out
[]
(assessment/metallic_fuel/EBRII/X441/analysis/group_C/x441_1_5D_C.i)
initial_fuel_density = 15800.0
[GlobalParams]
order = FIRST
family = LAGRANGE
energy_per_fission = 3.2e-11 # J/fission
displacements = 'disp_x'
[]
[Mesh]
coord_type = RZ
# Nominal Design Geometric Parameters (X441)
[layered1D_mesh]
type = Layered1DMeshGenerator
clad_thickness = ${clad_thickness}
pellet_outer_radius = ${pellet_outer_radius}
fuel_height = ${pellet_height}
plenum_height = ${clad_top_gap_height}
clad_gap_width = ${clad_gap_width}
# meshing parameters
clad_mesh_density = customize
pellet_mesh_density = customize
nx_p = 6
nx_c = 4
slices_per_block = 10
elem_type = EDGE2
[]
# mesh options
patch_update_strategy = auto
partitioner = centroid
centroid_partitioner_direction = y
[]
[UserObjects]
[pin_geometry]
type = Layered1DFuelPinGeometry
mesh_generator = layered1D_mesh
[]
[]
[Variables]
[temp]
initial_condition = 298
[]
[]
[AuxVariables]
[disp_y]
[]
[disp_z]
[]
# Aux variables for output
[gap_cond]
order = CONSTANT
family = MONOMIAL
[]
[coolant_htc]
order = CONSTANT
family = MONOMIAL
[]
[cumulative_damage_index]
order = CONSTANT
family = MONOMIAL
[]
[solid_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[gas_swell]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[volumetric_strain]
block = fuel
order = CONSTANT
family = MONOMIAL
[]
[hoop_stress]
order = CONSTANT
family = MONOMIAL
[]
[hoop_creep_strain]
order = CONSTANT
family = MONOMIAL
[]
[hoop_elastic_strain]
order = CONSTANT
family = MONOMIAL
[]
[total_hoop_strain]
order = CONSTANT
family = MONOMIAL
[]
[]
[Functions]
[power_history]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '0 44722 44722 0'
[]
[coolant_press_ramp]
type = PiecewiseLinear
x = '0 42000400'
y = '1.0 1.0'
[]
[coolant_temp_ramp]
type = PiecewiseLinear
x = '0 1e5 41990400 42000400'
y = '298.0 648.0 648.0 350.0'
[]
[axial_peaking_factors]
type = PowerPeakingFunction
fit = EBRII_ROW_4
pellet_length = ${pellet_height}
pellet_y_start = ${pellet_y_start}
[]
[clad_axial_pressure]
type = CladdingAxialPressureFunction
plenum_pressure = plenum_pressure
coolant_pressure = coolant_press_ramp
coolant_pressure_scaling_factor = 0.151e6
fuel_pin_geometry = pin_geometry
[]
[fuel_axial_pressure]
type = ParsedFunction
expression = plenum_pressure
symbol_names = plenum_pressure
symbol_values = plenum_pressure
[]
[]
[Physics/SolidMechanics/Layered1D]
[fuel]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = fuel_axial_pressure
generate_output = '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'
block = fuel
eigenstrain_names = 'fuel_thermal_strain gas_swelling_eigenstrain solid_swelling_eigenstrain'
mesh_generator = layered1D_mesh
[]
[clad]
strain = FINITE
add_variables = true
add_scalar_variables = true
out_of_plane_strain_name = strain_yy
fuel_pin_geometry = pin_geometry
out_of_plane_pressure_function = clad_axial_pressure
generate_output = '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'
block = clad
eigenstrain_names = 'clad_thermal_eigenstrain'
mesh_generator = layered1D_mesh
[]
[]
[Kernels]
# Define kernels for the various terms in the PDE system
[heat]
type = HeatConduction
variable = temp
[]
[heat_ie]
type = HeatConductionTimeDerivative
variable = temp
[]
[heat_source]
type = NeutronHeatSource
variable = temp
block = fuel
fission_rate = fission_rate
[]
[]
[AuxKernels]
[conductance]
type = MaterialRealAux
property = gap_conductance
variable = gap_cond
boundary = 10
[]
[gas_swell]
type = MaterialRealAux
variable = gas_swell
property = gas_swelling
execute_on = timestep_end
[]
[solid_swell]
type = MaterialRealAux
variable = solid_swell
property = solid_swelling
execute_on = timestep_end
[]
[volumetric_strain]
type = RankTwoScalarAux
rank_two_tensor = total_strain
variable = volumetric_strain
scalar_type = VolumetricStrain
execute_on = timestep_end
block = fuel
[]
[hoop_stress]
type = RankTwoAux
rank_two_tensor = stress
variable = hoop_stress
index_j = 2
index_i = 2
execute_on = timestep_end
[]
[hoop_creep_strain]
type = RankTwoAux
rank_two_tensor = creep_strain
variable = hoop_creep_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[hoop_elastic_strain]
type = RankTwoAux
rank_two_tensor = elastic_strain
variable = hoop_elastic_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[total_hoop_strain]
type = RankTwoAux
rank_two_tensor = total_strain
variable = total_hoop_strain
index_j = 2
index_i = 2
execute_on = timestep_end
block = clad
[]
[]
[Contact]
[pellet_clad_mechanical]
primary = 5
secondary = 10
penalty = 1e12
model = frictionless
formulation = kinematic
[]
[]
[ThermalContact]
[thermal_contact]
type = GapHeatTransfer
variable = temp
primary = 5
secondary = 10
quadrature = true
gap_conductivity = 61.0
min_gap = ${clad_gap_width}
[]
[]
[BCs]
[no_x_all]
type = DirichletBC
variable = disp_x
boundary = 12
value = 0.0
[]
[Pressure]
[coolantPressure]
boundary = '2'
function = coolant_press_ramp
[]
[]
[PlenumPressure]
[plenumPressure]
boundary = 9
initial_pressure = 0.084e6 # Pa
startup_time = 0
R = 8.3143
temperature = ave_temp_interior
volume = gas_volume
output = plenum_pressure
material_input = fis_gas_released
[]
[]
[]
[CoolantChannel]
[convective_clad_surface]
boundary = '2'
variable = temp
inlet_temperature = coolant_temp_ramp
inlet_pressure = coolant_press_ramp
inlet_massflux = 5261.5 # kg/m^2-sec
coolant_material = sodium
rod_diameter = 5.84e-3 # m
rod_pitch = 7.48e-3 # m (Pitch-to-diameter Ratio = 1.28)
linear_heat_rate = power_history
axial_power_profile = axial_peaking_factors
subchannel_geometry = triangular
[]
[]
[Materials]
[fission_rate]
type = UPuZrFissionRate
rod_linear_power = power_history
axial_power_profile = axial_peaking_factors
pellet_radius = ${pellet_outer_radius}
X_Zr = ${X_Zr}
X_Pu_function = ${X_Pu}
block = fuel
outputs = all
[]
[burnup]
type = UPuZrBurnup
initial_X_Zr = ${X_Zr}
initial_X_Pu = ${X_Pu}
density = ${initial_fuel_density}
block = fuel
outputs = all
[]
[fuel_elasticity_tensor]
type = UPuZrElasticityTensor
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
block = fuel
temperature = temp
[]
[fuel_inelastic_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'fuel_upuzrcreep'
block = fuel
[]
[fuel_upuzrcreep]
type = UPuZrCreepUpdate
block = fuel
temperature = temp
porosity = porosity
max_inelastic_increment = 1e-3
[]
[fuel_thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = fuel
thermal_expansion_coeff = 1.18e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = fuel_thermal_strain
[]
[gas_swelling]
type = UPuZrGaseousEigenstrain
temperature = temp
eigenstrain_name = gas_swelling_eigenstrain
initial_porosity = 0.0
bubble_number_density = 1e20
outputs = all
output_properties = 'porosity gaseous_porosity'
block = fuel
[]
[solid_swelling]
type = BurnupDependentEigenstrain
eigenstrain_name = solid_swelling_eigenstrain
block = fuel
swelling_name = 'solid_swelling'
[]
[metal_fuel_thermal]
type = UPuZrThermal
block = fuel
X_Zr = ${X_Zr}
X_Pu = ${X_Pu}
spheat_model = savage
thcond_model = lanl
porosity = porosity
temperature = temp
[]
[fuel_density]
type = StrainAdjustedDensity
block = fuel
strain_free_density = 15800
[]
[fission_gas_behavior]
type = UPuZrFissionGasRelease
block = fuel
fission_rate = fission_rate
[]
[clad_elasticity_tensor]
type = ComputeIsotropicElasticityTensor
youngs_modulus = 1.88e11
poissons_ratio = 0.236
block = clad
[]
[clad_stress]
type = ComputeMultipleInelasticStress
tangent_operator = nonlinear
inelastic_models = 'clad_ht9creep'
block = clad
[]
[fast_flux]
type = FastNeutronFlux
block = clad
factor = 2.47e19
[]
[clad_ht9creep]
type = HT9CreepUpdate
block = clad
temperature = temp
[]
[thermal_expansion]
type = ComputeThermalExpansionEigenstrain
block = clad
thermal_expansion_coeff = 1.2e-5
temperature = temp
stress_free_temperature = 295.0
eigenstrain_name = clad_thermal_eigenstrain
[]
[clad_thermal]
type = HT9Thermal
block = clad
temperature = temp
[]
[clad_density]
type = StrainAdjustedDensity
block = clad
strain_free_density = 7874.0
[]
[]
[Preconditioning]
[SMP]
type = SMP
full = true
[]
[]
[Executioner]
type = Transient
solve_type = 'PJFNK'
petsc_options = '-snes_ksp_ew'
petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -ksp_gmres_restart'
petsc_options_value = 'lu superlu_dist 51'
line_search = 'none'
l_max_its = 60
l_tol = 8e-3
nl_max_its = 40
nl_rel_tol = 5e-4
nl_abs_tol = 1e-7
end_time = 42000400
dtmin = 100
dtmax = 5e5
[Quadrature]
order = fifth
side_order = seventh
[]
[TimeStepper]
type = IterationAdaptiveDT
timestep_limiting_postprocessor = creep_timestep
dt = 1e2
time_t = '0 1e5 41990400 42000400'
time_dt = '1e2 1e2 1e2 1e2'
iteration_window = 4
optimal_iterations = 10
[]
[]
[Postprocessors]
[ave_temp_interior] # average temperature of the cladding interior and all pellet exteriors
type = LayeredSideAverageValuePostprocessor
boundary = 9
variable = temp
execute_on = 'initial linear'
fuel_pin_geometry = pin_geometry
[]
[clad_inner_vol] # volume inside of cladding
type = LayeredInternalVolumePostprocessor
boundary = 7
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[pellet_volume] # fuel pellet total volume
type = LayeredInternalVolumePostprocessor
boundary = 8
component = 0
fuel_pin_geometry = pin_geometry
out_of_plane_strain = strain_yy
execute_on = 'initial linear'
#outputs = exodus
[]
[avg_clad_temp] # average temperature of cladding interior
type = LayeredSideAverageValuePostprocessor
boundary = 7
variable = temp
fuel_pin_geometry = pin_geometry
execute_on = 'initial linear'
[]
[gas_volume]
type = LayeredInternalVolumePostprocessor
boundary = 9
execute_on = 'initial linear'
component = 0
out_of_plane_strain = strain_yy
fuel_pin_geometry = pin_geometry
addition = ${gas_addition}
[]
[flux_from_clad] # area integrated heat flux from the cladding
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 5
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[flux_from_fuel] # area integrated heat flux from the fuel
type = LayeredSideFluxIntegralPostprocessor
variable = temp
boundary = 10
diffusivity = thermal_conductivity
fuel_pin_geometry = pin_geometry
[]
[rod_total_power]
type = LayeredElementIntegralPowerPostprocessor
variable = temp
fission_rate = fission_rate
block = fuel
fuel_pin_geometry = pin_geometry
[]
[approx_FCT]
type = AverageNodalVariableValue
boundary = 12
variable = temp
[]
[max_approx_FCT]
type = TimeExtremeValue
value_type = max
postprocessor = approx_FCT
[]
[ave_FST]
type = SideAverageValue
boundary = 10
variable = temp
[]
[max_ave_FST]
type = TimeExtremeValue
value_type = max
postprocessor = ave_FST
[]
[ave_CIT]
type = SideAverageValue
boundary = 5
variable = temp
[]
[max_ave_CIT]
type = TimeExtremeValue
value_type = max
postprocessor = ave_CIT
[]
[peak_clad_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = clad
[]
[peak_fuel_temp]
type = ElementExtremeValue
variable = temp
value_type = max
block = fuel
[]
[max_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = max
block = fuel
[]
[min_hydro]
type = ElementExtremeValue
variable = hydrostatic_stress
value_type = min
block = fuel
[]
[peak_porosity]
type = ElementExtremeValue
variable = porosity
value_type = max
block = fuel
[]
[clad_fuel_gap]
type = NodalExtremeValue
variable = penetration
boundary = 10
[]
[max_cont_press]
type = NodalExtremeValue
variable = contact_pressure
boundary = 10
[]
[LHGR_W_per_cm]
type = FunctionValuePostprocessor
function = power_history
scale_factor = 0.01
[]
[average_burnup]
type = ElementAverageValue
block = fuel
variable = burnup
[]
[fis_gas_produced]
type = LayeredElementIntegralFisGasProducePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
[]
[fis_gas_released]
type = LayeredElementIntegralFisGasReleasePostprocessor
block = fuel
fuel_pin_geometry = pin_geometry
execute_on = 'initial timestep_end'
[]
[creep_timestep]
type = MaterialTimeStepPostprocessor
block = fuel
[]
[hydrostatic_stress]
type = ElementAverageValue
variable = hydrostatic_stress
execute_on = 'initial timestep_end'
block = fuel
[]
[solid_swelling]
type = ElementAverageValue
variable = solid_swell
block = fuel
[]
[gas_swelling]
type = ElementAverageValue
variable = gas_swell
block = fuel
[]
[volumetric_strain]
type = ElementAverageValue
variable = volumetric_strain
block = fuel
[]
[fission_rate]
type = ElementAverageValue
variable = fission_rate
block = fuel
[]
[porosity]
type = ElementAverageValue
variable = porosity
block = fuel
[]
[gaseous_porosity]
type = ElementAverageValue
variable = gaseous_porosity
block = fuel
[]
[fis_gas_percent]
type = FGRPercent
fission_gas_released = fis_gas_released
fission_gas_generated = fis_gas_produced
[]
[max_clad_hoop_creep]
type = ElementExtremeValue
value_type = max
block = clad
variable = hoop_creep_strain
[]
[max_total_hoop_strain]
type = ElementExtremeValue
value_type = max
block = clad
variable = total_hoop_strain
[]
[]
[VectorPostprocessors]
[clad_x_disp]
type = NodalValueSampler
variable = disp_x
boundary = 2
sort_by = y
outputs = 'vec1'
[]
[fuel_cl_temp]
type = NodalValueSampler
variable = temp
boundary = 12
sort_by = y
outputs = 'vec2'
[]
[fuel_surf_temp]
type = NodalValueSampler
variable = temp
boundary = 10
sort_by = y
outputs = 'vec3'
[]
[clad_inn_temp]
type = NodalValueSampler
variable = temp
boundary = 5
sort_by = y
outputs = 'vec4'
[]
[clad_out_temp]
type = NodalValueSampler
variable = temp
boundary = 2
sort_by = y
outputs = 'vec5'
[]
[fuel_surf_disp_x]
type = NodalValueSampler
variable = disp_x
boundary = 10
sort_by = y
outputs = 'vec7'
[]
[]
[PerformanceMetricOutputs]
outputs = performance_metrics_file
[]
[Outputs]
time_step_interval = 10
color = true
exodus = true
perf_graph = true
csv = true
sync_times = '1e3 5e3 1e4 5e4 1e5 5e6 1e6 5e6 1e7 2e7 3e7 4e7 41990400 41991000'
file_base = x441_${group_name}_1_5D
[out2]
type = CSV
file_base = x441_${group_name}_1_5D_out2
time_step_interval = 1
[]
[console]
type = Console
max_rows = 25
time_step_interval = 1
output_linear = true
[]
[chkfile]
type = CSV
file_base = x441_${group_name}_1_5D_chkfile
show = 'max_approx_FCT max_ave_FST max_ave_CIT average_burnup fis_gas_percent max_clad_hoop_creep max_total_hoop_strain'
execute_on = 'FINAL'
[]
[performance_metrics_file]
type = CSV
file_base = x441_${group_name}_1_5D_performance_metrics
show = 'simulation_alive_time number_linear_iterations number_nonlinear_iterations time_step_size total_linear_iterations total_nonlinear_iterations physical_memory_use number_dofs number_nonlinear_variables residual_compute_time jacobian_compute_time'
[]
[vec1]
type = CSV
file_base = x441_${group_name}_1_5D_vec1
execute_on = 'FINAL'
[]
[vec2]
type = CSV
file_base = x441_${group_name}_1_5D_vec2
execute_on = 'FINAL'
[]
[vec3]
type = CSV
file_base = x441_${group_name}_1_5D_vec3
execute_on = 'FINAL'
[]
[vec4]
type = CSV
file_base = x441_${group_name}_1_5D_vec4
execute_on = 'FINAL'
[]
[vec5]
type = CSV
file_base = x441_${group_name}_1_5D_vec5
execute_on = 'FINAL'
[]
[vec7]
type = CSV
file_base = x441_${group_name}_1_5D_vec7
execute_on = 'FINAL'
[]
[]
[Debug]
show_var_residual = 'disp_x temp'
show_var_residual_norms = true
[]