Publications

Abstract not provided.

Abstract not provided.

The strength of brittle porous media is of concern in numerous applications, for example, earth penetration, crater formation, and blast loading. Thus, it is of importance to possess techniques that allow for constitutive model calibration within the laboratory setting. The goal of the current work is to demonstrate an experimental technique allowing for strength assessment of porous media subjected to shock loading, which can be implemented into pressure-dependent yield surfaces within numerical simulation schemes. As a case study, the deviatoric response of distended α-SiO2 has been captured in a tamped Richtmyer-Meshkov instability (RMI) environment at a pressure regime of 4-10 GPa. Hydrocode simulations were used to interpret RMI experimental data, and a resulting pressure-dependent yield surface akin to the often employed modified Drucker-Prager model was calibrated. Simulations indicate that the resulting jet length generated by the RMI is sensitive to the porous media strength, thereby providing a feasible experimental platform capable of capturing the pressurized granular deviatoric response. Furthermore, in efforts to validate the RMI-calibrated strength model, a set of Mach-lens experiments was performed and simulated with the calibrated pressure-dependent yield surface. Excellent agreement between the resulting Mach-lens length in experiment and simulation provides additional confidence to the RMI yield-surface calibration scheme.

Mesoscale simulations of a LiF impactor colliding with a PMMA capsule containing aluminum powder (r00 = 1.5 g/cc) have been performed to investigate shock-induced melting in porous metals. Impact velocities of 1-2.5 km/s are chosen to coincide with in situ X-ray diffraction experiments, which provide direct evidence of shock-induced melting in aluminum powders. Mesoscale simulations show shock heating within the powder is highly nonuniform and melting remains incomplete over hundreds of nanoseconds behind the shock front despite equilibrium pressure-temperature states from continuum simulations lying above the experimental melt line. Such incomplete melting behavior is consistent with X-ray diffraction data obtained in experiment. For an impact velocity of ~1 km/s, mesoscale simulations predict re-solidification behind the shock front as high-temperature regions are cooled below the melt line. Reducing the grain size of the powder by a factor of two leads to a reduction in the time required to reach complete melt such that total melting of the powder may be observed experimentally for an impact velocity of 2.42 km/s.

Abstract not provided.

Abstract not provided.

A novel experimental methodology is presented to study the deviatoric response of powders in shock regimes. The powders are confined to a cylindrical wedge volume, and a projectile-driven shock wave with a sinusoidally varying front propagates through the powder. The perturbed shock wave exhibits a damping behavior due to irreversible processes of viscosity and strength (deviatoric) of the powder with propagation through increasing powder thicknesses. The inclined surface of the wedge is polished and coated to establish a diffuse surface suitable for reflecting incident laser light into a high-speed camera imaging at 5 MHz. Images of the contrast loss upon shock wave arrival at the observation surface are post-processed for qualitative and quantitative information. New data of shock damping behavior with parameters of perturbation wavelength and initial shock strength are presented for powders of copper, tantalum, and tungsten carbide as well as their mixtures. We present the first full-field images showing additional spatial disturbances on the perturbed shock front that appear dependent on particle material and morphology.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Understanding and quantifying the uncertainties in experimental results are crucial to properly interpreting simulations based on those results. While methods are reasonably well established for estimating those uncertainties in high-pressure shock experiments on homogeneous materials, it is much more difficult to treat relatively low-pressure experiments where shock rise times are significant and material strength is not negligible. Sample heterogeneity further complicates the issue, especially when that heterogeneity is not characterized in each sample. Here, we extend the Monte Carlo impedance matching approach used in high-pressure Z experiments to low-pressure experiments on heterogeneous porous materials. The approach incorporates uncertainties not only in the equation of state of the impedance matching standard but also those associated with its strength. In addition, we also examine approaches for determining material heterogeneity and evaluate its effect on the experimental results.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The shock Hugoniot for full-density and porous CeO2 was investigated in the liquid regime using ab initio molecular dynamics (AIMD) simulations with Erpenbeck's approach based on the Rankine-Hugoniot jump conditions. The phase space was sampled by carrying out NVT simulations for isotherms between 6000 and 100 000 K and densities ranging from ρ=2.5 to 20g/cm3. The impact of on-site Coulomb interaction corrections +U on the equation of state (EOS) obtained from AIMD simulations was assessed by direct comparison with results from standard density functional theory simulations. Classical molecular dynamics (CMD) simulations were also performed to model atomic-scale shock compression of larger porous CeO2 models. Results from AIMD and CMD compression simulations compare favorably with Z-machine shock data to 525 GPa and gas-gun data to 109 GPa for porous CeO2 samples. Using results from AIMD simulations, an accurate liquid-regime Mie-Grüneisen EOS was built for CeO2. In addition, a revised multiphase SESAME-Type EOS was constrained using AIMD results and experimental data generated in this work. This study demonstrates the necessity of acquiring data in the porous regime to increase the reliability of existing analytical EOS models.

Abstract not provided.

Abstract not provided.

CTH is an Eulerian hydrocode developed by Sandia National Laboratories (SNL) to solve a wide range of shock wave propagation and material deformation problems. Adaptive mesh refinement is also used to improve efficiency for problems with a wide range of spatial scales. The code has a history of running on a variety of computing platforms ranging from desktops to massively parallel distributed-data systems. For the Trinity Phase 2 Open Science campaign, CTH was used to study mesoscale simulations of the hypervelocity penetration of granular SiC powders. The simulations were compared to experimental data. A scaling study of CTH up to 8192 KNL nodes was also performed, and several improvements were made to the code to improve the scalability.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We have recently shown that the final density of silicon under shock compression is anomalously enhanced by introducing voids in the initial uncompressed material. Using molecular simulation, we also demonstrated a molecular mechanism for the effect, which is seen in a growing class of other similar materials. We have shown that this mechanism involves a premature local phase transition nucleated by local shear strain. At higher shock loads we show here that this transition becomes frustrated producing amorphous silicon.We also observe local melting below the equilibrium melt line for bulk silicon. Large-scale non-equilibrium molecular dynamics (NEMD) and Hugoniostat simulations of shock compressed porous silicon are used to study the mechanism. Final stress states and strength were characterized versus initial porosity and for various porosity microstructures.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A technique in which the evolution of a perturbation in a shock wave front is monitored as it travels through a sample is applied to granular materials. Although the approach was originally conceived as a way to measure the viscosity of the sample, here it is utilized as a means to probe the deviatoric strength of the material. Initial results for a tungsten carbide powder are presented that demonstrate the approach is viable. Simulations of the experiments using continuum and mesoscale modeling approaches are used to better understand the experiments. The best agreement with the limited experimental data is obtained for the mesoscale model, which has previously been shown to give good agreement with planar impact results. The continuum simulations indicate that the decay of the perturbation is controlled by material strength but is insensitive to the compaction response. Other sensitivities are assessed using the two modeling approaches. The simulations indicate that the configuration used in the preliminary experiments suffers from certain artifacts and should be modified to remove them. The limitations of the current instrumentation are discussed, and possible approaches to improve it are suggested.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Under shock compression, most porous materials exhibit lower densities for a given pressure than that of a full-dense sample of the same material. However, some porous materials exhibit an anomalous, or enhanced, densification under shock compression. We demonstrate a molecular mechanism that drives this behavior. We also present evidence from atomistic simulation that silicon belongs to this anomalous class of materials. Atomistic simulations indicate that local shear strain in the neighborhood of collapsing pores nucleates a local solid-solid phase transformation even when bulk pressures are below the thermodynamic phase transformation pressure. This metastable, local, and partial, solid-solid phase transformation, which accounts for the enhanced densification in silicon, is driven by the local stress state near the void, not equilibrium thermodynamics. This mechanism may also explain the phenomenon in other covalently bonded materials.

A nonlocal, ordinary peridynamic constitutive model is formulated to numerically simulate the pressure-dependent flow and fracture of heterogeneous, quasi-brittle ma- terials, such as concrete. Classical mechanics and traditional computational modeling methods do not accurately model the distributed fracture observed within this family of materials. The peridynamic horizon, or range of influence, provides a characteristic length to the continuum and limits localization of fracture. Scaling laws are derived to relate the parameters of peridynamic constitutive model to the parameters of the classical Drucker-Prager plasticity model. Thermodynamic analysis of associated and non-associated plastic flow is performed. An implicit integration algorithm is formu- lated to calculate the accumulated plastic bond extension and force state. The gov- erning equations are linearized and the simulation of the quasi-static compression of a cylinder is compared to the classical theory. A dissipation-based peridynamic bond failure criteria is implemented to model fracture and the splitting of a concrete cylinder is numerically simulated. Finally, calculation of the impact and spallation of a con- crete structure is performed to assess the suitability of the material and failure models for simulating concrete during dynamic loadings. The peridynamic model is found to accurately simulate the inelastic deformation and fracture behavior of concrete during compression, splitting, and dynamically induced spall. The work expands the types of materials that can be modeled using peridynamics. A multi-scale methodology for simulating concrete to be used in conjunction with the plasticity model is presented. The work was funded by LDRD 158806.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Density Functional Theory (DFT) based molecular dynamics has been established as a method capable of yielding high fidelity results for many materials at a wide range of pressures and temperatures and has recently been applied to complex polymers such as polyethylene, compounds such as ethane or CO2, and oxides such as MgO. We use this method to obtain a Grïneisen Γ and thereby build a Mie-Grüneisen equation of state (EOS) and a Rice-Walsh EOS for tantalum pentoxide (Ta2O5 or tantala) and compare to experimental data. The experimental data have initial densities (ρ00) of approximately 1.13, 3, and 7.4 g/cm 3 reduced from a crystalline of 8.36 g/cm3. We found that r becomes constant at higher temperatures and pressure, but is a function of both density and temperature at lower densities and temperatures. Finally, the Mie-Gruneisen EOS is adequate for modeling the slightly distended Hugoniot with an initial density of 7.4 g/cm3 however it is inadequate for the more porous Hugoniot, while the Rice-Walsh EOS combined with a P-λ crush model approximates the experimental data quite well. © Published under licence by IOP Publishing Ltd.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.