Engineering Papers⌕ Search

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.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a 5-dimensional (5D) continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular 2D velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirably large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. We describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma by using 1st- and 2nd-order elements with the pseudo-inverse method and comparing to the bilinear mapping.

Fusion Plasma↗
DOE OSTI2021

Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a five-dimensional continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular two-dimensional velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirably large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. Here we describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma using first- and second-order elements with the pseudo-inverse method and comparing with the bilinear mapping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Powering Post-COVID-19 Resilient Recovery with Green Stimulus

As the coronavirus disease (COVID-19) pandemic spread globally in early 2020, many governments enacted lockdown measures as an initial containment response. These lockdowns, while effective in slowing infection rates, have also had substantial economic consequences. A series of supply and demand shocks, workforce and supply chain disruptions, and other impacts of the pandemic and ensuing lockdowns have upended the livelihoods for much of the world's population. This quick read spotlights how green stimulus plans and other recovery measures can support near-term recovery and enable longer-term power system resilience against future threats. The term “green stimulus” encompasses fiscal measures (i.e., governmental tax and spending actions) that support short-term economic activity that enhances environmental and natural resource quality over a longer term. The deployment of energy efficiency (EE) and renewable energy (RE) technologies, as well as associated initiatives (e.g., grid modernization, electrification of transport, and so on) have been demonstrably effective in providing this kind of short-run economic boost while simultaneously strengthening the resilience of the power sector. Along with green stimulus packages, concurrent action on the regulatory side—for example, the development of green building codes, updating grid interconnection standards, creating programs and marketplaces for demand-side efficiency, and others—can be a force multiplier in ensuring that public spending achieves maximum impact in near-term recovery and long-term power sector transformation

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

An Overview of Behind-the-Meter Solar-Plus-Storage Regulatory Design: Approaches and Case Studies to Inform International Applications

Behind-the-meter energy storage systems paired with distributed photovoltaic (DPV) - with the capability to act as both generation and load - represent a unique and disruptive power sector technology capable of providing a range of important services to customers, utilities, and the broader power system. How should regulators, utilities, and policymakers manage the range of challenges and opportunities that increased behind-the-meter energy storage deployment will bring to the power system, in particular when these systems are paired with DPV? This report is intended to offer key regulatory considerations for facilitating DPV-plus-storage programs for retail customers; relevant cases from U.S. states are provided as examples of how novel regulatory issues related to behind-the-meter energy storage systems paired with DPV are being addressed in practice. At a high level, designing a regulatory framework that aligns DPV-plus-storage deployment with larger policy objectives requires thoughtful deliberation across a range of technical and economic issues. While this report attempts to segment many of these issues into distinct topics to enhance reader understanding, in reality, DPV-plus-storage regulatory issues are all closely integrated, and design decisions on a single aspect often have broader implications. With that in mind, this report outlines a series of steps that can be employed by regulators to approach DPV-plus-storage regulatory design.

14 SOLAR ENERGY↗
DOE OSTI2020

An Overview of Behind-the-Meter Solar-Plus-Storage Program Design: With Considerations for India

Indian consumers have been deploying behind-the-meter generation (predominantly diesel backup, and, more recently, photovoltaic) and storage systems (predominantly lead-acid and other kinds of batteries as uninterrupted power supplies) by the millions for decades) to address reliability issues within the Indian power system. The same energy storage systems could provide additional services to the consumer and distribution companies if properly regulated and designed from the outset to be grid interactive. Grid-connected distributed solar PV (DPV), or rooftop solar, has also seen wide deployment in India and features prominently in the Government of India’s plans for a transition to clean, reliable, and affordable energy for all. At the same time, many utilities and state governments, as well as the central government in India are currently funding-constrained for both operational and future capital expenditures in the power sector, and some perceive customer-sited resources as exacerbating existing financial challenges. In that context, behind-the-meter energy storage systems paired with distributed photovoltaic (DPV) - with the capability to act as both generation and load - represent a potentially unique and disruptive power sector technology capable of providing a range of important services to customers, utilities, and the broader power system in India. This report aims to offer a comprehensive, evidence-based approach to designing customer programs based on experience in the United States that can help regulators, utilities, and policymakers in India manage the range of challenges and opportunities that increased behind-the-meter energy storage deployment will bring to the power system, in particular when these systems are paired with DPV.

14 SOLAR ENERGY↗
DOE OSTI2024

code for "Improving the QBO forcing by resolved waves with vertical grid refinement in E3SMv2"

This is all code that used in the journal article "Improving the QBO forcing by resolved waves with vertical grid refinement in E3SMv2", including: Running and post processing E3SM: Run E3SM: run_E3SM.2023.scidac.MMF.amip.py Regrid: regrid.E3SM.py post processing, diagnostics, archive: run_post.2023.scidac.MMF.py TEM calculation, including filter TEM, E3SM: calculate.TEM.v2.py calculate.TEM.monthly_mean.v1.py TEM, ERA5: ERA5_calculate_TEM_90.py merge TEM data: merge_TEM_E3SM_ERA5.py TEM, filter, E3SM: filter_TEM_E3SM.py TEM, filter, ERA5: filter_TEM_ERA5.py WK wave analysis code from (https://github.com/brianpm/wavenumber_frequency): wavenumber_frequency_functions.py calculate the wave spectra so easy to plot: wk_calculation.py Make figures some functions for plot: plot.py code for drawing paper plots: figures_for_paper.ipynb

Hannah, WalterM↗
DOE OSTI2024

Towards exascale for wind energy simulations

We examine large-eddy-simulation modeling approaches and computational performance of two open-source computational fluid dynamics codes for the simulation of atmospheric boundary layer flows that are of direct relevance to wind energy production. The first code, NekRS, is a high-order, unstructured-grid, spectral element code. The second code, AMR-Wind, is a second-order, block-structured, finite-volume code with adaptive mesh refinement capabilities. The objective of this study is to co-develop these codes in order to improve model fidelity and performance for each. These features will be critical for running ABL-based applications such as wind farm analysis on advanced computing architectures. To this end, we investigate the performance of NekRS and AMR-Wind on the Oak Ridge Leadership Facility supercomputers Summit, using 4 to 800 nodes (24 to 4,800 NVIDIA V100 GPUs), and Crusher, the testbed for the Frontier exascale system, using 18 to 384 Graphics Compute Dies on AMD MI250X GPUs. We compare strong- and weak-scaling capabilities, linear solver performance, and time to solution. We also identify leading inhibitors to parallel scaling.

17 WIND ENERGY↗
DOE OSTI2022

INGRID: An interactive grid generator for 2D edge plasma modeling

A fusion boundary-plasma domain is defined by axisymmetric magnetic surfaces where the geometry is often complicated by the presence of one or more X-points; and modeling boundary plasmas usually relies on computational grids that account for the magnetic field geometry. The new grid generator INGRID (Interactive Grid Generator) presented here is a Python-based code for calculating grids for fusion boundary plasma modeling, for a variety of configurations with one or two X-points in the domain. INGRID first performs partitioning over the domain consisting of a small number of patches conforming to the magnetic field and wall geometry; then it generates a subgrid on each of the patches and joins them into a global grid. This domain partitioning strategy makes possible a uniform treatment of various configurations with one or two X-points in the domain. This includes single-null, double-null, and other configurations with two X-points in the domain. The INGRID design allows generating grids either interactively, via a parameterfile driven GUI, or using a non-interactive script-controlled workflow. Results of testing demonstrate that INGRID is a flexible, robust, and user-friendly grid-generation tool for fusion boundary-plasma modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

ExaWind: Exascale Predictive Wind Plant Flow Physics Modeling

The scientific goal of the ExaWind project is to advance our fundamental understanding of the flow physics governing whole wind plant performance, including wake formation, complex terrain impacts, and turbine-turbine-interaction effects. The primary application codes in the ExaWind environment are Nalu-Wind, an unstructured-grid computational fluid dynamics (CFD) code, AMR-Wind, a structured-grid CFD code, and OpenFAST, a whole-turbine simulation code. In this poster we present our progress towards simulating the ExaWind challenge problem, which is a predictive simulation of a wind farm with tens of megawatt-scale wind turbines dispersed over an area of 50 square kilometers.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2023

ExaWind: Then and Now

The scientific goal of the ExaWind project is to advance our fundamental understanding of the flow physics governing whole wind plant performance, including wake formation, complex terrain impacts, and turbine-turbine-interaction effects. The primary application codes in the ExaWind environment are Nalu-Wind, an unstructured-grid computational fluid dynamics (CFD) code, AMR-Wind, a structured-grid CFD code, and OpenFAST, a whole-turbine simulation code. In this poster we present the current status of the ExaWind software stack in the context of the modeling and simulation capabilities when the project started in 2016.

computational fluid dynamics↗
DOE OSTI2026

Feasibility Study on Implementing a Staggered-Grid Finite Volume Method for System Analysis Code Development Under the MOOSE Framework

Here, this work summarizes a feasibility study on testing numerical algorithms that are suitable and efficient for advanced system analysis code development under the mutli-physics framework, MOOSE. The key to the test bed is the implementation of high-order one-dimensional staggered-grid finite volume method (SG-FVM), and its direct interaction with the linear/nonlinear solver, PETSc. The test bed utilized a more flexible code structure to enable the finite volume method implementation and direct interacting with the solver package, instead of using the natively supported finite element method by the framework. Using a suite of selected test problems with different problem sizes and levels of complexity, the implemented SG-FVM demonstrated superior performance improvement against a direct finite element method implementation through MOOSE. On two computer systems, the speedup was observed to be significant, with at least one order of magnitude of solving time reduction. For a complex reactor model, transient simulation was performed using the newly developed finite volume method code, the results of which agree very well with the reference results from the finite element method code. Overall, this study demonstrates a successful feasibility study on the proposed numerical algorithms and software structure to support advanced system analysis tool development.

MOOSE↗
DOE OSTI2022

INGRID: An interactive grid generator for 2D edge plasma modeling

A fusion boundary-plasma domain is defined by axisymmetric magnetic surfaces where the geometry is often complicated by the presence of one or more X-points; and modeling boundary plasmas usually relies on computational grids that account for the magnetic field geometry. The new grid generator INGRID (Interactive Grid Generator) presented in this work is a Python-based code for calculating grids for fusion boundary plasma modeling, for a variety of configurations with one or two X-points in the domain. INGRID first performs partitioning over the domain consisting of a small number of patches conforming to the magnetic field and wall geometry; then it generates a subgrid on each of the patches and joins them into a global grid. This domain partitioning strategy makes possible a uniform treatment of various configurations with one or two X-points in the domain. This includes single-null, double-null, and other configurations with two X-points in the domain. The INGRID design allows generating grids either interactively, via a parameter-file driven GUI, or using a non-interactive script-controlled workflow. Results of testing demonstrate that INGRID is a flexible, robust, and user-friendly grid-generation tool for fusion boundary-plasma modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

EVs-at-RISC: A Secure and Resilient Interoperable SCM Control System Architecture for Electric Vehicle’s-at-Scale (Final Technical Report)

The EVs-at-RISC project was a five-year research, development, and demonstration initiative to create foundational tools for utility-scale fleet aggregation and Smart Charge Management (SCM) of Electric Vehicles (EV), Electric Vehicle Charging Infrastructure (EVCI), and related Distributed Energy Resources (DER). Rather than seeking to develop and demonstrate highly perfected SCM algorithms and control strategies, this project instead focused on creating foundational software solutions that enable unprecedented digital interoperability across the communications technologies and vendor platforms used to manage EV , EVCI, and DER, as well as existing energy management infrastructure operated by utilities, grid operators, and aggregators. This project then extends these novel interoperability capabilities to develop and deploy powerful middleware abstractions across grid edge networks and EVCI/DER fleet aggregations incorporating modern software tools and best practices, such as CI/CD, to bring the immense capabilities of infrastructure-as-code and policy-as-code to modern grid edge network environments. This addresses the foremost systemic issues preventing realization of any net operational benefits from scaled deployment of behind-the-meter EV, EVCI, and DER assets in electric power grids and markets today. The results of this approach and project unlock massive potential for new SCM capabilities to be easily prototyped, evaluated, and deployed at-scale within the existing grid edge network infrastructure and EVCI/DER technology ecosystem. The EVs-at-RISC project achieves this by extending Open Field Message Bus (OpenFMB), a conceptual model for digital interoperability and distributed intelligence in traditional front-of-meter utility SCADA networks, validating our hypothesis that OpenFMB could be similarly used to solve systemic digital interoperability issues in behind-the-meter environments and unlock real-world utility-scale SCM capabilities without requiring any new proprietary vendor solutions or significant infrastructure reconfiguration.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

pyRMG: A framework for high-throughput, large-cell DFT calculations on supercomputers

Exascale computing delivers the raw power to simulate ever larger and more chemically realistic systems, but realizing this potential requires codes that can efficiently use thousands of processors. Our real-space multigrid (RMG) density functional theory (DFT) code’s grid-decomposition approach scales nearly linearly with the number of graphics processing units (GPUs), even for simulations exceeding thousands of atoms. This scalability makes RMG a compelling tool for high-throughput DFT studies of materials that would otherwise be bottlenecked in other codes (for example, by global fast Fourier transforms in plane-wave DFT). However, the limited workflow infrastructure for RMG has thus far constrained its adoption to a small user community. In this work, we present pyRMG, a Python package designed to streamline the setup and execution of RMG DFT calculations. Built on the pymatgen and ASE (Atomic Simulation Environment) computational materials science Python packages, pyRMG automates input generation and convergence checking, and it integrates with modern job schedulers (e.g., Flux) on leadership-class platforms such as Frontier and Perlmutter. Here, we demonstrate pyRMG for a high-throughput study of strain effects in 2D 2L-Bi 2 Se 3 /2L-NbSe 2 heterostructures, which offers chemical insights into this system and shows that RMG-based workflows can converge with limited user intervention.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Rendezvous algorithms for large-scale modeling and simulation

Rendezvous algorithms encode a communication pattern that is useful when processors sending data do not know who the receiving processors should be, or vice versa. The idea is to define an intermediate decomposition where datums from different sending processors can ”rendezvous” to perform a computation, in a manner that both the senders and eventual receivers of the results can identify the appropriate rendezvous processor. Though they were originally designed for interpolating between overlaid grids with independent parallel decompositions (Plimpton et al., 2004), we have recently found rendezvous algorithms useful for a variety of operations in particle- or grid-based simulation codes when running large problems on large numbers of processors. In particular, we show they can perform well when a load-balanced intermediate decomposition is randomized and not spatial, requiring all-to-all communication to move data between processors. In this case rendezvous algorithms leverage the large bisection communication bandwidths which parallel machines provide. We describe how rendezvous algorithms work in a scientific computing context and give specific examples for molecular dynamics and Direct Simulation Monte Carlo codes which result in dramatic performance improvements versus simpler algorithms which do not scale as well. We explain how a generic rendezvous algorithm can be implemented, and also point out similarities with the MapReduce paradigm popularized by Google and Hadoop.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Computational Evaluation of Thermal Response of Open-Cell Foam With Circular Pore (Computional Evaluation of Thermal Response of Open-Cell Foam with Circular Pore)

The evaluation of effective material properties in heterogeneous materials (e.g., composites or multicomponent structures) critically relevant to a wide spectrum of applications, including nuclear power, electronic packaging, flame retardants, hypersonics, and gas turbine power. The work described in this paper is centered around the numerical assessment of the thermal behavior of porous materials obtained from finite element thermal modeling and simulation. Here, two-dimensional, steady state analyses were performed on unit cells with centered, circular pores using a second order accurate Galerkin finite element method (FEM). The effective thermal conductivities of the porous systems were examined, encompassing a range of porosities from 4.9% to 60.1%. The geometries of the models were generated based on ordered circular pores for each modeled porosity level. The system response quantity (SRQ) under investigation was the dimensionless effective thermal conductivity across the unit cell. The dimensionless effective thermal conductivity was compared across all simulated cases, producing a trend between porosity and effective thermal conductivity. In the presented investigation, the method of manufactured solutions (MMS) was used to perform code verification, and the grid convergence index (GCI) was employed to estimate discretization uncertainty as solution verification. Code verification concluded an approximately second order accurate Galerkin FEM solver. It was found that the introduction of porosity to the unit cell material structure reduces effective thermal conductivity, as anticipated. Numerical results obtained in this study are compared to an analytical solution and to a sample of empirical data.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dust Survival in Galactic Winds

This repository contains three-dimensional volumetric data from an Eulerian hydrodynamical simulation (conducted on a uniform Cartesian grid) generated by the Cholla hydrodynamics code. The datasets contain snapshots (full-grid, projections, and slices) in the HDF5 format of a multi-phase medium in which a hot, diffuse, dust-free background wind accelerates a cool, dense cloud of gas and dust. This scenario is intended to represent a supernova-driven galactic outflow, in which hot supernova winds are thought to accelerate cool interstellar medium material out of the galactic disk into the surrounding circumgalactic medium. There are three separate datasets for simulations corresponding to three cloud evolutionary scenarios: long-term cloud survival (surv), marginal cloud survival (disr), and cloud destruction (dest). Projection and slice images of the simulations are also included in this repository.

79 ASTRONOMY AND ASTROPHYSICS↗