Merge pull request #47 from sandialabs/schwarz_inclined_contact

schwarz inclined explicit contact bugfix + tests

examples/contact/explicit-dynamic/inclined-cubes/cube-test1-1.yaml

@@ -0,0 +1,44 @@
+type: single
+input mesh file: cube-test1-1.g
+output mesh file: cube-test1-1.e
+model:
+  type: solid mechanics
+  material:
+    blocks:
+      fine: elastic
+    elastic:
+      model: linear elastic
+      elastic modulus: 1.0e+09
+      Poisson's ratio: 0.25
+      density: 1000.0
+time integrator:
+  type: central difference
+  initial time: 0.0
+  final time: 1
+  time step: 0.01
+  CFL: 0.1
+  γ: 0.5
+initial conditions:
+  displacement:
+    - node set: nsall
+      component: x
+      function: "-1.0e-04 * 1.0"
+    - node set: nsall
+      component: y
+      function: "-1.0e-04 * 0.0"
+  velocity:
+    - node set: nsall
+      component: x
+      function: "100.0 * 1.0"
+    - node set: nsall
+      component: y
+      function: "100.0 * 0.0"
+boundary conditions:
+  Schwarz contact:
+    - side set: ssx+
+      source: "cube-test1-2.yaml"
+      source block: coarse
+      source side set: ssx-
+solver:
+  type: explicit solver
+  step: explicit

examples/contact/explicit-dynamic/inclined-cubes/cube-test1-2.yaml

@@ -0,0 +1,44 @@
+type: single
+input mesh file: cube-test1-2.g
+output mesh file: cube-test1-2.e
+model:
+  type: solid mechanics
+  material:
+    blocks:
+      coarse: elastic
+    elastic:
+      model: linear elastic
+      elastic modulus: 1.0e+09
+      Poisson's ratio: 0.25
+      density: 1000.0
+time integrator:
+  type: central difference
+  initial time: 0.0
+  final time: 1
+  time step: 0.01
+  CFL: 0.1
+  γ: 0.5
+initial conditions:
+  displacement:
+    - node set: nsall
+      component: x
+      function: "1.0e-04 * 1.0"
+    - node set: nsall
+      component: y
+      function: "1.0e-04 * 0.0"
+  velocity:
+    - node set: nsall
+      component: x
+      function: "-100.0 * 1.0"
+    - node set: nsall
+      component: y
+      function: "-100.0 * 0.0"
+boundary conditions:
+  Schwarz contact:
+    - side set: ssx-
+      source: "cube-test1-1.yaml"
+      source block: fine
+      source side set: ssx+
+solver:
+  type: explicit solver
+  step: explicit

examples/contact/explicit-dynamic/inclined-cubes/cube-test2-1.yaml

@@ -0,0 +1,44 @@
+type: single
+input mesh file: cube-test2-1.g
+output mesh file: cube-test2-1.e
+model:
+  type: solid mechanics
+  material:
+    blocks:
+      fine: elastic
+    elastic:
+      model: linear elastic
+      elastic modulus: 1.0e+09
+      Poisson's ratio: 0.25
+      density: 1000.0
+time integrator:
+  type: central difference
+  initial time: 0.0
+  final time: 1
+  time step: 0.01
+  CFL: 0.1
+  γ: 0.5
+initial conditions:
+  displacement:
+    - node set: nsall
+      component: x
+      function: "-1.0e-04 * 0.9238795325112867"
+    - node set: nsall
+      component: y
+      function: "-1.0e-04 * 0.3826834323650898"
+  velocity:
+    - node set: nsall
+      component: x
+      function: "100.0 * 0.9238795325112867"
+    - node set: nsall
+      component: y
+      function: "100.0 * 0.3826834323650898"
+boundary conditions:
+  Schwarz contact:
+    - side set: ssx+
+      source: "cube-test2-2.yaml"
+      source block: coarse
+      source side set: ssx-
+solver:
+  type: explicit solver
+  step: explicit

examples/contact/explicit-dynamic/inclined-cubes/cube-test2-2.yaml

@@ -0,0 +1,44 @@
+type: single
+input mesh file: cube-test2-2.g
+output mesh file: cube-test2-2.e
+model:
+  type: solid mechanics
+  material:
+    blocks:
+      coarse: elastic
+    elastic:
+      model: linear elastic
+      elastic modulus: 1.0e+09
+      Poisson's ratio: 0.25
+      density: 1000.0
+time integrator:
+  type: central difference
+  initial time: 0.0
+  final time: 1
+  time step: 0.01
+  CFL: 0.1
+  γ: 0.5
+initial conditions:
+  displacement:
+    - node set: nsall
+      component: x
+      function: "1.0e-04 * 0.9238795325112867"
+    - node set: nsall
+      component: y
+      function: "1.0e-04 * 0.3826834323650898"
+  velocity:
+    - node set: nsall
+      component: x
+      function: "-100.0 * 0.9238795325112867"
+    - node set: nsall
+      component: y
+      function: "-100.0 * 0.3826834323650898"
+boundary conditions:
+  Schwarz contact:
+    - side set: ssx-
+      source: "cube-test2-1.yaml"
+      source block: fine
+      source side set: ssx+
+solver:
+  type: explicit solver
+  step: explicit

examples/contact/explicit-dynamic/inclined-cubes/cube-test3-1.yaml

@@ -0,0 +1,44 @@
+type: single
+input mesh file: cube-test3-1.g
+output mesh file: cube-test3-1.e
+model:
+  type: solid mechanics
+  material:
+    blocks:
+      fine: elastic
+    elastic:
+      model: linear elastic
+      elastic modulus: 1.0e+09
+      Poisson's ratio: 0.25
+      density: 1000.0
+time integrator:
+  type: central difference
+  initial time: 0.0
+  final time: 1
+  time step: 0.01
+  CFL: 0.1
+  γ: 0.5
+initial conditions:
+  displacement:
+    - node set: nsall
+      component: x
+      function: "-1.0e-04 * 0.7071067811865476"
+    - node set: nsall
+      component: y
+      function: "-1.0e-04 * 0.7071067811865475"
+  velocity:
+    - node set: nsall
+      component: x
+      function: "100.0 * 0.7071067811865476"
+    - node set: nsall
+      component: y
+      function: "100.0 * 0.7071067811865475"
+boundary conditions:
+  Schwarz contact:
+    - side set: ssx+
+      source: "cube-test3-2.yaml"
+      source block: coarse
+      source side set: ssx-
+solver:
+  type: explicit solver
+  step: explicit

examples/contact/explicit-dynamic/inclined-cubes/cube-test3-2.yaml

@@ -0,0 +1,44 @@
+type: single
+input mesh file: cube-test3-2.g
+output mesh file: cube-test3-2.e
+model:
+  type: solid mechanics
+  material:
+    blocks:
+      coarse: elastic
+    elastic:
+      model: linear elastic
+      elastic modulus: 1.0e+09
+      Poisson's ratio: 0.25
+      density: 1000.0
+time integrator:
+  type: central difference
+  initial time: 0.0
+  final time: 1
+  time step: 0.01
+  CFL: 0.1
+  γ: 0.5
+initial conditions:
+  displacement:
+    - node set: nsall
+      component: x
+      function: "1.0e-04 * 0.7071067811865476"
+    - node set: nsall
+      component: y
+      function: "1.0e-04 * 0.7071067811865475"
+  velocity:
+    - node set: nsall
+      component: x
+      function: "-100.0 * 0.7071067811865476"
+    - node set: nsall
+      component: y
+      function: "-100.0 * 0.7071067811865475"
+boundary conditions:
+  Schwarz contact:
+    - side set: ssx-
+      source: "cube-test3-1.yaml"
+      source block: fine
+      source side set: ssx+
+solver:
+  type: explicit solver
+  step: explicit

examples/contact/explicit-dynamic/inclined-cubes/cube-test4-1.yaml

@@ -0,0 +1,44 @@
+type: single
+input mesh file: cube-test4-1.g
+output mesh file: cube-test4-1.e
+model:
+  type: solid mechanics
+  material:
+    blocks:
+      fine: elastic
+    elastic:
+      model: linear elastic
+      elastic modulus: 1.0e+09
+      Poisson's ratio: 0.25
+      density: 1000.0
+time integrator:
+  type: central difference
+  initial time: 0.0
+  final time: 1
+  time step: 0.01
+  CFL: 0.1
+  γ: 0.5
+initial conditions:
+  displacement:
+    - node set: nsall
+      component: x
+      function: "-1.0e-04 * 6.123233995736766e-17"
+    - node set: nsall
+      component: y
+      function: "-1.0e-04 * 1.0"
+  velocity:
+    - node set: nsall
+      component: x
+      function: "100.0 * 6.123233995736766e-17"
+    - node set: nsall
+      component: y
+      function: "100.0 * 1.0"
+boundary conditions:
+  Schwarz contact:
+    - side set: ssx+
+      source: "cube-test4-2.yaml"
+      source block: coarse
+      source side set: ssx-
+solver:
+  type: explicit solver
+  step: explicit

examples/contact/explicit-dynamic/inclined-cubes/cube-test4-2.yaml

@@ -0,0 +1,44 @@
+type: single
+input mesh file: cube-test4-2.g
+output mesh file: cube-test4-2.e
+model:
+  type: solid mechanics
+  material:
+    blocks:
+      coarse: elastic
+    elastic:
+      model: linear elastic
+      elastic modulus: 1.0e+09
+      Poisson's ratio: 0.25
+      density: 1000.0
+time integrator:
+  type: central difference
+  initial time: 0.0
+  final time: 1
+  time step: 0.01
+  CFL: 0.1
+  γ: 0.5
+initial conditions:
+  displacement:
+    - node set: nsall
+      component: x
+      function: "1.0e-04 * 6.123233995736766e-17"
+    - node set: nsall
+      component: y
+      function: "1.0e-04 * 1.0"
+  velocity:
+    - node set: nsall
+      component: x
+      function: "-100.0 * 6.123233995736766e-17"
+    - node set: nsall
+      component: y
+      function: "-100.0 * 1.0"
+boundary conditions:
+  Schwarz contact:
+    - side set: ssx-
+      source: "cube-test4-1.yaml"
+      source block: fine
+      source side set: ssx+
+solver:
+  type: explicit solver
+  step: explicit

examples/contact/explicit-dynamic/inclined-cubes/cubes-test1.yaml

@@ -0,0 +1,14 @@
+type: multi
+domains: ["cube-test1-1.yaml", "cube-test1-2.yaml"]
+Exodus output interval: 1
+CSV output interval: 0
+initial time: -1.0e-06
+time step: 1.e-06
+final time: .001
+same time step for domains: true
+relaxation parameter: 1.0
+naive stabilized: true
+minimum iterations: 1
+maximum iterations: 16
+relative tolerance: 1.0e-12
+absolute tolerance: 1.0e-08

examples/contact/explicit-dynamic/inclined-cubes/cubes-test2.yaml

@@ -0,0 +1,14 @@
+type: multi
+domains: ["cube-test2-1.yaml", "cube-test2-2.yaml"]
+Exodus output interval: 1
+CSV output interval: 0
+initial time: -1.0e-06
+time step: 1.e-06
+final time: .001
+same time step for domains: true
+relaxation parameter: 1.0
+naive stabilized: true
+minimum iterations: 1
+maximum iterations: 16
+relative tolerance: 1.0e-12
+absolute tolerance: 1.0e-08