Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local regression”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Distributed Monitoring of the R(sup 2) Statistic for Linear Regression

The problem of monitoring a multivariate linear regression model is relevant in studying the evolving relationship between a set of input variables (features) and one or more dependent target variables. This problem becomes challenging for large scale data in a distributed computing environment when only a subset of instances is available at individual nodes and the local data changes frequently. Data centralization and periodic model recomputation can add high overhead to tasks like anomaly detection in such dynamic settings. Therefore, the goal is to develop techniques for monitoring and updating the model over the union of all nodes data in a communication-efficient fashion. Correctness guarantees on such techniques are also often highly desirable, especially in safety-critical application scenarios. In this paper we develop DReMo a distributed algorithm with very low resource overhead, for monitoring the quality of a regression model in terms of its coefficient of determination (R2 statistic). When the nodes collectively determine that R2 has dropped below a fixed threshold, the linear regression model is recomputed via a network-wide convergecast and the updated model is broadcast back to all nodes. We show empirically, using both synthetic and real data, that our proposed method is highly communication-efficient and scalable, and also provide theoretical guarantees on correctness.

Bhaduri, Kanishka↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2021

Contribution of Meteorological Downscaling to Skill and Precision of Seasonal Drought Forecasts

Research in meteorological prediction on sub-seasonal to seasonal (S2S) timescales has seen growth in recent years. Concurrent with this, demand for seasonal drought forecasting has risen. While there is obvious synergy between these fields, S2S meteorological forecasting has typically focused on low resolution global models, while the development of drought can be sensitive to the local expression of weather anomalies and their interaction with local surface properties and processes. This suggests that downscaling might play an important role in the application of meteorological S2S forecasts to skillful forecasting of drought. Here, we apply the Generalized Analog Regression Downscaling (GARD) algorithm to downscale meteorological hindcasts from the NASA Goddard Earth Observing System (GEOS) global S2S forecast system. Downscaled meteorological fields are then applied to drive offline simulations with the Catchment Land Surface Model (CLSM) to forecast United States Drought Monitor (USDM) style drought indicators derived from simulated surface hydrology variables. We compare the representation of drought in these downscaled hindcasts to hindcasts that are not downscaled, using the North American Land Data Assimilation System Phase 2 (NLDAS-2) dataset as an observational reference. We find that downscaling using GARD improves hindcasts of temperature and temperature anomalies, but the results for precipitation are mixed and generally small. Overall, GARD downscaling led to improved hindcast skill for total drought across the Contiguous United States (CONUS), and improvements were greatest for extreme (D3) and exceptional (D4) drought categories.

Ryan A Zamora↗
NASA NTRS1994

Magnetic spectral signatures in the Earth's magnetosheath and plasma depletion layer

Correlations between plasma properties and magnetic fluctuations in the sub-solar magnetosheath downstream of a quasi-perpendicular shock have been found and indicate that mirror and ion cyclotronlike fluctuations correlate with the magnetosheath proper and plasma depletion layer, respectively (Anderson and Fueselier, 1993). We explore the entire range of magnetic spectral signatures observed from the Active Magnetospheric Particle Tracer Explorers/Charge Composition Explorer (AMPTE/CCE)spacecraft in the magnetosheath downstream of a quasi-perpendicular shock. The magnetic spectral signatures typically progress from predominantly compressional fluctuations,delta B(sub parallel)/delta B perpendicular to approximately 3, with F/F (sub p) less than 0.2 (F and F (sub p) are the wave frequency and proton gyrofrequency, respectively) to predominantly transverse fluctuations, delta B(sub parallel)/delta B perpendicular to approximately 0.3, extending up to F(sub p). The compressional fluctuations are characterized by anticorrelation between the field magnitude and electron density, n(sub e), and by a small compressibility, C(sub e) identically equal to (delta n(sub e)/n(sub e)) (exp 2) (B/delta B(sub parallel)) (exp 2) approximately 0.13, indicative of mirror waves. The spectral characteristics of the transverse fluctuations are in agreement with predictions of linear Vlasov theory for the H(+) and He(2+) cyclotron modes. The power spectra and local plasma parameters are found to vary in concert: mirror waves occur for beta(s ub parallel p) (beta (sub parallel p) identically = 2 mu(sub zero) n(sub p) kT (sub parallel p) / B(exp 2) approximately = 2, A(sub p) indentically = T(sub perpendicular to p)/T(sub parallel p) - 1 approximately = 0.4, whereas cyclotron waves occur for beta (sub parallel p) approximately = 0.2 and A(sub p) approximately = 2. The transition from mirror to cyclotron modes is predicted by linear theory. The spectral characteristics overlap for intermediate plasma parameters. The plasma observations are described by A(sub p) = 0.85 beta(sub parallel P) (exp - 0.48) with a log regression coefficient of -0.74. This inverse A(sub p) - beta(sub parallel p) correlation corresponds closely to the isocontours of maximum ion anisotropy instability growth, gamma (sub m)/omega(sub p) = 0.01, for the mirror and cyclotron modes. The agreement of observed properties and predictions of local theory suggests that the spectral signatures reflect the local plasma environment and that the anisotropy instabilities regulate A(sub p). We suggest that the spectral characteristics may provide a useful basis for ordering observations in the magnetosheath and that the A(sub p) - beta(sub parallel p) inverse correlation may be used as a beta-dependent upper limit on the proton anisotropy to represent kinetic effects.

Anderson, Brian J.↗
DOE OSTI2026

Revealing EDL-driven reduction mechanisms in binary, ternary, and quaternary fluorinated electrolytes via an integrated MD–DFT–ML framework

Accurately predicting solid electrolyte interphase (SEI) formation requires explicitly resolving the electric double layer (EDL) structure, which deviates significantly from that of the bulk electrolyte. Although an established molecular dynamics (MD) and Density Functional Theory (DFT) framework can model SEI formation by evaluating reduction reactions of local clusters in the EDL, it suffers from a combinatorial computational bottleneck. To overcome this limitation, we introduce a machine-learning-accelerated simulation workflow (MD–DFT–ML), integrating a gradient-boosted regression model trained on EDL composition data to efficiently predict reduction potentials. We apply this framework to seven fluorinated electrolytes comprising fluorinated anions, a fluorinated ester solvent, two types of diluent (ion-solvating ester vs. non-solvating ether), and an FEC additive. The analysis shows that the EDL selectively accumulates cation-binding species; consequently, the non–cation-binding ether diluent rarely enters the EDL and makes minimal contributions to SEI formation. DFT calculations on statistically representative EDL clusters provide reduction potentials and fluorine-release pathways, while the ML model, which substantially reduces the DFT workload, predicts cluster reduction energies with a mean absolute error of 0.1 eV. The combined MD–DFT–ML approach also quantifies contributions from different sources to LiF formation in the SEI. This methodology establishes a generalizable route for multiscale modeling electrolyte and interphase design for next-generation electrochemical energy-storage systems.

DFT-MD-ML workflow↗
DOE OSTI2025

A kinetic-based regularization method for data science applications

We propose a physics-based regularization technique for function learning, inspired by statistical mechanics. By drawing an analogy between optimizing the parameters of an interpolator and minimizing the energy of a system, we introduce corrections that impose constraints on the lower-order moments of the data distribution. This minimizes the discrepancy between the discrete and continuum representations of the data, in turn allowing to access more favorable energy landscapes, thus improving the accuracy of the interpolator. Our approach improves performance in both interpolation and regression tasks, even in high-dimensional spaces. Unlike traditional methods, it does not require empirical parameter tuning, making it particularly effective for handling noisy data. We also show that thanks to its local nature, the method offers computational and memory efficiency advantages over Radial Basis Function interpolators, especially for large datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A Comfort Model Simplification for Tight Integration with Grid Service Optimizations

Localized generation, storage, and intelligent power electronics are increasingly being integrated into building. The excess storage, thermal, and generation capacity can be used to develop new markets for ancillary services such as peak reduction, voltage and frequency support, or load flattening. Widespread participation in these markets will make the power grid more resilient. However, widespread participation in these markets will require a simple and intuitive control over the excess capacity by the building occupants. Most importantly for adopting these technologies, participation in these markets should not require sacrificing occupant comfort. Many of the proposed control and optimization technologies for supplying the ancillary services do not explicitly incorporate occupant comfort into their models and when they do, generally utilize a rigid constraint on the indoor air temperature set a priori. This reduces flexibility in using excess thermal capacity. It also cannot take into account uncertainties in local factors that affect thermal comfort such as clothing insulation or relative humidity. In this paper we present and justify some assumptions that simplify the standard predictors of thermal comfort. We also develop a regression model of the thermal comfort that is linear with respect to the indoor air temperature. Next we develop the concept of thermal comfort variation that eliminates the dependence of the thermal comfort on uncertain thermal comfort factors. Then, we present a quadratic program utilizing the thermal comfort variation in both the objective function and as a constraint. Finally, we show the results of a time-of-use cost optimization that utilizes the simplified thermal comfort model.

Melin, Alexander↗
NASA NTRSDate not supplied

Assessing Flooding Vulnerability to Assist High Water Intervention and Urban Planning Programs in the Charles River Watershed with NASA DEVELOP

The Charles River watershed intersects 35 municipalities within the Boston Metropolitan Area and has a total population of 1.2 million, making it one of the most densely populated watersheds in New England. In recent years, the watershed has observed higher rates of flood inundation, mainly due to increased development, extreme precipitation events, and increased surface runoff. As the frequency of flooding events increases and a changing climate poses an ongoing threat to local communities, governments and organizations in Massachusetts are in need of accurate flood risk assessments. NASA DEVELOP partnered with the Charles River Watershed Association, the Town of Natick’s Office of Sustainability, and the Massachusetts Audubon Society to assess flood vulnerability and susceptibility in the watershed. The team used Landsat 5 Thematic Mapper, Landsat 8 Operational Land Imager, Sentinel-1 C-Band Synthetic Aperture Radar, and Sentinel-2 MultiSpectral Instrument to assess the feasibility of identifying the extent of past flood events using remote sensing. After identifying images that overlapped with reported flood events, the team concluded that it was not feasible to use Earth observation data to detect localized flooding in the time available for this study. Instead, the Federal Emergency Management Agency (FEMA) 100-year floodplain was used as a proxy for areas where flooding may occur. The team used statistical regression analysis and validation and supervised classification to develop a flood susceptibility map, incorporating several flood conditioning factors. The susceptibility maps were calibrated to various thresholds, including two that highlight hypothetical flooding under more liberal and more conservative planning scenarios. These were overlaid with demographic and socioeconomic data to create flood vulnerability maps. The team’s flood susceptibility maps showed an improvement in capturing known flood events over the FEMA 100-year and 500-year floodplain maps. These results will be improved with the addition of stormwater drainage mapping and precipitation data. Results can be used to fill in the gaps to help the stakeholders understand their communities’ vulnerability and susceptibility to flooding and improve their preparedness plans.

Trista Brophy↗
NASA NTRSDate not supplied

A Climate and History Case Study of 18th- and 19th-Century Multidecadal Droughts in East Africa Using a new Tree-Ring Drought Atlas

Historians and paleoclimatologists have both identified the decades spanning the late-18th and early-19th centuries in many East African regions as a period of prolonged and severe drought. A challenge that emerges from both the historical evidence and the paleolimnological data that have been primarily used to characterize this drought period, however, concerns the dating of events and the characterization of their spatial character; both oral traditions and lake-core proxies cannot typically offer hydroclimate estimates with seasonal or annual resolution, nor can they provide gridded reconstructions on the order of 1° latitude and longitude. The new East African Drought Atlas (EADA) incorporates far-field dendroclimatic records and several local East African tree-ring chronologies to provide estimates of the self-calibrating Palmer Drought Severity Index (scPDSI) on a 0.5° latitude-longitude grid. We present the Point-by-Point Regression technique that was used to create the EADA, the rationale for the employed input data network, and the calibration and validation skill of the derived hydroclimatic field reconstruction. We use the EADA to characterize the timing of two decadal-scale droughts, centered over the Uganda and Kenya regions, which occurred in the late 18th century and early in the 19th century and were separated by a shorter, but pronounced wet period. We characterize the spatial patterns of the associated events and the dynamical causes that that the patterns imply. The reconstructed droughts are also compared to historical and paleolimnological data to evaluate the robustness of the reconstruction and the degree of agreement across the different sources of information.

Droughts↗
NASA NTRSDate not supplied

A Climate and History Case Study of 18th- and 19th-Century Multidecadal Droughts in East Africa Using a new Tree-Ring Drought Atlas

Historians and paleoclimatologists have both identified the decades spanning the late-18th and early-19th centuries in many East African regions as a period of prolonged and severe drought. A challenge that emerges from both the historical evidence and the paleolimnological data that have been primarily used to characterize this drought period, however, concerns the dating of events and the characterization of their spatial character; both oral traditions and lake-core proxies cannot typically offer hydroclimate estimates with seasonal or annual resolution, nor can they provide gridded reconstructions on the order of 1° latitude and longitude. The new East African Drought Atlas (EADA) incorporates far-field dendroclimatic records and several local East African tree-ring chronologies to provide estimates of the self-calibrating Palmer Drought Severity Index (scPDSI) on a 0.5° latitude-longitude grid. We present the Point-by-Point Regression technique that was used to create the EADA, the rationale for the employed input data network, and the calibration and validation skill of the derived hydroclimatic field reconstruction. We use the EADA to characterize the timing of two decadal-scale droughts, centered over the Uganda and Kenya regions, which occurred in the late 18th century and early in the 19th century and were separated by a shorter, but pronounced wet period. We characterize the spatial patterns of the associated events and the dynamical causes that that the patterns imply. The reconstructed droughts are also compared to historical and paleolimnological data to evaluate the robustness of the reconstruction and the degree of agreement across the different sources of information.

Droughts↗
DOE OSTI2025

LC-Opt: Benchmarking Reinforcement Learning and Agentic AI for End-to-End Liquid Cooling Optimization in Data Centers

Liquid cooling is critical for thermal management in high-density data centers with the rising AI workloads. However, machine learning-based controllers are essential to unlock greater energy efficiency and reliability, promoting sustainability. We present LC-Opt, a Sustainable Liquid Cooling (LC) benchmark environment, for reinforcement learning (RL) control strategies in energy-efficient liquid cooling of high-performance computing (HPC) systems. Built on the baseline of a high-fidelity digital twin of Oak Ridge National Lab's Frontier Supercomputer cooling system, LC-Opt provides detailed Modelica-based end-to-end models spanning site-level cooling towers to data center cabinets and server blade groups. RL agents optimize critical thermal controls like liquid supply temperature, flow rate, and granular valve actuation at the IT cabinet level, as well as cooling tower (CT) setpoints through a Gymnasium interface, with dynamic changes in workloads. This environment creates a multi-objective real-time optimization challenge balancing local thermal regulation and global energy efficiency, and also supports additional components like a heat recovery unit (HRU). We benchmark centralized and decentralized multi-agent RL approaches, demonstrate policy distillation into decision and regression trees for interpretable control, and explore LLM-based methods that explain control actions in natural language through an agentic mesh architecture designed to foster user trust and simplify system management. LC-Opt democratizes access to detailed, customizable liquid cooling models, enabling the ML community, operators, and vendors to develop sustainable data center liquid cooling control solutions.

Naug, Avisek [Hewlett Packard Enterprise]↗
DOE OSTI2022

MICCO: An Enhanced Multi-GPU Scheduling Framework for Many-Body Correlation Functions

Calculation of many-body correlation functions is one of the critical kernels utilized in many scientific computing areas, especially in Lattice Quantum Chromodynamics (Lattice QCD). It is formalized as a sum of a large number of contraction terms each of which can be represented by a graph consisting of vertices describing quarks inside a hadron node and edges designating quark propagations at specific time intervals. Due to its computation- and memory-intensive nature, real-world physics systems (e.g., multi-meson or multi-baryon systems) explored by Lattice QCD prefer to leverage multi-GPUs. Different from general graph processing, many-body correlation function calculations show two specific features: a large number of computation-/data-intensive kernels and frequently repeated appearances of original and intermediate data. The former results in expensive memory operations such as tensor movements and evictions. The latter offers data reuse opportunities to mitigate the data-intensive nature of many-body correlation function calculations. However, existing graph-based multi-GPU schedulers cannot capture these data-centric features, thus resulting in a sub-optimal performance for many-body correlation function calculations. To address this issue, this paper presents a multi-GPU scheduling framework, MICCO, to accelerate contractions for correlation functions particularly by taking the data dimension (e.g., data reuse and data eviction) into account. This work first performs a comprehensive study on the interplay of data reuse and load balance, and designs two new concepts: local reuse pattern and reuse bound to study the opportunity of achieving the optimal trade-off between them. Based on this study, MICCO proposes a heuristic scheduling algorithm and a machine-learning-based regression model to generate the optimal setting of reuse bounds. Specifically, MICCO is integrated into a real-world Lattice QCD system, Redstar, for the first time running on multiple GPUs. The evaluation demonstrates MICCO outperforms other state-of-art works, achieving up to 2.25× speedup in synthesized datasets, and 1.49× speedup in real-world correlation functions.

Wang, Qihan↗
DOE OSTI2021

Transcutaneous and End-Tidal CO 2 Measurements in Hypoxia and Hyperoxia

BACKGROUND: Transcutaneous measurement of carbon dioxide (CO 2 ) has been proposed for physiological monitoring of tactical jet aircrew because in some clinical settings it mirrors arterial CO 2 partial pressure (P a co 2 ). End-tidal monitoring in laboratory settings is known to give high-fidelity estimates of P a co 2 . METHODS: The correspondence between end-tidal (P ET co 2 ) and transcutaneous Pco 2 (tcPco 2 ) was examined in healthy volunteers under laboratory conditions of hyperoxia and hypoxia. Rest and exercise, skin heating and cooling, hyperventilation, and induced CO 2 retention were employed. RESULTS: Neither measure followed all known changes in P a co 2 and tcPco 2 changed when the skin temperature near the probe changed. Bland-Altman analysis showed significant nonzero slopes under most conditions. Regression analysis indicated that oxygen partial pressure (Po 2 ) in tissue measured as transcutaneous Po 2 (tcPo 2 ) is an important explanatory variable for tcPco 2 in addition to P ET co 2 , and that local skin temperature also has an effect. Additionally, absorption atelectasis from breathing 100% O 2 may cause P ET co 2 to deviate from P a co 2 . DISCUSSION: Even as a trend indicator for P a co 2 , tcPco 2 is not useful under conditions that resemble those in the highly dynamic tactical jet aircraft environment. P ET co 2 is also not a good indicator of CO 2 status in pilots who breathe nearly 100% O 2 . Shykoff BE, Lee LR, Gallo M, Griswold CA. Transcutaneous and end-tidal CO 2 measurements in hypoxia and hyperoxia . Aerosp Med Hum Perform. 2021; 92(11):864–872 .

Biophysics↗
NASA NTRS1992

Active and passive microwave measurements of soil moisture in FIFE

During the intensive field campaigns of the First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE) in May-October of 1987, several nearly simultaneous measurements were made with low-altitude flights of the L-band radiometer and C- and X-band scatterometers over two transects in the Konza Prairie Natural Research Area, some 8 km south of Manhattan, Kansas. These measurements showed that although the scatterometers were sensitive to soil moisture variations in most regions under the flight path, the L-band radiometer lost most of its sensitivity in regions unburned for many years. The correlation coefficient derived from the regression between the radar backscattering coefficient and the soil moisture was found to improve with the increase in antenna incidence angle. This is attributed to a steeper falloff of the backscattering coefficient as a function of local incidence at angles near nadir than at angles greater than 30 deg.

Wang, J. R.↗
NASA NTRS1995

Estimating the urban bias of surface shelter temperatures using upper-air and satellite data. Part 1: Development of models predicting surface shelter temperatures

Multiple regression techniques were used to predict surface shelter temperatures based on the time period 1986-89 using upper-air data from the European Centre for Medium-Range Weather Forecasts (ECMWF) to represent the background climate and site-specific data to represent the local landscape. Global monthly mean temperature models were developed using data from over 5000 stations available in the Global Historical Climate Network (GHCN). Monthly maximum, mean, and minimum temperature models for the United States were also developed using data from over 1000 stations available in the U.S. Cooperative (COOP) Network and comparative monthly mean temperature models were developed using over 1150 U.S. stations in the GHCN. Three-, six-, and full-variable models were developed for comparative purposes. Inferences about the variables selected for the various models were easier for the GHCN models, which displayed month-to-month consistency in which variables were selected, than for the COOP models, which were assigned a different list of variables for nearly every month. These and other results suggest that global calibration is preferred because data from the global spectrum of physical processes that control surface temperatures are incorporated in a global model. All of the models that were developed in this study validated relatively well, especially the global models. Recalibration of the models with validation data resulted in only slightly poorer regression statistics, indicating that the calibration list of variables was valid. Predictions using data from the validation dataset in the calibrated equation were better for the GHCN models, and the globally calibrated GHCN models generally provided better U.S. predictions than the U.S.-calibrated COOP models. Overall, the GHCN and COOP models explained approximately 64%-95% of the total variance of surface shelter temperatures, depending on the month and the number of model variables. In addition, root-mean-square errors (rmse's) were over 3 C for GHCN models and over 2 C for COOP models for winter months, and near 2 C for GHCN models and near 1.5 C for COOP models for summer months.

Epperson, David L.↗
NASA NTRS1986

The effect of local advection on the inference of soil moisture from thermal infrared radiances

The degree to which the GOES-VISSR infrared data can be used to infer area-averaged soil moisture is explored for a five-day case study period. Chosen variables are transformed and incorporated into a multiple linear regression. The actual observations, rather than a simplified model, are used to determine the relationship between soil moisture and GOES-IR radiance. It is shown that a depletion coefficient of 0.92 produces an index of ground truth which is best correlated with soil moisture as inferred from GOES thermal infrared data. When all individual daily soil estimates during the case study period are averaged at each point and compared to the average observed soil moisture, the data correlate at 0.85. This implies that the algorithm can distinguish at least four categories of soil moisture.

Wetzel, Peter J.↗
DOE OSTI2022

What drives the scatter of local star-forming galaxies in the BPT diagrams? A Machine Learning based analysis

ABSTRACT We investigate which physical properties are most predictive of the position of local star forming galaxies on the BPT diagrams, by means of different Machine Learning (ML) algorithms. Exploiting the large statistics from the Sloan Digital Sky Survey (SDSS), we define a framework in which the deviation of star-forming galaxies from their median sequence can be described in terms of the relative variations in a variety of observational parameters. We train artificial neural networks (ANN) and random forest (RF) trees to predict whether galaxies are offset above or below the sequence (via classification), and to estimate the exact magnitude of the offset itself (via regression). We find, with high significance, that parameters primarily associated to variations in the nitrogen-over-oxygen abundance ratio (N/O) are the most predictive for the [N ii]-BPT diagram, whereas properties related to star formation (like variations in SFR or EW(H α)) perform better in the [S ii]-BPT diagram. We interpret the former as a reflection of the N/O–O/H relationship for local galaxies, while the latter as primarily tracing the variation in the effective size of the S+ emitting region, which directly impacts the [S ii] emission lines. This analysis paves the way to assess to what extent the physics shaping local BPT diagrams is also responsible for the offsets seen in high redshift galaxies or, instead, whether a different framework or even different mechanisms need to be invoked.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Partner with a Third-Party Delivery Service or Not? A Prediction-and-Decision Tool for Restaurants Facing Takeout Demand Surges During a Pandemic

Amidst the COVID-19 pandemic, restaurants become more reliant on no-contact pick-up or delivery ways for serving customers. As a result, they need to make tactical planning decisions such as whether to partner with online platforms, to form their own delivery team, or both. In this paper, we develop an integrated prediction-decision model to analyze the profit of combining the two approaches and to decide the needed number of drivers under stochastic demand. We first use the susceptible-infected-recovered (SIR) model to forecast future infected cases in a given region and then construct an autoregressive-moving-average (ARMA) regression model to predict food-ordering demand. Using predicted demand samples, we formulate a stochastic integer program to optimize food delivery plans. We conduct numerical studies using COVID-19 data and food-ordering demand data collected from local restaurants in Nuevo Leon, Mexico, from April to October 2020, to show results for helping restaurants build contingency plans under rapid market changes. Our method can be used under unexpected demand surges, various infection/vaccination status, and demand patterns. Here, our results show that a restaurant can benefit from partnering with third-party delivery platforms when (i) the subscription fee is low, (ii) customers can flexibly decide whether to order from platforms or from restaurants directly, (iii) customers require more efficient delivery, (iv) average delivery distance is long, or (v) demand variance is high.

97 MATHEMATICS AND COMPUTING↗