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 523 records · Page 29

Evaluating Bay Area Methane Emission Inventory

As a regulatory agency, evaluating and improving estimates of methane (CH4) emissions from the San Francisco Bay Area is an area of interest to the Bay Area Air Quality Management District (BAAQMD). Currently, regional, state, and federal agencies generally estimate methane emissions using bottom-up inventory methods that rely on a combination of activity data, emission factors, biogeochemical models and other information. Recent atmospheric top-down measurement estimates of methane emissions for the US as a whole (e.g., Miller et al., 2013) and in California (e.g., Jeong et al., 2013; Peischl et al., 2013) have shown inventories underestimate total methane emissions by ~ 50% in many areas of California, including the SF Bay Area (Fairley and Fischer, 2015). The goal of this research is to provide information to help improve methane emission estimates for the San Francisco Bay Area. The research effort builds upon our previous work that produced methane emission maps for each of the major source sectors as part of the California Greenhouse Gas Emissions Measurement (CALGEM) project (http://calgem.lbl.gov/prior_emission.html; Jeong et al., 2012; Jeong et al., 2013; Jeong et al., 2014). Working with BAAQMD, we evaluate the existing inventory in light of recently published literature and revise the CALGEM CH4 emission maps to provide better specificity for BAAQMD. We also suggest further research that will improve emission estimates. To accomplish the goals, we reviewed the current BAAQMD inventory, and compared its method with those from the state inventory from the California Air Resources Board (CARB), the CALGEM inventory, and recent published literature. We also updated activity data (e.g., livestock statistics) to reflect recent changes and to better represent spatial information. Then, we produced spatially explicit CH4 emission estimates on the 1-km modeling grid used by BAAQMD. We present the detailed activity data, methods and derived emission maps by sector. In total, we estimate the anthropogenic emissions for BAAQMD to be 116.4 Gg (1 Gg = 109 g) CH4/yr, with a likely uncertainty of ~ 50% or more (e.g., NRC, 2010; US-EPA, 2015). Including the emissions from wetland (Jeong et al., 2013), the total CH4 emission estimate for BAAQMD is 120.1 Gg CH4/yr. Table 1 summarizes the estimated CH4 emissions for 2011 by sector. The sectors were categorized following those that are used in recent regional emission quantification studies (e.g., Jeong et al., 2013; Peischl et al., 2013; Wecht et al., 2014). However, we note that this result is marginally lower than the top-down estimate of 240 ± 60 Gg CH4/yr (at 95% confidence) reported by Fairley and Fischer (2015), suggesting some combination of systematic error in the top-down estimate, underestimation of emissions from known sources, or as yet unidentified sources may be present. With respect to the relative contributions from different source sectors, the CH4 emissions from the region are dominated by urban activities. Landfill emissions represent 53% of the District’s total emission followed by livestock (16%) and natural gas (15%). These three dominant sectors account for 84% of the total anthropogenic emission in BAAQMD. This suggests that mitigation efforts need to focus on these three sources. Figure 1 shows the gridded anthropogenic CH4 emissions on the BAAQMD’s 1-km grid. In general, the spatial pattern of emissions follows the density of population while strong point sources are also distributed in the rural areas of the District. Detailed methods and emissions for each sector and county are described in the following sections.

54 ENVIRONMENTAL SCIENCES↗

Economics of Solar with Storage for Municipal Buildings in the City of Orlando [Slides]

This analysis was conducted under the first round of the Solar Energy Innovation Network (SEIN), which is a program led by the National Renewable Energy Laboratory (NREL). The program assembles diverse teams of stakeholders to research solutions to real-world challenges associated with solar energy adoption. This analysis supported the efforts of the Orlando: Renewable and Resilient SEIN team, which consists of the City of Orlando, Orlando Utility Commission, University of Central Florida's Solar Energy Center, and the Greenlink Group. The City of Orlando aims to deploy solar and solar-plus-storage to support energy resiliency, environmental quality, and the continued affordability of the electricity supply. The team's efforts include identifying the total solar potential available city-wide, conducting more detailed analysis for municipal facilities and distribution grid modeling, and gathering insights on municipal-level solar policies and building guidelines to share with other municipalities. This analysis supports the team's work by analyzing the techno-economic potential of solar photovoltaics (PV) and lithium-ion battery energy storage at 14 municipal sites in the City of Orlando.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

EGSim - a C++ Toolkit for Analysis of Power Grid Systems

We describe the load flow formulation and the solution algorithms available in the Electric power Grid Simulator (EGSim) software toolkit. EGSim contains tools aimed at simulating static load flow solutions for electric power grids. It parses power grid models described in IEEE Common Data Format, and generates solutions for the bus voltages and voltage angles, and real and reactive power values through the transmission lines. The software, written in C++, implements both Gauss-Seidel and Newton solution methods. Example results for the 118 bus models and 300 bus models are also presented.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE AGM and NAERM FY23 Annual Accomplishments

This document summarizes FY2023 work performed by Lawrence Livermore National Laboratory (LLNL) in support of the DOE Office of Electricity’s (OE) Advanced Grid Modeling (AGM) program, North American Energy Resilience Model (NAERM), and Protective Relaying (PR) program.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Small-Scale Variability of Solar Radiation (S2VSR) Field Campaign Report

Clouds introduce significant variability in solar radiation, in particular for shallow cumulus and other convective cloud types. Variability at small spatio-temporal scales (specifically below 1 min and 1 km) is, however, not properly resolved by current satellite observations nor atmospheric models. Additionally, 3D photon transport is commonly neglected in analyses and causes significant deviations in 1D radiative transfer simulations at such small scales. A specific challenge is to account for the point-like nature of traditional surface radiative flux measurements in comparison to spatially averaged quantities such as satellite pixel radiances or fluxes representing model grid boxes.

54 ENVIRONMENTAL SCIENCES↗

Graph-Learning-Assisted State and Event Tracking for Solar-Penetrated Power Grids with Heterogeneous Data Sources

Unlike transmission systems, distribution systems do not typically contain sufficient metering to enable real-time state estimation. The lack of sufficient real-time measurements prohibits accurate and timely monitoring of the state of distribution systems. As a result, control and optimal operation of distribution systems, especially those containing large numbers of renewable generation units are not possible without proper data and information about the current state of the system. The main motivation of this project is to address this shortcoming by developing an approach which provides “predicted” real-time measurements so that they can be used to execute a distribution system state estimator. Thus, the objective of the project is to make the distribution systems fully observable, such that the hosting capacity for solar generation can be accurately estimated, and unnecessary solar curtailments can be avoided. In order to accomplish this goal, the project investigated the use of a grid-model-informed machine learning (ML) tool which integrates heterogeneous data streams obtained from AMI meters, SCADA as well as PMU measurements and created synchronous measurement snapshots for the state estimator (SE); and developed a hybrid robust SE which provides not only accurate state estimates but also real-time feedback for the ML model refinement.

14 SOLAR ENERGY↗

Data Center Market Report

The data center market is poised to explode in the coming decade due to undeniable drivers such as continued adoption of generative AI, increased data storage needs, and enterprise integration of AI in numerous industries [1] [2] [3]. Scalable power and increased computational capacity are at the forefront of considerations for hyperscalers, the major cloud service providers in this space. Lawrence Livermore National Laboratory is uniquely poised to help with informed decision making for data center market leaders during this phase of explosive expansion. National grid modeling expertise and cutting edge innovations in computer cooling systems place LLNL in an enviable position for creating economic impact in the data center industry by leveraging its expertise in these areas which can help the data center market keep up with growing demand.

97 MATHEMATICS AND COMPUTING↗

Advanced Computing Annual Report 2025 [Slides]

In Fiscal Year (FY) 2025, the National Laboratory of the Rockies (NLR) continued to advance computing as a cornerstone of energy innovation, expanding the Kestrel high-performance computing (HPC) system to 56 peak petaflops. This growth strengthened Kestrel's role as a national asset for applied energy research, enabling larger, more complex simulations and accelerating the integration of artificial intelligence (AI) methods across the laboratory's computing portfolio. In FY 2025, AI was a component of most projects running on Kestrel, underscoring its central role in modern energy science and engineering. Kestrel supported a broad and diverse set of 507 modeling and simulation projects, engaging 855 researchers across the U.S. Department of Energy's (DOE's) Office of Critical Minerals and Energy Innovation (CMEI) portfolio and other offices, as well as partners from industry, academia, and utilities. These efforts span critical materials discovery, energy systems modeling, grid modernization, advanced manufacturing, and other areas essential to strengthening U.S. energy security and competitiveness. Together, these collaborations produced 708 technical outputs, including 293 peer-reviewed publications, reflecting both the depth and impact of the science enabled by NLR's computing capabilities. This year's report highlights the growing importance and benefit of AI throughout NLR's research programs and features work by early career researchers who are helping shape the future of computing-enabled energy innovation. Explore these sections and the many project successes captured in the pages that follow.

97 MATHEMATICS AND COMPUTING↗

5-minute Wind Power Data based on WFIP2 WRF Simulation

The second Wind Forecast Improvement Project (WFIP2) was a public-private partnership funded by the U.S. Department of Energy and NOAA, aimed at enhancing the forecast skill of numerical weather prediction models for turbine-height winds in regions with complex terrain. An 18-month Weather Research and Forecasting (WRF) model simulation was conducted over the Pacific Northwest, with model outputs validated against observational data collected during WFIP2. Simulated wind speeds were used to estimate wind power generation using reV (the Renewable Energy Potential Model developed by NREL) at ten wind project sites. Two sets of results were produced: one using wind speeds extracted from the model grid cell at the project centroid, and another using wind speeds from the actual turbine locations. For each dataset, power output was calculated using both actual turbine-specific power curves and nine generic power curves to convert wind speed into power.

17 WIND ENERGY↗

Editorial: Resolving atmospheric flow in complex environments: recent experiments in terrain and forest canopies

The characterization of atmospheric flows in complex environments, which may include steep terrain slopes and heterogeneous vegetation and/or forest cover, is a long-standing challenge in boundary-layer meteorology. Atmospheric observations are complicated by the presence of transient, terrain-induced flow features, forest-canopy-atmosphere interactions, and atmospheric stability effects, not to mention the logistical hurdles involved with instrument deployment, data analysis, and quality control. Furthermore, challenges in atmospheric modeling arise due to numerical errors associated with complex terrain flows, as well as reliance on simplified parameterizations for unresolved processes such as turbulent mixing and land-surface or forest-canopy-atmosphere interactions. These modeling challenges are exacerbated in the so-called “gray zone,” wherein features of interest have length scales that are similar to the model grid spacing, or when the principal flow layer is smaller than the grid spacing (e.g., slope flows).

54 ENVIRONMENTAL SCIENCES↗

SPISEA: A Python-based Simple Stellar Population Synthesis Code for Star Clusters

We present Stellar Population Interface for Stellar Evolution and Atmospheres (SPISEA), an open-source Python package that simulates simple stellar populations. The strength of SPISEA is its modular interface which offers the user control of 13 input properties including (but not limited to) the initial mass function, stellar multiplicity, extinction law, and the metallicity-dependent stellar evolution and atmosphere model grids used. The user also has control over the initial–final mass relation in order to produce compact stellar remnants (black holes, neutron stars, and white dwarfs). We demonstrate several outputs produced by the code, including color–magnitude diagrams, HR-diagrams, luminosity functions, and mass functions. SPISEA is object-oriented and extensible, and we welcome contributions from the community. The code and documentation are available on GitHub (https://github.com/astropy/SPISEA) and ReadtheDocs (https://spisea.readthedocs.io/en/latest/), respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

Kilonova Luminosity Function Constraints Based on Zwicky Transient Facility Searches for 13 Neutron Star Merger Triggers during O3

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run (O3). We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization area of 4480 deg 2 , a median distance of 267 Mpc, and false-alarm rates ranging from 1.5 to 10 -25 yr -1 . The ZTF coverage in the g and r bands had a median enclosed probability of 39%, median depth of 20.8 mag, and median time lag between merger and the start of observations of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UltraViolet/Optical/InfraRed (UVOIR) photometric points, 64 OIR spectra, and three radio images using 17 different telescopes. We find no promising kilonovae (radioactivity-powered counterparts), and we show how to convert the upper limits to constrain the underlying kilonova luminosity function. Initially, we assume that all GW triggers are bona fide astrophysical events regardless of false-alarm rate and that kilonovae accompanying BNS and NSBH mergers are drawn from a common population; later, we relax these assumptions. Assuming that all kilonovae are at least as luminous as the discovery magnitude of GW170817 (-16.1 mag), we calculate that our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than-16.6 mag (the extrapolated peak magnitude of GW170817) and fade at a rate of 1 mag day -1 (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations based on the online classifications, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6 mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, no more than <57% (<89%) of putative kilonovae could be brighter than -16.6 mag assuming flat evolution (fading by 1 mag day -1 ) at the 90% confidence level. If we further take into account the online terrestrial probability for each GW trigger, we find that no more than <68% of putative kilonovae could be brighter than -16.6 mag. Comparing to model grids, we find that some kilonovae must have M ej < 0.03 M ⊙ , X lan > 10 -4 , or φ > 30° to be consistent with our limits. We look forward to searches in the fourth GW observing run; even 17 neutron star mergers with only 50% coverage to a depth of-16 mag would constrain the maximum fraction of bright kilonovae to <25%.

79 ASTRONOMY AND ASTROPHYSICS↗

The BTSbot-nearby Discovery of SN 2024jlf: Rapid, Autonomous Follow-up Probes Interaction in an 18.5 Mpc Type IIP Supernova

We present observations of the Type IIP supernova (SN) SN 2024jlf, including spectroscopy beginning just 0.7 days (∼17 hr) after first light. Rapid follow-up was enabled by the new BTSbot-nearby program, which involves autonomously triggering target-of-opportunity requests for new transients in Zwicky Transient Facility data that are coincident with nearby (D < 60 Mpc) galaxies and identified by the BTSbot machine learning model. Early photometry and nondetections shortly prior to first light show that SN 2024jlf initially brightened by >4 mag day −1 , quicker than ∼90% of Type II SNe. Early spectra reveal weak flash ionization features: narrow, short-lived (1.3 < τ[days] < 1.8) emission lines of Hα, He II , and C IV . Assuming a wind velocity of v w = 50 km s −1 , these properties indicate that the red supergiant progenitor exhibited enhanced mass loss in the last year before explosion. We constrain the mass-loss rate to $1{0}^{-4}\lt \dot{M}\,[{M}_{\odot }\,{\mathrm{yr}}^{-1}]\lt 1{0}^{-3}$ by matching observations to model grids from two independent radiative hydrodynamics codes. BTSbot-nearby automation minimizes spectroscopic follow-up latency, enabling the observation of ephemeral early-time phenomena exhibited by transients.

core-collapse supernovae↗

Molecular and physical characteristics of aerosol at a remote free troposphere site: implications for atmospheric aging

Aerosol properties are transformed by atmospheric processes during long-range transport and play a key role in the Earth's radiative balance. To understand the molecular and physical characteristics of free tropospheric aerosol, we studied samples collected at the Pico Mountain Observatory in the North Atlantic. The observatory is located in the marine free troposphere at 2225 m above sea level, on Pico Island in the Azores archipelago. The site is ideal for the study of long-range-transported free tropospheric aerosol with minimal local influence. Three aerosol samples with elevated organic carbon concentrations were selected for detailed analysis. FLEXPART retroplumes indicated that two of the samples were influenced by North American wildfire emissions transported in the free troposphere and one by North American outflow mainly transported within the marine boundary layer.Ultrahigh-resolution Fourier transform ion cyclotron resonance mass spectrometry was used to determine the detailed molecular composition of the samples. Thousands of molecular formulas were assigned to each of the individual samples. On average ~60 % of the molecular formulas contained only carbon, hydrogen, and oxygen atoms (CHO), ~30 %contained nitrogen (CHNO), and ~10 % contained sulfur (CHOS). The molecular formula compositions of the two wildfire-influenced aerosol samples transported mainly in the free troposphere had relatively low averageO/C ratios (0.48±0.13 and 0.45±0.11) despite the 7–10 days of transport time according to FLEXPART. In contrast, the molecular composition of the North American outflow transported mainly in the boundary layer had a higher average O/C ratio (0.57±0.17) with 3 days of transport time. To better understand the difference between free tropospheric transport and boundary layer transport, the meteorological conditions along the FLEXPART simulated transport pathways were extracted from the Global Forecast System analysis for the model grids. We used the extracted meteorological conditions and the observed molecular chemistry to predict the relative-humidity-dependent glass transition temperatures (T g ) of the aerosol components. Comparisons of the T g to the ambient temperature indicated that a majority of the organic aerosol components transported in the free troposphere were more viscous and therefore less susceptible to oxidation than the organic aerosol components transported inthe boundary layer. Although the number of observations is limited, the results suggest that biomass burning organic aerosol injected into the free troposphere is more persistent than organic aerosol in the boundary layer having broader implications for aerosol aging.

54 ENVIRONMENTAL SCIENCES↗

Global anthropogenic emissions (CAMS-GLOB-ANT) for the Copernicus Atmosphere Monitoring Service simulations of air quality forecasts and reanalyses

Anthropogenic emissions are the result of many different economic sectors, including transportation, power generation, industrial, residential and commercial activities, waste treatment and agricultural practices. Air quality models are used to forecast the atmospheric composition, analyze observations and reconstruct the chemical composition of the atmosphere during the previous decades. In order to drive these models, gridded emissions of all compounds need to be provided. This paper describes a new global inventory of emissions called CAMS-GLOB-ANT, developed as part of the Copernicus Atmosphere Monitoring Service (CAMS; https://doi.org/10.24380/eets-qd81, Soulie et al., 2023). The inventory provides monthly averages of the global emissions of 36 compounds, including the main air pollutants and greenhouse gases, at a spatial resolution of 0.1° × 0.1° in latitude and longitude, for 17 emission sectors. The methodology to generate the emissions for the 2000–2023 period is explained, and the datasets are analyzed and compared with publicly available global and regional inventories for selected world regions. Depending on the species and regions, good agreements as well as significant differences are highlighted, which can be further explained through an analysis of different sectors as shown in the figures in the Supplement.

54 ENVIRONMENTAL SCIENCES↗

DCMIP2016: a review of non-hydrostatic dynamical core design and intercomparison of participating models

Abstract. Atmospheric dynamical cores are a fundamental component of global atmospheric modeling systems and are responsible for capturing the dynamical behavior of the Earth's atmosphere via numerical integration of the Navier–Stokes equations. These systems have existed in one form or another for over half of a century, with the earliest discretizations having now evolved into a complex ecosystem of algorithms and computational strategies. In essence, no two dynamical cores are alike, and their individual successes suggest that no perfect model exists. To better understand modern dynamical cores, this paper aims to provide a comprehensive review of 11 non-hydrostatic dynamical cores, drawn from modeling centers and groups that participated in the 2016 Dynamical Core Model Intercomparison Project (DCMIP) workshop and summer school. This review includes a choice of model grid, variable placement, vertical coordinate, prognostic equations, temporal discretization, and the diffusion, stabilization, filters, and fixers employed by each system.

54 ENVIRONMENTAL SCIENCES↗

Medium Voltage Testbed for Comparing Advanced Power Line Sensors vs. Transformers with Electrical Grid Events

Electrical utilities have relied upon potential transformers (PTs) and current transformers (CTs) for very accurate metering and to provide reliable signals for protective relays. These devices measure phase voltages and currents and are commissioned by electrical engineers. PTs/CTs can detect and react to various electrical anomalies that could adversely affect electrical grid operations. Less expensive alternative sensing technologies offer the possibility of wider deployment, particularly in grids that employ distributed energy resources. In this work, the performance of an advanced medium-voltage sensor is compared with that of a reference PT and a CT and experimentally evaluated for different power grid scenarios on an advanced outdoor power line sensor testbed at the U.S. Department of Energy’s Oak Ridge National Laboratory. The sensor is based on a capacitive divider for voltage monitoring and a Rogowski coil with an integrator for current monitoring. This study simulated a power grid model based on a utility circuit at the Riverside EPB of Chattanooga. The simulation circuit was created with MATLAB/Simulink software and was integrated into an RT-LAB project to run with the OP4510 real-time simulator at the OPLST. During the tests, the real-time simulations were run for 40 s, and the signal to record the test event with the power meter was set at 30 s for the event-trigger circuit. The advanced outdoor power line sensor testbed has a real-time simulator that is used to generate transient scenarios (e.g., electrical faults, capacitor bank operation, and service restoration), while the analog signals are recorded by the same high-resolution power meter. The behaviors of analog signals, harmonic components, total harmonic distortion, and crest factors are assessed for this power line sensor and compared with those of the reference PT/CT because of the absence of testing standards for advanced outdoor power line sensors. The results showed that this OPLS technology responded identically to the PT and CT under all conditions.

Piesciorovsky, Emilio [ORNL]↗

Low-level circulations induced by an upper-level jet in a fine mesh numerical model

Attention is given to varying the initial conditions in the Penn State three-dimensional numerical model in the development of the low-level jet over the Ohio Valley during the 10 May 1973 severe weather event. The results of two forecasts are summarized: (1) the control forecast which interpolated the National Meteorological Center's operational analysis to the mesoscale model grid; and (2) modifying the control forecast by inserting artificial soundings and surface observations. Enhancements are made on the basis of subsynoptic scale analyses by Uccellini and surface data. The results indicate that the low-level jet was in part induced by pressure changes associated with deceleration and horizontal divergence of air in the exit region of the upper-level jet. The enhancement of the initial moisture field using surface reports of clouds and precipitation has a substantial effect on precipitation.

Burkhart, R. P.↗