Mixed Oxide Fuel Tutorial

[UserObjects<<<{"href": "../../syntax/UserObjects/index.html"}>>>]
  [pin_geometry]
    type = Layered1DFuelPinGeometry<<<{"description": "Computes LWR fuel pin geometry for 1D meshes by reading the input mesh. This object can be coupled to Burnup and other functions as an alternative to having the user supply parameters such as pellet radius and pellet-cladding gap.", "href": "../../source/userobject/Layered1DFuelPinGeometry.html"}>>>
     include_clad<<<{"description": "Whether to include the clad block"}>>> = false
    mesh_generator<<<{"description": "The name of the generator to use as the prefix for mesh meta data properties."}>>> = layered1D_mesh
  []
[]

[Mesh<<<{"href": "../../syntax/Mesh/index.html"}>>>]
  coord_type = RZ
  [layered1D_mesh]
    type = Layered1DMeshGenerator<<<{"description": "Creates an axisymmetric mesh composed of layers of 1-dimensional elements.", "href": "../../source/meshgenerators/Layered1DMeshGenerator.html"}>>>
    fuel_height<<<{"description": "Height of the fuel."}>>> = 0.1
    pellet_outer_radius<<<{"description": "Pellet outer radius.  If more than one given, number must match number of fuel blocks."}>>> = 0.00541
    include_clad<<<{"description": "Whether to include the clad block."}>>> = false
    pellet_bottom_coor<<<{"description": "Axial location of the start of the fuel stack."}>>> = 0.0
    pellet_mesh_density<<<{"description": "Sets the mesh density of the pellet (really_coarse, coarse, medium, fine, customize)."}>>> = customize
    nx_p<<<{"description": "Number of fuel elements in radial direction."}>>> = 200
    elem_type<<<{"description": "The type of element from libMesh to generate."}>>> = EDGE2
    slices_per_block<<<{"description": "For cases with only cladding, this is a single number providing the total number of slices for the cladding.  For cases with fuel, this is the number of slices per fuel block. Multiple values are provided if multiple fuel blocks are used."}>>> = 1
    include_plenum<<<{"description": "Whether to include the plenum."}>>> = false
  []
[]

[Variables<<<{"href": "../../syntax/Variables/index.html"}>>>]
  [temp]
    initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = 300
  []
  [pore]
    initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = 0.15
    scaling<<<{"description": "Specifies a scaling factor to apply to this variable"}>>> = 1e14
  []
[]

[AuxVariables<<<{"href": "../../syntax/AuxVariables/index.html"}>>>]
 [pore_speed_aux]
   order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = constant
   family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = monomial
 []
 [fission_rate_aux_variable]
   order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = first
   family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = lagrange
 []
 [fission_rate_aux_variable_mox]
   order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = first
   family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = lagrange
 []
 [grad_temp_x]
   order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = CONSTANT
   family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = MONOMIAL
 []
 [thermal_conductivity]
   order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = CONSTANT
   family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = MONOMIAL
   block = fuel
 []
[]

[Functions<<<{"href": "../../syntax/Functions/index.html"}>>>]
  [power_history1]
    type = PiecewiseLinear<<<{"description": "Linearly interpolates between pairs of x-y data", "href": "../../source/functions/PiecewiseLinear.html"}>>>
      x<<<{"description": "The abscissa values"}>>> = '0 10000'
      y<<<{"description": "The ordinate values"}>>> = '0 50000'
  []
  [fuel_surface_temp]
    type = PiecewiseLinear<<<{"description": "Linearly interpolates between pairs of x-y data", "href": "../../source/functions/PiecewiseLinear.html"}>>>
    data_file<<<{"description": "File holding CSV data"}>>> = fuel_surface_temp_bc.csv
    scale_factor<<<{"description": "Scale factor to be applied to the ordinate values"}>>> = 1
    format<<<{"description": "Format of csv data file that is in either in columns or rows"}>>> = columns
  []
[]

[Kernels<<<{"href": "../../syntax/Kernels/index.html"}>>>]

  [heat]         # gradient term in heat conduction equation
    type = HeatConduction<<<{"description": "Diffusive heat conduction term $-\\nabla\\cdot(k\\nabla T)$ of the thermal energy conservation equation", "href": "../../source/kernels/HeatConduction.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temp
  []
  [heat_ie]       # time term in heat conduction equation
    type = HeatConductionTimeDerivative<<<{"description": "Time derivative term $\\rho c_p \\frac{\\partial T}{\\partial t}$ of the thermal energy conservation equation.", "href": "../../source/kernels/HeatConductionTimeDerivative.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temp
  []
  [heat_source]  # source term in heat conduction equation
    type = NeutronHeatSource<<<{"description": "Compute heat generation due to fission.", "href": "../../source/kernels/NeutronHeatSource.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temp
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel # fission rate applied to the fuel (block 2) only
    fission_rate<<<{"description": "Coupled Fission Rate"}>>> = fission_rate_aux_variable_mox
  []

  [pore_continuity]
   type = MOXPoreContinuity<<<{"description": "MOX kernel used to simulate pore migration.", "href": "../../source/kernels/MOXPoreContinuity.html"}>>>
   variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = pore
   temperature<<<{"description": "Temperature."}>>> = temp
   debug<<<{"description": "Whether or not to debug"}>>> = 0
   alpha<<<{"description": "alpha"}>>> = 0.25
   beta<<<{"description": "beta"}>>> = 1
   heating_function<<<{"description": "function for ralp"}>>> = power_history1
  []
  [pore_diffusion]
   type = MOXPoreDiffusion<<<{"description": "MOX porosity diffusion kernel used with kernel MOXPoreContinuity.", "href": "../../source/kernels/MOXPoreDiffusion.html"}>>>
   variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = pore
   debug<<<{"description": "Whether or not to debug"}>>> = 0
#   nu = 3.25e-8 #seems to be THE value to use... result is super sensitive to this number
#   nu = 10e-10
   nu<<<{"description": "porosity diffusion coefficient"}>>> = 1e-12
   heating_function<<<{"description": "function for ralp"}>>> = power_history1
   v_upper<<<{"description": "v_upper"}>>> = 1e-12
   v_lower<<<{"description": "v_lower"}>>> = 1e-20
  # v_upper = 1
  # v_lower = 1
  []
  [poretimederivative]
   type = CoefTimeDerivative<<<{"description": "The time derivative operator with the weak form of $(\\psi_i, \\frac{\\partial u_h}{\\partial t})$.", "href": "../../source/kernels/CoefTimeDerivative.html"}>>>
   variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = pore
   Coefficient<<<{"description": "The coefficient for the time derivative kernel"}>>> = 1
   []
[]

[AuxKernels<<<{"href": "../../syntax/AuxKernels/index.html"}>>>]
  [pore_speed_aux]
    type = MaterialRealAux<<<{"description": "Outputs element volume-averaged material properties", "href": "../../source/auxkernels/MaterialRealAux.html"}>>>
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = pore_speed_aux
    property<<<{"description": "The material property name."}>>> = pore_velocity
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    execute_on<<<{"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."}>>> = 'initial timestep_end'
  []
  [fission_rate_aux_kernel]
    type = FissionRateGeneral<<<{"description": "Provides various implementations which allow the user to set the fission rate.", "href": "../../source/auxkernels/FissionRateGeneral.html"}>>>
    fission_rate_formulation<<<{"description": "The active fission rate formulation."}>>> = LWR
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = fission_rate_aux_variable
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    rod_ave_lin_pow<<<{"description": "The rod power history function."}>>> = power_history1
#    axial_power_profile = axial_peaking_factors     # using the axial power profile function defined above
    pellet_diameter<<<{"description": "Fuel pellet diameter in meters."}>>> = 0.01082
    pellet_inner_diameter<<<{"description": "Pellet inner diameter in meters for an annular pellet."}>>> = 0
    energy_per_fission<<<{"description": "Energy in Joules per fission."}>>> = 3.2e-11
    execute_on<<<{"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."}>>> = 'initial timestep_end'
    value<<<{"description": "If no function is provided, the value of the fission rate, otherwise, a scaling factor for the function."}>>> = 1.0
  []
  [fission_rate_aux_kernel_mox]
    type = FissionRateGeneral<<<{"description": "Provides various implementations which allow the user to set the fission rate.", "href": "../../source/auxkernels/FissionRateGeneral.html"}>>>
    fission_rate_formulation<<<{"description": "The active fission rate formulation."}>>> = MOX
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = fission_rate_aux_variable_mox
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    porosity<<<{"description": "Coupled porosity"}>>> = pore
    initial_porosity<<<{"description": "Initial porosity. Must be the same as the initial condition in the variable block."}>>> = 0.15
    rod_ave_lin_pow<<<{"description": "The rod power history function."}>>> = power_history1
    pellet_diameter<<<{"description": "Fuel pellet diameter in meters."}>>> = 0.01082
    pellet_inner_diameter<<<{"description": "Pellet inner diameter in meters for an annular pellet."}>>> = 0
    energy_per_fission<<<{"description": "Energy in Joules per fission."}>>> = 3.2e-11
    execute_on<<<{"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."}>>> = 'initial timestep_end'
    value<<<{"description": "If no function is provided, the value of the fission rate, otherwise, a scaling factor for the function."}>>> = 1.0
  []
  [grad_temp_x_aux]
    type = VariableGradientComponent<<<{"description": "Creates a field consisting of one component of the gradient of a coupled variable.", "href": "../../source/auxkernels/VariableGradientComponent.html"}>>>
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = grad_temp_x
    component<<<{"description": "The gradient component to compute"}>>> = x
    gradient_variable<<<{"description": "The variable from which to compute the gradient component"}>>> = temp
    execute_on<<<{"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."}>>> = 'initial timestep_end'
  []
  [ThermalConductivityAux]
    type = MaterialRealAux<<<{"description": "Outputs element volume-averaged material properties", "href": "../../source/auxkernels/MaterialRealAux.html"}>>>
     block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
     execute_on<<<{"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."}>>> = linear
     property<<<{"description": "The material property name."}>>> = thermal_conductivity
     variable<<<{"description": "The name of the variable that this object applies to"}>>> = thermal_conductivity
   []
[]

[BCs<<<{"href": "../../syntax/NuclearMaterials/BCs/index.html"}>>>]
[temp_outside] # pin pellets and clad along axis of symmetry (y)
  type = FunctionDirichletBC<<<{"description": "Imposes the essential boundary condition $u=g(t,\\vec{x})$, where $g$ is a (possibly) time and space-dependent MOOSE Function.", "href": "../../source/bcs/FunctionDirichletBC.html"}>>>
  variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temp
  boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = 10
  function<<<{"description": "The forcing function."}>>> = fuel_surface_temp
[]
[]

[Materials<<<{"href": "../../syntax/Materials/index.html"}>>>]
  [fuel_thermal]
    type = MAMOXThermal<<<{"description": "Computes the thermal conductivity for minor actinide fast MOX fuel.", "href": "../../source/materials/MAMOXThermal.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    temperature<<<{"description": "Coupled Temperature"}>>> = temp
    porosity<<<{"description": "Coupled Porosity from field variable"}>>> = pore
    porosity_limit<<<{"description": "Maximum allowed porosity for thermal conductivity calculation. Avoids k = 0 when porosity = 1"}>>> = 0.95
  []
  [density_block]
    type = GenericConstantMaterial<<<{"description": "Declares material properties based on names and values prescribed by input parameters.", "href": "../../source/materials/GenericConstantMaterial.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    prop_names<<<{"description": "The names of the properties this material will have"}>>> = density
    prop_values<<<{"description": "The values associated with the named properties"}>>> = 10662.0
  []
  [pore_velocity]
   type = MOXPoreVelocity<<<{"description": "Computes pore speed. Used with kernel [MOXPoreContinuity](/MOXPoreContinuity).", "href": "../../source/materials/MOXPoreVelocity.html"}>>>
   block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
   temperature<<<{"description": "Coupled Temperature"}>>> = temp
   limit<<<{"description": "limit for pore velocity in m/s"}>>> = 1e-3
#go back to this if necessary   scale_factor = 0.05
   scale_factor<<<{"description": "scale the velocity to account for uncertainty"}>>> = 0.1
  []
[]

[Preconditioning<<<{"href": "../../syntax/Preconditioning/index.html"}>>>]
  [SMP]
    type = SMP<<<{"description": "Single matrix preconditioner (SMP) builds a preconditioner using user defined off-diagonal parts of the Jacobian.", "href": "../../source/preconditioners/SingleMatrixPreconditioner.html"}>>>
    full<<<{"description": "Set to true if you want the full set of couplings between variables simply for convenience so you don't have to set every off_diag_row and off_diag_column combination."}>>> = true
  []
[]

[Dampers<<<{"href": "../../syntax/Dampers/index.html"}>>>]
  [limitT]
    type = MaxIncrement<<<{"description": "Limits a variable's update by some max fraction", "href": "../../source/dampers/MaxIncrement.html"}>>>
    max_increment<<<{"description": "The maximum newton increment for the variable."}>>> = 100.0
    variable<<<{"description": "The name of the variable that this damper operates on"}>>> = temp
  []
[]

[Executioner<<<{"href": "../../syntax/Executioner/index.html"}>>>]
  type = Transient
  solve_type = 'PJFNK'

  petsc_options = '-snes_ksp_ew'
  #petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'  # -mat_superlu_dist_fact'
  #petsc_options_value = 'lu superlu_dist' # SamePattern_SameRowPerm'
  petsc_options_iname = '-pc_type -sub_pc_type'
  petsc_options_value = 'asm lu'
  line_search = 'none'

  l_max_its = 50
  l_tol = 8e-3
  nl_max_its = 50
  nl_rel_tol = 1e-5
  nl_abs_tol = 1e-6
  end_time = 10000
  dtmax = 100
  dtmin = 0.25

  [TimeStepper<<<{"href": "../../syntax/Executioner/TimeStepper/index.html"}>>>]
    type = IterationAdaptiveDT
    dt = 1e1
    optimal_iterations = 15
    iteration_window = 2
    linear_iteration_ratio = 100
    growth_factor = 2
    cutback_factor = .5
    force_step_every_function_point = true
    timestep_limiting_function = power_history1
  []
[]

[Postprocessors<<<{"href": "../../syntax/Postprocessors/index.html"}>>>]
  [_dt]                     # time step
    type = TimestepSize<<<{"description": "Reports the timestep size", "href": "../../source/postprocessors/TimestepSize.html"}>>>
  []
  [z_nonlinear_its]           # number of nonlinear iterations at each timestep
    type = NumNonlinearIterations<<<{"description": "Outputs the number of nonlinear iterations", "href": "../../source/postprocessors/NumNonlinearIterations.html"}>>>
  []
  [a_run_time]               # average temperature of cladding interior
    type = PerfGraphData<<<{"description": "Retrieves performance information about a section from the PerfGraph.", "href": "../../source/postprocessors/PerfGraphData.html"}>>>
    section_name<<<{"description": "The name of the section to get data for"}>>> = Root
    data_type<<<{"description": "The type of data to retrieve for the section_name"}>>> = TOTAL
  []

  [ave_fuel_temp]
    type = ElementAverageValue<<<{"description": "Computes the volumetric average of a variable", "href": "../../source/postprocessors/ElementAverageValue.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    variable<<<{"description": "The name of the variable that this object operates on"}>>> = temp
  []
  [max_fuel_temp]
    type = NodalExtremeValue<<<{"description": "Finds either the min or max elemental value of a variable over the domain.", "href": "../../source/postprocessors/NodalExtremeValue.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    value_type<<<{"description": "Type of extreme value to return. 'max' returns the maximum value. 'min' returns the minimum value. 'max_abs' returns the maximum of the absolute value."}>>> = max
    variable<<<{"description": "The name of the variable that this postprocessor operates on"}>>> = temp
  []

  [ave_pore]
    type = ElementAverageValue<<<{"description": "Computes the volumetric average of a variable", "href": "../../source/postprocessors/ElementAverageValue.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    variable<<<{"description": "The name of the variable that this object operates on"}>>> = pore
  []
  [max_pore]
    type = NodalExtremeValue<<<{"description": "Finds either the min or max elemental value of a variable over the domain.", "href": "../../source/postprocessors/NodalExtremeValue.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    value_type<<<{"description": "Type of extreme value to return. 'max' returns the maximum value. 'min' returns the minimum value. 'max_abs' returns the maximum of the absolute value."}>>> = max
    variable<<<{"description": "The name of the variable that this postprocessor operates on"}>>> = pore
  []
  [min_pore]
    type = NodalExtremeValue<<<{"description": "Finds either the min or max elemental value of a variable over the domain.", "href": "../../source/postprocessors/NodalExtremeValue.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    value_type<<<{"description": "Type of extreme value to return. 'max' returns the maximum value. 'min' returns the minimum value. 'max_abs' returns the maximum of the absolute value."}>>> = min
    variable<<<{"description": "The name of the variable that this postprocessor operates on"}>>> = pore
  []
  [max_pore_speed]
    type = ElementExtremeValue<<<{"description": "Finds either the min or max elemental value of a variable over the domain.", "href": "../../source/postprocessors/ElementExtremeValue.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    value_type<<<{"description": "Type of extreme value to return. 'max' returns the maximum value. 'min' returns the minimum value. 'max_abs' returns the maximum of the absolute value."}>>> = max
    variable<<<{"description": "The name of the variable that this postprocessor operates on"}>>> = pore_speed_aux
  []
  [ave_fission_rate]
    type = ElementAverageValue<<<{"description": "Computes the volumetric average of a variable", "href": "../../source/postprocessors/ElementAverageValue.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    variable<<<{"description": "The name of the variable that this object operates on"}>>> = fission_rate_aux_variable
  []
  [rod_total_power]
    type = LayeredElementIntegralPowerPostprocessor<<<{"description": "Computes the power given the fission rate and energy per fission for Layered1D and Layered2D geometeries.", "href": "../../source/postprocessors/LayeredElementIntegralPowerPostprocessor.html"}>>>
    variable<<<{"description": "The name of the variable that this object operates on"}>>> = temp
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    fission_rate<<<{"description": "Coupled fission rate"}>>> = fission_rate_aux_variable
    fuel_pin_geometry<<<{"description": "Name of fuel pin geometry UserObject which is required for Layered1D or Layered2D calculations."}>>> = pin_geometry
  []
  [rod_total_power_mox]
    type = LayeredElementIntegralPowerPostprocessor<<<{"description": "Computes the power given the fission rate and energy per fission for Layered1D and Layered2D geometeries.", "href": "../../source/postprocessors/LayeredElementIntegralPowerPostprocessor.html"}>>>
    variable<<<{"description": "The name of the variable that this object operates on"}>>> = temp
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    fission_rate<<<{"description": "Coupled fission rate"}>>> = fission_rate_aux_variable_mox
    fuel_pin_geometry<<<{"description": "Name of fuel pin geometry UserObject which is required for Layered1D or Layered2D calculations."}>>> = pin_geometry
  []
  [rod_input_power]
    type = FunctionValuePostprocessor<<<{"description": "Computes the value of a supplied function at a single point (scalable)", "href": "../../source/postprocessors/FunctionValuePostprocessor.html"}>>>
    function<<<{"description": "The function which supplies the postprocessor value."}>>> = power_history1
    scale_factor<<<{"description": "A scale factor to be applied to the function"}>>> = 0.1 # rod height
  []
  [ave_themal_conductivity]
    type = ElementAverageValue<<<{"description": "Computes the volumetric average of a variable", "href": "../../source/postprocessors/ElementAverageValue.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    variable<<<{"description": "The name of the variable that this object operates on"}>>> = thermal_conductivity
  [] # end element average burnup
[]

#[VectorPostprocessors]
#  [line_value_vector_postprocessor_pore]
#    type = LineValueSampler
#    variable = pore
#    start_point = '0.0 0.05 0'
#    end_point = '0.0041 0.05 0'
#    num_points = 100
#    sort_by = x
#    execute_on = linear
#    outputs = stuff_v_rad
#    control_tags = a
#  []
#  [line_value_vector_postprocessor_gradT]
#    type = LineValueSampler
#    variable = grad_temp_x
#    start_point = '0.0 0.05 0'
#    end_point = '0.0041 0.05 0'
#    num_points = 100
#    sort_by = x
#    execute_on = linear
#    outputs = stuff_v_rad
#  []
#  [line_value_vector_postprocessor_pore_speed]
#    type = LineValueSampler
#    variable = pore_speed_aux
#    start_point = '0.0 0.05 0'
#    end_point = '0.0041 0.05 0'
#    num_points = 100
#    sort_by = x
#    execute_on = linear
#    outputs = stuff_v_rad
#  []
#  [line_value_vector_postprocessor_temp]
#    type = LineValueSampler
#    variable = temp
#    start_point = '0.0 0.05 0'
#    end_point = '0.0041 0.05 0'
#    num_points = 100
#    sort_by = x
#    execute_on = linear
#    outputs = stuff_v_rad
#  []
#[]

[Outputs<<<{"href": "../../syntax/Outputs/index.html"}>>>]
  perf_graph<<<{"description": "Enable printing of the performance graph to the screen (Console)"}>>> = true
  exodus<<<{"description": "Output the results using the default settings for Exodus output."}>>> = true
  csv<<<{"description": "Output the scalar variable and postprocessors to a *.csv file using the default CSV output."}>>> = false
  color<<<{"description": "Set to false to turn off all coloring in all outputs"}>>> = false
  [console]
    type = Console<<<{"description": "Object for screen output.", "href": "../../source/outputs/Console.html"}>>>
    max_rows<<<{"description": "The maximum number of postprocessor/scalar values displayed on screen during a timestep (set to 0 for unlimited)"}>>> = 25
    all_variable_norms<<<{"description": "If true, all variable norms will be printed after each solve"}>>> = true
  []
  #[stuff_v_rad]
  #  type = CSV
  #  execute_on = 'FINAL'
  #[]
[]

[Debug<<<{"href": "../../syntax/Debug/index.html"}>>>]
 show_var_residual_norms<<<{"description": "Print the residual norms of the individual solution variables at each nonlinear iteration"}>>> = true
[]