Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Grid computing”

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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Renewable energy integration and system operation challenge: control and optimization of millions of devices

The electric power infrastructure, originally designed and built on large-scale power plants, is evolving into a more resilient power generation and delivery system in which millions of smaller units of distributed energy generation resources units will be installed in sub-transmission and distribution networks. In order to control, manage and optimize the future grid, a hierarchical design is presented in this chapter which enables the distributed control on grid edge while inheriting the existing centralized control structure. This layered design of large-scale power system operation and control uses the following principle: reactive power control is treated as a primary control for voltage stability, and the real power control is primarily a grid-level control but can also be a supplementary control for voltage support in the case of insufficient reactive power control capacity. For the purpose of active control and operation at the distribution level, a recursive power network model is derived from nodal injection and branch power flow models. Based on the model, the proposed algorithms of hierarchical control, grid-edge inference and dynamic hosting allowance are developed and presented for multi-level controlled operation. And, a co-simulation architecture of integrated T&D system is presented to validate and demonstrate the feasibility and scalability of proposed algorithms.

Xu, Ying↗
DOE OSTI2025

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗
DOE OSTI2022

Physics makes the difference: Bayesian optimization and active learning via augmented Gaussian process

Abstract Both experimental and computational methods for the exploration of structure, functionality, and properties of materials often necessitate the search across broad parameter spaces to discover optimal experimental conditions and regions of interest in the image space or parameter space of computational models. The direct grid search of the parameter space tends to be extremely time-consuming, leading to the development of strategies balancing exploration of unknown parameter spaces and exploitation towards required performance metrics. However, classical Bayesian optimization (BO) strategies based on the Gaussian process (GP) do not readily allow for the incorporation of the known physical behaviors or past knowledge. Here we explore a hybrid optimization/exploration algorithm created by augmenting the standard GP with a structured probabilistic model of the expected system’s behavior. This approach balances the flexibility of the non-parametric GP approach with a rigid structure of physical knowledge encoded into the parametric model. The fully Bayesian treatment of the latter allows additional control over the optimization via the selection of priors for the model parameters. The method is demonstrated for a noisy version of a standard univariate test function used to evaluate optimization algorithms and further extended to physical lattice models. This methodology is expected to be universally suitable for injecting prior knowledge in the form of physical models and past data in the BO framework.

42 ENGINEERING↗
DOE OSTI2021

Seismic waves in medium with poroelastic/elastic interfaces: a two-dimensional P-SV finite-difference modelling

SUMMARY We present a new methodology of the finite-difference (FD) modelling of seismic wave propagation in a strongly heterogeneous medium composed of poroelastic (P) and (strictly) elastic (E) parts. The medium can include P/P, P/E and E/E material interfaces of arbitrary shapes. The poroelastic part can be with (i) zero resistive friction, (ii) non-zero constant resistive friction or (iii) JKD model of the frequency-dependent permeability and resistive friction. Our FD scheme is capable of subcell resolution: a material interface can have an arbitrary position in the spatial grid. The scheme keeps computational efficiency of the scheme for a smoothly and weakly heterogeneous medium (medium without material interfaces). Numerical tests against independent analytical, semi-analytical and spectral-element methods prove the efficiency and accuracy of our FD modelling. In numerical examples, we indicate effect of the P/E interfaces for the poroelastic medium with a constant resistive friction and medium with the JKD model of the frequency-dependent permeability and resistive friction. We address the 2-D P-SV problem. The approach can be readily extended to the 3-D problem.

58 GEOSCIENCES↗
DOE OSTI2026

Pixel-Resolved Long-Context Learning for Turbulence at Exascale: Resolving Small-scale Eddies Toward the Viscous Limit

Turbulence plays a crucial role in multiphysics applications, including aerodynamics, fusion, and combustion. Accurately capturing turbulence's multiscale characteristics is essential for reliable predictions of multiphysics interactions, but remains a grand challenge even for exascale supercomputers and advanced deep learning models. The extreme-resolution data required to represent turbulence, ranging from billions to trillions of grid points, pose prohibitive computational costs for models based on architectures like vision transformers. To address this challenge, we introduce a multiscale hierarchical Turbulence Transformer that reduces sequence length from billions to a few millions and a novel RingX sequence parallelism approach that enables scalable long-context learning. We perform scaling and science runs on the Frontier supercomputer. Our approach demonstrates excellent performance up to 1.1 EFLOPS on 32,768 AMD GPUs, with a scaling efficiency of 94\%. To our knowledge, this is the first AI model for turbulence that can capture small-scale eddies down to the dissipative range in three-dimensional turbulence at high Reynolds numbers.

Yin, Junqi [ORNL] (ORCID:0000000338435520)↗
DOE OSTI2023

Automating the Solar Interconnection Technical Evaluation Process: PREconfiguring and Controlling Inverter SEt-Points (PRECISE)

Utilities are receiving an increasing number of interconnection requests for distributed solar photovoltaic (PV) systems from their customers, requiring a solution to quickly and intelligently perform technical assessments of these requests and to manage these new renewable energy assets. In this paper, we present a stand-alone interconnection tool, PREconfiguring and Controlling Inverter SEt-points (PRECISETM), and how it overcomes the challenges of accelerating technical evaluations and leverages smart inverter functionality to provide local voltage support on an as-needed basis. The paper details the purpose and architecture of PRECISE, including tools, data source requirements, the model-based interconnection evaluation process, and technical results. PRECISE is an integrated software tool that customizes PV inverter settings for utilities by modeling and assessing impacts on local voltages and the need for advanced inverter functions (AIFs) (e.g., volt-VAR and volt-watt). PRECISE maximizes the use of readily available utility data sets to make fast online assessments of incoming PV interconnection requests and makes accept/reject recommendations along with custom settings for AIFs.

ENERGY PLANNING, POLICY, AND ECONOMY,SOLAR ENERGY↗
DOE OSTI2021

An Energy Service Interface for Distributed Energy Resources

Renewable energy resources, particularly wind and solar photovoltaic, are becoming significant contributors to electric power generation. These re-sources will contribute towards achieving sustainable electric power systems. However, renewable resources will dramatically increase the demand for flexible power system operations. This paper proposes an energy service interface that will allow aggregated distributed energy resources, such as residential loads and inverter-based systems, to participate in NERC-defined smart energy reliability services. Such cyber-physical systems will increase system flexibility by ensuring match between energy supply and energy demand.Aggregation and coordinated dispatch of millions of distributed energy resources will require development of large-scale computing networks. Several smart grid interface-enabling technologies, including IEEE 2030.5, Common Smart Inverter Profile, SunSpec Modbus, and CTA 2045, are discussed. Residential loads are categorized by their static and dynamic energy characteristics to identify services in which they can participate. The business model for the energy services interface as well as probabilistic modeling for resource estimation are highlighted as future considerations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Selecting Critical Scenarios of DER Adoption in Distribution Grids Using Bayesian Optimization

We develop a new methodology to select scenarios of DER adoption most critical for distribution grids. Anticipating risks of future voltage and line flow violations due to additional PV adopters is central for utility investment planning but continues to rely on deterministic or ad hoc scenario selection. We propose a highly efficient search framework based on multi-objective Bayesian Optimization. We treat underlying grid stress metrics as computationally expensive black-box functions, approximated via Gaussian Process surrogates and design an acquisition function based on probability of scenarios being Pareto-critical across a collection of line- and bus-based violation objectives. Our approach provides a statistical guarantee and offers an order of magnitude speed-up relative to a conservative exhaustive search. Case studies on realistic feeders with 200-400 buses demonstrate the effectiveness and accuracy of our approach.

Mulkin, Olivier↗
DOE OSTI2022

Game-Theoretic Strategies for Cyber-Physical Infrastructures Under Component Disruptions

In this work, networked infrastructures of recursively defined systems composed of discrete cyber and physical components are considered. The components of basic systems at the finest levels can be disrupted by cyber or physical means, and can be reinforced to survive at certain costs. A problem of ensuring the infrastructure performance is formulated as a game between a provider and an attacker, who probabilistically choose components to reinforce and attack, respectively. The disruptions of this infrastructure are characterized using the aggregate failure correlation function that specifies the conditional failure probability of the infrastructure given the failure of an individual system at that level. The survival probabilities of basic systems satisfy simple product-form, first-order differential equations expressed in terms of the multiplier functions. The utility functions of the provider and attacker are composed of the reward and cost terms, both expressed in terms of the component reinforcement and attack probabilities. The Nash equilibrium of this game is characterized, along with the sensitivity functions of the survival probabilities of basic systems that highlight their dependence individually on the cost-benefit terms, the correlation functions, and the multiplier functions. These results are illustrated using simplified models of a distributed cloud servers infrastructure, a 5G data network infrastructure, a high performance computing federation, and a smart energy grid infrastructure.

42 ENGINEERING↗
DOE OSTI2024

SDE_quark

This program provides a fixed-grid quadrature algorithm to compute integrals in the self-energy of the quark propagator within the Maris-Tandy model. The quark propagator both on t he spacelike real axis and at complex-valued momenta is determined from its Schwinger-Dyson equation (SDE). We first apply an iterative solver to find the quark propagator on the spacelike real axis. The propagator at complex-valued momenta is then computed from its self-energy based on this solution, where demanding integrals are encountered. In order to compute of these integrals, we apply customized variable transformations for the radial integral after subtracting the asymptotics. We subsequently apply an optional compound of quadrature rules for the angular integral.

Jia, Shaoyang↗
DOE OSTI2021

Energy Exascale Earth System Model v2.0

First release of version 2 of the Energy Exascale Earth System Model. The atmosphere component remains EAM. Major changes since version 1 include: all column-physics parameterizations are computed on a separate grid that has approximately half the number of points of the dynamics grid, a new nonhydrostatic dynamical core (running in hydrostatic mode) with semi-Lagrangian tracer transport, CLUBB updated from v1 to v2, a new convective trigger (dCAPE/ULL) based on the dynamic Convective Available Potential Energy (CAPE) (dCAPE) and the Unrestricted Launch Level (ULL) concepts is used in ZM. minimum cloud droplet number changed, gravity wave drag energy conservation fixed and new tunings used, dust emission size distribution changed to emit more coarse dust particles The land component is still ELM. Major changes since version 1 include: using SNICAR-AD for radiation in snow to match the sea-ice model and fixing bugs in snow compaction and water state calculation. The ocean component remains MPAS-ocean. Major change since version 1 include: Redi isopycnal mixing has been updated, tested, and tuned in combination with the Gent-McWilliams parameterization, a sign error was fixed in the 3rd-order flux routines, the mesh used in low-resolution coupled cases was modified and the time steps adjusted, new regionally refined meshes were created, one focused on North America and another on the Southern Ocean, including ice shelf cavities. The sea-ice component remains MPAS-seaice. Major changes since version 1 include: A new heat- and freshwater-conserving coupling of frazil ice, turning off of SSH filtering, addition of SNICAR-AD and snow grain aging. The land-ice component remains MPAS-Albany-landIce (MALI) and is a static ice sheet. There is more out-of-the-box support for cryosphere configurations including the new regionally refined configuration around Antarctica with ice shelf cavities. The river model is MOSART. The half-degree river mesh was redone so it no longer treats Black and Caspian seas as ocean. The coupler remains cpl7/MCT. Major changes since version 1 include handling of ice shelf melt fluxes (heat / freshwater exchange with the ocean), and data icebergs. All components allow regional refinement of their meshes and two separate example of refinement, one in and around North America and one around Antarctica and the Southern Ocean, are provided.

E3SM Project, DOE↗
DOE OSTI2021

Combinatorial Algorithms in Scientific Computing

We provide the final report for this grant, detailing the publications, software produced, students trained who have joined the DOE workforce, and the impact our work has had on computational mathematics and related disciplines.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Multi-Stage and Multi-Timescale Robust Co-Optimization Planning for Reliable and Sustainable Power Systems. Final Report

In this project, Clarkson University, in collaboration with Southern Methodist University and University of Pittsburgh, conducted the research to model, design, and implement a sophisticated generation and transmission co-optimization planning decision tool, called Multi-stage and Multi-timescale robust Co-Optimization Planning (MMCOP). The MMCOP decision tool intends to facilitate generation and transmission co-optimization planning of emerging power systems, while mitigating risks and uncertainties in both short-term operation dynamics and long-term policy and technology changes. Long-term power system planning aims at optimizing asset utilization by investing in a proper mix of various generation technologies and transmission lines to supply the future load growth. In particular, the Clean Power Plan (CPP), which is designed to combat climate change and reduce carbon emissions by setting a national limit on carbon pollution from power plants, may dramatically change the landscape of the power industry by further promoting clean energy and phasing out emissions-intensive generation technologies. In addition, novel non-wire alternatives (e.g., demand response (DR), distributed generation (DG), energy efficiency (EE), and smart grid technologies) and the computational complexity for large-scale systems significantly complicate the system planning procedure even further. However, existing conventional planning approaches neglect short-term variability and uncertainty of renewable energy, hourly chronological operation details, and physical nonlinear characteristics of the alternating current transmission network. As a result, existing conventional planning approaches may not work properly, and power systems reliability could be in jeopardy. In observing the limitations of existing conventional planning approaches and addressing new challenges of emerging power systems, the main scope of this project to develop co-optimization planning models within a multi-stage and multi-timescale framework. In particular, random contingencies, key uncertainty factors, and AC power flows are included to derive expansion plans while considering both long-term reliability and short-term flexibility.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

BACCHUS Analysis of Weak Lines in APOGEE Spectra (BAWLAS)

Elements with weak and blended spectral features in stellar spectra are challenging to measure and require specialized analysis methods to precisely measure their chemical abundances. In this work, we have created a catalog of approximately 120,000 giants with high signal-to-noise Apache Point Observatory Galactic Evolution Experiment (APOGEE) Data Release 17 (DR17) spectra, for which we explore weak and blended species to measure Na, P, S, V, Cu, Ce, and Nd abundances and 12 C/ 13 C isotopic ratios. We employ an updated version of the Brussels Automatic Code for Characterizing High-accuracy Spectra (BACCHUS) code to derive these abundances using the stellar parameters measured by APOGEE's DR17 Stellar Parameters and Chemical Abundances Pipeline, quality flagging to identify suspect spectral lines, and a prescription for upper limits. Combined, these allow us to provide our BACCHUS Analysis of Weak Lines in APOGEE Spectra catalog of precise chemical abundances for these weak and blended species, which agrees well with the literature and improves upon APOGEE abundances for these elements, some of which are unable to be measured with APOGEE's current, grid-based approach without computationally expensive expansions. This new catalog can be used alongside APOGEE and provides measurements for many scientific applications ranging from nuclear physics to Galactic chemical evolution and Milky Way population studies. To illustrate this we show some examples of uses for this catalog, such as showing that we observe stars with enhanced s-process abundances or that we can use the 12 C/ 13 C ratios to explore extra mixing along the red giant branch.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A Systematic Framework for Tuning Open-Source Multifunctional IBR Models To Emulate OEM Black-Box Fault Dynamics

This paper presents a systematic framework to tune a generic IBR EMT model to match with an OEM provided balckbox inverter model based on the fault current responses. The key learnings and findings are summarized as follows: The tunable key parameters include inner control loops and current limiters to align the fault current magnitude, sequence content, and phase trajectories with the OEM models across diverse fault type and locations. The tuned model's fidelity is validated through comparative analysis with an OEM blackbox model, assessing both the fault current response and the responses of multiple relay elements. The results demonstrate the tuned generic model can trigger relay decision logic that is identical or near identical to that of the OEM model, thus generating very good match model for fault studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1974

Development of inverse inviscid transonic solution methods

A numerical method suitable for the analysis and/or design of supercritical transonic airfoils is reported. In order to achieve accuracy, the method utilizes the full inviscid potential flow equations; and in order to remain simple it solves the problem in a stretched Cartesian grid system. The resulting computer program has several advantages over others of its type -- its use in either the direct analysis mode in which the airfoil shape is prescribed and the flow field and surface pressures are determined, or in the inverse mode in which the surface pressures are given and the airfoil shape and flow field are computed. Other advantages of the program include its use in a design program, the rotated finite difference scheme and its determination of the airfoil shape simultaneously with the flow field relaxation solution.

Carlson, L. A.↗
NASA NTRS1975

Recent improvements to BANDIT

The NASTRAN preprocessor BANDIT, which improves NASTRAN's computer efficiency by resequencing grid point labels for reduced matrix bandwidth is described. The addition of (1) the Gibbs-Poole-Stockmeyer (GPS) algorithm, and (2) the user option to reduce matrix profile rather than matrix bandwidth is also described. It is shown that, compared to the Cuthill-McKee algorithm on which BANDIT was originally based, GPS is faster and achieves similar results. Current capabilities and options of BANDIT are summarized.

Everstine, G. C.↗