Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “IMPLEMENTATION”

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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

VAN-DAMME: GPU-accelerated and symmetry-assisted quantum optimal control of multi-qubit systems

We present an open-source software package, VAN-DAMME (Versatile Approaches to Numerically Design, Accelerate, and Manipulate Magnetic Excitations), for massively-parallelized quantum optimal control (QOC) calculations of multi-qubit systems. To enable large QOC calculations, the VAN-DAMME software package utilizes symmetry-based techniques with custom GPU-enhanced algorithms. This combined approach allows for the simultaneous computation of hundreds of matrix exponential propagators that efficiently leverage the intra-GPU parallelism found in high-performance GPUs. In addition, to maximize the computational efficiency of the VAN-DAMME code, we carried out several extensive tests on data layout, computational complexity, memory requirements, and performance. These extensive analyses allowed us to develop computationally efficient approaches for evaluating complex-valued matrix exponential propagators based on Padé approximants. To assess the computational performance of our GPU-accelerated VAN-DAMME code, we carried out QOC calculations of systems containing 10 - 15 qubits, which showed that our GPU implementation is 18.4× faster than the corresponding CPU implementation. Our GPU-accelerated enhancements allow efficient calculations of multi-qubit systems, which can be used for the efficient implementation of QOC applications across multiple domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Record acceleration of the two-dimensional Ising model using a high-performance wafer-scale engine

The versatility and wide-ranging applicability of the Ising model, originally introduced to study phase transitions in magnetic materials, have made it a cornerstone in statistical physics and a valuable tool for evaluating the performance of emerging computer hardware. Here, we present a novel implementation of the two-dimensional Ising model on Cerebras Wafer-Scale Engine (WSE) – a revolutionary processor that is opening new frontiers in computing. In our deployment of the checkerboard algorithm, we optimized the Ising model to take advantage of the unique WSE architecture. Specifically, we employed a compressed bit representation storing 16 spins on each int16 word, and efficiently distributed the spins over the processing units enabling seamless weak scaling and limiting communications to only immediate neighboring units. Our implementation can handle up to 754 simulations in parallel, achieving an aggregate of over 61.8 trillion flip attempts per second for Ising models with up to 200 million spins. This represents a gain of up to 148 times over previously reported single-devices with a highly optimized implementation on NVIDIA V100 and up to 88 times in productivity compared to NVIDIA H100. Our findings highlight the significant potential of the WSE in scientific computing, particularly in the field of materials modeling.

Ising model↗
DOE OSTI2023

Evaluation of riparian enhancement actions in the Columbia River Basin

Riparian enhancement is a common restoration technique in the Columbia River Basin (CRB) and the Pacific Northwest. However, relatively few studies have evaluated its success and even fewer studies include long-term monitoring. Forty-one riparian planting projects located in the CRB, each with a paired treatment and control reach, were evaluated in the summer of 2018 and 2019 using an extensive post-treatment (EPT) design. At each reach, we quantified woody plant abundance, richness, diversity, and vegetation and canopy cover as riparian response variables, and measured terrace height and documented individuals with bud browse, deceased, and with predator protection as potential explanatory variables. Species richness, woody plant abundance (all height classes combined and within shrub height class), and the proportion of woody plants with bud browse or deceased were all higher in treatment than control reaches. However, no significant improvements were observed for any other riparian response or explanatory variables. Implementation techniques, including invasive species removal efforts, follow-up planting, watering, and the type of restoration project (floodplain restoration and planting, in-stream restoration and planting, or solely planting) were also considered as potential explanatory variables and were found to have significant impacts on restoration response for abundance and cover variables. Further, our results suggest that planting implementation methods, site level physical factors, and time since planting all influence the success of riparian planting projects and if not addressed in restoration design and implementation can contribute to a lack of detectable response.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Deploying advanced supervisory control strategies to small- and medium-sized commercial buildings: Case study and lessons learned

Deploying advanced supervisory control strategies (ASCSs) in small- and medium-sized commercial buildings (SMBs) is vital but faces two issues: (1) a lack of building control and communication infrastructure (BCCI) in SMBs and (2) the significant engineering efforts required to implement and configurate ASCSs. Despite hindering the large-scale adoption of ASCSs in SMBs, these issues have not been adequately explored in the literature, which tends to focus more on feasibility than scalability. This paper provides a comprehensive evaluation of these two issues through a case study of an occupied SMB in eastern Tennessee of the United States. Specifically, we design and implement a BCCI for the studied building with commercial off-the-shelf products to accommodate the needs for deploying ASCSs. Here, we then deploy two ASCSs—a rule-based setback control and a model predictive control (MPC)—with the BCCI and evaluate their performance throughout the summer of 2024. This study reveals that the main bottlenecks in establishing BCCIs for SMBs are the high initial cost (∼$56/m 2 ) and communication delays (up to 9 min). Additionally, the assessment of the two ASCSs indicates that the majority (∼75 %) of the engineering effort required for implementation and configuration is spent on model identification, debugging, and tuning.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

A fast and accurate modeling approach for water and steam thermodynamics with practical applications in district heating system simulation

In U.S. district heating (DH) systems, steam is the most common heat transport medium. Industry demand for new advanced modeling capabilities of complete steam DH systems is increasing; however, the existing models for water/steam thermodynamics are too slow for large system simulations because of computationally expensive algebraic loops that require the solution to nonlinear systems of equations. For practical applications, this work presents a novel split-medium approach that implements numerically efficient liquid water models alongside various water/steam models, breaking costly algebraic loops by decoupling mass and energy balance equations. New component models for steam DH systems are also presented. We implemented the models in the equation based Modelica language and evaluated accuracy and computing speed across multiple scales: from fundamental thermodynamic properties to complete districts featuring 10 to 200 buildings. Compared to district models with the IF97 water/steam model and equipment models from the Modelica Standard Library, the new implementation improves the scaling rate for large districts from cubic to quadratic with negligible compromise to accuracy. Additionally, for an annual simulation with 180 buildings, this translates to a computing time reduction from 33 to 1-1.5 h. These results are critically important for industry practitioners to simulate steam DH systems at large scales.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

City-level impacts of building tune-ups: Findings from Seattle's building tune-ups program

Many U.S. cities are implementing policies to reduce greenhouse gas (GHG) emissions of their buildings. These range from building energy benchmarking and disclosure to building performance standards (BPS) that require buildings to meet specific targets of energy use or emissions. The City of Seattle adopted a climate action plan in 2013 that set a goal of zero net GHG emissions in the road transportation, buildings, and waste sectors by 2050, with a number of near and long term actions. Seattle implemented mandatory building tune-ups in 2016, applying to commercial buildings larger than 50,000 sqft. Building tune-ups1 involve assessment and implementation of operational and maintenance (O + M) improvements to achieve energy and water efficiency, such as changes to thermostat set points or adjusting lighting or irrigation schedules. Seattle's tune-ups program covered 27 such improvements in HVAC, lighting, domestic hot water, and envelope systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2027

Nodal capacity expansion planning with flexible large-scale load siting

We propose explicitly incorporating large-scale load siting into a stochastic nodal power system capacity expansion planning model that concurrently co-optimizes generation, transmission, and storage expansion. The potential operational flexibility of some of these large loads is also taken into account by considering them as consisting of a set of tranches with different reliability requirements, which are modeled as a constraint on expected served energy across operational scenarios. We implement our model as a two-stage stochastic mixed-integer optimization problem with cross-scenario expectation constraints. To overcome the challenge of scalability, we build upon existing work to implement this model on a high performance computing platform and exploit scenario parallelization using an augmented Progressive Hedging Algorithm. The algorithm is implemented using the bounding features of mpisppy, which have shown to provide satisfactory provable optimality gaps despite the absence of theoretical guarantees of convergence. We test our approach and assess the value of this proactive planning framework on total system cost and reliability metrics using realistic testcases geographically assigned to San Diego and South Carolina, with datacenter and direct air capture facilities as large loads.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Go slow to go fast? A review of the impacts of permitting on large-scale solar project development

State and local permitting challenges could impede the ability of large-scale solar (LSS) to meet growing electricity demand in the United States. Here, we review research that explores LSS permitting and its impacts on the pace and scale of LSS project development. Research on LSS permitting is relatively scarce, such that we support our review with research in the context of wind permitting, where appropriate. Further, few studies attempt to rigorously quantify the effects of permitting on the pace and scale of LSS project development. The available evidence allows us to identify various hypotheses and identify gaps for further research. Our review suggests that differences in permitting policies, regulations, and ordinances explain relatively little variation in LSS permitting and development outcomes across jurisdictions, except where jurisdictions implement rules designed to impede LSS. The evidence suggests LSS permitting challenges largely accrue during the implementation of permitting processes. Recent research suggests that community opposition to project development is a key driver of LSS permitting challenges, given that project opponents often use permitting processes to translate opposition into legal action. We call on future researchers to more concretely describe the LSS permitting challenge and to identify the specific actors responsible for implementing solutions.

Community opposition↗
DOE OSTI2022

Demonstration of RELAP5-3D for transient analysis of a dual coolant lead lithium fusion blanket concept

Our research aims to build a multiphysics framework for transient analysis of the Dual- Coolant Lead-Lithium (DCLL) blanket design of the proposed Fusion Nuclear Science Facility (FNSF) for all Onticipated Operational Occurrences (AOOs) using Reduced Order Models (ROMs). This paper presents an assessment of using RELAP5-3D for transient thermal-hydraulic blanket analysis of the DCLL design. Preliminary validation studies for the RELAP5-3D properties were conducted using a representative vertical flow loop model and a comparison study using experimental data from the thermal convection corrosion loops at ORNL. These studies inspected the thermal hydraulic response of Lead-Lithium Eutectic (PbLi) systems within the RELAP5-3D code; mainly pressure drop and heat transfer, as compared to analytical calculations and experimental data. The Magnetohydrodynamic (MHD) pressure drop effect for liquid metal under the effects of a uniform magnetic field was implemented into our models using an equivalent forms loss coefficient based on the phasic momentum equations. We verified our implementation of MHD pressure drop within RELAP5-3D using a systematic approach consisting of two studies. Models were developed to perform comparison studies based on experimental data from the Argonne Liquid metal EXperiment (ALEX) facility and the Magnetohydrodynamic PbLi Experiment (MaPLE) facility. Using the developed validation basis, we built a simplified DCLL channel model within RELAP5-3D, inspected MHD pressure drop and performed transient temperature analysis. The DCLL channel model used radial heating data from the International Thermonuclear Experimental Reactor (ITER) which was normalized to the FNSF design. Here, a representative startup transient was performed using a relative power curve from the initial ramp up of a long-term power pulse of the Tore Supra reactor. This study demonstrated that RELAP5-3D ROMs are capable of accurate transient analysis of the DCLL blanket. We have since utilized MCNP to develop heating profiles based on the FNSF design and are working toward implementing them into our DCLL model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An exploration of online-simulation-driven portfolio scheduling in Workflow Management Systems

Workflow Management Systems used to automate the execution of scientific workflow applications on parallel and distributed computing platforms must make scheduling decisions at runtime. A large number of workflow scheduling algorithms have been proposed in the literature, but often these algorithms are evaluated based on simplifying assumptions that may not hold in practice. Furthermore, published algorithm evaluation and/or comparison results are necessarily only for a subset of all possible scenarios, and thus may not include scenarios relevant to particular use-cases. Consequently, it is difficult for Workflow Management Systems (WMSs) developers to decide which scheduling algorithm should be implemented. To obviate this difficulty, one possible approach is to implement a portfolio of scheduling algorithms and select the most effective algorithm at runtime. One method for performing this selection is to run an online simulation for each algorithm in the portfolio. The algorithm that leads to the best performance, in simulation, is selected for future use. The above simulation-driven portfolio scheduling (SDPS) approach has been proposed in a few parallel and distributed computing contexts. The main objective of this work is to evaluate the feasibility and potential merit of SDPS if implemented in WMSs. Here we perform this evaluation using simulated WMS executions, where the simulations are instantiated from real-world platform and workflow configurations. Our main finding is that SDPS is on par with or outperforms an approach in which a single algorithm is used, where this algorithm is the one that performs best on average across all our experimental scenarios. Furthermore, we find that SDPS remains an attractive proposition even in the presence of high levels of simulation error and for simulators with relatively low levels of sophistication. In many of our experimental scenarios we find that mitigating simulation error at runtime can further improve performance. Finally, we show that simulation overhead can be made sufficiently low for SDPS to be feasible in practice.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Process modeling, techno-economic assessment, and life cycle assessment of the electrochemical reduction of CO 2 : a review

The electrochemical reduction of CO 2 has emerged as a promising alternative to traditional fossil-based technologies for the synthesis of chemicals. Its industrial implementation could lead to a reduction in the carbon footprint of chemicals and the mitigation of climate change impacts caused by hard-to-decarbonize industrial applications, among other benefits. However, the current low technology readiness levels of such emerging technologies make it hard to predict their performance at industrial scales. During the past few years, researchers have developed diverse techniques to model and assess the electrochemical reduction of CO 2 toward its industrial implementation. The aim of this literature review is to provide a comprehensive overview of techno-economic and life cycle assessment methods and pave the way for future assessment approaches. First, we identify which modeling approaches have been conducted to extend analysis to the production scale. Next, we explore the metrics used to evaluate such systems, regarding technical, environmental, and economic aspects. Finally, we assess the challenges and research opportunities for the industrial implementation of CO 2 reduction via electrolysis.

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

Massively parallel transport sweeps on meshes with cyclic dependencies

When solving the first-order form of the linear Boltzmann equation, a common misconception is that the matrix-free computational method of “sweeping the mesh”, used in conjunction with the Discrete Ordinates method, is too complex or does not scale well enough to be implemented in modern high performance computing codes. This has led to considerable efforts in the development of matrix-based methods that are computationally expensive and is partly driven by the requirements placed on modern spatial discretizations. In particular, modern transport codes are required to support higher order elements, a concept that invariably adds a lot of complexity to sweeps because of the introduction of cyclic dependencies with curved mesh cells. In this article we will present a comprehensive implementation of sweeping, to a piecewise-linear DFEM spatial discretization with particular focus on handling cyclic dependencies and possible extensions to higher order spatial discretizations. We find that these methods are implemented in a new C++ simulation framework called Chi-Tech (). We present some typical simulation results with some performance aspects that one can expect during real world simulations, we also present a scaling study to >100k processes where Chi-Tech maintains greater than 80% efficiency solving a total of 87.7 trillion angular flux unknowns for a 116 group simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A deterministic verification strategy for electrostatic particle-in-cell algorithms in arbitrary spatial dimensions using the method of manufactured solutions

As simulations of kinetic plasmas continue to increase in scope and complexity, a rigorous and straightforward method for verifying particle-in-cell (PIC) implementations is necessary to ensure their correctness. Here, in this paper, we present a deterministic method for the rigorous verification of multidimensional, multispecies, electrostatic particle-in-cell codes based on the method of manufactured solutions. Specifically, we prove that rigorous verification is possible through the exclusive examination of errors of grid quantities (i.e., moments and/or fields), allowing for a very light-weight and non-intrusive implementation in existing PIC codes. This is a marked improvement over earlier PIC verification approaches (only demonstrated with one species in 1D-1V), which rely on the comparison of cumulative distribution functions, and are invasive on the code base, introduce additional stochastic noise, are significantly more computationally expensive, and lack rigorous convergence properties. Interestingly, we show that different grid quantities feature different rates of convergence with the number of particles and mesh size, impacting the verification process. These theoretical results are confirmed numerically with a multi-species 2D-2V particle-in-cell code, including the ability of the method to detect order reduction due to an incorrect implementation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. Here this algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Time-Domain Vortex Induced Vibration Modeling of Reference Dynamic Power Cable for the Gulf of Maine

Vortex-induced vibration (VIV) is a phenomenon known to decrease the fatigue life of dynamic power cables used in floating offshore wind systems through increased bending loading cycles. A variety of modeling solutions have been proposed to study this, but none enable fully-coupled simulations leveraging open-source tools. To address this, a time-domain VIV model has been successfully implemented into the open-source mooring dynamics model MoorDyn. The semi-empirical VIV model describes a lift force acting on segments of a flexible cylinder, making it well suited for implementation into MoorDyn. The new capability successfully predicted the frequencies and magnitudes of strain in both steady and oscillating flows when compared to the original validation results and matched peak spectral response when compared to experimental results, verifying successful implementation. MoorDyn with VIV was then leveraged to simulate a 15 MW floating turbine with a dynamic power cable in the Gulf of Maine using OpenFAST for six return periods of combined wind, wave, and current conditions. Increased curvature and tension fluctuations in higher flow speeds were observed when simulating VIV. Maximum tensions and curvatures also increased, with the 500 year conditions violating the curvature factor of safety of 2.0. These results highlight the importance of considering VIV when designing dynamic power cables for the Gulf of Maine, and indicate that cost effective mitigation strategies should be explored for the area. They also demonstrate the utility of this new modeling capability, which provides the first open-source tool for fully coupled time-domain floating offshore wind simulations with dynamic power cable VIV.

17 WIND ENERGY↗
DOE OSTI2022

Discrete dislocation dynamics for crystal RVEs. Part 1: Periodic network kinematics

A novel implementation of the dislocation flux boundary condition in discrete dislocation dynamics is presented. The continuity of the individual dislocation loops in a periodic representative crystal volume (RVE) is enforced across the boundary of the RVE with the help of a dual topological description for representing dislocation line kinematics in two equivalent spaces representing the deforming crystal, the RVE and the unbounded crystal spaces. The former describes the motion of the dislocations in the simulated crystal RVE whereas the latter represents the motion of dislocations in an infinite space containing all replicas of the RVE. A mapping between the two spaces forms the basis of the implementation of flux boundary condition. The implementation details are discussed in the context of statistical homogeneity of bulk crystals undergoing macroscopically homogeneous plastic deformation. In this case, the boundary nodes associated with dislocation segments bear no relevance in the motion of the dislocations. Finally, some test cases are presented and discussed to establish the proposed approach.

42 ENGINEERING↗
DOE OSTI2023

Complementing a continuum thermodynamic approach to constitutive modeling with symbolic regression

Here, an interpretable machine learning method, physics-informed genetic programming-based symbolic regression (P-GPSR), is integrated into a continuum thermodynamic approach to developing constitutive models. The proposed strategy for combining a thermodynamic analysis with P-GPSR is demonstrated by generating a yield function for an idealized material with voids, i.e., the Gurson yield function. First, a thermodynamic-based analysis is used to derive model requirements that are exploited in a custom P-GPSR implementation as fitness criteria or are strongly enforced in the solution. The P-GPSR implementation improved accuracy, generalizability, and training time compared to the same GPSR code without physics-informed fitness criteria. The yield function generated through the P-GPSR framework is in the form of a composite function that describes a class of materials and is characteristically more interpretable than GPSR-derived equations. The physical significance of the input functions learned by P-GPSR within the composite function is acquired from the thermodynamic analysis. Fundamental explanations of why the implemented P-GPSR capabilities improve results over a conventional GPSR algorithm are provided.

42 ENGINEERING↗