Some Numerical Considerations for Shock Simulations of Ceramics
Abstract
A comparison of the analytic solution for interface velocity at a planer impact interface with a momentum-conserving numerical contact algorithm suggests a guideline for setting the relative element sizes on either side of the interface. Scaling the element size by the sound speed is shown to produce an initial interface velocity consistent with the Rankine-Hugoniot solution. This introduces a smaller initial perturbation into the numerical simulation, and less noise is generated and propagated from the contact surface. The damping of the subsequent residual numerical noise by the traditional artificial bulk viscosity construct is also briefly examined. Damping by the linear artificial bulk viscosity compellingly reduces noise in elastic materials for coefficients greater than 0.1. A coefficient of 0.2 substantially reduces the noise, but viscous damping this strong will noticeably broaden the shock front.
Document Details
- Document Type
- Technical Report
- Publication Date
- Jun 10, 2020
- Accession Number
- AD1102564
Entities
People
- Richard Becker
Organizations
- United States Army Research Laboratory