Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PV System Modeling”

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 685 records · Page 38

Pricing Strategy of Electric Vehicle Aggregators Based on Locational Marginal Price to Minimize Photovoltaic (PV) Curtailment

The global climate crisis demands urgent action to mitigate global warming. Using renewable energy sources, such as solar and wind power, for electricity generation is crucial. This shift from centralized to distributed power systems, however, brings challenges, including voltage fluctuations and renewable energy curtailment. The rapid growth of the electric vehicle (EV) industry adds complexity, increasing overall electricity demand and straining the power supply during peak charging times. This paper proposes a scheduling strategy for EV aggregators to reduce renewable energy curtailment and stabilize grid operation by strategically scheduling EV charging. Using Multi -Agent Transport Simulation (MATSim), a traffic simulation tool, EV driving data in Denver, Colorado, USA, were modeled. The EV aggregator adjusts charging fees based on locational marginal prices, encouraging EVs to charge at different stations according to pricing. Simulations on an IEEE 33-bus system with distributed energy resources and EV charging stations validate the proposed algorithm, demonstrating its effectiveness in reducing curtailment by 12.55% and stabilizing grid operation.

33 ADVANCED PROPULSION SYSTEMS↗

Engineering improved ethylene production: Leveraging systems biology and adaptive laboratory evolution

Ethylene is a small hydrocarbon gas widely used in the chemical industry. Annual worldwide production currently exceeds 150 million tons, producing considerable amounts of CO 2 contributing to climate change. The need for a sustainable alternative is therefore imperative. Ethylene is natively produced by several different microorganisms, including Pseudomonas syringae pv. phaseolicola via a process catalyzed by the ethylene forming enzyme (EFE), subsequent heterologous expression of EFE has led to ethylene production in non-native bacterial hosts including E. coli and cyanobacteria. However, solubility of EFE and substrate availability remain rate limiting steps in biological ethylene production. We employed a combination of genome scale metabolic modelling, continuous fermentation, and protein evolution to enable the accelerated development of a high efficiency ethylene producing E. coli strain, yielding a 49-fold increase in production, the most significant improvement reported to date. Furthermore, we have clearly demonstrated that this increased yield resulted from metabolic adaptations that were uniquely linked to the EFE enzyme (WT vs mutant). Our findings provide a novel solution to deregulate metabolic bottlenecks in key pathways, which can be readily applied to address other engineering challenges.

59 BASIC BIOLOGICAL SCIENCES↗

Automated Shift Detection in Sensor-Based PV Power and Irradiance Time Series: Preprint

PV power and irradiance sensor-based measurements are prone to error, resulting in issues such as abrupt time series data shifts. These shifts, which are usually unintentional, may be caused by software or hardware configuration changes on a PV system, and do not reflect an actual change in overall system performance. Locating these shifts and segmenting the associated time series aids in more accurate future PV analysis. In this research, an offline changepoint detection (CPD) algorithm that automatically detects these abrupt data shifts in sensor-based time series is introduced. Data shift periods in 101 daily PV power and irradiance time series were labeled manually by two solar experts. These data streams represent sensor-based measurements, and display a variety of data shift behaviors. A changepoint detection algorithm was tuned using the 101 labeled data streams, with each model configuration's ability to detect labeled changepoints benchmarked using metrics such as F1-score, recall, and Rand Index. Best performing models on seasonality-corrected data streams include the Pruned Exact Linear (PELT) method, the Binary Segmentation method, and the Bottom-Up method, all scoring an average F1-score of 0.76 or greater at detecting labeled changepoints within a 30-day window for the labeled data sets. To promote further research in this space, we are releasing the labeled data shift sets on U.S. Department of Energy's (DOE) DuraMAT Data Hub, and the associated algorithm in the Python PVAnalytics package.

changepoint detection↗

CalTestBed - Delphire - Testing and Evaluation of Delphire Sentinel System (CRADA Final Report)

The Delphire Sentinel is a modular fire detection and communications system operating as a mobile field unit, with low voltage DC power supplied by onboard photovoltaics (PV) and batteries. The Sentinel addresses several aspects of fire detection, communications and data analysis. The Sentinel's mobility enables it to be rapidly deployed and operate independently of existing power and communications networks. The duration of independent operation depends critically on the energy consumption of the systems and performance of the onboard PV and battery. The purpose of this testing is to ascertain the power draw and energy consumption of the Delphire Sentinel prototype system under several operational states, including various data transfer packet sizes, transmission time and frequencies, and communication pathways (Wi-Fi, cellular, satellite) expected to be encountered in field deployments. It will also include procedures to test the ability of the Sentinel to operate for extended periods without loss of functionality. Based on results from energy and power measurements, and anticipated duty cycles in field deployments, we will model annual system autonomy (e.g. loss of load probability) for off-grid operation in representative locations.

47 OTHER INSTRUMENTATION↗

Grid Resiliency with a 100% Renewable Microgrid

San Diego Gas & Electric Company (SDG&E) installed America’s first and largest utility-scale microgrid in Borrego Springs in 2013. The first generation Borrego Springs Microgrid utilized diesel generators to form and stabilize the microgrid island, with support from grid-scale batteries and local solar photovoltaic (PV) generation. In this project, SDG&E in partnership with National Renewable Energy Laboratory (NREL) demonstrated through modeling, simulation and utility field testing that blackstart and islanding of the microgrid can be led with 100% renewable, inverter based resources (IBRs), to help reduce community reliance on conventional generation resources. Through equipment upgrades, grid-forming island leader capability was transitioned to a battery IBR instead of the Borrego Springs Microgrid diesel generators. A new microgrid controller was integrated to the microgrid and programmed to control and manage multiple energy storage systems. Synchrophasor and other power quality data verified autonomous, high-speed response of the IBRs through blackstart, islanding, and load step testing. Results of project field evaluations provide distribution systems operators (DSO) with increased confidence that renewable, IBR can replace traditional generators to blackstart and island microgrids and rapidly establish stable island frequency with rapid changes in peak power demand. Importantly, the project validated the integration feasibility of a distributed energy resource management system (DERMS) controller that manages multiple grid-forming and grid-following IBRs, establishing a standard design interface to reduce the complexity of integrating new DERs in the future and supporting replication by the industry. As a result of learnings in this project, SDG&E has implemented the microgrid controller strategy at multiple other microgrid sites, thereby validating the replicability of the solution. Hardware-in-the-loop (HIL) simulations including power and controller HIL hardware — along with electromagnetic transient (EMT) simulations of Borrego Springs Microgrid —informed adjustments to inverter parameters and were important to characterize the performance of the IBRs in relevant operating conditions before deployment. The EMT and HIL simulations of islanding the entire community are important contributions in providing confidence in IBR performance prior to future islanding of the community in the field. High-fidelity EMT and/or HIL simulation of IBRs can de-risk field operations, and its relevance and importance as a tool is increasing as distribution grids and microgrids become more complex and dynamic with an increasing proportion of renewable generation, distributed energy storage, and two-way power and energy flows.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Solar Thermal Energy Planner (STEP 1): A New Decision Support Tool for Solar Industrial Process Heat Applications

Solar thermal technologies are a promising technology to supply low-cost thermal energy to industrial processes, but there are often significant barriers to entry to industrial owners considering these technologies for their energy demands. To overcome this barrier and convey economic value to customers, NREL and Sandia National Laboratories developed Solar Thermal Energy Planner (STEP 1), a new web-based decision support tool for solar industrial process heat systems. At SolarPACES 2024, the STEP 1 tool was still under development; progress, methodologies, and a preliminary case study was presented. With the STEP 1 tool launch in May 2025, in this work, the initial version of the full public tool will be presented with demonstrations of its capabilities using a few case studies. First, the user's process heat needs such as location, process media (e.g., steam, air), process temperature, land availability, electricity and fuel costs, among other parameters. STEP 1 features a mapping interface that allows users to draw land and roof boundaries. The process media and temperature inform technology selection criteria modules that determine the appropriate solar thermal collection technologies, as well as congruent heat transfer media (e.g., hot water, oil, salt). Once the solar thermal technology selected, its nominal thermal production for the given site is characterized using NREL's System Advisor Model (SAM). Then, a modified version of NREL's REopt optimal sizing and dispatch optimization tool determines cost-optimal sizing. Within minutes, the user receives the results of the technoeconomics analysis, including the size and performance of the cost-optimal solar-plus-storage system. The cost of the system is compared to business-as-usual (e.g., an existing, standalone natural gas boiler). Users can download key results to store for sensitivity analyses. Examples of flat plate collector, parabolic trough, and molten salt tower applications with and without PV hybridization for different industrial facility types are presented in this work. The STEP 1 tool aims to reduce barriers to the adoption of solar heating solutions stemming from a lack of familiarity and technical background with solar system design options and costs among industry stakeholders.

14 SOLAR ENERGY↗

Automated Shift Detection in Sensor-Based PV Power and Irradiance Time Series

PV power and irradiance sensor-based measurements are prone to error, resulting in issues such as abrupt time series data shifts. These shifts, which are usually unintentional, may be caused by software or hardware configuration changes on a PV system, and do not reflect an actual change in overall system performance. Locating these shifts and segmenting the associated time series aids in more accurate future PV analysis. In this research, an offline changepoint detection (CPD) algorithm that automatically detects these abrupt data shifts in sensor-based time series is introduced. Data shift periods in 101 daily PV power and irradiance time series were labeled manually by two solar experts. These data streams represent sensor-based measurements, and display a variety of data shift behaviors. A changepoint detection algorithm was tuned using the 101 labeled data streams, with each model configuration's ability to detect labeled changepoints benchmarked using metrics such as F1-score, recall, and Rand Index. Best performing models on seasonality-corrected data streams include the Pruned Exact Linear (PELT) method, the Binary Segmentation method, and the Bottom-Up method, all scoring an average F1-score of 0.76 or greater at detecting labeled changepoints within a 30-day window for the labeled data sets. To promote further research in this space, we are releasing the labeled data shift sets on U.S. Department of Energy's (DOE) DuraMAT Data Hub, and the associated algorithm in the Python PVAnalytics package.

changepoint detection↗

PV Stormwater Management Research and Testing (PV-SMaRT) (Final Technical Report)

The objective of the Photovoltaic Stormwater Management Research and Testing (PV-SMaRT) project was to develop and disseminate research-based, solar-specific resources for estimating stormwater runoff at ground-mounted PV facilities and detail stormwater management and water quality best practices. The intended use of stormwater management and water quality best practices and stormwater runoff estimation resources is to reduce balance of system soft costs associated with stormwater infrastructure requirements and improve water quality through research-tested best practices. Currently, stormwater regulations and guidelines vary by jurisdiction and can lead to regulatory uncertainty and increased compliance costs for managing stormwater runoff at solar sites. To address these concerns, the NREL team and its partners (University of Minnesota, Great Plains Institute, and Fresh Energy) established and engaged an advisory Water Quality Task Force (WQTF); conducted field research on stormwater infiltration and runoff at five ground-mounted PV sites; validated a 3-D hydrologic model to predict water runoff and generate stormwater runoff coefficients for a range of site conditions and PV designs; developed a stormwater management and water quality best practices document; and engaged with local jurisdictions and other stakeholders to disseminate best practices, stormwater runoff coefficients, and other tools. Key outputs of this project were a PV-SMaRT Runoff Calculator developed by University of Minnesota, a webinar detailing project outcomes and how to use the PV-SMaRT runoff calculator, and a document on Best Practices for regulators and developers regarding stormwater runoff at PV sites.

14 SOLAR ENERGY↗

Montana Community-Scale Solar Strategy Project (Final Technical Report)

Grid-tied solar photovoltaic (PV) development began gradually in Montana, initially taking root in the early 2000s, largely in the form of small residential systems, rarely larger than two kilowatts. Utility support for solar technology through public benefit funding triggered the installation of the state’s first solar PV installations on community-owned facilities including schools, libraries, fire stations and local government buildings in the mid-2000s. By 2016, solar PV installations totaled approximately 28 megawatts of installed capacity, generating less than 0.05% of electricity sold in Montana. Montana’s solar fleet in 2016 included 17 megawatts of utility-scale solar farms, and three shared solar projects developed by rural electric cooperatives that ranged from 25 to 50 kilowatts in size. However, with only 3,000 households utilizing solar PV, either through an on-site installation or participation in a shared solar project, the potential benefits and reach of solar technology was limited. The Montana Solar Community Project (MSCP) was developed by the Montana Energy Office (MEO) at the Department of Environmental Quality in order to implement a “community-scale” solar energy strategy for the state that would cost-effectively expand Montanans’ access to solar PV. This project defines three types of community-scale solar projects: shared solar, community-sited solar, and group purchasing programs. The MSCP project conducted research and stakeholder engagement that helped to inform and develop model community-scale initiatives and an overarching community-scale solar strategy for Montana. Specifically, the project included a solar market assessment, stakeholder meetings across Montana, development of a menu of community-scale solar options, a study of solar potential on schools, community solar market research, and community-scale project implementation support. MSCP activities were successful in implementing solar projects and supporting solar development in communities across Montana. Project activities realized 280 kW of newly installed solar from two solarize programs, a shared solar array at a rural electric cooperative, and a solar array atop a Bozeman high school. The community-scale project implementation support and Solar on Schools components of the project resulted directly in the analysis and design of nearly 2 MW of distributed solar systems across the state. The Solar on Schools analysis alone produced solar designs for 25 schools, representing 920.8 kW of solar potential. Many of the deliverables from this project provide guidance on future development of community-scale solar in Montana, and resulted in a number of important conclusions about community-scale solar in Montana. There is strong interest in community-scale solar in Montana, as noted during stakeholder meetings, in response to solar market research surveys, and as seen by applications for implementation support grants provided by MEO. Project activities also suggest that Montanans demonstrated a desire and need for more information about solar energy technology, solar financing options, and the state’s policy and regulatory framework. Additionally, the project found that schools will be key partners for future community-scale solar programs. In addition to being excellent hosts for solar arrays, schools have the opportunity to use the installation of a solar array to develop educational materials and curricula to educate Montana’s future leaders about these technologies, which can be shared and replicated across the state. Project activities laid a clear path forward for solar energy in Montana. Certain activities have concrete next steps. For example, numerous solar arrays were designed for schools and other public buildings around the state. Many of those designs will be used as soon as funding becomes available. Other activities have less developed next steps, though lay the foundation for future work. The solar market assessment, solar market research, and community stakeholder meetings are good examples. These documents can be used as stepping stones to increase renewable energy education, increase stakeholder outreach, and support communities as they seek to develop community-scale solar projects.

14 SOLAR ENERGY↗

World Record Demonstration of > 30% Thermophotovoltaic Conversion Efficiency

Thermophotovoltaic (TPV) devices are solid-state heat engines that convert thermal radiation into electricity using semiconductor diodes, and have applications in energy storage, primary power conversion, and waste heat recovery. Higher temperature emitters and wider band gap photovoltaics (PVs) generally afford higher TPV efficiencies. High temperature emitters are accessible in thermal energy storage applications. We have demonstrated a world record 31% +/- 2% TPV conversion efficiency with a 0.9 cm2 GaAs-based PV device under a 2430 °C thermal emitter, producing an electrical output power of 2.23 W. Critical to the result was the cell’s high reflectance of photon energies below the device band gap. Unlike solar PV, for which sub-band gap (SBG) light is lost, a TPV cell can reflect and recycle SBG light to the thermal emitter. The demonstration was made on a custom-built measurement platform in which a ~100 cm2 graphite thermal emitter was heated under vacuum up to 2430 °C. The TPV efficiency was evaluated by comparing the measured electrical output power of the TPV with the calorimetrically measured total net power incident on the cell. The measured TPV efficiency as a function of thermal emitter temperature was corroborated by our full system modeling predictions. As far as the authors are aware, this is the highest TPV conversion efficiency ever measured, and device improvements should yield > 40% efficiency in the near future.

27 ARPA - Advanced Research Projects Agency-Energy↗

Generalized Analytical Estimation of Sensitivity Matrices in Unbalanced Distribution Networks

Fast and accurate estimation of sensitivity matrices is significant for the enhancement of distribution system modeling and automation. Analytical estimations have mainly focused on voltage magnitude sensitivity to active/reactive power injections for unbalanced networks with Wye-connected loads and neglecting DERs' smart inverter functionality. Hence, this paper enhances the scope of analytical estimation of sensitivity matrices for unbalanced networks with 1- Φ, 2- Φ, and 3- Φ Delta/Wye-connected loads, DERs with smart inverter functionality, and substation/line step-voltage regulators (SVR). A composite bus model comprising of DER, Delta- and Wye-connected load is proposed to represent a generic distribution bus, which can be simplified to load, PV, or voltage-controlled bus as required. Furthermore, the proposed matrix-based analytical method consolidates voltage magnitude and angle sensitivity to active/reactive power injection and tap-position of all SVRs into a single algorithm. Extensive case studies on IEEE and EPRI networks show the accuracy and wide scope of the proposed algorithm compared to the existing benchmark method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Debris-induced consequences on turbulence and vorticity in solar photovoltaic module-generated array wakes

Particle-laden flows in solar photovoltaic (PV) systems are inevitable, where wind-swept debris in open environments are carried by high winds and turbulence, coating panel surfaces or damaging structures. Particle deposition, or soiling, is a well-known issue for large-scale plants which rely on uninhibited solar rays for optimal production. But understanding the mechanisms leading to soiling requires a physical and fluid dynamics-centered focus, since turbulence dominates PV panel wakes and is also known to alter particle concentration and trajectories. This study presents an experimental campaign toward consequences of particle-laden flow between two model PV panels using time-resolved particle image velocimetry. The model array was subjected to varied particle volume fractions, including a tracer particle case and a water droplet case. Characterization of mean velocity, turbulence statistics, and mean kinetic energy within the single phase and, separately, particle phase flows showed modified features due to particle inertia. Images captured at a frequency of 1 kHz in the near wake of the upstream panel allow for a first experimental look at vorticity and convective velocity of vortex structures for single-phase and particle-phase flows which are crucial to debris transport and soiling in PV environments.

Energy & Fuels↗

Real-time Simulation Framework for Hardware-in-the-Loop Testing of Multi-port Autonomous Reconfigurable Solar Power Plant (MARS)

Multiport autonomous reconfigurable solar power plant (MARS) has been proposed for integrated development of photovoltaic (PV) and energy storage system (ESS) that can connect to high-voltage direct current (HVdc) and alternating current (ac) transmission grid. To de-risk the development of this complex integrated system that consists of hundreds to thousands of power electronics modules, a controller hardware-in-the-loop (cHIL) test setup will be extremely beneficial. The cHIL testing can be used for evaluation of modules as well as the hierarchical control system in MARS. With the unique configuration of power electronics modules in MARS, it becomes necessary to develop custom-designed real-time simulation models in the cHIL setup in absence of off-the-shelf models. In this paper, high-fidelity dynamic model of MARS, control algorithms at the lower level, and required communication algorithms are developed and optimized for real-time performance in the cHIL setup. Real-time experimental results from the cHIL are provided.

Dong, Zerui↗

Stratosphere-to-Troposphere Transport Revealed by Ground-based Lidar and Ozonesonde at a Midlatitude Site

This paper presents ozone structures measured by a ground-based ozone lidar and ozonesonde at Huntsville, Alabama, on 27-29 April 2010 originating from a stratosphere-to-troposphere transport event associated with a cutoff cyclone and tropopause fold. In this case, the tropopause reached 6 km and the stratospheric intrusion resulted in a 2-km thick elevated ozone layer with values between 70 and 85 ppbv descending from the 306-K to 298-K isentropic surface at a rate of 5 km day1. The potential temperature was provided by a collocated microwave profiling radiometer. We examine the corresponding meteorological fields and potential vorticity (PV) structures derived from the analysis data from the North American Mesoscale model. The 2-PVU (PV unit) surface, defined as the dynamic tropopause, is able to capture the variations of the ozone tropopause estimated from the ozonesonde and lidar measurements. The estimated ozone/PV ratio, from the measured ozone and model derived PV, for the mixing layer between the troposphere and stratosphere is approximately 41 ppbv/PVU with an uncertainty of approximately 33%. Within two days, the estimated mass of ozone irreversibly transported from the stratospheric into the troposphere is between 0.07 Tg (0.9 10(exp33) molecules) and 0.11 Tg (1.3 10(exp33) molecules) with an estimated uncertainty of 59%. Tropospheric ozone exhibited enormous variability due to the complicated mixing processes. Low ozone and large variability were observed in the mid-troposphere after the stratospheric intrusion due to the westerly advection including the transition from a cyclonic system to an anticyclonic system. This study using high temporal and vertical-resolution measurements suggests that, in this case, stratospheric air quickly lost its stratospheric characteristics once it is irreversibly mixed down into the troposphere.

Kuang, Shi↗

A Method for Estimating Time-Series PV Production Loss From Solar Tracking Failures

Solar tracking system failures can dramatically reduce photovoltaic (PV) system energy output through suboptimal irradiance capture and self-shading loss. Detecting tracker failures is not always straightforward; unlike other common PV system failures like inverter outages and blown fuses which take portions of a system completely offline, tracker failures only partially reduce output power and can go unnoticed as a result. Additionally, tracker failures can take one of several forms, each with their own loss characteristics. Here we present two methods of detecting tracker failure events from time-series production data and a method for estimating the associated production loss. Compared with existing detection methods, the proposed power-based detection method showed 0-17% reductions in Type I error rate depending on weather conditions. Here, the loss model estimated production loss with low error (mean bias error = -2.3 %, root-mean-squared error = 6%) in a rudimentary validation.

14 SOLAR ENERGY↗

Analysis of a Disturbance Event with Inverter-Based Resources Using EMT Simulations

Increasing penetration of inverter-based resources (IBRs) necessitates newer methods of planning and analysis of disturbances. The existing phasor-domain transient stability (TS) analysis may not capture the dynamics of IBRs during fault events. Here, in this paper, electromagnetic transient (EMT) simulations using high-fidelity detailed model of power grid and one of the affected photovoltaic (PV) plants during the Angeles Forest disturbance in 2018 are performed. In these simulations, the processes to develop EMT models of power grid from traditional phasor-domain TS data and PV plant from collected data are described. Thereafter, using these simulations, the response of the PV plant during the fault event in 2018 is replicated and a sensitivity analysis is performed. The sensitivity analysis consists of making changes to the components within the PV plant and in the power grid to evaluate the impact they have on the response observed by the PV plant during the fault event. This analysis provides an understanding of the components that impact the operation of a PV plant during fault events and provide guidance to system planners on the studies that need to be performed to maintain a reliable power grid as new IBR plants are integrated.

42 ENGINEERING↗

Inventory for Crystalline Silicon Module Recycling: Cooperative Research and Development (Final Report)

A critical challenge for the continued expansion of photovoltaics (PV) is to develop technically feasible, inexpensive and environmentally friendly practices for handling and recycling modules at the end of their usable life. The National Renewable Energy Laboratory (NREL) is requested by the Electric Power Research Institute (EPRI) to collect primary data regarding the environmental performance of currently operational PV module recycling facilities in Europe. Very little has been published regarding crystalline silicon (C-Si) module recycling. Thus, much effort will be needed in direct industry outreach, collection of information and other business intelligence strategies similar to NREL's approaches for developing cost models for PV manufacturing. The goal of this work effort is to produce a detailed inventory that accounts for physical (e.g., energy, water, materials) flows through each step of a C-Si recycling process. The inventory (a life cycle inventory, or LCI) shall be designed so that it can be extended to include an accounting of costs for each process step, inputs, etc. This work effort shall leverage prior LCI data collection NREL performed for the United States Department of Energy, Solar Energy Technologies Office, under the auspices of the U.S. contribution to International Energy Agency's Photovoltaics Power Systems (PVPS) Task 12 (Environmental Health and Safety), which SETO nominated NREL to chair. The primary purpose of this work effort is to augment the prior data collection to increase the sample size of manufacturers' primary data in the LCI.

14 SOLAR ENERGY↗

Storage Futures Study: Four Phases Framework and Modeling

This webinar, presented by Paul Denholm and Will Frazier, incorporates material from two publications of the NREL Storage Futures Study: "The Four Phases of Storage Deployment: A Framework for the Expanding Role of Storage in the U.S. Power System" (https://www.nrel.gov/docs/fy21osti/77480.pdf) and "Storage Futures Study: Economic Potential of Diurnal Storage in the U.S. Power Sector" (https://www.nrel.gov/docs/fy21osti/77449.pdf).

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗