Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “system partition”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56

Low-frequency variations in the atmospheric branch of the global hydrological cycle

The annual variation of the hydrological cycle is illustrated in terms of hemispheric-mean hydrological variables for the Northern and Southern Hemispheres, while the intraseasonal variations of the global hydrological cycle are illustrated with mean values over two hemispheres that form an east-west partition of the globe. This partition is defined by the 60 deg E-120 deg W great circle and was chosen so that the mean precipitation difference and the divergent water vapor transport between the two hemispheres was maximized. Two years (1979-80) of daily precipitation estimates from the Goddard Laboratory for Atmospheres an 14 years (1979-92) of upper-air data generated by the Global Data Assimilation System at the National Meteorological Center are used in making quantitative estimates of the annual and intraseasonal variations in the global hydrological cycle. The annual variations in hemispheric-mean precipitation (P-circumflex) and water vapor flux divergence (del(vector differential operator) dot Q-circumflex) for the Northern and Southern Hemispheres are comparable with amplitudes of about 0.5 approximately 0.7 mm/day. Both (P-circumflex) and (del(vector differential operator) dot Q-circumflex) vary annually in a coherent way in each hemisphere so that water vapor diverges from the winter hemisphere, where (P-circumflex) reaches its minimum, to the summer hemisphere, where (P-circumflex) attains its maximum. In fact, the hemispheric-mean divergence of water vapor flux changes sign during the annual cycle. Intraseasonal variations of hemispheric-mean precipitation mean P-tilde, evaporation mean E-tilde, and water vapor flux divergence mean del (vector differential operator) dot Q-tilde in the two hemisphres in the east-west direction are comparable with amplitudes of about 0.1 approximately 0.2 mm/day, although amplitudes in some cases exceed 0.3 mm/day. Hemispheric-mean precipitation mean P-tilde varies coherently in opposite phase for the two hemispheres, while mean del (vector differential operator) dot Q-tilde varies so that water vapor diverges from the hemisphere of maximum mean P-tilde to the hemisphere to the hemisphere of minimum mean P-tilde. Intraseasonal variations of mean P-tilde, mean E-tilde, and mean del (vector differential oprator) dot Q-tilde are in accord with the eastward propagation of the intraseasonal global divergent circulation.

Chen, Tsing-Chang↗

Dynamic Forms: Application to Aircraft Guidance - Part 2

The paper describes a method for guiding a dynamic system through a given set of points. The paradigm is a fully automatic aircraft subject to air traffic control (ATC). The ATC provides a sequence of waypoints through which the aircraft trajectory must pass. The waypoints typically specify time, position, and velocity. The guidance problem is to synthesize a system state trajectory that satisfies both the ATC and aircraft constraints. Complications arise because the controlled process is multidimensional, multiaxis, nonlinear, highly coupled, and the state space is not flat. In addition, there is a multitude of operating modes, which may number in the hundreds. Each such mode defines a distinct state space model of the process by specifying the state space coordinatization, the partition of the controls into active controls and configuration controls, and the output map. Furthermore, mode transitions are required to be smooth. The proposed guidance algorithm is based on the inversion of the pure feedback approximation, followed by correction for the effects of zero dynamics. The paper describes the structure and major modules of the algorithm, and the performance is illustrated by several example aircraft maneuvers.

Meyer, George↗

A Catchment-Based Approach to Modeling Land Surface Processes in a GCM: Model Structure - Part 1

A new strategy for modeling the land surface component of the climate system is described. The strategy is motivated by an arguable deficiency in most state-of-the-art land surface models (LSMs), namely the disproportionately higher emphasis given to the formulation of one-dimensional, vertical physics relative to the treatment of horizontal heterogeneity in surface properties -- particularly subgrid soil moisture variability and its effects on runoff generation. The new strategy calls for the partitioning of the continental surface into a mosaic of hydrologic catchments, delineated through analysis of high-resolution surface elevation data. The effective "grid" used for the land surface is therefore not specified by the overlying atmospheric grid. Within each catchment, the variability of soil moisture is related to characteristics of the topography and to three bulk soil moisture variables through a well-established model of catchment processes. This modeled variability allows the partitioning of the catchment into several areas representing distinct hydrological regimes, wherein distinct (regime-specific) evaporation and runoff parameterizations are applied. Care is taken to ensure that the deficiencies of the catchment model in regions of little to moderate topography are minimized.

Koster, Randal D.↗

Modeling the Joint Effects of Vegetation Characteristics and Soil Properties on Ecosystem Dynamics in a Panama Tropical Forest

Abstract In tropical forests, both vegetation characteristics and soil properties are important not only for controlling energy, water, and gas exchanges directly but also determining the competition among species, successional dynamics, forest structure and composition. However, the joint effects of the two factors have received limited attention in Earth system model development. Here we use a vegetation demographic model, the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM), ELM‐FATES, to explore how plant traits and soil properties affect tropical forest growth and composition concurrently. A large ensemble of simulations with perturbed vegetation and soil hydrological parameters is conducted at the Barro Colorado Island, Panama. The simulations are compared against observed carbon, energy, and water fluxes. We find that soil hydrological parameters, particularly the scaling exponent of the soil retention curve ( B sw ), play crucial roles in controlling forest diversity, with higher B sw values (>7) favoring late successional species in competition, and lower B sw values (1 ∼ 7) promoting the coexistence of early and late successional plants. Considering the additional impact of soil properties resolves a systematic bias of FATES in simulating sensible/latent heat partitioning with repercussion on water budget and plant coexistence. A greater fraction of deeper tree roots can help maintain the dry‐season soil moisture and plant gas exchange. As soil properties are as important as vegetation parameters in predicting tropical forest dynamics, more efforts are needed to improve parameterizations of soil functions and belowground processes and their interactions with aboveground vegetation dynamics.

54 ENVIRONMENTAL SCIENCES↗

Acetonitrile Destruction and Fate of Organics in the Reverse Osmosis System at the ETF

The Hanford Site Effluent Treatment Facility (ETF) currently treats aqueous waste streams that include condensates from the 242-A evaporator, leachate from the Environmental Restoration Disposal Facility (ERDF), as well as laboratory wastes and, in the future, will treat liquid effluents from the Hanford Tank Waste Treatment and Immobilization Plant (WTP) and Integrated Disposal Facility (IDF) leachate. Liquid effluents from the WTP will have significant concentrations of acetonitrile. Acetonitrile is formed by reaction of nitrates and sugar in the WTP low activity waste (LAW) melters and is prevalent in the submerged bed scrubber (SBS) and wet electrostatic precipitator (WESP) liquid effluents from WTP off-gas treatment. When these liquids are concentrated in the WTP Effluent Management Facility (EMF) evaporator in the direct feed low activity waste (DFLAW) flow-sheet, testing has shown that the majority of the acetonitrile partitions to the evaporator condensate. Since the evaporator condensate is directed to the ETF, this creates a potential issue with the ETF waste acceptance criteria. Consequently, there is a need to validate flow-sheet assumptions on the fate of acetonitrile and other organics within the ETF. The present plan includes the addition of a steam stripper to the ETF to remove acetonitrile. There is, therefore, also a need to determine a suitable method to destroy acetonitrile in the overhead condensate stream from the new steam stripper. Washington River Protection Solutions, LLC (WRPS) previously contracted with Atkins and the Vitreous State Laboratory (VSL) of The Catholic University of America (CUA) to perform development and testing work to evaluate potential methods for destruction of acetonitrile in WTP secondary liquid effluents. Based on the results of that work, WRPS requested that follow-on testing be conducted to further evaluate acetonitrile destruction in the steam stripper condensate using ultraviolet oxidation (UV/OX) with persulfate. WRPS also requested testing to assess the rejection rate of organics in the reverse osmosis (RO) system installed in the ETF. This report presents the results from testing to address those needs.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Efficient mapping algorithms for scheduling robot inverse dynamics computation on a multiprocessor system

Two efficient mapping algorithms for scheduling the robot inverse dynamics computation consisting of m computational modules with precedence relationship to be executed on a multiprocessor system consisting of p identical homogeneous processors with processor and communication costs to achieve minimum computation time are presented. An objective function is defined in terms of the sum of the processor finishing time and the interprocessor communication time. The minimax optimization is performed on the objective function to obtain the best mapping. This mapping problem can be formulated as a combination of the graph partitioning and the scheduling problems; both have been known to be NP-complete. Thus, to speed up the searching for a solution, two heuristic algorithms were proposed to obtain fast but suboptimal mapping solutions. The first algorithm utilizes the level and the communication intensity of the task modules to construct an ordered priority list of ready modules and the module assignment is performed by a weighted bipartite matching algorithm. For a near-optimal mapping solution, the problem can be solved by the heuristic algorithm with simulated annealing. These proposed optimization algorithms can solve various large-scale problems within a reasonable time. Computer simulations were performed to evaluate and verify the performance and the validity of the proposed mapping algorithms. Finally, experiments for computing the inverse dynamics of a six-jointed PUMA-like manipulator based on the Newton-Euler dynamic equations were implemented on an NCUBE/ten hypercube computer to verify the proposed mapping algorithms. Computer simulation and experimental results are compared and discussed.

Lee, C. S. G.↗

Small vs. Large Convective Cloud Objects from CERES Aqua Observations: Where are the Intraseasonal Variation Signals?

During inactive phases of Madden-Julian oscillation (MJO), there are plenty of deep but small convective systems and far fewer deep and large ones. During active phases of MJO, a manifestation of an increase in the occurrence of large and deep cloud clusters results from an amplification of large-scale motions by stronger convective heating. This study is designed to quantitatively examine the roles of small and large cloud clusters during the MJO life cycle. We analyze the cloud object data from Aqua CERES observations for tropical deep convective (DC) and cirrostratus (CS) cloud object types according to the real-time multivariate MJO index. The cloud object is a contiguous region of the earth with a single dominant cloud-system type. The size distributions, defined as the footprint numbers as a function of cloud object diameters, for particular MJO phases depart greatly from the combined (8-phase) distribution at large cloud-object diameters due to the reduced/increased numbers of cloud objects related to changes in the large-scale environments. The medium diameter corresponding to the combined distribution is determined and used to partition all cloud objects into "small" and "large" groups of a particular phase. The two groups corresponding to the combined distribution have nearly equal numbers of footprints. The medium diameters are 502 km for DC and 310 km for cirrostratus. The range of the variation between two extreme phases (typically, the most active and depressed phases) for the small group is 6-11% in terms of the numbers of cloud objects and the total footprint numbers. The corresponding range for the large group is 19-44%. In terms of the probability density functions of radiative and cloud physical properties, there are virtually no differences between the MJO phases for the small group, but there are significant differences for the large groups for both DC and CS types. These results suggest that the intreseasonal variation signals reside at the large cloud clusters while the small cloud clusters represent the background noises resulting from various types of the tropical waves with different wavenumbers and propagation directions/speeds.

Xu, Kuan-Man↗

A Survey of Multizone Airflow Simulation Software

H-Canyon simulates the building and process ventilation to ensure positive contamination control for any modifications as well as evaluating off-normal conditions. Simulations can also be used to identify air quality of personnel occupied areas. The models can be classified into two groups, single and multi-zone models. As the name states, a single zone model assumes that the building can be described by a single, well-mixed zone. Whereas large buildings have structures that must be described as multizone. Multizone airflow models deal with the difficulty of flows in a building by distinguishing the effects of internal flow restrictions. Several models have been developed that take internal partitions into account. CONTAM (currently used) is a multizone airflow program created by NIST to determine airflows including infiltration, exfiltration, wind pressures acting outside the building, buoyancy effects caused by the indoor and outdoor air temperature difference. COMIS is a multizone air infiltration, ventilation, and contaminant transport model that was developed at Lawrence Berkeley National Laboratory. COMIS deals with applied pressure distribution, airflow components, HVAC systems, schedules, and contaminant transport. EnergyPlus is designed by DOE. It supports flexible component level configuration of HVAC models and simulates sub-hourly time steps to control fast dynamic systems. The program features include providing accurate environmental performance data including energy consumption and carbon emissions, and reporting solar gains on surfaces. Heavent is used to model industrial ventilation systems to control airborne contaminants. It's designed based on the American Conference of Governmental Industrial Hygienists (ACGIH) Ventilation manual. DesignBuilder is the most recognized and advanced user interface to EnergyPlus. It gives access to all of the most frequently required simulation capabilities covering building fabric, thermal mass, glazing, shading renewables, and HVAC analysis. The objective of this study is to identify all relevant software available including commercially licensed, freely licensed, and open source, that can be used to simulate large, complex process facility ventilation systems. All discovered software tools are evaluated for salient features and compared to the software tool currently used to simulate the H-Canyon ventilation system. The design criteria of a ventilation system has two basic requirements: Supplying fresh air for the occupants and changing the air in the room sufficiently; therefore, the smells, fumes, and contaminants will be removed. In this study, first, whether the software simulates airflow (infiltration, natural ventilation, and interzonal) or considers them to be inputs is discussed. Second, if indoor contaminant levels are considered; and finally suggestions on how airflow analysis can be more appropriately accounted for when performing airflow simulation of a building. All of the airflow software reviewed for this study can account for constant infiltration and natural ventilation rates that are not affected by changes in indoor and outdoor conditions. Modelling infiltration in EnergyPlus or DesignBuilder are applicable for small buildings and not appropriate for larger buildings like the H-Canyon. The effect of wind on external pressures, infiltration, and natural and mechanical ventilation, can also be calculated using all of these software. These programs have control options to switch between natural to mechanical ventilation modes or both (hybrid). Interzonal airflow can be applied which means one zone can be specified to receive a fixed amount of air from another zone. The following ventilation modeling software is recommended for further testing: CONTAM, COMIS, Heavent. The next step in evaluating the software will consist of creating test facility simulations using evaluation copies of the software.

42 ENGINEERING↗

IRIS: A Performance-Portable Framework for Cross-Platform Heterogeneous Computing

From edge to exascale, computer architectures are becoming more heterogeneous and complex. The systems typically have fat nodes, with multicore CPUs and multiple hardware accelerators such as GPUs, FPGAs, and DSPs. This complexity is causing a crisis in programming systems and performance portability. Several programming systems are working to address these challenges, but the increasing architectural diversity is forcing software stacks and applications to be specialized for each architecture. As we show, all of these approaches critically depend on their software framework for discovery, execution, scheduling, and data orchestration. To address this challenge, we believe that a more agile and proactive software framework is essential to increase performance portability and improve user productivity. To this end, we have designed and implemented IRIS: a performance-portable framework for cross-platform heterogeneous computing. IRIS can discover available resources, manage multiple diverse programming platforms (e.g., CUDA, Hexagon, HIP, Level Zero, OpenCL, OpenMP) simultaneously in the same execution, respect data dependencies, orchestrate data movement proactively, and provide for user-configurable scheduling. To simplify data movement, IRIS introduces a shared virtual device memory with relaxed consistency among different heterogeneous devices. IRIS also adds an automatic kernel workload partitioning technique using the polyhedral model so that it can resize kernels for a wide range of devices. Our evaluation on three architectures, ranging from Qualcomm Snapdragon to a Summit supercomputer node, shows that IRIS improves portability across a wide range of diverse heterogeneous architectures with negligible overhead.

97 MATHEMATICS AND COMPUTING↗

Recent Advances of PyROS: A Pyomo Solver for Nonconvex Two-Stage Robust Optimization in Process Systems Engineering

The document presents recent algorithmic and implementation advances of the two-stage robust optimization (RO) solver PyROS, and a benchmarking study which demonstrates the utility of PyROS for two-stage RO problems. The advances include extensions of the scope of PyROS to models with uncertain variable bounds, improvements to the initializations of the subproblems used by the underlying cutting set algorithm, and extensions of the uncertainty set interfaces. The benchmarking study is performed on a library of over 8,500 instances, with variations in the nonlinearities, degree-of-freedom partitioning, uncertainty sets, and polynomial decision rule approximations. An amine-based CO2 capture case study is presented to demonstrate the utility of PyROS for large-scale process models. Overall, the results highlight the effectiveness of PyROS for obtaining robust solutions to optimization problems with uncertain equality constraints.

Sherman, Jason↗

Localization of proteins involved in the biogenesis and repair of the photosynthetic apparatus to thylakoid subdomains in Arabidopsis

Abstract Thylakoid membranes in chloroplasts and cyanobacteria harbor the multisubunit protein complexes that catalyze the light reactions of photosynthesis. In plant chloroplasts, the thylakoid membrane system comprises a highly organized network with several subcompartments that differ in composition and morphology: grana stacks, unstacked stromal lamellae, and grana margins at the interface between stacked and unstacked regions. The localization of components of the photosynthetic apparatus among these subcompartments has been well characterized. However, less is known about the localization of proteins involved in the biogenesis and repair of the photosynthetic apparatus, the partitioning of proteins between two recently resolved components of the traditional margin fraction (refined margins and curvature), and the effects of light on these features. In this study, we analyzed the partitioning of numerous thylakoid biogenesis and repair factors among grana, curvature, refined margin, and stromal lamellae fractions of Arabidopsis thylakoid membranes, comparing the results from illuminated and dark‐adapted plants. Several proteins previously shown to localize to a margin fraction partitioned in varying ways among the resolved curvature and refined margin fractions. For example, the ALB3 insertase and FtsH protease involved in photosystem II (PSII) repair were concentrated in the refined margin fraction, whereas TAT translocon subunits and proteins involved in early steps in photosystem assembly were concentrated in the curvature fraction. By contrast, two photosystem assembly factors that facilitate late assembly steps were depleted from the curvature fraction. The enrichment of the PSII subunit OE23/PsbP in the curvature fraction set it apart from other PSII subunits, supporting the previous conjecture that OE23/PsbP assists in PSII biogenesis and/or repair. The PSII assembly factor PAM68 partitioned differently among thylakoid fractions from dark‐adapted plants and illuminated plants and was the only analyzed protein to convincingly do so. These results demonstrate an unanticipated spatial heterogeneity of photosystem biogenesis and repair functions in thylakoid membranes and reveal the curvature fraction to be a focal point of early photosystem biogenesis.

Chotewutmontri, Prakitchai↗

Performance Results on CPU/GPU Exascale Architectures for OMEGA: The Ocean Model for E3SM Global Applications

The US Department of Energy (DOE) conducts climate simulations on some of the world’s largest supercomputers. These exascale machines use heterogeneous architectures with both CPUs and GPUs, and scientific codes must adapt to make full use of this computing power. Los Alamos National Lab is developing Omega: The Ocean Model for E3SM Global Applications, which is specifically designed for modern exascale computers. It uses external libraries that have been optimized for a variety of architectures to run on different supercomputers. Omega is an unstructured-mesh ocean model based on TRiSK numerical methods. It will be the new ocean component of the DOE’s Energy Exascale Earth System Model (E3SM). The algorithms in Omega follow those of the current ocean component, MPAS-Ocean, but it will be written in C++ rather than Fortran to take advantage of the Kokkos performance portability library. Omega spatial operators are written as Kokkos kernels to run efficiently on both CPUs and GPUs. Work on Omega began in 2023 with a new C++ framework for unstructured mesh partitioning, halo exchanges, parallel IO, and Kokkos interfaces. The current version, Omega-0, is being developed to solve the shallow water equations and at present includes all of the tendency terms but not time stepping. Here we share the results of Omega-0 verification and performance testing. Verification includes unit tests implemented with CTest as well as convergence tests in Polaris, an in-house python package with a large suite of test problems. Performance tests compare simulations conducted on CPUs versus GPUs and across different architectures: tests are run on Frontier, which has AMD “Optimized 3rd Gen EPYC” CPUs and AMD MI250X GPUs, as well as Perlmutter, which is composed of AMD EPYC 7763 CPUs and NVIDIA A100 GPUs.

58 GEOSCIENCES↗

Performance modeling and measurement of real-time multiprocessors with time-shared buses

A closed queueing network model is constructed to address workload effects on computer performance for a highly reliable unibus multiprocessor used in real-time control. The queueing model consists of multiserver nodes and a nonpreemptive priority queue. Use of this model requires partitioning the workload into task classes. The time average steady-state solution of the queueing model directly produces useful results that are necessary in performance evaluation. The model is experimentally justified with the Fault-Tolerant Multiprocessor (FTMP) located at the NASA AIRLAB. Extensive experiments are performed on FTMP with a synthetic workload generator (SWG) to directly measure performance parameters, such as processor idle time, system bus contention, and task processing times. These measurements determine values for parameters in the queueing model. Experimental and analytic results are then compared.

Woodbury, Michael H.↗

Interannual Variation in Phytoplankton Class-specific Primary Production at a Global Scale

Phytoplankton is responsible for over half of the net primary production on earth. The knowledge on the contribution of various phytoplankton groups to the total primary production is still poorly understood. Data from satellite observations suggest that for upwelling regions, photosynthetic rates by microplankton is higher than that of nanoplankton but that when the spatial extent is considered, the production by nanoplankton is comparable or even larger than microplankton. Here, we used the NASA Ocean Biogeochemical Model (NOBM) combined with remote sensing data via assimilation to evaluate the contribution of 4 phytoplankton groups to the total primary production. Globally, diatoms were the group that contributed the most to the total phytoplankton production (approx. 50%) followed by coccolithophores and chlorophytes. Primary production by diatoms was highest in high latitude (>45 deg) and in major upwelling systems (Equatorial Pacific and Benguela system). We assessed the effects of climate variability on the class-specific primary production using global (i.e. Multivariate El Nino Index, MEI) and 'regional' climate indices (e.g. Southern Annular Mode (SAM), Pacific Decadal Oscillation (PDO) and North Atlantic Oscillation (NAO)). Most interannual variability occurred in the Equatorial Pacific and was associated with climate variability. These results provide a modeling and data assimilation perspective to phytoplankton partitioning of primary production and contribute to our understanding of the dynamics of the carbon cycle in the oceans at a global scale.

OBM↗

QUEST: systematically approximating Quantum circuits for higher output fidelity

We present QUEST, a procedure to systematically generate approximations for quantum circuits to reduce their CNOT gate count. Our approach employs circuit partitioning for scalability with procedures to 1) reduce circuit length using approximate synthesis, 2) improve fidelity by running circuits that represent key samples in the approximation space, and 3) reason about approximation upper bound. Our evaluation results indicate that our approach of "dissimilar"approximations provides close fidelity to the original circuit. Overall, the results indicate that QUEST can reduce CNOT gate count by 30-80% on ideal systems and decrease the impact of noise on existing and near-future quantum systems.

Patel, Tirthak↗

Flux vector splitting and approximate Newton methods

In the present investigation, the basic approach is employed to view an iterative scheme as Newton's method or as a modified Newton's method. Attention is given to various modified Newton methods which can arise from differencing schemes for the Euler equations. Flux vector splitting is considered as the basic spatial differencing technique. This technique is based on the partition of a flux vector into groups which have certain properties. The Euler equations fluxes can be split into two groups, the first group having a flux Jacobian with all positive eigenvalues, and the second group having a flux Jacobian with all negative eigenvalues. Flux vector splitting based on a velocity-sound speed split is considered along with the use of numerical techniques to analyze nonlinear systems, and the steady Euler equations for quasi-one-dimensional flow in a nozzle. Results are given for steady flows with shocks.

Jespersen, D. C.↗

Water resources and hydrology of Mars

The surface of Mars has been extensively modified by a large variety of water erosional and depositional processes. Although liquid water is presently unstable on the planet's surface, in its cold, hyperarid climate, there is abundant geomorphological evidence of past fluvial valley development multiple episodes of catastrophic flooding, periglacial landforms, ice-related permafrost, lake deposits, eroded impact craters and possible glacial landforms throughout much of Mars' geological history. The amount of water required to form such features is estimated to be equivalent to a planet-wide layer approximately 50 meters deep. Some of this water undoubtedly was removed from the planet by atmospheric escape processes, but much probably remains in the subsurface of Mars. Jakosky summarized the present partitioning of water on Mars, expressed as an average global depth, as follows: in the polar caps, 30 meters; in the megaregolith, 500 to 1000 meters; structurally bound in clays, 10 meters; and in high latitude regolith, a few meters. However, most of this water is probably in the form of ice, except in anomalous areas of possible near surface liquid water, and in regions where hydrothermal systems are still active. The best locations for prospecting are those areas where water or ice is sufficiently concentrated at shallow enough depths to make it feasible to pump out or mine.

Baker, V. R.↗

Privacy-Preserving Control of Partitioned Energy Resources

Distributed energy resources are an increasingly important part of the electric grid. We examine the problem of partitioning a distributed energy resource among many users while providing privacy to them. In this model, clients can send requests to a server, the server can verify that the requests are valid and aggregate them, but it cannot see the actual values in the requests. Without privacy, each user is forced to reveal their daily schedule or energy use. Energy resources add a novel challenge that prior systems do not address: they require verifying limits on private power (a rate over time) and energy (a sum) values. Furthermore, the cryptographic mechanisms must run on embedded energy control systems. We describe Weft, a novel cryptographic system that verifies both power (rate) and energy (integral) constraints on private client values and aggregates them. The key insight behind the approach is to rely on additively homomorphic secret shares, which allows servers to compute sums from rates. We present 3 cryptographic proof systems with different system trade-off for embedded systems: bit-splitting proofs minimize memory use, sorting proofs minimize computation, and commitment proofs minimize network communication. Using bit-splitting proofs, it takes an IoT client using a CortexM microcontroller 4 minutes of compute time to privately control its share of an energy resource for a day at 20s granularity.

Laufer, Evan↗