Engineering topics
Herbold, E.
Publications and source records attributed to Herbold, E..
Pre-Shot Simulations of Near-Field and Far-Field Ground Motion for the Source Physics Experiment (SPE) Explosions at the Climax Stock, Nevada National Security Site: SPE3
The Source Physics Experiment (SPE) is designed to improve ground motion prediction capabilities for underground nuclear testing through the development of new physics-based models together with validation via new chemical explosion data. The SPE shots consist of different size and depth chemical explosions in a granite borehole at the Climax Stock on the Nevada National Security Site (NNSS). The third shot in the series (SPE3) will be the same size (1000 kg TNT equivalent) and depth (45 m canister centroid) as the previous shot (SPE2). The main purpose of the SPE3 is to investigate the effects of the damage zone caused by the SPE2 chemical explosion on the SPE3 explosion, keeping all other parameters (e.g. size, depth, explosive type) the same. This SPE3 pre-shot report addresses progress made in simulation of the expected wavefield in the both the near- and the far-field. In this report we consider the far-field to be where seismic wave motion can be characterized by elastic or visco-elastic approximations, in contrast to the near-field where other, often non-linear effects, play a major role. The current simulations build on previous work from the SPE1 and SPE2 shots. The simulations are intended to highlight differences between SPE2 and SPE3 using a variety to approximations intended to mimic the expected effect of the damage zone. The simplest physics of damage zone effects suggest changes in the source region material properties (density, velocity), and in the yield strength (coupling) of the rock, should both impact the far-field signal. We can illustrate this in a simplistic way for P-waves in Figure A using the empirical Mueller-Muphy (1971) far-field explosion model (MM). First if the SPE2 explosion-induced damage reduces velocities and densities without affecting the yield strength or coupling, then by energy conservation the amplitudes for SPE 3 would increase as shown in the top of Figure A. Similarly if the explosion induced damage increased density (say by compacting pore space) the amplitude would decrease. This impedance effect, where the product of the density, wave-velocity, and ground-velocity squared is conserved, is well known. The second and harder effect to quantify, are changes in the energy coupling or yield strength. In the simple MM model this is represented by the coupling factor between cavity radius and yield. If we assume the granite was slightly weakened by SPE2 such that the cavity radius increases slightly then the overall amplitudes of SPE3 decrease relative to SPE2 as shown in the bottom of Figure A. In reality both effects, may exist due to the damage zone, and furthermore we know they are interrelated in complex ways. For S-waves we expect even more complex effects due to motion on joints and faults. For this reason we are developing physics-based near field simulation models coupled to far-field propagation codes to predict the overall effects of changes in the source emplacement region on observed signals. Near-field simulations have been performed using the Eulerian code (GEODYN). An ensemble model, calibrated to match records from past large-scale nuclear events is used. Both 2D axisymmetric (no faults) and fully three-dimensional GEODYN calculations (with /without faults) were conducted in which the damage zones sustained by SPE2 are retained as the initial condition for the SPE3 calculations. The results indicate that the waveform peak amplitudes are slightly decayed since the ensemble model includes no softening mechanism. The pulse widths are reduced, indicating the energy is dissipated more by the damage zone as expected. The dominant motions calculated in all the scenarios are spherically symmetric and non-spherical components are small as seen for the SPE2 calculations. The arrival time of the SPE3 waveform at the boreholes is shorter than SPE2 due to compaction during the previous explosion. As shown in the SPE 1 and 2 observations shear waves are generated, indicating the potential contributions from the induced motion on distributed joints. The effects of joints are evaluated using the Lagrangian hydrocode (GEODYN-L), which was benchmarked against the GEODYN calculations for the ensemble model. A new granite model calibrated for samples from the SPE site is used. The new model includes effects of poroelasticity and water saturation and may be scaled allowing application to both small and large scales. Parametric tests were performed to evaluate the effect of joints within the granite and the results demonstrate that in situ joints introduce anisotropy and significant tangential motions.
Dynamic Behavior of Sand: Annual Report FY 11
Currently, design of earth-penetrating munitions relies heavily on empirical relationships to estimate behavior, making it difficult to design novel munitions or address novel target situations without expensive and time-consuming full-scale testing with relevant system and target characteristics. Enhancing design through numerical studies and modeling could help reduce the extent and duration of full-scale testing if the models have enough fidelity to capture all of the relevant parameters. This can be separated into three distinct problems: that of the penetrator structural and component response, that of the target response, and that of the coupling between the two. This project focuses on enhancing understanding of the target response, specifically granular geomaterials, where the temporal and spatial multi-scale nature of the material controls its response. As part of the overarching goal of developing computational capabilities to predict the performance of conventional earth-penetrating weapons, this project focuses specifically on developing new models and numerical capabilities for modeling sand response in ALE3D. There is general recognition that granular materials behave in a manner that defies conventional continuum approaches which rely on response locality and which degrade in the presence of strong response nonlinearities, localization, and phase gradients. There are many numerical tools available to address parts of the problem. However, to enhance modeling capability, this project is pursuing a bottom-up approach of building constitutive models from higher fidelity, smaller spatial scale simulations (rather than from macro-scale observations of physical behavior as is traditionally employed) that are being augmented to address the unique challenges of mesoscale modeling of dynamically loaded granular materials. Through understanding response and sensitivity at the grain-scale, it is expected that better reduced order representations of response can be formulated at the continuum scale as illustrated in Figure 1 and Figure 2. The final result of this project is to implement such reduced order models in the ALE3D material library for general use.