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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

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↗
DOE OSTI2025

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↗
DOE OSTI2024

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↗
DOE OSTI2020

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↗
DOE OSTI2020

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↗
DOE OSTI2025

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↗
DOE OSTI2018

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↗
DOE OSTI2024

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↗
DOE OSTI2017

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↗
DOE OSTI2025

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]↗
NASA NTRS1981

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.↗
NASA NTRS1984

Vectorization of a penalty function algorithm for well scheduling

In petroleum engineering, the oil production profiles of a reservoir can be simulated by using a finite gridded model. This profile is affected by the number and choice of wells which in turn is a result of various production limits and constraints including, for example, the economic minimum well spacing, the number of drilling rigs available and the time required to drill and complete a well. After a well is available it may be shut in because of excessive water or gas productions. In order to optimize the field performance a penalty function algorithm was developed for scheduling wells. For an example with some 343 wells and 15 different constraints, the scheduling routine vectorized for the CYBER 205 averaged 560 times faster performance than the scalar version.

Absar, I.↗
NASA NTRS1984

Analysis of airfoil transitional separation bubbles

A previously developed local inviscid-viscous interaction technique for the analysis of airfoil transitional separation bubbles, ALESEP (Airfoil Leading Edge Separation) has been modified to utilize a more accurate windward finite difference procedure in the reversed flow region, and a natural transition/turbulence model has been incorporated for the prediction of transition within the separation bubble. Numerous calculations and experimental comparisons are presented to demonstrate the effects of the windward differencing scheme and the natural transition/turbulence model. Grid sensitivity and convergence capabilities of this inviscid-viscous interaction technique are briefly addressed. A major contribution of this report is that with the use of windward differencing, a second, counter-rotating eddy has been found to exist in the wall layer of the primary separation bubble.

Davis, R. L.↗
NASA NTRS1984

Analysis of airfoil transitional separation bubbles

A previously developed local inviscid-viscous interaction technique for the analysis of airfoil transitional separation bubbles, ALESEP (Airfoil Leading Edge Separation), has been modified to utilize a more accurate windward finite difference procedure in the reversed flow region, and a natural transition turbulence model has been incorporated for the prediction of transition within the separation bubble. Numerous calculations and experimental comparisons are presented to demonstrate the effects of the windward differencing scheme and the natural transition turbulence model. Grid sensitivity and convergence capabilities of this inviscid-viscous interaction technique are briefly addressed. A major conclusion of this paper is that a second, counter-rotating eddy has been found to exist in the wall layer of the primary separation bubble with the use of windward differencing.

Davis, R. L.↗
NASA NTRS1985

Neutral parcel transport in the high latitude F-region

The NCAR Thermospheric General Circulation model (TGCM) is used together with a 'post-processor' diagnostic package to analyze neutral parcel transport in the high latitude F-region. The diagnostic package includes a feature that enables individual parcels of neutral gas to be tracked by interpolation through the TGCM model grid in time and space. Various trajectories have been calculated to illustrate the predicted loci of neutral gas parcels that encounter large ion-drag forces due to the ionospheric mapping of the magnetospheric convection electric field. The trajectories, calculated for a 'steady state', diurnally reproducible run of the TGCM, show that the cusp region plays an exceedingly important regulatory or nozzle role for the overall flow pattern and has great influence on the energization of the F-region polar cap. In particular, it is found that the trajectory of a 'typical' F-region neutral gas parcel passes through the cusp region at least once and sometimes twice in a 24 hour period due to the convergent geometry of the ion convection pattern.

Killeen, T. L.↗
NASA NTRS1985

The error associated with use of various forms of the divergence equation to diagnose geopotential and temperature

An experimental study was made of errors arising from diagnosing the geopotential and temperature for numerical weather prediction models by means of an inviscid divergence equation (DE). The experiment consisted of a simulation of a snow storm which occurred in 1979, with observational data being input for every point of the model grid and for the geopotential boundary condition, thereby eliminating errors in the horizontal and vertical analysis. The DE terms were recorded at each grid point, and were found to include synoptic-scale geopotential height errors up to 23 m and rms temperature errors exceeding 26 C in the PBL. The errors were minimal above 800 mb. The various unsuccessful attempts to reduce or eliminate the errors are described.

Modica, G. M.↗
NASA NTRS1989

Renormalization group analysis of turbulence

The objective is to understand and extend a recent theory of turbulence based on dynamic renormalization group (RNG) techniques. The application of RNG methods to hydrodynamic turbulence was explored most extensively by Yakhot and Orszag (1986). An eddy viscosity was calculated which was consistent with the Kolmogorov inertial range by systematic elimination of the small scales in the flow. Further, assumed smallness of the nonlinear terms in the redefined equations for the large scales results in predictions for important flow constants such as the Kolmogorov constant. It is emphasized that no adjustable parameters are needed. The parameterization of the small scales in a self-consistent manner has important implications for sub-grid modeling.

Smith, Leslie M.↗
NASA NTRS1991

Effects of engine emissions from high-speed civil transport aircraft: A two-dimensional modeling study, part 2

The AER two-dimensional chemistry-transport model is used to study the effect of supersonic and subsonic aircraft operation in the 2010 atmosphere on stratospheric ozone (O3). The results show that: (1) the calculated O3 response is smaller in the 2010 atmosphere compared to previous calculations performed in the 1980 atmosphere; (2) with the emissions provided, the calculated decrease in O3 column is less than 1 percent; and (3) the effect of model grid resolution on O3 response is small provided that the physics is not modified.

Ko, Malcolm K. W.↗