Grain Boundary interface PolycrystalVoronoi

[selarem]

Hello,
We are making CP models and using Polycrystalvoronoi for tesselation.
My question is how to get the GB interface name to use it with breakmeshbyblockgenerator to create a CZM at GB.
Thanks for any help,
Saber

[GiudGiud]

HellHe

Is the GB interface computed by the simulation?
If so you ll have to do this in two steps, saving the mesh in a checkpoint file then loading it using a FileMeshGenerator, then a BMBB can act on it

[ghost]

Hi, is it possible to obtain the interface from a test like this (attached below)? Because even in the output through paraview or other, it is possible to see the different values of the euler angles and elasticity tensor components but we cannot extrude each 'grain' as a block. If not possible, is there any other way to do so using Moose? Thanks

[Mesh]
    [top_half]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 120
    ny = 60
    ymin = 0
    ymax = 0.5
    elem_type = tri3
    boundary_id_offset = 0
    boundary_name_prefix = top_half
  []
  [top_stitch]
    type = BoundingBoxNodeSetGenerator
    input = top_half
    new_boundary = top_stitch
    bottom_left = '0.5 0 0'
    top_right = '1 0 0'
  []
  [bottom_half]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 120
    ny = 60
    ymin = -0.5
    ymax = 0
    elem_type = tri3
    boundary_id_offset = 5
    boundary_name_prefix = bottom_half
  []
  [bottom_stitch]
    type = BoundingBoxNodeSetGenerator
    input = bottom_half
    new_boundary = bottom_stitch
    bottom_left = '0.5 0 0'
    top_right = '1 0 0'
  []
  [stitch]
    type = StitchedMeshGenerator
    inputs = 'top_stitch bottom_stitch'
    stitch_boundaries_pairs = 'top_stitch bottom_stitch'
  []
   [./corner_node]
    type = ExtraNodesetGenerator
    new_boundary = BL
    coord = '0 0 0' 
    input = stitch
  [../]

  construct_side_list_from_node_list = true
[]


[GlobalParams]
  displacements = 'disp_x disp_y'
[]
[Variables]
  [./disp_x]
  [../]
  [./disp_y]
  [../]
[]

[Physics/SolidMechanics/QuasiStatic]
      [./All]
        add_variables = true
        strain = FINITE
        additional_generate_output = 'stress_xx  stress_xy  stress_yy'
        planar_formulation = PLANE_STRAIN
        save_in = 'resid_x resid_y'
        decomposition_method = HughesWinget
      [../]
[]

[AuxVariables]
  [./bounds_dummy]
  [../]
 [./C1111]
    order = CONSTANT
    family = MONOMIAL
  [../]
[./C1212]
    order = CONSTANT
    family = MONOMIAL
  [../]
[./C3131]
    order = CONSTANT
    family = MONOMIAL
  [../]
[./C1122]
    order = CONSTANT
    family = MONOMIAL
  [../]
[./euler1]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./euler2]
    order = CONSTANT
    family = MONOMIAL
  [../]
 [./euler3]
    order = CONSTANT
    family = MONOMIAL
  [../]
[./resid_x]
  [../]
  [./resid_y]
  [../]
[]

[BCs]
    [./bottomy]
    type = DirichletBC
    variable = disp_x
    boundary = BL
    value = 0.0
    preset = false
  [../]
  [./topx]
    type = FunctionDirichletBC
    variable = disp_y
    boundary = top_half_top
    function ='1e-5*t'
   preset = false
   [../]
   [./bottomx]
    type = DirichletBC
    variable = disp_y
    boundary = bottom_half_bottom
    value = 0.0
    preset = false
  [../]

[]

[AuxKernels]
[./C1111]
    type = RankFourAux
    variable = C1111
    rank_four_tensor = elasticity_tensor
    index_l = 0
    index_j = 0
    index_k = 0
    index_i = 0
    execute_on = timestep_end
  [../]
[./C1122]
    type = RankFourAux
    variable = C1122
    rank_four_tensor = elasticity_tensor
    index_l = 0
    index_j = 0
    index_k = 1
    index_i = 1
    execute_on = timestep_end
  [../]
[./C1212]
    type = RankFourAux
    variable = C1212
    rank_four_tensor = elasticity_tensor
    index_l = 0
    index_j = 1
    index_k = 0
    index_i = 1
    execute_on = timestep_end
  [../]
[./C3131]
    type = RankFourAux
    variable = C3131
    rank_four_tensor = elasticity_tensor
    index_l = 2
    index_j = 0
    index_k = 2
    index_i = 0
    execute_on = timestep_end
  [../]
[./euler1]
    type = MaterialRealVectorValueAux
    variable = euler1
    property = Euler_angles
    component = 0
    execute_on = timestep_end
  [../]
  [./euler2]
    type = MaterialRealVectorValueAux
    variable = euler2
    property = Euler_angles
    component = 1
    execute_on = timestep_end
  [../]
  [./euler3]
    type = MaterialRealVectorValueAux
    variable = euler3
    property = Euler_angles
    component = 2
    execute_on = timestep_end
[../]
[]


[Materials]
  [./elasticity_tensor]
    type = ComputeElasticityTensorCP
    C_ijkl = '48.7e3  11.9e3  11.9e3  48.7e3 11.9e3  48.7e3   13.1e3 13.1e3 13.1e3'
    fill_method = symmetric9
    read_prop_user_object = prop_read
  [../]
  [stress_copper]
    type = ComputeMultipleCrystalPlasticityStress
    crystal_plasticity_models = 'trial_xtalpl_copper'
    tan_mod_type = exact
    maximum_substep_iteration = 5
    rtol = 1e-6 # Constitutive stress residual relative tolerance
    maxiter_state_variable = 10 # Maximum number of iterations for stress update
    use_line_search = true
  []
  [trial_xtalpl_copper]
    type = CrystalPlasticityKalidindiUpdate
    number_slip_systems = 12
    slip_sys_file_name = 12input_slip_sys.txt
    ao = 2e-3
    gss_a = 2
    gss_initial = '8'
    h = 100
    r = 1.4
    xm = 0.1
    t_sat = 80
    []
[]

[UserObjects]
  [./prop_read]
    type = PropertyReadFile
    prop_file_name = 'input_euler.txt'
    # Enter file data as prop#1, prop#2, .., prop#nprop
    nprop = 3
    read_type = grain
    #ngrain = 20
    nvoronoi = 4
    use_random_voronoi = true
    rand_seed = 777
    rve_type = periodic
  [../]
[]


[Executioner]
 type = Transient
 petsc_options_iname = '-pc_type -pc_hypre_type -ksp_gmres_restart -snes_type'
 petsc_options_value = 'hypre boomeramg 31 vinewtonrsls'
 solve_type = 'newton'
 #scheme=crank-nicolson
 automatic_scaling = true
 dtmin = 5e-12
 dt = 0.5
 nl_abs_tol = 1e-6
 nl_rel_tol = 1e-6
 nl_max_its = 20
 end_time = 2e6

[]



[Postprocessors]
    [./top_x]
    type = NodalExtremeValue
    variable = disp_x
    boundary = top_half_top
    value_type = max
  [../]

[./top_y]
    type = NodalExtremeValue
    variable = disp_y
    boundary = top_half_top
    value_type = max
  [../]
  [./resid_x]
    type = NodalSum
    variable = resid_x
    boundary = bottom_half_bottom
  [../]
  [./resid_y]
    type = NodalSum
    variable = resid_y
    boundary = bottom_half_bottom
  [../]

[]

[Outputs]
  exodus = true
  time_step_interval = 2
  [./table]
    time_step_interval = 1
    type = CSV
   # delimiter = ' '
  [../]
[]

[GiudGiud]

I dont think as a single input it can currently be done.

If not possible, is there any other way to do so using Moose?

You can work to implement MeshModifiers that do the separation of the grain and the breaking-the-mesh

[ghost]

Are there any meshmodifiers on moose we can base ourselves on to implement this ?

[GiudGiud]

Yes see
https://mooseframework.inl.gov/source/meshmodifiers/SidesetAroundSubdomainUpdater.html
for creating / updating a sideset