Search NASA⌕ Search

SEARCH · Search NASA

Results for “Distributed space 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 505 records · Page 28

Quantifying the impact of sub-grid surface wind variability on sea salt and dust emissions in CAM5

This paper evaluates the impact of sub-grid variability of surface wind on sea salt and dust emissions in the Community Atmosphere Model version 5 (CAM5). The basic strategy is to calculate emission fluxes multiple times, using different wind speed samples of a Weibull probability distribution derived from model-predicted grid-box mean quantities. In order to derive the Weibull distribution, the sub-grid standard deviation of surface wind speed is estimated by taking into account four mechanisms: turbulence under neutral and stable conditions, dry convective eddies, moist convective eddies over the ocean, and air motions induced by mesoscale systems and fine-scale topography over land. The contributions of turbulence and dry convective eddy are parameterized using schemes from the literature. Wind variabilities caused by moist convective eddies and fine-scale topography are estimated using empirical relationships derived from an operational weather analysis data set at 15 km resolution. The estimated sub-grid standard deviations of surface wind speed agree well with reference results derived from 1 year of global weather analysis at 15 km resolution and from two regional model simulations with 3 km grid spacing. The wind-distribution-based emission calculations are implemented in CAM5. In terms of computational cost, the increase in total simulation time turns out to be less than 3 %. Simulations at 2° resolution indicate that sub-grid wind variability has relatively small impacts (about 7 % increase) on the global annual mean emission of sea salt aerosols, but considerable influence on the emission of dust. Among the considered mechanisms, dry convective eddies and mesoscale flows associated with topography are major causes of dust emission enhancement. With all the four mechanisms included and without additional adjustment of uncertain parameters in the model, the simulated global and annual mean dust emission increase by about 50 % compared to the default model. By tuning the globally constant dust emission scale factor, the global annual mean dust emission, aerosol optical depth, and top-of-atmosphere radiative fluxes can be adjusted to the level of the default model, but the frequency distribution of dust emission changes, with more contribution from weaker wind events and less contribution from stronger wind events. Lastly, in Africa and Asia, the overall frequencies of occurrence of dust emissions increase, and the seasonal variations are enhanced, while the geographical patterns of the emission frequency show little change.

58 GEOSCIENCES↗

Fullrmc, a rigid body reverse monte carlo modeling package enabled with machine learning and artificial intelligence

Here, a new Reverse Monte Carlo (RMC) package ‘fullrmc’ for atomic or rigid body and molecular, amorphous or crystalline materials is presented. fullrmc main purpose is to provide a fully modular, fast and flexible software, thoroughly documented, complex molecules enabled, written in a modern programming language (python, cython ,C and C++ when performance is needed) and complying to modern programming practices. fullrmc approach in solving an atomic or molecular structure is different from existing RMC algorithms and software. In a nutshell, traditional RMC methods and software randomly adjust atom positions until the whole system has the greatest consistency with a set of experimental data. In contrast, fullrmc applies smart moves endorsed with reinforcement machine learning to groups of atoms. While fullrmc allows running traditional RMC modelling, the uniqueness of this approach resides in its ability to customize grouping atoms in any convenient way with no additional programming efforts and to apply smart and more physically meaningful moves to the defined groups of atoms. Also fullrmc provides a unique way with almost no additional computational cost to recur a group’s selection, allowing the system to go out of local minimas by refining a group’s position or exploring through and beyond not allowed positions and energy barriers the unrestricted three dimensional space around a group.

97 MATHEMATICS AND COMPUTING↗

Predicting the Galactic population of free-floating planets from realistic initial conditions

ABSTRACT We present the first prediction for the mass distribution function of Galactic free-floating planets (FFPs) that aims to accurately include the relative contributions of multiple formation pathways and stellar populations. We derive our predicted distribution from dedicated simulations of planet birth, growth, migration, and ejection around circumbinary systems and extend these results to also include the contributions from single and wide binary systems. Our resulting FFP mass distribution shows several distinct features, including a strong peak at ${\sim} 8\, {\rm M}_{\oplus }$ arising from the transition between pebble and gas accretion regimes and a trough at ${\sim} 1\, {\rm M}_{\oplus }$ due to the shift in the dominant ejection process from planet–planet scattering to ejection through interactions with stars in circumbinary systems. We find that interactions with the central binary in close circumbinary systems are likely the dominant progenitor for FFPs more massive than Earth, leading to a steep power-law dependence in mass that agrees well with existing observations. In contrast, we find planet–planet scattering events in single and wide binary systems likely produce the majority of planets at Mars mass and below, resulting in a shallower power-law dependence. Our results suggest that existing extrapolations into the sub-terrestrial mass range may significantly overestimate the true FFP abundance. The features we predict in the mass distribution of FFPs will be detectable by upcoming space-based microlensing surveys and, if observed, will provide key insight into the origins of FFPs and the environments in which they form.

Coleman, Gavin A. L. (ORCID:0000000151118963)↗

Quantifying Risk in an Uncertain Future: The Evolution of Resource Adequacy

As our power grids transition toward a decarbonized energy mix, ensuring reliability and provision of grid services remains paramount. The power system has always been heavily influenced by the weather - extreme temperatures determine the timing of peak demand, winter cold snaps can limit natural gas supply, gas turbine reliability and output are affected by ambient conditions, and hydro output varies seasonally and annually. However, as the grid increasingly relies on variable renewable energy (VRE), like wind and solar, the attention to reliability and weather conditions is increasingly important. The implications of changing reliability are large. The Electric Reliability Council of Texas (ERCOT) rolling blackouts from earlier this year impacted millions of people across the state and could be seen from space (Figure 1).

hydraulic turbines↗

Reconstruction of Six-Dimensional Phase Space

A phase space is a mathematical representation of all possible physical states of a system. Particle beams at Fermilab exist within a six-dimensional (6D) phase space defined by three positional components, (x, y, z) and three momentum components, (px, py, pz). To reconstruct this space implies taking measurement data from detectors and mapping out particle behavior using computational methods. The beam detectors, however, are only able to detect spatial distribution among the events of the beam, therefore being limited to positional data. Also, due to the vast number of events in a particle beam, it is extremely difficult to analyze and differentiate every single one’s behavior. However, with Machine Learning (ML), which can distinguish between patterns and map out particle behavior more efficiently. We first used the particle beam software, G4beamline, to simulate a 10,000-event muon beam, adjusting parameters such as initial momentum magnitude (p¬0) and virtual detector position. Using ten virtual detectors, we analyzed p0 values such that minimum 9,990 events were analyzed by every detector. We then input the data from these beam simulations to a C++ program, that randomly selects 100 events, and creates a 2D histogram based on spatial distribution, detector position, and event intensity. This process is repeated 100 times to create 100 histograms per p0 value. These images were then input to a modified ResNet18 Convolutional Neural Network (CNN) for training, and to predict p0 from some unseen set of histograms. The model was accurate when trained on momentum increments of 5 MeV/c and provided with denser training samples around highly variable test values. These results displayed machine learning being able to accurately predict p0 from being trained on different particle behaviors.

Shirlee, Jermain [Fermilab]↗

Scrambling and quantum chaos indicators from long-time properties of operator distributions

Scrambling is a key concept in the analysis of nonequilibrium properties of quantum many-body systems. Most studies focus on its characterization via out-of-time-ordered correlation (OTOC) functions, particularly through the early-time decay of the OTOC. However, scrambling is a complex process which involves operator spreading and operator entanglement, and a full characterization requires one to access more refined information on the operator dynamics at several timescales. In this work we analyze operator scrambling by expanding the target operator in a complete basis and studying the structure of the expansion coefficients treated as a coarse-grained probability distribution in the space of operators. Here we study different features of this distribution, such as its mean, variance, and participation ratio, for the Ising model with longitudinal and transverse fields, kicked collective spin models, and random circuit models. We show that the long-time properties of the operator distribution display common features across these cases and discuss how these properties can be used as a proxy for the onset of quantum chaos. Finally, we discuss the connection with OTOCs and analyze the cost of probing the operator distribution experimentally using these correlation functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Imaging Wavepackets in Real-Space & Disentangling Ultrafast Solvation Dynamics with High-Energy Ultrafast X-ray Scattering

The microscopic information of solute and solute-solvent motions can be measured using ultrafast diffuse x-ray scattering, and compared to molecular dynamics simulations, as demonstrated in studies of solvated molecular systems done in X-ray free-electron lasers. However, for the typical photon energies used in these experiments, the photon momentum transfer range was limited to lower Q values where the signals arising from different parts of the studied system overlap significantly. Recently, this limitation was greatly alleviated by extending the photon energy at the Linac Coherent Light-Source (LCLS), significantly increasing the accessible momentum transfer range. This improvement is transformative for ultrafast diffuse x-ray scattering measurements, enabling for the first time to disentangle the solute scattering difference signal, that persists to higher Q values, from the bulk solvent and solute-solvent cross-terms difference signals. In addition, it enhances the prediction ability of modeling and simulation methods such as the hybrid QM/MM approach. The extended Q range also opens the way to resolve in real-space details regarding coherent wavepackets motions beyond their center-of-mass positions. Here, we present the first results on high-energy (18keV) diffuse x-ray scattering in solution, demonstrating high fidelity time-resolved scattering and analysis of the photoexcited model photocatalysts PtPOP ([Pt2(POP)4]4-), and IrDimen ([Ir2(dimen)4]2+) in several solvents. These complexes provide ideal systems for demonstrating the ability of high-Q ultrafast scattering and QM/MM simulations to decompose ultrafast solvation dynamics into specific changes in the solute-solvent pair distribution function with a particular focus on how electronically excited states change the interaction between the solvent and photo-catalytically active metal sites. We introduce an approach to further utilize the high-energy capability and demonstrate a single-shot ultra-wide-angle X-ray Scattering modality using two perpendicular detectors, spanning a scattering angle range of more than 100 degrees, allowing to extend the accessible momentum transfer range up to Q~14 Å-1. We develop a model-free method to invert the x-ray scattering signals and enable the recovery of multiple pair-density motions that happen simultaneously, allowing us to directly measure in real-space nuclear wavepacket motions. . [1] Natan, Adi. "Real-Space Inversion and Super-Resolution of Ultrafast X-ray Scattering using Natural Scattering Kernels." arXiv preprint arXiv:2107.05576 (2021)

Natan, Adi↗

Fast and Scalable Sparse Triangular Solver for Multi-GPU Based HPC Architectures

Designing efficient and scalable sparse linear algebra kernels on modern multi-GPU based HPC systems is a daunting task due to significant irregular memory references and workload imbalance across the GPUs. This is particularly the case for \textit{Sparse Triangular Solver (SpTRSV)} which introduces additional two-dimensional computation dependencies among subsequent computation steps. Dependency information is exchanged and shared among GPUs, thus warrant for efficient memory allocation, data partitioning, and workload distribution as well as fine-grained communication and synchronization support. In this work, we demonstrate that directly adopting unified memory can adversely affect the performance of SpTRSV on multi-GPU architectures, despite linking via fast interconnect like NVLinks and NVSwitches. Alternatively, we employ the latest NVSHMEM technology based on Partitioned Global Address Space programming model to enable efficient fine-grained communication and drastic synchronization overhead reduction. Furthermore, to handle workload imbalance, we propose a malleable task-pool execution model which can further enhance the utilization of GPUs. By applying these techniques, our experiments on the NVIDIA multi-GPU supernode V100-DGX-1 and DGX-2 systems demonstrate that our design can achieve on average 3.53x (up to 9.86x) speedup on a DGX-1 system and 3.66x (up to 9.64x) speedup on a DGX-2 system with 4-GPUs over the Unified-Memory design. The comprehensive sensitivity and scalability studies also show that the proposed zero-copy SpTRSV is able to fully utilize the computing and communication resources of the multi-GPU system.

Xie, Chenhao↗

Deep kernel methods learn better: from cards to process optimization

Abstract The ability of deep learning methods to perform classification and regression tasks relies heavily on their capacity to uncover manifolds in high-dimensional data spaces and project them into low-dimensional representation spaces. In this study, we investigate the structure and character of the manifolds generated by classical variational autoencoder (VAE) approaches and deep kernel learning (DKL). In the former case, the structure of the latent space is determined by the properties of the input data alone, while in the latter, the latent manifold forms as a result of an active learning process that balances the data distribution and target functionalities. We show that DKL with active learning can produce a more compact and smooth latent space which is more conducive to optimization compared to previously reported methods, such as the VAE. We demonstrate this behavior using a simple cards dataset and extend it to the optimization of domain-generated trajectories in physical systems. Our findings suggest that latent manifolds constructed through active learning have a more beneficial structure for optimization problems, especially in feature-rich target-poor scenarios that are common in domain sciences, such as materials synthesis, energy storage, and molecular discovery. The Jupyter Notebooks that encapsulate the complete analysis accompany the article.

97 MATHEMATICS AND COMPUTING↗

Pore-Scale Study on the Positive Feedback Between Stress and Porosity Caused by Pressure Solution in Porous Media

Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.

CT scanning↗

Emerging Energy Market Analysis Initiative, Methodological Framework

Planning and operations of the electric power sector are undergoing radical changes. Climate change mitigation efforts have forced rapid changes to the technology mix. Technologies like wind and solar have experienced rapid growth, while investment in fossil sources has peaked or is declining. These foundational changes are forcing changes to energy systems. Demand-side adoption of electrified technologies, including electric vehicles, is changing load profiles and opening up new avenues for consumer participation in the power systems. The implications of an evolving power system pertain to more than environmental and technical dimensions. Changes to the generation mix and its consequent upstream and downstream impacts such as fuel production have significant and highly concentrated consequences on economies and employment. Shifts towards distributed (or decentralized) generating assets offer the potential to reshape economic and employment opportunities associated with the energy sector across space and socioeconomic groups. The Emerging Energy Market Analysis (EMA) initiative aims to identify sustainable, regionally acceptable, and high-value energy solutions that are secure and equitable. Unlike short-term, least-cost choices that can narrowly account for traditional options, EMA’s focus on emerging energy markets recognizes that new or adapted practices and technologies can alter the frontier of solutions and advance a community’s social, economic, and natural pathways. Such change requires a more comprehensive analysis of societal input, resources, capabilities, and infrastructure. These considerations lay the foundation for community decision-making models that are responsive to community values as well as the history and drivers. The result is a community-based decision and engagement model that will be valuable to decisionmakers and developers of advanced and emerging energy solutions, seeking a social license to operate prior to project development.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Grid Planning for Building Electrification: A Report by the Energy Systems Integration Group’s Grid Planning for Building Electrification Task Force

The increased electrification of buildings across the United States is being driven by technical advancements, cost reductions for some building technologies, consumer preferences, and policy goals for decarbonization. However, the effects of this load growth on the electric distribution system are often only a minor consideration in policymaking and long-term planning studies. The long lead time and useful life of power system equipment means that the decisions taken today are expected to support our society well into the 2060s and beyond. Distribution planning stakeholders can take steps today to establish a grid foundation that captures the new challenges presented by building electrification. This report focuses on building electrification for the residential and commercial sectors, focusing most heavily on space heating, as this end use stands to most strongly impact demand.1 Across all end uses, building electrification could require somewhere between 10% and 70% more electricity generation capacity than exists today depending on technology adoption and energy efficiency.

24 POWER TRANSMISSION AND DISTRIBUTION↗

THE SPECTRAL ENERGY DISTRIBUTION OF THE COLDEST KNOWN BROWN DWARF

WISE J085510.83–071442.5 (hereafter WISE 0855–0714) is the coldest known brown dwarf (∼250 K) and the fourth-closest known system to the Sun (2.2 pc). It has been previously detected only in the J band and two mid-IR bands. To better measure its spectral energy distribution (SED), we have performed deep imaging of WISE 0855–0714 in six optical and near-IR bands with Gemini Observatory, the Very Large Telescope, and the Hubble Space Telescope . Five of the bands show detections, although one detection is marginal (S/N ∼ 3). We also have obtained two epochs of images with the Spitzer Space Telescope for use in refining the parallax of the brown dwarf. By combining astrometry from this work and previous studies, we have derived a parallax of 0.449 ± 0.008″ (2.23 ± 0.04 pc). We have compared our photometry for WISE 0855–0714 to data for known Y dwarfs and to the predictions of three suites of models by Saumon et al. and Morley et al. that are defined by the presence or absence of clouds and nonequilibrium chemistry. Our estimates of Y − J and J − H for WISE 0855–0714 are redder than colors of other Y dwarfs, confirming a predicted reversal of near-IR colors to redder values at temperatures below 300–400 K. In color–magnitude diagrams, no single suite of models provides a clearly superior match to the sequence formed by WISE 0855–0714 and other Y dwarfs. Instead, the best-fitting model changes from one diagram to the next. Similarly, all of the models have substantial differences from the SED of WISE 0855–0714. As a result, we are currently unable to constrain the presence of clouds or nonequilibrium chemistry in its atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigation of reactivity differences in cylinder arrays using various fill geometries with constant mass

The factors affecting criticality safety are well understood by the nuclear engineering industry. Array configurations present multiple challenges to the industry due to the competing effects that are present in array systems - most notably interaction/spacing, moderation, and reflection. However, even the fissile mass geometry has been shown to complicate array systems. This attribute is also seen in filled UF{sub 6} cylinders, thus making it important to understand the reactivity in cylinder arrays. The process of filling UF{sub 6} cylinders results in non-uniform cylinder fill geometries. Due to the non-uniform fill, it becomes necessary to model a specific UF{sub 6} mass, in cylinder arrays, in different ways to ensure understanding of the reactivity of the system. Some of the ways to model the UF{sub 6} include: spreading the UF{sub 6} mass across the cylinder, all UF{sub 6} on the bottom of the cylinder, UF{sub 6} consolidated in the middle of the cylinder, and UF{sub 6} consolidated on the outside of the cylinder. Spreading the UF{sub 6} mass across the cylinder is accomplished by adjusting the UF{sub 6} density to achieve the given UF{sub 6} mass, whereas all the other scenarios use a fixed UF{sub 6} density with varying thicknesses to achieve the given UF{sub 6} mass. A different modeling approach, which helps to capture mixture effects, is to model the UF{sub 6} mass using the nominal UF{sub 6} density along with the void space for a given mass and distribute the mixture throughout the entire cylinder. This paper compares the reactivity of a cylinder array at a constant mass value using three different modeling techniques - completely filling the cylinder using a reduced UF{sub 6} density, filling the cylinder with nominal UF{sub 6} density and void space throughout the cylinder (same amount of void space as in the last technique), and all UF{sub 6} in the bottom of the cylinder (all void space on top of UF{sub 6}). The results demonstrate that completely filling the cylinder produces the most reactive UF{sub 6} cylinder array system. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development and Application of a Risk Analysis Toolkit for Plant Resources Optimization

This report summarizes the R&D activities of the Risk Informed Asset Management (RIAM) project during fiscal year 2020 (FY20). This project focuses on the development of methods designed to optimize plant operations (e.g., maintenance/replacement schedule, optimal maintenance posture) provided system/component health/cost data. This project development lives in cooperation with the Plant Health Management (PHM) project which focuses on the development of methods that integrate component health data and propagate such information at the system level to evaluate most relevant sources of risk. This year’s activities for the RIAM project focused on the continuation of schedule optimization algorithms developed in FY19. While in FY19 we focused on both deterministic and stochastic capital budgeting methods, in FY20 we moved forward by implementing two versions of schedule optimization methods. The first one reformulates the capital budgeting problem in a distributionally robust form which allows the user to rely on data directly rather than proposing a distribution from the data itself. The second version reformulates the capital budgeting explicitly using risk measures as variables to maximize/minimize. Lastly, we focused on the development of methods designed to identify the optimal maintenance posture based on the Pareto Frontier analysis. Rather than performing a tradeoff analysis (i.e., identify the absolute best posture), the Pareto Frontier analysis performs a trade space exploration approach (i.e., identify value and costs of several postures and have the analyst perform the task of imposing desired value and cost constraints). This is performed by identifying maintenance postures that maximize value (e.g., system availability) and minimize operational costs, i.e., the Pareto frontier in a value-cost trade space.

97 MATHEMATICS AND COMPUTING↗

Investigation of boundary-layer wind predictions during nocturnal low-level jet events using the Weather Research and Forecasting model

Simulations of two periods featuring three consecutive low level jet (LLJ) events in the US Upper Great Plains during the autumn of 2011 were conducted to explore the impacts of various setup configurations and physical process models on simulated flow parameters within the lowest 200 m above the surface, using the Weather Research and Forecasting (WRF) model. Sensitivities of simulated flow parameters to the horizontal and vertical grid spacing, planetary boundary layer (PBL) and land surface model (LSM) physics options, were assessed. Data from a Light Detection and Ranging (lidar) system, deployed to the Weather Forecast Improvement Project (WFIP; Finley et al. 2013) were used to evaluate the accuracy of simulated wind speed and direction at 80 m above the surface, as well as their vertical distributions between 120 and 40 m, covering the typical span of contemporary tall wind turbines. All of the simulations qualitatively captured the overall diurnal cycle of wind speed and stratification, producing LLJs during each overnight period, however large discrepancies occurred at certain times for each simulation in relation to the observations. 54-member ensembles encompassing changes of the above discussed configuration parameters displayed a wide range of simulated vertical distributions of wind speed and direction, and potential temperature, reflecting highly variable representations of stratification during the weakly stable overnight conditions. Root mean square error (RMSE) statistics show that different ensemble members performed better and worse in various simulated parameters at different times, with no clearly superior configuration . Simulations using a PBL parameterization designed specifically for the stable conditions investigated herein provided superior overall simulations of wind speed at 80 m, demonstrating the efficacy of targeting improvements of physical process models in areas of known deficiencies. However, the considerable magnitudes of the RMSE values of even the best performing simulations indicate ample opportunities for further improvements.

Mirocha, Jeff D.↗

CWRF performance at downscaling China climate characteristics

The performance of the regional Climate-Weather Research and Forecasting model (CWRF) for downscaling China climate characteristics is evaluated using a 1980–2015 simulation at 30 km grid spacing driven by the ECMWF Interim reanalysis (ERI). It is shown that CWRF outperforms the popular Regional Climate Modeling system (RegCM4.6) in key features including monsoon rain bands, diurnal temperature ranges, surface winds, interannual precipitation and temperature anomalies, humidity couplings, and 95th percentile daily precipitation. Even compared with ERI, which assimilates surface observations, CWRF better represents the geographic distributions of seasonal mean climate and extreme precipitation. These results indicate that CWRF may significantly enhance China climate modeling capabilities.

Meteorology & Atmospheric Sciences↗

Modeling information flow in a computer processor with a multi-stage queuing model

In this paper, we introduce a nonlinear stochastic model to describe the propagation of information inside a computer processor. In this model, a computational task is divided into stages, and information can flow from one stage to another. The model is formulated as a spatially-extended, continuous-time Markov chain where space represents different stages. This model is equivalent to a spatially-extended version of the M/M/s queue. The main modeling feature is the throttling function which describes the processor slowdown when the amount of information falls below a certain threshold. We derive the stationary distribution for this stochastic model and develop a closure for a deterministic ODE system that approximates the evolution of the mean and variance of the stochastic model. In conclusion, we demonstrate the validity of the closure with numerical simulations.

97 MATHEMATICS AND COMPUTING↗