Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “load profile”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

COMPASO : A new halo finder for competitive assignment to spherical overdensities

We describe a new method (compaso) for identifying groups of particles in cosmological N-body simulations. compaso builds upon existing spherical overdensity (SO) algorithms by taking into consideration the tidal radius around a smaller halo before competitively assigning halo membership to the particles. In this way, the compaso finder allows for more effective deblending of haloes in close proximity as well as the formation of new haloes on the outskirts of larger ones. This halo-finding algorithm is used in the abacussummit suite of N-body simulations, designed to meet the cosmological simulation requirements of the Dark Energy Spectroscopic Instrument (DESI) survey. compaso is developed as a highly efficient on-the-fly group finder, which is crucial for enabling good load-balancing between the GPU and CPU and the creation of high-resolution merger trees. In this paper, we describe the halo-finding procedure and its particular implementation in abacus, accompanying it with a qualitative analysis of the finder. We test the robustness of the compaso catalogues before and after applying the cleaning method described in an accompanying paper and demonstrate its effectiveness by comparing it with other validation techniques. We then visualize the haloes and their density profiles, finding that they are well fit by the NFW formalism. Finally, we compare other properties such as radius–mass relationships and two-point correlation functions with that of another widely used halo finder, rockstar.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Modeling of ATR Fuel in DOE Standard Canisters with Helium Backfill

One pathway for road-ready and final disposition packaging configurations for the aluminum-clad spent nuclear fuel (ANSF) fuel is storage within helium backfilled sealed Department of Energy (DOE) standard canisters. The typical packaging configuration for the 15-foot DOE standard canisters places 10 advance test reactor (ATR) elements a basket, and three baskets are loaded within each DOE canister. During in-reactor operations and cooling pond storage conditions, oxyhydroxide layers form on the surface of the aluminum clad fuel. These layers produce hydrogen gas over time due to the fuel’s radiation field. As part of the packing procedure, the ATR fuel should be dried to remove any residual physio-/chemi- sorbed water bound to the surface. A 50-year CFD model of the DOE canister packaged with fuel was developed to provide a temperature profile for coupled chemical modeling of the conditions within the canister. The results of this modeling include results at fully saturated and fully dried fuel cladding conditions. In the associated experimental work, radiolysis experiments tests were completed in a helium environment, and G-values for the radiolytic production of hydrogen from the oxyhydroxide layers were provided. That reaction was coupled with the thermal profiles and gas-phase reactions to develop a 50-year model of the conditions within a sealed DOE canister with ATR fuel. For a nominal scenario of stored ATR fuel, after 50 years the model results give a 1.36 atm total pressure, 7% mole percent hydrogen, for the upper decay heat, 1.51 atm total pressure, 16% mole percent hydrogen, and for upper decay heat with undried fuel 2.6 atm total pressure, 15% mole percent hydrogen. No case modeled yields significant oxygen, and for the lower decay heat case that is modeled, hydrogen concentrations are under the 4% flammability limit after 50 years of storage. The modeled pressures for all cases modeled are below the pressure limit for the DOE standard sealed canister.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Shock Response of Polymethyl Methacrylate (PMMA) Under Explosive Loading

Abstract The dynamic response of polymethyl methacrylate (PMMA) is well understood for one-dimensional planar impact shocks, but limited research has been performed on the response of PMMA under spherical shock loading. In this work, the shock decay of an explosively-driven shock wave into PMMA was experimentally measured. PMMA cubes of various geometries were explosively loaded with an RP-80 detonator to produce the explosive shock wave. High-speed schlieren imaging was implemented to measure the explosively-driven shock wave velocity throughout the PMMA cubes. Photon Doppler velocimetry (PDV) was used to measure the particle velocity imparted by the shock wave at the surface of the cubes. The material shock response was studied at distances from 21.91 to 133.3 mm from the explosive source. The particle velocity history measured by PDV was compared to the wave profile visualized in the high-speed images. The shock wave pulse amplitude decreased with increased distance from the source. The conducted experiments extend the PMMA shock Hugoniot relating to the lower shock and particle velocity regime.

Torres, S. M. (ORCID:0009000287177206)↗
DOE OSTI2022

Do ambient shear and thermal stratification impact wind turbine tip-vortex breakdown?

Modern wind turbines experience uneven inflow conditions across the rotor, due to the ambient flow’s shear and thermal stratification. Such conditions alter the shape and length of turbine wakes and thus impact the loads and power generation of downstream turbines. To this end, understanding the spatial evolution of the individual wakes under different atmospheric conditions is key to controlling and optimising turbine arrays. With this numerical study we aim to obtain a better understanding of the fundamental physics governing the near-wake dynamics of wind turbines under shear and thermal stability, by examining their tip-vortex breakup mechanisms. Our approach considers scale-resolving simulations of a single turbine wake under a linear shear profile as well as the application of harmonic tip perturbations to trigger flow instabilities. For the subsequent analysis we use the proper orthogonal decomposition (POD) method to extract coherent structures from the flow, and we also calculate mean kinetic energy fluxes to quantify each coherent structure’s contribution to wake recovery. The wake’s helical spiral is found to hinder wake recovery for all studied ambient flow conditions, whereas the mutual inductance instability has positive MKE flux leading to an enhanced wake recovery. Finally, the ambient shear has the largest impact on the local MKE flux with respect to downstream location by changing the shape of the curve and location of extrema, whereas thermal stratification has only a minimal impact on the magnitude of the near-wake local MKE flux distribution.

17 WIND ENERGY↗
DOE OSTI2023

Distributed Data-Driven Optimization for Voltage Regulation in Distribution Systems

Here, this paper proposes a distributed data-driven optimization framework for voltage regulation in distribution systems. The recursive kernel regression and alternating direction method of multipliers (ADMM) are selected to cover the system learning and distributed optimization tasks. The proposed distributed data-driven framework is capable of having a rapid response to system or load changes while considering the operation optimality. Besides, the distributed algorithm parallels the computation tasks and reduces the computational expense of a single agent. To validate the performance of the proposed method, a hypothetical 7-Bus system and the IEEE 123-Bus system are selected to show the effectiveness of the proposed data-driven framework. According to the numerical study results, the proposed method offers great flexibility for selecting customized kernel models for different regions and can effectively improve the system voltage profile in a distributed manner.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Modelling and assessing the near-wake representation and turbulence behaviour of control-oriented wake models

Due to the interaction between the wake of an upstream turbine on a downstream turbine, power losses and increased fatigue loads occur. By yawing the upstream turbine with regard to the wind direction, one can potentially reduce the power losses of the downstream turbine and reduce the fatigue loads. The evolution of the wake depends on the pressure gradient within the near-wake region and the turbulent kinetic energy and must be incorporated in existing wake steering algorithms to provide an accurate estimation of the wake flow. This paper will show a first approach to implement a near-wake model and a turbulence model in the curled wake model within the controls-oriented framework FLORIS. The near-wake model is based on an analytical expression of the velocity profile to model the pressure gradient. Furthermore, two turbulence models are incorporated within the curled wake model based on a Gaussian-distribution and a mixing length formulation. The adapted curled wake model is then assessed with the measurement data acquired in the wind tunnel at ForWind – University of Oldenburg. The evaluation of the models show good agreement for the velocity deficit and representation of the near-wake region. Furthermore, the turbulent kinetic energy behaved as expected in comparison to other work, showing a ring of high turbulent kinetic energy at non-yawed condition which is deflected to a curled shape at large yaw angles with the turbulence model based on a mixing length formulation.

FLORIS↗
DOE OSTI2025

A Dynamic Pricing Method to Manage the Impact of EV Charging on the Grid Using RL

This work addresses the challenge of managing electrical vehicle (EV) charging loads on distribution feeders with the increase in deployment of fast charging stations. To mitigate the adverse impacts on feeder health, a novel dynamic grid-informed pricing approach is proposed. This approach leverages reinforcement learning (RL) to determine hourly charging prices based on real-time grid conditions. A synthetic environment was developed to train the reinforcement learning agent. A model of an IEEE 34-bus distribution feeder with EV charging stations has been developed in OpenDSS utilizing Caldera for realistic EV charging profiles. Test cases demonstrate that the dynamic pricing strategy achieves higher energy delivery to the EV end user at a lower cost compared to constant pricing methods, while lowering voltage deviations and congestion. This approach offers more granular price adjustments, responding dynamically to feeder conditions and potentially improving grid stability and efficiency. The communication architecture to implement this dynamic pricing method is described. This research contributes to the development of smart grid-informed charging solutions that can reduce the cost of charging to the end user while also helping the grid.

EV charging, dynamic pricing, grid-informed chargi↗
DOE OSTI2025

A dynamic pricing method to manage the impact of EV charging on the grid using RL

This work addresses the challenge of managing electrical vehicle (EV) charging loads on distribution feeders with the increase in deployment of fast charging stations. To mitigate the adverse impacts on feeder health, a novel dynamic grid-informed pricing approach is proposed. This approach leverages reinforcement learning (RL) to determine hourly charging prices based on real-time grid conditions. A synthetic environment was developed to train the reinforcement learning agent. A model of an IEEE 34-bus distribution feeder with EV charging stations has been developed in OpenDSS utilizing Caldera for realistic EV charging profiles. Test cases demonstrate that the dynamic pricing strategy achieves higher energy delivery to the EV end user at a lower cost compared to constant pricing methods, while lowering voltage deviations and congestion. This approach offers more granular price adjustments, responding dynamically to feeder conditions and potentially improving grid stability and efficiency. The communication architecture to implement this dynamic pricing method is described. This research contributes to the development of smart grid-informed charging solutions that can reduce the cost of charging to the end user while also helping the grid.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Field Study of 120-volt Heat Pump Water Heaters in the Big Easy

To meet long-term goals for reducing carbon emissions and lessen the impacts of climate change, the U.S. plans to decarbonize its building stock, which includes the replacement of fossil fuel-burning end uses with energy efficient electric alternatives. As part of this strategy, the replacement of fossil fuel water heaters with heat pump water heaters (HPWH) has the potential to avoid substantial carbon emissions. An estimated 10-15 million single-family homes with fossil fuel water heaters do not have the electrical panel capacity to install a 240-volt, 30-amp HPWH. For these homes, a technology recently introduced to the market, the 120-volt plug-in HPWH, can provide energy efficient electrification of water heating without an expensive panel and wiring retrofit. This paper presents the results of an ongoing 120-volt HPWH field study conducted in 17 homes in New Orleans, LA. Key installation scenarios are profiled for retrofitting from gas-fired water heaters to 120-volt HPWHs, accounting for space, air volume, air temperature, condensate drainage, and electrical. Each HPWH had its surrounding air temperature, relative humidity, inlet and outlet water temperature, flow, and energy consumption monitored. Using these data, hot water delivery and energy efficiency performance were analyzed based on home characteristics. Hot water run outs were investigated to understand how hot water usage, inlet water temperature, and surrounding air temperature impact the HPWH’s ability to meet load. From these results, best practices for 120-volt HPWH siting, sizing, and installation were developed. In addition, results are explored further to add insight for electrification policies and product development.

heat pump water heater, 120-volt, Residential buil↗
DOE OSTI2022

Microbial Tracking-2, a metagenomics analysis of bacteria and fungi onboard the International Space Station

The International Space Station (ISS) is a unique and complex built environment with the ISS surface microbiome originating from crew and cargo or from life support recirculation in an almost entirely closed system. The Microbial Tracking 1 (MT-1) project was the first ISS environmental surface study to report on the metagenome profiles without using whole-genome amplification. The study surveyed the microbial communities from eight surfaces over a 14-month period. The Microbial Tracking 2 (MT-2) project aimed to continue the work of MT-1, sampling an additional four flights from the same locations, over another 14 months. Eight surfaces across the ISS were sampled with sterile wipes and processed upon return to Earth. DNA extracted from the processed samples (and controls) were treated with propidium monoazide (PMA) to detect intact/viable cells or left untreated and to detect the total DNA population (free DNA/compromised cells/intact cells/viable cells). DNA extracted from PMA-treated and untreated samples were analyzed using shotgun metagenomics. Samples were cultured for bacteria and fungi to supplement the above results. Staphylococcus sp. and Malassezia sp. were the most represented bacterial and fungal species, respectively, on the ISS. Overall, the ISS surface microbiome was dominated by organisms associated with the human skin. Multi-dimensional scaling and differential abundance analysis showed significant temporal changes in the microbial population but no spatial differences. The ISS antimicrobial resistance gene profiles were however more stable over time, with no differences over the 5-year span of the MT-1 and MT-2 studies. Twenty-nine antimicrobial resistance genes were detected across all samples, with macrolide/lincosamide/streptogramin resistance being the most widespread. Metagenomic assembled genomes were reconstructed from the dataset, resulting in 82 MAGs. Functional assessment of the collective MAGs showed a propensity for amino acid utilization over carbohydrate metabolism. Co-occurrence analyses showed strong associations between bacterial and fungal genera. Culture analysis showed the microbial load to be on average 3.0 × 10 5 cfu/m 2 . Utilizing various metagenomics analyses and culture methods, we provided a comprehensive analysis of the ISS surface microbiome, showing microbial burden, bacterial and fungal species prevalence, changes in the microbiome, and resistome over time and space, as well as the functional capabilities and microbial interactions of this unique built microbiome. Data from this study may help to inform policies for future space missions to ensure an ISS surface microbiome that promotes astronaut health and spacecraft integrity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

THERMAL STRESS ANALYSIS OF A SPENT NUCLEAR FUEL CANISTER

The potential for chloride induced stress corrosion cracking (CISCC) in spent nuclear fuel dry storage canisters is a current topic of research by the US Department of Energy Spent Fuel and Waste Science and Technology program. One of the important prerequisites for CISCC is a tensile residual stress state. This study utilizes computational models to provide an initial analysis of thermal stresses in a generic spent nuclear fuel canister. A STAR-CCM+ thermal fluid model provides a temperature profile analysis of the canister based on four different ambient temperature conditions. An ANSYS APDL finite element model incorporates the temperature profiles to analyze the thermal stresses in the canister. The calculated thermal stress magnitudes are in the range of 10 MPa to 80 MPa, which could be significant for crack propagation through the canister wall. The maximum thermal stress is not a concern for the structural failure of the canister, but it is high enough that it is a significant contribution to the total stress state of the canister wall, which also includes residual stress from fabrication, internal pressure from helium cover gas, and potential transient mechanical loads, such as earthquakes. This paper describes the initial finite element thermal stress analysis that was completed in 2020, reports the initial findings, and identifies areas where model refinement is still needed to complete this work in the future.

Jensen, Ben J.↗
DOE OSTI2021

Flattening the Radial Temperature Profile across the Transformational Challenge Reactor Core

The Transformational Challenge Reactor (TCR) program is demonstrating an agile development approach to advanced nuclear reactor design, which has traditionally utilized a linear design process. In leveraging artificial intelligence, additive manufacturing, advanced materials, and cutting-edge modeling and simulation, the TCR program aims to minimize the high cost and lengthy deployment timelines now standard in the nuclear industry. Within a relatively short period of time, a robust and mature advanced gas-cooled reactor was iteratively designed under the TCR program, using these cutting-edge technologies. The TCR is a 3 MWt gas-cooled microreactor fueled with uranium nitride (UN) tristructural isotropic (TRISO) fuel particles. Though manufactured via traditional means, these UN TRISO particles are loaded into additively manufactured silicon carbide (SiC) cans [4]. Once loaded with TRISO particles, the SiC cans are densified using a chemical vapor infiltration process. The additively manufactured SiC enables significantly more freedom in the design of the fuel form than could ever be achieved using traditionally manufactured SiC. The helium coolant, pressurized to 5 MPa, enters the core at 300°C and nominally exits it at 500°C. Typically, the most thermally limiting components in any reactor design are the fuel assemblies in the core center. To provide a wide thermal margin in these central fuel assemblies, the flow may be biased toward the center of the core to more effectively cool these fuel assemblies with more power deposition and flatten the core’s radial temperature distribution. An analytical fluid model of the TCR core was developed to explore methods for biasing the flow away from the cooler outer fuel assemblies and towards the hotter inner ones. Higher-fidelity models developed in STAR-CCM+ 2020.3.1, a computational fluid dynamics code, were then utilized to verify the analytical model’s findings.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Cyclic bend over sheave performance of a steel tension member belt for wave energy conversion

Rotary power take-offs (PTOs) promise many benefits for commercial wave energy converters (WECs), including longer power strokes/reduced end-stop concerns, deployment advantages, and controllability. However, these winch-type PTOs can impose millions of bending cycles on the belt or rope element every year, which may lead to a failure mode known as Cyclic Bend Over Sheave, or “CBOS”. Synthetic rope, though attractive due to its high strength, low elongation, and durability, suffers from comparatively poor CBOS performance and struggles to reach a million load cycles without failing. Using belts instead of ropes for winching elements should improve CBOS performance due to having a smaller thickness that undergoes bending. Unfortunately, belts have their own host of problems. Maintaining even tension across the width of a belt subject to twist or fleet angles is a challenge, especially for WECs which can move in all 6 degrees of freedom in the ocean. This study sought to quantify the CBOS performance of a belt constructed of steel tension members in a polyurethane matrix containing longitudinal profiles. Belt samples were tested to failure or up to two million bend-unbend cycles and subject to a range of fleet angles. This paper summarizes the utilized testing methodology and the resulting belt performance.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Techno-Economic Analysis and Optimization of a Compressed-Air Energy Storage System Integrated with a Natural Gas Combined-Cycle Plant

To address the rising electricity demand and greenhouse gas concentration in the environment, considerable effort is being carried out across the globe on installing and operating renewable energy sources. However, the renewable energy production is affected by diurnal and seasonal variability. To ensure that the electric grid remains reliable and resilient even for the high penetration of renewables into the grid, various types of energy storage systems are being investigated. In this paper, a compressed-air energy storage (CAES) system integrated with a natural gas combined-cycle (NGCC) power plant is investigated where air is extracted from the gas turbine compressor or injected back into the gas turbine combustor when it is optimal to do so. First-principles dynamic models of the NGCC plant and CAES are developed along with the development of an economic model. The dynamic optimization of the integrated system is undertaken in the Python/Pyomo platform for maximizing the net present value (NPV). NPV optimization is undertaken for 14 regions/cases considering year-long locational marginal price (LMP) data with a 1 h interval. Design variables such as the storage capacity and storage pressure, as well as the operating variables such as the power plant load, air injection rate, and air extraction rate, are optimized. Results show that the integrated CAES system has a higher NPV than the NGCC-only system for all 14 regions, thus indicating the potential deployment of the integrated system under the assumption of the availability of caverns in close proximity to the NGCC plant. The levelized cost of storage is found to be in the range of 136–145 $/MWh. Roundtrip efficiency is found to be between 74.6–82.5%. A sensitivity study with respect to LMP shows that the LMP profile has a significant impact on the extent of air injection/extraction while capital expenditure reduction has a negligible effect.

25 ENERGY STORAGE↗
DOE OSTI2021

Multiobjective Dynamic Optimization for Optimal Load-Following of Natural Gas Combined Cycle Power Plants under Stress Constraints

With the increasing penetration of intermittent renewable energy sources into the electric grid, there is an associated need for conventional thermal power plants that are designed to operate at base-load conditions to cycle their load rapidly and frequently and operate under low-load conditions. Fast load-following operations lead to plant efficiency loss and the reduction of equipment life. A high-fidelity dynamic model of a natural gas combined cycle power plant, with rigorous equipment level submodels, is developed to capture the plant transient performance and the thermomechanical stress evolution at the high-pressure drum is computed to assess the drum life consumption. Stress evolution is modeled for the location of the edge at the drum/downcomer junction that experiences high circumferential stress amplitude during the fast load-following operation, leading to higher fatigue damage than locations considered under existing design standards. A dynamic optimization problem is solved for optimal load-following operation. Depending on the value of the stress constraint and the desired ramp rate, satisfying the desired stress constraint may be infeasible without relaxing the ramp rate. A multiobjective dynamic optimization problem is solved using a lexicographic approach that minimizes ramp rate relaxation and maximizes efficiency while satisfying stress constraints, avoiding spraying to saturation, and maintaining main steam and reheat steam temperature within bounds. It was observed that the optimal ramp rate can be highly nonlinear as opposed to the industry-standard linear ramp rate. Nonlinear optimal ramp-rate profiles not only help to avoid stress constraints that may be unavoidable by using the linear profile but also result in higher efficiency than the linear profile. Furthermore, the study shows that there is a strong tradeoff between the relaxation in the plant ramp rate and the time-average thermal efficiency.

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

EVI-EnSitePy (Electric Vehicle Infrastructure – Energy Estimation and Site Optimization Tool in Python) [EVI-X Modeling Suite] [SWR-25-07]

EVI-EnSitePy is a comprehensive agent-based tool designed for the analysis and design of high-power charging sites, encompassing a wide array of site agents including Electric Vehicles (EVs), chargers, energy storage units (ESS), renewable energy resources (DER), and loads. This versatile tool offers diverse functionalities and a modular modeling approach, allowing detailed configuration of agents based on power ratings, port numbers, energy capacities, demand requirements, charger interfaces, and flexibility to customize the tool for project specific goals. By simulating agent interactions and employing various metrics, EVI-EnSitePy enables the assessment of site performance, exploration of energy management systems (EMS), and implementation of innovative EV charging policies. Utilizing EV charge schedules and arrival states, the tool performs thorough charging site simulations, with outputs consisting of agent and site-level power profiles and statistical metrics. Employing a tree graph structure, EVI-EnSitePy supports nested site structures and power distribution modeling. The tool's ability to generate charging schedules deterministically or via stochastic analysis further enhances its versatility. Through its features and capabilities, EVI-EnSitePy offers a powerful platform for informed decision-making in the realm of high-power charging site design and operation.

Jackson, Derek [National Renewable Energy Laborato↗
DOE OSTI2025

Model Intercomparison of the Impacts of Varying Cloud Droplet–Nucleating Aerosols on the Life Cycle and Microphysics of Isolated Deep Convection

The microphysical impacts of aerosol particles on scattered isolated deep convective cells near Houston, Texas, on 19 June 2013, are examined using multiple cloud-resolving model (CRM) simulations initialized with vertical profiles of low and high concentrations of cloud droplet–nucleating aerosols. These simulations formed part of the Model Intercomparison Project (MIP) conducted by the Deep Convective Working Group of the Aerosol, Cloud, Precipitation and Climate (ACPC) initiative. Each CRM generated a field of convective cells representing those observed during the case study with varying degrees of accuracy. The Tracking and Object-Based Analysis of Clouds (tobac) cell-tracking algorithm was applied to each MIP CRM simulation to track relatively long-lived convective cells (20–60 min). Most of the CRMs produced similar aerosol loading impacts on the warm phase of tracked cell properties with reduced autoconversion and accretion growth of rain, increased cloud water, reduced rainfall, and reduced near-surface evaporation of rain. The sign of aerosol impacts on the warm-phase properties of the convective cells was also quite consistent over cell lifetimes with the greatest magnitude of influence in the first half of the life cycle in most CRMs. In contrast, the ice-phase response to aerosol loading was highly variable among CRMs and included increases or decreases in ice amounts at inconsistent stages of the cell life cycle and midlevel versus upper-level changes in ice. This intermodel variability in ice is indicative both of the complex indirect interactions between aerosols and ice-phase processes in deep convection and their associated parameterizations.

Aerosol-cloud interaction↗
DOE OSTI2021

Fast GPU-Based Generation of Large Graph Networks From Degree Distributions

Synthetically generated, large graph networks serve as useful proxies to real-world networks for many graph-based applications. The ability to generate such networks helps overcome several limitations of real-world networks regarding their number, availability, and access. Here, we present the design, implementation, and performance study of a novel network generator that can produce very large graph networks conforming to any desired degree distribution. The generator is designed and implemented for efficient execution on modern graphics processing units (GPUs). Given an array of desired vertex degrees and number of vertices for each desired degree, our algorithm generates the edges of a random graph that satisfies the input degree distribution. Multiple runtime variants are implemented and tested: 1) a uniform static work assignment using a fixed thread launch scheme, 2) a load-balanced static work assignment also with fixed thread launch but with cost-aware task-to-thread mapping, and 3) a dynamic scheme with multiple GPU kernels asynchronously launched from the CPU. The generation is tested on a range of popular networks such as Twitter and Facebook, representing different scales and skews in degree distributions. Results show that, using our algorithm on a single modern GPU (NVIDIA Volta V100), it is possible to generate large-scale graph networks at rates exceeding 50 billion edges per second for a 69 billion-edge network. GPU profiling confirms high utilization and low branching divergence of our implementation from small to large network sizes. For networks with scattered distributions, we provide a coarsening method that further increases the GPU-based generation speed by up to a factor of 4 on tested input networks with over 45 billion edges.

97 MATHEMATICS AND COMPUTING↗