Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive analytics”

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,171 records · Page 65

A Fabrication Strategy for Reconfigurable Millimeter–Scale Metamaterials

Rather than depending on material composition to primarily dictate performance metrics, metamaterials can leverage geometry to achieve specific properties of interest. For example, reconfigurable metamaterials have enabled programmable shape transformations, tunable mechanical properties, and energy absorption. While several methods exist to fabricate such structures, they often place severe restrictions on manufacturing materials, or require significant manual assembly. Moreover, these arrays are typically composed of unit cells that are either macro–scale or micro–scale in dimension. Here, the fabrication gap is bridged, and laminate manufacturing is used to develop a method for designing reconfigurable metamaterials at the millimeter–scale, that is compatible with a wide range of materials, and that requires minimal manual assembly. In addition to showing the versatility of this fabrication method, how the use of laminate manufacturing affects the behavior of these multi–component arrays is also characterized. To this end, a numerical model that captures the deformations exhibited by the structures is developed, and an analytic model that predicts the strain of the structure under compressive stress is built. Overall, this approach can be leveraged to develop millimeter–scale metamaterials for applications that require reconfigurable materials, such as in the design of tunable acoustics, photonic waveguides, and electromagnetic devices.

36 MATERIALS SCIENCE↗

Algebraic expressions for estimating the impact depths of a surface barrier over a homogeneous soil

Engineered surface barriers are used to isolate subsurface contaminants for effective long-term containment of municipal solid waste, other nonhazardous solid and liquid waste, hazardous and toxic wastes, and radioactive waste. The impact depths of a surface barrier are affected by the pre-barrier recharge rate and the properties of the soil beneath the barrier. In this paper, the pore-size-specific (PSS) water velocity is defined and an algebraic expression of PSS velocity is derived based on the stream tube concept and the Brooks and Corey hydraulic retention model. Algebraic expressions are developed to estimate drainage velocities and barrier impact depths after the emplacement of a surface barrier. Four impact depth terms are used to convey the protective effect: drainage front, average drainage, the location with 50% impact, and drainage tail. The drainage front depth is the deepest point at which the barrier has a detectable impact at a specific time (also called the near zero-impact depth). The average-impact depth is the depth at which average drainage occurs. At the 50% impact depth, the water flux rate is reduced by half because of the surface barrier. Lastly, the drainage tail depth (also called the full-impact depth) is the deepest depth at which the water conditions above it are in equilibrium with the barrier. The algebraic expressions show that the average-impact depth is no more than 1/3 of near zero-impact depth, while the 50% impact depth is slightly larger than 1/2 of the near zero-impact depth. The full-impact depth, depending on the final recharge rate from the surface barrier, is usually much smaller than the other impact depths. These differences lead to a very large transition zone beneath a surface barrier. Numerical simulations were conducted to replicate the same conditions. The numerical results corroborated the analytical models by predicting very similar water content profiles and near zero-, average-, 50%, and full-impact depths. The algebraic expressions provided in this paper are useful for quickly identifying sites where the depth of the existing contaminants could be beyond the protection of a surface barrier.

54 ENVIRONMENTAL SCIENCES↗

Floating platform effects on power generation in spar and semisubmersible wind turbines

The design and financing of commercial-scale floating offshore wind projects require a better understanding of how power generation differs between newer floating turbines and well-established fixed-bottom turbines. In floating turbines, platform mobility causes additional rotor motion that can change the time-averaged power generation. In this work, OpenFAST simulations examine the power generated by the National Renewable Energy Laboratory's 5-MW reference turbine mounted on the OC3-UMaine spar and OC4-DeepCWind semisubmersible floating platforms, subjected to extreme irregular waves and below-rated turbulent inflow wind from large-eddy simulations of a neutral atmospheric boundary layer. For these below-rated conditions, average power generation in floating turbines is most affected by two types of turbine displacements: an average rotor pitch angle that reduces power, caused by platform pitch; and rotor motion upwind-downwind that increases power, caused by platform surge and pitch. The relative balance between these two effects determines whether a floating platform causes power gains or losses compared to a fixed-bottom turbine; for example, the spar creates modest (3.1%–4.5%) power gains, whereas the semisubmersible creates insignificant (0.1%–0.2%) power gains for the simulated conditions. Furthermore, platform surge and pitch motions must be analyzed concurrently to fully capture power generation in floating turbines, which is not yet universal practice. Finally, a simple analytical model for predicting average power in floating turbines under below-rated wind speeds is proposed, incorporating effects from both the time-averaged pitch displacement and the dynamic upwind-downwind displacements.

17 WIND ENERGY↗

Stability study of aqueous foams under high-temperature and high-pressure conditions relevant to Enhanced Geothermal Systems (EGS)

In Enhanced Geothermal System (EGS) hydraulic fracturing is carried out by injecting cold water into deep Hot Dry Rocks (HDR) under carefully controlled conditions to create new or reopen existing fractures. Water-based fracturing fluids demonstrate some challenges including immense quantity of water usage, water sensitivity of the formations, water blocking, and lack of proppant carrying capacity and transportation. Thus, an alternative is to use foam-based fracturing fluid which offers potential advantage over conventional water-based fracturing fluid such as minimum water usage, reduced wellbore damage, high proppant carrying capacity, and less environmental damage. However, foams are complex mixture of gaseous phase and liquid phase which are thermodynamically unstable at downhole conditions, and their stability can decrease over time due to liquid drainage, bubble coarsening, and coalescence. This paper shows laboratory experiments executed to study foam stability at high temperature (200°C) and high pressure (6.9MPa) conditions which simulates the geothermal environment. Foam stability was characterized by half-life of foam, which is defined as the time taken by the foam to decreases by 50% of its original height due to drainage. In this paper, two types of gaseous phases, nitrogen (N 2 ) and carbon dioxide (CO 2 ) were investigated. Also, based on successful practice of foam-based fracturing fluid in oil and gas industries, four surfactants, including Alpha olefin sulfonate (AOS), Sodium dodecyl sulfonate (SDS), Tergitol™ (NP – 40), and Cetyltrimethylammonium chloride (CTAC) at optimum concentration of 1wt.% were tested for best stability performance. In addition, different stabilizing agents including guar gum, bentonite clay, crosslinker, silicon dioxide nanoparticles (SiO 2 ), graphene oxide (GO) were also studied. Experimental results showed that N 2 foams were more stable than CO 2 foams. It was observed that foam half-life decreased with the increase in temperature. Among all the surfactants, AOS foams showed the most promising thermal stability at high temperatures. Moreover, with the addition of stabilizing agents, foam's half-life was enhanced. Stabilizing agents such as crosslinker and GO dispersion showed the most stable foams with half-life recorded at 20min and 17min, respectively, at 200°C and 6.9MPa. Finally, pressure also showed a positive effect on foam stability; with increased pressure, foam half-life was increased. Based on the experimental data, analytical models for the effect of temperature and pressure were developed, considering foam degradation is a first-order kinetic reaction that linearly depends on the foam drainage mechanism. The effect of temperature on foam half-life was studied as an exponential decay model. In this model, foam half-life is a function of drainage rate constant (D A ) and activation energy (E a ) of the foam system. The effect of pressure on foam half-life was found to obey a power-law model where an increase in pressure showed an increase in foam half-life. Furthermore, a linear relation was studied for the effect of pressure on foam activation energy and drainage rate. Then the combined effects of temperature and pressure were studied, which yielded an analytical model to predict the foam stabilities in terms of half-life for different foam compositions. In conclusion, this research indicates that with an appropriate selection of surfactants and stabilizing agents, it is possible to obtain stable foams, which could replace conventional water fracturing fluid under EGS conditions.

15 GEOTHERMAL ENERGY↗

Generalized Analysis of Dynamic Flow Fouling on Heat Transfer Surfaces

Several energy applications involve flow of untreated heat transfer fluids that contain trace amounts of mineral salts, such as calcium carbonate and calcium sulfate, which show reverse solubility with temperature. The gradual accretion of these salts onto the heat transfer surfaces causes fouling that degrades thermal performance over time. While fouling of bare surfaces has been extensively studied, fouling of nonwetting superhydrophobic or liquid-infused surfaces is less understood, especially under dynamic flow conditions. This study presents a systematic experimental analysis to elucidate dynamic flow fouling of calcium carbonate and calcium sulfate on smooth, superhydrophobic and liquid-infused surfaces in a generalized manner. Here, copper tubes with modified inner surface wettability are fabricated and subjected to forced convection at different flow rates and fouling salt concentrations. Fouling characteristics of the two salts on the different surfaces are examined through microstructural analysis, and quantified in terms of asymptotic fouling resistance that measures the total extent of fouling, and a threshold time, which correlates to the onset of significant fouling. Using Taguchi design of experiments, the study reports, for the first time, closed-form analytical relationships for the asymptotic fouling thermal resistance and threshold time on the Reynolds number, dimensionless concentration of fouling agent and dimensionless infused liquid viscosity that represents the different surface types in a unified manner. It is shown that the analytical model accurately predicts the fouling characteristics in up to about 97% of the data, based on which contour maps of optimum surface designs for minimizing fouling resistance are presented. The results provide information on designing surfaces for reduced fouling or for estimating the fouling characteristics for a given surface in applications.

36 MATERIALS SCIENCE↗

On the spectral behaviour of the turbulence-driven power fluctuations of horizontal-axis turbines

In this article we consider the spectral behaviour of turbulence-driven power fluctuations for a single horizontal-axis turbine. To this end, a small-scale instrumented axial-flow hydrokinetic turbine model ( diameter=0.724 m ) is deployed in the long water flume situated in the laboratory facilities of IFREMER in Boulogne-sur-Mer, France, and synchronous measurements of the upstream velocity and the rotor are collected for different tip-speed ratios. The study confirms previous findings suggesting that the power spectra follow the velocity spectra behaviour in the large scales region and a steeper power law slope behaviour ( -11/3 ) over the inertial frequency sub-range. However, we show that both the amplitude of the power spectra and low-pass filtering effect over the inertial sub-range also depend on the rotor aero/hydrodynamics (e.g. dC L /dα) and the approaching flow deceleration and not solely on the rotational effects. In addition, we present a novel semi-analytical model to predict the dominant blade-passing frequency harmonics in the high-frequency regime using the rotationally sampled spectra technique. For all calculations, the distortion of incoming turbulence is taken into account.

42 ENGINEERING↗

A vortex sheet based analytical model of the curled wake behind yawed wind turbines

Motivated by the need for compact descriptions of the evolution of non-classical wakes behind yawed wind turbines, we develop an analytical model to predict the shape of curled wakes. Interest in such modelling arises due to the potential of wake steering as a strategy for mitigating power reduction and unsteady loading of downstream turbines in wind farms. We first estimate the distribution of the shed vorticity at the wake edge due to both yaw offset and rotating blades. By considering the wake edge as an ideally thin vortex sheet, we describe its evolution in time moving with the flow. Vortex sheet equations are solved using a power series expansion method, and an approximate solution for the wake shape is obtained. The vortex sheet time evolution is then mapped into a spatial evolution by using a convection velocity. Apart from the wake shape, the lateral deflection of the wake including ground effects is modelled. Our results show that there exists a universal solution for the shape of curled wakes if suitable dimensionless variables are employed. For the case of turbulent boundary layer inflow, the decay of vortex sheet circulation due to turbulent diffusion is included. Finally, we modify the Gaussian wake model by incorporating the predicted shape and deflection of the curled wake, so that we can calculate the wake profiles behind yawed turbines. Model predictions are validated against large-eddy simulations and laboratory experiments for turbines with various operating conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Unified Analytical Method to Quantify Three Types of Fast Frequency Response from Inverter-Based Resources

With more inverter-based resources (IBRs), our power systems have lower frequency nadirs following N-1 contingencies, and undesired under-frequency load shedding (UFLS) can occur. To address this challenge, IBRs can be programmed to provide at least three types of fast frequency response (FFR), e.g., step response, proportional response (P/f droop response), and derivative response (synthetic inertia). However, these heterogeneous FFR challenge the study of power system frequency dynamics. Thus, this paper develops an analytical frequency nadir prediction method that allows for the consideration of all three potential forms of FFR provided by IBRs. The proposed method provides fast and accurate frequency nadir estimation after N-1 generation tripping contingencies. Our method is grounded on the closed-form solution for the frequency nadir, which is solved from the second-order system frequency response model considering the governor dynamics and three types of FFR. The simulation results in the IEEE 39-bus system with different types of FFR demonstrate that the proposed method provides an accurate and fast prediction of the frequency nadir under various disturbances.

fast frequency response↗

Black hole-neutron star binaries near neutron star disruption limit in the mass regime of event GW230529

In May 2023, the LIGO-Virgo-KAGRA (LVK) Collaboration reported the likely black hole-neutron star (BHNS) merger GW230529_181500. The signal was observed with high significance in only one detector, limiting constraints on the black hole spin and motivating our study of disruption in this mass regime. That event is expected to be the merger of a 2.5–4.5 M ⊙ primary with a secondary compact object of mass between 1.2–2.0 M ⊙ . This makes it the first BHNS merger with a significant potential for the production of electromagnetic (EM) counterparts, and provides further evidence for compact objects existing within the suspected lower mass gap. To produce post-merger EM transients, the component of the black hole spin aligned with the orbital angular momentum must be sufficiently high, allowing the neutron star to be tidally disrupted. The disrupting BHNS binary may then eject a few percent of a solar mass of matter, leading to an observable kilonova driven by radioactive decays in ejecta, and/or a compact-binary gamma-ray burst (cbGRB) resulting from the formation of an accretion disk and relativistic jet. Determining which mergers lead to disruption of the neutron star is necessary to predict the prevalence of EM signals from BHNS mergers, yet most BHNS simulations so far have been performed far from the minimum spin required for tidal disruption. Here, we use the Spectral Einstein Code to explore the behavior of BHNS mergers in a mass range consistent with GW230529_181500 close to that critical spin, and compare our results against the mass remnant model currently used by the LVK Collaboration to predict the probability of tidal disruption. Furthermore, our numerical results reveal the emergence of non-zero accretion disks even below the predicted NS disruption limit, of low mass but capable of powering cbGRBs. Our results also demonstrate that the remnant mass model underpredicts the disk mass for the DD2 equation of state, while they are within expected modeling errors for SFHo. The disruption limit itself, however, is not found to significantly differ from the predictions of the analytical model, unless remnant masses M rem ≲ 0.001 M ⊙ prove interesting observationally. In all of our simulations, any kilonova signal would be dim and most likely dominated by post-merger disk outflows.

GW230529_181500↗

Solidification dynamics in metal additive manufacturing: analysis of model assumptions

Solidification characteristics in the meltpool drive the process-microstructure relationship which helps dictate the material properties of as-built parts in additive manufacturing; therefore, being able to accurately and quickly simulate the size, shape, and solidification characteristics in the melt pool is of great interest to the field. In this study, we investigate various important physical phenomena (dynamic material properties, fluid-flow, radiation and vaporization) which can either be included or neglected in a continuum finite volume model (FVM) and their effect on the solidification conditions. Additionally, since the simplest form of such a model (conduction only) has an analytic solution which is much faster, its viability is also considered. Since the inclusion of some of these physical phenomena will inherently change the net energy input as well as the amount of energy needed to achieve melting of a control volume, each set of included phenomena had an effective absorption efficiency which was calibrated to closely match the dimensions of the melt pool to that of the ground truth data. The ground truth data for this study was defined to be the output of the FVM which included all the physical phenomena (OF). This study then goes on to compare the effects on solidification conditions each of these calibrated models has. It was found that most of the change in solidification conditions comes from the inclusion of latent heat. A posterior correlation factor (PCF) is then introduced to enable an analytic model to predict similar solidification conditions to OF model.

36 MATERIALS SCIENCE↗

The effect of curvature induced broken potential vorticity conservation on drift wave turbulences

Abstract The correlation theory of turbulence suppression (Zhang and Mahajan 1993 Phys. Fluids B5 2000) by velocity shear was constructed by invoking the ansatz of potential vorticity conservation (PVC) that holds for relatively simple (slab) models of drift wave turbulence. It is, therefore, surprising that a detailed modern simulation of the H mode pedestal, using the gyrokinetic code GENE (Hatch et al 2018 Plasma Phys. Control. Fusion 60 084003), found ‘striking agreement’ with the predictions of the analytic model. To understand the reasons for this remarkable agreement, an extended theory that contains finite (magnetic) curvature, and which does not conserve potential vorticity, is developed and ‘solved’ by calculating an inhomogeneous Green function reflecting the fact that the new system has a potential vorticity source. It is, then, demonstrated that the effect of the broken PVC is insignificant for the normal operation parameters in tokamaks; the correction due to curvature is at the order of 2 L 0 / R , where L 0 is the scale length of the local gradient and R is the major radius. The excellent agreement between simulation and slab model 1993 theory is therefore not accidental; the latter can be applied with confidence to the tokamak pedestal.

Physics↗

Photon phase-space dynamics in a plasma wakefield accelerator

Frequency up-shifting of laser light in a beam-driven plasma wakefield has the potential to provide high-intensity sources of short wavelength radiation. Simulations have demonstrated that a laser pulse can undergo large frequency shifts, limited only by the drive beam energy, when the plasma density is tailored to match the accelerating phase of the wake to the group velocity of the pulse. Here, we study the dynamical evolution of photons in the 1D1P phase-space vicinity of the plasma wake-phase matching condition. Numerical calculations using a photon kinetic model are validated by direct comparison with 1D and quasi-3D particle-in-cell simulations. These calculations form the basis of a linear theory of the photon dynamics which reveals several important results, including scalings for the properties of the witness pulse and a self-similar solution for the photon phase-space dynamics. One prediction of the analytic theory is that the pulse can be compressed indefinitely, though the amount of compression would ultimately be limited by practical constraints. These results suggest that photon acceleration can provide a novel source of sub-femtosecond, short wavelength radiation.

XUV generation↗

Comparison between MAST-U conventional and Super-X configurations through SOLPS-ITER modelling

Abstract MAST-U has recently started operating with a Super-X divertor, designed to increase total flux expansion and neutral trapping, both predicted through simple analytic models and SOLPS calculations to reduce the plasma and impurity density detachment thresholds. In this study, utilising the SOLPS-ITER code, we are quantifying the possible gain allowed by the MAST-U Super-X and neutral baffling geometry, in terms of access to detachment. We show that a significant reduction of the upstream density detachment threshold (up to a factor 1.6) could be achieved in MAST-U, for the Super-X, as opposed to conventional divertor geometry, mainly through an increased total flux expansion, neutral trapping being found very similar between the different configurations. We also show that variations of the strike-point angle are complex to interpret in such a tightly baffled geometry, and that a case in which the target normal points more towards the separatrix does not necessarily imply a lower detachment threshold. As in previous calculations for TCV, we quantify the role of neutral effects through developing and applying a quantitative definition of neutral trapping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Infernal instabilities in negative-triangularity plasmas with negative central shear

Abstract A systematic numerical investigation is carried out to understand magnetohydrodynamic stability of the ideal infernal-kink instability in tokamak plasmas with both negative triangularity (neg-D) shaping and negative central shear for the equilibrium safety factor profile. The latter is motivated by the desire to form the internal transport barrier in the neg-D configuration, which is known to have difficulty in forming the edge transport barrier. The infernal-kink mode is generally found to be more unstable in neg-D plasmas as compared to their positive D-shaped (pos-D) counterpart. This is mainly due to less favorable (or even unfavorable) average magnetic curvature near the radial location of the minimum safety factor ( q min ) as compared to the pos-D configuration. The larger Shafranov shift associated with the neg-D shape helps the mode stabilization but is not sufficient to overcome the destabilizing effect due to bad curvature. Strong poloidal mode coupling due to plasma shaping (toroidicity, elongation, triangularity, etc.) helps explain the slight shift with respect to that predicted by the analytic theory of the peak location of the computed mode growth versus q min .

Zhang, X. M.↗

Heavy flavor jet substructure at lepton colliders

We provide a detailed analysis of event-shape observables, namely the energy correlation function and jet angularity, for heavy-flavor jets produced in electron-positron collisions, focusing on quantum chromodynamics (QCD) interactions. Using modern jet substructure techniques, we investigate the dead-cone effect, where QCD radiation is suppressed around a heavy quark within an angle proportional to its mass. Our analysis achieves next-to-leading logarithmic accuracy, combined with partial fixed-order contributions, to improve the description of the transition near the dead-cone threshold. To ensure a comprehensive perspective, we compare our analytical results with predictions from the pythia, herwig, and sherpa Monte Carlo simulations at past and future lepton colliders.

Perturbative QCD↗

Buckling-induced interaction between circular inclusions in an infinite thin plate

Design of slender artificial materials and morphogenesis of thin biological tissues typically involve stimulation of isolated regions (inclusions) in the growing body. These inclusions apply internal stresses on their surrounding areas that are ultimately relaxed by out-of-plane deformation (buckling). Here, we utilize the Föppl-von Kármán model to analyze the interaction between two circular inclusions in an infinite plate that their centers are separated a distance of 2ℓ. In particular, we investigate a region in phase space where buckling occurs at a narrow transition layer of length ℓ D around the radius of the inclusion, $R (ℓ_D \ll R)$. We show that the latter length scale defines two regions within the system, the close separation region, $ℓ – R \sim ℓ_D$, where the transition layers of the two inclusions approximately coalesce, and the far separation region, $ℓ – R \gg ℓ_D$. While the interaction energy decays exponentially in the latter region, $E_{\text{int}} ∝ e^{– (ℓ – R)/ℓ_D}$, it presents nonmonotonic behavior in the former region. While this exponential decay is predicted by our analytical analysis and agrees with the numerical observations, the close separation region is treated only numerically. In particular, we utilize the numerical investigation to explore two different scenarios within the final configuration: The first where the two inclusions buckle in the same direction (up-up solution) and the second where the two inclusions buckle in opposite directions (up-down solution). We show that the up-down solution is always energetically favorable over the up-up solution. In addition, we point to a curious symmetry breaking within the up-down scenario; we show that this solution becomes asymmetric in the close separation region.

42 ENGINEERING↗

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Learning Global Proliferation Expertise Evolution Using AI-Driven Analytics and Public Information

Detecting and anticipating global proliferation expertise and capability evolution from unstructured, noisy, and incomplete public data streams is a highly desired, but extremely challenging task. Here, in this article, we present our pioneering data-driven approach to support the non-proliferation mission to detect and explain the evolution of proliferation expertise and capability development globally from terabytes of publicly available information (PAI), focusing on our knowledge extraction pipeline and descriptive analytics. We first discuss how we fuse nine open-source data streams, including multilingual data, to convert 4 TB of unstructured data to structured knowledge and encode dynamically evolving proliferation expertise representations—content and context graphs. For this, we rely on natural language processing (NLP) and deep learning (DL) models to perform information extraction, topic modeling, and distributed text representation (aka embedding) learning. We then present interactive, usable, and explainable descriptive analytics to refine domain knowledge and present it in a human-understandable form. Finally, we introduce future work avenues that will leverage our dynamic knowledge representations and descriptive analytics to enable predictive and prescriptive inferences to achieve real-time domain understanding and contextual reasoning about global proliferation expertise and capability evolution.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗