Publications

A method is described for applying a sequence of peridynamic models with different length scales concurrently to subregions of a body. The method allows the smallest length scale, and therefore greatest spatial resolution, to be focused on evolving defects such as cracks. The peridynamic horizon in each of the models is half of that of the next model in the sequence. The boundary conditions on each model are provided by the solution predicted by the model above it. Material property characterization for each model is derived by coarse-graining the more detailed resolution in the model below it. Implementation of the multiscale method in the PDMS code is described. Examples of crack growth modeling illustrate the ability of the method to reproduce the main features of crack growth seen in a model with uniformly small resolution. Comparison of the multiscale model results with XFEM and cohesive elements is also given for a crack growth problem.

Surface effects are critical to the accurate simulation of electromagnetics (EM) as current tends to concentrate near material surfaces. Sandia EM applications, which include exploding bridge wires for detonator design, electromagnetic launch of flyer plates for material testing and gun design, lightning blast-through for weapon safety, electromagnetic armor, and magnetic flux compression generators, all require accurate resolution of surface effects. These applications operate in a large deformation regime, where body-fitted meshes are impractical and multimaterial elements are the only feasible option. State-of-the-art methods use various mixture models to approximate the multi-physics of these elements. The empirical nature of these models can significantly compromise the accuracy of the simulation in this very important surface region. We propose to substantially improve the predictive capability of electromagnetic simulations by removing the need for empirical mixture models at material surfaces. We do this by developing an eXtended Finite Element Method (XFEM) and an associated Conformal Decomposition Finite Element Method (CDFEM) which satisfy the physically required compatibility conditions at material interfaces. We demonstrate the effectiveness of these methods for diffusion and diffusion-like problems on node, edge and face elements in 2D and 3D. We also present preliminary work on h -hierarchical elements and remap algorithms.