Subdomain do not exist

[Emmanuel-324]

Below is a snippet of my mesh block

[Mesh]
   [file]
      type = FileMeshGenerator
      file = my_kks_noload_twophase_out.e-s180
      use_for_exodus_restart = true
      add_subdomain_ids = '1'
      add_subdomain_names = file
   []
[]

Below is the error I get I have tried adding the block id but it does not work

The following error occurred in the MooseMesh 'mesh' of type MeshGeneratorMesh.check_inp(moose_srv)
Requested subdomain 1 does not exist.

[GiudGiud]

Hello

The parameter is not at the right place.

Add_subdomain_ should be under the Mesh block directly

[GiudGiud]

in that input you have several blocks that specify block = 1

if you do not need a block of id 1, then delete all those

[Emmanuel-324]

I get this error when I take out the block 1 under the Physics/Solid Mechanics block

No SolidMechanics/QuasiStatic action can be block unrestricted if more than one SolidMechanics/QuasiStatic action is specified.

[GiudGiud]

you have to take out the entire block, not just block=1

[Emmanuel-324]

So I am trying to use the crystal plasticity model to simulate an already generated kks model having two parameters eta1 and eta2. I use to use the monolithic crystal plasticity model in the MOOSE modules to simulate as below as shown below

[Mesh]
  [file]
    type = FileMeshGenerator
    file = kks_multiphase_largescale_large_out.e-s183
    use_for_exodus_restart = true
  []
[]

[Variables]
  [disp_x]
  []
  [disp_y]
  []

  # order parameter 0
  [eta0]
    initial_from_file_var = eta1
  []
  # order parameter 1
  [eta1]
    initial_from_file_var = eta2
  []
  [eta2]
    initial_from_file_var = eta3
  []
[]

[AuxVariables]
  [vonmises]
    order = CONSTANT
    family = MONOMIAL
  []
  [stress_xx]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
[]

[AuxKernels]
  [vonmises]
    type = RankTwoScalarAux
    rank_two_tensor = stress
    variable = vonmises
    scalar_type = VonMisesStress
    execute_on = timestep_end
    #   block = 0
  []
  [stress_xx]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_xx
    index_j = 0
    index_i = 0
    execute_on = timestep_end
    block = 0
  []
[]

[BCs]
  [symmy]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0
  []
  [symmx]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0
  []
  [tdisp]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = right
    function = '0.1*t'
  []
[]

[Materials]
  [elasticity_tensor_phase0]
    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
    base_name = phase0
    block = 0
  []
  [crysp_phase0]
    type = FiniteStrainCPSlipRateRes
    slip_sys_file_name = input_slip_sys.txt
    nss = 12
    num_slip_sys_flowrate_props = 2 #Number of properties in a slip system
    flowprops = '1 4 1e-6 0.05 5 8 1e-6 0.05 9 12 1e-6 0.05'
    hprops = '1.0 215 78.6 200.0 1.75'
    gprops = '1 4 78.6 5 8 78.6 9 12 78.6'
    slip_incr_tol = 1
    maximum_substep_iteration = 12
    rtol = 1e-9
    abs_tol = 1e-10
    base_name = phase0
  []
  [strain_phase0]
    type = ComputeFiniteStrain
    displacements = 'disp_x disp_y'
    base_name = phase0
    #    eigenstrain_names = eigenstrain2
  []

  [elasticity_tensor_phase1]
    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
    euler_angle_1 = 0.0
    euler_angle_2 = 45.0
    euler_angle_3 = 0.0
    base_name = phase1
  []
  [crysp_phase1]
    type = FiniteStrainCPSlipRateRes
    slip_sys_file_name = input_slip_sys.txt
    nss = 12
    num_slip_sys_flowrate_props = 2 #Number of properties in a slip system
    flowprops = '1 4 1e-6 0.05 5 8 1e-6 0.05 9 12 1e-6 0.05'
    hprops = '1.0 250 118.6 230.0 1.75'
    gprops = '1 4 118.6 5 8 118.6 9 12 118.6'
    slip_incr_tol = 1
    maximum_substep_iteration = 12
    rtol = 1e-9
    abs_tol = 1e-10
    base_name = phase1
  []
  [strain_phase1]
    type = ComputeFiniteStrain
    displacements = 'disp_x disp_y'
    base_name = phase1
  []

  [elasticity_tensor_phase2]
    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
    euler_angle_1 = 0.0
    euler_angle_2 = 45.0
    euler_angle_3 = 30.0
    base_name = phase2
  []
  [crysp_phase2]
    type = FiniteStrainCPSlipRateRes
    slip_sys_file_name = input_slip_sys.txt
    nss = 12
    num_slip_sys_flowrate_props = 2 #Number of properties in a slip system
    flowprops = '1 4 1e-6 0.05 5 8 1e-6 0.05 9 12 1e-6 0.05'
    hprops = '1.0 250 118.6 230.0 1.75'
    gprops = '1 4 118.6 5 8 118.6 9 12 118.6'
    slip_incr_tol = 1
    maximum_substep_iteration = 12
    rtol = 1e-9
    abs_tol = 1e-10
    base_name = phase2
  []
  [strain_phase2]
    type = ComputeFiniteStrain
    displacements = 'disp_x disp_y'
    base_name = phase2
  []

  # Switching functions for each phase
  [h0]
    type = SwitchingFunctionMultiPhaseMaterial
    phase_etas = eta0
    all_etas = 'eta0 eta1 eta2'
    h_name = h0
  []
  [h1]
    type = SwitchingFunctionMultiPhaseMaterial
    phase_etas = eta1
    all_etas = 'eta0 eta1 eta2'
    h_name = h1
  []
  [h2]
    type = SwitchingFunctionMultiPhaseMaterial
    phase_etas = eta2
    all_etas = 'eta0 eta1 eta2'
    h_name = h2
  []

  # Generate the global stress from the phase stresses
  [global_stress]
    type = MultiPhaseStressMaterial
    phase_base = 'phase0 phase1 phase2'
    h = 'h0     h1     h2'
  []
[]

[Kernels]
  [eta0_dt]
    type = TimeDerivative
    variable = eta0
  []
  [eta1_dt]
    type = TimeDerivative
    variable = eta1
  []
  [eta2_dt]
    type = TimeDerivative
    variable = eta2
  []

  [TensorMechanics]
    displacements = 'disp_x disp_y'
  []
[]

[Postprocessors]
  [vonmises]
    type = ElementAverageValue
    variable = vonmises
    block = 'ANY_BLOCK_ID 0'
  []
  [stress_xx]
    type = ElementAverageValue
    variable = stress_xx
    block = 'ANY_BLOCK_ID 0'
  []
[]

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

[Executioner]
  type = Transient
  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -snes_type'
  petsc_options_value = 'lu            superlu_dist          vinewtonrsls'

  #  petsc_options_iname = '-pc_type -pc_asm_overlap -sub_pc_type -ksp_type -ksp_gmres_restart'
  #  petsc_options_value = ' asm      2              lu            gmres     200'
  l_max_its = 20
  nl_max_its = 10
  nl_rel_tol = 1.0e-8
  nl_abs_tol = 1.0e-9

  end_time = 300

  [TimeStepper]
    type = IterationAdaptiveDT
    dt = 5e-3
    cutback_factor = 0.75
    growth_factor = 1.2
    optimal_iterations = 20
  []
[]

[Outputs]
  exodus = true
  csv = true
[]

The input script above is now deprecated and I have been trying to find the suitable kernel in the stress_update_method(modules/solid_mechanics/test/tests/crystal_plasticity) to do the same but I run into a lot materials property errors.

Is there one in there you can suggest to help with my issue where it can easily take the two or more order parameters. Note that is what I have been trying to do with bicrytsal but I get the domain issues in the already generated mesh.

[Emmanuel-324]

Thanks, my model works now.