Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Measuring Flexibility through Reduction Potential

As electrical loads, electric vehicles (EVs) are distinguished in many cases by their substantial flexibility. Harnessing this flexibility, however, requires a precise characterization of its timing and degree. Current approaches to flexibility characterization tend to be limited either by the complexity of their representation or their inability to capture aggregate flexibility across multiple vehicles. In this paper, we present the reduction potential matrix, a novel metric capturing fleet-level EV load flexibility while being both straightforward to calculate and intuitive to interpret. We use the metric to quantify flexibility for two operationally distinct commercial groups: freight vehicles and transit buses. While we find that both vehicle groups have substantial flexibility, its timing, magnitude, and duration varies across the groups. We describe how this variability translates into variation in the each group's role as a grid resource, and conclude with a discussion of how system planners, fleet operators, and other stakeholders can use the novel metric to cultivate symbiosis between vehicles and the grid.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Small-Signal Stability Constrained Optimal Power Flow of Inverter-Dominated Power Systems with Flexible Operation Mode Selection

Given the intermittence and low inertia nature of inverter-based resources (IBRs), modern power systems with high penetration of IBRs challenge the conventional optimal power flow (OPF) analysis and the system may experience unexpected failures if stability constraints are not incorporated. This study proposes a small-signal stability-constrained OPF (SSSC-OPF) with flexible operation mode selection between grid-forming (GFM) and grid-following (GFL) modes for IBRs to address these challenges. The approach aims to maintain system stability with a sufficient stability margin while minimizing operation costs. The effectiveness of the proposed method is validated through extensive case studies on the IEEE 14-bus system. The results demonstrate that the proposed method is able to support system-level power flow analysis, reduce generation costs, and ensure stability under various disturbances.

grid-following↗
DOE OSTI2026

Optimal Co-Design of Integrated Thermal-Electrical Networks and Control Systems for Grid-interactive Efficient District (GED) Energy Systems

This project advances a unified, open-source framework for the optimal co-design of thermal, electrical, and control systems in grid-interactive efficient districts (GEDs). As communities integrate growing levels of distributed energy resources, traditional approaches that model thermal and electrical networks independently lead to reduced efficiency, limited flexibility, and missed opportunities for coordinated operation. To address these challenges, the research team developed a comprehensive suite of physics-based models, control algorithms, and software tools that enable holistic simulation, optimization, and demonstration of district-scale energy systems.

14 SOLAR ENERGY↗
DOE OSTI2026

Dynamic Transmission Line Switching Amid Wildfire-Prone Weather Under Decision-Dependent Uncertainty

During dry and windy seasons, environmental conditions significantly increase the risk of wildfires, exposing power grids to disruptions caused by transmission line failures. Wildfire propagation exacerbates grid vulnerability, potentially leading to prolonged power outages. To address this challenge, we propose a multistage optimization model that dynamically adjusts transmission grid topology in response to wildfire propagation, aiming to develop an optimal response policy. By accounting for decision-dependent uncertainty, where line survival probabilities depend on usage, we employ distributionally robust optimization to model uncertainty in line survival distributions. We adapt the stochastic nested decomposition algorithm and derive a deterministic upper bound for its finite convergence. To enhance computational efficiency, we exploit the Lagrangian dual problem structure for a faster generation of Lagrangian cuts. Using realistic data from the California transmission grid, we demonstrate the superior performance of dynamic response policies against two-stage alternatives through a comprehensive case study. In addition, after solving the multistage formulation, we construct easy-to-implement policies that significantly reduce computational burden while maintaining good performance in real-time deployment. History: Accepted by Russell Bent, Area Editor for Network Optimization: Algorithms and Applications. Funding: This work was supported by the U.S. Department of Energy, Office of Electricity [Grant DE-AC02-05CH11231]. The work of R. Jiang was supported in part by the U.S. National Science Foundation, Division of Electrical, Communications and Cyber Systems [Grant ECCS-1845980] and the U.S. Air Force Office of Scientific Research [Grant FA9550-23-1-0323]. Supplemental Material: The software that supports the findings of this study is available within the paper and its Supplemental Information ( https://pubsonline.informs.org/doi/suppl/10.1287/ijoc.2025.1210 ) as well as from the IJOC GitHub software repository ( https://github.com/INFORMSJoC/2025.1210 ). The complete IJOC Software and Data Repository is available at https://informsjoc.github.io/ .

Estrada-Garcia, Juan-Alberto↗
NASA NTRS1982

Multilayer, Front-Contact Grid for Solar Cells

Proposed multilayer, front-contact grid structure for solar cells optimizes collection of photogenerated current with minimum power losses. It is constructed of several layers of conducting grids. With multilayer concept, peak efficiency can occur at higher output-power levels. Because of this, higher solar concentrations can be applied to solar-cell arrays.

Milnes, A. G.↗
DOE OSTI2026

SODAs: sparse optimization for the discovery of differential and algebraic equations

Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.

DAE↗
NASA NTRS2003

Computational Analysis of Cryogenic Flow Through a Control Valve

The initial efforts to develop the capability to model valves used in rocket engine component testing at Stennis Space Center are documented. An axisymmetric model of a control valve with LN2 as the working fluid was developed. The goal was to predict the effect of change in the plug/sear region of the valve prior to testing. The valve flow coefficient was predicted for a range of plug positions. Verification of the calculations was carried out to quantify the uncertainty in the numerical answer. The modeled results compared well qualitatively to experimental trends. Additionally, insights into the flow processes in the valve were obtained. Benefits from the verification process included the ability to use coarser grids and insight into ways to reduce computational time by using double precision accuracy and non-integer grid ratios. Future valve modeling activities will include shape optimization of the valve/seat region and dynamic grid modeling.

Danes, Russell↗
NASA NTRS1995

Rapid Airplane Parametric Input Design (RAPID)

RAPID is a methodology and software system to define a class of airplane configurations and directly evaluate surface grids, volume grids, and grid sensitivity on and about the configurations. A distinguishing characteristic which separates RAPID from other airplane surface modellers is that the output grids and grid sensitivity are directly applicable in CFD analysis. A small set of design parameters and grid control parameters govern the process which is incorporated into interactive software for 'real time' visual analysis and into batch software for the application of optimization technology. The computed surface grids and volume grids are suitable for a wide range of Computational Fluid Dynamics (CFD) simulation. The general airplane configuration has wing, fuselage, horizontal tail, and vertical tail components. The double-delta wing and tail components are manifested by solving a fourth order partial differential equation (PDE) subject to Dirichlet and Neumann boundary conditions. The design parameters are incorporated into the boundary conditions and therefore govern the shapes of the surfaces. The PDE solution yields a smooth transition between boundaries. Surface grids suitable for CFD calculation are created by establishing an H-type topology about the configuration and incorporating grid spacing functions in the PDE equation for the lifting components and the fuselage definition equations. User specified grid parameters govern the location and degree of grid concentration. A two-block volume grid about a configuration is calculated using the Control Point Form (CPF) technique. The interactive software, which runs on Silicon Graphics IRIS workstations, allows design parameters to be continuously varied and the resulting surface grid to be observed in real time. The batch software computes both the surface and volume grids and also computes the sensitivity of the output grid with respect to the input design parameters by applying the precompiler tool ADIFOR to the grid generation program. The output of ADIFOR is a new source code containing the old code plus expressions for derivatives of specified dependent variables (grid coordinates) with respect to specified independent variables (design parameters). The RAPID methodology and software provide a means of rapidly defining numerical prototypes, grids, and grid sensitivity of a class of airplane configurations. This technology and software is highly useful for CFD research for preliminary design and optimization processes.

Smith, Robert E.↗
DOE OSTI2024

Electrolyte and Cutoff Potential Effects on Cycle Life of Li4Ti5O12/LiNi0.9Mn0.1O2 Batteries for Behind-the-Meter Storage Applications

Behind-the-Meter Storage (BTMS) is a stationary battery energy storage system that is connected to the electrical distribution system on the customer's side of the utility's service meter. BTMS systems are used to store electrical energy from the grid as well as inconstant, renewable energy, such as local solar and wind generation. A successful BTMS system will allow the customer to pair their energy generation and storage to optimize electrical consumption from the grid, improving reliability and minimizing cost. For BTMS applications, batteries must be designed and optimized with different set of criteria from other leading segments of the Li-ion battery market, like transportation, due the system being stationary and proximal to the residential or commercial building it's benefitting. BTMS applications prioritize safety, cost (low/no-critical materials), reliability (20-year calendar life), and durability (10,000 cycle life), while having the ability to (minimally) compromise energy density and rate capability. Lithium titanate (Li4Ti5O12-, LTO) is a promising anode candidate for BTMS applications due to its high safety and capacity retention, while maintaining a reasonable 160 mAhg-1 reversable capacity and composition of relatively abundant materials. (1) Specifically, LTO has a high working voltage which helps to prevent Li dendrite formation, improving safety. Furthermore, LTO also has negligible lithiation-based volume change, leading to less mechanical pulverization, or loss of active material, upon cycling. For the cathode, materials with little or no Co are of high interest due to the high cost and low abundance of Co. LiMn2O4 (LMO) has been paired with LTO for BTMS applications in the past due to its safety, low cost (abundancy), and reasonably high operating voltage. (2-4) However, the low capacity of LMO limits energy density and specific energy. While not the highest priority for BTMS applications, increasing energy density will enable deployment in space constrained BTMS applications and decrease total cost. LiNi0.9Mn0.1O2 (LN-MO) is a recently developed material with promise due to its high operating voltage and relatively low price. (5) However, Ni-rich layered oxides, including LNMO, tend to struggle with capacity retention during high-voltage cycling due to mechanical pulverization, irreversible phase transitions, and unstable solid-electrolyte interphase. The study presented here focuses on building an understanding of how electrolyte solvent and varied cutoff potentials will impact the cycle life of LTO/LN-MO cells. Specifically, a comparison is provided between ethylene carbonate (EC), ethyl methyl carbonate (EMC), fluoroethylene carbonate (FEC), and Gen2 electrolyte solvents with 1M Lithium hexafluorophosphate (LiPF6) salt, cycling to two upper termination potentials, 2.6V and 2.7V. Electrochemical testing and diagnostics (e.g., differential capacity analysis, area specific impedance, constant voltage hold, and rate capability) and post-mortem characterization will be used to understand the aging behavior and failure mechanisms of the 8 cell combinations (four electrolytes and two voltage cutoffs). Cells with FEC electrolyte showed a lower initial capacity compared to cells with Gen2, EMC, and EC cycling at both voltages; however, the cells with FEC showed consistent trends in capacity retention with 2.6V and 2.7V termination potentials, while the cells with the other electrolytes showed much higher rates of capacity loss when cycling to the higher voltage. These results indicate that FEC may play a role in improving durability of high-voltage, Ni-rich electrode systems for use in high-cycle applications, such as BTMS.

electrolyte↗
DOE OSTI2025

MTUQ: a framework for estimating moment tensors, point forces, and their uncertainties

SUMMARY We introduce MTUQ, an open-source Python package for seismic source estimation and uncertainty quantification, emphasizing flexibility and operational scalability. MTUQ provides MPI-parallelized grid search and global optimization capabilities, compatibility with 1-D and 3-D Green’s function database formats, customizable data processing, C-accelerated waveform and first-motion polarity misfit functions, and utilities for plotting seismic waveforms and visualizing misfit and likelihood surfaces. Applicability to a range of full- and constrained-moment tensor, point force, and centroid inversion problems is possible via a documented application programming interface, accompanied by example scripts and integration tests. We demonstrate the software using three different types of seismic events: (1) a 2009 intraslab earthquake near Anchorage, Alaska; (2) an episode of the 2021 Barry Arm landslide in Alaska; and (3) the 2017 Democratic People’s Republic of Korea underground nuclear test. With these events, we illustrate the well-known complementary character of body waves, surface waves, and polarities for constraining source parameters. We also convey the distinct misfit patterns that arise from each individual data type, the importance of uncertainty quantification for detecting multimodal or otherwise poorly constrained solutions, and the software’s flexible, modular design.

58 GEOSCIENCES↗
DOE OSTI2024

Adaptive Online Model Update Algorithm for Predictive Control in Networked Systems

In this article, we introduce an adaptive on-line model update algorithm designed for predictive control applications in networked systems, particularly focusing on power distribution systems. Unlike traditional methods that depend on historical data for offline model identification, our approach utilizes real-time data for continuous model updates. This method integrates seamlessly with existing online control and optimization algorithms and provides timely updates in response to real-time changes. This methodology offers significant advantages, including a reduction in the communication network bandwidth requirements by minimizing the data exchanged at each iteration and enabling the model to adapt after disturbances. Furthermore, our algorithm is tailored for non-linear convex models, enhancing its applicability to practical scenarios. The efficacy of the proposed method is validated through a numerical study, demonstrating improved control performance using a synthetic IEEE test case.

data-driven model predictive control↗
DOE OSTI2025

Comparison of nested geometry treatments within GPU-based Monte Carlo neutron transport simulations of fission reactors

Monte Carlo (MC) neutron transport provides detailed estimates of radiological quantities within fission reactors. This involves tracking individual neutrons through a computational geometry. CPU-based MC codes use multiple polymorphic tracker types with different tracking algorithms to exploit the repeated configurations of reactors, but virtual function calls have high overhead on the GPU. The Shift MC code was modified to support GPU-based tracking with three strategies: dynamic polymorphism with virtual functions, static polymorphism, and a single tracker type with tree-based acceleration. On the Frontier supercomputer these methods achieve 77.8%, 91.2%, and 83.4%, respectively, of the tracking rate obtained using a specialized tracker optimized for rectilinear-grid-based reactors. This indicates that all three methods are suitable for typical reactor problems in which tracking does not dominate runtime. The flexibility of the single tracker method is highlighted with a hexagonal-grid microreactor problem, performed without hexagonal-grid-specific tracking routines, providing a 2.19× speedup over CPU execution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Recovery Simulator and Analysis Formulation: Mathematical Framework for Enhanced Resilience and Resource Allocation

This report introduces recovery simulator and analysis (RSA), a framework aimed at enhancing the resilience of electrical grids post-disruption. The RSA model leverages an optimization problem formulation that focuses on maximizing the load served (or optionally customers served) through a coordinated and cooptimized recovery of non-black start generation, transmission lines, feeders and substations subject to labor budget constraints. By integrating advanced linear programming techniques, the simulator selects efficient reocovery pathways, optimizing both short-term and long-term grid recovery strategies. The mathematical framework guides decision-making through a comprehensive evaluation of potential recovery actions, factoring in the trade-offs between labor constraints and load (or optionally customer) restoration efficacy. This enables grid operators to simulate diverse outage scenarios and delineate optimal recovery pathways, thereby prioritizing critical repair tasks and ensuring resource allocation is both economical and effective. The intended use case of RSA is to allow planners to explore many recovery scenarios quickly and determine assets most critical across a wide range of scenarios, and therefore strong candidates for hardening or additional investment. RSA might also be used in an operational setting, following a single event, for exploring efficient recovery pathways.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1984

18-percent efficient terrestrial silicon solar cells

Silicon solar cells are described which operate at energy conversion efficiencies in excess of 18 percent under standard terrestrial test conditions (AM1.5, 100 mW/sq cm, 28 C). These are believed to be the most efficient silicon cells reported to date. The high efficiency is a result of the combination of high open-circuit voltage due to the careful attention paid to passivation of the top surface of the cell; high fill factors due to the high open-circuit voltage and low parasitic resistance losses; and high short-circuit current due to the use of shallow diffusions, a low grid coverage, and an optimized double-layer antireflection coating.

Blakers, A. W.↗
NASA NTRS1984

High-efficiency silicon solar cells

Silicon solar cells are described which operate at energy conversion efficiencies independently measured at 18.7 percent under standard terrestrial test conditions (AM1.5, 100 mW/sq cm, 28 C). These are apparently the most efficient silicon cells fabricated to date. The high-efficiency results from a combination of high open-circuit voltage due to the careful attention paid to the passivation of the top surface of the cell, high fill factor due to the high open-circuit voltage and low parasitic resistance losses, and high short-circuit current density due to the use of shallow diffusions, a low grid coverage, and an optimized double layer antireflection coating.

Green, M. A.↗
NASA NTRS1989

Computational Fluid Dynamics Program at NASA Ames Research Center

The Computational Fluid Dynamics (CFD) Program at NASA Ames Research Center is reviewed and discussed. The technical elements of the CFD Program are listed and briefly discussed. These elements include algorithm research, research and pilot code development, scientific visualization, advanced surface representation, volume grid generation, and numerical optimization. Next, the discipline of CFD is briefly discussed and related to other areas of research at NASA Ames including experimental fluid dynamics, computer science research, computational chemistry, and numerical aerodynamic simulation. These areas combine with CFD to form a larger area of research, which might collectively be called computational technology. The ultimate goal of computational technology research at NASA Ames is to increase the physical understanding of the world in which we live, solve problems of national importance, and increase the technical capabilities of the aerospace community. Next, the major programs at NASA Ames that either use CFD technology or perform research in CFD are listed and discussed. Briefly, this list includes turbulent/transition physics and modeling, high-speed real gas flows, interdisciplinary research, turbomachinery demonstration computations, complete aircraft aerodynamics, rotorcraft applications, powered lift flows, high alpha flows, multiple body aerodynamics, and incompressible flow applications. Some of the individual problems actively being worked in each of these areas is listed to help define the breadth or extent of CFD involvement in each of these major programs. State-of-the-art examples of various CFD applications are presented to highlight most of these areas. The main emphasis of this portion of the presentation is on examples which will not otherwise be treated at this conference by the individual presentations. Finally, a list of principal current limitations and expected future directions is given.

Holst, Terry L.↗
NASA NTRS1998

Acoustics and Aeroperformance of Nozzles With Screwdriver Shaped and Axisymmetric Plugs

The recent experimental and numerical tests of corrugated nozzles have shown some acoustic and thrust benefits relative to traditional round nozzles. For example, a Bluebell nozzle which was obtained by 3D nozzle design incorporating a corrugated cross section nozzle shape with a sinusoidal lip line nozzle edge, can provide an acoustic benefit up to 4dB with about a 1% thrust augmentation. In references, this effect was explained as being the result of the corrugated design producing more efficient mixing of the exhausted jet with ambient air. Based on this argument, the authors have proposed the application of this concept for a centerbody (plug) which can form several vortices downstream from the centerbody. Several different corrugated designs are proposed and described in detail in this paper. The main design is a Screwdriver shaped centerbody or plug (SCR) which was tested experimentally and numerically. The acoustic tests were conducted in the anechoic chamber of the Central AeroHydrodynamics Institute (TsAGI, Moscow) under Civilian Research and Development Foundation (CRDF) grant. These experiments have shown an essential acoustic benefit of about 10-13% with the application of the co-annular nozzles by comparison with the reference round nozzle with the same mass flow rate. However, the expected acoustic benefits with the application of the 4-petal Screwdriver shaped centerbody were not obtained by comparison with the reference axisymmetric centerbody (CON) having the same length and the same cross section areas at the same distance from the nozzle throat. Moreover, for some angles (Theta = 60 deg and 90 deg) noise increase was observed (about 1-3%). These tests will be continued with the goal of obtaining better acoustic results. In particular, acoustic characteristics are hoped to be improved by moving t lie centerbody into the nozzle and using penetrable walls for the SCR and/or for the main nozzle. Preliminary results for such approach are very promising. Aeroperformance effects were analyzed numerically. The nozzle thrust, calculations were based on a full Navier--Stokes equations solver (NSE), and both full and marching Euler codes: CFL3D, CRAFT, and Krayko-Godunov . Grid preparation and its optimization were conducted using GRIDGEN and our own codes. The general conclusion of this numerical analysis is some thrust loss with the application of SCR design (about 1 - 1.5%). But, again, some constructive features of SCR design give some promising perspectives for its application in aviation and domestic industries.

Gilinsky, M.↗
NASA NTRS1995

CFD-Based Design Optimization Tool Developed for Subsonic Inlet

The traditional approach to the design of engine inlets for commercial transport aircraft is a tedious process that ends with a less-than-optimum design. With the advent of high-speed computers and the availability of more accurate and reliable computational fluid dynamics (CFD) solvers, numerical optimization processes can effectively be used to design an aerodynamic inlet lip that enhances engine performance. The designers' experience at Boeing Corporation showed that for a peak Mach number on the inlet surface beyond some upper limit, the performance of the engine degrades excessively. Thus, our objective was to optimize efficiency (minimize the peak Mach number) at maximum cruise without compromising performance at other operating conditions. Using a CFD code NPARC, the NASA Lewis Research Center, in collaboration with Boeing, developed an integrated procedure at Lewis to find the optimum shape of a subsonic inlet lip and a numerical optimization code, ADS. We used a GRAPE-based three-dimensional grid generator to help automate the optimization procedure. The inlet lip shape at the crown and the keel was described as a superellipse, and the superellipse exponents and radii ratios were considered as design variables. Three operating conditions: cruise, takeoff, and rolling takeoff, were considered in this study. Three-dimensional Euler computations were carried out to obtain the flow field. At the initial design, the peak Mach numbers for maximum cruise, takeoff, and rolling takeoff conditions were 0.88, 1.772, and 1.61, respectively. The acceptable upper limits on the takeoff and rolling takeoff Mach numbers were 1.55 and 1.45. Since the initial design provided by Boeing was found to be optimum with respect to the maximum cruise condition, the sum of the peak Mach numbers at takeoff and rolling takeoff were minimized in the current study while the maximum cruise Mach number was constrained to be close to that at the existing design. With this objective, the optimum design satisfied the upper limits at takeoff and rolling takeoff while retaining the desirable cruise performance. Further studies are being conducted to include static and cross-wind operating conditions in the design optimization procedure. This work was carried out in collaboration with Dr. E.S. Reddy of NYMA, Inc.

Source record↗