Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPLOSIVE”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,477 records · Page 82

Powder Diffraction Crystallography : 98 Years as Plan B?

A historical essay written in 2014 to accompany a talk in the Transactions Symposium "100 Years of Crystallography" at the 2014 American Crystallographic Association Annual Meeting (talk at https://history.amercrystalassn.org/h-brian-toby-videos). Abstract: Powder diffraction patterns were used for structure solution almost contemporaneously with the Bragg’s initial single-crystal determinations. Structures have continued to be determined from powders, in the subsequent century. There has been a nearly explosive growth in this technique since the 1970’s. This was initially driven by neutron diffraction, which prompted the invention of Rietveld analysis, but synchrotrons have gone on to have at least as large an impact. There is no apparent upper limit to the scale of systems that can be structurally studied with powder diffraction; novel protein structures have been determined. When short-range ordering or disordered regions of a material are of interest, the Fourier transform of a powder diffraction pattern, the pair distribution function can be modeled.

Toby, Brian H.↗

Neutronics Modeling of the pulsed plasma rocket reactor using rattlesnake

In the pulsed plasma rocket (PPR) reactor, neutron induced fission processes are utilized to implement a series of pulsed micro-explosions of very high power and eject plasma as a propellant. More specifically, projectiles (bullets) composed of moderated uranium are sent through the chamber of an unmoderated uranium barrel. By inducing rapid fission within the bullets, a plasma can be generated with appropriate delivery of neutrons. This study aims to assess the neutronics performance of the PPR reactor during normal operations using the MOOSE-based Rattlesnake code through the evaluation of the impact of the movement of the fuel bullet and the rotation of control drums on the criticality of the system. The Monte Carlo (MC) code Serpent 2 was employed to generate material-based cross sections for use in Rattlesnake and the reference neutronics solution. Cubit was used to generate the mesh for the Rattlesnake model. Parametric studies were conducted to evaluate the best approaches for cross section and mesh generation to ensure accurate results from Rattlesnake. As part of the verification process, the eigenvalue results of the system at various fuel bullet positions were obtained using Rattlesnake and compared with the reference solutions. The acceptable differences show that the Rattlesnake model with appropriate cross section and mesh generation procedure is a sufficiently accurate approximation of the continuous energy (CE) MC model for micro-sized reactors like the PPR reactor. Next, the impact of the rotation of the control drums on the system criticality was evaluated using the verified model. It was shown that the high-fidelity simulation using the deterministic code Rattlesnake can produce sufficiently accurate results for the evaluation of the reactor's neutronics performance in different phases of the normal operation with acceptable computational cost.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Performing Numerical Analysis of Cybersecurity Options Using Dynamic Risk Analysis Tool EMRALD

Cyberattacks can have many different attack paths, durations, and goals. There are also many different mitigation options involving hardware, software, and/or humans. Considering a cyber threat should involve defense-in-depth methods and a quantitative or numerical evaluation of overall effectiveness against dynamic, time-dependent attacks to make cost and risk-informed decisions. Typical cyberattack modeling methods only provide a qualitative evaluation. The main areas of cybersecurity are confidentiality, integrity, and availability. For companies with cyber-physical systems such as advanced nuclear reactors, cyber-related safety is a requirement set by North American Electric Reliability and the U.S. Nuclear Regulatory Commission. They are also concerned about availability or reliability as a business case. As cyber threats are evolving to a business-for-hire structure, more attacks may focus on disrupting business success and reliability, causing financial and economic stability risk. Companies want to know business reliability and recovery from those threats, and that requires modeling physical behavior of the targets. Dynamic-state-based and Markov-based modeling provides a method for better cyber scenario modeling with different tools having issues such as state-base explosion. Dynamic modeling enables time and conditional features not found in other numerical evaluation methods. EMRALD (Event Modeling Risk Assessment using Lined Diagrams) is a dynamic risk analysis modeling and simulation tool and has features that reduce modeling issues. It has been used to model different time-dependent events including plant behavior and operator procedures. As a general modeling tool, EMRALD can also be used to model cyberattack scenarios with varying mitigation options and quantify effectiveness, producing numerical data for risk-informed decisions. This paper uses EMRALD to demonstrate that dynamic numerical risk analysis can be used for cyber threat modeling to provide insights for design decision-making and optimize defense strategies. Keywords: cyber modeling; cyber-physical systems; numerical cyber modeling

97 - MATHEMATICS AND COMPUTING↗

Three-beam rotational coherent anti-Stokes Raman spectroscopy thermometry in scattering environments

Three-beam rotational coherent anti-Stokes Raman scattering (CARS) measurements performed in highly scattering environments are susceptible to contamination by two-beam CARS signals generated by the pump–probe and Stokes–probe interactions at the measurement volume. If this occurs, differences in the Raman excitation bandwidth between the two-beam and three-beam CARS signals can add significant errors to the spectral analysis. This interference to the best of our knowledge has not been acknowledged in previous three-beam rotational CARS experiments, but may introduce measurement errors up to 25% depending on the temperature, amount of scattering, and differences between the two-beam and three-beam Raman excitation bandwidths. Here, the presence of two-beam CARS signal contamination was experimentally verified using a femtosecond–picosecond rotational CARS instrument in two scattering environments: (1) a fireball generated by a laboratory-scale explosion that contained particulate matter, metal fragments, and soot, and (2) a flow of air and small liquid droplets. A polarization scheme is presented to overcome this interference. By rotating the pump and Stokes polarizations +55° and –55° from the probe, respectively, the two-beam and three-beam CARS signals are orthogonally polarized and can be separated using a polarization analyzer. Using this polarization arrangement, the Raman-resonant three-beam CARS signal amplitude is reduced by a factor of 2.3 compared to the case where all polarizations are parallel. This method is successfully demonstrated in both scattering environments. A theoretical model is presented, and the temperature measurement error is studied for different experimental conditions. The criteria for when this interference may be present are discussed.

42 ENGINEERING↗

Time-stretch spectroscopy for fast infrared absorption spectra of acetylene and hydroxyl radicals during combustion

We have developed a diagnostic that uses time-domain spectroscopy to measure transient infrared absorption spectra in gases. Using a time-stretch Fourier transform approach, we can determine pressure, temperature, and gas concentrations with sub-microsecond time resolution for over two milliseconds. We demonstrate high-resolution (0.015 nm), time-resolved spectral measurements in an acetylene-oxygen gas mixture undergoing combustion. Within a 5 µs period during the reaction, the acetylene line intensities decrease substantially, and new spectra appear that are consistent with the hydroxyl (OH) radical, a common by-product in the combustion, deflagration, and detonation of fuels and explosives. Post-reaction pressures and temperatures were estimated from the OH spectra. The technique measures spectra from 1520 to 1620 nm using fiber optics, photodetectors, and digitizers. No cameras or spectrometers are required.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Ultraviolet digital holographic microscopy (DHM) of micron-scale particles from shocked Sn ejecta

A cloud of very fast, O (km/s), and very fine, O (µm), particles may be ejected when a strong shock impacts and possibly melts the free surface of a solid metal. To quantify these dynamics, this work develops an ultraviolet, long-working distance, two-pulse Digital Holographic Microscopy (DHM) configuration and is the first to replace film recording with digital sensors for this challenging application. A proposed multi-iteration DHM processing algorithm is demonstrated for automated measures of the sizes, velocities, and three-dimensional positions of non-spherical particles. Ejecta as small as 2 µm diameter are successfully tracked, while uncertainty simulations indicate that particle size distributions are accurately quantified for diameters ≥4 µm. These techniques are demonstrated on three explosively driven experiments. Measured ejecta size and velocity statistics are shown to be consistent with prior film-based recording, while also revealing spatial variations in velocities and 3D positions that have yet to be widely investigated. Having eliminated time-consuming analog film processing, the methodologies proposed here are expected to significantly accelerate future experimental investigation of ejecta physics.

47 OTHER INSTRUMENTATION↗

Building epsilon near zero materials from layered uniaxial metamaterials

Recently, there has been an explosion of activity in the fields of optics and photonics with the advent of fabrication techniques which enable the design of metamaterials which possess properties not encountered in the natural world. In this work, we are concerned with zero permittivity materials and a new scheme to design metamaterials for which all components of the dielectric tensor are approximately zero. Our approach involves the alternate layering of many, very thin, slices of two constituent metamaterials, a uniaxial layered medium and a uniaxial nanowire array. With a simple optimization strategy we demonstrate a candidate configuration which very nearly satisfies our design goal of zero permittivity.

Gray, Stephen K. [Argonne National Laboratory (ANL↗

Dose estimates and their uncertainties for use in epidemiological studies of radiation-exposed populations in the Russian Southern Urals

Many residents of the Russian Southern Urals were exposed to radioactive environmental pollution created by the operations of the Mayak Production Association in the mid- 20 th century. There were two major releases: the discharge of about 1x10 17 Bq of liquid waste into the Techa River between 1949 and 1959; and the atmospheric release of 7.4 * 10 16 Bq as a result an explosion in the radioactive waste-storage facility in 1957. The releases into the Techa River resulted in the exposure of more than 30,000 people who lived in riverside villages between 1950 and 1961. The 1957 accident contaminated a larger area with the highest exposure levels in an area that is called the East Urals Radioactive Trace (EURT). Current epidemiologic studies of the exposed populations are based on dose estimates obtained using a Monte-Carlo dosimetry system (TRDS-2016MC) that provides multiple realizations of the annual doses for each cohort member. These dose realizations provide a central estimate of the individual dose and information on the uncertainty of these dose estimates. In addition, the correlation of individual annual doses over realizations provides important information on shared uncertainties that can be used to assess the impact of shared dose uncertainties on risk estimate uncertainty.This paper considers dose uncertainties in the TRDS-2016MC. Individual doses from external and internal radiation sources were reconstructed for 48,036 people based on environmental contamination patterns, residential histories, individual 90 Sr body-burden measurements and dietary intakes. Dietary intake of 90 Sr resulted in doses accumulated in active bone marrow (or simply, marrow) that were an order of magnitude greater than those in soft tissues. About 84% of the marrow dose and 50% of the stomach dose was associated with internal exposures. The lognormal distribution is well-fitted to the individual dose realizations, which, therefore, could be expressed and easily operated in terms of geometric mean (GM) and geometric standard deviation (GSD). Cohort average GM for marrow and stomach cumulative doses are 0.21 and 0.03 Gy, respectively. Cohort average dose uncertainties in terms of GSD are as follows: for marrow it is 2.93 (90%CI: 2.02–4.34); for stomach and the other non-calcified tissues it is 2.32 (90% CI: 1.78–2.9).

61 RADIATION PROTECTION AND DOSIMETRY↗

Solid-state ionics: The key to the discovery and domination of lithium batteries: some learnings from β-alumina and titanium disulfide

We report solid state ionics, the study of fast ion transport in solids, expanded explosively after the discovery of sodium ion transport in b-alumina 50 years ago and has revolutionized energy storage. Lithium ion batteries have come from a dream with titanium disulfide to enabling the communications revolution and are enabling renewable energy. Much can be learnt from these early studies, that will allow the attainment of a fossil-free society and a sustainable environment.

36 MATERIALS SCIENCE↗

Luminescent silica microagglomerates, synthesis, and environmental testing

Abstract Tracking mass through harsh environments requires surrogate particles that withstand the event and endure until sampling. Silica-covered quantum dots have been shown to withstand a range of environmental pHs from months to years; in this work they are shown to endure in anticipated local environments. Two methods of particle synthesis were employed to produce luminescent silica with particle diameters 0.1–4 μm. These tracer particles scale for mass production, tolerate harsh environments, and endure in debris. They could be deployed in places such as chemical explosions, industrial processes, geologic test beds, oil and gas fields, nuclear reactors, and geothermal plants to track mass under harsh conditions. Graphical abstract

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Local-Distance Seismic Event Relocation and Relative Magnitude Estimation, Applications to Mining Related Seismicity in the Powder River Basin, Wyoming

Recent efforts to characterize small ( M w < 3 ) seismic events at local distances have become more important because of the increased observation of human-triggered and induced seismicity and the need to advance nuclear explosion monitoring capabilities. The signals generated by low-magnitude seismic sources necessitate the use of nearby short-period observations, which are sensitive to local geological heterogeneity. Local to near-regional distance ( < 300 km ) surface and shear waves can dominate short-period observations from small, shallow seismic sources. In this work, we utilize these observations to estimate precise, relative locations and magnitudes of ~ 700 industrial mining events in Wyoming, using nearly 360,000 observations. The precise, relative location estimates (with formal location uncertainty estimates of less than 1 km) collapse a diffuse collection of mining events into discrete clusters associated with individual blasting operations. We also invert the cross-correlation amplitudes to estimate precise, relative moment magnitude estimates, which help validate and identify disparities in the event sizes reported by regional network catalogs. Joint use of multiple phases allows for the inclusion of more seismic events due to the increase in the number of observations. In some cases, using a single phase allowed us to relocate only 50% of the original reported seismic events within a cluster. Combining shear- and surface-wave phases increased the number of events to above 90% of the original events, allowing us to characterize a broader range of event sizes, source to station distances, and event distributions. This analysis takes a step toward making a fuller characterization of small industrial seismic events observed at local distances.

58 GEOSCIENCES↗

A New Decade in Seismoacoustics (2010–2022)

Several sources of interest often generate both low‐frequency acoustic and seismic signals due to energy propagation through the atmosphere and the solid Earth. Seismic and acoustic observations are associated with a wide range of sources, including earthquakes, volcanoes, bolides, chemical and nuclear explosions, ocean noise, and others. The fusion of seismic and acoustic observations contributes to a better understanding of the source, both in terms of constraining source location and physics, as well as the seismic to acoustic coupling of energy. Here in this review, we summarize progress in seismoacoustic data processing, including recent developments in open‐source data availability, low‐cost seismic and acoustic sensors, and large‐scale deployments of collocated sensors from 2010 to 2022. Similarly, we outline the recent advancements in modeling efforts for both source characteristics and propagation dynamics. Finally, we highlight the advantages of fusing multiphenomenological signals, focusing on current and future techniques to improve source detection, localization, and characterization efforts. This review aims to serve as a reference for seismologists, acousticians, and others within the growing field of seismoacoustics and multiphenomenology research.

58 GEOSCIENCES↗

SKS Polarization Anomalies Due to the Coriolis Force

The Earth’s Coriolis force has been well–known to impact surface waves and normal modes, which is essential to accurately interpret these waves. However, the Coriolis force on body waves has been assumed to be negligible and mostly ignored. It has been previously shown that the Coriolis force impacts polarizations of shear waves, whereas the wavefronts remain unaffected. We expand on the potential influences of Earth’s Coriolis force on shear–wave polarization measurements by conducting 3D numerical simulations for elastic waves generated by earthquake and explosive sources in a radially symmetric, and 3D mantle and crustal models. The Coriolis force can produce polarization anomalies of mantle shear waves up to 7° and core phases, such as SKS and SKKS, up to 4°. Uncorrected shear–wave polarizations due to the Coriolis force can cause an additional source of error (5°–10° in fast direction, and 0.2–0.3 s delay time depending on the method and seismic phase), inaccurate interpretation of station misalignments, and imprecise estimates of the core–mantle boundary topography. We show how to correct for the Coriolis force on teleseismic shear waves using 1D ray tracing for well–isolated phases. Here, we recommend the use of full waveform simulations to accurately account for earthquake sources parameters, poorly isolated phases that could include interfering phase arrivals within the measurement time window, and the effect of the Coriolis force on the polarizations of shear waves.

58 GEOSCIENCES↗

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure of the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

Here we present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. The new model enables fully 3D simulations of complete three-component waveforms and provides improved fits that were not possible with previous models. We inverted over 32,000 waveforms from 192 earthquakes recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region with highly uneven coverage. Inversion iterations proceeded from the period band 50–100 s in six stages and 54 total iterations reducing the minimum period to 30 s. Our final model, MESWA, improves waveform fits compared to the starting and other models for both the data used in the inversion and an independent validation set of 66 events. Restitution tests indicate that the model resolves features in the central part of the model to depths of about 150 km. The new model reveals tectonic features imaged by other studies and methods but in a new holistic model of anisotropic shear and compressional wavespeeds (V S and V P , respectively) covering a larger domain with smaller scale length and amplified features. Examples include low crustal V S in the Tethyan belt and low mantle VS following divergent (Gulf of Aden, Red Sea) and transform (Dead Sea fault) margins of the Arabian plate. Low V S is imaged below Cenozoic volcanic centers of the Mecca–Madina–Nafud Line, Arabian Peninsula, and the Türkiye–Iran border region. Elevated V S tracks Makran subduction under southeast Iran with near vertical dip. MESWA could be used as a starting model for further improvements, say, using waveforms from in-country seismic networks that are not currently openly available and/or smaller-scale studies targeting a shorter period. The model could be used to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗

Improved Regional Moment Tensor Inversion for Moderately Large Earthquakes in the Western United States Using a 3D Earth Model Based on Full Waveform Tomography

The nature of seismic sources for moderately large (moment magnitude, M w 5.0–6.5) events are commonly characterized by their moment tensor (MT) solutions and obtained by inversion of regional distance (200–1600 km) long‐period (20–50 s) waveforms. Regional MT estimates are often calculated from average plane‐layered, one‐dimensional (1D) velocity models. However, 1D model calculations can produce misfits in the arrival times and waveform shapes that introduce errors, particularly at longer distances or for shorter periods, which are necessary for analyzing lower magnitude events. Approximate Earth models (e.g., 1D) representing broad areas may be inadequate, particularly in the crust and uppermost mantle of tectonically complex regions. In this study, we show how a three‐dimensional (3D) Earth model obtained from full waveform inversion tomography can improve waveform fits and decrease phase errors. We developed a platform and workflow to perform routine 3D MT inversions and inverted MTs for 25 earthquakes in the western United States and seven nuclear explosions using an average 1D and a recent 3D Earth model, WUS256 (Rodgers et al., 2022). Using the 3D model improves waveform fits (variance reduction and phase time shifts) compared with the 1D model, and the 3D MT solutions are stable across large distances. This study shows that 3D models obtained from full waveform tomography can improve MTs and source characterization especially at far regional distances (>800 km).

Geosciences↗

Assessing the Expansion of Ground-Motion Sensing Capability in Smart Cities via Internet Fiber-Optic Infrastructure

Monitoring ground motion in smart cities can improve the public safety by providing critical insights on natural and anthropogenic hazards, for example, earthquakes, landslides, explosions, infrastructure failures, and so forth. Although seismic activity is typically measured using dedicated point sensors (e.g., geophones and accelerometers), techniques such as distributed acoustic sensing have demonstrated the utility of using fiber-optic cable to detect seismic activity over comparable distances. In this article, we present the results of a study that quantifies the expansion in an area monitored for low-amplitude ground-motion events by augmenting existing point sensors with the internet fiber-optic cable infrastructure. Here we begin by describing our methodology, which utilizes geospatial data on point sensors and internet optical fiber deployed in metropolitan statistical areas (MSAs) in the United States. We extend these data to identify the area that can be monitored by (1) considering the observed seismic noise data in target locations, (2) applying the model from Wilson et al. (2021) to understand the potential coverage area gains using optical fiber sensing, and (3) optimizing the selection of fiber segments to maximize coverage and minimize deployment costs. We implement our methodology in ArcGIS to assess the additional area that can be monitored for low-amplitude ground-motion events (i.e., magnitude >0.5) by utilizing internet fiber-optic cables in the 100 most populous MSAs in the United States. We find that the addition of internet fiber-based sensors in MSAs would increase the area monitored on average by over an order of magnitude from 1% to 12%, if the subset of fiber cable segments that maximize coverage and minimize deployment costs is chosen even if only 20% of all fibers are used.

58 GEOSCIENCES↗

A Deep Learning-Aided Workflow for Decoding the Stress Regime of Southern Nevada

The Rock Valley fault zone in southern Nevada has a notable history of seismic activity and is the site of a future direct comparison experiment of explosion and earthquake sources. This study aims to gain insight into regional tectonic processes by leveraging recent advances in seismic monitoring capabilities to elucidate the local stress regime. A crucial step in this investigation is the accurate determination of P-wave first-motion polarities, which play a vital role in resolving earthquake focal mechanisms of small earthquakes. Here, we deploy a deep learning-based method for automatic determination of first-motion polarities to vastly expand the polarity dataset beyond what has been reviewed by human analysts. By the integrating P-wave polarities with new measurements of S/P amplitude ratios, we obtain robust focal mechanism estimates for 1306 earthquakes with a local magnitude of 1 and above occurring between 2010 and 2023 in southern Nevada. We then use the focal mechanism catalog to examine the regional stress orientation, confirming an overall trans-tensional stress regime with smaller scale complexities illuminated by individual earthquake sequences. These findings demonstrate how detailed analyses of small earthquakes can provide fundamental information for understanding earthquake processes in the region and inform future experiments at the Nevada National Security Site.

58 GEOSCIENCES↗

The Value of Hyperparameter Optimization in Phase-Picking Neural Networks

The effectiveness of neural networks for picking seismic phase arrival times has been demonstrated through several case studies, and seismic monitoring programs are starting to adopt the technology into their workflows. However, published models were designed and trained using rather arbitrary choices of hyperparameters, limiting their performance. In this study, we use phase picks from both routine and template-matching analyses from multiple regions (Ridgecrest, California; Kilauea, Hawaii; Yellowstone, Wyoming–Montana–Idaho) to test a hyperparameter optimization scheme for phase-picking neural networks and to evaluate their performance. We show that a published model, namely PhaseNet (Zhu and Beroza, 2019), can be simplified and improved with reasonable effort and there are preferred choices of hyperparameters that increase the performance. We also show that models optimized based on the arrival times reported in routine event catalogs consistently perform well when picking arrival times of smaller events, which is crucial for many tasks from microseismicity to explosion monitoring.

58 GEOSCIENCES↗