fv_material_quadrature.i

#
# Parsed material properties depend on the physical location of the element
# This requires the initialization of the quadrature in the FVFlux loop

[Mesh]
  type = GeneratedMesh
  dim = 2
  xmin = 0
  xmax = 1
  ymin = 0
  ymax = 1
  nx = 4
  ny = 4
  elem_type = QUAD9
[]

[Functions]
  [linear_x]
    type = ParsedFunction
    value = 'x'
  []
[]

[Variables]
  [u]
    family = MONOMIAL
    order = CONSTANT
    fv = true
  []
[]

[FVKernels]
  [diff]
    type = FVDiffusion
    variable = u
    coeff = k1
  []

  [r]
    type = FVReaction
    variable = u
  []
[]

[FVBCs]
  [all]
    type = FVDirichletBC
    variable = u
    boundary = 'left right bottom top'
    value = 1
  []
[]

[Materials]
  [k1]
    type = ADGenericFunctionMaterial
    prop_names = 'k1'
    prop_values = linear_x
    block = 0
  []
[]

[Executioner]
  type = Steady
  solve_type = NEWTON

  petsc_options_iname = '-pc_type'
  petsc_options_value = 'lu'
[]

[Problem]
  kernel_coverage_check = false
[]

[Outputs]
  execute_on = 'timestep_end'
  exodus = true
[]