Material Property '...', requested by '...' is not defined on block 0

[mhelel]

I'm working with the Phase Field module and used the KKS formulation. After just copying the code from #https://github.com/idaholab/moose/blob/next/modules/phase_field/examples/kim-kim-suzuki/kks_example_noflux.i I wanted to add a Material Class that defines a small anisotropy:

C file:

registerMooseObject("otterApp", AnisotropicEpsilon);

InputParameters AnisotropicEpsilon::validParams() {
    InputParameters params = Material::validParams();

    params.addRequiredParam<Real>("alpha", "Anisotropy strength");
    params.addParam<Real>("eps0", 1.0, "Base epsilon");

    params.addCoupledVar("eta", "Order parameter to compute derivatives from");

    return params;
}

AnisotropicEpsilon::AnisotropicEpsilon(const InputParameters & parameters) : 
    Material(parameters), 
    _epsilon(declareProperty<Real>("epsilon")),
    _epsilon_sq(declareProperty<Real>("epsilon_sq")),
    _alpha(getParam<Real>("alpha")),
    _eps0(getParam<Real>("eps0")),
    _eta_grad(coupledGradient("eta")) 
    
    {}

void AnisotropicEpsilon::computeQpProperties() {
    RealVectorValue A(1.0, 0.0, 0.0);

    const Gradient & grad_eta = _eta_grad[_qp];

    Real dot = grad_eta * A;

    Real norm_grad = grad_eta.norm();
    Real norm_A = A.norm();

    if (norm_grad > 1e-12 && norm_A > 1e-12) {
        Real cos_theta = dot / (norm_grad * norm_A);
        cos_theta = std::max(Real(-1.0), std::min(Real(1.0), cos_theta));
        Real theta = std::acos(cos_theta);

        // For example, store epsilon depending on theta
        _epsilon[_qp] = _eps0 * (1.0 + _alpha * std::cos(4 * theta));
    } else {
        _epsilon[_qp] = _eps0;
    }

    _epsilon_sq[_qp] = _epsilon[_qp] * _epsilon[_qp];


}

H-File:

#include "Material.h"

class AnisotropicEpsilon : public Material {
    public: 
        AnisotropicEpsilon(const InputParameters & parameters);

        static InputParameters validParams();
    
    protected:
        virtual void computeQpProperties() override;

        /// member variable to hold the computed epsilon
        MaterialProperty<Real> & _epsilon;
        MaterialProperty<Real> & _epsilon_sq;

        /// A place to hold the coupled variable gradient for calculating the convection velocity property
        const VariableGradient & _eta_grad;

        Real _alpha;
        Real _eps0;

};

With the input file beeing:

[Mesh]
  [gen]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 10
    ny = 10
    nz = 0
    xmin = -10
    xmax = 10
    ymin = -10
    ymax = 10
    zmin = 0
    zmax = 0
    elem_type = QUAD4
  []
  [./Partitioner]
    type = LibmeshPartitioner
    partitioner = PARMETIS
  [../]
[]

[AuxVariables]
  [./Fglobal]
    order = CONSTANT
    family = MONOMIAL
  [../]
[]

[Variables]

  [./eta] # order parameter
    order = FIRST
    family = LAGRANGE
  [../]

  [./c] # solute concentration
    order = FIRST
    family = LAGRANGE
  [../]

  [./w] # chemical potential
    order = FIRST
    family = LAGRANGE
  [../]

  [./cl] # Liquid phase solute concentration
    order = FIRST
    family = LAGRANGE
    initial_condition = 0.1
  [../]
  
  [./cs] # Solid phase solute concentration
    order = FIRST
    family = LAGRANGE
    initial_condition = 0.9
  [../]
[]

[Functions]
  [./ic_func_eta]
    type = ParsedFunction
    expression = '0.5*(1.0-tanh((sqrt(x^2+y^2)-3)/sqrt(2.0)))'
  [../]
  [./ic_func_c]
    type = ParsedFunction
    expression = '0.6*0.5*(1.0-tanh((sqrt(x^2+y^2)-3)/sqrt(2.0)))+0.3'
  [../]
[]

[ICs]
  [./eta]
    variable = eta
    type = FunctionIC
    function = ic_func_eta
  [../]
  [./c]
    variable = c
    type = FunctionIC
    function = ic_func_c
  [../]
[]

[Materials]
  # Free energy of the liquid
  [./fl]
    type = DerivativeParsedMaterial
    property_name = fl
    coupled_variables = 'cl'
    expression = '(0.1-cl)^2'
  [../]

  # Free energy of the solid
  [./fs]
    type = DerivativeParsedMaterial
    property_name = fs
    coupled_variables = 'cs'
    expression = '(0.9-cs)^2'
  [../]

  # h(eta)
  [./h_eta]
    type = SwitchingFunctionMaterial
    h_order = HIGH
    eta = eta
  [../]

  # g(eta)
  [./g_eta]
    type = BarrierFunctionMaterial
    g_order = SIMPLE
    eta = eta
  [../]

  [./anisotropic_eps]
    type = AnisotropicEpsilon
    eta = eta
    alpha = 0.05
    eps0 = 1.0
  [../]

  # constant properties
  [./constants]
    type = GenericConstantMaterial
    prop_names  = 'M   L  '
    prop_values = '0.7 0.7'
  [../]

[]


[Kernels]
  active = 'PhaseConc ChemPotSolute CHBulk ACBulkF ACBulkC ACInterface dcdt detadt ckernel'

  # enforce c = (1-h(eta))*cl + h(eta)*cs
  [./PhaseConc]
    type = KKSPhaseConcentration
    ca       = cl
    variable = cs
    c        = c
    eta      = eta
  [../]

  # enforce pointwise equality of chemical potentials
  [./ChemPotSolute]
    type = KKSPhaseChemicalPotential
    variable = cl
    cb       = cs
    fa_name  = fl
    fb_name  = fs
  [../]

  #
  # Cahn-Hilliard Equation
  #
  [./CHBulk]
    type = KKSSplitCHCRes
    variable = c
    ca       = cl
    fa_name  = fl
    w        = w
  [../]

  [./dcdt]
    type = CoupledTimeDerivative
    variable = w
    v = c
  [../]
  [./ckernel]
    type = SplitCHWRes
    mob_name = M
    variable = w
  [../]

  #
  # Allen-Cahn Equation
  #
  [./ACBulkF]
    type = KKSACBulkF
    variable = eta
    fa_name  = fl
    fb_name  = fs
    w        = 1.0
    coupled_variables = 'cl cs'
  [../]
  [./ACBulkC]
    type = KKSACBulkC
    variable = eta
    ca       = cl
    cb       = cs
    fa_name  = fl
  [../]

  [./ACInterface]
    type = ACInterface
    variable = eta
    kappa_name = anisotropic_eps.epsilon_sq
  [../]
  
  [./detadt]
    type = TimeDerivative
    variable = eta
  [../]
[]

[AuxKernels]
  [./GlobalFreeEnergy]
    variable = Fglobal
    type = KKSGlobalFreeEnergy
    fa_name = fl
    fb_name = fs
    w = 1.0
  [../]
[]

[Executioner]
  type = Transient
  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -sub_pc_type -sub_pc_factor_shift_type'
  petsc_options_value = 'asm      ilu          nonzero'

  l_max_its = 100
  nl_max_its = 100

  num_steps = 50
  dt = 0.1
[]


# Precondition using handcoded off-diagonal terms
[Preconditioning]
  [./full]
    type = SMP
    full = true
  [../]
[]


[VectorPostprocessors]
  [./c]
    type =  LineValueSampler
    start_point = '-10 0 0'
    end_point = '10 0 0'
    variable = c
    num_points = 151
    sort_by =  id
    execute_on = timestep_end
  [../]
  [./eta]
    type =  LineValueSampler
    start_point = '-10 0 0'
    end_point = '10 0 0'
    variable = eta
    num_points = 151
    sort_by =  id
    execute_on = timestep_end
  [../]
[]

[Adaptivity]
  max_h_level = 3
  interval = 1
  marker = combo
  initial_marker = combo
  initial_steps = 3

  [./Markers]
    [./combo]
      type = ComboMarker
      markers = 'mark_eta_grad'
    [../]

    [./value_marker]
      type = ValueRangeMarker
      lower_bound = 0.1   # adjust for eta ∈ [0,1]
      upper_bound = 0.9
      buffer_size = 0.1
      variable = eta
    [../]

    [./mark_eta_grad]
      type = ErrorFractionMarker
      coarsen = 0.2
      refine = 0.6
      indicator = ind_eta_grad
    [../]
  [../]

  [./Indicators]
    [./ind_eta_grad]
      type = GradientJumpIndicator
      variable = eta
    [../]
  [../]
[]


[Outputs]
  exodus = false
  [./Soln]
    type = Exodus
    file_base = 'problems/ACCHAniso/ACCH_out'
    execute_on = 'INITIAL TIMESTEP_END FINAL'
    time_step_interval = 5
  [../]
[]

But now I get the error:

*** ERROR ***
The following error occurred in the Problem 'MOOSE Problem' of type FEProblem.

Material property 'anisotropic_eps.epsilon_sq', requested by 'ACInterface' is not defined on block 0

I don't have any idea, why this would not work. I already tried deleting the [./Partitioner] Block. I also tried adding block = 0 in the [./anisotropic_eps] block. But it still did not work

[GiudGiud]

Hello

you don’t need to prefix property names with the material name.

kappa_name = anisotropic_eps.epsilon_sq
could just be
kappa_name = epsilon_sq

It is fairly common to add prefixes to material property but those need to be a parameter to the material

[mhelel]

Thanks this solved it