An error occurs when using both MeshModifiers and Contact simultaneously

[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
The following error occurred when using both meshmodifiers and Contact simultaneously. Is there any good solution to this problem?

*** ERROR ***
/home/pdd/projects8/moose/modules/combined/problems/0410/test1.i:108.3:
The following occurred in the Action 'c1' of type ContactAction.

Contact requires updated coordinates.  Use the 'displacements = ...' parameter in the Mesh block.

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
  boundary_restricted_elem_integrity_check = false
[]

new_domain_blocks = '2001 2002'

[Mesh]
  [mesh]
    type = GeneratedMeshGenerator
    dim = 3
    xmin = 0
    xmax = 10
    ymin = 0
    ymax = 10
    zmin = 0
    zmax = 10
    nx = 10
    ny = 10
    nz = 10
  []
  [box1]
    input = mesh
    type = SubdomainBoundingBoxGenerator
    block_id = 1
    bottom_left = '0 0 5'
    top_right = '10 10 10'
  []
  [breakmesh]
    input = box1
    type = BreakMeshByBlockGenerator
    block_pairs = '0 1'
  []

  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}'
  []
[]

# ===== Tensor mechanics for all active domain blocks =====
[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [all]
        strain = SMALL
        incremental = true
        block = '${new_domain_blocks}'
        add_variables = false
        generate_output = 'stress_yy stress_zz min_principal_stress max_principal_stress'
      []
    []
  []
[]

# ===== Gravity =====
[Kernels]
  [Gravity1]
    type = Gravity
    block = '${new_domain_blocks}'
    variable = disp_y
    value = -9.8
  []
[]

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

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

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

[Contact]
  [c1]
    primary = '6'
    secondary = '7'
    model = coulomb
    formulation = mortar
    c_normal = 5e8
    c_tangential = 5e7
    friction_coefficient = 0.5
    # tangential_lm_scaling = 1.0e-6
    # normal_lm_scaling = 1.0e-6
  []
[]

# ===== Materials (linear-elastic to keep it simple) =====
[Materials]
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    block = '2001 2002'
    youngs_modulus = 30e9
    poissons_ratio = 0.167
  []
  [stress]
    type = ComputeFiniteStrainElasticStress
    block = '2001 2002'
  []

  [density]
    type = GenericConstantMaterial
    block = '2001 2002'
    prop_names = density
    prop_values = 2400.0
  []
[]

[MeshModifiers]
  [m1]
    type = TimedSubdomainModifier
    times = '1e7 2e7'
    blocks_from = '0 1'
    blocks_to = '2001 2002'
    execute_on = 'TIMESTEP_BEGIN'
  []
[]

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

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

  end_time = 3e7
  dt = 1e7

  solve_type = 'PJFNK'
  petsc_options = '-snes_converged_reason'

  petsc_options_iname = '-ksp_gmres_restart -pc_type -pc_hypre_type -pc_hypre_boomeramg_max_iter'
  petsc_options_value = '  201               hypre    boomeramg      10'

  nl_abs_tol = 1E-5
  nl_max_its = 400

  l_tol = 1E-8
  l_max_its = 200
[]

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

When I set 'use_displaced_mesh = true', the following problem occurred.

Time Step 0, time = 0

Time Step 1, time = 1e+07, dt = 1e+07

Performing automatic scaling calculation

Segmentation fault (core dumped)

[GiudGiud]

Hello

When I set 'use_displaced_mesh = true', the following problem occurred.

you'll want to set this on a per object basis. The mesh modifiers have likely not been coded for the displaced mesh

Thanks for reporting. i ll create an issue

[GiudGiud]

For the seg fault I am hitting this in devel mode

A libMesh DegenerateMap was raised during FEProblemBase::computeUserObjectsInternal
Assertion `comp < this->n_comp(s,var)' failed.
comp = 0
this->n_comp(s,var) = 0


Stack frames: 26
0: 0   libmesh_devel.0.dylib               0x000000010980b8d8 libMesh::print_trace(std::__1::basic_ostream<char, std::__1::char_traits<char>>&) + 1008
1: 1   libmesh_devel.0.dylib               0x00000001098085b0 libMesh::MacroFunctions::report_error(char const*, int, char const*, char const*, std::__1::basic_ostream<char, std::__1::char_traits<char>>&) + 280
2: 2   libmoose-devel.0.dylib              0x00000001040779e0 libMesh::DofObject::dof_number(unsigned int, unsigned int, unsigned int) const + 856
3: 3   libmoose-devel.0.dylib              0x0000000105540358 Assembly::computeFaceMap(libMesh::Elem const&, unsigned int, std::__1::vector<double, std::__1::allocator<double>> const&) + 5728
4: 4   libmoose-devel.0.dylib              0x000000010553db98 Assembly::computeADFace(libMesh::Elem const&, unsigned int) + 248
5: 5   libmoose-devel.0.dylib              0x000000010554af70 Assembly::reinitElemFaceRef(libMesh::Elem const*, unsigned int, double, std::__1::vector<libMesh::Point, std::__1::allocator<libMesh::Point>> const*, std::__1::vector<double, std::__

you'll want to check the block restriction of variables and objects acting on them. Something is not restricted correctly

[GiudGiud]

The displacements are block restricted.
I dont see subdomain 2002 in the domain.

Could you keep the displacements block restricted to all the 3D domains?
And you'll want to check other contact inputs to see if they should be defined on lowerD (2D) domains as well. I suspect we only want them on 3D domains

[PEI0214]

Hi GiudGiud,
I appreciate you pointing out the subdomain issue! To explain the setup:

• Subdomains 2001 and 2002 are newly created, initially empty blocks, which will be activated sequentially to simulate the step-by-step filling process of the embankment.
• The block restriction on displacement variables is intentional: it ensures the solver only computes displacements on the already filled, active subdomains, avoiding calculations on empty/undeveloped regions which would cause numerical issues.

[GiudGiud]

The block restriction on displacement variables is intentional: it ensures the solver only computes displacements on the already filled, active subdomains, avoiding calculations on empty/undeveloped regions which would cause numerical issues.

you can use a NullKernel on the empty regions to avoid numerical issues.
It will have the added advantage to get around the issue you are currently hitting.
On the lowerD regions i think it is best not to define displacements as well.