Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “design constraints”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Technical and Economic Evaluation of the First Ever Polymer Flood Field Pilot to Enhance the Recovery of Heavy Oils on Alaska's North Slope via Machine Assisted History Matching

Polymer flooding has become globally established as a potential enhanced oil recovery method for heavy oils. To determine whether this technology may be useful in developing the substantial heavy oil resources on the Alaska North Slope, a polymer flood field pilot commenced at the Milne Point Unit in August 2018. This study seeks to evaluate the results of the field pilot on a technical and economic basis. A reservoir simulation model is constructed and calibrated to predict the oil recovery performance of the pilot through machine-assisted reservoir simulation techniques. To replicate the early water breakthrough observed during waterflooding, transmissibility contrasts are introduced into the simulation model, forcing viscous fingering effects. In the ensuing polymer flood, these transmissibility contrasts are reduced to replicate the restoration of injection conformance during polymer flooding. Transmissibility contrasts are later reinstated to replicate fracture overextension interpreted in one of the producing wells. The calibrated simulation models produced at each stage of the history matching process are used to forecast oil recovery. These forecasts are used as input for economic analysis, incremental to waterflooding expectations. The simulation forecasts indicate that polymer flooding significantly increases the heavy oil production for this field pilot compared to waterflooding alone, yielding attractive project economics. However, meaningful variations between simulation scenarios demonstrate that a simulation model is only valid for prediction if flow behavior in the reservoir remains consistent with that observed during the history matched period. Critically, this means that a simulation model calibrated for waterflooding may not fully capture the technical and economic benefits of an enhanced oil recovery process such as polymer flooding. Subsequently, the simulation model and economic model are used in conjunction to conduct a sensitivity analysis for polymer flood design parameters, from which recommendations are provided for both the continued operation of the current field pilot and future polymer flood designs. The results demonstrate that a higher polymer concentration can be injected due to the development of fractures in the reservoir. The throughput rate should remain high without exceeding operating constraints. A calculated point-forward polymer utilization parameter demonstrates the decreasing efficiency of the polymer flood at later times in the pattern life. Future projects will benefit from starting polymer injection earlier in the pattern life. A pattern with tighter horizontal well spacing will observe a greater incremental benefit from polymer flooding.

Keith, Cody↗
DOE OSTI2025

Nondestructive Evaluation of Concrete: Elastic Property Imaging Through Full-Waveform Inversion

Concrete is a vital material in construction—especially in the nuclear industry, where it is used in critical structures such as containment vessels. Over time, concrete can degrade due to harsh operational and environmental conditions, necessitating that its elastic properties be accurately evaluated to ensure structural integrity and safety. Traditional nondestructive evaluation methods such as ultrasound-based techniques often rely on simplifying assumptions that may not hold true for concrete. This paper presents an advanced ultrasound-based method that uses elastic full-waveform inversion (EFWI) to create detailed images of concrete’s mechanical properties. By accurately modeling wave behaviors such as scattering and reflection, we aim to overcome the limitations of conventional ultrasonic-based methods. In this work, the imaging problem involved reconstructing the various elastic properties of a heterogenous concrete block with three steel rebars embedded in it. The ultrasonic measurements were synthetically generated from multiple sources and receivers, and the reconstruction process was performed using a gradient-based optimization algorithm. Our approach leveraged EFWI to reconstruct high-resolution images of the pressure wave speed, shear wave speed, and density. Multiple misfit functions—including L2-norm, cross-correlation (CC), and L1-norm—combined with total variation (TV) regularization and parameter constraints using a Sigmoid function—were explored for the reconstruction. The results demonstrated that using the L1-norm misfit function in conjunction with TV regularization and Sigmoid constraints significantly improved the reconstruction quality in comparison to traditional methods. This approach provided clearer images with fewer artifacts and better captured background heterogeneity. Our findings highlight that, when properly designed, EFWI carries great potential for providing comprehensive, more accurate, and more reliable assessments of concrete conditions, as is crucial for the maintenance and safety of nuclear power plant structures.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Optimizing Information Automation Using a New Method Based on System-Theoretic Process Analysis: Tool Development and Method Evaluation

This report is an update to a prior report that describes progress and findings for a program of research supporting the design and optimization of information automation systems for nuclear power plants. Much of the domestic nuclear fleet is currently focused on modernizing technologies and processes, including transitioning toward digitalization in the control room and throughout the plant, along with a greater interest in the use of automation, artificial intelligence, robotics, and other emerging technologies. While there are significant opportunities to apply these technologies toward greater plant safety, efficiency, and overall cost-effectiveness, optimizing their design and avoiding potential safety and performance risks depends on ensuring that human performance-related organizational and technical design issues are identified and addressed early in the design process. This report describes modeling tools and techniques, based on sociotechnical systems theory, to support these design goals and their application in the current research effort. The report is primarily intended for senior nuclear energy stakeholders, including regulators, corporate management, and senior plant management. We have developed and employed a method to design an optimized information automation ecosystem (IAE) based on the systems-theoretic constructs underlying sociotechnical systems theory in general and the Systems-Theoretic Accident Modeling and Processes (STAMP) approach in particular. We argue that an IAE can be modeled as an interactive information control system whose behavior can be understood in terms of dynamic control, feedback, and communication relationships amongst the system’s technical and organizational components. We have employed two STAMP-based tools in this effort. The first is Causal Analysis based on STAMP (CAST), an accident and incident analysis technique that was used to examine a performance- and safety-related incident at an industry partner’s plant involving the unintentional activation of an emergency diesel generator. This analysis provided insight into the behavior of the plant’s current information control structure within the context of a specific, significant event. The second tool is Systems Theoretic Process Analysis (STPA) which is a proactive risk analysis tool used to examine existing and potential, planned sociotechnical systems. STPA was used to identify risk factors in the current design of a generic nuclear power plant (NPP) preventive maintenance system. Our analyses focused on identifying near-term system improvements and longer-term design requirements for an optimized IAE system. CAST analysis findings indicate an important underlying contributor to the incident under investigation, and a significant risk to information automation system performance, was perceived time and schedule pressure, which exposed weaknesses in interdepartmental coordination between and within responsible plant organizations and challenged the resilience of established plant processes, until a human caused the eventual event. These findings are discussed in terms of their risk to overall system performance and their implications for information automation system resilience and brittleness. STPA findings exposed several areas of concern in the design of current preventive maintenance systems. We also present two preliminary information automation models. The proactive issue resolution (PIR) model is a test case of an information automation concept with significant near-term potential for application and subsequent reduction in significant plant events. The IAE model is a more general representation of a broader, plantwide information automation system and represents an end-state vision for our work. From our results, we have generated an initial set of preliminary system-level requirements and safety constraints for these models. We have also focused on early development of easy to learn, easy to use “transportable” tools for sociotechnical systems analysis. We intend these to be used by NPP personnel as a means of gaining reliable and relatively quick insight into (1) sociotechnical systems factors impacting incidents and accidents, (2) potential sociotechnical risk factors in existing or planned system designs, and (3) potential weaknesses in a system’s safety and/or information control structure. We conclude the report with a set of summary recommendations, a discussion of planned and potential follow-on research and development, and a draft list of system-level requirements and safety constraints for optimized information automation systems.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

Achieving high ethylene yield in non-oxidative ethane dehydrogenation

Steam cracking of ethane, a non-catalytic thermochemical process, remains the dominant means of ethylene production. The severe reaction conditions and energy expenditure involved in this process incentivize the search for alternative reaction pathways and reactor designs which maximize ethylene yield while minimizing cost and energy input. Herein, we report a comparison of catalytic and non-catalytic non-oxidative dehydrogenation of ethane. We achieve ethylene yields as high as 67 % with an open tube quartz reactor without the use of a catalyst at residence times ~4 s. The open tube reactor design promotes simplicity, low cost, and negligible coke formation. Pristine quartz tubes were most effective, since coke formation was detected when defects were introduced by scratching the surface of the quartz. Surprisingly, the addition of solids to the quartz tube, such as quartz sand, alumina powder, or even Pt-based intermetallic catalysts, led to lower ethylene yield. Pt alloy catalysts are effective at lower temperatures, such as at 575 °C, but conversion is limited due to thermodynamic constraints. When operated at industrially relevant temperatures, such as 700 °C and above, these catalysts were not stable in our tests, causing ethylene yield to drop below that of the open tube. Furthermore, these results suggest that future research on non-oxidative dehydrogenation should be directed at optimizing reactor designs to improve the conversion of ethane to ethylene, since this approach shows promise for decentralized production of ethylene from natural gas deposits.

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

A strategy for automated core design to increase economic viability and minimize fuel fragmentation, relocation, and dispersal susceptibility in high-burnup cores

The nuclear industry aims to increase the cycle length of pressurized water reactors from 18 to 24 months to increase power plant capacity factors and economic viability. These cycle length extensions will inherently require fuel rods to exceed the current peak rod average burnup limit of 62 GWd/MTU. A chief concern of operating beyond the current burnup limit is the fuel fragmentation, relocation, and dispersal (FFRD) phenomenon in which pulverized fuel fragments can axially relocate and escape through a burst in the cladding formed during a loss-of-coolant accident. In this work, we demonstrate an approach for automating core design employing an optimization tool based on a penalty-free, parallel simulated annealing algorithm to produce pressurized water reactor core designs with two different optimization objectives. The two objectives were to produce core designs with (1) mitigated FFRD susceptibility while achieving 24-month cycle lengths (2) maximum cycle length with no regard for the likelihood of FFRD. Batch size was considered in tandem with both cases to maximize economic viability. The PARCS nodal model was the primary reactor physics tool used in the optimizations and used nuclear cross sections calculated with 2D Polaris lattice physics models. Reactor performance and safety characteristics of the optimized cores were verified using high-fidelity Virtual Environment for Reactor Applications models. The core designs produced by the optimization tool are compared with each other and to a high-burnup core design produced and analyzed in previous works to highlight the fuel management strategies that may enhance high-burnup reactor safety and economic viability. The optimized cores satisfied their respective objective functions, producing a maximum cycle length of 720 effective full-power days in one core design and one that may reduce FFRD susceptibility by up to 50% based on the first-order approximation to FFRD risk formulated in this work. The optimized cores met most constraints but exceeded the hot channel factor limit, especially in FFRD cases where fresh fuel carried more power. Furthermore, this highlights the need for future lattice-level optimizations and broader assembly options.

Cycle length↗
DOE OSTI2025

Distribution System Blackstart and Restoration Using DERs and Dynamically Formed Microgrids

Extreme weather events have led to long-duration outages in the distribution system (DS), necessitating novel approaches to blackstart and restore the system. Existing blackstart solutions utilize blackstart units to establish multiple microgrids (MGs), sequentially energize non-blackstart units, and restore loads. However, these approaches often result in isolated MGs. In DERs-aided blackstart, the continuous operation of these MGs is limited by the finite energy capacity of commonly used blackstart units like battery energy storage (BES)-based gridforming inverters (GFMIs). To address this issue, this article proposes a holistic blackstart and restoration framework that incorporates synchronization between dynamic MGs and the entire DS with the transmission grid (TG). To support synchronization, we leveraged virtual synchronous generator-based control for GFMIs to estimate their frequency response to load pick-up events using only initial/final quasi-steady-state points. Subsequently, a synchronization switching condition is developed to model synchronizing switches, aligning them seamlessly with a linearized branch flow problem. Finally, we designed a bottomup blackstart and restoration framework that considers the switching structure of the DS, energizing/synchronizing switches, DERs with grid-following inverters, and BES-based GFMIs with frequency security constraints. In conclusion, the proposed framework is validated in IEEE-123-bus system, considering cases with two and four GFMIs under various TG recovery instants.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

A Non-Cooperative Game-Based Distributed Beam Scheduling Framework for 5G Millimeter-Wave Cellular Networks

Here, this paper studies the problem of distributed beam scheduling for 5G millimeter-Wave (mm-Wave) cellular networks where base stations (BSs) belonging to different operators share the same spectrum without centralized coordination among them. Our goal is to design efficient distributed scheduling algorithms to maximize the network utility, which is a function of the achieved throughput by the user equipment (UEs), subject to the average and instantaneous power consumption constraints of the BSs. We propose a Media Access Control (MAC) and a power allocation/adaptation mechanism utilizing the Lyapunov stochastic optimization framework and non-cooperative games. In particular, we first decompose the original utility maximization problem into two sub-optimization problems for each time frame, which are a convex optimization problem and a non-convex optimization problem, respectively. By formulating the distributed scheduling problem as a non-cooperative game where each BS is a player attempting to optimize its own utility, we provide a distributed solution to the non-convex sub-optimization problem via finding the Nash Equilibrium (NE) of the game whose weights are determined optimally by the Lyapunov optimization framework. Finally, we conduct simulation under various network settings to show the effectiveness of the proposed game-based beam scheduling algorithm in comparison to that of several reference schemes.

42 ENGINEERING↗
DOE OSTI2022

Hydrogen Energy Storage Integrated with a Combined Cycle Plant

A project is being developed that will build upon the existing infrastructure and resources at the Intermountain Power Project (IPP) site to provide reliable, dispatchable energy and to support the transmission of renewable energy resources while transitioning to an economical green energy future. The concept study depicted in this report outlines a techno-economic optimization to fulfill the demand for 30% vol hydrogen co-firing in the IPP 840 MW advanced class combined cycle power plant. In an initial step, a site assessment concluded the site has sufficient land available to co-locate a hydrogen production and storage facility. The team evaluated and defined a scalable concept that considered technology characteristics, including input and output models to be used for optimization purposes. The concept for the hydrogen production and storage system integrates multiple technologies, to determine system size and scalable approach, for each of the technologies evaluated, the team defined component and subcomponent sizes, minimum and maximum capacity, modularity, component utility consumption (electric, water), component flexibility and servicing, layout, and technology status, as well as technology alternatives. For hydrogen generation, the project considers Siemens Energy’s Silyzer-300 (S300) technology, a 17.5 MW modular Proton Exchange Membrane (PEM) electrolyzer. For the S300 configuration, the team determined that three S300 arrays, or approximately 1,000 kg/hr, per block would yield a compact block design. This configuration results in a fairly wide and flexible arrangement that fits well into the spaces available at the site. Therefore, the overall design approach is based on multiple identical blocks of 3 arrays to minimize engineering cost and optimize constructability. In parallel, a transmission screening study was conducted to determine any potential transmission constraints from the energy sources that could feed the hydrogen production equipment. The study results show that minimum transmission constraints would be encountered to deliver 400 MW renewable generation from southern California, or south-central Wyoming. At last, the techno-economic analysis concluded that a scenario that uses solar and wind power yields the lowest levelized cost of hydrogen (LCOH 2 ) production and the lowest cost per tonne of CO 2 reduced. In this optimized scenario, the hydrogen production plant was determined as 6,201 kg/hr and the hydrogen storage (underground cavern) was determined as 4,600 tonnes. The resulting capacity factor for the hydrogen production plant was 66.33% with 8,745 operating hours in one year. This techno-economic analysis provided various options for integrating hydrogen storage at the Intermountain Power Plant site to co-fire the CCPP units. The results provide insightful data about the magnitude of capacity needed and the economics of producing hydrogen and reducing CO 2 emissions.

08 HYDROGEN↗
DOE OSTI2026

Deployable UHV Pump

This project presents an external motorized actuation system for a deployable ultra-high vacuum (UHV) pump, enabling internal motion without compromising vacuum conditions. All active components remain external due to environmental and space constraints. A motor-driven mechanical feedthrough transfers motion into the pump, supported by a modular, adjustable mounting system that maintains alignment and integrates with existing hardware. CAD modeling and iterative design were used to refine geometry and ensure proper fit. The final design reliably transfers motion while maintaining alignment and structural integrity. Its adjustability improves installation and maintenance, demonstrating a practical solution for actuation in UHV systems.

Remington, Austin [Northern Illinois U.]↗
DOE OSTI2020

VALVE TIMING AND AREA PROFILE SELECTION FOR HYDRAULIC PUMPS AND MOTORS

A significant amount of research has been conducted to select valve timing and area profiles that create efficient and quiet hydraulic pumps and motors. Numerous active valve architectures have been modeled and optimized, but the rationale for the final solution is often unclear. The solution is usually highly dependent on the modeled valve geometry constraints and the duty cycle of the pump or motor for which the valve was optimized. This paper first presents an ideal valve area profile and what compromises are made to design for a working pressure and speed range. The ideal valve area profile is resolved using a piston-cylinder model with a pressure dependent bulk modulus and variable valve openings. The piston-cylinder is first modeled as a pump with passive (check) valves and a specified pressure drop across the valves to determine the ideal area profile. A generalization of the ideal timing is demonstrated for positive pressure differential and positive rotation direction, also known as the first quadrant. The effects on valve timing when varying pressure are shown for the first quadrant, but the trends can be extrapolated for all quadrants of operation. The ideal valve profile is implemented as fixed valve timing, meaning the valve area is only a function of the timing angle, in a pump-motor. The pump-motor is simulated in one rotation direction through a pressure range. Performance is high in pumping operation, but when the pressure differential is reversed, cylinder pressure spikes ensue. Two strategies to modify an ideal valve area profile are presented: timing grooves and a pressure shifted valve timing. Timing grooves reduce pressure spikes and cavitation in the cylinder but generally increase throttling losses. A pressure shifted valve timing has lower throttling energy losses, making it the favored solution.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Closeout Report of BNL LDRD 23-050 A Second EIC Detector: Physics Case and Conceptual Design

This document is the closeout report for LDRD 23-050, a type-A LDRD project awarded in FY2022 under the title “ A Second EIC Detector: Physics Case and Conceptual Design ”. The project was motivated by the strong interest within the EIC community in a second general-purpose detector and interaction region, and by the recognition that such a detector is essential to fully exploit the scientific potential of the EIC over its multi-decade lifetime. The key goals of the LDRD were to (i) strengthen the case for a second EIC detector, building on the arguments already articulated in the community Yellow Report, (ii) provide a realistic detector concept that is complementary to the current project detector, ePIC, in terms of physics reach, precision, and control of systematics, and (iii) broaden the overall physics program of the EIC facility. Since a possible second detector is expected to be realized with a delay of several years relative to the first detector, the project explicitly aimed at identifying technologies that are not yet sufficiently mature for ePIC but could be deployed on the later timescale of a second detector, thereby providing genuine complementarity and room for innovation. LDRD 23-050 provided support for two postdoctoral researchers and partial support for several BNL staff members, and the effort ramped up in mid-2023. Work was carried out in close collaboration with colleagues from the EIC project, the Cold QCD and STAR groups at BNL, and with substantial input from members of the EIC User Group and the ePIC Collaboration as reflected in the author list. Over the period from October 2022 through September 2025, this team developed and refined a physics program tailored to a second detector, translated the resulting physics requirements into detector-level performance targets, and explored multiple conceptual detector layouts capable of meeting these targets within the constraints of the second interaction region. The expected outcome, as stated in the original proposal, was a document detailing the physics potential and requirements of a second EIC detector, accompanied by a comprehensive conceptual detector design and an outline of the remaining R&D needed to bring the relevant technologies to maturity. The present report summarizes the progress that has been achieved toward these goals. It consolidates the physics studies, detector concepts, and technology assessments developed under this LDRD and places them in the broader context of worldwide detector R&D. While the evolving priorities of the EIC project and the substantial community effort required to design and construct ePIC naturally limited the scope of EIC-wide activities toward a second detector during this period, the work documented here is intended to provide a foundation and a point of reference for future efforts. Whatever form a second EIC detector may ultimately take, we hope that this report will serve as a useful guide for colleagues who continue to advance this program in the near- and mid-term future.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Measuring the robustness of predictive probability for early stopping in two-group comparisons

We report physical experiments are often expensive and time-consuming. Test engineers must certify the compatibility of aircraft and their weapon systems before they can be deployed in the field, but the testing required is time consuming, expensive, and resource limited. Adopting Bayesian adaptive designs is a promising way to borrow from the successes seen in the clinical trials domain. The use of predictive probability (PP) to stop testing early and make faster decisions is particularly appealing given the aforementioned constraints. Given the high-consequence nature of the tests performed in the national security space, a strong understanding of new methods is required before being deployed. Although PP has been thoroughly studied for binary data, there is less work with continuous data, where many reliability studies are interested in certifying the specification limits of components. A simulation study evaluating the robustness of this approach indicates early stopping based on PP is reasonably robust to minor assumption violations, especially when only a few interim analyses are conducted. The simulation study also compares PP to conditional power, showing its relative strengths and weaknesses. A post-hoc analysis exploring whether release requirements of a weapon system from an aircraft are within specification with desired reliability resulted in stopping the experiment early and saving 33% of the experimental runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Review of Remote Welding and Nondestructive Examination Technologies for the DOE Standard Canister

The U.S. Department of Energy (DOE) manages a wide variety of spent nuclear fuel (SNF) that poses a unique management challenge. To help address this challenge, the DOE Standard Canister (DOESC), designed to remain sealed during handling, storage, transportation, and disposal, was conceptualized as a standardized containment vessel to accommodate DOE-managed SNF. Since 1999, several welding and examination processes have been independently developed for the DOESC’s closure welds. However, neither the DOESC nor these processes have been realized in an operational capacity. This review paper seeks to present and compare previously developed DOESC closure weld, nondestructive examination, and repair processes and technologies. Specific processes developed for the Idaho Spent Fuel Facility, in preparation for the Yucca Mountain geological repository, and the recent Road-Ready Demonstration Project are discussed. Further, specific focus is given to how different operating constraints and the American Society of Mechanical Engineers Boiler and Pressure Vessel Code (BPVC) have driven certain welding and nondestructive examination requirements. Historical DOESC welding and examination strategies are assessed against current regulatory and BPVC requirements. The comparison of welding processes, technologies, and DOESC designs presented in this review paper will inform further construction efforts for other commercial and DOE-managed SNF containments, including the DOESC.

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

Are light curve classification metrics good proxies for SN Ia cosmological constraining power?

Context. When selecting a light curve classifier for use as part of a photometric supernova Ia (SN Ia) cosmological analysis, it is common to make decisions based on metrics of classification performance, such as the contamination within the photometrically classified SN Ia sample, rather than a measure of cosmological constraining power. If the former is an appropriate proxy for the latter, this practice would eliminate the computational expense of a full cosmology forecast in the analysis pipeline design process. Aims. This study tests the assumption that light curve classification metrics are an appropriate proxy for cosmology metrics. Methods. We emulated photometric SN Ia cosmology light curve samples with controlled contamination rates of individual contaminant classes and evaluated each of them under a set of classification metrics. We then derived cosmological parameter constraints from all samples under two common analysis approaches and quantified the impact of contamination by each contaminant class on the resulting cosmological parameter estimates. Results. We observe that cosmology metrics are sensitive to both the contamination rate and the class of the contaminating population, whereas the classification metrics are shown to be insensitive to the latter. Conclusions. Based on these findings, we discourage any exclusive reliance on light curve classification-based metrics for analysis design decisions, which (counterintuitively) include but are not limited to the classifier choice. Instead, we recommend optimising science analysis pipeline design choices using a metric of the information gained about the physical parameters of interest.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Feedforward equilibrium trajectory optimization with GSPulse

One of the common tasks required for designing new plasma scenarios or evaluating capabilities of a tokamak is to design the desired equilibria using a Grad-Shafranov (GS) equilibrium solver. However, most standard equilibrium solvers are time-independent and do not include dynamic effects such as plasma current flux consumption, induced vessel currents, or voltage constraints. Another class of tools, plasma equilibrium evolution simulators, do include time-dependent effects. These are generally structured to solve the forward problem of evolving the plasma equilibrium given feedback-controlled voltages. In this work, we introduce GSPulse, a novel algorithm for equilibrium trajectory optimization, that is more akin to a pulse planner than a pulse simulator. GSPulse includes time-dependent effects and solves the inverse problem: given a user-specified set of target equilibrium shapes, as well as limits on the coil currents and voltages, the optimizer returns trajectories of the voltages, currents, and achievable equilibria. This task is useful for scoping performance of a tokamak and exploring the space of achievable pulses. The computed equilibria satisfy both Grad-Shafranov force balance and axisymmetric circuit dynamics. The optimization is performed by restructuring the free-boundary equilibrium evolution equations into a form where it is computationally efficient to optimize the entire dynamic sequence. GSPulse can solve for hundreds of equilibria simultaneously within a few minutes. GSPulse has been validated against NSTX-U and MAST-U experiments and against SPARC feedback control simulations, and is being used to perform scenario design for SPARC. The computed trajectories can be used as feedforward inputs that are connected to the feedback controller to inform and improve feedback performance. The code for GSPulse is available open-source at github.com/jwai-cfs/GSPulse_public.

equilibrium↗
DOE OSTI2023

BICEP Array: 150 GHz Detector Module Development

The Background Imaging of Cosmic Extragalactic Polarization (BICEP)/Keck (BK) collaboration is currently leading the quest for the highest-sensitivity measurements of the polarized cosmic microwave background (CMB) anisotropies on a degree scale with a series of cryogenic telescopes, of which BICEP Array (BA) is the latest Stage-3 upgrade with a total of ~ 32,000 detectors. The instrument comprises 4 receivers spanning 30–270 GHz, with the low-frequency 30/40 GHz deployed to the South Pole Station in late 2019. The full complement of receivers is forecast to set the most stringent constraints on the tensor-to-scalar ratio r. Building on these advances, the overarching small-aperture telescope concept is already being used as the reference for further Stage-4 experiment design. This paper describes the development of the BICEP Array 150 GHz detector module and its fabrication requirements, with highlights on the high-density time division multiplexing (TDM) design of the cryogenic circuit boards. The low-impedance wiring required between the detectors and the frst stage of superconducting quantum interference device amplifers is crucial to maintaining a stable bias current on the detectors. Here, a novel multilayer FR4 Printed Circuit Board with superconducting traces, capable of reading out up to 648 detectors, is detailed along with its validation tests. An ultra-high-density TDM detector module concept we developed for a CMB-S4-like experiment that allows up to 1920 detectors to be read out is also presented. TDM has been chosen as the detector readout technology for the Cosmic Microwave Background Stage-4 (CMB-S4) experiment based on its proven low-noise performance, predictable costs, and overall maturity of the architecture. The heritage for TDM is rooted in mm- and sub-mm-wave experiments dating back 20 years and has since evolved to support a multiplexing factor of 64x in Stage-3 experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Pyomo.GDP: an ecosystem for logic based modeling and optimization development

We present three core principles for engineering-oriented integrated modeling and optimization tool sets—intuitive modeling contexts, systematic computer-aided reformulations, and flexible solution strategies—and describe how new developments in Pyomo.GDP for Generalized Disjunctive Programming (GDP) advance this vision. We describe a new logical expression system implementation for Pyomo.GDP allowing for a more intuitive description of logical propositions. The logical expression system supports automated reformulation of these logical constraints to linear constraints. We also describe two new logic-based global optimization solver implementations built on Pyomo.GDP that exploit logical structure to avoid “zero-flow” numerical difficulties that arise in nonlinear network design problems when nodes or streams disappear. These new solvers also demonstrate the capability to link to external libraries for expanded functionality within an integrated implementation. We present these new solvers in the context of a flexible array of solution paths available to GDP models. Finally, we present results on a new library of GDP models demonstrating the value of multiple solution approaches.

42 ENGINEERING↗
DOE OSTI2026

Machine learning-enhanced MPC for demand flexibility in small commercial buildings: An experimental study

Small- and medium-sized commercial buildings (SMCBs) represent the majority of U.S. commercial building stock and a significant share of peak electricity demand, yet they often lack centralized building automation systems, representing a significant untapped resource for urban energy management. This infrastructure gap makes advanced control implementation challenging, limiting the potential for widespread demand flexibility. Model Predictive Control (MPC) has shown strong potential for load shifting, peak demand reduction, and cost savings, but its effectiveness is hindered by unmeasured disturbances such as internal heat gains. This paper presents a Hybrid MPC framework that integrates a physics-based gray-box building thermal model, identified using a lumped disturbance (LD) approach, with a machine learning (ML) model for forecasting unmeasured disturbances. The hybrid approach is designed for buildings with multiple individually controlled heat pump and thermostat pairs, common in SMCBs, and aims to optimize coordinated scheduling of multiple heat pumps under dynamic electricity pricing while respecting comfort constraints. The methodology is validated through both simulations of case study buildings and experimental studies at a highly-instrumented test facility. Simulation results show that the Hybrid MPC achieves substantial load shifting and peak demand reduction, approaching the performance of an ideal MPC with perfect disturbance knowledge, and outperforming a conventional MPC without disturbance forecasting. In experiments, the Hybrid MPC reduced daily HVAC energy costs by 8.7%, peak-price time load (load shifting) by 41.7%, and peak demand by 29.2% compared to baseline control, demonstrating comparable benefits to the 11.6% cost savings, 42.9% load shifting, and 23.2% peak reduction of the ideal MPC. These results demonstrate that the proposed hybrid modeling approach can significantly improve MPC performance in real-world SMCB applications without requiring additional disturbance measurements.

Demand Flexibility↗