SEARCH · Engineering Papers

Results for “grid codes”

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.

SourceResearch connectionFrom yearThrough year

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

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2026

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part I: Lagrangian and indirect Euler AMR algorithms

Many applications of physics and engineering involve wide ranges of time and spatial scales. The numerical simulation of localized small scales such as shock waves and material interfaces requires a large number of computational cells in these regions. For these applications, Lagrangian and Arbitrary-Lagrangian-Eulerian (ALE) related methods are engaging since the moving mesh feature naturally brings mesh cells on shock discontinuities and material interfaces are carefully captured. In addition, Adaptive-Mesh-Refinement (AMR) strategies aim to optimize computational resources by concentrating finer mesh cells only in areas of interest while using coarser cells elsewhere. A key but challenging AMR requirement consists in efficiently distributing the computational effort to achieve high accuracy without the prohibitive computational costs associated with uniformly fine grids. Here, in this document, the coupling of the p4est AMR library with a cell-centered Lagrangian scheme is presented with the goal to perform reliable 3D Lagrangian-AMR and indirect Euler-AMR multi-material simulations. In particular, it is shown that starting from a 3D indirect ALE code, the memory management and load balancing requirements can be delegated to an external library (here the p4est library) to unlock ALE-AMR capabilities. First, we present a strategy to transcribe the octant-based connectivity of the 3D AMR framework with that of an unstructured mesh of polygonal cells used in Lagrangian hydrodynamics. Then, we show how refinement and coarsening operations must be adapted to the particular Lagrangian framework to ensure the conservation of volume during those steps. Finally, several numerical test cases are presented that demonstrate the capabilities of the Lagrangian-AMR and indirect Euler-AMR algorithms.

3D cell-centered Lagrangian numerical scheme↗

DOE OSTI2025

Decoding Resilience by Modeling Outage and Restoration Processes in Distribution Grids: A Pittsburgh Case Study

Climate-induced extreme weather events, such as floods and heatwaves, pose significant challenges to the resilience of urban power distribution grids. This paper examines the outage and restoration dynamics of Pittsburgh's power grid during the flood event in April 2024 and three heatwaves in June, July, and August 2024. We introduce a comprehensive modeling framework that integrates outage and restoration processes, enabling the quantification of resilience metrics, including total customer-hours of power outage, maximum residual values, and restoration durations across various ZIP codes. Our analysis highlights spatial disparities in outage impacts and restoration efficiencies, with ZIP Code 15222 experiencing the highest cumulative disruptions, while ZIP Code 15217 shows lower susceptibility. By linking statistical trends with weather events, this study underscores the critical need for targeted infrastructure upgrades and advanced restoration strategies. Furthermore, the proposed framework offers valuable insights for planning and managing resilient power systems in the face of increasing climate stress.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

Particle-in-cell Simulations of Relativistic Magnetic Reconnection with Advanced Maxwell Solver Algorithms

Abstract Relativistic magnetic reconnection is a nonideal plasma process that is a source of nonthermal particle acceleration in many high-energy astrophysical systems. Particle-in-cell (PIC) methods are commonly used for simulating reconnection from first principles. While much progress has been made in understanding the physics of reconnection, especially in 2D, the adoption of advanced algorithms and numerical techniques for efficiently modeling such systems has been limited. With the GPU-accelerated PIC code WarpX, we explore the accuracy and potential performance benefits of two advanced Maxwell solver algorithms: a nonstandard finite-difference scheme (CKC) and an ultrahigh-order pseudo-spectral method (PSATD). We find that, for the relativistic reconnection problem, CKC and PSATD qualitatively and quantitatively match the standard Yee-grid finite-difference method. CKC and PSATD both admit a time step that is 40% longer than that of Yee, resulting in a ∼40% faster time to solution for CKC, but no performance benefit for PSATD when using a current deposition scheme that satisfies Gauss’s law. Relaxing this constraint maintains accuracy and yields a 30% speedup. Unlike Yee and CKC, PSATD is numerically stable at any time step, allowing for a larger time step than with the finite-difference methods. We found that increasing the time step 2.4–3 times over the standard Yee step still yields accurate results, but it only translates to modest performance improvements over CKC, due to the current deposition scheme used with PSATD. Further optimization of this scheme will likely improve the effective performance of PSATD.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

County Level Annual Population Projections for SSP3 and SSP5

This dataset consist of annual (2020-2100) county-level population projections for the United States (U.S.) for Shared Socioeconomic Pathways (SSPs) 3 and 5. The original decadal state-level data is used as an input to the gridded population data model to downscale the state-level projections for each SSP to a 1 km grid. The 1 km data is then aggregated to the county-level using the 2020 TIGER/Line county shapefiles from the U.S. Census Bureau. The two data files are shared as .csv files with the following structure: Rows = Data for individual counties which are identified using their Federal Information Processing Standard (FIPS) code. The FIPS code is stored in the first column. Columns = Each year from 2020-2100 is an individual column. For easier reference, the state name associated with each row is stored in the last column. Values in each cell are the downscaled estimated population for that county-year combination. Values are fractional due to the weighting scheme used in the downscaling. The paper and model source code cited in the "Related Works" section describes the initial source of the population projections.

FIPS↗

DOE OSTI2025

Verification of a Fluid-Based Plasma-Edge Model Within the Multiphysics Object-Oriented Simulation Environment (MOOSE) Framework

As the goal of achieving fusion power on the grid comes closer to fruition, fully coupled multiphysics models of fusion devices will be crucial. These models must incorporate the interconnected phenomena of these devices, including plasma physics, neutronics, first wall interactions, and tritium transport. Currently, there are two main approaches to developing these platforms: (1) loosely coupled, where one couples existing codes and solvers together through input and output parameters and data, and (2) tightly coupled, where one develops the necessary models within a singular, integrated framework. This work focuses on the latter approach for magnetically confined fusion devices by developing a fluid-based plasma-edge model within the Multiphysics Object Oriented Simulation Environment (MOOSE) Framework. This effort is coordinated with other efforts to develop, test, demonstrate, and deploy fusion relevant multiphysics capabilities including electromagnetics, particle-in-cell plasma, tritium transport, and fusion blanket design. This new model is an expansion of the MOOSE-based plasma application, Zapdos, which was originally formulated to model low-temperature, non-magnetized plasma processes. Verification studies have been conducted using newly developed magnetic plasma capabilities. These involved convergence analyses utilizing the method of manufactured solutions to verify new operators and case studies. A modular approach was taken here to demonstrate increasingly complicated simulation scenarios, which included a singular fluid with uniform magnetic field case, a singular fluid with spatially varying magnetic field case, and a coupled multifluid case.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

REopt Lite User Manual

REopt Lite evaluates the economic viability of grid-connected solar photovoltaics, wind, combined heat and power, and electric and thermal storage at commercial and small industrial sites. It allows building owners to identify the system sizes and dispatch strategies that minimize the site’s life cycle cost of energy. REopt Lite also estimates the amount of time on-site generation and storage can sustain the site's critical load during a grid outage and allows the user the choice of optimizing for energy resilience. It is primarily used to inform project development decisions and to support research on the factors that drive project feasibility for market development and policy analysis. It is available through a web interface, application programing interface, and open-source code. This user manual provides an overview of the model, including its capabilities and typical applications; inputs and outputs; economic calculations; technology descriptions; and model parameters, variables, and equations. The model is highly flexible and is continually evolving to meet the needs of each analysis. Therefore, this report is not an exhaustive description of all capabilities, but rather a summary of the core components of the model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

ROOT’s RNTuple I/O Subsystem: The Path to Production

The RNTuple I/O subsystem is ROOT’s future event data file format and access API. It is driven by the expected data volume increase at upcoming HEP experiments, e.g. at the HL-LHC, and recent opportunities in the storage hardware and software landscape such as NVMe drives and distributed object stores. RNTuple is a redesign of the TTree binary format and API and has shown to deliver substantially faster data throughput and better data compression both compared to TTree and to industry standard formats. In order to let HENP computing workflows benefit from RNTuple’s superior performance, however, the I/O stack needs to connect efficiently to the rest of the ecosystem, from grid storage to (distributed) analysis frameworks to (multithreaded) experiment frameworks for reconstruction and ntuple derivation. With the RNTuple binary format soon arriving at its first production release, we present RNTuple’s feature set, integration efforts, and its performance impact on the time-to-solution. We show the latest performance figures of RDataFrame analysis code of realistic complexity, comparing RNTuple and TTree as data sources. We discuss RNTuple’s approach to functionality critical to the HENP I/O (such as multithreaded writes, fast data merging, schema evolution) and we provide an outlook on the road to its use in production.

Blomer, Jakob↗

DOE OSTI2021

Adaptive Control Algorithm to Adjust Settings in Photovoltaic Inverters for Electric Grid Cybersecurity (DERAC) v0.1

Volt-VAR and Volt-Watt functionality in photovoltaic (PV) smart inverters provide mechanisms to ensure system voltage magnitudes and power factors remain within acceptable limits. However, these control functions can become unstable, introducing oscillations in system voltages when not appropriately configured or maliciously altered during a cyberattack. In the event that Volt-VAR and Volt-Watt control functions in a portion of PV smart inverters in a distribution grid are unstable, the proposed adaptation scheme utilizes the remaining and stably-behaving PV smart inverters and other Distributed Energy Resources to mitigate the effect of the instability in real-time. The adaptation mechanism is entirely decentralized, model-free, communication-free, and requires virtually no external configuration. This repository provides code to simulate the algorithm in experiments on the IEEE 37 node test feeder.

Arnold, Daniel↗

DOE OSTI2025

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

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

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

(U) Marker Lagrangian Simulations using Geometry from CAD Assembly Files [Slides]

The Eulerian code PAGOSA now contains an option for the forward time integration of any user selected material in the Lagrangian frame of reference. This option is accomplished using a Multifield version of Brackbill’s FLIP plus Sulsky’s MPM. FLIP+MPM is a Marker Lagrangian scheme that is not to be confused with well known Grid Lagrangian schemes. This Marker Lagrangian method allows materials to undergo completely arbitrary deformation, up to and including fracture and separation. Using this option makes PAGOSA into a mixed frame simulation tool such that fluids can be integrated in the Eulerian frame and solids integrated in the Lagrangian frame. Boundaries between the fields are said to be ’immersed’ in the computational domain, and their nature is formulated using Multifield Theory. Each Marker in a field of FLIP+MPM material tracks the thermodynamic state of a small (but finite) piece of mass. Initiation of the Markers using the geometry and material information from CAD Assembly files is the subject of the talk. The ease and speed with which the dynamics of the assembly can be studied is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

A Structural Design Approach Tailored for the Rapid Preliminary Design of Microreactor Components

High-temperature microreactors can play a role in developing reliable, portable energy sources for off-grid remote locations, microgrid concepts, and industrial process heat. Portability and passive safety criteria tend to skew microreactor structural component designs toward complex geometries, high thermal stresses, and design bases with large numbers of startup/shutdown cycles. Current design rules, as typified by Section III of the American Society of Mechanical Engineers (ASME) Boiler & Pressure Vessel Code, are less than optimal for these conditions, particularly for preliminary component designs where developers need to rapidly consider a large number of potential component configurations. This paper presents a design method targeted toward rapid, efficient evaluation of preliminary component designs using modern finite element analysis. The new method retains key connections with the ASME Code rules and design data while streamlining the design approach. This paper presents the design method, several verification examples illustrating the similarities and differences between the new method and the current ASME rules, and the application of the new approach to the evaluation of a test article mimicking key features of a heat pipe–cooled microreactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

PBX9502 PAGOSA/SURF Calibration for 100 µm Grid Size

There are two methods to calibrate the SURF reactive burn explosives model.The first is to calibrate using Pop plot data. Because there are only two main parameters in SURF, called A and B, there is just enough flexibility in SURF to determine the Pop plot line, i.e., equivalent information to a slope and a single point on the line. However, due to the linearity of the Pop plot being based on logarithmic scales, a slight deviation off the Pop plot (though the fit to the Pop plot might be very good overall) can make considerable difference when comparing with the velocity profile data. The second calibration method thusly focuses strongly on matching entire velocity profile data. This is to say that individual velocity profiles extracted from embedded gages (in the case of PBX 9502, we refer to R. Gustavsen’s experimental data) may fit poorly, while the overall fit to the Pop plot may seem reasonable. In addition, there are other factors that could produce slightly different calibrations for SURF, while leading to an overall acceptable fit to data. One of those may be the physics code used and here we use Pagosa. Another is mesh size dependence. A third one stems from fitting a reactive burn model being done in conjunction with an equation of state (EOS) of the reactants and one for the products – different EOS can result in slightly different parameterizations. Finally, there are several secondary material parameters in SURF, and the overall fit will depend on the entire parameterization. The goal of this report is to clarify some of our previously made statements about SURF calibration, as well as to redo a systematic calibration using 100 µm grid resolution.

42 ENGINEERING↗

DOE OSTI2023

Geothermal Play Fairway Analysis, Part 2: GIS methodology

Play Fairway Analysis (PFA) in geothermal exploration originates from a systematic methodology developed within the petroleum industry and is based on a geologic, geophysical, and hydrologic framework of identified geothermal systems. We tailored this methodology to study the geothermal resource potential of the Snake River Plain and surrounding region, but it can be adapted to other geothermal resource settings. We adapted the PFA approach to geothermal resource exploration by cataloging the critical elements controlling exploitable hydrothermal systems, establishing risk matrices that evaluate these elements in terms of both probability of success and level of knowledge, and building a code-based ‘processing model’ to process results. A geographic information system was used to compile a range of different data types, which we refer to as elements (e.g., faults, vents, heat flow, etc.), with distinct characteristics and measures of confidence. Discontinuous discrete data (points, lines, or polygons) for each element were transformed into continuous interpretive 2D grid surfaces called evidence layers. Because different data types have varying uncertainties, most evidence layers have an accompanying confidence layer which reflects spatial variations in these uncertainties. Risk layers, as defined here, are the product of evidence and confidence layers, and are the building blocks used to construct Common Risk Segment (CRS) maps for heat, permeability, and seal, using a weighted sum for permeability and heat, but a different approach with seal. CRS maps quantify the variable risk associated with each of these critical components. In a final step, the three CRS maps were combined into a Composite Common Risk Segment (CCRS) map, using a modified weighted sum, for results that reveal favorable areas for geothermal exploration. Additional maps are also presented that do not mix contributions from evidence and confidence (to allow an isolated view of evidence and confidence), as well as maps that calculate favorability using the product of components instead of a weighted sum (to highlight where all components are present). Our approach helped to identify areas of high geothermal favorability in the western and central Snake River Plain during the first phase of study and helped identify more precise local drilling targets during the second phase of work. By identifying favorable areas, this methodology can help to reduce uncertainty in geothermal energy exploration and development.

15 GEOTHERMAL ENERGY↗

DOE OSTI2026

Principles and Options for Designing Battery Energy Storage Zoning Ordinances

Deployment of battery energy storage (BESS) systems, both standalone and as part of hybrid systems paired with generation, has rapidly increased in the United States in recent years as utilities and communities have deployed storage to improve electric grid reliability and act as a cost-effective alternative to larger grid infrastructure. The modular nature of BESS technologies means systems may be built near other existing land uses, creating the potential for conflicts with neighboring landowners that can be managed and mitigated through zoning and permitting requirements established by local jurisdictions. While many cities and counties have adopted ordinances specific to BESS into their local zoning codes, these ordinances vary significantly in their requirements and level of detail. Meanwhile, many other jurisdictions, including those home to proposed or existing BESS projects, lack any specific language related to BESS in their zoning codes. Local planning and zoning officials have limited capacity and may lack the familiarity with BESS technologies needed to develop ordinances or otherwise make reasonable zoning decisions that balance safety, community impacts, and other goals. The resulting uncertainty at the local zoning level has led developers to withdraw projects in some areas and has spurred moratoria or bans on energy storage projects in others. This report intends to provide practical resources for practitioners interested in reasonable and effective local regulation of battery energy storage. It does not present a model zoning ordinance, but rather is intended to complement model ordinances developed by others by providing additional context and analysis regarding the structure of energy storage zoning ordinances and the decision points for local officials. Zoning ordinances at the city, town, and county level across the U.S. were surveyed alongside two template model ordinances to identify common elements and options for regulating the zoning and siting of BESS. Common elements identified and analyzed include definitions and general requirements, including cutoffs or tiers used to apply regulations to different system sizes and the permitted zones where jurisdictions allow BESS to be sited; visual, noise, and aesthetic requirements, including property line setbacks, fencing and visual screening, noise, and lighting requirements; and safety and planning requirements, such as site plans, decommissioning plans or funds, and requirements for access by emergency services. The report also summarizes some of the more unique regulations, including those that place additional restrictions on BESS at the local level.

25 ENERGY STORAGE↗

DOE OSTI2021

Atmospheric Transmission Spectroscopy of Hot Jupiter KELT-10b using Synthetic Telluric Correction Software

High-resolution spectroscopic visible data were obtained with the Ultraviolet and Visible Echelle Spectrograph on the Very Large Telescope. Our goal was to analyze the data in an effort to detect the presence of sodium in the atmosphere of hot Jupiter exoplanet KELT-10b, as well as characterize the orbit of the planet via the Rossiter-McLaughlin effect. Eighty spectra were collected during a single transit of KELT-10b. After standard spectroscopic calibration using ESO-Reflex, the synthetic telluric modeling software molecfit was applied to remove terrestrial atmospheric effects, and to refine the wavelength calibration. Sodium is recognized by its characteristic absorption doublet located at 5895.924 and 5889.951 Å, which can be seen in the planet atmosphere transmission spectrum and through excess absorption during the transit. The radial velocity of the host star was analyzed by measuring the average shift of absorption features from spectrum to spectrum. Our results indicate a sodium detection in the planet transmission spectrum with a line contrast of 0.66% and 0.43% ± 0.09% for the sodium DII and DI lines, respectively. Excess absorption measurements agree to within one half combined standard deviation between the planet transmission spectrum (0.143% ± 0.020%, a 7σ detection) and during the time series (0.124% ± 0.034%, a 3.6σ detection) in a band 1.25 Å wide. The wavelength grid corrections provided by molecfit were insufficient to determine radial velocities and measure the Rossiter-McLaughlin effect.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

ATEQ: Adaptive toroidal equilibrium code

A radially adaptive numerical scheme is developed to solve the Grad–Shafranov equation for axisymmetric magnetohydrodynamic equilibrium. Here, a decomposition with independent solutions is employed in the radial direction, and Fourier decomposition is used in the poloidal direction. The independent solutions are then obtained using an adaptive shooting scheme together with the multi-region matching technique in the radial direction. Accordingly, the adaptive toroidal equilibrium (ATEQ) code is constructed for axisymmetric equilibrium studies. The adaptive numerical scheme in the radial direction improves considerably the accuracy of the equilibrium solution. The decomposition with independent solutions effectively reduces the matrix size in solving the magnetohydrodynamic equilibrium problem. The reduction of the matrix size is about an order of magnitude as compared with the conventional radially grid-based numerical schemes. Also, in this ATEQ numerical scheme, no matter how accuracy in the radial direction is imposed, the size of matrices basically does not change. The small matrix size scheme gives ATEQ more flexibility to address the requirement of the number of Fourier components in the poloidal direction in tough equilibrium problems. These two unique features, the adaptive shooting and small matrix size, make ATEQ useful to improve tokamak equilibrium solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

High fidelity actuator line data from 9 turbine wind farm simulations using ExaWind

This data was generated with the ExaWind code suite (https://github.com/Exawind) to investigate the performance of different Active Wake Mixing turbine control in a wind farm situated in a stable atmospheric boundary layer. All cases correspond to a 3x3 wind farm in a 10km x 10km domain using a total mesh size that varied between 1.6 X 10^9 to 1.85 X 10^9 grid cells. The simulations were run across 1800-2000 GPUs on Frontier. The case description and data generation process is fully documented in Yalla, G. R., Brown, K., Cheung, L., Houck, D., deVelder, N., and Balaji, J. (2025). "Estimating annual energy production of wake mixing control strategies including comparisons to wake steering." Wind Energy Sciences (https://doi.org/10.5194/wes-2025-250).

17 WIND ENERGY↗

DOE OSTI2024

Energy Exascale Computational Fluid Dynamics Simulations With the Spectral Element Method

Development and application of the open-source GPU-based fluid-thermal simulation code, NekRS, are described. Time advancement is based on an efficient kth-order accurate timesplit formulation coupled with scalable iterative solvers. Spatial discretization is based on the high-order spectral element method (SEM), which affords the use of fast, low-memory, matrix-free operator evaluation. Further, recent developments include support for nonconforming meshes using overset grids and for GPU-based Lagrangian particle tracking. Results of large-eddy simulations of atmospheric boundary layers for wind-energy applications as well as extensive nuclear energy applications are presented.

42 ENGINEERING↗