error when use CoupledVarThresholdElementSubdomainModifier and mesh adaptivity at the same time

[PEI0214]

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
When I used CoupledVarThresholdElementSubdomainModifier and mesh adaptivity at the same time , appeared the following error:

libMesh terminating:
ExodusII_IO_Helper: block id 2 not found in block_ids.
[0] ../src/mesh/exodusII_io_helper.C, line 3266, compiled Aug 20 2024 at 14:09:52

Abort(1) on node 0 (rank 0 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

The input file is as follows:

test.i

# primary model units (m | s | Pa | kg --> N/m²)
[Problem]
  kernel_coverage_check = false
  material_coverage_check = false
  boundary_restricted_node_integrity_check = false
[]

new_domain_blocks = '1001	1002'

[Mesh]
  [BaseMesh]
    type = GeneratedMeshGenerator
    dim = 3
    nx = 20
    ny = 20
    nz = 2
    xmin = -5
    xmax = +5
    ymin = -5
    ymax = +5
    zmin = -2
    zmax = +2
  []
  [Box1]
    type = SubdomainBoundingBoxGenerator
    input = "BaseMesh"
    block_id = 1
    location = "INSIDE"
    bottom_left = "-10 -10 -2"
    top_right   = "+10 +10 +2"
  []
  [Box2]
    type = SubdomainBoundingBoxGenerator
    input = "Box1"
    block_id = 2
    location = "INSIDE"
    bottom_left = "-10 -10 0"
    top_right   = "+10 +10 +2"
  []

  add_subdomain_ids = ${new_domain_blocks}

[]

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

[Variables]
  [disp_x]
    block = '${new_domain_blocks}'
  []
  [disp_y]
    block = '${new_domain_blocks}'
  []
  [disp_z]
    block = '${new_domain_blocks}'
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  [./all]
    strain = finite
    add_variables = true
    block = '${new_domain_blocks}'
    generate_output = 'stress_xx stress_zz strain_xx strain_zz'
  [../]
[]

# ===== Gravity =====
[Kernels]
  [Gravity]
    type = Gravity
    block = '${new_domain_blocks}'
    use_displaced_mesh = false
    variable = disp_z
    value = -10
  []
[]

# ===== Boundary Conditions =====
[BCs]
  [archor_x]
    type = DirichletBC
    boundary = 'back'
    variable = disp_x
    value = 0
  []

  [archor_y]
    type = DirichletBC
    boundary = 'back'
    variable = disp_y
    value = 0
  []

  [archor_z]
    type = DirichletBC
    boundary = 'back'
    variable = disp_z
    value = 0
  []
[]

# ===== Materials (linear-elastic to keep it simple) =====
[Materials]
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    block = '${new_domain_blocks}'
    youngs_modulus = 50E6 # 50 MPa
    poissons_ratio = 0.3
  []

  [stress]
    type = ComputeFiniteStrainElasticStress
    block = '${new_domain_blocks}'
  []

  [density]
    type = GenericConstantMaterial
    block = '${new_domain_blocks}'
    prop_names = density
    prop_values = 2000
  []
[]

# ===== Block Activation/Deactivation =====
[Functions]
  [Box_ActivationStateFunction]
    type = PiecewiseConstant
    direction = left
    xy_data = '0   2
               100 2
               200 1
               300 1
               400 0
               800 0'
  []
[]

[AuxVariables]
  [Box_Activation_AuxVariable]
    family = MONOMIAL
    order = CONSTANT
  []
[]
# link Box2_ActivationStateFunction and Box2_Activation_AuxVariable
[AuxKernels]
  [Box_Activation_FunctionAux]
    type = FunctionAux
    variable = 'Box_Activation_AuxVariable'
    function = 'Box_ActivationStateFunction'
    execute_on = 'INITIAL TIMESTEP_BEGIN'
  []
[]
# move Box2 between subdomains
[MeshModifiers]
  [Box1_SubdomainModifier]
    type = CoupledVarThresholdElementSubdomainModifier
    coupled_var = 'Box_Activation_AuxVariable'
    criterion_type = ABOVE
    threshold = 1.5
    block = 1
    subdomain_id = 1
    complement_subdomain_id = 1001
    execute_on = 'INITIAL TIMESTEP_BEGIN'
  []
  [Box2_SubdomainModifier]
    type = CoupledVarThresholdElementSubdomainModifier
    coupled_var = 'Box_Activation_AuxVariable'
    criterion_type = ABOVE
    threshold = 0.5
    block = 2
    subdomain_id = 2
    complement_subdomain_id = 1002
    execute_on = 'INITIAL TIMESTEP_BEGIN'
  []
[]

[Adaptivity]
  marker = marker
  max_h_level = 2
  initial_steps = 1
  stop_time = 200
  [Markers]
    [marker]
      type = UniformMarker
      mark = REFINE
      block = '2'
    []
  []
[]

# ===== Executioner =====
[Executioner]
  type = Transient
  automatic_scaling = true

  end_time = 800
  dt = 100

  solve_type = 'NEWTON'
  petsc_options = '-snes_converged_reason'
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = ' lu       mumps'

  nl_abs_tol = 1E-2
  nl_max_its = 400

  l_tol = 1E-8
  l_max_its = 200
[]

[Outputs]
  [out]
    type = Exodus
  []
  perf_graph = true
[]

[Debug]
  show_var_residual_norms = true
  # show_material_props = true
[]

When mesh adaptivity is not used, the simulation runs normally. CoupledVarThresholdElementSubdomainModifier and mesh adaptivity seems to conflicts. Is there a way to solve this problem?

[GiudGiud]

Hello

This seems to come from the exodus output. It's a known issue
#28005

[PEI0214]

Hello, is there a solution to this problem now?

[jmeier]

Dear @PEI0214, I dont know of any solution but there is an corresponding issue over there with libmesh: libMesh/libmesh#3949