Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random Systems”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Contrasting Community Assembly Forces Drive Microbial Structural and Potential Functional Responses to Precipitation in an Incipient Soil System

Microbial communities in incipient soil systems serve as the only biotic force shaping landscape evolution. However, the underlying ecological forces shaping microbial community structure and function are inadequately understood. We used amplicon sequencing to determine microbial taxonomic assembly and metagenome sequencing to evaluate microbial functional assembly in incipient basaltic soil subjected to precipitation. Community composition was stratified with soil depth in the pre-precipitation samples, with surficial communities maintaining their distinct structure and diversity after precipitation, while the deeper soil samples appeared to become more uniform. The structural community assembly remained deterministic in pre- and post-precipitation periods, with homogenous selection being dominant. Metagenome analysis revealed that carbon and nitrogen functional potential was assembled stochastically. Sub-populations putatively involved in the nitrogen cycle and carbon fixation experienced counteracting assembly pressures at the deepest depths, suggesting the communities may functionally assemble to respond to short-term environmental fluctuations and impact the landscape-scale response to perturbations. We propose that contrasting assembly forces impact microbial structure and potential function in an incipient landscape; in situ landscape characteristics (here homogenous parent material) drive community structure assembly, while short-term environmental fluctuations (here precipitation) shape environmental variations that are random in the soil depth profile and drive stochastic sub-population functional dynamics.

16S amplicon sequencing↗
DOE OSTI2020

Defect evolution in Ni and solid-solution alloys of NiFe and NiFeCoCr under ion irradiation at 16 and 300 K

Single-phase concentrated solid-solution alloys (SP-CSAs) have shown unique chemical complexity at the levels of electrons and atoms, and their defect evolution is expected to be different from conventional dilute alloys. Single crystals of Ni, NiFe and NiFeCoCr are chosen as model systems to understand the chemical complexity on defect formation and damage accumulation in SP-CSAs under ion irradiation. The high-quality crystals were irradiated at 16 and 300 K to different ion fluences, to form irradiated region with little to heavy damages. Additionally, the ion-induced damage was determined using Rutherford backscattering spectrometry technique along a channeling direction (RBS/C) and the level of lattice damage in irradiated Ni and SP-CSAs was quantified from Monte Carlo (MC) simulations. The results are interpreted using the Multi Step Damage Accumulation model to reveal material damage accumulation kinetics. Key findings of the study are that in case of room temperature irradiations the damage level measured for complex alloys at the highest irradiation fluence of 2 × 10 15 cm -2 (~3 dpa) is significantly higher than that obtained for pure nickel samples and suggest two-step damage accumulation process with a defect transformation taking place at a fluence of about 1.5 × 10 15 cm -2 . Furthermore, structural and damage kinetic differences clearly imply that, with increasing degree of chemical complexity and high solid-solution strengthening effects from Ni to NiFe and to NiFeCoCr, the enhanced lattice stiffness resists to randomization of atomic configurations and inhibits the growth of extended defects.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mean-field control for efficient mixing of energy loads

In this work, we pose an engineering challenge of controlling an ensemble of energy devices via coordinated, implementation-light, and randomized on/off switching as a problem in nonequilibrium statistical mechanics. We show that mean-field control with nonlinear feedback on the cumulative consumption, assumed available to the aggregator via direct physical measurements of the energy flow, allows the ensemble to recover from its use in the demand response regime, i.e., transition to a statistical steady state, significantly faster than in the case of the fixed feedback. Moreover when the nonlinearity is sufficiently strong, one observes the phenomenon of “super-relaxation,” where the total instantaneous energy consumption of the ensemble transitions to the steady state much faster than the underlying probability distribution of the devices over their state space, while also leaving almost no devices outside of the comfort zone.

25 ENERGY STORAGE↗
DOE OSTI2021

Redefining Resource Adequacy for Modern Power Systems: A Report of the Redefining Resource Adequacy Task Force

Today's rapidly increasing levels of wind, solar, storage, and load flexibility require the industry to rethink reliability planning and resource adequacy methods for modern power systems. Periods with a risk of shortfall often no longer coincide with peak demand - reliability risks are less about peak load and more about the daily setting of the sun, extended cloud cover, wind speeds, cold snaps, and heat waves. In addition, demand is increasingly flexible. Key resources are time-sensitive, as batteries need time to recharge and electricity customers can only be asked to provide demand response for just so long. And reliability failures are often correlated - with one another and with the weather. Two driving factors require the industry to reconsider its analytical approach for resource adequacy: (1) Chronological grid operations: The increasing importance of variable renewable resources (such as wind and solar) and of energy-limited resources (storage and demand response) make it essential to understand the full year of chronological operation of the grid. Specific attention must be paid to hourly, seasonal, and inter-annual resource variability. The sequence of the variability is key, as energy-limited resources such as batteries or demand response require either a preceding period or subsequent period of high production to be useful for grid reliability. (2) Correlated events: Historically, resource adequacy analysis focused on shortfalls caused by random, discrete mechanical failures of large generating units. In contrast, shortfalls today are often caused by multiple, correlated events caused by common weather patterns. Resource adequacy analysis must increasingly shift its focus to these correlated events. The redesign of resource adequacy methods will benefit from a set of guiding principles to better allow for sharing of insights and best practices, interregional resource coordination, and a smoother regulatory process for resource procurement. The objective of this report is to move this redesign forward. It provides an overview of key drivers changing the way resource adequacy needs to be evaluated, identifies shortcomings of conventional approaches, and outlines first principles for practitioners to consider as they adapt their approaches. The central message is: what got us here won't get us there.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Spatial Mapping of Riverbed Grain-Size Distribution Using Machine Learning

Recent alluvial sediments in riverbeds play a significant role in controlling hydrologic exchange flows (HEFs) in river systems. The alluvial layer is usually associated with strong heterogeneity in physical properties (e.g., permeability and hydraulic conductivity), which affects local HEFs and therefore biogeochemical processes. The spatial distribution of these physical properties needs to be determined to inform the numerical models used to reveal the realistic hydro-biogeochemical behaviors. Such information can be obtained based on the intrinsic link between sediment grain-size distribution and hydraulic properties where sediment texture information is available. However, grain-size measurements are usually spatially sparse and do not have adequate coverage and resolution, particularly for a relatively large domain such as the Hanford Reach of the Columbia River. In this paper, we adopted machine learning (ML) approaches for categorizing and mapping the spatial distributions of riverbed substrate grain size and filling in missing areas of substrate data using the ML models along the reach. Such ML models for substrate size mapping were trained at 13,372 locations using measured substrate sizes along with observed and simulated attributes, including bathymetric attributes (e.g., elevation, slope, and aspect ratio) from LIDAR and bathymetric surveys, and hydrodynamic properties (e.g., water depth, velocity, shear stress, and their statistical moments). An ensemble bagging-based ML technique, Random Forest, was adopted to identify the most influential factors as predictors to develop the predictive models with over-fitting issues addressed. The models were evaluated with respect to each individual substrate size class and the lumped group, and then used to generate the final substrate size maps covering all the grid cells in the numerical modeling domain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Thermalization of nuclear matter in heavy-ion collisions at Fermi energies

We analyze the time evolution of the kinetic properties of nuclear matter produced in heavy-ion collisions at Fermi energies. The collision system is simulated using Constrained Molecular Dynamics (CoMD) transport calculations whose output is the isospin, position, and momentum of the nucleons. Focusing on central 35 A·MeV 40 Ca+ 40 Ca collisions we utilize this information to extract localized momentum distributions in volume elements of 8 fm 3 and time steps of 5 fm/$c$. We then parameterize the single-particle momentum distributions with thermally motivated fit functions in the local rest frame of each cell. While the transverse-momentum distributions are well reproduced by thermal ones, the longitudinal ones carry a marked imprint of the initial nuclear motion which we capture by introducing a centroid motion into our fit functions. In particular, we find that Fermi distributions yield significantly better fits than Boltzmann ones, a consequence of the Pauli blocking implemented in CoMD. From the fits we extract the time dependence of the thermodynamic and collective properties of the excited nuclear medium. We find that the transverse temperature gradually rises to about 6 MeV, which is accompanied by a dissipation of the initial centroid motion of the incoming nuclei which vanishes at about 100 fm/c after initial impact. In conclusion, we are therefore able to track the transition of beam energy into random kinetic energy for nucleons, suggesting a three-dimensional equilibration of energy in the late stages of the collision.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Construction of Women’s All-Around Speed Skating Event Performance Prediction Model and Competition Strategy Analysis Based on Machine Learning Algorithms

Introduction Accurately predicting the competitive performance of elite athletes is an essential prerequisite for formulating competitive strategies. Women’s all-around speed skating event consists of four individual subevents, and the competition system is complex and challenging to make accurate predictions on their performance. Objective The present study aims to explore the feasibility and effectiveness of machine learning algorithms for predicting the performance of women’s all-around speed skating event and provide effective training and competition strategies. Methods The data, consisting of 16 seasons of world-class women’s all-around speed skating competition results, used in the present study came from the International Skating Union (ISU). According to the competition rules, distinct features are filtered using lasso regression, and a 5,000 m race model and a medal model are built using a fivefold cross-validation method. Results The results showed that the support vector machine model was the most stable among the 5,000 m race and the medal models, with the highest AUC (0.86, 0.81, respectively). Furthermore, 3,000 m points are the main characteristic factors that decide whether an athlete can qualify for the final. The 11th lap of the 5,000 m, the second lap of the 500 m, and the fourth lap of the 1,500 m are the main characteristic factors that affect the athlete’s ability to win medals. Conclusion Compared with logistic regression, random forest, K-nearest neighbor, naive Bayes, neural network, support vector machine is a more viable algorithm to establish the performance prediction model of women’s all-around speed skating event; excellent performance in the 3,000 m event can facilitate athletes to advance to the final, and athletes with outstanding performance in the 500 m event are more likely competitive for medals.

Liu, Meng↗
DOE OSTI2021

Lattice collective modes from a continuum model of magic-angle twisted bilayer graphene

Here, we show that the insulating states of magic-angle twisted bilayer graphene support a series of collective modes corresponding to local particle-hole excitations on triangular lattice sites. Our theory is based on a continuum model of the magic angle flat bands. When the system is insulating at moiré band filling ν = -3, our calculations show that the ground state supports seven low-energy modes that lie well below the charge gap throughout the moiré Brillouin zone, one of which couples strongly to THz photons. The low-energy collective modes are faithfully described by a model with a local SU(8) degree of freedom in each moiré unit cell that we identify as the direct product of spin, valley, and an orbital pseudospin. Apart from spin and valley-wave modes, the collective mode spectrum includes a low-energy intraflavor exciton mode associated with transitions between flat valence and conduction band orbitals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A Comparison of Three Types of Computer-Based Procedures: An Experiment Using the Rancor Microworld Simulator

The nuclear power industry has historically used paper-based procedures, but a shift towards computer-based procedures (CBPs) has the potential to reduce human errors, alleviate mental workload, and improve work performance. Twenty-seven participants were randomly assigned to one of three CBP types using the Rancor Microworld Simulator, and each performed two different types of operational scenarios (startup and loss of feedwater). The three CBP types varied in levels of digitalization. It was hypothesized that there would be less favorable impressions of the most basic procedures, and that these would demonstrate lower usability than types with greater digitalization. Overall, our predictions were partially supported with some interesting caveats, specifically with some performance benefits for CBPs that provided indicators but not embedded controls. We discuss our findings in terms of optimal levels of digitalization/automation within nuclear operations and suggest pathways for future directions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Building heights and urban canopy parameters for urban modeling

GLObal Building heights for Urban Studies (UT-GLOBUS) is a random forest model based framework that provides a level-of-detail-1 (LoD-1) building height dataset. The primary objective of UT-GLOBUS is not to precisely predict the height and footprint of individual buildings, but rather to offer a functional framework for generating building level information using open-source datasets for modeling applications. Specifically, UT-GLOBUS is tailored to meet the requirements of deriving urban canopy parameters (UCPs) for the multi-layer model within the Weather Research and Forecasting (WRF) model and building heights for the SOLWEIG and SUEWS model. Building-level data is accessible in vector file format (GeoPackage: .gpkg), which can be converted into raster file format (geoTIFF). The vector files employ the Universal Transverse Mercator (UTM) projection. The vector files are compatible with GIS platforms like QGIS and ArcGIS, and can be imported for analysis using programming languages such as Python. We are also providing UCPs required by the multi-layer urban model in the urban WRF in binary file format. Additionally, we provide the urban fractions calculated using ESA world cover dataset (https://esa-worldcover.org/en) for WRF model in binary file format. These binary files can be directly incorporated into the WRF pre-processing system (WPS).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

An increase in marine heatwaves without significant changes in surface ocean temperature variability

Marine heatwaves (MHWs)—extremely warm, persistent sea surface temperature (SST) anomalies causing substantial ecological and economic consequences—have increased worldwide in recent decades. Concurrent increases in global temperatures suggest that climate change impacted MHW occurrences, beyond random changes arising from natural internal variability. Moreover, the long-term SST warming trend was not constant but instead had more rapid warming in recent decades. Here we show that this nonlinear trend can—on its own—appear to increase SST variance and hence MHW frequency. Using a Linear Inverse Model to separate climate change contributions to SST means and internal variability, both in observations and CMIP6 historical simulations, we find that most MHW increases resulted from regional mean climate trends that alone increased the probability of SSTs exceeding a MHW threshold. Our results suggest the need to carefully attribute global warming-induced changes in climate extremes, which may not always reflect underlying changes in variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗
DOE OSTI2020

Modeling the role of local crystallographic correlations in microstructures of Ti-6Al-4V using a correlated structure visco-plastic self-consistent polycrystal plasticity formulation

This paper presents a multi-level crystal plasticity-based simulation framework for modeling mechanical response and microstructure evolution of Ti-6Al-4V with α-lath/lamellar microstructures. Here, the model is a correlated structure visco-plastic self-consistent (CS-VPSC) formulation linking three scales: a single crystals micro-scale, a lath/lamellar colony meso-scale, and a lath/lamellar aggregate macro-scale. A selected hardening law for the evolution of critical resolved shear stress per slip system used in CS-VPSC is phenomenological. However, it adjusts the resistances of basal and prismatic slip systems based on the geometry of slip transfer between adjacent lamellae. Consistent with experimental evidences, the resolved shear stress on the pyramidal slip planes is dependent not only on the stress in the direction of slip but also on the two orthogonal shear stress components and the three normal stress components (non-Schmid effects). Electron backscatter diffraction (EBSD) data in conjunction with a procedure relying on α→β phase transformation is used to construct paired variants of α-laths/lamellae satisfying their local crystallographic correlations. The procedure fits volume fractions of individual laths/lamellae with the experimental EBSD data and selects a distribution of habit planes between adjacent variants with respect to the loading direction. The simulation framework is applied to interpret the deformation behavior in tension and compression along two sample directions of Ti-6Al-4V fabricated via laser powder bed fusion. Moreover, the model is used to simulate texture evolution during rolling of the material to large strains. It is demonstrated that the model is capable of predicting plastic anisotropy/asymmetry and the concomitant texture evolution. While the model reveals a significant effect of habit plane inclination with respect to the loading direction on yield stress, the comparison of the data and model predictions shows that a random distribution of habit planes fits the flow response. It is further inferred that the tension-compression asymmetry arises from the non-Schmid effects.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Solar Influencer Next Door: Predicting Low-Income Solar Referrals and Leads

Increasing the adoption of solar among low-to-moderate income (LMI) households remains an important policy goal because of its promise to simultaneously reduce energy burden and support the just distribution of benefits of renewable energy. However, scaling LMI solar remains challenging due to affordability and access issues. Most existing LMI adoption has occurred under public-funded programs, highlighting the importance of increasing the cost-effectiveness of these programs at scale. We develop a new household-level data set on LMI solar lead acquisition, referrals, and adoption to understand the processes through which LMI solar uptake has occurred in California. Then, we develop models to predict two sub-mechanisms in the solar adoption process: whether an otherwise qualified lead becomes "lost" i.e. non-responsive to outreach and, for existing clients, whether they refer solar to others. For the program analyzed, participants received their solar system at no cost, which deemphasizes economic drivers of solar adoption and could differ from other program experiences. Both models substantially improved the accuracy of prediction relative to a baseline. Overall, we find that peer effects and solar economics are important to predicting referrals, and household demographic factors in lead loss prediction. Finally, we find that referrals are both the highest quality and largest source of LMI solar leads, providing a promising mechanism to expand LMI programs further.

customer acquisition costs↗
DOE OSTI2020

Data-driven Alignment of the Tracker for the Mu2e Experiment

An alignment of the Mu2e Tracker is investigated using the Mu2e Offline Software and Millepede-II. Simulated misalignments of Tracker planes were inserted between Monte-Carlo event generation and before reconstruction using randomly generated $(x,y,z)$-translations and rotations of the rigid-body planes in the detector model. Cosmic ray muons, with no magnetic field present in the Detector Solenoid, will follow a straight-line trajectory. A maximum-likelihood fit is used to reconstruct a track and determine its parameters. Tracks traversing at least four planes, with ten registered straw hits, and a maximum time residual of $20$ ns are selected to be used in an alignment sample. The cuts were chosen to discard tracks with low information content, such as those with high levels of multiple scattering or those with low numbers of measurements contributing to a single track fit. Millepede-II produces the resulting alignment constants according to a simultaneous least-squares fit of the alignment parameters and the track parameters for all input tracks. The improved alignment constants are used to repeat the process and recover an increased number of suitable tracks, further improving the estimate of the alignment constants in subsequent iterations. A performance study of the process was conducted on a misaligned Tracker in the case of translation-only, and translation plus rotational misalignments. In both cases, the correct alignment constants were recovered to within statistical error after three global fit iterations. Additionally, a continuous integration system was developed for the \texttt{Mu2e/Offline} GitHub repository, introducing a convenient automated testing process for all proposed changes (``pull requests"). The tests developed included a full build of the software, followed by tests of the validation routines, geometry overlap checking, and multi-threading mode. Additionally, code quality checks and low-statistics physics comparison checks a re also available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Multistage mixed precision iterative refinement

Abstract Low precision arithmetic, in particular half precision (16‐bit) floating point arithmetic, is now available in commercial hardware. Using lower precision can offer significant savings in computation and communication costs with proportional savings in energy. Motivated by this, there has been a renewed interest in mixed precision iterative refinement schemes for solving linear systems , and new variants of GMRES‐based iterative refinement have been developed. Each particular variant with a given combination of precisions leads to different condition number‐based constraints for convergence of the backward and forward errors, and each has different performance costs. The constraints for convergence given in the literature are, as an artifact of the analyses, often overly strict in practice, and thus could lead a user to select a more expensive variant when a less expensive one would have sufficed. In this work, we develop a multistage mixed precision iterative refinement solver which aims to combine existing mixed precision approaches to balance performance and accuracy and improve usability. For a user‐specified initial combination of precisions, the algorithm begins with the least expensive approach and convergence is monitored via inexpensive computations with quantities produced during the iteration. If slow convergence or divergence is detected using particular stopping criteria, the algorithm switches to use a more expensive, but more reliable variant. A novel aspect of our approach is that, unlike existing implementations, our algorithm first attempts to use “stronger” GMRES‐based solvers for the solution update before resorting to increasing the precision(s). In some scenarios, this can avoid the need to refactorize the matrix in higher precision. We perform extensive numerical experiments on a variety of random dense problems and problems from real applications which confirm the benefits of the multistage approach.

Oktay, Eda↗
DOE OSTI2022

Gas generation and intramolecular isotope study in laboratory pyrolysis of the Springfield coal from the Illinois Basin

Position-specific (PS) isotopes of propane have been proposed as a potential geochemical tool to decipher various geological processes (e.g., thermal cracking, biodegradation, H exchange) in natural reservoirs. The limited studies have been conducted on the PS isotopes of propane from the pyrolysis gases from marine shales, and natural gases sourced from lacustrine and marine kerogens, but little is known on gases produced from the humic kerogen. Here, this study investigated the PS δ 13 C of propane in the closed-system pyrolysis of the Springfield coal, Illinois Basin, Indiana, at 310 to 470 °C (Easy %Ro: 0.76 to 3.07). The C kinetic isotope effect (KIE) of CH 4 produced in both this study and previous low-temperature pyrolysis of the same coal indicates the cleavage of C—O bonds is the main generation pathway at the early kerogen cracking stage, followed by the breakdown of alkyl groups. At the wet-gas cracking stage, C 3 H 8 production from thermally stable compounds has a significant influence on the bulk and position-specific C KIE in the pyrolysis of marine Woodford kerogen and Springfield coal. According to the PS δ 13 C of propane, the central site is likely more enriched in 13 C and the δ 13 C of the terminal site is relatively heterogeneous within the propyl group attached to different functional groups of the gas-prone kerogen. Our findings based on the pyrolysis experiments and natural gas samples indicate thermal cracking and biodegradation appear to alter the δ 13 C cen values more significantly than the δ 13 C ter values of propane. The larger magnitude of ΔC c-t in the kMC simulations (Peterson et al., 2018) compared with those from the marine shale and coal possibly implies the non-random distributions of 13 C of propane precursors in the kerogens. As a new dimension of intramolecular isotopic information of propane, the PS δ 13 C values can contribute to fingerprinting the gas origins and identifying the various geological processes (e.g., kerogen cracking, wet-gas cracking, microbial activities) in sedimentary basins.

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

Cooperative Routing for an Air–Ground Vehicle Team—Exact Algorithm, Transformation Method, and Heuristics

This paper considers a cooperative vehicle routing problem for an intelligence, surveillance, and reconnaissance mission in the presence of communication constraints between the vehicles. The proposed framework uses a ground vehicle and an unmanned aerial vehicle (UAV) that travel cooperatively and visit a set of targets while satisfying the communication constraints. The problem is formulated as a mixed-integer linear program, and a branch-and-cut algorithm is developed to solve the problem to optimality. Furthermore, a transformation method and a heuristic are also developed for the problem. The effectiveness of all the algorithms is corroborated through extensive computational experiments on several randomly generated instances. Furthermore, this paper is motivated by an intelligence, surveillance, and reconnaissance mission involving a single unmanned aerial vehicle (UAV) and a ground vehicle, where the vehicles must coordinate their activity in the presence of communication constraints. The combination of a small UAV and a ground vehicle is an ideal platform for such missions, since small UAVs can fly at low altitudes and can avoid obstacles or threats that would be problematic for the ground vehicle alone. This paper addresses the coordinated routing problem involving these two vehicles and presents an algorithm to obtain an optimal solution for this problem, fast heuristics to obtain good feasible solutions, and also a transformation method to transform any instance of this cooperative vehicle routing problem to an instance of the one-in-a-set traveling salesman problem.

33 ADVANCED PROPULSION SYSTEMS↗