Search NASA⌕ Search

SEARCH · Search NASA

Results for “developable planes”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

PV Rooftop Database

The National Renewable Energy Laboratory's (NREL) Photovoltaic (PV) Rooftop Database (PVRDB) is a lidar-derived, geospatially-resolved dataset of suitable roof surfaces and their PV technical potential for 128 metropolitan regions in the United States. The PVRDB data are organized by city and year of lidar collection. Five geospatial layers are available for each city and year: 1) the raster extent of the lidar collection, 2) buildings identified from the lidar data, 3) suitable developable planes for each building, 4) aspect values of the developable planes, and 5) the technical potential estimates of the developable planes.

Array↗
DOE OSTI2021

Embedded training of neural-network subgrid-scale turbulence models

We report that the weights of a deep neural-network model are optimized in conjunction with the governing flow equations to provide a model for subgrid-scale stresses in a temporally developing plane turbulent jet at Reynolds number Re 0 = 6000 . The objective function for training is first based on the instantaneous filtered velocity fields from a corresponding direct numerical simulation, and the training is by a stochastic gradient descent method, which uses the adjoint Navier-Stokes equations to provide the end-to-end sensitivities of the model weights to the velocity fields. In-sample and out-of-sample testing on multiple dual-jet configurations show that its required mesh density in each coordinate direction for prediction of mean flow, Reynolds stresses, and spectra is half that needed by the dynamic Smagorinsky model for comparable accuracy. The same neural-network model trained directly to match filtered subgrid-scale stresses, without the constraint of being embedded within the flow equations during the training, fails to provide a qualitatively correct prediction. The coupled formulation is generalized to train based only on mean-flow and Reynolds stresses, which are more readily available in experiments. The mean-flow training provides a robust model, which is important, though a somewhat less accurate prediction for the same coarse meshes, as might be anticipated due to the reduced information available for training in this case. The anticipated advantage of the formulation is that the inclusion of resolved physics in the training increases its capacity to extrapolate. This is assessed for the case of passive scalar transport, for which it outperforms established models due to improved mixing predictions.

42 ENGINEERING↗
DOE OSTI2022

Maximum strength and dislocation patterning in multi–principal element alloys

Multi–principal element alloys (MPEAs) containing three or more components in high concentrations render a tunable chemical short-range order (SRO). Leveraging large-scale atomistic simulations, we probe the limit of Hall-Petch strengthening and deformation mechanisms in a model CrCoNi alloy and unravel chemical ordering effects. The presence of SRO appreciably increases the maximum strength and lowers the propensity for faulting and structure transformation, accompanied by intensification of planar slip and strain localization. Deformation grains exhibit notably different microstructures and dislocation patterns that prominently depend on their crystallographic orientation and the number of active slip planes. Grain of single-planar slip attains the highest volume fraction of deformation-induced structure transformation, and grain with double-slip planes develops the densest dislocation network. These results advancing the fundamental understanding of deformation mechanisms and dislocation patterning in MPEAs suggest a mechanistic strategy for tuning mechanical behavior through simultaneously tailoring grain texture and local chemical order.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Generating porosity during olivine carbonation via dissolution channels and expansion cracks

Abstract. The olivine carbonation reaction, in which carbon dioxide is chemically incorporated to form carbonate, is central to the emerging carbon sequestration method using ultramafic rocks. The rate of this retrograde metamorphic reaction is controlled, in part, by the available reactive surface area: as the solid volume increases during carbonation, the feasibility of this method ultimately depends on the maintenance of porosity and the creation of new reactive surfaces. We conducted in situ dynamic X-ray microtomography and nanotomography experiments to image and quantify the porosity generation during olivine carbonation. We designed a sample setup that included a thick-walled cup (made of porous olivine aggregates with a mean grain size of either ∼ 5 or ∼ 80 µm) filled with loose olivine sands with grain sizes of 100–500 µm. The whole sample assembly was reacted with a NaHCO 3 aqueous solution at 200 °C, under a constant confining pressure of 13 MPa and a pore pressure of 10 MPa. Using synchrotron-based X-ray microtomography, the three-dimensional (3-D) pore structure evolution of the carbonating olivine cup was documented until the olivine aggregates became disintegrated. The dynamic microtomography data show a volume reduction in olivine at the beginning of the reaction, indicating a vigorous dissolution process consistent with the disequilibrium reaction kinetics. In the olivine cup with a grain size of ∼ 80 µm (coarse-grained cup), dissolution planes developed within 30 h, before any precipitation was observed. In the experiment with the olivine cup of ∼ 5 µm mean grain size (fine-grained cup), idiomorphic magnesite crystals were observed on the surface of the olivine sands. The magnesite shows a near-constant growth throughout the experiment, suggesting that the reaction is self-sustained. Large fractures were generated as the reaction proceeded and eventually disintegrated the aggregate after 140 h. Detailed analysis show that these are expansion cracks caused by the volume mismatch in the cup walls, between the expanding interior and the near-surface which keeps a nearly constant volume. Nanotomography images of the reacted olivine cup reveal pervasive etch pits and wormholes in the olivine grains. We interpret this perforation of the solids to provide continuous fluid access, which is likely key to the complete carbonation observed in nature. Reactions proceeding through the formation of nano- to micron-scale dissolution channels provide a viable microscale mechanism in carbon sequestration practices. For the natural peridotite carbonation, a coupled mechanism of dissolution and reaction-induced fracturing should account for the observed self-sustainability of the reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Effect of annealing ambient on anisotropic retraction of film edges during solid-state dewetting of thin single crystal films

During solid-state dewetting of thin single crystal films, film edges retract at a rate that is strongly dependent on their crystallographic orientations. Edges with kinetically stable in-plane orientations remain straight as they retract, while those with other in-plane orientations develop in-plane facets as they retract. Kinetically stable edges have retraction rates that are lower than edges with other orientations and thus determine the shape of the natural holes that form during solid-state dewetting. In this paper, measurements of the retraction rates of kinetically stable edges for single crystal (110) and (100) Ni films on MgO are presented. Relative retraction rates of kinetically stable edges with different crystallographic orientations are observed to change under different annealing conditions, and this accordingly changes the initial shapes of growing holes. The surfaces of (110) and (100) films were also characterized using low energy electron diffraction, and different surface reconstructions were observed under different ambient conditions. The observed surface structures were found to correlate with the observed changes in the relative retraction rates of the kinetically stable edges.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Flow Development in Radial Plane of Rotating Detonation Engine Integrated with Aerospike

In this study, axial and radial velocity components are measured downstream of a rotating detonation engine (RDE) integrated with an aerospike by using particle image velocimetry (PIV) at 100 kHz. The RDE is operated at high pressures obtained by restricting the RDE exit with a converging nozzle, which also directs the flow radially toward the aerospike. Reactants, methane fuel, and enriched air (67% oxygen and 33% nitrogen) are supplied from separate plenums at ambient temperature. PIV data are presented to investigate how pressure ratios within the RDE channel affect the flow evolution across the aerospike. The RDE is shown to operate stably and consistently in a single wave mode for all six test runs. The study highlights large temporal and spatial variations in both axial and radial flow velocities at the nozzle throat that persist downstream across the aerospike. Each point in the flowfield oscillates at the frequency recorded inside the RDE channel. Overall, the temporally and spatially varying aerospike flowfield is far from the ideal case of a uniform flow at its exit, and it points toward the need to condition the flow within the RDE channel to produce more uniform conditions at the nozzle throat.

Engineering↗
DOE OSTI2019

“The development of a through-plane reactive excitation technique for detection of pinholes in membrane-containing MEA sub-assemblies”

As markets for polymer-electrolyte-membrane fuel cells (PEMFC) expand, there continues to be a need for the development and implementation of in-line quality inspection techniques to enable high-volume manufacturing methods for membrane-electrode assemblies (MEAs). These techniques mitigate the potential impact of poor quality on the cost of MEA components by real-time identification of process-induced irregularities, which can result from the manufacturing process or subsequent handling of the material during assembly into a full MEA. A particular type of irregularity known to impact MEA lifetime is loss of membrane integrity, e.g. via a pinhole. In this work, we focus on development of a technique to detect pinholes in membrane-containing MEA sub-assemblies - as opposed to in just the membrane - using a reactive excitation strategy coupled with infrared thermography for detection. A specialized device is introduced and, using various MEA sub-assemblies, is shown to enable detection of pinholes in stationary samples. Multi-physics modeling is utilized to understand and predict the impact of various design and operational parameters of this technique. The feasibility of the technique is demonstrated through both modeling and in-situ results. For example, a detectable (≥1 degree C) thermal response for a membrane pinhole of approximately 90 um in diameter is measured in a fraction of a second. The observed detection time is considered to be appropriate for further exploration of the technique in an in-line configuration.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Technical Performance of Refractory Liners for Molten Chloride Salt Thermal Energy Storage Systems

A chloride-based molten-salt system that uses a ternary blend of MgCl2/KCl/NaCl is investigated to provide higher temperature thermal energy storage capability. Despite higher thermal stability, molten chlorides present several unique challenges, including the design of internal refractory-ceramic liners to prevent the corrosion and thermal stress of alloy tank shells. This work discusses issues and potential solutions related to containment of molten chloride salt, specifically the optimization of the refractory material at the molten salt interface (hot face). The down-selected hot face candidate limits permeation of salt through the material and forms a highly stable secondary surface phase in equilibrium with the molten salt. A mortar is created using the corrosion resistant hot face brick. Brick and mortar composite are subjected to mechanical stress/strain analysis, in order to calculate composite material properties and better inform thermomechanical models. The U.S. Department of Energy Generation 3 (DOE Gen3) program seeks to develop higher efficiency CSP plants that can provide cost-competitive, flexible power in the U.S. electric grid. The proposed Gen3 Liquid Pathway CSP plant closely resembles the configuration of current nitrate salt power towers with two-tank storage (Gen2). The differences between Gen2 and Gen3 include the types of compatible materials used in salt storage tank construction. Stainless steel loses strength at Gen3 temperatures, and although nickel superalloys would be capable of withstanding sustained high temperatures, these materials are prohibitively expensive at scale. Uninsulated tank shells also pose a significant risk as common steels are highly susceptible to chemical attack from molten chloride salt. To address these concerns, refractory-ceramic based containment materials are proposed to line the inside of the hot and cold storage tanks. In doing so, stainless or carbon steel shells may be used in construction depending on the level of insulation provided. The composition of the internal liner requires careful consideration to maximize the efficacy of multiple parameters including corrosion resistance, strength at operating temperature, durability, and cost. This is particularly true for the material at the interface with the salt, known as the "hot face", which is responsible for protecting the insulating layers between the tank shell and the hot face brick layer. The molten salt in this system is superheated over 300 °C above its freeze temperature. Therefore, unlike other industrial processes which use refractory-lined vessels, it is not expected that a freeze plane will develop in the hot face. Therefore, the hot face must be designed to withstand chemical corrosion and inhibit permeation of molten salt into the insulating layers. A down selection was performed to identify a hot face candidate best equipped to maintain thermal, mechanical, and chemical integrity when exposed to molten salt over extended periods of time. Long-duration chemical capability experiments were conducted with the down selected hot face refractory fully immersed in molten chloride salt for up to 3000 hours. The average salt penetration does not exceed 100 microns. When extrapolated to 20 and 30 years of continuous exposure, the expected salt penetration depth is approximately 2.0 mm and 2.9 mm, respectively. A magnesium-rich secondary phases develops at the salt/refractory interface. X-ray diffraction identifies the material as forsterite (Mg2SiO4), which is reported to form synthetically in molten chloride salt solutions. These results suggest the selected hot face will adequately inhibit salt permeation. While there is optimism that the hot face brick will inhibit salt permeation, mortar joints are typically the weakest point of a refractory brick liner. From a thermochemical perspective, differing thermal expansion coefficients may result in the mortar and brick to grow independent of each other, creating gaps through which molten salt can penetrate. To address this issue, NREL has developed an in-house mortar composed of the down selected hot face brick that has been shown to be compatible with the salt. Compressive stress/strain analyses have been performed on the brick/ mortar composites to generate stress/strain curves. Modulus of elasticity and Poisson's ratio of the composite may be calculated from the stress/strain curves, in order to provide more representative data to finite element mechanical models for accurate approximation of stress on the tank shell and the amount of thermal expansion expected within the tank liner.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2026

Equilibrium and non-equilibrium effects in high pressure phase transformations of carbon

The behavior of carbon in the range 1–100 GPa and 1–10 kK is central to problems in planetary interiors, inertial confinement fusion targets, and high-pressure synthesis of carbon-based materials, but experiments in this regime are difficult and often provide only indirect constraints on phase behavior. As a result, phase boundary loci, structure, and limits of metastability at high pressure remain uncertain. In this work, machine-learning enhanced atomistic simulations are used to address this knowledge gap. We determine the melt line up to 100 GPa, the graphite-diamond phase boundary up to the melt line, and analyze structure of the coexisting phases. We show that the coexisting liquid evolves smoothly with pressure without evidence for a first-order liquid–liquid transition. Orientation-resolved graphite melting simulations indicate that basal-plane interfaces develop a dewetting layer and undergo layer-by-layer melting, producing kinetic hysteresis and an apparent orientation dependence of the melt line. Non-equilibrium quenches from the melt are used to construct a kinetically limiting graphite–diamond phase boundary for rapid quenches from above the melt line, and show that graphite is metastable at pressures of up to ≈ 25 GPa. These results provide bounds on equilibrium and metastable behavior in carbon relevant for interpreting high-pressure experiments and for designing synthesis pathways to specific carbon microstructures.

Lyu, Yanjun [Department of Materials Science and E↗
DOE OSTI2016

SU-C-209-05: Monte Carlo Model of a Prototype Backscatter X-Ray (BSX) Imager for Projective and Selective Object-Plane Imaging

Purpose: To develop a Monte Carlo N-Particle (MCNP) model for the validation of a prototype backscatter x-ray (BSX) imager, and optimization of BSX technology for medical applications, including selective object-plane imaging. Methods: BSX is an emerging technology that represents an alternative to conventional computed tomography (CT) and projective digital radiography (DR). It employs detectors located on the same side as the incident x-ray source, making use of backscatter and avoiding ring geometry to enclose the imaging object. Current BSX imagers suffer from low spatial resolution. A MCNP model was designed to replicate a BSX prototype used for flaw detection in industrial materials. This prototype consisted of a 1.5mm diameter 60kVp pencil beam surrounded by a ring of four 5.0cm diameter NaI scintillation detectors. The imaging phantom consisted of a 2.9cm thick aluminum plate with five 0.6cm diameter holes drilled halfway. The experimental image was created using a raster scanning motion (in 1.5mm increments). Results: A qualitative comparison between the physical and simulated images showed very good agreement with 1.5mm spatial resolution in plane perpendicular to incident x-ray beam. The MCNP model developed the concept of radiography by selective plane detection (RSPD) for BSX, whereby specific object planes can be imaged by varying kVp. 10keV increments in mean x-ray energy yielded 4mm thick slice resolution in the phantom. Image resolution in the MCNP model can be further increased by increasing the number of detectors, and decreasing raster step size. Conclusion: MCNP modelling was used to validate a prototype BSX imager and introduce the RSPD concept, allowing for selective object-plane imaging. There was very good visual agreement between the experimental and MCNP imaging. Beyond optimizing system parameters for the existing prototype, new geometries can be investigated for volumetric image acquisition in medical applications. This material is based upon work supported under an Integrated University Program Graduate Fellowship sponsored by the Department of Energy Office of Nuclear Energy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Continuous Test and Transition Infrastructure for Quantum Networking for Science Complex

We propose a network architecture with a separate quantum dataplane and a conventional control plane: (a) Quantum Data Plane: The Quantum Data plane consists of links of dark fibers connecting quantum switches and repeaters, which in turn, connect to quantum computers, memory and sensors. Since the current reach is limited to local areas, it will begin as a collection of site networks at laboratories, (b) Conventional Control Plane: The Control plane provides management access to quantum devices for configuration and provisioning via control nodes with firewall and encryption capabilities. Continuous Test and Transition Infrastructure: We propose an infrastructure with a control plane connecting multiple sites via encrypted tunnels over conventional networks consisting of the following: (a) Site Quantum Networks: Individual site networks supported by their fiber plants connect laboratories that house and connect to their quantum devices for testing and interoperability.(b) Site Control Planes: Sites are connected over individual control planes with control hosts with Software Defined Networking capabilities, which can be peered with other networks via ESnet.(c) Progressive Expansion: Initially, sites will develop individual data planes with their specific quantum devices and fiber connections, and mechanisms to interface with conventional networks. They progressively expand and interconnect, under a wide-area ecosystem of peered control planes for device testing, interoperability development and roll off into production environments.

Rao, Nageswara S.↗
DOE OSTI2018

Rhenium-Doped and Stabilized MoS 2 Atomic Layers with Basal-Plane Catalytic Activity

The development of stable and efficient hydrogen evolution reaction (HER) catalysts is essential for the production of hydrogen as a clean energy resource. A combination of experiment and theory demonstrates that the normally inert basal planes of 2D layers of MoS 2 can be made highly catalytically active for the HER when alloyed with rhenium (Re). The presence of Re at the ≈50% level converts the material to a stable distorted tetragonal (DT) structure that shows enhanced HER activity as compared to most of the MoS 2 -based catalysts reported in the literature. More importantly, this new alloy catalyst shows much better stability over time and cycling than lithiated 1T-MoS 2 . Density functional theory calculations find that the role of Re is only to stabilize the DT structure, while catalysis occurs primarily in local Mo-rich DT configurations, where the HER catalytic activity is very close to that in Pt. The study provides a new strategy to improve the overall HER performance of MoS 2 -based materials via chemical doping.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multiscale analysis of large twist ferroelectricity and swirling dislocations in bilayer hexagonal boron nitride

With its atomically thin structure and intrinsic ferroelectric properties, heterodeformed bilayer hexagonal boron nitride (hBN) has gained prominence in next-generation non-volatile memory applications. However, studies to date have focused almost exclusively on small-twist bilayer hBN, leaving the question of whether ferroelectricity can persist under small heterostrain and large heterodeformation entirely unexplored. In this work, we establish the crystallographic origin of ferroelectricity in bilayer hBN configurations heterodeformed relative to high-symmetry configurations such as AA-stacking and 21.786789° twisted configurations (Σ7), using Smith normal form bicrystallography. We then demonstrate out-of-plane ferroelectricity in bilayer hBN across configurations vicinal to both the AA and Σ7 stackings. Atomistic simulations reveal that AA-vicinal systems support ferroelectricity under both small twist and small strain, with polarization switching in the latter governed by the deformation of swirling dislocations rather than the straight interface dislocations seen in the former. For Σ7-vicinal systems, where existing interatomic potentials underperform particularly under extreme out-of-plane compression, we develop a density-functional-theory-informed continuum framework—the bicrystallography-informed frame-invariant multiscale (BFIM) model, which captures out-of-plane ferroelectricity in heterodeformed configurations vicinal to Σ7 stacking. Interface dislocations in these large heterodeformed bilayer configurations exhibit markedly smaller Burgers vectors compared to interface dislocations in small-twist and small-strain bilayer hBN. The BFIM model reproduces experimental results and provides a powerful, computationally efficient framework for predicting ferroelectricity in large-unit-cell heterostructures where atomistic simulations are prohibitively expensive.

Ahmed, Md Tusher [Univ. of Illinois at Urbana-Cham↗
DOE OSTI2016

SU-F-T-259: GPR Tables for the Estimation of Mid-Plane Dose Using EPID

Purpose: To develop a simple method for estimating the mid-plane dose (MPD) of a patient using Electronic Portal imaging Device (EPID). Methods: A Varian TrueBeam with aSi100 EPID was used in this study. The EPID images were acquired for a 30 cm × 30 cm homogeneous slab phantom and a 30 cm diameter 20 cm thick cylindrical phantom in the continuous dosimetry mode. The acquired EPID images in XIM format were imported into in-house MATLAB program for the data analysis. First, the dosimetric characteristics of EPID were studied for dose-response linearity, dose-rate dependence, and field size dependence. Next, the average pixels values of the EPID images were correlated with the MPD measured by an ionisation chamber for various thicknesses of the slab phantom (8 cm – 30 cm) and for various square field sizes (3×3 cm{sup 2} – 25×25 cm{sup 2} at the isocenter). Look-up tables called as GPR tables were then generated for both SSD and SAD setup by taking the ratio of MPD measured by the ionisation chamber and the corresponding EPID pixel values. The accuracy of the GPR tables was evaluated by varying the field size, phantom thickness, and wedge angles with the slab and cylindrical phantoms. Results: The dose response of EPID was linear from 20 MU to 300 MU. The EPID response for different dose rates from 40 MU/min to 600 MU/min was within ±1%. The difference in the doses from the GPR tables and the doses measured by the ionization chambers were within 2% for slab phantoms, and 3% for the cylindrical phantom for various field sizes, phantom thickness, and wedge angles. Conclusion: GPR tables are a ready reckoner for in-vivo dosimetry and it can be used to quickly estimate the MPD value from the EPID images with an accuracy of ±3% for common clinical treatment. project work funded by Union for International cancer control(UICC) under ICRETT fellowship.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Measurement of energy resolution with the NEXT-White silicon photomultipliers

The NEXT-White detector, a high-pressure gaseous xenon time projection chamber, demonstrated the excellence of this technology for future neutrinoless double beta decay searches using photomultiplier tubes (PMTs) to measure energy and silicon photomultipliers (SiPMs) to extract topology information. This analysis uses $^{83m}$Kr data from the NEXT-White detector to measure and understand the energy resolution that can be obtained with the SiPMs, rather than with PMTs. The energy resolution obtained of (10.9 ± 0.6)%, full-width half-maximum, is slightly larger than predicted based on the photon statistics resulting from very low light detection coverage of the SiPM plane in the NEXT-White detector. The difference in the predicted and measured resolution is attributed to poor corrections, which are expected to be improved with larger statistics. Furthermore, the noise of the SiPMs is shown to not be a dominant factor in the energy resolution and may be negligible when noise subtraction is applied appropriately, for high-energy events or larger SiPM coverage detectors. These results, which are extrapolated to estimate the response of large coverage SiPM planes, are promising for the development of future, SiPM-only, readout planes that can offer imaging and achieve similar energy resolution to that previously demonstrated with PMTs.[graphic not available: see fulltext]

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Aimsgb: An algorithm and open-source python library to generate periodic grain boundary structures

An algorithm implemented in an open-source python library was developed for building periodic grain boundary models in a universal fashion. The software framework, aimsgb, aims to generate tilt and twist grain boundaries from an input cubic or non-cubic crystal structure for ab-initio and classical atomistic simulation. It can output a coincidence site lattice (CSL) grain boundary for a cubic input structure and a non-CSL grain boundary for a non-cubic input structure. This framework has two useful features: (i) it can calculate all the CSL matrices for generating CSL from a given Sigma (Σ) value and rotation axis, allowing the users to build the specific CSL and grain boundary models; (ii) it provides a convenient command line tool to enable high-throughput generation of tilt and twist grain boundaries by assigning an input crystal structure, Σ value, rotation axis, and grain boundary plane. Here, the developed algorithm in the open-source python library is expected to facilitate studies of grain boundary in materials science. The software framework is available on the website: aimsgb.org.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Optimization of battery swapping infrastructure for e-commerce drone delivery

Drone delivery is widely-researched to alter the current e-commerce delivery convention for providing short delivery lead times. Yet, flight range of drones, constrained by the available battery technology, sets a milestone toward realizing the sole-drone delivery. To tackle with the flight range limitation, locating automated battery swapping machines (ABSM) have been proposed and a few studies modeled the problem. Using the ABSMs, drones can be loaded with fully-charged batteries along the route to a demand location. Here, we introduce a mixed-integer nonlinear program to model the problem. The objective of the program is to optimally select ABSM locations, determine the delivery-mode choices (drone-only, truck-only, and mixed delivery) of demand locations, find drone delivery routes, and approximate the baseline requirements for the number of drones and batteries needed. The program minimizes the overall delivery system costs including: ABSM, delivery, drone ownership, battery inventory, and service congestion. A cutting-plane method is developed to find exact solutions in finite iterations. Computational experiments showed that the method quickly yields the optimal solution to instances with less than 60 ABSM candidates and 20 demand locations. A case study shows that the optimal drone delivery infrastructure can save almost 20% cost compared to the conventional truck-only delivery. Sensitivity analyses were conducted to reveal the impact of key parameters in the decision-making and found that a decrease in the ABSM and drone costs highly affect the system cost.

25 ENERGY STORAGE↗