Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “system simulation”

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 OSTI2025

Privacy Preservation from High-Performance Computing to Autonomous Science [Industrial and Governmental Activities]

High-Performance Computing (HPC) and Leadership-Class Supercomputing are driving forces behind scientific advancements, enabling researchers to tackle complex challenges in physics, chemistry, biology, and engineering. These systems power vast simulations and data analyses, fueling discoveries in fields ranging from materials science to climate modeling. However, their use often involves processing sensitive data—such as proprietary industry simulations, biomedical records, and national security computations—posing significant privacy concerns. In conclusion, this issue is amplified in collaborative environments like Department of Energy (DOE) user facilities, where HPC resources are shared across institutions to foster innovation.

Kotevska, Olivera [Oak Ridge National Laboratory (↗
DOE OSTI2020

Virtualized Logical Qubits: A 2.5D Architecture for Error-Corrected Quantum Computing

Current, near-term quantum devices have shown great progress in the last several years culminating recently with a demonstration of quantum supremacy. In the medium-term, however, quantum machines will need to transition to greater reliability through error correction, likely through promising techniques like surface codes which are well suited for near-term devices with limited qubit connectivity. We discover quantum memory, particularly resonant cavities with transmon qubits arranged in a 2.5D architecture, can efficiently implement surface codes with substantial hardware savings and performance/fidelity gains. Specifically, we virtualize logical qubits by storing them in layers of qubit memories connected to each transmon. Surprisingly, distributing each logical qubit across many memories has a minimal impact on fault tolerance and results in substantially more efficient operations. Our design permits fast transversal application of CNOT operations between logical qubits sharing the same physical address (same set of cavities) which are 6x faster than standard lattice surgery CNOTs. We develop a novel embedding which saves approximately 10x in transmons with another 2x savings from an additional optimization for compactness. Although qubit virtualization pays a 10x penalty in serialization, advantages in the transversal CNOT and in area efficiency result in fault-tolerance and performance comparable to conventional 2D transmon-only architectures. Our simulations show our system can achieve fault tolerance comparable to conventional two-dimensional grids while saving substantial hardware. Furthermore, our architecture can produce magic states at 1.22x the baseline rate given a fixed number of transmon qubits. Here, this is a critical benchmark for future fault-tolerant quantum computers as magic states are essential and machines will spend the majority of their resources continuously producing them. This architecture substantially reduces the hardware requirements for fault-tolerant quantum computing and puts within reach a proof-of-concept experimental demonstration of around 10 logical qubits, requiring only 11 transmons and 9 attached cavities in total.

quantum computing↗
DOE OSTI2020

A Model-Free Frequency Control Approach for Diesel-Wind Powered Microgrids

Islanded microgrids usually consist of diesel generators and renewable energy sources (RES) to reduce the operating cost. Such microgrids have shown values for powering remote locations, but introduced unique challenges for frequency control of the grid due to the variability of renewable energy as well as the decoupled design of converter interface in RES. So, supportive control in RES becomes mandatory in such networks. This paper proposes a new dynamic control strategy, based on model-free control (MFC) approach, to support the primary frequency control of such islanded microgrids. The practical values of MFC have been discussed in various domains due to its control capability without modeling procedure and its efficient computation. In particular, we consider the diesel-wind system and the proposed MFC strategy is utilized as an online controller, which computes supplementary control signals for the rotor-side converter of the available wind turbine generators (WTGs). The calculated input signals allow WTGs to respond to frequency variations and improve the frequency response of the system. The controller is implemented and verified using the modified IEEE 33-bus full nonlinear model of the three-phase diesel-wind system in Simulink. Simulation results show the enhanced frequency response from the MFC strategy.

Park, Byungkwon↗
DOE OSTI2021

Trap-Assisted Dopant Compensation Prevents Shunting in poly-Si Passivating Interdigitated Back Contact Silicon Solar Cells

Interdigitated back contact (IBC) solar cells achieve the highest efficiencies of single-junction architectures, but complicated patterning of the rear fingers and spreading of dopants during processing inhibit their mainstream adoption due to concerns of shunting between the IBC fingers. One method of simplifying patterning at the rear is by using contact masks combined with plasma-enhanced chemical vapor deposition (PECVD) or ion implantation. However, the intrinsic isolation region becomes contaminated during high-temperature annealing by lateral diffusion of dopants and during masked PECVD by spreading of dopant radicals through region between the mask and the substrate. Despite this contamination, we show through scanning spreading resistance microscopy and Kelvin probe force microscopy that a ~20 µm wide compensating region exists with high enough resistivity to prevent shunting. We model this p-i-n poly-Si system using two simulation models: a simple resistor model considering only the capture of charge carriers by trap defects in poly-Si to reduce the conductivity, and a more refined 1-dimensional finite element model using Poisson’s equation, drift-diffusion equations, and recombination of carriers. Using this model, we show that high defect density significantly decreases the current across the region between the p- and n-type fingers, preventing shunting.

atom probe↗
DOE OSTI2020

Four-Switch Buck-Boost Converter Based on Model Predictive Control with Smooth Mode Transition Capability

Four-switch buck-boost converter supports both voltage step-up and step-down functionalities, but it suffers from mode transfer challenge that would need reliable mode detection when designed to operate in multi different modes. In this paper, a novel control method based on model predictive current control is proposed for such converter with inherent smooth mode transfer capability without extra design on mode detection and transfer scheme. Modulator and mode detection are replaced by an optimization process through cost function. This largely simplifies the design and makes it easily implemented. Therefore, seamless transfer between buck and boost modes is achieved with many other system level benefits. Simulation and experimental results are provided to verify the effectiveness of proposed method for the four-switch buck-boost converter.

42 ENGINEERING↗
DOE OSTI2024

Buyers Collusion in Incentivized Forwarding Networks: A Multi-Agent Reinforcement Learning Study

We present the issue of monetarily incentivized forwarding in a multi-hop mesh network architecture from an economic perspective. It is anticipated that credit-incentivized forwarding and relaying will be a simple method of exchanging transmission power and spectrum for connectivity. However, gateways and forwarding nodes, like any other free market, may create an oligopolistic market for the users they serve. In this study, a coalition scheme between buyers aims to address price control by gateways or nodes closer to gateways. In a Stackelberg competition game, buyer agents (users) and sellers (gateways) make decisions using reinforcement learning (RL), with decentralized Deep Q-Networks to buy and sell forwarding resources. We allow communication links between the buyers with a limited messaging space, without defining a collusion mechanism. The idea is to demonstrate that through messaging, and RL tacit collusion can emerge between agents in a decentralized setup. The multi-agent reinforcement learning (MARL) system is presented and analyzed from a machine-learning perspective. Moreover, MARL dynamics are discussed via mean field analysis to better understand divergence causes and make implementation recommendations for such systems. Finally, the simulation results show the results of coordination among the users.

42 ENGINEERING↗
DOE OSTI2026

Specific ion effects on ion transport in charged polymer membranes

Selective separation of like-charged ions is a central challenge in applications such as critical mineral recovery. Electrochemical membrane–based separations offer promising pathways to address this need, but a limited fundamental understanding of ion transport in charged polymer membranes hampers the development of highly selective materials. This study elucidates the role of specific ion effects (SIEs) in ion transport through charged polymer membranes by integrating experimental measurements of ion mobility and in situ ion/ion and ion/water interactions in model charged polymer membranes with results from molecular dynamics simulations on analogous systems. We demonstrate that solvent-mediated ion interactions drive pronounced SIEs within the studied membranes, with the ion softness, i.e., the malleability of ion hydration shells, from hard/soft acid/base (HSAB) theory emerging as a key predictor of transport properties. HSAB theory also explains the observed mechanism of solvent-mediated ion interactions. Our findings offer a mechanistic framework for designing membranes with tailored ion selectivity, potentially enabling efficient separations of chemically similar ions.

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

GPU acceleration of Swendsen–Wang dynamics

When simulating a lattice system near its critical temperature, local algorithms for modeling the system’s evolution can introduce very large autocorrelation times into sampled data. Here, this critical slowing down places restrictions on the analysis that can be completed in a timely manner of the behavior of systems around the critical point. Because it is often desirable to study such systems around this point, a new algorithm must be introduced. Therefore, we turn to cluster algorithms, such as the Swendsen–Wang algorithm and the Wolff clustering algorithm. They incorporate global updates which generate new lattice configurations with little correlation to previous states, even near the critical point. We look to accelerate the rate at which these algorithm are capable of running by implementing and benchmarking a parallel implementation of each algorithm designed to run on GPUs under NVIDIA’s CUDA framework. A 17 and 90 fold increase in the computational rate was, respectively, experienced when measured against the equivalent algorithm implemented in serial code.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Automatic Differentiation of C++ Codes on Emerging Manycore Architectures with Sacado

Automatic differentiation (AD) is a well-known technique for evaluating analytic derivatives of calculations implemented on a computer, with numerous software tools available for incorporating AD technology into complex applications. However, a growing challenge for AD is the efficient differentiation of parallel computations implemented on emerging manycore computing architectures such as multicore CPUs, GPUs, and accelerators as these devices become more pervasive. In this work, we explore forward mode, operator overloading-based differentiation of C++ codes on these architectures using the widely available Sacado AD software package. In particular, we leverage Kokkos, a C++ tool providing APIs for implementing parallel computations that is portable to a wide variety of emerging architectures. Here we describe the challenges that arise when differentiating code for these architectures using Kokkos, and two approaches for overcoming them that ensure optimal memory access patterns as well as expose additional dimensions of fine-grained parallelism in the derivative calculation. We describe the results of several computational experiments that demonstrate the performance of the approach on a few contemporary CPU and GPU architectures. We then conclude with applications of these techniques to the simulation of discretized systems of partial differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Neurons to Algorithms

A framework for modeling and simulating biological neural systems

Warrender, Christina↗
DOE OSTI2020

ddcMD.os

ddcMD is a general purpose molecular dynamics (MD) code that supports MPI parallelism. MD codes are used for simulation of particles systems and capture all the many-body effects of the underlying particle potential that defines the physical system. Though MD can be used to model systems from the subatomic to astrological length scales ddcMD is mainly focused on the atomic scale length scale, In this release of ddcMD the support will be mainly for systems using the coarse-grain Martini potential, a particle potential for biological systems.

Glosli, JamesN↗
DOE OSTI2023

PFLOTRAN 5

PFLOTRAN leverages massively parallel, high performance computing to simulate large-scale non-isothermal multiphase flow, multicomponent reactive transport and electrical resistivity tomography (ERT) problems in the subsurface environment. Researchers have employed PFLOTRAN to simulate these Earth system processes on leadership class supercomputers for over two decades. The code is designed to predict the future estate of environmental systems and better inform stakeholders in the regulatory decision making process (e.g., fate of contaminants, long-term stewardship for nuclear waste, impact of climate change, etc.). A diverse team of scientists oversees PFLOTRAN development and maintenance under an open-source licensing agreement and manages contributions from an international community of researchers.

Hammond, Glenn↗
DOE OSTI2024

pnnl/HNSRigidBodySimulator

HNS Rigid Body Simulatoris a tool for performing rigid-body dynamics simulations of nanoscale systems (e.g., macromolecular assemblies) in Python.

Marina, Sofia↗
DOE OSTI2026

The U.S. DOE ARM User Facility Establishes a New Site for Studies of Land–Aerosol–Cloud Interactions in the Southeastern United States

The U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility has established a new site in the Bankhead National Forest (BNF) in northern Alabama that will gather data on how clouds, the land surface, and aerosols interact at a hierarchy of scales important to understanding and simulating the Earth system. Starting its operations in October 2024, the BNF site provides a multi-year opportunity for scientists to unravel complex land-atmosphere interactions. A suite of ground-based sensors, elevated tower-based instrumentation, and aerial facilities will enable scientists to investigate those interactions from within the canopy to the clouds. The Southeastern United States was recommended by the DOE ARM and its collaborators in the broader community as an important region to address their common scientific questions given the region’s abundant surface-forced convective clouds and mesoscale convective systems that pose ongoing challenges in Earth system models. The region is also home to significant terrain complexity and land-use heterogeneity that will unleash new understanding of anthropogenic and biogenic aerosol processes, boundary-layer aerosol-cloud interactions, and the interactions between the terrestrial ecosystem and coupled aerosol-cloud-radiation processes.

Aerosols↗
DOE OSTI2025

Forest Breeze–Cold Pool Interactions Drive Convective Organization over Heterogeneous Vegetation

Heterogeneous landscapes can influence the development of convection through the generation of thermally driven mesoscale circulations. To assess the impacts of these circulations and their interaction with sea breezes, we simulated convection in an idealized coastal environment using the Regional Atmospheric Modeling System (RAMS). We compared simulations with striped patterns of surface vegetation to those of uniform vegetation to identify the importance of vegetation heterogeneity in impacting convective development. Under dry soil conditions representative of those during the Tracking Aerosol Convection Interactions Experiment (TRACER) and Experiment of Sea Breeze Convection, Aerosols, Precipitation, and Environment (ESCAPE) campaigns in June 2022, we found that these vegetation-induced circulations, referred to in the literature as “forest breezes,” are more important than the sea breeze in determining the location of convection initiation. Convection and precipitation are also found to be favored over forests and suppressed over pasture and suburban landscapes as a result of greater surface sensible heat flux over the forest. Our findings also indicate that forest breezes are important for initiating convection along the boundaries of the forest, but that cold pools may play a key role in propagating the forest breezes toward the center of the forest stripe. In our simulations, the collisions of these breezes in the center of the forest stripe lead to uplift and strong convection there; however, a different width of the forest stripe would alter when the forest breezes collide or whether they collide at all. The presence of these cold pools may therefore impact the “ideal stripe width,” the width of each vegetation stripe which maximizes domain-wide precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A numerical approach for detecting switch-like bistability in mass action chemical reaction networks with conservation laws

Theoretical analysis of signaling pathways can provide a substantial amount of insight into their function. One particular area of research considers signaling pathways capable of assuming two or more stable states given the same amount of signaling ligand. This phenomenon of bistability can give rise to switch-like behavior, a mechanism that governs cellular decision making. Investigation of whether or not a signaling pathway can confer bistability and switch-like behavior, without knowledge of specific kinetic rate constant values, is a mathematically challenging problem. Recently a technique based on optimization has been introduced, which is capable of finding example parameter values that confer switch-like behavior for a given pathway. Although this approach has made it possible to analyze moderately sized pathways, it is limited to reaction networks that presume a uniterminal structure. It is this limited structure we address by developing a general technique that applies to any mass action reaction network with conservation laws. In this paper we developed a generalized method for detecting switch-like bistable behavior in any mass action reaction network with conservation laws. The method involves (1) construction of a constrained optimization problem using the determinant of the Jacobian of the underlying rate equations, (2) minimization of the objective function to search for conditions resulting in a zero eigenvalue, (3) computation of a confidence level that describes if the global minimum has been found and (4) evaluation of optimization values, using either numerical continuation or directly simulating the ODE system, to verify that a bistability region exists. The generalized method has been tested on three motifs known to be capable of bistability. We have developed a variation of an optimization-based method for the discovery of bistability, which is not limited to uniterminal chemical reaction networks. Successful completion of the method provides an S-shaped bifurcation diagram, which indicates that the network acts as a bistable switch for the given optimization parameters.

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

A fast particle-based approach for calibrating a 3-D model of the Antarctic ice sheet

We consider the scientifically challenging and policy-relevant task of understanding the past and projecting the future dynamics of the Antarctic ice sheet. The Antarctic ice sheet has shown a highly nonlinear threshold response to past climate forcings. Triggering such a threshold response through anthropogenic greenhouse gas emissions would drive drastic and potentially fast sea level rise with important implications for coastal flood risks. Previous studies have combined information from ice sheet models and observations to calibrate model parameters. These studies have broken important new ground but have either adopted simple ice sheet models or have limited the number of parameters to allow for the use of more complex models. These limitations are largely due to the computational challenges posed by calibration as models become more computationally intensive or when the number of parameters increases. Here, we propose a method to alleviate this problem: a fast sequential Monte Carlo method that takes advantage of the massive parallelization afforded by modern high-performance computing systems. We use simulated examples to demonstrate how our sample-based approach provides accurate approximations to the posterior distributions of the calibrated parameters. The drastic reduction in computational times enables us to provide new insights into important scientific questions, for example, the impact of Pliocene era data and prior parameter information on sea level projections. These studies would be computationally prohibitive with other computational approaches for calibration such as Markov chain Monte Carlo or emulation-based methods. We also find considerable differences in the distributions of sea level projections when we account for a larger number of uncertain parameters. For example, based on the same ice sheet model and data set, the 99th percentile of the Antarctic ice sheet contribution to sea level rise in 2300 increases from 6.5 m to 13.1 m when we increase the number of calibrated parameters from three to 11. With previous calibration methods, it would be challenging to go beyond five parameters. Here, this work provides an important next step toward improving the uncertainty quantification of complex, computationally intensive and decision-relevant models.

54 ENVIRONMENTAL SCIENCES↗