AHeaD Test Cases in Norma.jl

Irina Tezaur

​

Sandia National Laboratories, Livermore, CA, USA

​

SAND2025-14004PE

November 6, 2025

PRESENTED BY

Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

2

General Notes

Clamped Problem

Figure above: clamped problem solution with symmetric Gaussian IC.

Figure below: clamped problem solution with rounded square IC.

Cuboid Problem

Notched Cylinder Problem

• Geometry: See right. Larger radius = 0.032m. Height=0.048m
• Initial condition: 0 disp, 0 velo
• Boundary conditions:
• 0 x-disp in positive x (symmetry BC)
• 0 y-disp in positive y (symmetry BC)
• 0 z-disp at bottom (symmetry BC)
• z-disp = 0.0064 * (0.5 - 0.5 * cos(pi * t)) at top (pull from top)

​

• Notes:
• Can run linear or nonlinear.
• More contrived version of tension specimen (next slide)
• Had extreme difficulty converging non-overlapping
• Quasistatic variant considered in the following paper:
• https://doi.org/10.1016/j.cma.2017.02.006

Tension Specimen Problem

• Geometry: See right. Larger radius = 0.00635m. Height = 0.04445m.
• Initial condition: 0 disp, 0 velo
• Boundary conditions:
• 0 x-disp in positive x (symmetry BC)
• 0 y-disp in positive y (symmetry BC)
• 0 z-disp at bottom (symmetry BC)
• z-disp = 0.005 * 0.5 * (1.0 - cos(pi*t)) at top (pull from top)
• Notes:
• Single domain not in repo but can be set up (ask Irina)
• Can run linear or nonlinear.
• Will exhibit necking if run with inelastic material model (not yet available in Norma but in progress)
• Variant considered in the following papers:
• https://doi.org/10.1002/nme.6982
• In-prep Schwarz + OpInf paper

Torsion Problem

Laser Weld Problem

• Geometry: See right. Dimensions are 0.1cm x 0.3cm x 0.1cm
• Initial condition: 0 disp, 0 velo
• Boundary conditions:
• Pulling in pos/neg y at +/- 0.1 * (0.5 - 0.5 * cos(pi * t)
• Clamped or symmetry BCs
• Notes:
• Studied more quasistatically than dynamically
• Symmetry BC case seems easier/more robust
• Not all cases have symmetry and clamped BCs in repo, but can be set up
• Difficulties converging non-overlap
• Only >2 subdomain problem set up in Norma currently (there are 3 subdomains)
• Issues converging non-overlapping Schwarz.
• Variant of this problem considered in the following paper:
• https://doi.org/10.1002/nme.6982

​

Bolted Joint Problem

• Geometry: See right. Plate is 0.127m x 0.127m. Height = 0.07m.
• Initial condition: 0 disp, 0 velo
• Boundary conditions:
• 0 disp at base plate
• x-disp = 4.0 * 0.5 * (1.0 - cos(pi*t / (2.0e-3))) on top side of “parts” domain (see arrow in image at right).
• Notes:
• Non-overlapping version implies different physics: bolts are “sticky”. Still looking at how to avoid this.
• Variant of this problem considered in the following paper:
• https://doi.org/10.1002/nme.6982
• In-prep Schwarz + OpInf paper

​

​

​

Plate Problem

Bracket Problem

Pressure Vessel Problem

12

• Geometry: See right. Height = 0.0837794 m, radius = 0.04445 m.
• Initial condition: 0 disp, 0 velo.
• Boundary conditions:
• Symmetry BC in x and y
• Disp z = 0 at 1 point at top
• Pressure (Neumann) condition inside can: p = -937.5e3 * 1500 * (0.5 - 0.5 * cos(pi * t / 1500.0))
• Notes:
• Need to implement pressure NBC for ROM in order to set up this problem for ROM couplings.
• DD set up for Schwarz is arbitrary and should be modified to get better demonstration of coupling.
• Variant of this problem considered in the following paper:
•  https://doi.org/10.1002/nme.7056

​