Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed optimization”

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 271 records · Page 15

A Hierarchical Local Electricity Market for a DER-rich Grid Edge

We report with increasing penetration of distributed energy resources (DER) in the distribution system, it is critical to design market structures that enable smooth integration of DERs. A hierarchical local electricity market (LEM) structure is proposed in this paper with a secondary market (SM) at the lower level representing secondary feeders and a primary market (PM) at the upper level, representing primary feeders, in order to effectively use DERs to increase grid efficiency and resilience. The lower level SM enforces budget, power balance, and flexibility constraints and accounts for costs related to consumers, such as their disutility, flexibility limits, and commitment reliability, while the upper level PM enforces grid physics constraints such as power balance and capacity limits, and also minimizes line losses. The hierarchical LEM is extensively evaluated using a modified IEEE-123 bus with high DER penetration, with each primary feeder consisting of at least three secondary feeders. Data from a GridLAB-D model is used to emulate realistic power injections and load profiles over the course of 24 hours. The performance of the LEM is illustrated by delineating the family of power-injection profiles across the primary and secondary feeders as well as corresponding local electricity tariffs that vary across the distribution grid. Through numerical simulations, the hierarchical LEM is shown to improve the efficiency of the market in terms of lowering overall costs, including both the distribution-level locational marginal prices (d-LMP) as well as retail tariffs paid by customers. Together, it represents an overall framework for a Distribution System Operator (DSO) who can provide the oversight for the entire LEM.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Spatial biases of information influence global estimates of soil respiration: how can we improve global predictions?

Soil respiration (Rs), the efflux of CO2 from soils to the atmosphere, is a major component of the terrestrial carbon cycle, but is poorly constrained from regional to global scales. The global soil respiration database (SRDB) is a compilation of in-situ Rs observations from around the globe that has been consistently updated with new measurements over the past decade. It is unclear whether the addition of data to new versions of the SRDB has produced better-constrained global Rs estimates. We compared two versions of the SRDB (v3.0 n=5173 and v5.0 n=10366) to determine how additional data influenced the global sum and uncertainty of Rs as predicted by a quantile regression forest model. A model parameterized using SRDBv3 yielded a global Rs sum of 88.6 Pg C yr-1, and associated uncertainty of 29.9 (mean absolute error) and 57.9 (standard deviation) Pg C yr-1, whereas parameterization using SRDBv5 yielded 96.5 Pg C yr-1 and associated uncertainty of 30.2 (mean average error) and 73.4 (standard deviation) Pg C yr-1. Empirically estimated global heterotrophic respiration (Rh) from v3 and v5 were 49.9-50.2 (mean 50.1) and 53.3-53.5 (mean 53.4) Pg C yr-1, respectively. Interestingly, SRDBv5’s inclusion of new data from underrepresented regions (e.g., Asia, Africa, South America) thus resulted in overall higher model uncertainty. The SRDBv5 is still biased towards northern latitudes and temperate zones, so we tested that an optimized global distribution of Rs measurements resulted in a global sum of 96.4 ± 21.4 Pg C yr-1 and lower model uncertainty. These results support current global estimates of Rs but highlight spatial inconsistencies that influence model parameterization and interpretation that should be addressed to improve regional- to global-scale predictions.

Stell, Emma↗

Retrospective analysis of biochemical limitations to photosynthesis in 49 species: C 4 crops appear still adapted to pre‐industrial atmospheric [ CO 2 ]

Abstract Leaf CO 2 uptake ( A ) in C 4 photosynthesis is limited by the maximum apparent rate of PEPc carboxylation ( V pmax ) at low intercellular [CO 2 ] ( c i ) with a sharp transition to a c i ‐saturated rate ( V max ) due to co‐limitation by ribulose‐1:5‐bisphosphate carboxylase/oxygenase (Rubisco) and regeneration of PEP. The response of A to c i has been widely used to determine these two parameters. V max and V pmax depend on different enzymes but draw on a shared pool of leaf resources, such that resource distribution is optimized, and A maximized, when V max and V pmax are co‐limiting. We collected published A / c i curves in 49 C 4 species and assessed variation in photosynthetic traits between phylogenetic groups, and as a function of atmospheric [CO 2 ]. The balance of V max ‐ V pmax varied among evolutionary lineages and C 4 subtypes. Operating A was strongly V max ‐limited, such that re‐allocation of resources from V pmax towards V max was predicted to improve A by 12% in C 4 crops. This would not require additional inputs but rather altered partitioning of existing leaf nutrients, resulting in increased water and nutrient‐use efficiency. Optimal partitioning was achieved only in plants grown at pre‐industrial atmospheric [CO 2 ], suggesting C 4 crops have not adjusted to the rapid increase in atmospheric [CO 2 ] of the past few decades.

Pignon, Charles P.↗

Understanding Platinum Ionomer Interface Properties of Polymer Electrolyte Fuel Cells

A well-designed cathode catalyst layer with optimal ionomer distribution is critical to minimizing amount of Platinum (Pt) content in polymer electrolyte fuel cells (PEFCs). The impact of Pt loading, ionomer content and carbon support types on the catalyst/ionomer interface were also investigated at dry and wet conditions. Higher Pt loadings resulted in higher double layer capacity (C dl ) and similar electrochemical surface area (ECSA) due to well dispersed ionic phase material. Higher ionomer content resulted in higher ionic conductivity but also showed similar SO 3 − group coverage. High surface area (HSA) carbon support had larger ECSA and C dl at both dry and wet conditions, as less agglomerated Pt was well dispersed in the meso-pores of the support. Lower SO 3 − group coverages were observed for HSA carbon than for Vulcan carbon due to Pt particles being buried within the porous HSA carbon support. The effect of cell conditioning and voltage recovery on the PEFC cathode catalyst layer was shown to have minimal impact on SO 3 − group coverage despite a decrease in C dl and ECSA due to the size increase of Pt particles. At dry condition, a significant increases in SO 3 − group coverage were observed for all MEAs due to higher adsorptivity of ionomer in dry conditions.

25 ENERGY STORAGE↗

Modeling Synergistic Fuel Cell Membrane Degradation with Mitigating Effects of Cerium

During operation, polymer-electrolyte-membrane (PEM) fuel cells undergo mechanical and chemical degradation mechanisms, which behave synergistically and lead to accelerated membrane degradation over time. This study builds upon previous modeling work on mechanical degradation as described by a pinhole in the membrane and the effects of cerium on chemical degradation.1 By combining these two models, analysis can be carried out on the coupled degradation methods and how the mitigation effects of cerium disrupt the degradation cycle. The mechanical model represents a pinhole in the membrane using an effective void fraction, which allows for increased gas crossover through the membrane and impacts the membrane transport and mechanical properties. A microkinetic model for the chemical degradation is included in the model, including attack of the membrane polymer by hydroxyl radicals as well as the mitigation reaction for quenching of hydroxyl radicals with cerium. A concentrated solution theory approach is used to model the transport of cerium ions throughout the cell.2 The model results show how the location of cerium in the cell can be used to prevent chemical degradation. In addition, the model shows how cerium slows down the rate of pinhole growth by reducing the gas crossover and membrane thinning rate. Finally, the model can be used to optimize the distribution of cerium in the membrane and catalyst layers by balancing trade-offs between lowering degradation rates and decreasing fuel cell performance.

25 ENERGY STORAGE↗

Knowledge Oriented Graph Unified Transformer (KOGUT) v0.1

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

Joachimiak, Marcin [Lawrence Berkeley National Lab↗

OASIS: Optimal Analysis-Specific Importance Sampling for event generation

We propose a technique called Optimal Analysis-Specific Importance Sampling (OASIS) to reduce the number of simulated events required for a high-energy experimental analysis to reach a target sensitivity. We provide recipes to obtain the optimal sampling distributions which preferentially focus the event generation on the regions of phase space with high utility to the experimental analyses. OASIS leads to a conservation of resources at all stages of the Monte Carlo pipeline, including full-detector simulation, and is complementary to approaches which seek to speed-up the simulation pipeline.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An MILP-Based Distributed Energy Management for Coordination of Networked Microgrids

An MILP-based distributed energy management for the coordination of networked microgrids is proposed in this paper. Multiple microgrids and the utility grid are coordinated through iteratively adjusted price signals. Based on the price signals received, the microgrid controllers (MCs) and distribution management system (DMS) update their schedules separately. Then, the price signals are updated according to the generation–load mismatch and distributed to MCs and DMS for the next iteration. The iteration continues until the generation–load mismatch is small enough, i.e., the generation and load are balanced under agreed price signals. Through the proposed distributed energy management, various microgrids and the utility grid with different economic, resilient, emission and socio-economic objectives are coordinated with generation–load balance guaranteed and the microgrid customers’ privacy preserved. In particular, a piecewise linearization technique is employed to approximate the augmented Lagrange term in the alternating direction method of multipliers (ADMM) algorithm. Thus, the subproblems are transformed into mixed integer linear programming (MILP) problems and efficiently solved by open-source MILP solvers, which would accelerate the adoption and deployment of microgrids and promote clean energy. The proposed MILP-based distributed energy management is demonstrated through various case studies on a networked microgrids test system with three microgrids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Carbon Dioxide Pressure and Catalyst Quantity Dependencies in Artificial Photosynthesis of Hydrocarbon Chains on Nanostructured Co/CoO Surfaces

In this study, we investigated the influence of pressure and the quantity of Co/CoO catalyst on an artificial photosynthesis process that converts CO 2 and H 2 O into hydrocarbons (C n H 2n+2 , where n ≤ 18). The adsorption of CO 2 and H 2 O on Co/CoO surfaces proved to be pivotal in this photo-catalytic reaction. Photoexcited carbon dioxide and water molecules ((CO 2 )* and (H 2 O)*) generated by illuminating the catalyst surface led to the formation of alkene hydrocarbon molecules with carbon numbers following an approximate Poisson distribution. The optimal pressure was found to be 0.40 MPa. Pressure less than 0.40 MPa resulted in low CO 2 adsorption, impeding excitation for photosynthesis. At greater pressure, oil/wax accumulation on Co/CoO surfaces hindered CO 2 adsorption, limiting further photosynthesis reactions. The average number of carbon atoms in the hydrocarbons and hydrocarbon yield were correlated. The amount of Co/CoO was also found to affect the hydrocarbon yield. Our study contributes to the understanding of Co/CoO-catalyzed photosynthesis and suggests that an open-flow system could potentially enhance the productivity of long-chain hydrocarbons.

59 BASIC BIOLOGICAL SCIENCES↗

Validation of the Scientific Program for the Dark Energy Spectroscopic Instrument

The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^{2}$ over 5 yr to constrain the cosmic expansion history through precise measurements of baryon acoustic oscillations (BAO). The scientific program for DESI was evaluated during a 5 month survey validation (SV) campaign before beginning full operations. This program produced deep spectra of tens of thousands of objects from each of the stellar Milky Way Survey (MWS), Bright Galaxy Survey (BGS), luminous red galaxy (LRG), emission line galaxy (ELG), and quasar target classes. These SV spectra were used to optimize redshift distributions, characterize exposure times, determine calibration procedures, and assess observational overheads for the 5 yr program. In this paper, we present the final target selection algorithms, redshift distributions, and projected cosmology constraints resulting from those studies. We also present a One-Percent Survey conducted at the conclusion of SV covering 140 deg$^{2}$ using the final target selection algorithms with exposures of a depth typical of the main survey. The SV indicates that DESI will be able to complete the full 14,000 deg$^{2}$ program with spectroscopically confirmed targets from the MWS, BGS, LRG, ELG, and quasar programs with total sample sizes of 7.2, 13.8, 7.46, 15.7, and 2.87 million, respectively. These samples will allow exploration of the Milky Way halo, clustering on all scales, and BAO measurements with a statistical precision of 0.28% over the redshift interval z < 1.1, 0.39% over the redshift interval 1.1 < z < 1.9, and 0.46% over the redshift interval 1.9 < z < 3.5.

79 ASTRONOMY AND ASTROPHYSICS↗

Emergency Generation Asset Positioning for Resilient Transmission Grid Operation: Preprint

In this paper we study a problem of hurricane emergency preparedness via placement of emergency generation assets prior to its strike. We present a two-stage stochastic model for choosing locations and quantities of emergency storage to help support the power grid through a hurricane event. The expectation of losses in our two-stage model is estimated using the sample average approximation. We construct damage scenarios for sample average approximation using WIND Toolkit meteorological data and fragility curves of various electric grid components. We demonstrate the efficacy of our two-stage planning model by simulating operations during Hurricane Dolly on the 2000-bus transmission test system. Our model, coupled with our scenario selection strategy, is effective at mitigating loss of load when compared to a model without emergency generation assets placed prior to an extreme event.

emergency preparedness↗

Summary of a Workshop on Integrating Energy Equity into Power System Planning and Operation

On November 29 and 30, 2023, the Energy Systems Integration Group convened a one-and-a-half-day workshop in Arlington, Virginia, to improve knowledge and practices for integrating energy equity considerations into power system planning and operation. Workshop participants were invited to discuss approaches for integrating energy equity aspects into models and tools and to explore possible pathways for enhanced technical practices. Participants included utility planning and modeling practitioners, social science and energy justice academics, representatives from energy and environmental justice nongovernmental organizations, non-profits, consumer advocates, and representatives from the U.S. Department of Energy, national labs, and other government agencies. This report summarizes the discussion points and outcomes of the workshop, which prioritized the technical aspects of incorporating the challenges and solutions identified by electricity-futures stakeholders and communities. These technical activities include the fundamental grid structure and hardware, detailed modeling approaches that benchmark energy system performance and justify future investments, and engineering designs and limitations for potential solutions to meet energy equity goals in communities.

energy equity↗

Sub-microsecond Transformers for Jet Tagging on FPGAs

We present the first sub-microsecond transformer implementation on an FPGA achieving competitive performance for state-of-the-art high-energy physics benchmarks. Transformers have shown exceptional performance on multiple tasks in modern machine learning applications, including jet tagging at the CERN Large Hadron Collider (LHC). However, their computational complexity prohibits use in real-time applications, such as the hardware trigger system of the collider experiments up until now. In this work, we demonstrate the first application of transformers for jet tagging on FPGAs, achieving $\mathcal{O}(100)$ nanosecond latency with superior performance compared to alternative baseline models. We leverage high-granularity quantization and distributed arithmetic optimization to fit the entire transformer model on a single FPGA, achieving the required throughput and latency. Furthermore, we add multi-head attention and linear attention support to hls4ml, making our work accessible to the broader fast machine learning community. This work advances the next-generation trigger systems for the High Luminosity LHC, enabling the use of transformers for real-time applications in high-energy physics and beyond.

Laatu, Lauri [Imperial Coll., London]↗

Distribution Market-Clearing and Pricing Considering Coordination of DSOs and ISO: An EPEC Approach

Distribution-level electricity market provides a platform for trading energy and grid services from large-amount of small-scale distributed energy resources (DERs) located on distribution grids. The behavior of the DERs in the distribution electricity market may ultimately impact the market-clearing and locational marginal prices (LMPs) in the wholesale market. This paper proposes a bi-level optimization model for distribution market clearing and distribution locational marginal pricing (DLMP) considering the interactions between distribution and transmission wholesale markets. In the proposed model, the upper-level model represents the distribution system operator (DSO) market-clearing and the lower-level model represents the wholesale market-clearing by the independent system operator (ISO). The LMP at the substation will impact the DER dispatch and power demands of the DSO as well as the DLMP. In turn, the power demands of the DSO will further impact the ISO market and its LMPs. The equilibrium problem with equilibrium constraints (EPEC) approach is applied to find the equilibria of multiple DSOs and the ISO. The EPEC problem is transformed into a single-level mixed-integer convex problem in order to allow for efficient solving. Finally, the effectiveness of the proposed model and solution method are demonstrated through case studies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Distribution Feeder-Scale Fast Frequency Response via Optimal Coordination of Net-load Resources Part II: Large-Scale Demonstration

This work is the second of a two-part series in which we develop and experimentally demonstrate a hierarchical control solution for optimally coordinating thousands of deferrable loads and distributed energy resources (DERs) to provide fast frequency response (FFR) from an entire distribution feeder. In Part I, we developed and proved practical algorithms for fast, cost-based optimal dispatch and for determining the optimal amount of headroom to operate solar inverters with to support FFR dispatch while minimizing opportunity cost. Simulation results in Part I demonstrated the advantages of the hierarchical dispatch approach in being able to maintain fast solution times needed for FFR even when the problem size increases. In Part II, we implement the algorithms developed in Part I in a novel, large-scale power hardware-in-the-loop experiment including embedded controllers and more than 100 powered appliance loads and DER connected to a simulated real-world distribution system with more than 10,000 controlled devices. Experimental results from multiple scenarios confirm that the optimal FFR dispatch approach scales well and can optimally coordinate more than 10,000 net-load resources across a distribution network while achieving hardware response times within 500 ms, which is not possible using state-of-the-art optimal coordination approaches.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The Challenges of Modeling Distributed Energy Resources (DERs) as Blackstart Resources and for Volt-VAR Optimality

Modeling, testing and instituting Distributed Energy Resources (DERs) as Blackstart Resoucrces presents several challenges due to the fundamental differences between traditional black start resources (e.g., large generators) versus DERs like solar PV systems with battery storage. This paper addresses some of these key challenges. These include intermittency of coordinating DERs to provide continuous power during a black start event, especially during extended periods of cloud cover or when battery energy storage is depleted. This paper additionally addresses the collapsing voltage and stability control challenges specific to maintaining bulk power system stability during black start synchronization These physical and engineering limitations require careful engineering design, modeling and engineering to ensure that DERs can support critical loads and substations during black start events. A variety of additional challenges also exist. Additionally, there are challenges associated with feeder location and low voltage secondary system impacts on DER functions and settings. We compare typical functions and settings for DERs for power factor control and correction. We also demonstrate how voltag control via Volta-VAR power factor correction can be done

Mukherjee, Srijib↗

Real-Time Optimization and Control of Next-Generation Distribution Infrastructure

This presentation highlights the key developments of the ARAP-e NODES project, including the innovative real-time optimal power flow (RT-OPF) algorithm for distributed energy resource (DER) management, trip planning for T+D coordination, and extensive validation of the RT-OPF algorithm in NREL ESIF laboratory with 100+ physical hardware devices, controller-hardware-in-the-loop test at Southern California Edison (SCE), and two field demonstrations.

hardware-in-the-loop↗