Search NASA⌕ Search

SEARCH · Search NASA

Results for “IMPACT PARAMETER”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Typical and Unusual Properties of Magnetic Clouds during the WIND Era

A list of 33 magnetic clouds as identified in WIND magnetic field and plasma data has been compiled. The intervals for these events are provided as part of NASA/GSFC, WIND-MFI's Website under the URL http://lepmfi.qsfc.nasa.gov/mfi/mag_cloud publ.html#table The period covered in this study is from early 1995 to November 1998 which primarily occurs in the quiet part of the solar cycle. A force free, cylindrically symmetric, magnetic field model has been applied to the field data in 1-hour averaged form for all of these events (except one small event where 10 min avg's were used) and the resulting fit-parameters examined. Each event was provided a semi-quantitatively determined quality factor (excellent, good or poor). A set of 28 good or better cases, spanning a surprisingly large range of values for its various properties, was used for further analysis. These properties are, for example, durations, attitudes, sizes, asymmetries, axial field strengths, speeds, and relative impact parameters. They will be displayed and analyzed, along with some related derived quantities, with emphasis on typical vs unusual properties and on the magnetic fields magnetic clouds' relationships to the Sun and to upstream interplanetary shocks, where possible. For example, it is remarkable how narrowly distributed the speeds of these clouds are, and the overall average speed (390 techniques km/s) is less than that normally quoted for the average solar wind speed (420 km/s) despite the fact that many of these clouds are d"drivers" of interplanetary shocks. On average, a cloud appears to be a little less symmetric when the spacecraft is able to pass close to the cloud's axis as compared to a farther out passage. The average longitude and latitude (in GSE) of the axes of the clouds are 85 degrees and 8 degrees, respectively, with standard deviations near 40 degrees. Also, the half=yearly averaged axial magnetic flux has approximately tripled. almost monotonically, from about 6 to 17 X 10(exp 29) Mx over the first 3.5 years of consideration, but with a large uncertainty on each of the half-year estimates, because of small sampling. If true,this finding implies an approximate tripling of the events' solar fluxes over this period as it goes into solar maximum.

Lepping, R. P.↗
DOE OSTI2024

Analytical model for the motion and interaction of two-dimensional active nematic defects

Here, we develop an approximate, analytical model for the velocity of defects in active nematics by combining recent results for the velocity of topological defects in nematic liquid crystals with the flow field generated from individual defects in active nematics. Importantly, our model takes into account the long-range interactions between defects that result from the flows they produce as well as the orientational coupling between defects inherent in nematics. Our work complements previous studies of active nematic defect motion by introducing a linear approximation that allows us to treat defect interactions as two-body interactions and incorporates the hydrodynamic screening length as a tuning parameter. We show that the model can analytically predict bound states between two +1/2 winding number defects, effective attraction between two –1/2 defects, and the scaling of a critical unbinding length between ±1/2 defects with activity. The model also gives predictions for the trajectories of defects, such as the scattering of +1/2 defects by –1/2 defects at a critical impact parameter that depends on activity. In the presence of circular confinement, the model predicts a braiding motion for three +1/2 defects that was recently seen in experiments, as well as stable and ergodic trajectories for four or more defects.

36 MATERIALS SCIENCE↗
NASA NTRS1997

The Age of the Surface of Venus

Impact craters on Venus appear to be uniformly and randomly scattered over a once, but no longer, geologically active planet. To first approximation, the planet shows a single surface of a single age. Here we use Monte Carlo cratering simulations to estimate the age of the surface of Venus. The simulations are based on the present populations of Earth-approaching asteroids, Jupiter-family, Halley-family, and long period comets; they use standard Schmidt-Housen crater scalings in the gravity regime; and they describe interaction with the atmosphere using a semi-analytic 'pancake' model that is calibrated to detailed numerical simulations of impactors striking Venus. The lunar and terrestrial cratering records are also simulated. Both of these records suffer from poor statistics. The Moon has few young large craters and fewer still whose ages are known, and the record is biased because small craters tend to look old and large craters tend to look young. The craters of the Earth provide the only reliable ages, but these craters are few, eroded, of uncertain diameter, and statistically incomplete. Together the three cratering records can be inverted to constrain the flux of impacting bodies, crater diameters given impact parameters, and the calibration of atmospheric interactions. The surface age of Venus that results is relatively young. Alternatively, we can use our best estimates for these three input parameters to derive a best estimate for the age of the surface of Venus. Our tentative conclusions are that comets are unimportant, that the lunar and terrestrial crater records are both subject to strong biases, that there is no strong evidence for an increasing cratering flux in recent years, and that that the nominal age of the surface of Venus is about 600 Ma, although the uncertainty is about a factor of two. The chief difference between our estimate and earlier, somewhat younger estimates is that we find that the venusian atmosphere is less permeable to impacting bodies than supposed by earlier studies. An older surface increases the likelihood that Venus is dead.

Zahnle, K. J.↗
DOE OSTI2021

Modeling airborne pathogen transport and transmission risks of SARS-CoV-2

An integrated modeling approach has been developed to better understand the relative impacts of different expiratory and environmental factors on airborne pathogen transport and transmission, motivated by the recent COVID-19 pandemic. Computational fluid dynamics (CFD) modeling was used to simulate spatial-temporal aerosol concentrations and quantified risks of exposure as a function of separation distance, exposure duration, environmental conditions (e.g., airflow/ventilation), and face coverings. The CFD results were combined with infectivity models to determine probability of infection, which is a function of the spatial-temporal aerosol concentrations, viral load, infectivity rate, viral viability, lung-deposition probability, and inhalation rate. Uncertainty distributions were determined for these parameters from the literature. Probabilistic analyses were performed to determine cumulative distributions of infection probabilities and to determine the most important parameters impacting transmission. This modeling approach has relevance to both pathogen and pollutant dispersion from expelled aerosol plumes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Discovery of a Damped Lyα Absorber Originating in a Spectacular Interacting Dwarf Galaxy Pair at z = 0.026

We present the discovery of neutral gas detected in both damped Lyα absorption (DLA) and H I 21 cm emission outside of the stellar body of a galaxy, the first such detection in the literature. A joint analysis between the Cosmic Ultraviolet Baryon Survey and the MeerKAT Absorption Line Survey reveals an H I bridge connecting two interacting dwarf galaxies (log (M star /M ⊙ ) = 8.5 ± 0.2) that host a z = 0.026 DLA with log[N(H I )/cm –2 ] = 20.60 ± 0.05 toward the QSO J2339–5523 (z QSO = 1.35). At impact parameters of d = 6 and 33 kpc, the dwarf galaxies have no companions more luminous than ≈0.05L * within at least Δv = ±300 km s –1 and d ≈ 350 kpc. The H I 21 cm emission is spatially coincident with the DLA at the 2σ–3σ level per spectral channel over several adjacent beams. However, H I 21 cm absorption is not detected against the radio-bright QSO; if the background UV and radio sources are spatially aligned, the gas is either warm or clumpy (with a spin temperature to covering factor ratio T s /f c > 1880 K). Observations with VLT-MUSE demonstrate that the α-element abundance of the ionized interstellar medium (ISM) is consistent with the DLA (≈10% solar), suggesting that the neutral gas envelope is perturbed ISM gas. This study showcases the impact of dwarf–dwarf interactions on the physical and chemical state of neutral gas outside of star-forming regions. In the SKA era, joint UV and H I 21 cm analyses will be critical for connecting the cosmic neutral gas content to galaxy environments.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2023

Additive Manufacturing of Oxide Dispersion Strengthened Multi Principle Element Alloys for Future Aerospace Applications

Oxide Dispersion Strengthened (ODS) materials have long been of interest for their high temperature applications, and additive manufacturing enables their manufacturing viability. The ODS multi-principle element alloy NiCoCr was prepared using powder metallurgy techniques, additively manufactured, and evaluated for its processingmicrostructure-property relationships. The high temperature foundations of nickel-base superalloys and ODS materials were combined with the manufacturing advantages of 3D printing and the chemical simplicity of NiCoCr to inspire this work, which was divided into powder and printed material assessments. The project was achieved through multiple iterative project loops to assess the processing parameters’ impact on the microstructure and mechanical properties of the feedstock powder and printed material. The powder investigations (Chapter 3) focused on understanding the oxide coating that formed on the metal powder following acoustic mixing. Time of Flight Secondary Ion Mass Spectrometry was used to semi-quantitatively assess the amount of yttrium on the surface of the mixed powders, and indicated that a combination of higher mixing condition energy and moderate mixing time resulted in the most oxide coating on the NiCoCr powder. The results were supported by a qualitative assessment of scanning electron images of coated powder particles. Following mixing, the ODS NiCoCr was consolidated by Laser Powder Bed Fusion. The evaluations of the printed material (Chapter 4) frst considered screening experiments including Archimedes’ density, porosity, and grain size and number metrics from electron backscatter diffraction data. After the ideal additive manufacturing parameters were identifed, both the oxide homogeneity and yield strength were discussed for the idealized printed material. Overall, the project suggests that the combined use of qualitative or semi-quantitative powder surface analysis with Archimedes’ density analyses can be a valid high-throughput technique which can lead to process optimization of Laser Powder Bed Fusion additively manufactured ODS material.

Laura G Wilson↗
DOE OSTI2020

Orbital Refinement and Stellar Properties for the HD 9446, HD 43691, and HD 179079 Planetary Systems

The Transit Ephemeris Refinement and Monitoring Survey is a project that aims to detect transits of intermediate-long period planets by refining orbital parameters of the known radial velocity planets using additional data from ground-based telescopes, calculating a revised transit ephemeris for the planet, then monitoring the planet host star during the predicted transit window. Here we present the results from three systems that had high probabilities of transiting planets: HD 9446 b and c, HD 43691 b, and HD 179079 b. We provide new radial velocity (RV) measurements that are then used to improve the orbital solution for the known planets. We search the RV data for indications of additional planets in orbit and find that HD 9446 shows a strong linear trend of 4.8σ. Using the newly refined planet orbital solutions, which include a new best-fit solution for the orbital period of HD 9446 c, and an improved transit ephemerides, we found no evidence of transiting planets in the photometry for each system. Transits of HD 9446 b can be ruled out completely and transits HD 9446 c and HD 43691 b can be ruled out for impact parameters up to b = 0.5778 and b = 0.898, respectively, due to gaps in the photometry. A transit of HD 179079 b cannot be ruled out, however, due to the relatively small size of this planet compared to the large star and thus low signal to noise. We determine properties of the three host stars through spectroscopic analysis and find through photometric analysis that HD 9446 exhibits periodic variability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Moments of axial-vector GPD from lattice QCD: quark helicity, orbital angular momentum, and spin-orbit correlation

In this work, we present a lattice QCD calculation of the Mellin moments of the twist-2 axial-vector generalized parton distribution (GPD), $\overset{\sim }{H}\left(x,\xi, t\right)$ , at zero skewness, ξ, with multiple values of the momentum transfer, t. Our analysis employs the short-distance factorization framework on ratio-scheme renormalized quasi-GPD matrix elements. The calculations are based on an N f = 2 + 1 + 1 twisted mass fermions ensemble with clover improvement, a lattice spacing of a = 0.093 fm, and a pion mass of m π = 260 MeV. We consider both the iso-vector and iso-scalar cases, utilizing next-to-leading-order perturbative matching while omitting the disconnected contributions and gluon mixing in the iso-scalar case. For the first time, we determine the Mellin moments of $\overset{\sim }{H}$ up to the fifth order. From these moments, we discuss the quark helicity and orbital angular momentum contributions to the nucleon spin, as well as the spin-orbit correlations of the quarks. Additionally, we perform a Fourier transform over the momentum transfer, which allows us to explore the spin structure in the impact-parameter space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Three-dimensional imaging of pion using lattice QCD: generalized parton distributions

In this work, we report a lattice calculation of x-dependent valence pion generalized parton distributions (GPDs) at zero skewness with multiple values of the momentum transfer −t. The calculations are based on an N f = 2 + 1 gauge ensemble of highly improved staggered quarks with Wilson-Clover valence fermion. The lattice spacing is 0.04 fm, and the pion valence mass is tuned to be 300 MeV. We determine the Lorentz-invariant amplitudes of the quasi-GPD matrix elements for both symmetric and asymmetric momenta transfers with similar values and show the equivalence of both frames. Then, focusing on the asymmetric frame, we utilize a hybrid scheme to renormalize the quasi-GPD matrix elements obtained from the lattice calculations. After the Fourier transforms, the quasi-GPDs are then matched to the light-cone GPDs within the framework of large momentum effective theory with improved matching, including the next-to-next-to-leading order perturbative corrections, and leading renormalon and renormalization group resummations. We also present the 3-dimensional image of the pion in impact-parameter space through the Fourier transform of the momentum transfer −t.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Gravitational causality and the self-stress of photons

We study causality in gravitational systems beyond the classical limit. Using on-shell methods, we consider the 1-loop corrections from charged particles to the photon energy-momentum tensor — the self-stress — that controls the quantum interaction between two on-shell photons and one off-shell graviton. The self-stress determines in turn the phase shift and time delay in the scattering of photons against a spectator particle of any spin in the eikonal regime. We show that the sign of the β-function associated to the running gauge coupling is related to the sign of time delay at small impact parameter. Our results show that, at first post-Minkowskian order, asymptotic causality, where the time delay experienced by any particle must be positive, is respected quantum mechanically. Contrasted with asymptotic causality, we explore a local notion of causality, where the time delay is longer than the one of gravitons, which is seemingly violated by quantum effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Test of lepton flavour universality in W-boson decays into electrons and τ-leptons using pp collisions at s = 13 TeV with the ATLAS detector

A measurement of the ratio of the branching fractions, Rτ/e = B(W → τν)/B(W → eν), is performed using a sample of W bosons originating from top-quark decays to final states containing τ-leptons or electrons. This measurement uses pp collisions at s$$ \sqrt{s} $$ = 13 TeV, collected by the ATLAS experiment at the Large Hadron Collider during Run 2, corresponding to an integrated luminosity of 140 fb−1. The W → τντ (with τ → eνeντ) and W → eνe decays are distinguished using the differences in the impact parameter distributions and transverse momentum spectra of the electrons. The measured ratio of branching fractions Rτ/e = 0.975 ± 0.012 (stat.) ± 0.020 (syst.), is consistent with the Standard Model assumption of lepton flavour universality in W-boson decays.

Aad, G↗
DOE OSTI2023

Dielectron production at midrapidity at low transverse momentum in peripheral and semi-peripheral Pb–Pb collisions at ${\sqrt{s}}_{\textrm{NN}}$ = 5.02 TeV

The first measurement of the e + e – pair production at low lepton pair transverse momentum (p T,ee ) and low invariant mass (m ee ) in non-central Pb–Pb collisions at ${\sqrt{s}}_{\textrm{NN}}$ = 5.02 TeV at the LHC is presented. The dielectron production is studied with the ALICE detector at midrapidity (|η e | < 0.8) as a function of invariant mass (0.4 ≤ m ee < 2.7 GeV/c 2 ) in the 50–70% and 70–90% centrality classes for p T,ee < 0.1 GeV/c, and as a function of p T,ee in three m ee intervals in the most peripheral Pb–Pb collisions. Below a p T,ee of 0.1 GeV/c, a clear excess of e + e – pairs is found compared to the expectations from known hadronic sources and predictions of thermal radiation from the medium. The m ee excess spectra are reproduced, within uncertainties, by different predictions of the photon–photon production of dielectrons, where the photons originate from the extremely strong electromagnetic fields generated by the highly Lorentz-contracted Pb nuclei. Lowest-order quantum electrodynamic (QED) calculations, as well as a model that takes into account the impact-parameter dependence of the average transverse momentum of the photons, also provide a good description of the p T,ee spectra. The measured $\sqrt{\left\langle {p}_{\textrm{T},\textrm{ee}}^2\right\rangle }$ of the excess p T,ee spectrum in peripheral Pb–Pb collisions is found to be comparable to the values observed previously at RHIC in a similar phase-space region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Search for higgsinos in compressed mass spectra using low-momentum tracks in pp collisions at s=13 TeV with the ATLAS detector

This paper presents two searches for the electroweak production of higgsinos with compressed mass spectra using 140 fb−1 of s=13$$ \sqrt{s}=13 $$ TeV proton-proton collision data collected by the ATLAS experiment at the Large Hadron Collider. Events are required to feature an energetic jet, large missing transverse momentum, and at least one low-momentum charged particle that serves as a candidate higgsino decay product. In the first search, targeting higgsino mass splittings in the range of 0.3–1 GeV, the higgsinos are expected to predominantly decay into pions that are identified as low-momentum charged particles with large transverse impact parameters due to the long higgsino lifetime (cτ ≈ ?(0.1–10 mm)), and neural networks are used to discriminate between signal and background processes. The second search targets larger mass splittings in the range of 1–3 GeV, where the higgsinos are expected to decay promptly into low-momentum leptons, one of which is identified by dedicated low-momentum electron or muon taggers based on neural networks utilising tracking and calorimeter information. No significant excess above the Standard Model prediction is observed in either search and the results are interpreted within simplified models, to set lower limits on the masses of the higgsino-like charginos and neutralinos. Together, these searches exclude chargino masses below 126 GeV at 95% confidence level for mass splittings between the chargino and lightest neutralino in the range of 0.3–2 GeV. This represents the first ATLAS constraints in a portion of this parameter space and surpasses the limits previously set by other experiments.

Aad, G↗
DOE OSTI2021

Sharp boundaries for the swampland

We reconsider the problem of bounding higher derivative couplings in consistent weakly coupled gravitational theories, starting from general assumptions about analyticity and Regge growth of the S-matrix. Higher derivative couplings are expected to be of order one in the units of the UV cutoff. Our approach justifies this expectation and allows to prove precise bounds on the order one coefficients. Our main tool are dispersive sum rules for the S-matrix. We overcome the difficulties presented by the graviton pole by measuring couplings at small impact parameter, rather than in the forward limit. We illustrate the method in theories containing a massless scalar coupled to gravity, and in theories with maximal supersymmetry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Measurement of the top quark pair production cross section in PbPb collisions at $\sqrt{s_{\mathrm{NN}}}=5.36$ TeV

The inclusive cross section for top quark pair ($\mathrm{t}\overline{\mathrm{t}}$) production in lead-lead (PbPb) collisions is reported for the first time at a center-of-mass energy per nucleon pair of 5.36 TeV. The analysis uses data corresponding to an integrated luminosity of 1.58 nb −1 collected by the CMS experiment at the CERN LHC in 2023. The $\mathrm{t}\overline{\mathrm{t}}$ production cross section, ${\sigma}_{\mathrm{t}\overline{\mathrm{t}}}={3.42}_{-0.51}^{+0.54}{\left(\mathrm{stat}\right)}_{-0.43}^{+0.50}\left(\mathrm{syst}\right)$ μb, is measured in dilepton final states using a fit to a multivariate discriminator that combines the decay electron and muon kinematic properties with the multiplicity of bottom quark jets. The result is consistent with perturbative quantum chromodynamics calculations at next-to-next-to-leading order (NNLO) accuracy employing several nuclear parton distribution functions. In addition, the Drell–Yan production cross section (σ DY ) for dilepton masses above 10 GeV and the ratio of $\mathrm{t}\overline{\mathrm{t}}$ to DY cross sections $\left({R}_{\mathrm{t}\overline{\mathrm{t}}/\mathrm{DY}}\right)$ are found to be compatible with the NNLO predictions. The observables ${\sigma}_{\mathrm{t}\overline{\mathrm{t}}}$, σ DY , and ${R}_{\mathrm{t}\overline{\mathrm{t}}/\mathrm{DY}}$ are measured separately for central and semicentral PbPb collisions to investigate for the first time the dependence of top quark production on the collision impact parameter.

Heavy Ion Experiments↗
DOE OSTI2025

Study of $\langle {p}_{\text{T}}\rangle$ and its higher moments, and extraction of the speed of sound in Pb-Pb collisions with ALICE

Ultrarelativistic heavy-ion collisions produce a state of hot and dense strongly interacting QCD matter called quark-gluon plasma (QGP). On an event-by-event basis, the volume of the QGP in ultracentral collisions is mostly constant, while its total entropy can vary significantly due to quantum fluctuations, leading to variations in the temperature of the system. Exploiting this unique feature of ultracentral collisions allows for the interpretation of the correlation of the mean transverse momentum ($\langle$p T $\rangle$) of produced charged hadrons and the number of charged hadrons as a measure for the speed of sound, c s . This speed is related to the rate at which compression waves travel in the QGP and is determined by fitting the relative increase in $\langle$p T $\rangle$ with respect to the relative change in the average charged-particle density ($\langle$dN ch /dη$\rangle$) measured at mid-rapidity. This study reports the event-average $\langle$p T $\rangle$ of charged particles as well as the variance, skewness, and kurtosis of the event-by-event transverse momentum per charged particle ([p T ]) distribution in ultracentral Pb-Pb collisions at a center-of-mass energy of 5.02 TeV per nucleon pair using the ALICE detector. Different centrality estimators based on charged-particle multiplicity or the transverse energy of the event are used to select ultracentral collisions. By ensuring a pseudorapidity gap between the region used to define the centrality and the region used to perform the measurement, the influence of biases and their potential effects on the rise of the mean transverse momentum is tested. The measured c$^{2}_{s}$ is found to strongly depend on the exploited centrality estimator and ranges between 0.1146±0.0028 (stat.)±0.0065 (syst.) and 0.4374±0.0006 (stat.)±0.0184 (syst.) in natural units. The self-normalized variance shows a steep decrease towards ultracentral collisions, while the self-normalized skewness variables show a maximum, followed by a fast decrease. These non-Gaussian features are understood in terms of the vanishing of the impact-parameter fluctuations contributing to the event-to-event [p T ] distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Quantifying spatiotemporal variability in occupant exposure to an indoor airborne contaminant with an uncertain source location

Well-mixed zone models are often employed to compute indoor air quality and occupant exposures. While effective, a potential downside to assuming instantaneous, perfect mixing is underpredicting exposures to high intermittent concentrations within a room. When such cases are of concern, more spatially resolved models, like computational-fluid dynamics methods, are used for some or all of the zones. But, these models have higher computational costs and require more input information. A preferred compromise would be to continue with a multi-zone modeling approach for all rooms, but with a better assessment of the spatial variability within a room. Here to do so, we present a quantitative method for estimating a room's spatiotemporal variability, based on influential room parameters. Our proposed method disaggregates variability into the variability in a room's average concentration, and the spatial variability within the room relative to that average. This enables a detailed assessment of how variability in particular room parameters impacts the uncertain occupant exposures. To demonstrate the utility of this method, we simulate contaminant dispersion for a variety of possible source locations. We compute breathing-zone exposure during the releasing (source is active) and decaying (source is removed) periods. Using CFD methods, we found after a 30 minutes release the average standard deviation in the spatial distribution of exposure was approximately 28% of the source average exposure, whereas variability in the different average exposures was lower, only 10% of the total average. We also find that although uncertainty in the source location leads to variability in the average magnitude of transient exposure, it does not have a particularly large influence on the spatial distribution during the decaying period, or on the average contaminant removal rate. By systematically characterizing a room's average concentration, its variability, and the spatial variability within the room important insights can be gained as to how much uncertainty is introduced into occupant exposure predictions by assuming a uniform in-room contaminant concentration. We discuss how the results of these characterizations can improve our understanding of the uncertainty in occupant exposures relative to well-mixed models.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Designing a drone delivery network with automated battery swapping machines

Drones are projected to alter last-mile delivery, but their short travel range is a concern. In this study, we propose a drone delivery network design using automated battery swapping machines (ABSMs) to extend ranges. The design minimizes the long-term delivery costs, including ABSM investment, drone ownership, and cost of the delivery time, and locates ABSMs to serve a set of customers. We build a mixed-integer nonlinear program that captures the nonlinear waiting time of drones at ABSMs. To solve the problem, we create an exact solution algorithm that finds the globally optimal solution using a derivative-supported cutting-plane method. To validate the applicability of our program, we conduct a case study on the Chicago Metropolitan area using cost data from leading ABSM manufacturer and geographical data from the planning and operations language for agent-based regional integrated simulation (more commonly known as POLARIS). A sensitivity analysis identifies that ABSM service times and costs are the key parameters impacting the long-term adoption of drone delivery.

25 ENERGY STORAGE↗