Gaining a proper understanding of how Earth structure and other near-source properties affect estimates of explosion yield is important to the nonproliferation mission. The yields of explosion sources are often based on seismic moment or waveform amplitudes. Quantifying how the seismic waveforms or estimates of the source characteristics derived from those waveforms are influenced by natural or man-made structures within the near-source region, where the wavefield behaves nonlinearly, is required to understand the full range of uncertainty in those yield estimates. We simulate tamped chemical explosions using a nonlinear, shock physics code and couple the ground motions beyond the elastic radius to a linear elastic, full waveform seismic simulation algorithm through 3D media. In order to isolate the effects of simple small-scale 3D structures on the seismic wavefield and linear seismic source estimates, we embed spheres and cylinders close to the fully- tamped source location within an otherwise homogenous half-space. The 3 m diameters spheres, given their small size compared to the predominate wavelengths investigated, not surprisingly are virtually invisible with only negligible perturbations to the far-field waveforms and resultant seismic source time functions. Similarly, the 11 m diameter basalt sphere has a larger, but still relatively minor impact on the wavefield. However, the 11 m diameter air-filled sphere has the largest impact on both waveforms and the estimated seismic moment of any of the investigated cases with a reduction of ~25% compared to the tamped moment. This significant reduction is likely due in large part to the cavity collapsing from the shock instead of being solely due to diffraction effects . Although the cylinders have the same diameters as the 3 m spheres, their length of interaction with the wavefield produces noticeable changes to the seismic waveforms and estimated source terms with reductions in the peak seismic moment on the order of 10%. Both the cylinders and 11 m diameter spheres generate strong shear waves that appear to emanate from body force sources.
We used the CTH shock physics code to simulate the explosion of an 18-t chemical explosive at a depth of 250 m. We used the CTH in the two-dimensional axisymmetric (cylindrical) geometry (2DC) and most simulations included fully tamped explosions in wet tuff. Our study focused on parametric studies of three of the traditional strength models available in CTH, namely, geologic-yield, elastic perfectly-plastic von Mises, and Johnson-Cook strength (flow stress) models. We processed CTH results through a code that generates Reduced Displacement Potential (RDP) histories for each simulation. Since RDP is the solution of the linear wave equation in spherical coordinates, it is mainly valid at far-enough distance from the explosion the elastic radius. Among various parameters examined, we found the yield strength to have the greatest effect on the resulting RDP, where the peak RDP reduces almost linearly in log-log space as the yield strength increases. Moreover, an underground chemical explosion results in a cavity whose final diameter is inversely proportional to the material yield strength, i.e., as the material's yield strength increases the resulting final cavity radius decreases. Additionally, we found the choice of explosive material (COMP-C4 versus COMP-B) has minor effects on the peak RDP, where denser COMP-C4 shows higher peak RDP than the less dense COMP-B by a factor of ~1.1. In addition to wet tuff, we studied explosions in dry tuff, salt, and basalt, for a single strength model and yield strength value. We found wet tuff has the highest peak RDP value, followed by dry tuff, salt, and basalt. 2DC simulations of explosions in 11 m radius spherical, hemispherical, and cylindrical cavities showed the RDP signals have much lower magnitude than tamped explosions, where the cavity explosions mimicked nearly decoupled explosions.
Seismic source modeling allows researchers both to simulate how a source that induces seismic waves interacts with the Earth to produce observed seismograms and, inversely, to infer what the time histories, sizes, and force distributions were for a seismic source given observed seismograms. In this report, we discuss improvements made in FY21 to our software as applies to both the forward and inverse seismic source modeling problems. For the forward portion of the problem, we have added the ability to use full 3-D nonlinear simulations by implementing 3-D time varying boundary conditions within Sandia’s linear seismic code Parelasti. Secondly, on the inverse source modeling side, we have developed software that allows us to invert seismic gradiometer-derived observations in conjunction with standard translational motion seismic data to infer properties of the source that may improve characterization in certain circumstances. First, we describe the basic theory behind each software enhancement and then demonstrate the software in action with some simple examples.
Underground explosions nonlinearly deform the surrounding earth material and can interact with the free surface to produce spall. However, at typical seismological observation distances the seismic wavefield can be accurately modeled using linear approximations. Although nonlinear algorithms can accurately simulate very near field ground motions, they are computationally expensive and potentially unnecessary for far field wave simulations. Conversely, linearized seismic wave propagation codes are orders of magnitude faster computationally and can accurately simulate the wavefield out to typical observational distances. Thus, devising a means of approximating a nonlinear source in terms of a linear equivalent source would be advantageous both for scenario modeling and for interpretation of seismic source models that are based on linear, far-field approximations. This allows fast linear seismic modeling that still incorporates many features of the nonlinear source mechanics built into the simulation results so that one can have many of the advantages of both types of simulations without the computational cost of the nonlinear computation. In this report we first show the computational advantage of using linear equivalent models, and then discuss how the near-source (within the nonlinear wavefield regime) environment affects linear source equivalents and how well we can fit seismic wavefields derived from nonlinear sources.
A conceptual model of the MIU site in central Japan, was developed to predict the groundwater system response to pumping. The study area consisted of a fairly large three-dimensional domain, having the size 4.24 x 6 x 3 km{sup 3} with three different geological units, upper and lower fractured zones and a single fault unit. The resulting computational model comprised of 702,204 finite difference cells with variable grid spacing. Both steady-state and transient simulations were completed to evaluate the influence of two different surface boundary conditions: fixed head and no flow. Steady state results were used for particle tracking and also serving as the initial conditions (i.e., starting heads) for the transient simulations. Results of the steady state simulations indicate the significance of the choice of surface (i.e., upper) boundary conditions and its effect on the groundwater flow patterns along the base of the upper fractured zone. Steady state particle tracking results illustrate that all particles exit the top of the model in areas where groundwater discharges to the Hiyoshi and Toki rivers. Particle travel times range from 3.6 x 10{sup 7} sec (i.e., {approx}1.1 years) to 4.4 x 10{sup 10} sec (i.e., {approx}1394 years). For the transient simulations, two pumping zones one above and another one below the fault are considered. For both cases, the pumping period extends for 14 days followed by an additional 36 days of recovery. For the pumping rates used, the maximum drawdown is quite small (ranging from a few centimeters to a few meters) and thus, pumping does not severely impact the groundwater flow system. The range of drawdown values produced by pumping below the fault are generally much less sensitive to the choice of the boundary condition than are the drawdowns resulted from the pumping zone above the fault.
We consider the use of a hypodiffusive governing equation (HDE) for the porous-continuum modeling of gravity-driven fingers (GDF) as occur in initially dry, highly nonlinear, and hysteretic porous media. In addition to the capillary and gravity terms within the traditional Richards equation, the HDE contains a hypodiffusive term that models an experimentally observed hold-back-pile-up (HBPU) effect and thus imparts nonmonotonicity at the wetting front. In its dimensionless form the HDE contains the dimensionless hypodiffusion number, NHD. As NHD increases, one-dimensional (1D) numerical solutions transition from monotonic to nonmonotonic. Considering the experimentally observed controls on GDF occurrence, as either the initial moisture content and applied flux increase or the material nonlinearity decreases, solutions undergo the required transition back to monotonic. Additional tests for horizontal imbibition and capillary rise show the HDE to yield the required monotonie response but display sharper fronts for NHD > 0. Finally, two-dimensional (2D) numerical solutions illustrate that in parameter space where the 1D HDE yields nonmonotonicity, in 2D it forms nonmonotonic GDF.
The authors consider the ability of the numerical solution of Richards equation to model gravity-driven fingers. Although gravity-driven fingers can be easily simulated using a partial downwind averaging method, they find the fingers are purely artificial, generated by the combined effects of truncation error induced oscillations and capillary hysteresis. Since Richards equation can only yield a monotonic solution for standard constitutive relations and constant flux boundary conditions, it is not the valid governing equation to model gravity-driven fingers, and therefore is also suspect for unsaturated flow in initially dry, highly nonlinear, and hysteretic media where these fingers occur. However, analysis of truncation error at the wetting front for the partial downwind method suggests the required mathematical behavior of a more comprehensive and physically based modeling approach for this region of parameter space.