Problems reading material property in InterfaceKernel #28908

bmurgas-LANL · 2024-10-22T21:56:32Z

bmurgas-LANL
Oct 22, 2024

Check these boxes if you have followed the posting rules.

Q&A General is the most appropriate section for my question
I have consulted the posting Guidelines on the Discussions front page
I have searched the Discussions forum and my question has not been asked before
I have searched the MOOSE website and the documentation does not answer my question
I have formatted my post following the posting guidelines (screenshots as a last resort, triple back quotes around pasted text)

Question

Hello Moose community,

I am using a crystal plasticity model where I define material properties that I need to use inside a modified version of ConservativeAdvection kernel and InterfaceDiffusion interface kernel.
When I read the property in the ConservativeAdvection kernel, the material properties are correct and have the same value computed in the crystal plasticity model.
However, when I read the material properties in the InterfaceDiffusion interface kernel the variables are set to zero. Do you know why this can happen?

This is the validParameters and constructor of the "ConservativeDiffusion" kernel:

InputParameters
ConservativeAdvectionDislocation::validParams()
{
  InputParameters params = Kernel::validParams();
  params.addClassDescription(
      "Conservative form of $\\nabla \\cdot \\vec{v} u$ which in its weak "
      "form is given by: $(-\\nabla \\psi_i, \\vec{v} u)$. "
      "Velocity is defined as a material property computed in DiscoFluxCPStressUpdate");
  MooseEnum upwinding_type("none full", "none");
  params.addParam<MooseEnum>("upwinding_type",
                             upwinding_type,
                             "Type of upwinding used.  None: Typically results in overshoots and "
                             "undershoots, but numerical diffusion is minimized.  Full: Overshoots "
                             "and undershoots are avoided, but numerical diffusion is large");
  MooseEnum dislocation_character("edge screw", "edge");
  params.addRequiredParam<MooseEnum>(
      "dislocation_character", dislocation_character, "Character of dislocation");
  MooseEnum dislocation_sign("positive negative", "positive");
  params.addRequiredParam<MooseEnum>("dislocation_sign", dislocation_sign, "Sign of dislocation");
  params.addRequiredParam<int>("slip_system_index",
                               "Slip system index to get slip direction"
                               "FCC: 1 to 12.");
  return params;
}

ConservativeAdvectionDislocation::ConservativeAdvectionDislocation(
    const InputParameters & parameters)
  : Kernel(parameters),
    _dislocationcharacter(
        getParam<MooseEnum>("dislocation_character").getEnum<DislocationCharacter>()),
    _dislocationsign(getParam<MooseEnum>("dislocation_sign").getEnum<DislocationSign>()),
    _slip_direction_edge(getMaterialProperty<std::vector<RealVectorValue>>("slip_direction_edge")),
    _slip_plane_normalboth(
        getMaterialProperty<std::vector<RealVectorValue>>("slip_plane_normalboth")),
    _slip_system_index(getParam<int>("slip_system_index")),
    _dislo_velocity_CP_edge(getMaterialProperty<std::vector<Real>>("dislo_velocity_edge")),
    _dislo_velocity_CP_screw(getMaterialProperty<std::vector<Real>>("dislo_velocity_screw")),
    _upwinding(getParam<MooseEnum>("upwinding_type").getEnum<UpwindingType>()),
    _u_nodal(_var.dofValues()),
    _upwind_node(0),
    _dtotal_mass_out(0)
{
}

This is the validParameters and constructor of the "InterfaceDiffusion" interface kernel:

InputParameters
InterfaceDiffusionDislocation::validParams()
{
  InputParameters params = InterfaceKernel::validParams();
  params.addParam<Real>("density_critical", 1.0,"Critical density beyond which there will be dislocation transfer across Grain Boundary"); 
  params.addParam<Real>("tau_critical", 0.0,"Critical resolved shear stres beyond which there will be dislocation transfer across Grain Boundary");   
  MooseEnum dislocation_character("edge screw", "edge");
  params.addParam<MooseEnum>("dislocation_character", dislocation_character, "Character of dislocation");
  MooseEnum dislocation_sign("positive negative", "positive");
  params.addParam<MooseEnum>("dislocation_sign", dislocation_sign, "Sign of dislocation");
  params.addRequiredParam<int>("slip_system_index",
                               "Slip system index to get slip direction"
                               "FCC: 1 to 12.");
  params.addRequiredParam<int>("slip_system_index_neighbor",
                               "Slip system index to get slip direction of the neighbor"
                               "FCC: 1 to 12.");
  params.addClassDescription(
      "The kernel is utilized to establish flux equivalence on an interface for dislocation" 
      "assuming a divergence free velocity.");
  return params;
}

InterfaceDiffusionDislocation::InterfaceDiffusionDislocation(const InputParameters & parameters)
  : InterfaceKernel(parameters),
    _density_critical(getParam<Real>("density_critical")),
    _tau_critical(getParam<Real>("tau_critical")),
    _slip_system_index(getParam<int>("slip_system_index")),
    _slip_system_index_neighbor(getParam<int>("slip_system_index_neighbor")),
    _dislo_velocity_CP_edge(getMaterialProperty<std::vector<Real>>("dislo_velocity_edge")),
    _dislo_velocity_CP_edge_neighbor(getNeighborMaterialProperty<std::vector<Real>>("dislo_velocity_edge")),
    _dislocationcharacter(getParam<MooseEnum>("dislocation_character").getEnum<DislocationCharacter>()),
    _dislocationsign(getParam<MooseEnum>("dislocation_sign").getEnum<DislocationSign>()), 
    _slip_direction_edge(getMaterialProperty<std::vector<RealVectorValue>>("slip_direction_edge")),
    _slip_plane_normalboth(getMaterialProperty<std::vector<RealVectorValue>>("slip_plane_normalboth")),
    _slip_direction_edge_neighbor(getNeighborMaterialProperty<std::vector<RealVectorValue>>("slip_direction_edge")),
    _slip_plane_normalboth_neighbor(getNeighborMaterialProperty<std::vector<RealVectorValue>>("slip_plane_normalboth")),
    _tau(getMaterialProperty<std::vector<Real>>("applied_shear_stress")),
    _slip_resistance(getMaterialProperty<std::vector<Real>>("slip_resistance"))
{
}

In the interface kernel, the only variable that is correct is _slip_resistance. The other variables are always set to zero. Can you please help me with this.

Thanks a lor for your help!

Best,
Brayan

GiudGiud · 2024-10-23T14:00:04Z

GiudGiud
Oct 23, 2024
Collaborator

Hello

You should share your input and the header as well for the interface kernel.

The interface kernel is set to act on an internal boundary between two subdomains right?

2 replies

bmurgas-LANL Oct 23, 2024
Author

Hi @GiudGiud,
Thank you for the reply.
This is the repository I am working on: https://github.com/bmurgas-LANL/discopanther
I declare the property in the file DiscoFluxCPStressUpdate.C
For simplicity, I created a copy of the interfacediffusion kernel named interfacediffusiontest in https://github.com/bmurgas-LANL/discopanther/blob/main/src/interfacekernels/InterfaceDiffusionTest.C
In this modified interfaceDiffusion kernel I just tried to read the _slip_direction_edge material property created in DiscoFluxCPStressUpdate.C .
This is the source file:

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "InterfaceDiffusionTest.h"

registerMooseObject("MooseApp", InterfaceDiffusionTest);

InputParameters
InterfaceDiffusionTest::validParams()
{
  InputParameters params = InterfaceKernel::validParams();
  params.addParam<MaterialPropertyName>("D", "D", "The diffusion coefficient.");
  params.addParam<MaterialPropertyName>(
      "D_neighbor", "D_neighbor", "The neighboring diffusion coefficient.");
  params.addClassDescription(
      "The kernel is utilized to establish flux equivalence on an interface for variables.");
  return params;
}

InterfaceDiffusionTest::InterfaceDiffusionTest(const InputParameters & parameters)
  : InterfaceKernel(parameters),
    _D(getMaterialProperty<Real>("D")),
    _D_neighbor(getNeighborMaterialProperty<Real>("D_neighbor")),
    _slip_direction_edge(getMaterialProperty<std::vector<RealVectorValue>>("slip_direction_edge")),
    _slip_plane_normalboth(
        getMaterialProperty<std::vector<RealVectorValue>>("slip_plane_normalboth"))
{
}

Real
InterfaceDiffusionTest::computeQpResidual(Moose::DGResidualType type)
{
  Real r = 0;

  switch (type)
  {
    case Moose::Element:
      r = _test[_i][_qp] * -_D_neighbor[_qp] * _grad_neighbor_value[_qp] * _normals[_qp];
      break;

    case Moose::Neighbor:
      r = _test_neighbor[_i][_qp] * _D[_qp] * _grad_u[_qp] * _normals[_qp];
      break;
  }

  return r;
}

Real
InterfaceDiffusionTest::computeQpJacobian(Moose::DGJacobianType type)
{
  Real jac = 0;

  switch (type)
  {
    case Moose::ElementElement:
    case Moose::NeighborNeighbor:
      break;

    case Moose::NeighborElement:
      jac = _test_neighbor[_i][_qp] * _D[_qp] * _grad_phi[_j][_qp] * _normals[_qp];
      break;

    case Moose::ElementNeighbor:
      jac = _test[_i][_qp] * -_D_neighbor[_qp] * _grad_phi_neighbor[_j][_qp] * _normals[_qp];
      break;
  }

  return jac;
}

And this is the header file:

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "InterfaceKernel.h"

/**
 * DG kernel for interfacing diffusion between two variables on adjacent blocks
 */
class InterfaceDiffusionTest : public InterfaceKernel
{
public:
  static InputParameters validParams();

  InterfaceDiffusionTest(const InputParameters & parameters);

protected:
  virtual Real computeQpResidual(Moose::DGResidualType type) override;
  virtual Real computeQpJacobian(Moose::DGJacobianType type) override;

  const MaterialProperty<Real> & _D;
  const MaterialProperty<Real> & _D_neighbor;

  // dislocation velocities
  // slip direction and normal element
  const MaterialProperty<std::vector<RealVectorValue>> & _slip_direction_edge;
  const MaterialProperty<std::vector<RealVectorValue>> & _slip_plane_normalboth;
};

In this case, the two material properties are not correct. The value is always zero.

Thanks again for your help.

Best,
Brayan

bmurgas-LANL Oct 23, 2024
Author

I forgot the input file:

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh] 
  [./nepermesh]
    type = FileMeshGenerator
    file = '/Users/bmurgas/workspace/meshfiles/bar/px3-bar/px3_bar_hex8_cl16e-3.e'
  [../]
  [create_sideset]
    type = SideSetsFromNodeSetsGenerator
    input = nepermesh
  []

  # [./x0_modifier]
  #   type = BoundingBoxNodeSetGenerator
  #   input = nepermesh
  #   new_boundary = x0
  #   top_right = '0.1 1.1 0.35'
  #   bottom_left = '-0.1 -0.1 -0.1'
  # [../]
  # [./y0_modifier]
  #   type = BoundingBoxNodeSetGenerator
  #   input = x0_modifier
  #   new_boundary = y0
  #   top_right = '1.1 0.1 0.35'
  #   bottom_left = '-0.1 -0.1 -0.1'
  # [../]
  # [./z0_modifier]
  #   type = BoundingBoxNodeSetGenerator
  #   input = y0_modifier
  #   new_boundary = z0
  #   top_right = '1.1 1.1 0.1'
  #   bottom_left = '-0.1 -0.1 -0.1'
  # [../]
  # [./x1_modifier]
  #   type = BoundingBoxNodeSetGenerator
  #   input = z0_modifier
  #   new_boundary = x1
  #   top_right = '1.1 1.1 0.35'
  #   bottom_left = '0.9 -0.1 -0.1'
  # [../]
  # [./y1_modifier]
  #   type = BoundingBoxNodeSetGenerator
  #   input = x1_modifier
  #   new_boundary = y1
  #   top_right = '1.1 1.1 0.35'
  #   bottom_left = '-0.1 0.9 -0.1'
  # [../]
  displacements = 'disp_x disp_y disp_z'
[]

# [Debug]
#   #show_execution_order = 'ALWAYS'
#   show_material_props = true
#   show_var_residual_norms = true
# []

[Outputs]
    file_base = method_Discoflux
    csv = true
    exodus = true
    console = true
[]

[UserObjects]
  [./prop_read]
    type = PropertyReadFile
    prop_file_name = '/Users/bmurgas/workspace/auxfiles/px3-bar-sp10-sp01.inp' #sxori-0-0-0.inp' #px3-bar-sp10-sp01.inp'
    nprop = 3
    read_type = block
    nblock= 3
  [../]
[]

[Variables]
  [./disp_x]
    order = FIRST
    family = LAGRANGE
  [../]

  [./disp_y]
    order = FIRST
    family = LAGRANGE
  [../]

  [./disp_z]
    order = FIRST
    family = LAGRANGE
  [../]
  # DISLOCATION DENSITIES
  ## POSITIVE EDGE
  [./DD_EdgePositive_1]
    order = FIRST
    family = LAGRANGE
    # [./InitialCondition]
    #   type = FunctionIC
    #   function = init_rho_edge_pos_1
    #   # block = 1
    # [../]
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgePositive_2]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgePositive_3]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgePositive_4]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgePositive_5]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgePositive_6]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgePositive_7]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgePositive_8]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgePositive_9]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgePositive_10]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgePositive_11]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgePositive_12]
    order = FIRST
    family = LAGRANGE
    # [./InitialCondition]
    #   type = ConstantIC
    #   value = 0.5
    # [../]
    [./InitialCondition]
      type = FunctionIC
      function = init_rho_edge_pos_1
      # block = 1
    [../]
  [../]
  ## NEGATIVE EDGE
  [./DD_EdgeNegative_1]
    order = FIRST
    family = LAGRANGE
    # [./InitialCondition]
    #   type = FunctionIC
    #   function = init_rho_edge_neg_1
    #   # block = 1
    # [../]
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
      # block = 1
    [../]
  [../]
  [./DD_EdgeNegative_2]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgeNegative_3]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgeNegative_4]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgeNegative_5]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgeNegative_6]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgeNegative_7]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgeNegative_8]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgeNegative_9]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgeNegative_10]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_EdgeNegative_11]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC 
      value = 0.5
    [../]
  [../]
  [./DD_EdgeNegative_12]
    order = FIRST
    family = LAGRANGE
    # [./InitialCondition]
    #   type = ConstantIC
    #   value = 0.5
    # [../]
    [./InitialCondition]
      type = FunctionIC
      function = init_rho_edge_neg_1
      # block = 1
    [../]
  [../]
  # POSITIVE SCREW
  [./DD_ScrewPositive_1]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewPositive_2]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewPositive_3]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewPositive_4]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewPositive_5]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewPositive_6]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewPositive_7]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewPositive_8]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewPositive_9]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewPositive_10]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewPositive_11]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewPositive_12]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  ## NEGATIVE Screw
  [./DD_ScrewNegative_1]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewNegative_2]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewNegative_3]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewNegative_4]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewNegative_5]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewNegative_6]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewNegative_7]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewNegative_8]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewNegative_9]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewNegative_10]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewNegative_11]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
  [./DD_ScrewNegative_12]
    order = FIRST
    family = LAGRANGE
    [./InitialCondition]
      type = ConstantIC
      value = 0.5
    [../]
  [../]
[]

[Functions]
  [./init_rho_edge_pos_1]
    type = ParsedFunction
    value = 'if(x>=0.24,1.0,0.0)*if(x<=0.40,1.0,0.0)*if(z>=0.08,1.0,0.0)*if(z<=0.24,1.0,0.0)*0.4+0.1'
  [../]
  [./init_rho_edge_neg_1]
    type = ParsedFunction
	  value = 'if(x>=0.24,1.0,0.0)*if(x<=0.40,1.0,0.0)*if(z>=0.08,1.0,0.0)*if(z<=0.24,1.0,0.0)*0.4+0.1'
  [../]
  [./disp_load]
    type = ParsedFunction
    value = '0.1*t'
  [../]
[]

[Kernels]
  [./dot_DD_EdgePositive_1]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgePositive_1
  [../]
  # [./Advection_DD_EdgePositive_1]
  #   type = ConservativeAdvectionDislocation
  #   variable = DD_EdgePositive_1
  #   upwinding_type = full
  #   dislocation_character = edge
  #   dislocation_sign = positive
  #   slip_system_index = 1
  # [../]
  [./dot_DD_EdgePositive_2]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgePositive_2
  [../]
  #   [./Advection_DD_EdgePositive_2]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgePositive_2
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = positive
  #     slip_system_index = 2
  #   [../]
  [./dot_DD_EdgePositive_3]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgePositive_3
  [../]
  #   [./Advection_DD_EdgePositive_3]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgePositive_3
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = positive
  #     slip_system_index = 3
  #   [../]
  [./dot_DD_EdgePositive_4]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgePositive_4
  [../]
  #   [./Advection_DD_EdgePositive_4]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgePositive_4
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = positive
  #     slip_system_index = 4
  #   [../]
  [./dot_DD_EdgePositive_5]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgePositive_5
  [../]
  #   [./Advection_DD_EdgePositive_5]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgePositive_5
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = positive
  #     slip_system_index = 5
  #   [../]
  [./dot_DD_EdgePositive_6]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgePositive_6
  [../]
  #   [./Advection_DD_EdgePositive_6]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgePositive_6
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = positive
  #     slip_system_index = 6
  #   [../]
  [./dot_DD_EdgePositive_7]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgePositive_7
  [../]
  #   [./Advection_DD_EdgePositive_7]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgePositive_7
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = positive
  #     slip_system_index = 7
  #   [../]
  [./dot_DD_EdgePositive_8]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgePositive_8
  [../]
  #   [./Advection_DD_EdgePositive_8]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgePositive_8
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = positive
  #     slip_system_index = 8
  #   [../]
  [./dot_DD_EdgePositive_9]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgePositive_9
  [../]
  #   [./Advection_DD_EdgePositive_9]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgePositive_9
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = positive
  #     slip_system_index = 9
  #   [../]
  [./dot_DD_EdgePositive_10]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgePositive_10
  [../]
  #   [./Advection_DD_EdgePositive_10]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgePositive_10
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = positive
  #     slip_system_index = 10
  #   [../]
  [./dot_DD_EdgePositive_11]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgePositive_11
  [../]
  #   [./Advection_DD_EdgePositive_11]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgePositive_11
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = positive
  #     slip_system_index = 11
  #   [../]
  [./dot_DD_EdgePositive_12]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgePositive_12
  [../]
    [./Advection_DD_EdgePositive_12]
      type = ConservativeAdvectionDislocation
      variable = DD_EdgePositive_12
      upwinding_type = full
      dislocation_character = edge
      dislocation_sign = positive
      slip_system_index = 12
    [../]
  
  [./dot_DD_EdgeNegative_1]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgeNegative_1
  [../]
    # [./Advection_DD_EdgeNegative_1]
    #   type = ConservativeAdvectionDislocation
    #   variable = DD_EdgeNegative_1
    #   upwinding_type = full
    #   dislocation_character = edge
    #   dislocation_sign = negative
    #   slip_system_index = 1
    # [../]
  [./dot_DD_EdgeNegative_2]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgeNegative_2
  [../]
  #   [./Advection_DD_EdgeNegative_2]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgeNegative_2
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = negative
  #     slip_system_index = 2
  #   [../]
  [./dot_DD_EdgeNegative_3]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgeNegative_3
  [../]
  #   [./Advection_DD_EdgeNegative_3]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgeNegative_3
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = negative
  #     slip_system_index = 3
  #   [../]
  [./dot_DD_EdgeNegative_4]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgeNegative_4
  [../]
  #   [./Advection_DD_EdgeNegative_4]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgeNegative_4
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = negative
  #     slip_system_index = 4
  #   [../]
  [./dot_DD_EdgeNegative_5]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgeNegative_5
  [../]
  #   [./Advection_DD_EdgeNegative_5]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgeNegative_5
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = negative
  #     slip_system_index = 5
  #   [../]
  [./dot_DD_EdgeNegative_6]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgeNegative_6
  [../]
  #   [./Advection_DD_EdgeNegative_6]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgeNegative_6
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = negative
  #     slip_system_index = 6
  #   [../]
  [./dot_DD_EdgeNegative_7]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgeNegative_7
  [../]
  #   [./Advection_DD_EdgeNegative_7]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgeNegative_7
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = negative
  #     slip_system_index = 7
  #   [../]
  [./dot_DD_EdgeNegative_8]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgeNegative_8
  [../]
  #   [./Advection_DD_EdgeNegative_8]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgeNegative_8
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = negative
  #     slip_system_index = 8
  #   [../]
  [./dot_DD_EdgeNegative_9]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgeNegative_9
  [../]
  #   [./Advection_DD_EdgeNegative_9]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgeNegative_9
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = negative
  #     slip_system_index = 9
  #   [../]
  [./dot_DD_EdgeNegative_10]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgeNegative_10
  [../]
  #   [./Advection_DD_EdgeNegative_10]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgeNegative_10
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = negative
  #     slip_system_index = 10
  #   [../]
  [./dot_DD_EdgeNegative_11]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgeNegative_11
  [../]
  #   [./Advection_DD_EdgeNegative_11]
  #     type = ConservativeAdvectionDislocation
  #     variable = DD_EdgeNegative_11
  #     upwinding_type = full
  #     dislocation_character = edge
  #     dislocation_sign = negative
  #     slip_system_index = 11
  #   [../]
  [./dot_DD_EdgeNegative_12]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_EdgeNegative_12
  [../]
    [./Advection_DD_EdgeNegative_12]
      type = ConservativeAdvectionDislocation
      variable = DD_EdgeNegative_12
      upwinding_type = full
      dislocation_character = edge
      dislocation_sign = negative
      slip_system_index = 12
    [../]
  
  
  [./dot_DD_ScrewPositive_1]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewPositive_1
  [../]
  [./dot_DD_ScrewPositive_2]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewPositive_2
  [../]
  [./dot_DD_ScrewPositive_3]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewPositive_3
  [../]
  [./dot_DD_ScrewPositive_4]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewPositive_4
  [../]
  [./dot_DD_ScrewPositive_5]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewPositive_5
  [../]
  [./dot_DD_ScrewPositive_6]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewPositive_6
  [../]
  [./dot_DD_ScrewPositive_7]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewPositive_7
  [../]
  [./dot_DD_ScrewPositive_8]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewPositive_8
  [../]
  [./dot_DD_ScrewPositive_9]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewPositive_9
  [../]
  [./dot_DD_ScrewPositive_10]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewPositive_10
  [../]
  [./dot_DD_ScrewPositive_11]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewPositive_11
  [../]
  [./dot_DD_ScrewPositive_12]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewPositive_12
  [../]
  
  [./dot_DD_ScrewNegative_1]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewNegative_1
  [../]
  [./dot_DD_ScrewNegative_2]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewNegative_2
  [../]
  [./dot_DD_ScrewNegative_3]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewNegative_3
  [../]
  [./dot_DD_ScrewNegative_4]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewNegative_4
  [../]
  [./dot_DD_ScrewNegative_5]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewNegative_5
  [../]
  [./dot_DD_ScrewNegative_6]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewNegative_6
  [../]
  [./dot_DD_ScrewNegative_7]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewNegative_7
  [../]
  [./dot_DD_ScrewNegative_8]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewNegative_8
  [../]
  [./dot_DD_ScrewNegative_9]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewNegative_9
  [../]
  [./dot_DD_ScrewNegative_10]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewNegative_10
  [../]
  [./dot_DD_ScrewNegative_11]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewNegative_11
  [../]
  [./dot_DD_ScrewNegative_12]
    # type = TimeDerivative  
    type = MassLumpedTimeDerivative
    variable = DD_ScrewNegative_12
  [../]
[]

[InterfaceKernels]
  [./Intf_Avd_DD_EdgePositive_12_1]
    type = InterfaceDiffusionTest
    variable = DD_EdgePositive_12
    neighbor_var = DD_EdgePositive_1
    boundary = surface_GB
    D = 1.0
    D_neighbor = 1.0
  [../]
  # [./Intf_Avd_DD_EdgePositive_12_1]
  #   type = InterfaceDiffusionDislocation
  #   variable = DD_EdgePositive_12
  #   slip_system_index = 12
  #   neighbor_var = DD_EdgePositive_1
  #   slip_system_index_neighbor = 1
  #   dislocation_character = edge
  #   dislocation_sign = positive
  #   boundary = surface_GB
  #   density_critical = 0.2
  #   tau_critical = 2
  # [../]
  [./Intf_Avd_DD_EdgePositive_12_2]
    type = InterfaceDiffusionDislocation
    variable = DD_EdgePositive_12
    slip_system_index = 12
    neighbor_var = DD_EdgePositive_2
    slip_system_index_neighbor = 2
    dislocation_character = edge
    dislocation_sign = positive
    boundary = surface_GB
    density_critical = 0.2
    tau_critical = 2
  [../]
  [./Intf_Avd_DD_EdgePositive_12_3]
    type = InterfaceDiffusionDislocation
    variable = DD_EdgePositive_12
    slip_system_index = 12
    neighbor_var = DD_EdgePositive_3
    slip_system_index_neighbor = 3
    dislocation_character = edge
    dislocation_sign = positive
    boundary = surface_GB
    density_critical = 0.2
    tau_critical = 2
  [../]
  [./Intf_Avd_DD_EdgePositive_12_4]
    type = InterfaceDiffusionDislocation
    variable = DD_EdgePositive_12
    slip_system_index = 12
    neighbor_var = DD_EdgePositive_4
    slip_system_index_neighbor = 4
    dislocation_character = edge
    dislocation_sign = positive
    boundary = surface_GB
    density_critical = 0.2
    tau_critical = 2
  [../]
  [./Intf_Avd_DD_EdgePositive_12_5]
    type = InterfaceDiffusionDislocation
    variable = DD_EdgePositive_12
    slip_system_index = 12
    neighbor_var = DD_EdgePositive_5
    slip_system_index_neighbor = 5
    dislocation_character = edge
    dislocation_sign = positive
    boundary = surface_GB
    density_critical = 0.2
    tau_critical = 2
  [../]
  [./Intf_Avd_DD_EdgePositive_12_6]
    type = InterfaceDiffusionDislocation
    variable = DD_EdgePositive_12
    slip_system_index = 12
    neighbor_var = DD_EdgePositive_6
    slip_system_index_neighbor = 6
    dislocation_character = edge
    dislocation_sign = positive
    boundary = surface_GB
    density_critical = 0.2
    tau_critical = 2
  [../]
  [./Intf_Avd_DD_EdgePositive_12_7]
    type = InterfaceDiffusionDislocation
    variable = DD_EdgePositive_12
    slip_system_index = 12
    neighbor_var = DD_EdgePositive_7
    slip_system_index_neighbor = 7
    dislocation_character = edge
    dislocation_sign = positive
    boundary = surface_GB
    density_critical = 0.2
    tau_critical = 2
  [../]
  [./Intf_Avd_DD_EdgePositive_12_8]
    type = InterfaceDiffusionDislocation
    variable = DD_EdgePositive_12
    slip_system_index = 12
    neighbor_var = DD_EdgePositive_8
    slip_system_index_neighbor = 8
    dislocation_character = edge
    dislocation_sign = positive
    boundary = surface_GB
    density_critical = 0.2
    tau_critical = 2
  [../]
  [./Intf_Avd_DD_EdgePositive_12_9]
    type = InterfaceDiffusionDislocation
    variable = DD_EdgePositive_12
    slip_system_index = 12
    neighbor_var = DD_EdgePositive_9
    slip_system_index_neighbor = 9
    dislocation_character = edge
    dislocation_sign = positive
    boundary = surface_GB
    density_critical = 0.2
    tau_critical = 2
  [../]
  [./Intf_Avd_DD_EdgePositive_12_10]
    type = InterfaceDiffusionDislocation
    variable = DD_EdgePositive_12
    slip_system_index = 12
    neighbor_var = DD_EdgePositive_10
    slip_system_index_neighbor = 10
    dislocation_character = edge
    dislocation_sign = positive
    boundary = surface_GB
    density_critical = 0.2
    tau_critical = 2
  [../]
  [./Intf_Avd_DD_EdgePositive_12_11]
    type = InterfaceDiffusionDislocation
    variable = DD_EdgePositive_12
    slip_system_index = 12
    neighbor_var = DD_EdgePositive_11
    slip_system_index_neighbor = 11
    dislocation_character = edge
    dislocation_sign = positive
    boundary = surface_GB
    density_critical = 0.2
    tau_critical = 2
  [../]
  [./Intf_Avd_DD_EdgePositive_12_12]
    type = InterfaceDiffusionDislocation
    variable = DD_EdgePositive_12
    slip_system_index = 12
    neighbor_var = DD_EdgePositive_12
    slip_system_index_neighbor = 12
    dislocation_character = edge
    dislocation_sign = positive
    boundary = surface_GB
    density_critical = 0.2
    tau_critical = 2
  [../]
  # [./Intf_Avd_DD_EdgePositive_1_12]
  #   type = InterfaceDiffusionDislocation
  #   variable = DD_EdgePositive_1
  #   slip_system_index = 1
  #   neighbor_var = DD_EdgePositive_12
  #   slip_system_index_neighbor = 12
  #   dislocation_character = edge
  #   dislocation_sign = positive
  #   boundary = surface_GB
  #   density_critical = 0.2
  #   tau_critical = 2
  # [../]
  #   [./Intf_Avd_DD_EdgePositive_2_12]
  #     type = InterfaceDiffusionDislocation
  #     variable = DD_EdgePositive_2
  #     slip_system_index = 2
  #     neighbor_var = DD_EdgePositive_12
  #     slip_system_index_neighbor = 12
  #     dislocation_character = edge
  #     dislocation_sign = positive
  #     boundary = surface_GB
  #     density_critical = 0.2
  #     tau_critical = 2
  #   [../]
  #     [./Intf_Avd_DD_EdgePositive_3_12]
  #       type = InterfaceDiffusionDislocation
  #       variable = DD_EdgePositive_3
  #       slip_system_index = 3
  #       neighbor_var = DD_EdgePositive_12
  #       slip_system_index_neighbor = 12
  #       dislocation_character = edge
  #       dislocation_sign = positive
  #       boundary = surface_GB
  #       density_critical = 0.2
  #       tau_critical = 2
  #     [../]
  #       [./Intf_Avd_DD_EdgePositive_4_12]
  #         type = InterfaceDiffusionDislocation
  #         variable = DD_EdgePositive_4
  #         slip_system_index = 4
  #         neighbor_var = DD_EdgePositive_12
  #         slip_system_index_neighbor = 12
  #         dislocation_character = edge
  #         dislocation_sign = positive
  #         boundary = surface_GB
  #         density_critical = 0.2
  #         tau_critical = 2
  #       [../]
  #         [./Intf_Avd_DD_EdgePositive_5_12]
  #           type = InterfaceDiffusionDislocation
  #           variable = DD_EdgePositive_5
  #           slip_system_index = 5
  #           neighbor_var = DD_EdgePositive_12
  #           slip_system_index_neighbor = 12
  #           dislocation_character = edge
  #           dislocation_sign = positive
  #           boundary = surface_GB
  #           density_critical = 0.2
  #           tau_critical = 2
  #         [../]
  #           [./Intf_Avd_DD_EdgePositive_6_12]
  #             type = InterfaceDiffusionDislocation
  #             variable = DD_EdgePositive_6
  #             slip_system_index = 6
  #             neighbor_var = DD_EdgePositive_12
  #             slip_system_index_neighbor = 12
  #             dislocation_character = edge
  #             dislocation_sign = positive
  #             boundary = surface_GB
  #             density_critical = 0.2
  #             tau_critical = 2
  #           [../]
  #             [./Intf_Avd_DD_EdgePositive_7_12]
  #               type = InterfaceDiffusionDislocation
  #               variable = DD_EdgePositive_7
  #               slip_system_index = 7
  #               neighbor_var = DD_EdgePositive_12
  #               slip_system_index_neighbor = 12
  #               dislocation_character = edge
  #               dislocation_sign = positive
  #               boundary = surface_GB
  #               density_critical = 0.2
  #               tau_critical = 2
  #             [../]
  #               [./Intf_Avd_DD_EdgePositive_8_12]
  #                 type = InterfaceDiffusionDislocation
  #                 variable = DD_EdgePositive_8
  #                 slip_system_index = 8
  #                 neighbor_var = DD_EdgePositive_12
  #                 slip_system_index_neighbor = 12
  #                 dislocation_character = edge
  #                 dislocation_sign = positive
  #                 boundary = surface_GB
  #                 density_critical = 0.2
  #                 tau_critical = 2
  #               [../]
  #                 [./Intf_Avd_DD_EdgePositive_9_12]
  #                   type = InterfaceDiffusionDislocation
  #                   variable = DD_EdgePositive_9
  #                   slip_system_index = 9
  #                   neighbor_var = DD_EdgePositive_12
  #                   slip_system_index_neighbor = 12
  #                   dislocation_character = edge
  #                   dislocation_sign = positive
  #                   boundary = surface_GB
  #                   density_critical = 0.2
  #                   tau_critical = 2
  #                 [../]
  #                   [./Intf_Avd_DD_EdgePositive_10_12]
  #                     type = InterfaceDiffusionDislocation
  #                     variable = DD_EdgePositive_10
  #                     slip_system_index = 10
  #                     neighbor_var = DD_EdgePositive_12
  #                     slip_system_index_neighbor = 12
  #                     dislocation_character = edge
  #                     dislocation_sign = positive
  #                     boundary = surface_GB
  #                     density_critical = 0.2
  #                     tau_critical = 2
  #                   [../]
  #                     [./Intf_Avd_DD_EdgePositive_11_12]
  #                       type = InterfaceDiffusionDislocation
  #                       variable = DD_EdgePositive_11
  #                       slip_system_index = 11
  #                       neighbor_var = DD_EdgePositive_12
  #                       slip_system_index_neighbor = 12
  #                       dislocation_character = edge
  #                       dislocation_sign = positive
  #                       boundary = surface_GB
  #                       density_critical = 0.2
  #                       tau_critical = 2
  #                     [../]
[]

[AuxVariables]
  [euler_angle_1]
    order = CONSTANT
    family = MONOMIAL
  []
  [euler_angle_2]
    order = CONSTANT
    family = MONOMIAL
  []
  [euler_angle_3]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[AuxKernels]
  [euler_angle_1]
    type = MaterialRealVectorValueAux
    variable = euler_angle_1
    property = updated_Euler_angle
    component = 0
    execute_on = 'initial final'
  []
  [euler_angle_2]
    type = MaterialRealVectorValueAux
    variable = euler_angle_2
    property = updated_Euler_angle
    component = 1
    execute_on = 'initial final'
  []
  [euler_angle_3]
    type = MaterialRealVectorValueAux
    variable = euler_angle_3
    property = updated_Euler_angle
    component = 2
    execute_on = 'initial final'
  []
[]

[BCs]
  ##############################################
  # X : Left and Right
  # Y : Bottom and Top
  # Z : Back and Front
  ##############################################
  [./BC_X]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = 'x0'
    # boundary = left
    function = 0.0
  [../]
    [./BC_Y]
    type = FunctionDirichletBC
    variable = disp_y
    boundary = 'y0'
    # boundary = bottom
    function = 0.0
  [../]
    [./BC_Z]
    type = FunctionDirichletBC
    variable = disp_z
    boundary = 'z0'
    # boundary = back
    function = 0.0
  [../]

  [./BC_Loading]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = 'x0 x1 '
    # boundary = top
    function = '0.001*x*t'
  [../]

  # [./BC_X0]
  #   type = FunctionDirichletBC
  #   variable = disp_x
  #   boundary = 'x0 y0 z0'
  #   function = 0.0
  # [../]
  # [./BC_Y0]
  #   type = FunctionDirichletBC
  #   variable = disp_y
  #   boundary = 'x0 y0 z0'
  #   function = 0.0
  # [../]
  # [./BC_Z0]
  #   type = FunctionDirichletBC
  #   variable = disp_z
  #   boundary = 'x0 y0 z0'
  #   function = 0.0
  # [../]

  # [./Periodic]
  #   [./rho_edge_positive_1]
  #     variable = 'DD_EdgePositive_1 DD_EdgePositive_2 DD_EdgePositive_3 DD_EdgePositive_4 DD_EdgePositive_5 DD_EdgePositive_6 DD_EdgePositive_7 DD_EdgePositive_8 DD_EdgePositive_9 DD_EdgePositive_10 DD_EdgePositive_11 DD_EdgePositive_12'
  #     auto_direction = 'x y z'
  #   [../]
  #   [./rho_edge_negative]
  #     variable = 'DD_EdgeNegative_1 DD_EdgeNegative_2 DD_EdgeNegative_3 DD_EdgeNegative_4 DD_EdgeNegative_5 DD_EdgeNegative_6 DD_EdgeNegative_7 DD_EdgeNegative_8 DD_EdgeNegative_9 DD_EdgeNegative_10 DD_EdgeNegative_11 DD_EdgeNegative_12'
  #     auto_direction = 'x y z'
  #   [../]
  # [../] 
  
[]

[Physics/SolidMechanics/QuasiStatic/all]
  strain = FINITE
  add_variables = true
  new_system = true
  formulation = total
  incremental = true
[]

[Materials]
  [./elasticity_tensor]
    type = ComputeElasticityTensorCP
    C_ijkl = '1.684e5 1.214e5 1.214e5 1.684e5 1.214e5 1.684e5 0.754e5 0.754e5 0.754e5'
    fill_method = symmetric9
    read_prop_user_object = prop_read
  [../]
  [stress_brass]
    type = ComputeMultipleCrystalPlasticityStress
    crystal_plasticity_models = 'CP_DiscoFlux'
    tan_mod_type = exact
  []
  [./CP_DiscoFlux]
    type = DiscoFluxCPStressUpdate
    number_slip_systems = 12
    slip_sys_file_name = /Users/bmurgas/workspace/auxfiles/input_slip_sys.inp
    lattice_friction = 15
    Coeff_hardening = 0.85
    slip_increment_tolerance = 2.0e-2 #2.0e-2 10.0
    DD_EdgePositive_1 = DD_EdgePositive_1
    DD_EdgePositive_2 = DD_EdgePositive_2
    DD_EdgePositive_3 = DD_EdgePositive_3
    DD_EdgePositive_4 = DD_EdgePositive_4
    DD_EdgePositive_5 = DD_EdgePositive_5
    DD_EdgePositive_6 = DD_EdgePositive_6
    DD_EdgePositive_7 = DD_EdgePositive_7
    DD_EdgePositive_8 = DD_EdgePositive_8
    DD_EdgePositive_9 = DD_EdgePositive_9
    DD_EdgePositive_10 = DD_EdgePositive_10
    DD_EdgePositive_11 = DD_EdgePositive_11
    DD_EdgePositive_12 = DD_EdgePositive_12
    DD_EdgeNegative_1 = DD_EdgeNegative_1
    DD_EdgeNegative_2 = DD_EdgeNegative_2
    DD_EdgeNegative_3 = DD_EdgeNegative_3
    DD_EdgeNegative_4 = DD_EdgeNegative_4
    DD_EdgeNegative_5 = DD_EdgeNegative_5
    DD_EdgeNegative_6 = DD_EdgeNegative_6
    DD_EdgeNegative_7 = DD_EdgeNegative_7
    DD_EdgeNegative_8 = DD_EdgeNegative_8
    DD_EdgeNegative_9 = DD_EdgeNegative_9
    DD_EdgeNegative_10 = DD_EdgeNegative_10
    DD_EdgeNegative_11 = DD_EdgeNegative_11
    DD_EdgeNegative_12 = DD_EdgeNegative_12
    DD_ScrewPositive_1 = DD_ScrewPositive_1
    DD_ScrewPositive_2 = DD_ScrewPositive_2
    DD_ScrewPositive_3 = DD_ScrewPositive_3
    DD_ScrewPositive_4 = DD_ScrewPositive_4
    DD_ScrewPositive_5 = DD_ScrewPositive_5
    DD_ScrewPositive_6 = DD_ScrewPositive_6
    DD_ScrewPositive_7 = DD_ScrewPositive_7
    DD_ScrewPositive_8 = DD_ScrewPositive_8
    DD_ScrewPositive_9 = DD_ScrewPositive_9
    DD_ScrewPositive_10 = DD_ScrewPositive_10
    DD_ScrewPositive_11 = DD_ScrewPositive_11
    DD_ScrewPositive_12 = DD_ScrewPositive_12
    DD_ScrewNegative_1 = DD_ScrewNegative_1
    DD_ScrewNegative_2= DD_ScrewNegative_2
    DD_ScrewNegative_3 = DD_ScrewNegative_3
    DD_ScrewNegative_4 = DD_ScrewNegative_4
    DD_ScrewNegative_5 = DD_ScrewNegative_5
    DD_ScrewNegative_6 = DD_ScrewNegative_6
    DD_ScrewNegative_7 = DD_ScrewNegative_7
    DD_ScrewNegative_8 = DD_ScrewNegative_8
    DD_ScrewNegative_9 = DD_ScrewNegative_9
    DD_ScrewNegative_10 = DD_ScrewNegative_10
    DD_ScrewNegative_11 = DD_ScrewNegative_11
    DD_ScrewNegative_12 = DD_ScrewNegative_12
    print_state_variable_convergence_error_messages = true
  [../]
  [updated_euler_angle]
    type = ComputeUpdatedEulerAngle
    radian_to_degree = true
  []
[]

[Postprocessors]
  [euler_angle_1] 
    type = ElementAverageValue
    variable = euler_angle_1
  []
  [euler_angle_2]
    type = ElementAverageValue
    variable = euler_angle_2
  []
  [euler_angle_3]
    type = ElementAverageValue
    variable = euler_angle_3
  []
[]

[GlobalParams]
    dislo_density_initial = 1.0e+05
    C_multi =  0.0  #8.5e-06
    C_trap =   0.0  #5.5e-03
    C_m_ann =  0.0  #0.5
    C_im_ann = 0.0 #0.5
    burgers_vector_mag = 2.52e-07
[]

[Preconditioning]
  [./smp]
    type = SMP
    full = true
  [../]
[]

[Executioner]
  type = Transient

  # solve_type = 'PJFNK'
  # line_search = 'none'
  # petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -snes_type'
  # petsc_options_value = 'lu       superlu_dist                  vinewtonrsls'
  # nl_abs_tol = 1e-5 
  # nl_max_its = 10
  # l_abs_tol = 1e-5
  # l_max_its = 20

  # solver used in cpam
  solve_type = 'PJFNK'
  petsc_options = '-snes_ksp_ew'
  petsc_options_iname = '-pc_type -pc_hypre_type -ksp_gmres_restart'
  petsc_options_value = 'hypre    boomeramg          31'
  line_search = 'none'
  l_max_its = 50
  nl_max_its = 50
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-6
  l_tol = 1e-8

  dtmax = 0.2
  dtmin = 0.000001
  end_time = 100 #100
  [./TimeStepper]
    type = ConstantDT
    dt = 0.2
    growth_factor = 1.01
  [../]
[]

[AuxVariables]
  [./pk2_1]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./pk2_2]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./pk2_3]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./pk2_4]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./pk2_5]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./pk2_6]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./E_1]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./E_2]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./E_3]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./E_4]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./E_5]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./E_6]
    order = CONSTANT
    family = MONOMIAL
  [../]
[]

[AuxKernels]
  [./pk2_1] 
   type = RankTwoAux
   variable = pk2_1
   rank_two_tensor = stress
   index_j = 0
   index_i = 0
  [../]
  [./pk2_2] 
    type = RankTwoAux
    variable = pk2_2
    rank_two_tensor = stress
    index_j = 1
    index_i = 1
  [../]
  [./pk2_3] 
   type = RankTwoAux
   variable = pk2_3
   rank_two_tensor = stress
   index_j = 2
   index_i = 2
  [../]
  [./pk2_4] 
    type = RankTwoAux
    variable = pk2_4
    rank_two_tensor = stress
    index_j = 0
    index_i = 1
  [../]
  [./pk2_5] 
    type = RankTwoAux
    variable = pk2_5
    rank_two_tensor = stress
    index_j = 0
    index_i = 2
  [../]
  [./pk2_6] 
    type = RankTwoAux
    variable = pk2_6
    rank_two_tensor = stress
    index_j = 1
    index_i = 2
  [../]
  [./E_1]
   type = RankTwoAux
   variable = E_1
   rank_two_tensor = total_strain
   index_j = 0
   index_i = 0
  [../]
  [./E_2] 
    type = RankTwoAux
    variable = E_2
    rank_two_tensor = total_strain
    index_j = 1
    index_i = 1
  [../]
  [./E_3] 
   type = RankTwoAux
   variable = E_3
   rank_two_tensor = total_strain
   index_j = 2
   index_i = 2
  [../]
  [./E_4] 
    type = RankTwoAux
    variable = E_4
    rank_two_tensor = total_strain
    index_j = 0
    index_i = 1
  [../]
  [./E_5] 
    type = RankTwoAux
    variable = E_5
    rank_two_tensor = total_strain
    index_j = 0
    index_i = 2
  [../]
  [./E_6] 
    type = RankTwoAux
    variable = E_6
    rank_two_tensor = total_strain
    index_j = 1
    index_i = 2
  [../]
[]

[Postprocessors]
  [./pk2_1]
    type = ElementAverageValue
    variable = pk2_1
  [../]
  [./pk2_2]
    type = ElementAverageValue
    variable = pk2_2
  [../]
  [./pk2_3]
    type = ElementAverageValue
    variable = pk2_3
  [../]
  [./pk2_4]
    type = ElementAverageValue
    variable = pk2_4
  [../]
  [./pk2_5]
    type = ElementAverageValue
    variable = pk2_5
  [../]
  [./pk2_6]
    type = ElementAverageValue
    variable = pk2_6
  [../]
  [./E_1]
    type = ElementAverageValue
    variable = E_1
  [../]
  [./E_2]
    type = ElementAverageValue
    variable = E_2
  [../]
  [./E_3]
    type = ElementAverageValue
    variable = E_3
  [../]
  [./E_4]
    type = ElementAverageValue
    variable = E_4
  [../]
  [./E_5]
    type = ElementAverageValue
    variable = E_5
  [../]
  [./E_6]
    type = ElementAverageValue
    variable = E_6
  [../]
[]

# [AuxVariables]
#   [./DD_EP_01]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EP_02]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EP_03]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EP_04]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EP_05]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EP_06]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EP_07]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EP_08]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EP_09]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EP_10]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EP_11]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EP_12]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
# []


# [AuxKernels]
#   [./DD_EP_01]
#    type = MaterialRealAux
#    variable = DD_EP_01
#    property = DD_EdgePositive_1
#   [../]
#   [./DD_EP_02]
#    type = MaterialRealAux
#    variable = DD_EP_02
#    property = DD_EdgePositive_2  
#   [../]
#   [./DD_EP_03]
#    type = MaterialRealAux
#    variable = DD_EP_03
#    property = DD_EdgePositive_3
#   [../]
#   [./DD_EP_04]
#    type = MaterialRealAux
#    variable = DD_EP_04
#    property = DD_EdgePositive_4
#   [../]
#   [./DD_EP_05]
#    type = MaterialRealAux
#    variable = DD_EP_05
#    property = DD_EdgePositive_5
#   [../]
#   [./DD_EP_06]
#    type = MaterialRealAux
#    variable = DD_EP_06
#    property = DD_EdgePositive_6
#   [../]
#   [./DD_EP_07]
#    type = MaterialRealAux
#    variable = DD_EP_07
#    property = DD_EdgePositive_7
#   [../]
#   [./DD_EP_08]
#    type = MaterialRealAux
#    variable = DD_EP_08
#    property = DD_EdgePositive_8
#   [../]
#   [./DD_EP_09]
#    type = MaterialRealAux
#    variable = DD_EP_09
#    property = DD_EdgePositive_9
#   [../]
#   [./DD_EP_10]
#    type = MaterialRealAux
#    variable = DD_EP_10
#    property = DD_EdgePositive_10
#   [../]
#   [./DD_EP_11]
#    type = MaterialRealAux
#    variable = DD_EP_11
#    property = DD_EdgePositive_11
#   [../]
#   [./DD_EP_12]
#    type = MaterialRealAux
#    variable = DD_EP_12
#    property = DD_EdgePositive_12
#   [../]
# []
  

# [AuxVariables]
#   [./DD_EN_01]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EN_02]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EN_03]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EN_04]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EN_05]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EN_06]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EN_07]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EN_08]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EN_09]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EN_10]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EN_11]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_EN_12]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
# []


# [AuxKernels]
#   [./DD_EN_01]
#    type = MaterialRealAux
#    variable = DD_EN_01
#    property = DD_EdgeNegative_1
#   [../]
#   [./DD_EN_02]
#    type = MaterialRealAux
#    variable = DD_EN_02
#    property = DD_EdgeNegative_2
#   [../]
#   [./DD_EN_03]
#    type = MaterialRealAux
#    variable = DD_EN_03
#    property = DD_EdgeNegative_3
#   [../]
#   [./DD_EN_04]
#    type = MaterialRealAux
#    variable = DD_EN_04
#    property = DD_EdgeNegative_4
#   [../]
#   [./DD_EN_05]
#    type = MaterialRealAux
#    variable = DD_EN_05
#    property = DD_EdgeNegative_5
#   [../]
#   [./DD_EN_06]
#    type = MaterialRealAux
#    variable = DD_EN_06
#    property = DD_EdgeNegative_6
#   [../]
#   [./DD_EN_07]
#    type = MaterialRealAux
#    variable = DD_EN_07
#    property = DD_EdgeNegative_7
#   [../]
#   [./DD_EN_08]
#    type = MaterialRealAux
#    variable = DD_EN_08
#    property = DD_EdgeNegative_8
#   [../]
#   [./DD_EN_09]
#    type = MaterialRealAux
#    variable = DD_EN_09
#    property = DD_EdgeNegative_9
#   [../]
#   [./DD_EN_10]
#    type = MaterialRealAux
#    variable = DD_EN_10
#    property = DD_EdgeNegative_10
#   [../]
#   [./DD_EN_11]
#    type = MaterialRealAux
#    variable = DD_EN_11
#    property = DD_EdgeNegative_11
#   [../]
#   [./DD_EN_12]
#    type = MaterialRealAux
#    variable = DD_EN_12 
#    property = DD_EdgeNegative_12
#   [../]
# []

# [AuxVariables]
#   [./DD_SP_01]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SP_02]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SP_03]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SP_04]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SP_05]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SP_06]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SP_07]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SP_08]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SP_09]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SP_10]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SP_11]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SP_12]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
# []


# [AuxKernels]
#   [./DD_SP_01]
#    type = MaterialRealAux
#    variable = DD_SP_01
#    property = DD_ScrewPositive_1
#   [../]
#   [./DD_SP_02]
#    type = MaterialRealAux
#    variable = DD_SP_02
#    property = DD_ScrewPositive_2
#   [../]
#   [./DD_SP_03]
#    type = MaterialRealAux
#    variable = DD_SP_03
#    property = DD_ScrewPositive_3
#   [../]
#   [./DD_SP_04]
#    type = MaterialRealAux
#    variable = DD_SP_04
#    property = DD_ScrewPositive_4
#   [../]
#   [./DD_SP_05]
#    type = MaterialRealAux
#    variable = DD_SP_05
#    property = DD_ScrewPositive_5
#   [../]
#   [./DD_SP_06]
#    type = MaterialRealAux
#    variable = DD_SP_06
#    property = DD_ScrewPositive_6
#   [../]
#   [./DD_SP_07]
#    type = MaterialRealAux
#    variable = DD_SP_07
#    property = DD_ScrewPositive_7
#   [../]
#   [./DD_SP_08]
#    type = MaterialRealAux
#    variable = DD_SP_08
#    property = DD_ScrewPositive_8
#   [../]
#   [./DD_SP_09]
#    type = MaterialRealAux
#    variable = DD_SP_09
#    property = DD_ScrewPositive_9
#   [../]
#   [./DD_SP_10]
#    type = MaterialRealAux
#    variable = DD_SP_10 
#    property = DD_ScrewPositive_10
#   [../]
#   [./DD_SP_11]
#    type = MaterialRealAux
#    variable = DD_SP_11 
#    property = DD_ScrewPositive_11
#   [../]
#   [./DD_SP_12]
#    type = MaterialRealAux
#    variable = DD_SP_12
#    property = DD_ScrewPositive_12
#   [../]
# []
  

# [AuxVariables]
#   [./DD_SN_01]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SN_02]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SN_03]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SN_04]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SN_05]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SN_06]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SN_07]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SN_08]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SN_09]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SN_10]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SN_11]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
#   [./DD_SN_12]
#    order = CONSTANT
#    family = MONOMIAL
#   [../]
# []


# [AuxKernels]
#   [./DD_SN_01]
#    type = MaterialRealAux
#    variable = DD_SN_01
#    property = DD_ScrewNegative_1
#   [../]
#   [./DD_SN_02]
#    type = MaterialRealAux
#    variable = DD_SN_02
#    property = DD_ScrewNegative_2
#   [../]
#   [./DD_SN_03]
#    type = MaterialRealAux
#    variable = DD_SN_03
#    property = DD_ScrewNegative_3
#   [../]
#   [./DD_SN_04]
#    type = MaterialRealAux
#    variable = DD_SN_04
#    property = DD_ScrewNegative_4
#   [../]
#   [./DD_SN_05]
#    type = MaterialRealAux
#    variable = DD_SN_05
#    property = DD_ScrewNegative_5
#   [../]
#   [./DD_SN_06]
#    type = MaterialRealAux
#    variable = DD_SN_06
#    property = DD_ScrewNegative_6
#   [../]
#   [./DD_SN_07]
#    type = MaterialRealAux
#    variable = DD_SN_07
#    property = DD_ScrewNegative_7
#   [../]
#   [./DD_SN_08]
#    type = MaterialRealAux
#    variable = DD_SN_08
#    property = DD_ScrewNegative_8
#   [../]
#   [./DD_SN_09]
#    type = MaterialRealAux
#    variable = DD_SN_09
#    property = DD_ScrewNegative_9
#   [../]
#   [./DD_SN_10]
#    type = MaterialRealAux
#    variable = DD_SN_10
#    property = DD_ScrewNegative_10
#   [../]
#   [./DD_SN_11]
#    type = MaterialRealAux
#    variable = DD_SN_11
#    property = DD_ScrewNegative_11
#   [../]
#   [./DD_SN_12]
#    type = MaterialRealAux #ElementAverageValue
#    variable = DD_SN_12
#    property = DD_ScrewNegative_12
#   [../]
# []

[AuxVariables]
  [./shear_stress_01]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./shear_resistance_01]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./back_stress_01]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./DD_mobile_01]
   order = CONSTANT
   family = MONOMIAL
  [../]
[]


[AuxKernels]
  [./shear_stress_01]
   type = MaterialStdVectorAux
   variable = shear_stress_01
   property = applied_shear_stress  
   index = 0
  [../]
  [./back_stress_01]
   type = MaterialStdVectorAux
   variable = back_stress_01
   property = back_stress
   index = 0
  [../]
  [./shear_resistance_01]
    type = MaterialStdVectorAux
    variable = shear_resistance_01
    property = slip_resistance  
    index = 0
  [../]
  [./DD_mobile_01]
   type = MaterialStdVectorAux
   variable = DD_mobile_01
   property = dislocation_immobile
   index = 0
  [../]
[]

[AuxVariables]
  [./shear_stress_02]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./shear_resistance_02]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./back_stress_02]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./DD_mobile_02]
   order = CONSTANT
   family = MONOMIAL
  [../]
[]

[AuxKernels]
[./shear_stress_02]
   type = MaterialStdVectorAux
   variable = shear_stress_02
   property = applied_shear_stress  
   index = 1
  [../]
[./back_stress_02]
   type = MaterialStdVectorAux
   variable = back_stress_02
   property = back_stress
   index = 1
  [../]
  [./shear_resistance_02]
    type = MaterialStdVectorAux
    variable = shear_resistance_02
    property = slip_resistance  
    index = 1
  [../]
[./DD_mobile_02]
   type = MaterialStdVectorAux
   variable = DD_mobile_02
   property = dislocation_immobile
   index = 1
  [../]
[]

[AuxVariables]
  [./shear_stress_03]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./shear_resistance_03]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./back_stress_03]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./DD_mobile_03]
   order = CONSTANT
   family = MONOMIAL
  [../]
[]

[AuxKernels]
[./shear_stress_03]
   type = MaterialStdVectorAux
   variable = shear_stress_03
   property = applied_shear_stress  
   index = 2
  [../]
[./back_stress_03]
   type = MaterialStdVectorAux
   variable = back_stress_03
   property = back_stress
   index = 2
  [../]
  [./shear_resistance_03]
    type = MaterialStdVectorAux
    variable = shear_resistance_03
    property = slip_resistance  
    index = 2
  [../]
[./DD_mobile_03]
   type = MaterialStdVectorAux
   variable = DD_mobile_03
   property = dislocation_immobile
   index = 2
  [../]
[]

[AuxVariables]
  [./shear_stress_04]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./shear_resistance_04]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./back_stress_04]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./DD_mobile_04]
   order = CONSTANT
   family = MONOMIAL
  [../]
[]

[AuxKernels]
[./shear_stress_04]
   type = MaterialStdVectorAux
   variable = shear_stress_04
   property = applied_shear_stress  
   index = 3
  [../]
[./back_stress_04]
   type = MaterialStdVectorAux
   variable = back_stress_04
   property = back_stress
   index = 3
  [../]
  [./shear_resistance_04]
    type = MaterialStdVectorAux
    variable = shear_resistance_04
    property = slip_resistance  
    index = 3
  [../]
[./DD_mobile_04]
   type = MaterialStdVectorAux
   variable = DD_mobile_04
   property = dislocation_immobile
   index = 3
  [../]
[]

[AuxVariables]
  [./shear_stress_05]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./shear_resistance_05]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./back_stress_05]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./DD_mobile_05]
   order = CONSTANT
   family = MONOMIAL
  [../]
[]

[AuxKernels]
[./shear_stress_05]
   type = MaterialStdVectorAux
   variable = shear_stress_05
   property = applied_shear_stress  
   index = 4
  [../]
[./back_stress_05]
   type = MaterialStdVectorAux
   variable = back_stress_05
   property = back_stress
   index = 4
  [../]
  [./shear_resistance_05]
    type = MaterialStdVectorAux
    variable = shear_resistance_05
    property = slip_resistance  
    index = 4
  [../]
[./DD_mobile_05]
   type = MaterialStdVectorAux
   variable = DD_mobile_05
   property = dislocation_immobile
   index = 4
  [../]
[]

[AuxVariables]
  [./shear_stress_06]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./shear_resistance_06]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./back_stress_06]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./DD_mobile_06]
   order = CONSTANT
   family = MONOMIAL
  [../]
[]

[AuxKernels]
[./shear_stress_06]
   type = MaterialStdVectorAux
   variable = shear_stress_06
   property = applied_shear_stress  
   index = 5
  [../]
[./back_stress_06]
   type = MaterialStdVectorAux
   variable = back_stress_06
   property = back_stress
   index = 5
  [../]
    [./shear_resistance_06]
      type = MaterialStdVectorAux
      variable = shear_resistance_06
      property = slip_resistance
      index = 5
    [../]
[./DD_mobile_06]
   type = MaterialStdVectorAux
   variable = DD_mobile_06
   property = dislocation_immobile
   index = 5
  [../]
[]

[AuxVariables]
  [./shear_stress_07]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./shear_resistance_07]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./back_stress_07]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./DD_mobile_07]
   order = CONSTANT
   family = MONOMIAL
  [../]
[]

[AuxKernels]
[./shear_stress_07]
   type = MaterialStdVectorAux
   variable = shear_stress_07
   property = applied_shear_stress  
   index = 6
  [../]
[./back_stress_07]
   type = MaterialStdVectorAux
   variable = back_stress_07
   property = back_stress
   index = 6
  [../]
    [./shear_resistance_07]
      type = MaterialStdVectorAux
      variable = shear_resistance_07
      property = slip_resistance
      index = 6
    [../]
[./DD_mobile_07]
   type = MaterialStdVectorAux
   variable = DD_mobile_07
   property = dislocation_immobile
   index = 6
  [../]
[]


[AuxVariables]
  [./shear_stress_08]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./shear_resistance_08]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./back_stress_08]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./DD_mobile_08]
   order = CONSTANT
   family = MONOMIAL
  [../]
[]

[AuxKernels]
[./shear_stress_08]
   type = MaterialStdVectorAux
   variable = shear_stress_08
   property = applied_shear_stress  
   index = 7
  [../]
[./back_stress_08]
   type = MaterialStdVectorAux
   variable = back_stress_08
   property = back_stress
   index = 7
  [../]
    [./shear_resistance_08]
      type = MaterialStdVectorAux
      variable = shear_resistance_08
      property = slip_resistance
      index = 7
    [../]
[./DD_mobile_08]
   type = MaterialStdVectorAux
   variable = DD_mobile_08
   property = dislocation_immobile
   index = 7
  [../]
[]



[AuxVariables]
  [./shear_stress_09]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./shear_resistance_09]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./back_stress_09]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./DD_mobile_09]
   order = CONSTANT
   family = MONOMIAL
  [../]
[]

[AuxKernels]
[./shear_stress_09]
   type = MaterialStdVectorAux
   variable = shear_stress_09
   property = applied_shear_stress  
   index = 8
  [../]
[./back_stress_09]
   type = MaterialStdVectorAux
   variable = back_stress_09
   property = back_stress
   index = 8
  [../]
    [./shear_resistance_09]
      type = MaterialStdVectorAux
      variable = shear_resistance_09
      property = slip_resistance
      index = 8
    [../]
[./DD_mobile_09]
   type = MaterialStdVectorAux
   variable = DD_mobile_09
   property = dislocation_immobile
   index = 8
  [../]
[]


[AuxVariables]
  [./shear_stress_10]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./shear_resistance_10]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./back_stress_10]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./DD_mobile_10]
   order = CONSTANT
   family = MONOMIAL
  [../]
[]

[AuxKernels]
[./shear_stress_10]
   type = MaterialStdVectorAux
   variable = shear_stress_10
   property = applied_shear_stress  
   index = 9
  [../]
[./back_stress_10]
   type = MaterialStdVectorAux
   variable = back_stress_10
   property = back_stress
   index = 9
  [../]
    [./shear_resistance_10]
      type = MaterialStdVectorAux
      variable = shear_resistance_10
      property = slip_resistance
      index = 9
    [../]
[./DD_mobile_10]
   type = MaterialStdVectorAux
   variable = DD_mobile_10
   property = dislocation_immobile
   index = 9
  [../]
[]

[AuxVariables]
  [./shear_stress_11]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./shear_resistance_11]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./back_stress_11]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./DD_mobile_11]
   order = CONSTANT
   family = MONOMIAL
  [../]
[]

[AuxKernels]
[./shear_stress_11]
   type = MaterialStdVectorAux
   variable = shear_stress_11
   property = applied_shear_stress  
   index = 10
  [../]
[./back_stress_11]
   type = MaterialStdVectorAux
   variable = back_stress_11
   property = back_stress
   index = 10
  [../]
    [./shear_resistance_11]
      type = MaterialStdVectorAux
      variable = shear_resistance_11
      property = slip_resistance
      index = 10
    [../]
[./DD_mobile_11]
   type = MaterialStdVectorAux
   variable = DD_mobile_11
   property = dislocation_immobile
   index = 10
  [../]
[]

[AuxVariables]
  [./shear_stress_12]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./shear_resistance_12]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./back_stress_12]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./DD_mobile_12]
   order = CONSTANT
   family = MONOMIAL
  [../]
[]

[AuxKernels]
[./shear_stress_12]
   type = MaterialStdVectorAux
   variable = shear_stress_12
   property = applied_shear_stress  
   index = 11
  [../]
[./back_stress_12]
   type = MaterialStdVectorAux
   variable = back_stress_12
   property = back_stress
   index = 11
  [../]
    [./shear_resistance_12]
      type = MaterialStdVectorAux
      variable = shear_resistance_12
      property = slip_resistance
      index = 11
    [../]
[./DD_mobile_12]
   type = MaterialStdVectorAux
   variable = DD_mobile_12
   property = dislocation_immobile
   index = 11
  [../]
[]

[Materials]
  [./compute_stress_wrapper]
    type = ComputeLagrangianWrappedStress
  [../]
[]

[AuxVariables]
  [./GND_density_01]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./GND_density_02]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./GND_density_03]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./GND_density_04]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./GND_density_05]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./GND_density_06]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./GND_density_07]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./GND_density_08]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./GND_density_09]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./GND_density_10]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./GND_density_11]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./GND_density_12]
    order = CONSTANT
    family = MONOMIAL
  [../]
[]

[AuxKernels]
[./GND_density_01]
    type = MaterialStdVectorAux
    variable = GND_density_01
    property = kappa  
    index = 0
  [../]
[./GND_density_02]
    type = MaterialStdVectorAux
    variable = GND_density_02
    property = kappa  
    index = 1
  [../]
[./GND_density_03]
    type = MaterialStdVectorAux
    variable = GND_density_03
    property = kappa  
    index = 2
  [../]
[./GND_density_04]
    type = MaterialStdVectorAux
    variable = GND_density_04
    property = kappa  
    index = 3
  [../]
[./GND_density_05]
    type = MaterialStdVectorAux
    variable = GND_density_05
    property = kappa  
    index = 4
  [../]
[./GND_density_06]
    type = MaterialStdVectorAux
    variable = GND_density_06
    property = kappa  
    index = 5
  [../]
[./GND_density_07]
    type = MaterialStdVectorAux
    variable = GND_density_07
    property = kappa  
    index = 6
  [../]
[./GND_density_08]
    type = MaterialStdVectorAux
    variable = GND_density_08
    property = kappa  
    index = 7
  [../]
[./GND_density_09]
    type = MaterialStdVectorAux
    variable = GND_density_09
    property = kappa  
    index = 8
  [../]
[./GND_density_10]
    type = MaterialStdVectorAux
    variable = GND_density_10
    property = kappa  
    index = 9
  [../]
[./GND_density_11]
    type = MaterialStdVectorAux
    variable = GND_density_11
    property = kappa  
    index = 10
  [../]
[./GND_density_12]
    type = MaterialStdVectorAux
    variable = GND_density_12
    property = kappa  
    index = 11
  [../]
[]

[AuxVariables]
  [./slip_increment_01]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./slip_increment_02]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./slip_increment_03]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./slip_increment_04]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./slip_increment_05]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./slip_increment_06]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./slip_increment_07]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./slip_increment_08]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./slip_increment_09]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./slip_increment_10]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./slip_increment_11]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./slip_increment_12]
    order = CONSTANT
    family = MONOMIAL
  [../]
[]

[AuxKernels]
[./slip_increment_01]
    type = MaterialStdVectorAux
    variable = slip_increment_01
    property = slip_increment  
    index = 0
  [../]
[./slip_increment_02]
    type = MaterialStdVectorAux
    variable = slip_increment_02
    property = slip_increment  
    index = 1
  [../]
[./slip_increment_03]
    type = MaterialStdVectorAux
    variable = slip_increment_03
    property = slip_increment  
    index = 2
  [../]
[./slip_increment_04]
    type = MaterialStdVectorAux
    variable = slip_increment_04
    property = slip_increment  
    index = 3
  [../]
[./slip_increment_05]
    type = MaterialStdVectorAux
    variable = slip_increment_05
    property = slip_increment  
    index = 4
  [../]
[./slip_increment_06]
    type = MaterialStdVectorAux
    variable = slip_increment_06
    property = slip_increment  
    index = 5
  [../]
[./slip_increment_07]
    type = MaterialStdVectorAux
    variable = slip_increment_07
    property = slip_increment  
    index = 6
  [../]
[./slip_increment_08]
    type = MaterialStdVectorAux
    variable = slip_increment_08
    property = slip_increment  
    index = 7
  [../]
[./slip_increment_09]
    type = MaterialStdVectorAux
    variable = slip_increment_09
    property = slip_increment  
    index = 8
  [../]
[./slip_increment_10]
    type = MaterialStdVectorAux
    variable = slip_increment_10
    property = slip_increment  
    index = 9
  [../]
[./slip_increment_11]
    type = MaterialStdVectorAux
    variable = slip_increment_11
    property = slip_increment  
    index = 10
  [../]
[./slip_increment_12]
    type = MaterialStdVectorAux
    variable = slip_increment_12
    property = slip_increment  
    index = 11
  [../]
[]

[AuxVariables]
  [./dde_vel_01]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./dde_vel_02]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./dde_vel_03]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./dde_vel_04]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./dde_vel_05]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./dde_vel_06]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./dde_vel_07]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./dde_vel_08]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./dde_vel_09]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./dde_vel_10]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./dde_vel_11]
   order = CONSTANT
   family = MONOMIAL
  [../]
  [./dde_vel_12]
    order = CONSTANT
    family = MONOMIAL
  [../]
[]

[AuxKernels]
[./dde_vel_01]
    type = MaterialStdVectorAux
    variable = dde_vel_01
    property = dislo_velocity_edge  
    index = 0
  [../]
[./dde_vel_02]
    type = MaterialStdVectorAux
    variable = dde_vel_02
    property = dislo_velocity_edge  
    index = 1
  [../]
[./dde_vel_03]
    type = MaterialStdVectorAux
    variable = dde_vel_03
    property = dislo_velocity_edge  
    index = 2
  [../]
[./dde_vel_04]
    type = MaterialStdVectorAux
    variable = dde_vel_04
    property = dislo_velocity_edge  
    index = 3
  [../]
[./dde_vel_05]
    type = MaterialStdVectorAux
    variable = dde_vel_05
    property = dislo_velocity_edge  
    index = 4
  [../]
[./dde_vel_06]
    type = MaterialStdVectorAux
    variable = dde_vel_06
    property = dislo_velocity_edge  
    index = 5
  [../]
[./dde_vel_07]
    type = MaterialStdVectorAux
    variable = dde_vel_07
    property = dislo_velocity_edge  
    index = 6
  [../]
[./dde_vel_08]
    type = MaterialStdVectorAux
    variable = dde_vel_08
    property = dislo_velocity_edge  
    index = 7
  [../]
[./dde_vel_09]
    type = MaterialStdVectorAux
    variable = dde_vel_09
    property = dislo_velocity_edge  
    index = 8
  [../]
[./dde_vel_10]
    type = MaterialStdVectorAux
    variable = dde_vel_10
    property = dislo_velocity_edge  
    index = 9
  [../]
[./dde_vel_11]
    type = MaterialStdVectorAux
    variable = dde_vel_11
    property = dislo_velocity_edge  
    index = 10
  [../]
[./dde_vel_12]
    type = MaterialStdVectorAux
    variable = dde_vel_12
    property = dislo_velocity_edge  
    index = 11
  [../]
[]

bmurgas-LANL · 2024-10-31T19:35:12Z

bmurgas-LANL
Oct 31, 2024
Author

Hi @GiudGiud,

I have been trying to find the problem with the Material Properties in the InterfaceDiffusion kernel.

The values were set to zero because I used the wrong values when initializing the variables in DiscoFluxCPStressUpdate::initQpStatefulProperties().
I initialize the variables with the correct values and I can read them in the ConservativeAdvection and InterfaceDiffusion kernels. However, the Material Properties are evolving every time step in the crystal plasticity kernel, when reading the updated values I get the correct values in the ConservativeAdvection kernel but I keep getting the initial values in the InterfaceDiffusion.
I suspect this error is coming from the fact that Moose creates three copies of the Material Properties (https://mooseframework.inl.gov/syntax/Materials/index.html). After reading this, I started to check the memory address of the same variable in the three kernels: DiscoFluxCPStressUpdate, ConservativeAdvection and InterfaceDiffusion.
This is the result from lldb:
In crystal plasticity

(lldb) frame var -L _dislo_velocity_edge
0x00007fd001ec6470: (MaterialProperty<std::vector<double> > &) _dislo_velocity_edge = 0x000060000123c090: {
0x000060000123c090:   MaterialPropertyBase<std::__1::vector<double, std::__1::allocator<double> >, false> = {
0x000060000123c090:     PropertyValue = {
0x000060000123c098:       _id = 23
    }
0x000060000123c0a0:     _value = {
0x000060000123c0a0:       _data_ptr = size=12 {
0x000060000123c0a0:         pointer = 0x00006000027715f8 size=12
      }
0x000060000123c0a8:       _data = 0x00006000027715f8 size=12
0x000060000123c0b0:       _size = 8
0x000060000123c0b4:       _allocated_size = 8
    }
  }
}

In interfacediffusiondislocation

(lldb) frame var -L _dislo_velocity_CP_edge
0x00007fd001f81ff0: (const MaterialProperty<std::vector<double> > &) _dislo_velocity_CP_edge = 0x000060000123e970: {
0x000060000123e970:   MaterialPropertyBase<std::__1::vector<double, std::__1::allocator<double> >, false> = {
0x000060000123e970:     PropertyValue = {
0x000060000123e978:       _id = 23
    }
0x000060000123e980:     _value = {
0x000060000123e980:       _data_ptr = size=12 {
0x000060000123e980:         pointer = 0x0000600003765188 size=12
      }
0x000060000123e988:       _data = 0x0000600003765188 size=12
0x000060000123e990:       _size = 4
0x000060000123e994:       _allocated_size = 4
    }
  }
}

And in conservativeadvectiondislocation

(lldb) frame var -L _dislo_velocity_CP_edge
0x00007fd001f57e30: (const MaterialProperty<std::vector<double> > &) _dislo_velocity_CP_edge = 0x000060000123c090: {
0x000060000123c090:   MaterialPropertyBase<std::__1::vector<double, std::__1::allocator<double> >, false> = {
0x000060000123c090:     PropertyValue = {
0x000060000123c098:       _id = 23
    }
0x000060000123c0a0:     _value = {
0x000060000123c0a0:       _data_ptr = size=12 {
0x000060000123c0a0:         pointer = 0x00006000027715f8 size=12
      }
0x000060000123c0a8:       _data = 0x00006000027715f8 size=12
0x000060000123c0b0:       _size = 8
0x000060000123c0b4:       _allocated_size = 8
    }
  }
}

As you can see the memory address of the variable in the crystal plasticity (DiscoFluxCPStressUpdate) and the advection kernel (ConservativeAdvection) is the same. And the address of the interface diffusion kernel is different. Does it means that in the interface kernel system we are using the element face material properties?
I wonder is there is an issue with Moose when the material properties are being updated in the crystal plasticity model and the elemental material variables don't update.
I found this discussion #19062 and found it helpful. I'll try to create a Material Object to see if when calling the function copmputeQpProperties both the element and element face material properties are updated with the correct values.
Thank you for your help.
Brayan

3 replies

GiudGiud Nov 1, 2024
Collaborator

hi

sorry for the delay. I ll get to this this weekend. I tried to find a resource to help earlier but failed.

GiudGiud Nov 2, 2024
Collaborator

I can't run the input file, I need the px3_bar_hex8_cl16e-3.e mesh and the px3-bar-sp10-sp01.inp data file.

I ll make a MWE to see if I can reproduce this without those

GiudGiud Nov 3, 2024
Collaborator

ok I made do with sxori-0-0-0.inp instead of px3-bar-sp10-sp01.inp.

I am not getting uniformly 0 values on my case. So it could have something to do with your mesh. Is it possible for you to share that with me?
Alternatively, maybe I am not checking the materials the same way as you. How did you determine slip_resistance /... were uniformly 0 in the interface kernels ?

And the address of the interface diffusion kernel is different. Does it means that in the interface kernel system we are using the element face material properties?

Yes that is expected. The interface kernels request "boundary material data" from their inheritance on TwoMaterialPropertyInterface.

InputParameters
TwoMaterialPropertyInterface::validParams()
{

  // Objects inheriting from TwoMaterialPropertyInterface rely on Boundary MaterialData
  InputParameters params = MaterialPropertyInterface::validParams();
  params.set<Moose::MaterialDataType>("_material_data_type") = Moose::BOUNDARY_MATERIAL_DATA;
  return params;
}

Which under the hood maps to boundary material data

MaterialData &
FEProblemBase::getMaterialData(Moose::MaterialDataType type, const THREAD_ID tid)
{
  switch (type)
  {
    case Moose::BLOCK_MATERIAL_DATA:
      return _material_props.getMaterialData(tid);
    case Moose::NEIGHBOR_MATERIAL_DATA:
      return _neighbor_material_props.getMaterialData(tid);
    case Moose::BOUNDARY_MATERIAL_DATA:
    case Moose::FACE_MATERIAL_DATA:
    case Moose::INTERFACE_MATERIAL_DATA:
      return _bnd_material_props.getMaterialData(tid);
  }

So it's a different address than the kernels which use block material data, which is from regular materials.

bmurgas-LANL · 2025-11-13T21:46:42Z

bmurgas-LANL
Nov 13, 2025
Author

Hi @GiudGiud,

We found out after a while that the material properties were not correct because in the Crystal Plasticity class (ComputeMultipleCrystalPlasticityStress.C) there is this line at the beginning of the updateStress function:

// Does not support face/boundary material property calculation
  if (isBoundaryMaterial())
  return;

I commented this and the Crystal Plasticity code is working but I'm not sure if that'll generate trouble in the convergence later. I'll share my input in case I see that.

Thanks again for your help!
Brayan

8 replies

GiudGiud Nov 19, 2025
Collaborator

do you know of an alternative method to access the materialProperties in the BCs and InterfaceKernels systems?

getMaterialsPorperty is all I know, and when called on an object executing on a boundary it gets you the face Qp material property values.
Are you thinking of using element properties instead?
Or lagging by a step (by projecting onto an auxvariable) ?

bmurgas-LANL Nov 26, 2025
Author

Hi @GiudGiud,

I think element properties work too. I read other discussions and found this one: #27879
Do you think Functor Materials would work for this case?

I tried to create a Material object but I have the same issue, the values of the Material Property are not defined in the boundary/interface elements.
Also, created a Interface Material object and got this error:

*** ERROR ***
/Users/bmurgas/workspace/simus/discoflux/cu_bx_exp/micropillar-base/sx-base/base-1-debug/micropillar-base-1-sx.i:2956.3:
The following occurred in the MaterialBase 'DislocVelocity_neighbor' of type DislocationVelocity.

MaterialBase::subdomainSetup in Material'DislocVelocity_neighbor' needs to be implemented

Thanks for your help,
Brayan

GiudGiud Nov 26, 2025
Collaborator

I think element properties work too.

my question wuld be, which quadrature point from the element would we use for evaluating the material property when at the interface?
Seems like we need to use an element side quadrature point.

Do you think Functor Materials would work for this case?

FunctorMaterials are a whole different system for materials. You will need to adapt any kernel consuming these properties for full compatibility. If you are open to writing the relevant calls in the interface kernel to use a functor material, then please feeel free to try

Also, created a Interface Material object and got this error:

You mean you created a derived class from interface material and got an error, was that error from the interface kernel or from using it in a different kernel?

You can implement subdomain_setup. It's likely just doing nothing in your case.

void
NewMaterial::subdomainSetup() override {}
```

bmurgas-LANL Dec 2, 2025
Author

my question wuld be, which quadrature point from the element would we use for evaluating the material property when at the interface?
Seems like we need to use an element side quadrature point.

I need the value of a quadrature point close to the interface. I am using hex 8 elements so the values may change at the 4 quadrature points.

You mean you created a derived class from interface material and got an error, was that error from the interface kernel or from using it in a different kernel?
You can implement subdomain_setup. It's likely just doing nothing in your case.

I tried to do something similar to SideSetHeatTransferMaterial (https://mooseframework.inl.gov/source/materials/SideSetHeatTransferMaterial.html). The properties are computed in SideSetHeatTransferMaterial and then used in the SideSetHeatTransferKernel.
I first tried to use the Material parameters
InputParameters params = Material::validParams();
and
Material(parameters),
But I have the same issue that the properties were zero in the InterfaceKernel and BCs. So I noticed SideSetHeatTransferMaterial would be something similar to what I need in order to use the material properties in the interface kernel or bcs but it looks like it's not working.

I have a question regarding this line of code

// Does not support face/boundary material property calculation
  if (isBoundaryMaterial())
  return;

If I bypass this condition, does it mean the residuals and jacobians are considered twice: One for the regular materials and other for the interface materials?

I appreciate your help!
Thank you

GiudGiud Dec 2, 2025
Collaborator

I need the value of a quadrature point close to the interface. I am using hex 8 elements so the values may change at the 4 quadrature points.

Right. I think there is indeed an element Qp that is closer to the side Qp being considered. It's still something we have not implemented, preferring to compute the value of the material property at the side Qp.

But I have the same issue that the properties were zero in the InterfaceKernel and BCs

In an interface kernel, this does not seem right.
In the BC I could see it, with InterfaceMaterial not being computed outside of interfaces

If I bypass this condition, does it mean the residuals and jacobians are considered twice: One for the regular materials and other for the interface materials?

The residual (and Jacobian) are computed in an element loop across the entire mesh, where for each element, first the element Qps are considered, and then its sides (for DG, for boundary conditions, interface kernels).
The material properties for element Qps are computed before the element contributions (from kernels), and for side Qps before the contributions from the sides

Problems reading material property in InterfaceKernel #28908

Uh oh!

bmurgas-LANL Oct 22, 2024

Check these boxes if you have followed the posting rules.

Question

Replies: 3 comments · 13 replies

Uh oh!

GiudGiud Oct 23, 2024 Collaborator

Uh oh!

bmurgas-LANL Oct 23, 2024 Author

Uh oh!

bmurgas-LANL Oct 23, 2024 Author

Uh oh!

bmurgas-LANL Oct 31, 2024 Author

Uh oh!

GiudGiud Nov 1, 2024 Collaborator

Uh oh!

Uh oh!

GiudGiud Nov 2, 2024 Collaborator

Uh oh!

Uh oh!

GiudGiud Nov 3, 2024 Collaborator

Uh oh!

Uh oh!

bmurgas-LANL Nov 13, 2025 Author

Uh oh!

GiudGiud Nov 19, 2025 Collaborator

Uh oh!

bmurgas-LANL Nov 26, 2025 Author

Uh oh!

GiudGiud Nov 26, 2025 Collaborator

Uh oh!

bmurgas-LANL Dec 2, 2025 Author

Uh oh!

GiudGiud Dec 2, 2025 Collaborator

bmurgas-LANL
Oct 22, 2024

Replies: 3 comments 13 replies

GiudGiud
Oct 23, 2024
Collaborator

bmurgas-LANL Oct 23, 2024
Author

bmurgas-LANL Oct 23, 2024
Author

bmurgas-LANL
Oct 31, 2024
Author

GiudGiud Nov 1, 2024
Collaborator

GiudGiud Nov 2, 2024
Collaborator

GiudGiud Nov 3, 2024
Collaborator

bmurgas-LANL
Nov 13, 2025
Author

GiudGiud Nov 19, 2025
Collaborator

bmurgas-LANL Nov 26, 2025
Author

GiudGiud Nov 26, 2025
Collaborator

bmurgas-LANL Dec 2, 2025
Author

GiudGiud Dec 2, 2025
Collaborator