SEARCH · Engineering Papers

Results for “stochastic transition”

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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Theoretical Study of the Single Electron Optical Stochastic Cooling Experiment at the IOTA Facility (Final Technical Report)

The objectives of this proposal were to explore fundamental aspects of radiation from individual and small groups of electrons in the optical stochastic cooling experiment conducted at Fermi National Laboratory's IOTA facility. Professor Swapan Chattopadhyay collaborated directly with the IOTA team on their initial experimental endeavors. Concurrently, a team from Berkeley, including Prof. Wurtele, Dr. Charman, and Dr. Greggory Penn from LBNL, and Prof. Chattapadhyay, pursued a theoretical exploration to identify potential subsequent experiments for IOTA. A wide range of possible experiments were considered, but they do not allow for tests for new quantum mechanical phenomena at the precision level of the current IOTA apparatus. However, the experimentally observed optical stochastic cooling of a single electron in IOTA validated the relevance of the classical electromagnetic field nature of the continuous electron-radiation interaction (as opposed to discrete instantaneous photon emissions and associated 'classical particle trajectory jumps') responsible for cooling and pointed towards innovative methods of experimental techniques to measure 'turn-by-turn' field correlations. Such measurements have the promise of bridging the transition from classical fields building up to discrete quantum transitions and enhance proper theoretical understanding.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1996

Uncertainty and Intelligence in Computational Stochastic Mechanics

Classical structural reliability assessment techniques are based on precise and crisp (sharp) definitions of failure and non-failure (survival) of a structure in meeting a set of strength, function and serviceability criteria. These definitions are provided in the form of performance functions and limit state equations. Thus, the criteria provide a dichotomous definition of what real physical situations represent, in the form of abrupt change from structural survival to failure. However, based on observing the failure and survival of real structures according to the serviceability and strength criteria, the transition from a survival state to a failure state and from serviceability criteria to strength criteria are continuous and gradual rather than crisp and abrupt. That is, an entire spectrum of damage or failure levels (grades) is observed during the transition to total collapse. In the process, serviceability criteria are gradually violated with monotonically increasing level of violation, and progressively lead into the strength criteria violation. Classical structural reliability methods correctly and adequately include the ambiguity sources of uncertainty (physical randomness, statistical and modeling uncertainty) by varying amounts. However, they are unable to adequately incorporate the presence of a damage spectrum, and do not consider in their mathematical framework any sources of uncertainty of the vagueness type. Vagueness can be attributed to sources of fuzziness, unclearness, indistinctiveness, sharplessness and grayness; whereas ambiguity can be attributed to nonspecificity, one-to-many relations, variety, generality, diversity and divergence. Using the nomenclature of structural reliability, vagueness and ambiguity can be accounted for in the form of realistic delineation of structural damage based on subjective judgment of engineers. For situations that require decisions under uncertainty with cost/benefit objectives, the risk of failure should depend on the underlying level of damage and the uncertainties associated with its definition. A mathematical model for structural reliability assessment that includes both ambiguity and vagueness types of uncertainty was suggested to result in the likelihood of failure over a damage spectrum. The resulting structural reliability estimates properly represent the continuous transition from serviceability to strength limit states over the ultimate time exposure of the structure. In this section, a structural reliability assessment method based on a fuzzy definition of failure is suggested to meet these practical needs. A failure definition can be developed to indicate the relationship between failure level and structural response. In this fuzzy model, a subjective index is introduced to represent all levels of damage (or failure). This index can be interpreted as either a measure of failure level or a measure of a degree of belief in the occurrence of some performance condition (e.g., failure). The index allows expressing the transition state between complete survival and complete failure for some structural response based on subjective evaluation and judgment.

Ayyub, Bilal M.↗

DOE OSTI2025

Changing-look Active Galactic Nuclei from the Dark Energy Spectroscopic Instrument. II. Statistical Properties from the First Data Release

We present the identification of changing-look active galactic nuclei (CL-AGNs) from the Dark Energy Spectroscopic Instrument First Data Release and Sloan Digital Sky Survey Data Release 16 at z ≤ 0.9. To confirm the CL-AGNs, we utilize spectral flux calibration assessment via an [O III ]-based calibration, pseudophotometry examination, and visual inspection. This rigorous selection process allows us to compile a statistical catalog of 561 CL-AGNs, encompassing 527 Hβ, 149 Hα, and 129 Mg II CL behaviors. In this sample, we find (1) a 283:278 ratio of turn-on to turn-off CL-AGNs. (2) The median Eddington ratio for CL-AGNs in the dim state is approximately λ Edd ∼ 0.01. (3) A strong correlation between the change in the luminosity of the broad emission lines (BELs) and variation in the continuum luminosity, with Mg II and Hβ displaying similar responses during CL phases. (4) The Baldwin–Phillips–Terlevich diagram for CL-AGNs shows no statistical difference from the general AGN catalog. (5) Five CL-AGNs are associated with asymmetrical mid-infrared flares, possibly linked to tidal disruption events. Given the large CL-AGN sample and the stochastic sampling of spectra, we propose that some CL phenomena are inherently due to typical AGN variability during low accretion rates, particularly for CL phenomenon only occurring on one BEL. Finally, we introduce a monotonically dimming CL phase for objects characterized by a gradual decline over decades in the light curve and the complete disappearance of entire BELs in faint spectra, indicative of a real transition in the accretion disk.

accretion↗

DOE OSTI2025

Changing-look Active Galactic Nuclei from the Dark Energy Spectroscopic Instrument. II. Statistical Properties from the First Data Release

We present the identification of changing-look active galactic nuclei (CL-AGNs) from the Dark Energy Spectroscopic Instrument First Data Release and Sloan Digital Sky Survey Data Release 16 at z≤ 0.9. To confirm the CL-AGNs, we utilize spectral flux calibration assessment via an [O iii]-based calibration, pseudophotometry examination, and visual inspection. This rigorous selection process allows us to compile a statistical catalog of 561 CL-AGNs, encompassing 527 Hβ, 149 Hα, and 129 Mg ii CL behaviors. In this sample, we find (1) a 283:278 ratio of turn-on to turn-off CL-AGNs. (2) The median Eddington ratio for CL-AGNs in the dim state is approximately λ$_{Edd}$ ∼ 0.01. (3) A strong correlation between the change in the luminosity of the broad emission lines (BELs) and variation in the continuum luminosity, with Mg ii and Hβ displaying similar responses during CL phases. (4) The Baldwin–Phillips–Terlevich diagram for CL-AGNs shows no statistical difference from the general AGN catalog. (5) Five CL-AGNs are associated with asymmetrical mid-infrared flares, possibly linked to tidal disruption events. Given the large CL-AGN sample and the stochastic sampling of spectra, we propose that some CL phenomena are inherently due to typical AGN variability during low accretion rates, particularly for CL phenomenon only occurring on one BEL. Finally, we introduce a monotonically dimming CL phase for objects characterized by a gradual decline over decades in the light curve and the complete disappearance of entire BELs in faint spectra, indicative of a real transition in the accretion disk.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Exponential Escape Rate of Filamentary Incubation in Mott Spiking Neurons

Mott materials such as vanadium oxides, when subject to a strong applied voltage, present an inhomogeneous insulator-to-metal transition with formation of metallic filaments within the insulating bulk. This property is enabling the development of compact and power-efficient neuromorphic devices known as Mott neurons. However, the nature of the transition has not been fully understood yet, as it may be attributed to different effects, including Joule self-heating and hot-carrier injection. Moreover, the experimental determination of the threshold voltage needed to induce the transition has proven to be challenging, as the transition becomes increasingly unpredictable when the threshold is approached. The physical understanding of these issues would not only deepen our understanding of Mott insulators, but would also be an important step toward the realization of neuromorphic devices based on such materials. In this work we use numerical simulations based on the Mott resistor network model to study the nature of the filament incubation and formation process. Here, we show that both electronic and thermal effects, in the form of current density focusing and Joule self-heating, respectively, contribute to the filamentary incubation and growth. Remarkably, we find that the percolation of the metallic filaments near the threshold is intrinsically stochastic, qualitatively similar to the familiar Arrhenius activated behavior and to the stochastic firing of biological neurons. More precisely, we characterize the filament percolation as a Poisson point process, which has the same probability distribution as mathematical models of neuronal firing with an exponential escape rate. Finally, we support the numerical simulation results by performing experiments in VO 2 that are in agreement with the exponential escape rate behavior. Thus, we establish a functionality of Mott insulators that opens a path toward implementing neuromorphic hardware with quantum materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

AI-Assisted Conceptual Development of a Pre-Geometric Cosmological Model - An Exercise in AI-Assisted Conceptual Framework Generation, Paper III: Cosmological Structure and Predictions

This paper develops the cosmological consequences of the replication-driven cosmogenesis framework introduced in Paper I and the emergent geometric structure established in Paper II. After the replication epoch freezes out, the coherent sector occupies a finite spectral band and contains a population of excited states. The relaxation of these excited coherent configurations does not produce coherent radiation; instead, all released energy flows into the incoherent substrate, where the randomizer acts as a rapid phase-scrambling mechanism. This process generates an effectively thermal radiation bath, providing a natural reheating mechanism that requires neither inflaton oscillations nor scalar-field potentials, and can be contrasted with standard scenarios of nonperturbative reheating dynamics. Subsequent symmetry-breaking transitions in the coherent vacuum inject additional radiation, yielding a multi-stage thermal history with well-defined energy transfers. We derive the effective equations of state for each component—the cosmological vacuum, the coherent vacuum, and the radiation bath—and show how their interplay produces an FRW-like expansion. The discrete sequence of coherent-state relaxations imprints a distinctive multi-peaked stochastic gravitational-wave background, whose spectral structure reflects the underlying hierarchy of coherent frequencies. Potential observational signatures in the LISA and mid-band frequency ranges are highlighted, providing concrete avenues to test this replication-based cosmological framework in the context of standard cosmological gravitational-wave backgrounds and LISA-oriented forecasts.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Co-Location of Cellulosic Bioethanol and Alcohol-to-Jet (ATJ) Production Facilities for Targeted Scale-Up of Sustainable Aviation Fuel (SAF) Production

Achieving aerospace industry net-zero emissions by 2050 requires rapid scaling of sustainable aviation fuel (SAF) production. Leveraging existing infrastructure, proven technologies like Alcohol-to-Jet (ATJ), and low carbon intensity (CI) feedstocks (e.g., switchgrass and miscanthus) can support this transition and help achieve near-term emissions reduction targets. This study evaluates the implications of lignocellulosic ethanol biorefinery siting and integration with petroleum refineries to produce SAF across 1000 sites randomly sampled from areas suitable for perennial grasses in the U.S. rainfed region. To better understand the logistics of material transport and handoffs, we integrated models of biomass harvest, transport, ethanol, and ATJ production in a stochastic framework based on Monte Carlo simulations to characterize SAF minimum selling price (MSP) and carbon intensity (CI), considering site-specific parameters (e.g., feedstock production, transportation, taxes, incentives). The results indicate trade-offs between MSP and CI across locations, with median MSP ranging from 7.9 to 12.8 USD·gal −1 and CI from −9.7 to 39.4 gCO 2 e·MJ −1 . Despite high estimated decarbonization costs (580 USD·tonCO 2 e −1 ), our results indicate that site-specific deployment of ATJ with low-CI feedstocks can improve sustainability outcomes. The framework provides a systematic approach to assess cost and sustainability trade-offs across locations, considering the end-to-end supply chain and supporting an informed investment in SAF production.

09 BIOMASS FUELS↗

NASA NTRS1991

Generalized covariance analysis for partially autonomous deep space missions

A new covariance analysis method is presented that is suitable for the evaluation of multiple impulsive controllers acting on some stochastic process x. The method accommodates batch and sequential estimators with equal ease and accounts for time-delay effects in a natural manner. The formalism is developed in terms of a generalized state vector that is formed from the system state vector x, augmented by various fixed epoch estimates, and a data vector formed from discrete time observations of the system. Recursions are developed for time transition, measurement incorporation, and impulsive control updating of the generalized covariance matrix. Means of limiting the dimensional growth of the generalized state vector via the processes of estimator epoch adjustment and measurement vector deflation are described and the application of numerically stable matrix factorization methods to the generalized covariance recursions is outlined. The method is applied to the Magellan spacecraft to demonstrate the capability of ground-based optimal estimation and control of gyro/star scanner misalignment.

Boone, Jack N.↗

NASA NTRS2004

Computation Molecular Kinetics Model of HZE Induced Cell Cycle Arrest

Cell culture models play an important role in understanding the biological effectiveness of space radiation. High energy and charge (HZE) ions produce prolonged cell cycle arrests at the G1/S and G2/M transition points in the cell cycle. A detailed description of these phenomena is needed to integrate knowledge of the expression of DNA damage in surviving cells, including the determination of relative effectiveness factors between different types of radiation that produce differential types of DNA damage and arrest durations. We have developed a hierarchical kinetics model that tracks the distribution of cells in various cell phase compartments (early G1, late G1, S, G2, and M), however with transition rates that are controlled by rate-limiting steps in the kinetics of cyclin-cdk's interactions with their families of transcription factors and inhibitor molecules. The coupling of damaged DNA molecules to the downstream cyclin-cdk inhibitors is achieved through a description of the DNA-PK and ATM signaling pathways. For HZE irradiations we describe preliminary results, which introduce simulation of the stochastic nature of the number of direct particle traversals per cell in the modulation of cyclin-cdk and cell cycle population kinetics. Comparison of the model to data for fibroblast cells irradiated photons or HZE ions are described.

Cucinotta, Francis A.↗

NASA NTRS2009

Alternative Representations of Convective Processes in the NASA GEOS-5 AGCM

The gap in explicit resolution of phenomena between global climate models and cloud resolving models is shrinking at a steady pace with global integrations of several km in resolution now practical for at least some time scales. In moving toward finer resolution the long standing problem of convective parameterization is being examined along with the assumption of convective quasi-equilibrium and how this can be reconciled with the stochastic and intermittent nature of convection. In this context we examine the nature of parameterized convection in the NASA Goddard Earth Observing System (GEOS-5) Atmospheric General Circulation Model. Our analysis uses both coarse (2.5 degree) and fine scale 0.25 degree spatial resolution integrations. Two basic formulations are compared: the default option is the Relaxed Arakawa-Schubert (RAS) scheme which invokes a sequence of linearly entraining plumes and quasi-equilibrium closure. An optional modification of this method (the "Stochastic Tokioka" constraint) places a random lower limit on plume entrainment. An alternative representation is the Kain-Fritsch parameterization which was originally developed for mesoscale numerical modeling strategies. Here entrainment is determined by a crude buoyancy sorting approach that allows the plume spectrum to be more responsive to ambient vertical stratification of moisture. Diagnostics of the model behavior are referenced to recent observational evidence of continuous phase transition behavior. In particular we examine the relationship between column water vapor and probablility of convective presence and intensity. Sensitivity of the statistics of convective behavior to parcel mixing/entrainment formulations and parcel initial thermodynamics are considered. Observational statistics from A-Train and TRMM sensors provide validation of the model integrations.

Robertson, Franklin↗

DOE OSTI2023

Measurement-induced entanglement and teleportation on a noisy quantum processor

Measurement has a special role in quantum theory: by collapsing the wavefunction, it can enable phenomena such as teleportation and thereby alter the ‘arrow of time’ that constrains unitary evolution. When integrated in many-body dynamics, measurements can lead to emergent patterns of quantum information in space–time that go beyond the established paradigms for characterizing phases, either in or out of equilibrium. For present-day noisy intermediate-scale quantum (NISQ) processors, the experimental realization of such physics can be problematic because of hardware limitations and the stochastic nature of quantum measurement. Here we address these experimental challenges and study measurement-induced quantum information phases on up to 70 superconducting qubits. By leveraging the interchangeability of space and time, we use a duality mapping to avoid mid-circuit measurement and access different manifestations of the underlying phases, from entanglement scaling to measurement-induced teleportation. We obtain finite-sized signatures of a phase transition with a decoding protocol that correlates the experimental measurement with classical simulation data. The phases display remarkably different sensitivity to noise, and we use this disparity to turn an inherent hardware limitation into a useful diagnostic. Our work demonstrates an approach to realizing measurement-induced physics at scales that are at the limits of current NISQ processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Impact of confined geometries on hopping and trapping of motile bacteria in porous media

We use a random walk particle-tracking (RWPT) approach to elucidate the impact of porous media confinement and cell-cell interactions on bacterial transport. The model employs stochastic alternating motility states consisting of hopping movement and trapping reorientation. The stochastic motility patterns are defined based on direct visualization of individual trajectory data. We validate our model against experimental data, at single-cell resolution, of bacterial E. coli motion in three-dimensional confined porous media. Results show that the model is able to efficiently simulate the spreading dynamics of motile bacteria as it captures the impact of cell-cell interaction and pore confinement, which marks the transition to a late-time subdiffusive regime. Furthermore, the model is able to qualitatively reproduce the observed directional persistence. Lastly, our RWPT model constitutes a meshless simple method which is easy to implement and does not invoke ad hoc assumptions but represents the basis for a multiscale approach to the study of bacterial dispersal in porous systems.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1987

Implementation of variable time step stochastic dynamics for electronically inelastic gas-surface collisions

A variable time step algorithm has been implemented for solving the stochastic equations of motion for gas-surface collisions. It has been tested for a simple model of electronically inelastic collisions with an insulator surface in which the phonon manifold acts as a heat bath and electronic states are localized. In addition to reproducing the accurate nuclear dynamics of the surface atoms, numerical calculations have shown the algorithm to yield accurate ensemble averages of physical observables such as electronic transition probabilities and total energy loss of the gas atom to the surface. This new algorithm offers a gain in efficieny of up to an order of magnitude compared to fixed time step integration.

Garrett, Bruce C.↗

DOE OSTI2022

Spectroscopy of particle couplings with gravitational waves

We discuss the possibility to measure particle couplings with stochastic gravitational wave backgrounds (SGWBs). Under certain circumstances a sequence of peaks of different amplitude and frequency — a stairway — emerges in a SGWB spectrum, with each peak probing a different coupling. The detection of such signature opens the possibility to reconstruct couplings (spectroscopy) of particle species involved in high energy phenomena generating SGWBs. Stairwaylike signatures may arise in causally produced backgrounds in the early Universe, e.g., from preheating or first order phase transitions. As a proof of principle we study a preheating scenario with an inflaton $\phi$ coupled to multiple daughter fields {$χ_j$} with different coupling strengths. As a clear stairway signature is imprinted in the SGWB spectrum, we reconstruct the relevant couplings with various detectors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Computational Math Problems for a Clean Energy Future

Cutting edge computational mathematics are ubiquitous in renewable energy research. Problems in resilient and reliable electric grid operations, infrastructure planning, wind farm yaw control, and more demand sophisticated and scalable computational tools that enable the transition of renewable energy technologies from proof of concept to deployment into our energy system. The mission of the Computational Science Center at NREL is to lead the lab's efforts to solve energy challenges using high-performance computing (HPC), computational science, applied mathematics, scientific data management, visualization, and informatics. In this poster, we provide a short overview of three areas of computational mathematics research at NREL: wind power scenario generation for stochastic grid operations and infrastructure planning, improved rational function approximations for electromagnetic transients codes, and wind farm yaw control using a combination of the Alternating Direction Method of Multipliers (ADMM) and reinforcement learning (RL). Increasing penetrations of renewable energy into power grids motivate the investigation of new approaches to characterizing uncertainty for five-minute economic dispatch problems. Similarly, as the penetration of distributed energy resources on power grids increases, it becomes important to revisit our methods of modelling transient phenomena, i.e. electromagnetic transients programs. Finally, the combination of ADMM and RL for wind farm yaw control presented here can potentially increase the efficiency of the deployed distributed controllers by orders of magnitude.

ADMM↗

DOE OSTI2025

Fuel reid vapor pressure level and ethanol content on stochastic preignition, effects at steady and unsteady engine operation

The present work investigates relations between fuel Reid vapor pressure (RVP) and biofuel (ethanol) content on stochastic preignition (SPI) at both sustained steady-state engine operation and following load transients. This work stems from in-field observations that automotive original equipment manufacturers have observed consistent seasonal increases in United States customer drivability complaints and warranty claims during September and October where SPI is suspected to be responsible. The seasonal timing of these events coincides with the United States seasonal fuel property changeover initiating on September 15 each year, where fuel RVP increases. To explore potential linkage between fuel RVP and SPI the present study employs engine SPI experiments coupled with laboratory spray measurements of fuels with RVPs of 8, 12, and 16 psi in both E10 (10% ethanol) and E25 (25% ethanol) fuels. Engine results are partitioned into fuel RVP and ethanol content effects on SPI in steady-state, sustained high-load engine operation and unsteady-state low- to high-load transitions, where off-engine spray vessel patternation and tip penetration results help to elucidate the observed fuel effects on SPI. A boosted direct-injected, spark-ignition engine was fueled with three market relevant E10 and E25 fuels with RVPs of 8, 12, and 16 to characterize the interplay between winter fuels and abnormal combustion behavior. The steady-state work shows that for high-load, steady-state engine operation, SPI is directly linked to fuel retention, which was found to be dependent on fuel distillation. The unsteady-state engine operation work shows that following low-to high-load transitions, SPI can occur from a memory of fuel property effects at low-load operation. Specifically, the fuel RVP effect on fuel spray collapse at low loads was found to correlate with SPI with a more than 95% confidence interval following low- to high-engine-load transitions. Results suggest that fuel-wall impingement at low-load operation could carry over into high-load transitions and generate SPI events following low- to high-load transitions.

09 BIOMASS FUELS↗

DOE OSTI2021

Gravitational waves from breaking of an extra U (1) in SO (10) grand unification

Abstract In a class of gauged $U(1)$ extended Standard Models (SMs), the breaking of the $U(1)$ symmetry is not only a source for Majorana masses of right-handed (RH) neutrinos crucial for the seesaw mechanism, but also a source of stochastic gravitational wave (GW) background. Such $U(1)$ extended models are well-motivated from the viewpoint of grand unification. In this paper, we discuss a successful ultraviolet completion of a $U(1)$ extended SM by an $SO(10)$ grand unified model through an intermediate step of $SU(5) \times U(1)$ unification. With a parameter set that is compatible with the $SO(10)$ grand unification, we find that a first-order phase transition associated with the $U(1)$ symmetry breaking can be strong enough to generate GWs with a detectable size of amplitude. We also find that the resultant GW amplitude reduces and its peak frequency becomes higher as the RH neutrino masses increase.

Physics↗

DOE OSTI2021

The effects of electrically exploding gold bridgewires into inert and explosive powder beds

Abstract The particle velocity created in beds of both low-density inert sugar and explosive PETN as a function of distance from an exploding bridgewire was measured using optical velocimetry and a silvered PMMA window. As expected, more violent bridge-bursts (from a greater-stored-energy capacitive discharge unit) resulted in greater particle velocities and a better supported compaction wave in sugar. In all cases, ramp waves, not shocks, were observed in the inert sugar. Large window velocities were observed for very powerful bursts (up to 270 m/s), but bursts required for stochastic detonator operation conditions resulted in sugar/PMMA window velocities of only 8–10 m/s 0.85 mm from the bridge location. In contrast, after a distance of only 0.65 mm, a building shock wave was observed in PETN under both threshold and reliable firing conditions. Subsequently a hot-spot-driven shock-to-detonation (SDT) process was observed prior to full detonation. The measured buildup process accounts for $$\approx $$ ≈ 66% of the so-called excess transit time (ETT) between the observed and theoretical total function time for the particular exploding-bridge-wire (EBW) detonator studied. The remainder must occur in the powerful output pellet region. In contrast to a common understanding, the ETT is found to be a weak function of the discharge energy. Thus, the operation of the detonator after a bridge-burst energy-to-powder reaction transition process is found to be hot-spot-driven SDT in both the low- and high-density pellets.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗