Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Structured Grids”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The Future of Vehicle Grid Integration: Harnessing the Flexibility of EV Charging

This document lays out a shared vision for a beneficial, EV-integrated future where EVs are safely and securely connected, reliably served, and harmonized with the electric grid. It was developed as part of the U.S. Department of Energy’s (DOE) EVGrid Assist initiative. Stakeholder input was gathered through individual and collective conversations, including eight listening sessions with more than 100 participants representing utilities of different sizes and operating structures, manufacturers of vehicles and chargers, national associations, standards organizations, Tribes, fleet managers, consumer advocates, charging network operators, community-based organizations, utility regulators, consultants, vendors, labor, and environmental justice organizations. This document focuses on the future of electric on-road U.S. transportation, specifically the integration of light-duty vehicles (LDV) and medium- and heavy-duty vehicles (MHDV) and their charging infrastructure with the electric grid. However, many of the insights here may apply to electrifying other transportation modes, as well as other distributed energy resources (DER).

EVGrid Assist initiative, EV-integrated future, el↗
DOE OSTI2020

REopt Lite Overview & Training Exercise

This training exercise provides users with an introduction to and hands-on, interactive exploration of REopt Lite's capabilities. REopt Lite is a free, publicly available techno-economic optimization web tool for distributed energy systems, developed at the National Renewable Energy Laboratory (NREL). REopt Lite helps organizations evaluate the economic viability of grid-connected solar photovoltaics (PV), wind turbines, and battery storage; identify system sizes and battery dispatch strategies to minimize energy costs; and estimate how long a system can sustain critical load during a grid outage. The model is formulated as a mixed-integer linear program based in an underlying application programming interface (API) that is also free and publicly available. This training activity is structured as a group exercise. Participants split into eight groups and each group is assigned a different hypothetical site to model and assess the opportunity for solar PV + battery storage. Groups work together to develop results for their site and then re-convene to compare and discuss results, inputs/drivers of the analysis, and other factors impacting the decision-making process for behind-the-meter solar PV and battery storage.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

The June 2012 North American Derecho: A Testbed for Evaluating Regional and Global Climate Modeling Systems at Cloud‐Resolving Scales

Abstract In this paper, we introduce a testbed for evaluating and comparing climate modeling systems at cloud resolving scales using hindcasts of the June 2012 North American derecho. To demonstrate its utility for model intercomparison, the testbed is applied to two models: the regionally‐refined Simple Cloud‐Resolving E3SM Atmosphere Model (SCREAM) at 6.5, 3.25 and 1.625 km grid spacing and the Weather Research and Forecasting (WRF) model with 3.2 and 1.6 km grid spacing. We find the simulation results to be highly sensitive to the initial conditions (ICs), initialization time, and model configurations, with ICs from the Rapid Refresh producing the best simulation. Significant improvement is identified in both models as horizontal grid spacing is refined. While a propagation delay of approximately 2 hr is found in both models, SCREAM at 1.625 km simulates the observed bow echo structure of the derecho well and predicts strong surface gusts that exceed 30 m/s. In comparison, WRF has difficulty producing surface wind over 25 m/s, with wind gusts in WRF 42%–46% lower than in SCREAM. However, WRF has a lower bias in simulating cloud top temperature and extent, but overestimates precipitation intensity. Both models reproduce the observed outgoing longwave radiation spatial patterns well (Pearson correlation >0.88), but, compared with NEXRAD observations, simulate generally larger areas of composite radar reflectivity >40 dBZ and underestimate the precipitating area by ∼47%.

2012 North American derecho↗
DOE OSTI2025

Advancing Dynamic Modeling of Grid-Connected PV Inverter Using Bi-LSTM-Based AI Model

Power electronic converters (PECs) are widely used in modern power systems to facilitate the interconnection between various AC or DC sources and loads. Because of the extensive integration, the power system has grown into a more dynamic system in which the dynamics of the PECs must be adequately modeled. The paper presents a new bidirectional long short-term memory (Bi-LSTM) method for evaluating grid-connected inverter-based resources (IBR) dynamics. The method is tested using real hardware data from a grid -connected commercial inverter in laboratory experiments. Results show the Bi-LSTM model accurately reproduces the detailed IBR model's dynamics, even when the internal structure is unknown and parameters are unknown, preventing the disclosure of the manufacturer's confidential data.

Subedi, Sunil [ORNL] (ORCID:000000034069090X)↗
DOE OSTI2024

Impact of Reordering on the LU Factorization Performance of Bordered Block-Diagonal Sparse Matrix

Power engineers rely on computer-based simulation tools to assess grid performance and ensure security. At the core of these tools are solvers for sparse linear equations. When transformed into a bordered block-diagonal (BBD) structure, part of the sparse linear equation solving can be parallelized. This work focuses on using the Schur-complement-based method for LU factorization on BBD matrices, specifically, Jacobian matrices from large-scale systems. Our findings show that the natural ordering method outperforms the default ordering method in computational performance for each block of the BBD matrix. This observation is validated using synthetic 25k-bus and 70k-bus cases, showing a speedup of up to 38% when using natural ordering without permutation. Additionally, the impact of the number of partitions is studied, and the result shows that computational performance improves with more, smaller partitions in the BBD matrices.

BBD matrix↗
DOE OSTI2023

Vacant niches help predict invasion risk by birds

Abstract Aim If communities have a ceiling determining the number of species that can coexist, then the ability of alien species to establish at any given location should be related to the quantity of vacant niches available. We developed a new approach to estimate the extent to which niches are vacant and then explored the relationship between vacant niches and alien species. Location Global with focal tests in Europe and North America. Taxon Birds. Methods Drawing on a global classification of trophic structures for birds, we calculated a ‘vacant niche ratio’ metric to quantify the expected level of saturation (i.e. number of vacant niches) for each 1° × 1° grid cell globally, based on the difference between expected and observed numbers of bird species for all trophic guilds. Next, we used random forests to examine if the presence of plant‐invertivore (whose food source represents plants, seeds, fruits and invertebrates) and granivore alien bird species was associated with the vacant niche ratio across well‐sampled regions. Results In Europe, we found a significant relationship between alien species and vacant niches, with greater numbers of alien species being found in communities that offered greater numbers of vacant niches overall, and across habitat types. In North America, we found no significant relationship between plant‐invertivores and vacant niches, while for granivores, we found that areas with fewer vacant niches had greater numbers of alien species, especially in forests. However, vacant niches alone correctly predict 69% of the presence of alien bird species when combining both regions and trophic guilds. Main Conclusions Most regions of the world have unsaturated bird communities, with the level of saturation in communities varying within regions and trophic guilds. We found that although often‐neglected, vacant niches are likely to be, at least partially, related with the successful establishment of alien bird species.

González‐del‐Pliego, Pamela↗
DOE OSTI2024

In situ characterization of metastable Pb3O5 and Pb2O3 phases during thermal decomposition of PbO2 to PbO

Nonstoichiometric lead oxides play a key role in the formation and cycling of the positive electrodes in a lead acid battery. These phases have been linked to the underutilization of the positive active material but also play a key role in the battery’s cycle life, providing inter-particle adhesion and the connection to the underlying lead grid. Similar phases have previously been identified by mass loss or color change during thermal annealing of PbO2 to PbO, suggesting that at least two intermediate PbOx phases exist. Using multiple in situ analysis techniques (powder diffraction, x-ray absorption, x-ray photoelectron spectroscopy) and ex situ nuclear magnetic resonance measurements, the structural conversion and changes in the lead oxidation states were identified during this process. Isolation of the PbOx phases enabled confirmation of Pb3O5 and Pb2O3 by diffraction and the first 207Pb NMR measurement of these intermediates.

Kinnibrugh, Tiffany L.↗
DOE OSTI2023

Control Design of Passive Grid-Forming Inverters in Port-Hamiltonian Framework

This article presents a modified dispatchable virtual oscillator control approach for achieving the passivity of gridforming inverters (GFMs), without assuming constant voltage and constant frequency. The proposed control framework utilizes the Port-Hamiltonian (PH) based structure that mimics the behaviors of coupled harmonic oscillators, along with an energy ‘pumpingor-damping’ block and the Control by Interconnection (CbI) technique, to render the inverter passive. Once passivity is achieved, transient stability of the system will be guaranteed. The proposed control framework is composed of three loops: an outer power dispatching loop that generates the voltage and frequency references, a virtual oscillator loop that emulates the spontaneous synchronization of oscillators, and an inductor current loop that maintains lossless interconnection in PH systems. In conclusion, the study shows that the proposed control approach ensures the passivity of GFMs, facilitating the transient stability design of multi-inverter systems, as interconnections of passive systems remain passive and stable.

42 ENGINEERING↗
DOE OSTI2020

A High-Energy NASICON-Type Cathode Material for Na-Ion Batteries

Over the last decade, Na-ion batteries have been extensively studied as low-cost alternatives to Li-ion batteries for large-scale grid storage applications; however, the development of high-energy positive electrodes remains a major challenge. Materials with a polyanionic framework, such as Na superionic conductor (NASICON)-structured cathodes with formula Na x M 2 (PO 4 ) 3 , have attracted considerable attention because of their stable 3D crystal structure and high operating potential. Herein, a novel NASICON-type compound, Na 4 MnCr(PO 4 ) 3 , is reported as a promising cathode material for Na-ion batteries that deliver a high specific capacity of 130 mAh g –1 during discharge utilizing high-voltage Mn 2+/3+ (3.5 V), Mn 3+/4+ (4.0 V), and Cr 3+/4+ (4.35 V) transition metal redox. In addition, Na 4 MnCr(PO 4 ) 3 exhibits a high rate capability (97 mAh g –1 at 5 C) and excellent all-temperature performance. In situ X-ray diffraction and synchrotron X-ray diffraction analyses reveal reversible structural evolution for both charge and discharge.

25 ENERGY STORAGE↗
DOE OSTI2025

On the numerical sensitivity of cellular automata grain structure predictions to large thermal gradients and cooling rates

Cellular automata (CA) models of as-solidified grain structure, originally developed and applied to casting, have become a common means of predicting grain structure resulting from Additive Manufacturing (AM) processes. The majority of these models are based on the decentered octahedron approach, which attempts to correct for the effect of grid anisotropy on the prediction of competitive solidification of dendritic grains. However, AM solidification occurs under cooling rates ($\dot{T}$) and thermal gradients (G) that are orders of magnitude larger than those encountered in casting, and no systematic investigation on the effect of the CA model cell size (Δx) and time step (Δt) on AM microstructure predictions has been performed. Here, in this study, such an investigation is first performed via simulation of individual grains of various crystallographic orientations with a fixed, unidirectional G, showing that CA prediction of the steady-state undercooling matched the expected values based on the interfacial response function at small G and deviated from the expected values at large G. Simulation of competitive growth of multiple grains showed a weakening of the predicted texture as G and Δx became large. Simulation of solidification under AM conditions, where G and $\dot{T}$ vary spatially across the melt pools, showed that not only does grain selection weaken and deviate from expectations at large Δx, but grains with crystallographic $\langle$100$\rangle$ aligned with the grid directions are more adversely affected by the temperature field discontinuities than grains with other crystallographic orientations. Despite the fact that the exact grain competition results depended on Δt, the overall texture development was notably less sensitive to Δt than Δx, provided that a reasonable value of Δt is selected based on the ratio of Δx to the maximum local solidification velocity in the simulation domain. Finally, from the directional solidification and AM simulation results, an analysis of computational cost compared to simulation resolution is performed based on an equation derived to quantify the relatively inaccuracy in grain selection based on the model and temperature field inputs. From this analysis, it is concluded that there is a need for algorithmic improvements to improve CA grain competition accuracy for large G processing conditions as sufficiently small Δx to resolve the necessary competition is intractable for many AM processing conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations

There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Scalable bottom-up synthesis of Co-Ni–doped graphene

Introducing heteroatoms into graphene is a powerful strategy to modulate its catalytic, electronic, and magnetic properties. At variance with the cases of nitrogen (N)– and boron (B)–doped graphene, a scalable method for incorporating transition metal atoms in the carbon (C) mesh is currently lacking, limiting the applicative interest of model system studies. This work presents a during-growth synthesis enabling the incorporation of cobalt (Co) alongside nickel (Ni) atoms in graphene on a Ni(111) substrate. Single atoms are covalently stabilized within graphene double vacancies, with a Co load ranging from 0.07 to 0.22% relative to C atoms, controllable by synthesis parameters. Structural characterization involves variable-temperature scanning tunneling microscopy and ab initio calculations. The Co- and Ni-codoped layer is transferred onto a transmission electron microscopy grid, confirming stability through scanning transmission electron microscopy and electron energy loss spectroscopy. This method holds promise for applications in spintronics, gas sensing, electrochemistry and catalysis, and potential extension to graphene incorporation of similar metals.

Science & Technology - Other Topics↗
DOE OSTI2021

Optimal toolpath generation system and method for additively manufactured composite materials

The present disclosure relates to a system for optimizing an additive manufacturing (AM) process. The system may use a printing component for using a material to form a component in a layer-by-layer process. An electronic controller may control movement of one of the printing component or a substrate on which the component is being formed, in a manner to optimize a toolpath for the printing component as each layer of the component is formed. Optimization operations are performed using an optimization subsystem which enables manufacturability constraints, as well as optimized toolpaths for each layer of the component, to be defined using contours of a level-set function. The level-set function may be used to define the optimized toolpaths within a fixed, predetermined grid, and the optimized toolpaths then used to generate suitable code for controlling movement of the printing component relative to the substrate to create the part or structure in a layer-by-layer process, using the optimized toolpaths.

Lewicki, James↗
DOE OSTI2025

Resilient Distributed Frequency Regulation of Renewable Generators under Communication Interruptions

Modern power systems (MPSs) face significant challenges due to the high penetration of renewable energy sources (RESs) and new types of loads such as electric vehicles (EVs). Traditional load frequency control (LFC) methods struggle with the intermittent, stochastic nature of RESs, the near-zero inertia of power-electronics-based generators, and the mobility of controllable loads and battery systems. This paper introduces a novel resilient distributed frequency regulation method to address these issues. The proposed method employs a state space model to represent the dynamic behavior of participating power sources while accounting for stochastic switching processes to model structural and parameter variations caused by disruptions such as generator connection/disconnection, communication interruptions, and physical faults. By integrating these dynamic and stochastic components, the method treats power grids as a comprehensive stochastic hybrid system. Our method enhances conventional frequency control by incorporating local stability control, neighborhood control decoupling, and coordination feedback. Theoretical analyses establish the stability, convergence, and resilience of the proposed method, and its effectiveness is validated through case studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Prospects for Astrobiology and Technosignature Searches with the Vera C. Rubin Observatory Legacy Survey of Space and Time

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will map sources in multiband colour--variability space. We present a prototype coherence-based framework for astrobiology and technosignature searches, in which candidates are treated as structured departures from natural astrophysical manifolds rather than isolated photometric outliers. We illustrate the framework with three simulated cases: five Kuiper Belt Object (KBO) surface/activity states, a grid of 649 synthetic exoplanet spectra with vegetation-red-edge-like (VRE) perturbations, and 500 synthetic multiband light curves, each projected into LSST-like observable space and analysed through colour geometry, chromatic variability, and cross-band coherence. Key results include a full-colour Mahalanobis distance $D\approx5.1$ for the weak-coma KBO state (${\sim}5σ$ in the five-dimensional colour vector), an indicative VRE coherence threshold at $f_{\rm crit}\approx0.13$, and an idealised stacking forecast reaching $5σ$ under optimistic assumptions. We show, using a small Gaia~DR3 stellar sample, that stellar colour and photometric stability may inform the prioritisation of Galactic regions for applying such coherence diagnostics.

Kovačević, Andjelka B. [Belgrade U.] (ORCID:000000↗
DOE OSTI2026

Automated ICRF heating surrogate modeling via machine learning

This work introduces automated machine learning workflows that address critical bottlenecks in surrogate model development for Ion Cyclotron Range of Frequencies (ICRF) heating applications. The automated framework includes data analysis tools that transform raw datasets into actionable insights in seconds, replacing weeks of manual exploratory effort and ensuring consistent, reproducible dataset characterization. By integrating advanced hyperparameter optimization (HPO) methods including Bayesian optimization via BoTorch and Tree-structured Parzen Estimators (TPE), the framework significantly reduces model development time from weeks to hours, decreasing computational cost and required expertise, while enabling high-accuracy surrogate models. Compared to traditional hyperparameter scanning (HPS) techniques such as methodical, randomized, and grid searches, HPO methods achieve superior convergence and predictive performance, even when compared to already well-tuned reference models. On NSTX High Harmonic Fast Wave (HHFW) heating datasets, both Random Forest Regressor (RFR) and neural network surrogates demonstrate improved accuracy, achieving R 2 values beyond 0.97 and 0.98, respectively. The results show that while HPO gains are modest for robust architectures like RFR, they become essential for more sensitive models such as neural networks, highlighting the trade-offs across optimization strategies. Through automated workflows that eliminate manual hyperparameter tuning and require minimal ML expertise, this work enables widespread adoption of high-fidelity surrogate models across the fusion community for real-time plasma control, uncertainty quantification, rapid experimental scenario development, and integrated system optimization.

Sanchez-Villar, Alvaro [Princeton Plasma Physics L↗
DOE OSTI2022

Volumetric Rendering on Wavelet-Based Adaptive Grid

Numerical modeling of physical phenomena frequently involves processes across a wide range of spatial and temporal scales. In the last two decades, the advancements in wavelet-based numerical methodologies to solve partial differential equations, combined with the unique properties of wavelet analysis to resolve localized structures of the solution on dynamically adaptive computational meshes, make it feasible to perform large-scale numerical simulations of a variety of physical systems on a dynamically adaptive computational mesh that changes both in space and time. Volumetric visualization of the solution is an essential part of scientific computing, yet the existing volumetric visualization techniques do not take full advantage of multi-resolution wavelet analysis and are not fully tailored for visualization of a compressed solution on the wavelet-based adaptive computational mesh. Our objective is to explore the alternatives for the visualization of time-dependent data on space-time varying adaptive mesh using volume rendering while capitalizing on the available sparse data representation. Two alternative formulations are explored. The first one is based on volumetric ray casting of multi-scale datasets in wavelet space. Rather than working with the wavelets at the finest possible resolution, a partial inverse wavelet transform is performed as a preprocessing step to obtain scaling functions on a uniform grid at a user-prescribed resolution. As a result, a solution in physical space is represented by a superposition of scaling functions on a coarse regular grid and wavelets on an adaptive mesh. An efficient and accurate ray casting algorithm is based just on these coarse scaling functions. Additional details are added during the ray tracing by taking an appropriate number of wavelets into account based on support overlap with the interpolation point, wavelet coefficient magnitude, and other characteristics, such as opacity accumulation (front to back ordering) and deviation from frontal viewing direction. The second approach is based on complementing of wavelet-based adaptive mesh to the traditional Adaptive Mesh Refinement (AMR) mesh. Both algorithms are illustrated and compared to the existing volume visualization software for Rayleigh-Benard thermal convection and electron density data sets in terms of rendering time and visual quality for different data compression of both wavelet-based and AMR adaptive meshes.

Vezolainen, Alexei V.↗