Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Comparative framework”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Examining Metal Identity and Proximity Effects on Acetylene Hydrogenation with Azolate-Based MOFs

Liquid organic hydrogen carriers (LOHCs) are an attractive fuel source due to their compatibility with existing transportation methods and ease of use. However, they suffer from sluggish (de)hydrogenation kinetics. One promising platform for developing next-generation catalysts is metal–organic frameworks (MOFs), which can enable systematic interrogation into the influence of metal identity and spatial arrangement. In this study, the effect of the coordination environment was investigated using Ni- and Co-based azolate MOFs: MFU-4l-OH (M x Zn 5–x (OH) 4 (BTDD) 3 ; x = 4 for M = Co and x = 3 for M = Ni, H 2 BTDD = bis(1H-1,2,3-triazolo[4,5-b][4′,5′-i])dibenzo[1,4]dioxin), composed of single-site nodes, and M(OH) 2 BBTA (M = Ni, Co; H 2 BBTA = 1H,5H-benzo(1,2-d:4,5-d’)bistriazole), composed of extended chain-type nodes. The catalysts were characterized by isotherms, powder X-ray diffraction (PXRD), scanning electron microscopy (SEM), inductively coupled plasma-optical emission spectroscopy (ICP-OES), and X-ray photoelectron spectroscopy (XPS) analysis. Acetylene hydrogenation activity under steady state conditions (150 °C, 1:1 C 2 H 2 :H 2 ) revealed higher turnover frequencies (TOFs) up 1.17 × 10 –2 mol C 2 H 2 mol −1 Ni min −1 and 1.98 × 10 –3 mol C 2 H 2 mol −1 Co min −1 for Ni-MFU-4l-OH and Co-MFU-4l-OH, respectively, compared to their BBTA analogues. However, the Co-based MOFs, particularly Co 2 (OH) 2 -BBTA, exhibited greater selectivity (up to 19%) for the fully hydrogenated ethane product. Isosteric heat of adsorption (Q st ) measurements for ethylene and ethane revealed that the BBTA framework had stronger binding to the products than MFU-4l. Furthermore, these findings demonstrate that metal identity and coordination environment may modulate acetylene hydrogenation performance, leading to design principles for tuning LOHC hydrogenation catalysts.

catalysts↗
DOE OSTI2025

Joint modelling of astrophysical systematics for cosmology with LSST cosmic shear

ABSTRACT We present a novel framework for jointly modelling the weak lensing source galaxy redshift distribution and the intrinsic alignment (IA) of galaxies through a shared luminosity function (LF). In the context of a Rubin Observatory’s Legacy Survey of Space and Time (LSST) Year 1 and Year 10 cosmic shear analysis, we show that our novel approach produces cosmological parameter constraints which are comparable to standard methods, while offering more physical insight into IA and selection effects. We clarify the relationship between individual parameters of a Schechter LF and the redshift distribution of a magnitude-limited sample, showing the consequences of marginalizing over these parameters when modelling IAs in standard cosmic shear analyses. We explore the impact of the shape of the LF on the cosmic shear data vector, and we outline the potential of this method to naturally model selection functions in redshift distribution estimation. Although this work focuses on LSST cosmic shear, the proposed joint modelling framework is broadly applicable to weak lensing surveys.

Šarčević, Nikolina (ORCID:0000000173016415)↗
DOE OSTI2023

Measuring the giant radio galaxy length distribution with the LoTSS

Many massive galaxies launch jets from the accretion disk of their central black hole, but only ~10 3 instances are known in which the associated outflows form giant radio galaxies (GRGs, or giants): luminous structures of megaparsec extent that consist of atomic nuclei, relativistic electrons, and magnetic fields. Large samples are imperative to understanding the enigmatic growth of giants, and recent systematic searches in homogeneous surveys constitute a promising development. For the first time, it is possible to perform meaningful precision statistics with GRG lengths, but a framework to do so is missing. We measured the intrinsic GRG length distribution by combining a novel statistical framework with a LOFAR Two-metre Sky Survey (LoTSS) sample of freshly discovered giants. In turn, this allowed us to answer an array of questions on giants. For example, we can now assess how rare a 5 Mpc giant is compared with one of 1 Mpc, and how much larger – given a projected length – the corresponding intrinsic length is expected to be. Notably, we can now also infer the GRG number density in the Local Universe. We assumed the intrinsic GRG length distribution to be Paretian (i.e. of power-law form) with tail index ξ, and predicted the observed distribution by modelling projection and selection effects. To infer ξ, we also systematically searched the LoTSS for hitherto unknown giants and compiled the largest catalogue of giants to date. We show that if intrinsic GRG lengths are Pareto distributed with index ξ, then projected GRG lengths are also Pareto distributed with index ξ. Selection effects induce curvature in the observed projected GRG length distribution: angular length selection flattens it towards the lower end, while surface brightness selection steepens it towards the higher end. We explicitly derived a GRG’s posterior over intrinsic lengths given its projected length, laying bare the ξ dependence. We also discovered 2060 giants within LoTSS DR2 pipeline products; our sample more than doubles the known population. Spectacular discoveries include the largest, second-largest, and fourth-largest GRG known (l p = 5.1 Mpc, l p = 5.0 Mpc, and l p = 4.8 Mpc), the largest GRG known hosted by a spiral galaxy (l p = 2.5 Mpc), and the largest secure GRG known beyond redshift 1 (l p = 3.9 Mpc). We increase the number of known giants whose angular length exceeds that of the Moon from 10 to 23; among the discoveries is the angularly largest known radio galaxy in the Northern Sky, which is also the angularly largest known GRG ($\phi$ = 2°). Combining theory and data, we determined that intrinsic GRG lengths are well described by a Pareto distribution, and measured the index ξ = –3.5 ± 0.5. This implies that, given its projected length, a GRG’s intrinsic length is expected to be just 15% larger. Finally, we determined the comoving number density of giants in the Local Universe to be n GRG = 5 ± 2(100 Mpc) –3 . We developed a practical mathematical framework that elucidates the statistics of giant radio galaxy lengths. Through a LoTSS search, we also discovered 2060 new giants. By combining both advances, we determined that intrinsic GRG lengths are well described by a Pareto distribution with index ξ = –3.5 ± 0.5, and that giants are truly rare in a cosmological sense: most clusters and filaments of the Cosmic Web are not currently home to a giant. Thus, our work yields new observational constraints for analytical models and simulations featuring radio galaxy growth.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Post synthetic amine functionalization of MOF-808 for CO 2 sorption via ligand exchange and Michael addition

Post synthetic modification of metal organic frameworks presents a viable route for amine functionalization which can significantly enhance CO 2 sorption capacity. We present a facile means of amine incorporation using limited synthetic steps and low-cost reagents that results in a high density of primary alkyl amines distributed through a zirconium-based metal organic framework (MOF). Both the MOF synthesis and the post synthetic modification take place under aqueous conditions and result in strongly bound molecular amines available for sorbate interaction throughout the MOF pores. Furthermore, this amine incorporation protocol results in a significantly increased CO 2 capacity compared with the unmodified MOF-808. Specifically, CO 2 isotherms collected at 298 K for the unmodified MOF-808 show uptakes of 0.06, 0.2, and 1.2 mmol/g at 4, 15, and 100 kPa, respectively, which can be compared with 0.3, 0.7, and 2.5 mmol/g for the glycine grafted MOF-808 and 0.5, 0.9, and 2.3 mmol/g for ethylenediamine grafted MOF-808.

Amine sorbent↗
DOE OSTI2022

Beyond carbon flux partitioning: Carbon allocation and nonstructural carbon dynamics inferred from continuous fluxes

Carbon (C) allocation and nonstructural carbon (NSC) dynamics play essential roles in plant growth and survival under stress and disturbance. However, quantitative understanding of these processes remains limited. Here, in this work, we propose a framework where we connect commonly measured carbon cycle components (eddy covariance fluxes of canopy CO 2 exchange, soil CO 2 efflux, and allometry-based biomass and net primary production) by a simple mass balance model to derive ecosystem-level NSC dynamics (NSC i ), C translocation (dC i ), and the biomass production efficiency (BPE i ) in above- and belowground plant (i = agp and bgp) compartments. We applied this framework to two long-term monitored loblolly pine (Pinus taeda) plantations of different ages in North Carolina and characterized the variations of NSC and allocation in years under normal and drought conditions. The results indicated that the young stand did not have net NSC flux at the annual scale, whereas the mature stand stored a near-constant proportion of new assimilates as NSC every year under normal conditions, which was comparable in magnitude to new structural growth. Roots consumed NSC in drought and stored a significant amount of NSC post drought. The above- and belowground dC i and BPE i varied more from year to year in the young stand and approached a relatively stable pattern in the mature stand. The belowground BPE bgp differed the most between the young and mature stands and was most responsive to drought. With the internal C dynamics quantified, this framework may also improve biomass production estimation, which reveals the variations resulting from droughts. Overall, these quantified ecosystem-scale dynamics were consistent with existing evidence from tree-based manipulative experiments and measurements and demonstrated that combining the continuous fluxes as proposed here can provide additional information about plant internal C dynamics. Given that it is based on broadly available flux data, the proposed framework is promising to improve the allocation algorithms in ecosystem C cycle models and offers new insights into observed variability in soil–plant–climate interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A variational phase-field model For ductile fracture with coalescence dissipation

A novel phase-field model for ductile fracture is presented. The model is developed within a consistent variational framework in the context of finite-deformation kinematics. A novel coalescence dissipation introduces a new coupling mechanism between plasticity and fracture by degrading the fracture toughness as the equivalent plastic strain increases. The proposed model is compared with a recent alternative where plasticity and fracture are strongly coupled. Several representative numerical examples motivate specific modeling choices. In particular, a linear crack geometric function provides an “unperturbed” ductile response prior to crack initiation, and Lorentz-type degradation functions ensure that the critical fracture strength remains independent of the phase-field regularization length. In addition, the response of the model is demonstrated to converge with a vanishing phase-field regularization length. Additionally, the model is then applied to calibrate and simulate a three-point bending experiment of an aluminum alloy specimen with a complex geometry. The effect of the proposed coalescence dissipation coupling on simulations of the experiment is first investigated in a two-dimensional plane strain setting. The calibrated model is then applied to a three-dimensional calculation, where the calculated load-deflection curves and the crack trajectory show excellent agreement with experimental observations. Finally, the model is applied to simulate crack nucleation and growth in a specimen from a recent Sandia Fracture Challenge.

42 ENGINEERING↗
DOE OSTI2023

A deep learning model for automatic analysis of cavities in irradiated materials

Transmission electron microscopy (TEM) is a commonly used technique in materials science for defect investigation. Quantitative analysis of defects is important for understanding the properties of a material, but manual analysis of TEM micrographs can be time-consuming and prone to error, especially when the defects have irregular shapes rather than spherical shapes. Many existing methods or deep learning models do not handle a wide range of sizes for the same object type within a single image. In this work, we present a framework that enables users to train an instance segmentation model called Mask R- CNN on any microstructure dataset, perform multi-detection on the same image at different scales, and obtain properties (e.g., size, area) of the objects based on the desired shape (e.g., circle, ellipse, rectangle). Additionally, we have developed a parallel detection module that uses multiple GPUs to increase the efficiency of the object detection process. We demonstrate the capabilities of our framework using a set of TEM images of cavities with different shapes, size distributions, and background contrasts. Finally, we show that the performance of our model in terms of density, size, and swelling of the cavities is comparable to the human average and that our model achieves the highest recall value compared to existing methods due to the use of image multi-rescaling.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Predictive process mapping for laser powder bed fusion: A review of existing analytical solutions

One of the main challenges in the laser powder bed fusion (LPBF) process is making dense and defect-free components. These porosity defects are dependent upon the melt pool geometry and the processing conditions. Power-velocity (PV) processing maps can aid in visualizing the effects of LPBF processing variables and mapping different defect regimes such as lack-of-fusion, under-melting, balling, and keyholing. This work presents an assessment of existing analytical equations and models that provide an estimate of the melt pool geometry as a function of material properties. The melt pool equations are then combined with defect criteria to provide a quick approximation of the PV processing maps for a variety of materials. Finally, the predictions of these processing maps are compared with experimental data from the literature. Here, the predictive processing maps can be computed quickly and can be coupled with dimensionless numbers and high-throughput (HT) experiments for validation. The present work provides a boundary framework for designing the optimal processing parameters for new metals and alloys based on existing analytical solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Half-cell electrode assessments of a crossover-tolerant direct methanol fuel cell with a platinum group metal-free cathode

In this work, a platinum group metal- free (PGM-free) Fe-N-C oxygen reduction reaction (ORR) catalyst derived from metal organic framework precursors (Fe-MOF) is evaluated at the cathode side of a direct methanol fuel cell (DMFC). A wide range of methanol concentrations were used, and the performance of the DMFC using air was evaluated and compared to a commercial Pt/C cathode. We conducted our tests with a custom, integrated H 2 reference electrode that separately provided the anode and cathode overpotentials. This allowed us to perform a Tafel analysis on both electrodes and understand the dominant polarizations involved at the PGM-free cathode and the PtRu/C anode. Overall, we achieved a high cell performance using a Fe-N-C based catalyst under liquid methanol and air operation with a peak power density of 111 mW/cm 2 , where the cathode presented a markedly larger contribution to the voltage loss versus the anode. Additionally, we used nanoscale X-ray computed tomography (nano-CT) to image the Fe-N-C cathode and investigate the ionomer and pore size distributions allowing us to identify areas of improvement for catalyst layer fabrication.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Bulk properties and multi-particle correlations in large and small systems

Charged particle production is calculated in a hybrid framework consisting of the IP-Glasma initial state, Music viscous relativistic fluid dynamics, and the UrQMD microscopic hadronic cascade. Using one set of parameters, we compute observables for a large variety of collision systems. Furthermore, we compare to experimental data in p+p, p+Pb, Xe+Xe, Au+Au and Pb+Pb collisions at various energies, and make predictions for potential O+O runs at RHIC and LHC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A physiologically based pharmacokinetic (PBPK) model to align dosimetry of the isobutyl metabolic series in rats and humans

Here, we developed a physiologically based pharmacokinetic (PBPK) model in rats and humans for the isobutyl metabolic series including isobutyl acetate, isobutanol, isobutyraldehyde, and isobutyric acid. Chemical manufactures routinely use these compounds as solvents, for chemical synthesis, as potential biofuels, de-icing fluids, and additives for food and/or fragrance in consumer products. Human exposure to isobutyl compounds can occur through inhalation or oral routes. We previously developed a PBPK model for the propyl metabolic series and utilized it as a framework to create the isobutyl PBPK model due to the chemical similarities between the two series. To support model development, we measured in vitro metabolism of isobutyl acetate in rat and human blood and liver S9 fractions. Compared to rats, humans demonstrated faster isobutyl acetate hydrolysis in liver S9 fractions, while the hydrolysis rates in blood were similar between the two species. We used concentrations of isobutyl compounds measured in air and blood from rats exposed to isobutyl acetate and isobutanol as well as other published data to further parameterize the model. Following exposure to either isobutyl acetate or isobutanol, we observed isobutanol concentrations highest among the isobutyl compounds in the blood of rats. Overall, the model accurately predicts measured time course concentrations of isobutyl acetate, isobutanol, and isobutyric acid in available data in rats and humans. Sensitivity analyses identified alveolar ventilation rates, isobutyl metabolism rates, and cardiac output as the most sensitive parameters affecting concentrations of isobutyl compounds in blood. The isobutyl PBPK model enables comparisons of internal dose metrics across various isobutyl compound exposures and species and allows for calculation of equivalent external exposures that result in the same dose metric. Regulators can employ this PBPK model to predict and align internal dose metrics of isobutyl compounds for risk assessment purposes.

2-methyl-1-propanol↗
DOE OSTI2025

Multi-modal effects on indirect noise induced by turbulent entropy fields

Planar entropy waves are commonly assumed for predicting indirect combustion noise. However, the non-planar and turbulent nature of flows found in most practical combustors challenges this assumption. In the present paper, we examine the indirect noise generated by non-planar and turbulent entropy fields through subsonic nozzles. Firstly, we introduce a new transfer function framework that accounts for the contribution of non-planar Fourier modes of the entropy field to the indirect noise spectra. When applied to a turbulent flow field, this method demonstrates a significant improvement in spectral predictions compared with a conventional approach that only considers the planar mode. Secondly, simulations show that non-planar Fourier modes become significant above a threshold frequency ƒ thresh , found in the mid- to high-frequency range. This contribution of non-planar modes is explained by two-dimensional shear effects that distort the entropy waves. A scaling relation that uses residence times along streamlines is developed for ƒ thresh , showing good agreement with simulation results. Finally, we show that the indirect noise from non-planar entropy modes found in aviation combustors can be significant at frequencies below 1 kHz, which might be relevant in situations of thermo-acoustic instabilities coupled to indirect noise.

42 ENGINEERING↗
DOE OSTI2024

Conical Intersection Accessibility Dictates Brightness in Red Fluorescent Proteins

Red fluorescent protein (RFP) variants are highly sought after for in-vivo imaging since longer wavelengths improve depth and contrast in fluorescence imaging. However, the lower energy emission wavelength usually correlates with a lower fluorescent quantum yield compared to their green emitting counterparts. To guide the rational design of bright variants, we have theoretically assessed two variants (mScarlet and mRouge) which are reported to have very different brightness. Using an α-CASSCF QM/MM framework (chromophore and all protein residues within 6 Å of it in the QM region, for a total of more than 450 QM atoms), we identify key points on the ground and first excited state potential energy surfaces. The brighter variant mScarlet has a rigid scaffold, and the chromophore stays largely planar on the ground state. The dimmer variant mRouge shows more flexibility and can accommodate a pre-twisted chromophore conformation which provides easier access to conical intersections. Notably, the main difference between the variants lies in the intersection seam regions, which appear largely inaccessible in mScarlet but partially accessible in mRouge. This observation is mainly related with changes in the cavity charge distribution, the hydrogen-bonding network involving the chromophore and a key ARG/THR mutation (which changes both charge and steric hindrance).

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

Neural network potentials with effective charge separation for non-equilibrium dynamics of ionic solids: a ZnO case study

Developing neural network potentials (NNPs) accurate under non-equilibrium dynamics is challenging, as such systems require extensive sampling beyond equilibrium phases. Here we construct high-fidelity NNPs for zinc oxide (ZnO), a polymorphic ionic solid, using density functional theory (DFT) reference data. To efficiently capture transitional configurations, we combine enhanced-sampling molecular dynamics with empirical potentials, data distillation, and pretraining on short-range atomic energies (A-Train), followed by transfer learning with DFT-relabeled datasets. This hierarchical approach improves transferability across polymorphs and stress states. We further introduce effective charge separation, treating long-range Coulombic terms analytically while short-range residual interactions are learned by the NNP. The optimal effective charges fall in the range 0.5–1.0 q e , consistent with dielectric-screened values derived from formal charges but distinct from Bader estimates. Motivated by this observation, we propose a simple data-driven protocol in which effective charges are optimized by comparing DFT reference energies with explicit Coulomb calculations, without additional NNP training. This strategy improves accuracy and transferability in DFT-level predictions of energies, forces, and stress. Together, these results provide a practical charge-selection framework for robust NNP development in ionic solids, enabling reliable simulation of polymorphic phase transformations and non-equilibrium dynamics.

Chemistry↗
DOE OSTI2020

Synthesis and structural characterization of Ca 12 Ga 14 O 33

Ca 12 Ga 14 O 33 was successfully synthesized using a wet chemistry technique to promote the homogenous mixing of the Ca and Ga cations. Rietveld refinements on X-ray and neutron powder diffraction data confirm that the compound is isostructural to Ca 12 Al 14 O 33 , however, with a significantly larger lattice parameter allowing for the cages that result from the framework arrangement to expand. In naturally occurring Ca 12 Al 14 O 33 , the mineral mayenite, these cages are occupied by O 2- anions, however, experimental studies exchanging the O 2- anions with other anions has led to a host of applications, depending on the caged anion. The functional nature of the structure, where framework distortions coupled with cage occupants, are correlated to electronic band structure and modifications to the framework could lead to interesting physical properties. The phase evolution was tracked using thermogravimetric analysis and high temperature X-ray diffraction and showed a lower formation temperature for the Ca 12 Ga 14 O 33 analogue compared to Ca 12 Al 14 O 33 synthesized using the same wet chemistry technique. Analyzing both X-ray and neutron powder diffraction using the Rietveld method with two different starting models results in one structural model, with one Ca position and the caged O on a 24d special position, being preferred.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phase transformation path in Aluminum under ramp compression; simulation and experimental study

Abstract We present a framework based on non-equilibrium molecular dynamics (NEMD) to reproduce the phase transformation event of Aluminum under ramp compression loading. The simulated stress-density response, virtual x-ray diffraction patterns, and structure analysis are compared against the previously observed experimental laser-driven ramp compression in-situ x-ray diffraction data. The NEMD simulations show the solid–solid phase transitions are consistent to experimental observations with a close-packed face-centered cubic (fcc) (111), hexagonal close-packed (hcp) structure (002), and body-centered cubic bcc (110) planes remaining parallel. The atomic-level analysis of NEMD simulations identifiy the exact phase transformation pathway happening via Bain transformation while the previous in situ x-ray diffraction data did not provide sufficient information for deducing the exact phase transformation path.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HPDR: High-Performance Portable Scientific Data Reduction Framework

The rapid growth in scientific data generation is outpacing advancements in computing systems necessary for efficient storage, transfer, and analysis, particularly in the context of exascale computing. With the deployment of first-generation exascale computing systems and next-generation experimental facilities, this gap is widening and necessitates effective data reduction techniques to manage enormous data volumes. Over the past decade, various data reduction methods, including lossless compression, error-controlled lossy compression, and data refactoring, have been developed to accelerate I/O in scientific workflows. Despite significant reductions in data volume, these methods introduce considerable computational overhead, which can become the new bottleneck in data processing. To mitigate this, GPU-accelerated data reduction algorithms have been introduced. However, challenges remain in their integration into exascale workflows, including limited portability across different GPU architectures, substantial memory transfer overhead, and reduced scalability on dense multi-GPU systems. To address these challenges, we propose HPDR, a high-performance and portable data reduction framework. HPDR is designed to enable the execution of state-of-the-art reduction algorithms across diverse processor architectures while reducing memory transfer overhead to 2.3 % of the original, resulting in up to 3.5× faster throughput compared to existing solutions. It also achieves up to 96% of the theoretical speedup in multi-GPU settings. In addition, evaluations on accelerating I/O operations at scale up to 1,024 nodes of the Frontier supercomputer demonstrate that HPDR can achieve up to 103 TB/s reduction throughput, providing up to 4× acceleration in parallel I/O performance compared to existing data reduction routines. This work highlights the potential of HPDR to significantly enhance data reduction efficiency in exascale computing environments.

Chen, Jieyang [University of Oregon]↗
DOE OSTI2026

A Hierarchical Optimization Method for Electric Vertical Takeoff and Landing Aircraft Network Design

Electric vertical takeoff and landing aircraft (eVTOLs) are expected to serve urban air mobility in a station-to-station configuration, which makes the optimal network design of eVTOL stations a critical question to explore. Existing approaches often face limitations, such as the inability to interact station locations with demand or difficulty in finding the optimal solution for large study regions. Here, this paper first proposes a mathematical model to generate optimal eVTOL station locations while considering associated potential eVTOL demand, and then proposes a heuristic algorithm, Hierarchical Optimization MEthod (HOME), to efficiently solve the model. With a case study of Southern California, HOME was compared to 1) directly solving the original integer linear programming-based network design problem, and 2) employing the widely used genetic algorithm. Results suggest that HOME can find optimal solutions with limited computational resources. The proposed framework powered by HOME provides a computationally efficient way to support urban air mobility planning.

97 MATHEMATICS AND COMPUTING↗