Contact Software Design Description

Introduction

The SDD for Contact is a description of the software product and serves as guideline that describes the architecture of the system and all of its parts.

Dependencies

The Contact application is developed using MOOSE and is based on various modules, as such the SDD for Contact is dependent upon the following documents.

• Framework System Design Description

Requirements Cross Reference

The following is a list of all design documents and the associated requirements for Contact.

• contact: FEProblemBase
• F4.1.1The system shall use grid sequencing in order to improve the performance of the nonlinear solve in a frictional contact problem

  Specification: bouncing-block-contact/grid-sequencing:grid_sequencing

  Design: FEProblemBase

  Issue(s): #14166

  20
• F4.1.2The system shall be able to produce the same result when not using grid sequencing

  Specification: bouncing-block-contact/grid-sequencing:no_grid_sequencing

  Design: FEProblemBase

  Issue(s): #14166

  20

• contact: Constraints System
• F4.1.3We will solve the frictionless bouncing block problem with nodal constraint enforcement, mortar application of forces, and min NCP function

  Specification: bouncing-block-contact:frictionless-nodal-min-lm-mortar-disp-coarse

  Design: Constraints System

  Issue(s): #13080

  20
• F4.1.4We will solve the frictionless bouncing block problem with nodal constraint enforcement, mortar application of forces, and fb NCP function

  Specification: bouncing-block-contact:frictionless-nodal-fb-lm-mortar-disp-coarse

  Design: Constraints System

  Issue(s): #13080

  20
• F4.1.5We will solve the frictionless bouncing block problem with nodal constraint enforcement, nodal application of forces, and min NCP function

  Specification: bouncing-block-contact:frictionless-nodal-min-lm-nodal-disp-coarse

  Design: Constraints System

  Issue(s): #13080

  20
• F4.1.6We will solve the frictionless bouncing block problem with nodal constraint enforcement, nodal application of forces, and fb NCP function

  Specification: bouncing-block-contact:frictionless-nodal-fb-lm-nodal-disp-coarse

  Design: Constraints System

  Issue(s): #13080

  20
• F4.1.7We will solve the frictionless bouncing block problem with mortar constraint enforcement, mortar application of forces, and min NCP function

  Specification: bouncing-block-contact:frictionless-mortar-min-lm-mortar-disp-coarse

  Design: Constraints System

  Issue(s): #13080

  20
• F4.1.8We will solve the frictionless bouncing block problem with mortar constraint enforcement, mortar application of forces, and fb NCP function

  Specification: bouncing-block-contact:frictionless-mortar-fb-lm-mortar-disp-coarse

  Design: Constraints System

  Issue(s): #13080

  20
• F4.1.9We will solve the frictional bouncing block problem with mortar constraint enforcement, mortar application of forces, and min NCP function

  Specification: bouncing-block-contact:frictional-mortar-min-lm-mortar-disp-coarse

  Design: Constraints System

  Issue(s): #13080

  20
• F4.1.10We will solve the frictional bouncing block problem with mortar constraint enforcement, mortar application of forces, and fb NCP function

  Specification: bouncing-block-contact:frictional-mortar-fb-lm-mortar-disp-coarse

  Design: Constraints System

  Issue(s): #13080

  20
• F4.1.11We will solve the frictional bouncing block problem with nodal constraint enforcement, mortar application of forces, and min NCP function

  Specification: bouncing-block-contact:frictional-nodal-min-lm-mortar-disp-coarse

  Design: Constraints System

  Issue(s): #13080

  14
• F4.1.12We will solve the frictional bouncing block problem with nodal constraint enforcment for the normal LM using min NCP, mortar constraint enforcement of the tangential LM with fb, and mortar application of forces

  Specification: bouncing-block-contact:frictional-nodal-min-normal-lm-mortar-fb-tangential-lm-mortar-disp-coarse

  Design: Constraints System

  Issue(s): #13080

  20
• F4.1.13We will solve the frictional bouncing block problem with nodal constraint enforcment for the normal LM using min NCP, mortar constraint enforcement of the tangential LM with fb, and mortar application of forces using the ContactAction. The result shall be identical to that obtained by directly specifying individual objects.

  Specification: bouncing-block-contact:frictional-nodal-min-normal-lm-mortar-fb-tangential-lm-mortar-action

  Design: Constraints System

  Issue(s): #13080

  20
• F4.1.14The system shall be able to solve frictionless mechanical contact using a reduced active nonlinear function set scheme (RANFS) in conjunction with a node-face geometric discretization. The RANFS scheme shall be
  1. nonsingular both with bounds projection and
  2. without bounds projection and be
  3. solvable with amg both with bounds projection
  4. and without bounds projection.
  5. The system's RANFS scheme shall have a perfect Jacobian for mechanical contact that only has one non-zero normal component
  6. The system shall be able to detect when a slave node is ping-ponging back and forth between different master faces and consequently tie the locations of the slave and corresponding master node using Lagrange Multipliers corresponding to equality constraints, e.g. more RANFS
  7. The system shall be able to solve a smaller model of the full ping-ponging problem

  Specification: bouncing-block-contact:reduced_active_set

  Design: Constraints System

  Issue(s): #13080

  2020202020220
• F4.1.15Using a RANFS scheme with Lagrange multipliers corresponding to equality constraints the system shall be able to
  1. tie nodes together and
  2. have a perfect Jacobian

  Specification: bouncing-block-contact:tied_node

  Design: Constraints System

  Issue(s): #13080

  2020
• F4.2.1We shall be able to do mechanical contact with the new mortar system. It's beautiful.

  Specification: mechanical-small-problem:mortar-mech

  Design: Constraints System

  Issue(s): #13080#13676

  20
• F4.2.2We shall be able to do mechanical contact with mortar for the displacements and node enforcement of the constraints.

  Specification: mechanical-small-problem:frictionless-nodal-lm-mortar-disp

  Design: Constraints System

  Issue(s): #13080#13676

  20
• F4.2.3We shall be able to do mechanical contact with mortar for the displacements and node enforcement of the constraints while using the Contact action block. The results should be identical to those gained by specifying individual mortar constraints, variables, and mesh objects.

  Specification: mechanical-small-problem:frictionless-nodal-lm-mortar-disp-action

  Design: Constraints System

  Issue(s): #13080#13676

  20

[grid_sequencing]
  type = Exodiff
  input = grid-sequencing.i
  exodiff = grid-sequencing_exo.e
  requirement = 'The system shall use grid sequencing in order to improve the performance of the nonlinear solve in a frictional contact problem'
  cli_args = 'Postprocessors/active=contact'
  # Fix below when I understand duplicate element centroids
  map = False
  partial = True
  mesh_mode = replicated # Fix when mortar segment mesh is distributed
[]

[no_grid_sequencing]
  type = Exodiff
  input = grid-sequencing.i
  exodiff = no-grid-sequencing_exo.e
  requirement = 'The system shall be able to produce the same result when not using grid sequencing'
  cli_args = 'Postprocessors/active=contact Mesh/uniform_refine=4 Executioner/num_grids=1 Outputs/file_base=no-grid-sequencing_exo'
  # Fix below when I understand duplicate element centroids
  map = False
  partial = True
  mesh_mode = replicated # Fix when mortar segment mesh is distributed
  rel_err = 2e-5
[]

[frictionless-nodal-min-lm-mortar-disp-coarse]
  type = 'Exodiff'
  input = 'frictionless-nodal-min-lm-mortar-disp.i'
  exodiff = 'frictionless-nodal-min-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictionless bouncing block problem with nodal constraint enforcement, mortar application of forces, and min NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
  allow_test_objects = true
[]

[frictionless-nodal-fb-lm-mortar-disp-coarse]
  type = 'Exodiff'
  input = 'frictionless-nodal-fb-lm-mortar-disp.i'
  exodiff = 'frictionless-nodal-fb-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictionless bouncing block problem with nodal constraint enforcement, mortar application of forces, and fb NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
  allow_test_objects = true
[]

[frictionless-nodal-min-lm-nodal-disp-coarse]
  type = 'Exodiff'
  input = 'frictionless-nodal-min-lm-nodal-disp.i'
  exodiff = 'frictionless-nodal-min-lm-nodal-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictionless bouncing block problem with nodal constraint enforcement, nodal application of forces, and min NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
  allow_test_objects = true
[]

[frictionless-nodal-fb-lm-nodal-disp-coarse]
  type = 'Exodiff'
  input = 'frictionless-nodal-fb-lm-nodal-disp.i'
  exodiff = 'frictionless-nodal-fb-lm-nodal-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictionless bouncing block problem with nodal constraint enforcement, nodal application of forces, and fb NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
  allow_test_objects = true
[]

[frictionless-mortar-min-lm-mortar-disp-coarse]
  allow_test_objects = true
  type = 'Exodiff'
  input = 'frictionless-mortar-min-lm-mortar-disp.i'
  exodiff = 'frictionless-mortar-min-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictionless bouncing block problem with mortar constraint enforcement, mortar application of forces, and min NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
[]

[frictionless-mortar-fb-lm-mortar-disp-coarse]
  type = 'Exodiff'
  input = 'frictionless-mortar-fb-lm-mortar-disp.i'
  exodiff = 'frictionless-mortar-fb-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictionless bouncing block problem with mortar constraint enforcement, mortar application of forces, and fb NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
  allow_test_objects = true
[]

[frictional-mortar-min-lm-mortar-disp-coarse]
  type = 'Exodiff'
  input = 'frictional-mortar-min-lm-mortar-disp.i'
  exodiff = 'frictional-mortar-min-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictional bouncing block problem with mortar constraint enforcement, mortar application of forces, and min NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
  max_parallel = 1
  allow_test_objects = true
[]

[frictional-mortar-fb-lm-mortar-disp-coarse]
  type = 'Exodiff'
  input = 'frictional-mortar-fb-lm-mortar-disp.i'
  exodiff = 'frictional-mortar-fb-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictional bouncing block problem with mortar constraint enforcement, mortar application of forces, and fb NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active=''"
  allow_test_objects = true
[]

[frictional-nodal-min-lm-mortar-disp-coarse]
  type = 'Exodiff'
  input = 'frictional-nodal-min-lm-mortar-disp.i'
  exodiff = 'frictional-nodal-min-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictional bouncing block problem with nodal constraint enforcement, mortar application of forces, and min NCP function'
  cli_args = "Mesh/file=long-bottom-block-1elem-blocks-coarse.e Executioner/end_time=50 Postprocessors/active='' Outputs/exodus/sync_only=true --redirect-output output.log"
  allow_warnings = True
  max_parallel = 1
  method = 'OPT'
  allow_test_objects = true
[]

[frictional-nodal-min-normal-lm-mortar-fb-tangential-lm-mortar-disp-coarse]
  type = 'Exodiff'
  input = 'frictional-nodal-min-normal-lm-mortar-fb-tangential-lm-mortar-disp.i'
  exodiff = 'frictional-nodal-min-normal-lm-mortar-fb-tangential-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictional bouncing block problem with nodal constraint enforcment for the normal LM using min NCP, mortar constraint enforcement of the tangential LM with fb, and mortar application of forces'
  cli_args = "Executioner/end_time=50 Postprocessors/active=''"
  allow_test_objects = true
[]

[frictional-nodal-min-normal-lm-mortar-fb-tangential-lm-mortar-action]
  type = 'Exodiff'
  input = 'frictional-nodal-min-normal-lm-mortar-fb-tangential-lm-mortar-action.i'
  exodiff = 'frictional-nodal-min-normal-lm-mortar-fb-tangential-lm-mortar-action_out.e'
  mesh_mode = 'REPLICATED'
  requirement = 'We will solve the frictional bouncing block problem with nodal constraint enforcment for the normal LM using min NCP, mortar constraint enforcement of the tangential LM with fb, and mortar application of forces using the ContactAction. The result shall be identical to that obtained by directly specifying individual objects.'
  cli_args = "Executioner/end_time=50 Postprocessors/active=''"
[]

[reduced_active_set]
  requirement = 'The system shall be able to solve frictionless mechanical contact using a reduced active nonlinear function set scheme (RANFS) in conjunction with a node-face geometric discretization. The RANFS scheme shall be'
  [ranfs-project-svd]
    type = 'Exodiff'
    input = 'frictionless-nodal-reduced-active-set.i'
    exodiff = 'ranfs-project-svd.e'
    absent_out = '[1-9]+[0-9]* of 16 singular values'
    expect_out = '0 of 16 singular values'
    detail = 'nonsingular both with bounds projection and'
    cli_args = 'Outputs/file_base=ranfs-project-svd -pc_type svd -pc_svd_monitor'
  []
  [ranfs-dont-project-svd]
    type = 'Exodiff'
    input = 'frictionless-nodal-reduced-active-set.i'
    exodiff = 'ranfs-dont-project-svd.e'
    detail = 'without bounds projection and be'
    cli_args = 'Outputs/file_base=ranfs-dont-project-svd Executioner/line_search=none -pc_type svd -pc_svd_monitor'
  []
  [ranfs-project-amg]
    type = 'Exodiff'
    input = 'frictionless-nodal-reduced-active-set.i'
    exodiff = 'ranfs-project-amg.e'
    detail = 'solvable with amg both with bounds projection'
    cli_args = 'Outputs/file_base=ranfs-project-amg'
    expect_out = '[bB][oO][oO][mM][eE][rR][aA][mM][gG]'
  []
  [ranfs-dont-project-amg]
    type = 'Exodiff'
    input = 'frictionless-nodal-reduced-active-set.i'
    exodiff = 'ranfs-dont-project-amg.e'
    detail = 'and without bounds projection.'
    cli_args = 'Outputs/file_base=ranfs-dont-project-amg Executioner/line_search=none'
    expect_out = '[bB][oO][oO][mM][eE][rR][aA][mM][gG]'
  []
  [ranfs-jacobian]
    type = 'PetscJacobianTester'
    input = 'frictionless-nodal-reduced-active-set.i'
    difference_tol = 1e-7
    cli_args = "Outputs/active='' Executioner/line_search=none"
    expect_out = '[bB][oO][oO][mM][eE][rR][aA][mM][gG]'
    detail = "The system's RANFS scheme shall have a perfect Jacobian for mechanical contact that only has one non-zero normal component"
  []
  [full-bouncing-block]
    type = 'Exodiff'
    input = 'bouncing-block-ranfs.i'
    exodiff = 'bouncing-block-ranfs_exo.e'
    detail = 'The system shall be able to detect when a slave node is ping-ponging back and forth between different master faces and consequently tie the locations of the slave and corresponding master node using Lagrange Multipliers corresponding to equality constraints, e.g. more RANFS'
    cli_args = "Postprocessors/active=''"
    custom_cmp = 'bouncing-block-ranfs.cmp'
    heavy = true
  []
  [bouncing-block]
    type = 'Exodiff'
    input = 'bouncing-block-ranfs.i'
    exodiff = 'bouncing-block-coarse.e'
    detail = 'The system shall be able to solve a smaller model of the full ping-ponging problem'
    cli_args = "Postprocessors/active='' Mesh/file=long-bottom-block-no-lower-d-coarse.e Executioner/end_time=100 Outputs/file_base=bouncing-block-coarse"
  []
[]

[tied_node]
  requirement = 'Using a RANFS scheme with Lagrange multipliers corresponding to equality constraints the system shall be able to '
  [physics]
    type = 'Exodiff'
    input = tied-nodes.i
    exodiff = 'tied-nodes_exo.e'
    absent_out = '[1-9]+[0-9]* of 16 singular values'
    expect_out = '0 of 16 singular values'
    detail = 'tie nodes together and'
    allow_test_objects = True
    cli_args = '-pc_type svd -pc_svd_monitor'
  []
  [jac]
    type = 'PetscJacobianTester'
    input = tied-nodes.i
    cli_args = "Outputs/active=''"
    detail = 'have a perfect Jacobian'
    allow_test_objects = True
  []
[]

[mortar-mech]
  type = 'Exodiff'
  input = 'mortar-mech.i'
  exodiff = 'mortar-mech_out.e'
  mesh_mode = 'REPLICATED'
  requirement = "We shall be able to do mechanical contact with the new mortar system. It's beautiful."
  # Need to set an absolute zero here because there is a slight difference in results between
  # running with FD and PJFNK at the time step when we have 'resting contact' e.g. there is no
  # gap, but also no force
  abs_zero = '2e-8'
  allow_test_objects = true
[]

[frictionless-nodal-lm-mortar-disp]
  type = 'Exodiff'
  input = 'frictionless-nodal-lm-mortar-disp.i'
  exodiff = 'frictionless-nodal-lm-mortar-disp_out.e'
  mesh_mode = 'REPLICATED'
  requirement = "We shall be able to do mechanical contact with mortar for the displacements and node enforcement of the constraints."
  allow_test_objects = true
[]

[frictionless-nodal-lm-mortar-disp-action]
  type = 'Exodiff'
  input = 'frictionless-nodal-lm-mortar-disp-action.i'
  exodiff = 'frictionless-nodal-lm-mortar-disp-action_out.e'
  mesh_mode = 'REPLICATED'
  requirement = "We shall be able to do mechanical contact with mortar for the displacements and node enforcement of the constraints while using the Contact action block. The results should be identical to those gained by specifying individual mortar constraints, variables, and mesh objects."
[]