Add ahead tests

examples/ahead/nonoverlap/bracket/bracket-tet-tet2.jou

@@ -0,0 +1,81 @@
+reset
+#left part of Omega 
+brick x 0.2 y 0.2 z 0.01
+#movie volume 1 -2 -2 0
+create Cylinder height 1 radius 0.075
+#move Volume 2 midpoint location curve 9  include_merged 
+subtract volume 2 from volume 1 
+brick x 0.01 y 0.2 z 0.2
+#move Volume 3 midpoint location curve 8  include_merged 
+volume 3 move -0.1 0 0
+unite volume 1 volume 3
+
+webcut volume 1 with plane xplane offset 0.0
+
+
+delete volume 1 
+
+block 1 volume all
+block 1 name 'left'
+
+sideset 1 surface 32 33
+sideset 1 name "left-frontSS"
+
+nodeset 1 surface 16
+nodeset 1 name "left-backNS"
+nodeset 2 surface 32 33
+nodeset 2 name "left-frontNS"
+nodeset 3 volume all
+nodeset 3 name "left-nsall"
+
+volume all scheme Tetmesh
+volume 4 size auto factor 7
+#volume all size auto factor 10
+block 1 element type TETRA4
+mesh volume all
+#refine volume all numsplit 1 
+export mesh "bracket-1.g" overwrite 
+
+reset
+#right part of Omega 
+brick x 0.2 y 0.2 z 0.01
+#movie volume 1 -2 -2 0
+create Cylinder height 1 radius 0.075
+#move Volume 2 midpoint location curve 9  include_merged 
+subtract volume 2 from volume 1 
+brick x 0.01 y 0.2 z 0.2
+#move Volume 3 midpoint location curve 8  include_merged 
+volume 3 move -0.1 0 0
+unite volume 1 volume 3
+
+webcut volume 1 with plane xplane offset 0.0
+delete volume 4 
+
+block 1 volume all
+block 1 name 'right'
+
+sideset 1 surface 26 25
+sideset 1 name "right-backSS"
+
+nodeset 1 surface 26 25
+nodeset 1 name "right-backNS"
+nodeset 2 surface 6 
+nodeset 2 name "right-frontNS"
+nodeset 3 volume all
+nodeset 3 name "right-nsall"
+
+volume all scheme Tetmesh
+volume all size auto factor 7
+#volume all size auto factor 14
+block 1 element type TETRA4
+mesh volume all 
+refine volume all numsplit 1 
+export mesh "bracket-2.g" overwrite 
+
+
+
+
+
+
+
+

examples/ahead/nonoverlap/bracket/dynamic/bracket-1.yaml

@@ -0,0 +1,40 @@
+type: single
+input mesh file: ../bracket-1.g
+output mesh file: bracket-1.e
+model:
+  type: solid mechanics
+  material:
+    blocks:
+      left: hyperelastic
+    hyperelastic:
+      model: neohookean
+      elastic modulus: 200.0e+09
+      Poisson's ratio: 0.25
+      density: 7800.0
+time integrator:
+  type: Newmark
+  β: 0.49
+  γ: 0.9
+boundary conditions:
+  Dirichlet:
+    - node set: left-backNS
+      component: x
+      function: "0.0"
+    - node set: left-backNS
+      component: y
+      function: "0.0"
+    - node set: left-backNS
+      component: z
+      function: "0.0"
+  Schwarz nonoverlap:
+    - side set: left-frontSS
+      source: bracket-2
+      source block: right
+      source side set: right-backSS
+solver:
+  type: Hessian minimizer
+  step: full Newton
+  minimum iterations: 1
+  maximum iterations: 16
+  relative tolerance: 1.0e-10
+  absolute tolerance: 1.0e-06

examples/ahead/nonoverlap/bracket/dynamic/bracket-2.yaml

@@ -0,0 +1,41 @@
+type: single
+input mesh file: ../bracket-2.g
+output mesh file: bracket-2.e
+model:
+  type: solid mechanics
+  material:
+    blocks:
+      right: hyperelastic
+    hyperelastic:
+      model: neohookean
+      elastic modulus: 200.0e+09
+      Poisson's ratio: 0.25
+      density: 7800.0
+time integrator:
+  type: Newmark
+  β: 0.49
+  γ: 0.9
+initial conditions:
+  velocity:
+    - node set: right-nsall
+      component: x
+      function: "0.0"
+    - node set: right-nsall
+      component: y
+      function: "0.0"
+    - node set: right-nsall
+      component: z
+      function: "ifelse(x > 0.099, 100.0, 0.0)"
+boundary conditions:
+  Schwarz nonoverlap:
+    - side set: right-backSS
+      source: bracket-1
+      source block: left
+      source side set: left-frontSS
+solver:
+  type: Hessian minimizer
+  step: full Newton
+  minimum iterations: 1
+  maximum iterations: 16
+  relative tolerance: 1.0e-10
+  absolute tolerance: 1.0e-06

examples/ahead/nonoverlap/bracket/dynamic/bracket.yaml

@@ -0,0 +1,13 @@
+type: multi
+domains: ["bracket-1.yaml", "bracket-2.yaml"]
+Exodus output interval: 1
+CSV output interval: 0
+initial time: 0.0
+final time: 2.5e-3
+#final time: 20.0e-3 #original final time
+time step: 1.0e-6
+minimum iterations: 1
+maximum iterations: 128
+relative tolerance: 1.0e-06
+absolute tolerance: 1.0e-06
+relaxation parameter: 0.2

examples/ahead/nonoverlap/plate/dynamic/plate-1.yaml

@@ -13,8 +13,8 @@ model:
       density: 7800.0
 time integrator:
   type: Newmark
-  β: 0.25
-  γ: 0.5
+  β: 0.49
+  γ: 0.9
 boundary conditions:
   Dirichlet:
     - node set: nsx-

examples/ahead/nonoverlap/plate/dynamic/plate-2.yaml

@@ -13,8 +13,8 @@ model:
       density: 7800.0
 time integrator:
   type: Newmark
-  β: 0.25
-  γ: 0.5
+  β: 0.49
+  γ: 0.9
 initial conditions:
   velocity:
     - node set: nsall

examples/ahead/nonoverlap/plate/dynamic/plate.yaml

@@ -9,4 +9,4 @@ minimum iterations: 1
 maximum iterations: 128
 relative tolerance: 1.0e-12
 absolute tolerance: 1.0e-08
-relaxation parameter: 0.4 
+relaxation parameter: 0.1 

examples/ahead/nonoverlap/plate/plate.jou

@@ -1,7 +1,7 @@
 ${side = 0.2}
 ${height = 0.03}
 ${overlap_vol_frac = 0.0}
-${h_fine = 0.01}
+${h_fine = 0.05}
 ${h_coarse = 0.01}
 
 ${area = side * height}

examples/ahead/overlap/bracket/dynamic/bracket-1.yaml

@@ -13,8 +13,8 @@ model:
       density: 7800.0
 time integrator:
   type: Newmark
-  β: 0.25
-  γ: 0.5
+  β: 0.49
+  γ: 0.9
 boundary conditions:
   Dirichlet:
     - node set: left-backNS

examples/ahead/overlap/bracket/dynamic/bracket-2.yaml

@@ -13,8 +13,8 @@ model:
       density: 7800.0
 time integrator:
   type: Newmark
-  β: 0.25
-  γ: 0.5
+  β: 0.49
+  γ: 0.9
 initial conditions:
   velocity:
     - node set: right-nsall

examples/ahead/overlap/bracket/dynamic/bracket.yaml

@@ -7,6 +7,6 @@ final time: 2.5e-3
 #final time: 20.0e-3 #original final time
 time step: 1.0e-6
 minimum iterations: 1
-maximum iterations: 16
-relative tolerance: 1.0e-12
-absolute tolerance: 1.0e-08
+maximum iterations: 128
+relative tolerance: 1.0e-10
+absolute tolerance: 1.0e-06

examples/ahead/overlap/laser-weld/dynamic-opinf-fom/clamped/holder-0.yaml

@@ -24,7 +24,7 @@ boundary conditions:
       function: "0.0"
     - node set: nodeset_negative_y
       component: y
-      function: "-0.1 * t"
+      function: "-0.1 * (0.5 - 0.5 * cos(pi * t))"
     - node set: nodeset_negative_y
       component: z
       function: "0.0"

examples/ahead/overlap/laser-weld/dynamic-opinf-fom/clamped/holder-1.yaml

@@ -24,7 +24,7 @@ boundary conditions:
       function: "0.0"
     - node set: nodeset_positive_y
       component: y
-      function: "0.1 * t"
+      function: "0.1 * (0.5 - 0.5 * cos(pi * t))"
     - node set: nodeset_positive_y
       component: z
       function: "0.0"

examples/ahead/overlap/laser-weld/dynamic-opinf-rom/clamped/gauge.yaml

@@ -18,6 +18,13 @@ time integrator:
   β: 0.25
   γ: 0.5
 boundary conditions:
+  Dirichlet:
+    - node set: nodeset_negative_y
+      component: x
+      function: "0.0"
+    - node set: nodeset_negative_y
+      component: z
+      function: "0.0"
   Schwarz overlap:
     - side set: sideset_negative_y
       source: holder-0

examples/ahead/overlap/laser-weld/dynamic-opinf-rom/clamped/holder-0.yaml

@@ -17,15 +17,15 @@ time integrator:
   γ: 0.5
 boundary conditions:
   Dirichlet:
-    - node set: nodeset_negative_y
-      component: x
-      function: "0.0"
+    #- node set: nodeset_negative_y
+    #  component: x
+    #  function: "0.0"
     - node set: nodeset_negative_y
       component: y
-      function: "-0.1 * t"
-    - node set: nodeset_negative_y
-      component: z
-      function: "0.0"
+      function: "-0.1 * (0.5 - 0.5 * cos(pi * t))"
+        #- node set: nodeset_negative_y
+        #component: z
+        #function: "0.0"
   Schwarz overlap:
     - side set: sideset_positive_y
       source: gauge

examples/ahead/overlap/laser-weld/dynamic-opinf-rom/clamped/holder-1.yaml

@@ -18,15 +18,15 @@ time integrator:
   γ: 0.5
 boundary conditions:
   Dirichlet:
-    - node set: nodeset_positive_y
-      component: x
-      function: "0.0"
+    #- node set: nodeset_positive_y
+    #  component: x
+    #  function: "0.0"
     - node set: nodeset_positive_y
       component: y
-      function: "0.1 * t"
-    - node set: nodeset_positive_y
-      component: z
-      function: "0.0"
+      function: "0.1 * (0.5 - 0.5 * cos(pi * t))"
+    #- node set: nodeset_positive_y
+    #component: z
+    #function: "0.0"
   Schwarz overlap:
     - side set: sideset_negative_y
       source: gauge

examples/ahead/overlap/laser-weld/dynamic-opinf-rom/clamped/laser-weld.yaml

@@ -3,8 +3,8 @@ domains: ["holder-0.yaml", "gauge.yaml", "holder-1.yaml"]
 Exodus output interval: 1
 initial time: 0.0
 final time: 0.1
-time step: 0.001
+time step: 0.0001
 minimum iterations: 1
 maximum iterations: 64
-relative tolerance: 1.0e-10
-absolute tolerance: 1.0e-06
+relative tolerance: 1.0e-07
+absolute tolerance: 1.0e-05

examples/ahead/overlap/notched-cylinder/dynamic-opinf-fom/notched-cylinder-1.yaml

@@ -1,5 +1,5 @@
 type: single
-input mesh file: notched-cylinder-1.g
+input mesh file: ../notched-cylinder-1.g
 output mesh file: notched-cylinder-1.e
 CSV output interval: 1
 CSV write sidesets: true

examples/ahead/overlap/notched-cylinder/dynamic-opinf-fom/notched-cylinder-2.yaml

@@ -1,5 +1,5 @@
 type: single
-input mesh file: notched-cylinder-2.g
+input mesh file: ../notched-cylinder-2.g
 output mesh file: notched-cylinder-2.e
 CSV output interval: 1
 CSV write sidesets: true
@@ -37,6 +37,7 @@ boundary conditions:
       source: notched-cylinder-1
       source side set: +Z_bottomSS
       source block: fine
+      search tolerance: 5.0e-02
 solver:
   type: Hessian minimizer
   step: full Newton