Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sparse data”

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 OSTI2018

Comparison of sparse domain approaches for 4D SPECT dynamic image reconstruction

Purpose: Dynamic imaging (DI) provides additional diagnostic information in emission tomography in comparison to conventional static imaging at the cost of being computationally more challenging. Dynamic single photon emission computed tomography (SPECT) reconstruction is particularly difficult because of the limitations in the sampling geometry present in most existing scanners. We have developed an algorithm Spline Initialized Factor Analysis of Dynamic Structures (SIFADS) that is a matrix factorization method for reconstructing the dynamics of tracers in tissues and blood directly from the projections in dynamic cardiac SPECT, without first resorting to any 3D reconstruction. Methods: SIFADS is different from "pure" factor analysis in dynamic structures (FADS) in that it employs a dedicated spline-based pre-initialization. In this paper, we analyze the convergence properties of SIFADS and FADS using multiple metrics. The performances of the two approaches are evaluated for numerically simulated data and for real dynamic SPECT data from canine and human subjects. Results: For SIFADS, metrics analyzed for reconstruction algorithm convergence show better features of the metric curves vs iterations. In addition, SIAFDS provides better tissue segmentations than that from pure FADS. Measured computational times are also typically better for SIFADS implementations than those with pure FADS. Conclusion: The analysis supports the utility of the pre-initialization of a factorization algorithm for better dynamic SPECT image reconstruction.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Knowledge Oriented Graph Unified Transformer (KOGUT) v0.1

KOGUT — Knowledge Oriented Graph Unified Transformer KOGUT implements the Relational Graph Transformer (RelGT) architecture for knowledge graph link prediction in biological domains, with a primary focus on microbial growth media prediction. While the original RelGT (arXiv:2505.10960) targets relational tables, time series, and multi-table databases, KOGUT adapts this architecture for heterogeneous biological knowledge graphs, providing first-in-class AI predictive models for microbial cultivation. Key Adaptations Beyond Original RelGT: - Knowledge Graph Focus: Applied to biological KGs with semantic node types (taxa, chemicals, media, phenotypes, environments) versus generic relational database tables, trained on the KG-Microbe knowledge graph (1.3M entities, 2.9M edges, 24 relation types). - Multimodal Node Encoding: Integrates node labels, categories, descriptions, and synonyms from KG metadata through learned embedding layers—adapting relational column features to graph node attributes with textual semantics. - Extended K-Hop Subgraph Strategy: Optimized neighborhood sampling (3-hop default, configurable up to 200 nodes) tuned for sparse biological networks, building on the original local-global attention framework with biological relation preservation. - Biolink Predicate Preservation: Type-specific transformations for 24 biological edge semantics (occurs_in, consumes, produces, has_phenotype, subclass_of) beyond standard relational foreign keys, enabling multi-relation link prediction. - Inductive Learning Support: Enables zero-shot predictions for novel taxa through feature-based embeddings (temperature, oxygen requirements, gram stain, cell shape), extending the original transductive relational benchmark scope to uncultured microorganisms. CheapSOTA Performance Optimizations (This Distribution): - VQ-EMA Centroid Attention: Vector quantization with exponential moving average for improved global context modeling (+5-10% MRR improvement). - HDF5 Precomputed Data Loading: One-time preprocessing of k-hop subgraphs to eliminate redundant graph traversals (2-5× training speedup). - Distributed Data Parallel Training: Multi-GPU support for scaling to larger knowledge graphs (tested on 4× NVIDIA A100 GPUs at NERSC Perlmutter). - Mixed Precision Training: Automatic mixed precision (AMP) for memory efficiency and faster training. Advantages Over Standard Knowledge Graph Embedding Models: Combines RelGT's proven multi-element tokenization (features, type, hop, structure) with graph-native biological representations, enabling interpretable link prediction across heterogeneous entities that standard embedding models (TransE, RotatE, ComplEx) and table-based transformers cannot directly model. Achieves near-perfect performance on microbial growth media prediction (MRR: 0.9966, Precision@1: 0.9932, Hit@10: 1.0000) while maintaining explainability through attention-based reasoning over biological pathways. Training Data: - KG-Microbe merged knowledge graph: 1,379,337 nodes, 2,960,472 edges - 24 biological relation types including taxonomic hierarchies, metabolic interactions, phenotype associations, and environmental relationships - Primary prediction task: Growth media suitability for microbial taxa (biolink:occurs_in, 50K edges) - Multi-relation capability: Predicts links for any of the 24 relation types, including chemical consumption/production, phenotype associations, and taxonomic classification Citation: Original RelGT Architecture: Dwivedi et al., "Relational Graph Transformer", arXiv:2505.10960, 2025 KOGUT Implementation: Knowledge Oriented Graph Unified Transformer for Microbial Growth Media Prediction Developed at Lawrence Berkeley National Laboratory (LBNL) Trained on NERSC Perlmutter supercomputer

Joachimiak, Marcin [Lawrence Berkeley National Lab↗
DOE OSTI2025

Evidence for a Single Holocene Paleoseismic Event on the Pajarito Fault, Northern New Mexico

Low-slip rate fault systems tend to be less studied than their high-slip rate counterparts, and paleoseismic techniques used to study them may pose challenges in interpretation that differ from high-slip rate systems. A good example of this is the Pajarito fault system (PFS), a normal fault complex within the Rio Grande rift. Despite numerous previous paleoseismic trenching studies conducted on the PFS between 1990 and 2003, considerable uncertainty remains regarding its Holocene paleoseismic history, particularly for the primary Pajarito fault (PF). To further clarify the PF paleoseismic history, we present data from paleoseismic investigations of 6 trenches at 3 distinct locations along the PF. Though the totality of the age and structural data obtained in this study is complex and not entirely consistent with any one interpretation, a single Holocene paleoearthquake occurring younger than ∼1,600 to 2,300 kcal yr BP is the simplest interpretation. It is possible that the PF records two Holocene events, with a penultimate event 6.9–2.4 kcal yr BP event and the aforementioned most recent event (MRE) between 2.3 and 1.6 kcal yr BP. However, only a single wall of one trench, out of a total of 12 walls in our 6 trenches, provides evidence supporting that interpretation. This study finds evidence of a single late Holocene paleoseismic event on the PF and sparse evidence for 2 Holocene paleoseismic events on the PF and highlights the benefits of logging multiple trench walls to better understand the complexity that results from this low-slip rate, low-deposition-rate fault system.

58 GEOSCIENCES↗
DOE OSTI2019

Aggressive Mammary Cancers Lacking Lymphocytic Infiltration Arise in Irradiated Mice and Can Be Prevented by Dietary Intervention

Because the incidence of breast cancer increases decades after ionizing radiation exposure, aging has been implicated in the evolution of the tumor microenvironment and tumor progression. In this paper, we investigated radiation-induced carcinogenesis using a model in which the mammary glands of 10-month-old BALB/c mice were transplanted with Trp53-null mammary tissue 3 days after exposure to low doses of sparsely ionizing γ-radiation or densely ionizing particle radiation. Mammary transplants in aged, irradiated hosts gave rise to significantly more tumors that grew more rapidly than those in sham-irradiated mice, with the most pronounced effects seen in mice irradiated with densely ionizing particle radiation. Tumor transcriptomes identified a characteristic immune signature of these aggressive cancers. Consistent with this, fast-growing tumors exhibited an immunosuppressive tumor microenvironment with few infiltrating lymphocytes, abundant immunosuppressive myeloid cells, and high COX-2 and TGFβ. Only irradiated hosts gave rise to tumors lacking cytotoxic CD8 + lymphocytes (defined here as immune desert), which also occurred in younger irradiated hosts. These data suggest that host irradiation may promote immunosuppression. To test this, young chimera mice were fed chow containing a honeybee-derived compound with anti-inflammatory and immunomodulatory properties, caffeic acid phenethyl ester (CAPE). CAPE prevented the detrimental effects of host irradiation on tumor growth rate, immune signature, and immunosuppression. These data indicated that low-dose radiation, particularly densely ionizing exposure of aged mice, promoted more aggressive cancers by suppressing antitumor immunity. Dietary intervention with a nontoxic immunomodulatory agent could prevent systemic effects of radiation that fuel carcinogenesis, supporting the potential of this strategy for cancer prevention.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

A search for a new gauge boson A'

In the Standard Model, gauge bosons mediate the strong, weak, and electromagnetic forces. New forces could have escaped detection only if their mediators are either heavier than order(TeV) or weakly coupled to charged matter. New vector bosons with small coupling {alpha}' arise naturally from a small kinetic mixing with the photon and have received considerable attention as an explanation of various dark matter related anomalies. Such particles can be produced in electron-nucleus fixed-target scattering and then decay to e + e -+ pairs. New light vector bosons and their associated forces are a common feature of Standard Model extensions, but existing constraints are remarkably sparse. The APEX experiment will search for a new vector boson A' with coupling α'/α fs > 6 × 10 -8 to electrons in the mass range 65MeV < mass A' < 550MeV. The experiment will study e + e - production off an electron beam incident on a high-Z target in Hall A at Jefferson Lab. The e - and e + will be detected in the High Resolution Spectrometers (HRSs). The invariant mass spectrum of the e + e - pairs will be scanned for a narrow resonance corresponding to the mass of the A'. A test run for the APEX experiment was held in the summer of 2010. Using the test run data, an A' search was performed in the mass range 175-250 MeV. The search found no evidence for an A' → e + e - reaction, and set an upper limit of {alpha}'/{alpha}{sub fs} ~ 10 -6 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Direct retrieval of isoprene from satellite-based infrared measurements

Isoprene is the atmosphere’s most important non-methane organic compound, with key impacts on atmospheric oxidation, ozone, and organic aerosols. In-situ isoprene measurements are sparse, and satellite-based constraints have employed an indirect approach using its oxidation product formaldehyde, which is affected by non-isoprene sources plus uncertainty and spatial smearing in the isoprene-formaldehyde relationship. Direct global isoprene measurements are therefore needed to better understand its sources, sinks, and atmospheric impacts. Here we show that the isoprene spectral signatures are detectable from space using the satellite-borne Cross-track Infrared Sounder (CrIS), develop a full-physics retrieval methodology for quantifying isoprene abundances from these spectral features, and apply the algorithm to CrIS measurements over Amazonia. The results are consistent with model output and in-situ data, and establish the feasibility of direct global space-based isoprene measurements. Finally, we demonstrate the potential for combining space-based measurements of isoprene and formaldehyde to constrain atmospheric oxidation over isoprene source regions.

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

A projection method for particle resampling

Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase-space due to their better scalability than continuum approaches with respect to dimension. Complex processes collectively referred to as particle noise hamper long time simulations with particle methods. One approach to address this problem is particle mesh adaptivity, or remapping, known as particle resampling and remeshing. Here, this work introduces a resampling method that projects particles to and from a (finite element) function space. The method is simple, using standard sparse linear algebra and finite element techniques, and it preserves all moments up to the order of a polynomial represented exactly by the continuum function space. It is distinguished from most other mesh-based methods in that new particle positions and number are decoupled from the mesh, allowing particle and continuum meshes to be adapted relatively independently. While this work is developed with structured particle and continuum phase-space grids on 1X + 1V Vlasov-Poisson models of Landau damping and two-stream instability, the method is well-suited to unstructured grids. Stable long time dynamics are demonstrated up to time T = 500. Reproducibility artifacts and data are publicly available.

Kinetic methods↗
DOE OSTI2025

Accurate and uncertainty-aware multi-task prediction of HEA properties using prior-guided deep Gaussian processes

Surrogate modeling techniques have become indispensable in accelerating the discovery and optimization of high-entropy alloys (HEAs), especially when integrating computational predictions with sparse experimental observations. This study systematically evaluates the training and testing performance of four prominent surrogate models—conventional Gaussian processes (cGP), Deep Gaussian processes (DGP), encoder-decoder neural networks for multi-output regression and eXtreme Gradient Boosting (XGBoost)—applied to a hybrid dataset of experimental and computational properties of the 8-component HEA system Al-Co-Cr-Cu-Fe-Mn-Ni-V. We specifically assess their capabilities in predicting correlated material properties, including yield strength, hardness, modulus, ultimate tensile strength, elongation, and average hardness under dynamic/quasi-static conditions, alongside auxiliary computational properties. The comparison highlights the strengths of hierarchical deep modeling approaches in handling heteroscedastic, heterotopic, and incomplete data commonly encountered in materials science. Our findings illustrate that combined surrogate models such as DGPs infused with machine-learned priors outperform other surrogates by effectively capturing inter-property correlations and by assimilating prior knowledge. This enhanced predictive accuracy positions the combined surrogate models as powerful tools for robust and data-efficient materials design.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sparse Bayesian mass mapping with uncertainties: hypothesis testing of structure

ABSTRACT A crucial aspect of mass mapping, via weak lensing, is quantification of the uncertainty introduced during the reconstruction process. Properly accounting for these errors has been largely ignored to date. We present a new method to reconstruct maximum a posteriori (MAP) convergence maps by formulating an unconstrained Bayesian inference problem with Laplace-type l1-norm sparsity-promoting priors, which we solve via convex optimization. Approaching mass mapping in this manner allows us to exploit recent developments in probability concentration theory to infer theoretically conservative uncertainties for our MAP reconstructions, without relying on assumptions of Gaussianity. For the first time, these methods allow us to perform hypothesis testing of structure, from which it is possible to distinguish between physical objects and artefacts of the reconstruction. Here, we present this new formalism, and demonstrate the method on simulations, before applying the developed formalism to two observational data sets of the Abell 520 cluster. Initial reconstructions of the Abell 520 catalogues reported the detection of an anomalous ‘dark core’ – an overdense region with no optical counterpart – which was taken to be evidence for self-interacting dark matter. In our Bayesian framework, it is found that neither Abell 520 data set can conclusively determine the physicality of such dark cores at $99{{\ \rm per\ cent}}$ confidence. However, in both cases the recovered MAP estimators are consistent with both sets of data.

Price, M. A.↗
DOE OSTI2016

MAPPING THE GAS TURBULENCE IN THE COMA CLUSTER: PREDICTIONS FOR ASTRO-H

Astro-H will be able for the first time to map gas velocities and detect turbulence in galaxy clusters. One of the best targets for turbulence studies is the Coma cluster, due to its proximity, absence of a cool core, and lack of a central active galactic nucleus. To determine what constraints Astro-H will be able to place on the Coma velocity field, we construct simulated maps of the projected gas velocity and compute the second-order structure function, an analog of the velocity power spectrum. We vary the injection scale, dissipation scale, slope, and normalization of the turbulent power spectrum, and apply measurement errors and finite sampling to the velocity field. We find that even with sparse coverage of the cluster, Astro-H will be able to measure the Mach number and the injection scale of the turbulent power spectrum—the quantities determining the energy flux down the turbulent cascade and the diffusion rate for everything that is advected by the gas (metals, cosmic rays, etc.). Astro-H will not be sensitive to the dissipation scale or the slope of the power spectrum in its inertial range, unless they are outside physically motivated intervals. We give the expected confidence intervals for the injection scale and the normalization of the power spectrum for a number of possible pointing configurations, combining the structure function and velocity dispersion data. Importantly, we also determine that measurement errors on the line shift will bias the velocity structure function upward, and show how to correct this bias.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Simulations of Fractal Star Cluster Formation. I. New Insights for Measuring Mass Segregation of Star Clusters with Substructure

We compare the existent methods, including the minimum spanning tree based method and the local stellar density based method, in measuring mass segregation of star clusters. We find that the minimum spanning tree method reflects more the compactness, which represents the global spatial distribution of massive stars, while the local stellar density method reflects more the crowdedness, which provides the local gravitational potential information. It is suggested to measure the local and the global mass segregation simultaneously. We also develop a hybrid method that takes both aspects into account. This hybrid method balances the local and the global mass segregation in the sense that the predominant one is either caused by dynamical evolution or purely accidental, especially when such information is unknown a priori. In addition, we test our prescriptions with numerical models and show the impact of binaries in estimating the mass segregation value. As an application, we use these methods on the Orion Nebula Cluster (ONC) observations and the Taurus cluster. We find that the ONC is significantly mass segregated down to the 20th most massive stars. In contrast, the massive stars of the Taurus cluster are sparsely distributed in many different subclusters, showing a low degree of compactness. The massive stars of Taurus are also found to be distributed in the high-density region of the subclusters, showing significant mass segregation at subcluster scales. Meanwhile, we also apply these methods to discuss the possible mechanisms of the dynamical evolution of the simulated substructured star clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Dynamic global vegetation models underestimate net CO 2 flux mean and inter-annual variability in dryland ecosystems

Despite their sparse vegetation, dryland regions exert a huge influence over global biogeochemical cycles because they cover more than 40% of the world surface (Schimel 2010 Science 327 418–9). It is thought that drylands dominate the inter-annual variability (IAV) and long-term trend in the global carbon (C) cycle (Poulter et al 2014 Nature 509 600–3, Ahlstrom et al 2015 Science 348 895–9, Zhang et al 2018 Glob. Change Biol. 24 3954–68). Projections of the global land C sink therefore rely on accurate representation of dryland C cycle processes; however, the dynamic global vegetation models (DGVMs) used in future projections have rarely been evaluated against dryland C flux data. Here, we carried out an evaluation of 14 DGVMs (TRENDY v7) against net ecosystem exchange (NEE) data from 12 dryland flux sites in the southwestern US encompassing a range of ecosystem types (forests, shrub- and grasslands). We find that all the models underestimate both mean annual C uptake/release as well as the magnitude of NEE IAV, suggesting that improvements in representing dryland regions may improve global C cycle projections. Across all models, the sensitivity and timing of ecosystem C uptake to plant available moisture was at fault. Spring biases in gross primary production (GPP) dominate the underestimate of mean annual NEE, whereas models' lack of GPP response to water availability in both spring and summer monsoon are responsible for inability to capture NEE IAV. Errors in GPP moisture sensitivity at high elevation forested sites were more prominent during the spring, while errors at the low elevation shrub and grass-dominated sites were more important during the monsoon. We propose a range of hypotheses for why model GPP does not respond sufficiently to changing water availability that can serve as a guide for future dryland DGVM developments. Our analysis suggests that improvements in modeling C cycle processes across more than a quarter of the Earth's land surface could be achieved by addressing the moisture sensitivity of dryland C uptake.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

AmeriFlux CA-SCC Scotty Creek Landscape

This is the AmeriFlux version of the carbon flux data for the site CA-SCC Scotty Creek Landscape. Site Description - The Scotty Creek flux tower is located in an organic-rich boreal forest-wetland landscape about 50 km south of Fort Simpson in the Taiga Plains of the Mackenzie watershed. The tower was installed in 2013 and operates an open-path EC system year-round running on solar power only. Flux footprints contain about 50 % forested peat plateaus and 50 % wetlands (i.e., collapse-scar bogs). The forests are underlain by permafrost, while the treeless wetlands are permafrost-free. The tower itself is located on a forested peat plateau. Black spruce tree density on plateaus is sparse and the mean canopy height is ca. 5 m.

Sonnentag, Oliver↗
DOE OSTI2020

AWB-GCN: A Graph Convolutional Network Accelerator with Runtime Workload Rebalancing

The recent development of deep learning has been mostly focusing on Euclidean data, such as images, videos, audios, etc. However, most real-world information and relation are often expressed as graphs. To efficiently learn from graph data, graph convolutional networks (GCNs) emerge as a promising approach, showing advantages in several practical applications such as social network analysis, knowledge discovery, 3D modeling, motion capturing, etc. Real-world graphs are usually extremely large and imbalanced, posting significant performance demand and design challenges on the hardware dedicated for GCN inference. In this paper, we propose an architecture design called UW-GCN to accelerate graph convolutional network inference. To tackle the major performance bottleneck from workload imbalance, we propose dynamic neighborhood stealing and remote chunk shuffling techniques, relying on hardware flexibility to achieve hardware auto-tuning under negligible area or delay overhead. Specifically, UW-GCN is able to smartly profile the sparse graph pattern while continuously adjusting the workload distribution via routing reconfiguration among parallel processing elements (PEs). The ideal configuration is then reused in the remaining iterations. To the best of our knowledge, this is the first accelerator design particularly for GCN and the first work relying on hardware auto-tuning, which is normally based on software, to achieve near-optimal workload balance in processing sparse structures.

Geng, Tong↗
DOE OSTI2024

Radiation induced athermal diffusivity in uranium mononitride

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuel considered for light water reactors and advanced reactor designs. Properties like self- and fission gas diffusivity need to be better understood, given that they influence key fuel performance phenomena such as fission gas swelling and release. In particular, the radiation induced athermal (D 3 ) diffusivity remains challenging to accurately predict and has only been sparsely characterized in UN, despite its importance as it likely governs diffusion at the low temperatures this high-thermal-conductivity fuel form may operate. Molecular Dynamics simulations are used to estimate the mean square displacement induced by a primary knock-on atom (PKA) with a given kinetic energy. These results are combined with the PKA energy distributions obtained from binary collision approximation calculations to obtain the displacement due to a particular fission fragment. Finally, this is combined with experimental fission fragment yields to determine the displacement due to an average fission event and, thus, express the athermal diffusivity as a function of the fission rate density. These results are in excellent agreement with available experimental data. In conclusion, a particular importance is given to the understanding and the quantification of the variability of these results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Bayesian differential programming for robust systems identification under uncertainty

This paper presents a machine learning framework for Bayesian systems identification from noisy, sparse and irregular observations of nonlinear dynamical systems. The proposed method takes advantage of recent developments in differentiable programming to propagate gradient information through ordinary differential equation solvers and perform Bayesian inference with respect to unknown model parameters using Hamiltonian Monte Carlo sampling. This allows an efficient inference of the posterior distributions over plausible models with quantified uncertainty, while the use of sparsity-promoting priors enables the discovery of interpretable and parsimonious representations for the underlying latent dynamics. A series of numerical studies is presented to demonstrate the effectiveness of the proposed methods, including nonlinear oscillators, predator–prey systems and examples from systems biology. Taken together, our findings put forth a flexible and robust workflow for data-driven model discovery under uncertainty. All codes and data accompanying this article are available at https://bit.ly/34FOJMj .

Science & Technology - Other Topics↗
DOE OSTI2025

AmeriFlux US-TLR Timberlake Observatory for Wetland Restoration (TOWeR)

This is the AmeriFlux version of the carbon flux data for the site US-TLR Timberlake Observatory for Wetland Restoration (TOWeR). Site Description - This tower is located in Timberlake Forest, a restored forested wetland located 4 km from Albermarle Sound on the North Carolina coast. The forest to the south of the tower had historically been ditched, drained and converted to agricultural land, and subsequently reforested with Black Gum and Cypress trees and naturally re-wetted in 2006. This southern region has a bottomland forest ecosystem which saw a large-scale succession event of pine trees following the reforestation efforts. The northern part of the forest had been ditched, but never drained or converted for other land uses. It is much wetter and more sparsely vegetated when compared to the southern region, and has a mixed bottomland forest and swamp ecosystem.

Rey-Sanchez, Camilo [North Carolina State Universi↗
DOE OSTI2021

Improving the performance of first- and last-mile mobility services through transit coordination, real-time demand prediction, advanced reservations, and trip prioritization

Socio-demographic trends and recent economic development patterns have resulted in travel behavior changes that call for more flexible and accessible public transit options. Because flexible transit services vary in scope, size, and service type, new data-informed methods are useful to optimize services based on the specific needs of local communities and riders. In this study, real-world demand and vehicle trajectory data were used to evaluate and optimize system performance for an existing first-mile–last-mile (FMLM) service in Robinson Township, PA. A general FMLM model for arbitrary demand and service supply was then developed to quantify system performance—both travel time costs and day-to-day reliability—for various operational polices considering spatio-temporal demand variation and transportation network dynamics. Heuristics were used for optimal real-time vehicle routing in sizable real-world networks accommodating various service types and scopes. In this case study, total user costs were reduced by 18.6% when rides were coordinated with mainline fixed-route transit. Predictive routing strategies were shown to marginally improve system performance under sparse and variable spatio-temporal demand. The case study also highlights potentially large travel time and user reliability improvements—reductions of 51% and 53.8%, respectively—when trip requests were made in advance of their desired pickup time. Finally, we show that travel time reliability can be improved for time-inflexible trips with trip prioritization without increasing total user costs. These results were stable to changes in demand density.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗