Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid modeling”

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

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

At least 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Near-UV absorption in very cool DA white dwarfs

The atmospheres of very cool, hydrogen-rich white dwarfs (WDs) (T{sub eff} < 6000 K) are challenging to model because of the increased complexity of the equation of state, chemical equilibrium, and opacity sources in a low-temperature, weakly ionized dense gas. In particular, many models that assume relatively simple models for the broadening of atomic levels and mostly ideal gas physics overestimate the flux in the blue part of their spectra. A solution to this problem that has met with some success is that additional opacity at short wavelengths comes for the extreme broadening of the Lyman α line of atomic H by collisions primarily with H{sub 2}. For the purpose of validating this model more rigorously, we acquired Hubble Space Telescope STIS spectra of eight very cool WDs (five DA and three DC stars). Combined with their known parallaxes, BVRIJHK, and Spitzer IRAC photometry, we analyze their entire spectral energy distribution (from 0.24 to 9.3 μm) with a large grid of model atmospheres and synthetic spectra. We find that the red wing of the Lyman α line reproduces the rapidly decreasing near-UV flux of these very cool stars very well. We determine better constrained values of T{sub eff} and gravity as well as upper limits to the helium abundance in their atmospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2025

CGSim: A Simulation Framework for Large Scale Distributed Computing Environment

Large-scale distributed computing infrastructures such as the Worldwide LHC Computing Grid (WLCG) require comprehensive simulation tools for evaluating performance, testing new algorithms, and optimizing resource allocation strategies. However, existing simulators suffer from limited scalability, hardwired algorithms, lack of real-time monitoring, and inability to generate datasets suitable for modern machine learning approaches. We present CGSim, a simulation framework for large-scale distributed computing environments that addresses these limitations. Built upon the validated SimGrid simulation framework, CGSim provides high-level abstractions for modeling heterogeneous grid environments while maintaining accuracy and scalability. Key features include a modular plugin mechanism for testing custom workflow scheduling and data movement policies, interactive real-time visualization dashboards, and automatic generation of event-level datasets suitable for AI-assisted performance modeling. We demonstrate CGSim’s capabilities through a comprehensive evaluation using production ATLAS PanDA workloads, showing significant calibration accuracy improvements across WLCG computing sites. Scalability experiments show near-linear scaling for multi-site simulations, with distributed workloads achieving 6 × better performance compared to single-site execution. The framework enables researchers to simulate WLCG-scale infrastructures with hundreds of sites and thousands of concurrent jobs within practical time budget constraints on commodity hardware.

Vatsavai, Sairam Sri [Brookhaven National Laborato↗
NASA NTRS1991

A numerical modeling study of the interaction between the tides and the circulation forced by high-latitude plasma convection

A spectral, time-varying thermospheric general circulation model has been used to study the nonlinear interaction at high latitudes between the tides propagating into the thermosphere from below and the circulation induced by magnetospheric forcing and in situ solar heating. The model is discrete in the vertical with 27 layers spaced by half a scale height. In the horizontal, the fields are expanded in a series of spherical harmonics using a triangular truncation at wave number 31, equivalent to a homogeneous global resolution with a minimum wavelength of 1270 km. A hypothetical uniform grid point model would require a horizontal spacing of 417 km to describe the same minimum wavelength. In the high-latitude F region the tides affect the dusk vortex of the neutral flow very little, but the dawn vortex is either suppressed or amplified dependent upon the universal time and tidal phase. In the E region neutral flow, both the dusk and dawn vortices are shifted in local time by the tides, again as a function of universal time and tidal phase. At dusk a nonlinear amplification of the sunward winds occurs for certain combinations of parameters, and at dawn the winds may be completely suppressed. Below 120 km altitude the magnetospheric forcing creates a single cyclonic vortex which is also sensitive to the high-latitude tidal structure.

Mikkelsen, I. S.↗
DOE OSTI2024

Dependencies of Simulated Convective Cell and System Growth Biases on Atmospheric Instability and Model Resolution

Abstract This study evaluates convective cell properties and their relationships with convective and stratiform rainfall within a season‐long convection‐permitting weather research and forecasting simulation over central Argentina using radar, satellite, and radiosonde measurements from the RELAMPAGO‐CACTI field campaign. The simulation slightly underestimates radar‐estimated rainfall over the ∼3.5‐month evaluation period but underestimates stratiform rainfall by 46% and overestimates convective rainfall by 43%. As convective available potential energy (CAPE) increases, the convective rainfall overestimation decreases, but the stratiform rainfall underestimation increases such that the contribution of convective to total rainfall remains constantly high biased by ∼26%. Overestimated convective rainfall arises from the simulation generating 2.6 times more precipitating convective cells (14,299) than observed by radar (5,662) despite similar observed and simulated cell growth processes, with relatively wide cells contributing mostly to excessive convective rainfall. Relatively shallow cells, typically reaching heights of 4–7 km, contribute most to the cell number bias. This cell number bias increases as CAPE decreases, potentially because cells and their updrafts become narrower and more under‐resolved as CAPE decreases. The gross overproduction of precipitating shallow cells leads to overly efficient precipitation and inadequate detrainment of ice aloft, thereby diminishing the formation of robust stratiform rainfall regions. Decreasing model horizontal grid spacing from 3 to 1 or 0.333 km for low (<300 J kg −1 ) and high CAPE (>1,000 J kg −1 ) cases results in minimal change to cell number, depth, and convective‐to‐stratiform partitioning biases. This suggests that improving prediction of these convective properties depends on factors beyond solely increasing model resolution.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2019

Utilization of Airborne Observations to Assess Model Parametrizations of Critical RH Profiles in the Arctic Ocean

Observational paucity is a reality in the Arctic Ocean. This is especially true for near-surface variables such as temperature, moisture, heat fluxes and BL clouds. As a result, modeling has become one of the major avenues for understanding current and future Arctic trends. Reanalyses are frequently used to force global ocean circulation and sea-ice models. But in northern high latitudes, model integrations and reanalyses are known to have large uncertainties in temperature and humidity profiles, and in boundary layer cloudiness. These are common sources of error in the surface radiative budget terms.An important way to diagnose these biases spatially and temporally is by using satellite remote sensing data. However, remotely-sensed observations also have large uncertainties, especially in near-surface temperature and relative humidity profiles.In situ observational studies are important in bridging our knowledge gap in regions such as the Arctic Ocean. Here, we utilize airborne and ship observations from the ARISE, ACME-IV, and ASCOS campaigns to construct critical relative humidity (RH) profiles over the Beaufort Sea. Such profiles are used as parameterization inputs in the NASA GOES global model to derive the total water condensate in a model grid-box, which determines the cloud fraction. Currently, the critical RH profiles are derived by global AIRS data, relaying mostly on mid-latitude regions, which are not necessarily relevant to the Arctic.We derive campaign-wide mean, standard deviation, and critical RH values, for grid size of 50x50 km and altitude bins between 50 to 400 m, covering both open-ocean and sea-ice covered regions. We compare profiles over open ocean and sea-ice, and look at correlations between the observed critical RH values and water condensate (by cloud number concentration) from observations versus the modeled ones. We then input our calculated values of minimal critical RH values into a set of GEOS single column model (SCM) simulations over the ARISE and ASCOS regions and compare the differences between the predicted values of cloud liquid water path (LWP), ice water path (IWP) and surface fluxes with the observed ones under the range of input parameterizations. Finally, we discuss the implications on surface radiative budget predictions in this region.

Segal-Rozenhaimer, Michal↗
DOE OSTI2024

Impact of Systematic Modeling Uncertainties on Kilonova Property Estimation

The precise atomic structure and therefore the wavelength-dependent opacities of lanthanides are highly uncertain. This uncertainty introduces systematic errors in modeling transients like kilonovae and estimating key properties such as mass, characteristic velocity, and heavy metal content. Here, we quantify how atomic data from across the literature as well as choices of thermalization efficiency of r-process radioactive decay heating impact the light curve and spectra of kilonovae. Specifically, we analyze the spectra of a grid of models produced by the radiative transfer code Sedona that span the expected range of kilonova properties to identify regions with the highest systematic uncertainty. Our findings indicate that differences in atomic data have a substantial impact on estimates of lanthanide mass fraction, spanning approximately 1 order of magnitude for lanthanide-rich ejecta, and demonstrate the difficulty in precisely measuring the lanthanide fraction in lanthanide-poor ejecta. Mass estimates vary typically by 25%–40% for differing atomic data. Similarly, the choice of thermalization efficiency can affect mass estimates by 20%–50%. Observational properties such as color and decay rate are highly model dependent. Velocity estimation, when fitting solely based on the light curve, can have a typical error of ∼100%. Atomic data of light r-process elements can strongly affect blue emission. Even for well-observed events like GW170817, the total lanthanide production estimated using different atomic data sets can vary by a factor of ∼6.

Gravitational wave sources↗
DOE OSTI2026

Challenges and Opportunities for Electric Utility Modeling and Asset Valuation Frameworks: Case Study on Valuing New Pumped Storage Hydropower

Asset valuation by electric utilities is becoming increasingly difficult in the rapidly changing electric sector. Rapid deployment of variable generation and inverter-based storage systems along with uncertain demand growth, climate, policies, and other factors create a challenging environment for understanding the value proposition of a new potential asset. This report describes an effort between the Tennessee Valley Authority (TVA) and three U.S. Department of Energy laboratories to perform a detailed review of utility modeling and analysis practices for asset valuation and identify challenges and opportunities for advancing its methods into the future. It focuses on a case study of new potential pumped storage hydropower (PSH) because of growing interest in new PSH capacity to provide energy balancing, firm capacity, and a range of ancillary services. Staff from the DOE labs conducted systematic interviews about current practices in capacity expansion modeling, production-cost modeling, hydrological modeling, and transmission stability modeling while also discussing how scenario analysis is conducted and how models and data are integrated. The effort resulted in a set of model, integration, and scenario recommendations that could be valuable to TVA, other utilities, system operators, and other stakeholders conducting integrated grid analysis. Individual model recommendations suggest exploring computational tradeoffs with detail and resolution across spatiotemporal structure, supply- and demand-side details, transmission overlays, market interactions, and ancillary services. Automated processes to pass data between models and conduct larger scenario suites could also enhance valuation practices by enabling a more consistent study of asset value across a broader range of uncertain future grid conditions where PSH could be particularly valuable. TVA and other industry stakeholders can learn from and adapt applied research-grade methods developed by DOE laboratories and other research institutions to improve decision making and accelerate progress towards a reliable, economic, sustainable energy system.

13 HYDRO ENERGY↗
NASA NTRS1994

Navier-Stokes Solutions of Rotor and Rotor-Body Flows

This paper will review the advances made recently in the Navier-Stokes CFD methods to simulate aerodynamics and aeroacoustics of helicopter rotors and rotor-body flows. Although a complete flowfield simulation of full helicopter is currently not feasible with these methods, impressive gains have been made in analyzing individual components of this complex problem in a very detailed manner. The use of the state-of-the-art numerical algorithms in solution methods, in conjunction with powerful supercomputers, like the Cray-2, have enabled noticeable progress to be made in modeling viscous-inviscid interactions, blade-vortex interactions, tip-vortex: simulation and wake effects, as well as high speed impulsive noise in hover and forward flight for isolated rotor blades. This paper will critically evaluate the presently available Euler and Navier-Stokes methods, both finite-difference and finite volume methods using structured and unstructured grids for helicopter applications for accuracy, suitability, and computational efficiency. The review will also include the recent progress made using overset grids to model rotor-body flows. All the material for this review will be drawn from the published material shown below.

Srinivasan, G. R.↗
DOE OSTI2014

Pair instability supernovae of very massive population III stars

Numerical studies of primordial star formation suggest that the first stars in the universe may have been very massive. Stellar models indicate that non-rotating Population III stars with initial masses of 140-260 M {sub ☉} die as highly energetic pair-instability supernovae. We present new two-dimensional simulations of primordial pair-instability supernovae done with the CASTRO code. Our simulations begin at earlier times than previous multidimensional models, at the onset of core contraction, to capture any dynamical instabilities that may be seeded by core contraction and explosive burning. Such instabilities could enhance explosive yields by mixing hot ash with fuel, thereby accelerating nuclear burning, and affect the spectra of the supernova by dredging up heavy elements from greater depths in the star at early times. Our grid of models includes both blue supergiants and red supergiants over the range in progenitor mass expected for these events. We find that fluid instabilities driven by oxygen and helium burning arise at the upper and lower boundaries of the oxygen shell ∼20-100 s after core bounce. Instabilities driven by burning freeze out after the SN shock exits the helium core. As the shock later propagates through the hydrogen envelope, a strong reverse shock forms that drives the growth of Rayleigh-Taylor instabilities. In red supergiant progenitors, the amplitudes of these instabilities are sufficient to mix the supernova ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

The Second Wind Forecast Improvement Project (WFIP 2): General Overview

In 2015 the U.S. Department of Energy initiated a four-year study, the second Wind Forecast Improvement Project (WFIP 2), to improve the representation of boundary-layer physics and related processes in mesoscale models for better wind and wind power forecasts. This goal challenges numerical weather prediction models in complex terrain due to inherent assumptions underlying their boundary-layer parameterizations. The WFIP 2 effort involved the wind industry, universities, NOAA, and DOE's national laboratories in an integrated observational and modeling study. Observations spanned 18-months to assure a full annual cycle of continuously recorded observations from remote-sensing and in situ measurement systems. The study area comprised the Columbia Basin of eastern Washington and Oregon, containing more than 4 GW of installed wind capacity; nests of observational systems captured important atmospheric scales from mesoscale to NWP sub-grid scale. Model improvements targeted NOAA's High-Resolution Rapid Refresh (HRRR) model to facilitate transfer of improvements to NWS operational forecast models. This paper describes the general WFIP 2 scope and objectives, the particular scientific challenges of improving wind forecasts in complex terrain, early successes of the project, and an integrated approach to archiving observations and model output. It provides an introduction for a set of more detailed BAMS papers addressing WFIP 2 observational science, modeling challenges and solutions, practical transfer of boundary layer advances to industry use, and advances in coupling improved mesoscale models to microscale models that can represent interactions between wind plants and the atmosphere.

WFIP2, boundary layer, complex terrain, forecastin↗
NASA NTRSDate not supplied

The System Modeling and Analysis of Resiliency in STEReO (SMARt-STEReO)

Wildfire emergency response has remained rooted in relatively low-tech solutions for coordination between ground and aerial assets. These low-tech solutions are robust for the remote environments in which wildfires are usually fought, but limit strategic cross-organizational support and the ability to deploy and effectively utilize aerial assets. As aircraft become more advanced and new technology, including drones, become available to firefighters, a new, more modern method of asset coordination is needed. NASA is working on a project called ‘Scalable Traffic Management for Emergency Response Operations’ (STEReO) to integrate unmanned aerial systems (UAS)and UAS traffic management (UTM)into wildfire response. STEReO’s goals include simplifying the coordination of aerial assets, improving the existing UAS framework, and increasing the role of additional autonomous systems to reduce human risk and to increase system resilience. This paper describes the development of the ‘System Modeling and Analysis of Resiliency in STEReO’ (SMARt-STEReO) project, which aims to model wildfire response and to quantify the additional system resilience that STEReO technology provides firefighters. This paper verifies SMARt-STEReO and defines its scope; it includes experimental and statistical analysis of the impact that the addition of UAS has on both performance metrics and also on performance resiliency response to a given fault. SMARt-STEReO is a grid-based model of fire propagation that incorporates varying crew responses. Through the use of a Python package called ‘fmdtools’, the model easily allows for the addition of faults to the system. These faults allow analysts to investigate various response parameters. Factors including terrain, fuel type and wind speed can be modified to affect the fire propagation; additionally, the number of ground crews, engines, fixed wing aircraft, helicopters, and UAS can be changed to affect the crew response. The communication lines between actors mimic those used in real life situations. This paper explains the development of SMARt-STEReO including background research, verification and validation, and preliminary experimental analysis of system resilience to both a minor and major fault in systems with and without UAS.

Resiliency↗
NASA NTRSDate not supplied

Improving Simulation of Irrigation in a Global Land Surface Model

Studies have shown that irrigated agriculture has substantial effects on local to regional hydrologic and atmospheric processes. These effects include increased evapotranspiration and near surface humidity, decreased air and ground temperature, augmented downwind precipitation, enhanced recharge, depletion of source water stocks, and overall changes to the regional water budget. However, developers of global offline and coupled modeling systems have been slow to incorporate irrigation. This presentation will describe progress towards simulating irrigation worldwide in the Land Information System (LIS) and Global Land Data Assimilation System (GLDAS). We have developed and refined the ability to simulate sprinkler irrigation, furrow irrigation, paddy irrigation, and drip irrigation. We have reconciled multiple global and national irrigation, land cover, and crop type databases in an effort to (1) distinguish irrigated and rainfed agriculture, (2) simulate the proper distribution of crops (as sub-grid tiles) in each model grid pixel, and (3) apply the most likely irrigation method in each tile based on the crop type and country or state level information on the relative usage of sprinkler, flood, and drip irrigation. Using agricultural water withdrawal data, as available, for calibration and evaluation, we show that the implemented approach is promising but that much fine tuning remains to be done. In particular, irrigation model parameters that produce good results in one location often produce less satisfying results in others.

Matthew Rodell↗
DOE OSTI2018

Sensitivities of Simulated Satellite Views of Clouds to Subgrid-Scale Overlap and Condensate Heterogeneity

Satellite simulators are often used to account for limitations in satellite retrievals of cloud properties in comparisons between models and satellite observations. The purpose of this framework is to enable more robust evaluation of model cloud properties, so that differences between models and observations can more confidently be attributed to model errors. A critical step in this process is accounting for the difference between the spatial scales at which cloud properties are retrieved with those at which clouds are simulated in global models. In this study, we create a series of sensitivity tests using 4-km global model output from the multiscale modeling framework to evaluate the sensitivity of simulated satellite retrievals to common assumptions about cloud and precipitation overlap and condensate variability used in climate models whose grid spacing is many tens to hundreds of kilometers. We find the simulated retrievals are sensitive to these assumptions. Using maximum-random overlap with homogeneous cloud and precipitation condensate leads to errors in Multiangle Imaging Spectroradiometer and International Satellite Cloud Climatology Project-simulated cloud cover and in CloudSat-simulated radar reflectivity that are significant compared to typical differences between the model simulations and observations. A more realistic treatment of unresolved clouds and precipitation is shown to substantially reduce these errors. The sensitivity to these assumptions underscores the need for the adoption of more realistic subcolumn treatments in models and the need for consistency among subcolumn assumptions between models and simulators to ensure that simulator-diagnosed errors are consistent with the model formulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Multiscale Electricity Modeling for Evaluating Carbon Capture and Sequestration Technologies (MEME-CCS)

This effort employs and adapts an existing, rigorous multiscale electricity modeling platform at the National Renewable Energy Laboratory (NREL) to evaluate carbon capture and sequestration (CCS) and negative emissions technologies (NET) from the ARPA-E FLECCS program. NREL's modeling platform includes the Regional Energy Deployment System (ReEDS) electric sector capacity expansion model, which projects future electricity generation mixes at sub-state-level resolution that are downscaled to the unit-level to enable hourly, zonal or nodal electricity production cost modeling in the PLEXOS model. The resulting hourly price data from PLEXOS is provided to technology developers under the ARPA-E FLECCS program to enable technology-specific economic analysis. The ReEDS-PLEXOS modeling suite is well-established for examining electric sector futures with high renewable energy penetrations, energy storage, electrification, and distributed generation. The key advancement proposed herein utilizes collaboration with CCS experts at the University of Wyoming and the FLECCS teams to create innovative methods for representing CCS and NET in the ReEDS and PLEXOS models. Expanded technology options, new operational parameterizations, and detailed data defining CO2 capture, transportation, and storage systems are being integrated into these models to allow an unprecedented combination of scope and resolution for exploring the future of CCS and NET. The resulting capabilities will take advantage of high-performance computing resources to permit wide-ranging scenario analysis to assess CCS and NET deployment under alternative CO2 prices, fossil fuel prices, electricity demand growth, and other electric sector characteristics. Final outcomes will include publicly available hourly grid operation and price data for any U.S. region of interest along with open-access capacity expansion tools for evaluating CCS/NET systems. These products will help an emerging CCS/NET industry in the United States by providing economics-driven guidance to technology developers while informing policy and investment decisions in the public and private sectors.

air capture↗
DOE OSTI2015

THE ABUNDANCES OF LIGHT NEUTRON-CAPTURE ELEMENTS IN PLANETARY NEBULAE. III. THE IMPACT OF NEW ATOMIC DATA ON NEBULAR SELENIUM AND KRYPTON ABUNDANCE DETERMINATIONS

The detection of neutron(n)-capture elements in several planetary nebulae (PNe) has provided a new means of investigating s-process nucleosynthesis in low-mass stars. However, a lack of atomic data has inhibited accurate trans-iron element abundance determinations in astrophysical nebulae. Recently, photoionization (PI) and recombination data were determined for Se and Kr, the two most widely detected n-capture elements in nebular spectra. We have incorporated these new data into the photoionization code Cloudy. To test the atomic data, numerical models were computed for 15 PNe that exhibit emission lines from multiple Kr ions. We found systematic discrepancies between the predicted and observed emission lines that are most likely caused by inaccurate PI and recombination data. These discrepancies were removed by adjusting the Kr{sup +}–Kr{sup 3+} PI cross sections within their cited uncertainties and the dielectronic recombination rate coefficients by slightly larger amounts. From grids of models spanning the physical conditions encountered in PNe, we derive new, broadly applicable ionization correction factor (ICF) formulae for calculating Se and Kr elemental abundances. The ICFs were applied to our previous survey of near-infrared [Kr iii] and [Se iv] emission lines in 120 PNe. The revised Se and Kr abundances are 0.1–0.3 dex lower than former estimates, with average values of [Se/(O, Ar)] = 0.12 ± 0.27 and [Kr/(O, Ar)] = 0.82 ± 0.29, but correlations previously found between their abundances and other nebular and stellar properties are unaffected. We also find a tendency for high-velocity PNe that can be associated with the Galactic thick disk to exhibit larger s-process enrichments than low-velocity PNe belonging to the thin-disk population.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

cerf: A Python package to evaluate the feasibility and costs of power plant siting for alternative futures

Long-term electric power sector planning and capacity expansion is a key area of interest to stakeholders across a wide range of organizations because it helps in making informed decisions about investments in infrastructure within the context of potential future vulnerabilities under various natural and human stressors. Future power plant siting costs will depend on a number of factors including the characteristics of the electricity capacity expansion and electricity demand (e.g., fuel mix of future electric power capacity, and the magnitude and geographic distribution of electricity demand growth) as well as the geographic location of power plants. Electricity technology capacity expansion plans modeled to represent alternate future conditions meeting a set of scenario assumptions are traditionally compared against historical trends which may not be consistent with current and future conditions. We present the `cerf` Python package (a.k.a., the Capacity Expansion Regional Feasibility model) which helps evaluate the feasibility and structure of future, scenario-driven electricity capacity expansion plans by siting power plants in areas that have been deemed the least cost option while considering dynamic future conditions. We can use `cerf` to gain insight to research topics such as: 1) under what conditions future projected electricity expansion plans from models such as GCAM-USA are possible to achieve, 2) where and which on-the-ground barriers to siting (e.g., protected areas, cooling water availability) may influence our ability to achieve certain expansions, and 3) how electricity infrastructure build-outs and value may evolve into the future when considering locational marginal pricing (LMP) based on the supply and demand of electricity from a grid operations model.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Enabling Extended-Term Simulation of Power Systems with High PV Penetration. Final Report

This is the final Technical Report for DOE-SETO Project Award # DE-EE0036461. The goal of this project is to advance the understanding of the grid impact of high penetration of photovoltaic (PV) generation by developing novel numerical methods to solve the differential algebraic equations (DAEs) that define power systems. This will overcome the limitations of current software packages – namely that they only consider fast dynamics over brief time periods. The work presented in this final project report covers results over the entire period of the project. This includes results on model development, code development for the PST repository, datasets in the PST repository, algorithm development and results from variable time-step simulations, development and results from multirate simulations, and sensitivity analysis of key parameter in variable time-step methods. In addition, this report discusses project outreach activities to stakeholders, and a summary of project products. Also covered in this final report is the writing of two conference papers (one of which has already been accepted) and a journal paper. In addition, the updating of two inverter models (both grid forming and grid following) to be compatible with the latest version of PST software is discussed.

14 SOLAR ENERGY↗