Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed 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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Development of an Encoding Method on an Co-simulation Platform for Mitigating the Impact of Unreliable Communication

This report presents a hardware-in-the-loop (HIL) based modeling approach for simulating impacts of unreliable communication on the performance of centralized volt-var control and for developing an encoding method to mitigate the impacts. First, an asynchronous real-time HIL simulation platform is introduced to enable multi-rate co-simulation of a distribution system with many inverter-based distributed energy resources (DERs). The distribution system is modeled by milliseconds phasor-based models and the DERs are modeled by micro-seconds power electronic models. Communication connections between a centralized volt-var controller (modeled externally to the HIL testbed) and smart inverters are built by implementing Modbus links and the Long Term Evolution network. On this co-simulation platform, an enhanced, augmented Lagrangian multiplier based encoded data recovery (EALM-EDR) algorithm for mitigating the impact of unreliable communication is developed and validated. Simulation results demonstrate the efficacy of using the HIL-based co-simulation platform as a power grid digital twin for developing algorithms that coordinate a large number of heterogeneous control systems through wired and wireless communication links.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Simulating Magnetized Electron Cooling for EIC with JSPEC

We present a possible electron cooling configuration for the proposed Electron Ion Collider (EIC) facility, developed using a Nelder-Mead Simplex optimization procedure built into JSPEC, an electron cooling code developed at Jefferson Lab. We show the time evolution of the emittance of the ion beam in the presence of this cooler evaluated assuming the ion distribution remains Gaussian. We also show that bi-gaussian distributions emerge in simulations of ion macro-particles. We show how intra-beam scattering can be treated with a core-tail model in simulations of ion macro-particles. The Sirepo/JSPEC* and Sirepo/Jupyter** apps will be presented, with instructions enabling the community to reproduce our simulations.

Coleman, S. J.↗
DOE OSTI2026

MixPI: Mixed-time slicing path integral software for quantized molecular dynamics simulations

We introduce the MixPI software to implement path integral molecular dynamics (PIMD) simulations for the study of condensed phase systems where nuclear quantum effects (NQEs) are important. In contrast to existing PIMD simulation software, MixPI enables the implementation of mixed quantum–classical path integral simulations where only a subset of system degrees of freedom (dofs) are treated quantum mechanically in an extended phase space while the remaining dofs are described classically. We expect this software to be particularly useful for simulations of electron and proton transfer in condensed phase systems, as well as for the study of biological and material systems where only a handful of dofs contribute significantly to the observed NQEs. We demonstrate the use of MixPI in two different systems. The first is a simple water model where we implement a set of mixed quantum–classical simulations to compute average energy and radial distribution functions. We use these simulations to benchmark the effectiveness of MixPI and to demonstrate how it enables systematic investigation into the origin of observed NQEs. We then compute radial distribution functions for a system where MixPI is essential: a solvated metal (M 2+ ) cation described using an explicit quantized electron localized on an M 3+ ion in water.

chemical physics↗
DOE OSTI2015

One-dimensional particle-in-cell simulations of electrostatic Bernstein waves in plasmas with kappa velocity distributions

Electrostatic Bernstein waves that propagate exactly perpendicularly to a static magnetic field in an electron-ion plasma are investigated using one-and-two-halves dimensional particle-in-cell simulations. An ion-to-electron mass ratio of m{sub i}/m{sub e} = 100 is used, allowing sufficient separation of the electron and ion time scales while still accounting for the ion dynamics without resorting to exceptionally long simulation run times. As a consequence of the mass ratio used, both the high frequency electron Bernstein wave and the lower frequency ion Bernstein wave are resolved within a single simulation run. The simulations presented here use isotropic three-dimensional kappa velocity distributions as well as the widely used Maxwellian velocity distribution, and the results from using each of these velocity distributions are analysed and compared. The behaviour of the Bernstein waves is found to be significantly dependent on the spectral index, κ, of the kappa distribution in all frequency domains of the Bernstein waves. In both the Maxwellian and kappa cases, spectral analysis of the electric field (wave) intensities, as a function of ω and k, show very good agreement between the simulation results and the linear dispersion relation for Bernstein waves. This agreement serves to validate the simulation techniques used, as well as the theory of Bernstein waves in plasmas with a kappa velocity distribution. The intensity of the field fluctuations in the simulations containing an abundance of superthermal particles, i.e., where the plasma has a kappa velocity distribution with a low kappa index, is slightly higher compared to the simulations of plasmas with higher kappa values. The plasmas with low kappa values also exhibit a broader region in frequency space of high intensity field fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Influence of microwave driver coupling design on plasma density at Testbench for Ion sources Plasma Studies, a 2.45 GHz Electron Cyclotron Resonance Plasma Reactor

A comparative study of two microwave driver systems (preliminary and optimized) for a 2.45 GHz hydrogen Electron Cyclotron Resonance plasma generator has been conducted. The influence on plasma behavior and parameters of stationary electric field distribution in vacuum, i.e., just before breakdown, along all the microwave excitation system is analyzed. 3D simulations of resonant stationary electric field distributions, 2D simulations of external magnetic field mapping, experimental measurements of incoming and reflected power, and electron temperature and density along the plasma chamber axis have been carried out. By using these tools, an optimized set of plasma chamber and microwave coupler has been designed paying special attention to the optimization of stationary electric field value in the center of the plasma chamber. This system shows a strong stability on plasma behavior allowing a wider range of operational parameters and even sustaining low density plasma formation without external magnetic field. In addition, the optimized system shows the capability to produce values of plasma density four times higher than the preliminary as a consequence of a deeper penetration of the magnetic resonance surface in relative high electric field zone by keeping plasma stability. The increment of the amount of resonance surface embedded in the plasma under high electric field is suggested as a key factor.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2011

Patient-specific radiation dose and cancer risk estimation in CT: Part I. Development and validation of a Monte Carlo program

Purpose: Radiation-dose awareness and optimization in CT can greatly benefit from a dose-reporting system that provides dose and risk estimates specific to each patient and each CT examination. As the first step toward patient-specific dose and risk estimation, this article aimed to develop a method for accurately assessing radiation dose from CT examinations. Methods: A Monte Carlo program was developed to model a CT system (LightSpeed VCT, GE Healthcare). The geometry of the system, the energy spectra of the x-ray source, the three-dimensional geometry of the bowtie filters, and the trajectories of source motions during axial and helical scans were explicitly modeled. To validate the accuracy of the program, a cylindrical phantom was built to enable dose measurements at seven different radial distances from its central axis. Simulated radial dose distributions in the cylindrical phantom were validated against ion chamber measurements for single axial scans at all combinations of tube potential and bowtie filter settings. The accuracy of the program was further validated using two anthropomorphic phantoms (a pediatric one-year-old phantom and an adult female phantom). Computer models of the two phantoms were created based on their CT data and were voxelized for input into the Monte Carlo program. Simulated dose at various organ locations was compared against measurements made with thermoluminescent dosimetry chips for both single axial and helical scans. Results: For the cylindrical phantom, simulations differed from measurements by -4.8% to 2.2%. For the two anthropomorphic phantoms, the discrepancies between simulations and measurements ranged between (-8.1%, 8.1%) and (-17.2%, 13.0%) for the single axial scans and the helical scans, respectively. Conclusions: The authors developed an accurate Monte Carlo program for assessing radiation dose from CT examinations. When combined with computer models of actual patients, the program can provide accurate dose estimates for specific patients.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Primary Frequency Response with Aggregated DERs

Power networks have to withstand a variety of disturbances that affect system frequency, and the problem is compounded with the increasing integration of intermittent renewable generation. Following a large-signal generation or load disturbance, system frequency is arrested leveraging primary frequency control provided by governor action in synchronous generators. In this work, we propose a framework for distributed energy resources (DERs) deployed in distribution networks to provide (supplemental) primary frequency response. Particularly, we demonstrate how power-frequency droop slopes for individual DERs can be designed so that the distribution feeder presents a guaranteed frequency-regulation characteristic at the feeder head. Furthermore, the droop slopes are engineered such that injections of individual DERs conform to a well-defined fairness objective that does not penalize them for their location on the distribution feeder. Time-domain simulations for an illustrative network composed of a combined transmission network and distribution network with frequency-responsive DERs are provided to validate the approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

IGMS: An Integrated ISO-to-Appliance Scale Grid Modeling System

This paper describes the Integrated Grid Modeling System (IGMS), a novel electric power system modeling platform for integrated transmission-distribution analysis that co-simulates off-the-shelf tools on high performance computing (HPC) platforms to offer unprecedented resolution from ISO markets down to appliances and other end uses. Specifically, the system simultaneously models hundreds or thousands of distribution systems in co-simulation with detailed Independent System Operator (ISO) markets and AGC-level reserve deployment. IGMS uses a new MPI-based hierarchical co-simulation framework to connect existing sub-domain models. Our initial efforts integrate opensource tools for wholesale markets (FESTIV), bulk AC power flow (MATPOWER), and full-featured distribution systems including physics-based end-use and distributed generation models (many instances of GridLAB-D[TM]). The modular IGMS framework enables tool substitution and additions for multi-domain analyses. This paper describes the IGMS tool, characterizes its performance, and demonstrates the impacts of the coupled simulations for analyzing high-penetration solar PV and price responsive load scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Coupled Electro-Thermal Analysis of Permanent Magnet Synchronous Motor for Electric Vehicles

Permanent magnet synchronous motors (PMSM) are extensively used in electric vehicles. However, high internal heat generation and inefficient heat dissipation often limit the operational reliability, and longevity of the PMSM. Therefore, proper quantification of heat generation in electric motor and advanced embedded motor cooling techniques remain topics of immense interest. In order to accurately predict electro-magnetic performance, i.e., efficiency, component-wise heat losses and the corresponding temperature distribution of a jacket cooled machine, this paper presents a two-way iteratively coupled electro-thermal modeling framework. Finite element based software Motor-CAD has been utilized for electro-magnetic performance calculation of BMW i3 PMSM. Finite volume based computational fluid dynamics/heat transfer (CFD/HT) software ANSYS® FLUENT® has been employed to simulate the temperature distribution of the PMSM, using the electro-magnetic losses as heat input. Computed heat losses, stator, winding, rotor, and magnet temperatures are utilized as coupling parameters between the electro-magnetic and thermal models. A conventional one-way coupling algorithm has also been developed and compared to the newly proposed two-way coupling algorithm. Numerical results confirm that at high current density, one-way coupling algorithm significantly over-predicts the motor temperature compared to the two-way algorithm. A comprehensive analysis has been carried out to characterize the influences of current density, speed, and forced convection heat transfer coefficient on the heat losses, overall efficiency, and maximum temperature of the PMSM. Lastly, an efficiency map has been interpreted from the coupled electro-magnetic simulation.

42 ENGINEERING↗
DOE OSTI2016

Joule heating induced thermomigration failure in un-powered microbumps due to thermal crosstalk in 2.5D IC technology

Thermal-crosstalk induced thermomigration failure in un-powered microbumps has been found in 2.5D integrated circuit (IC) circuit. In 2.5D IC, a Si interposer was used between a polymer substrate and a device chip which has transistors. The interposer has no transistors. If transistors are added to the interposer chip, it becomes 3D IC. In our test structure, there are two Si chips placed horizontally on a Si interposer. The vertical connections between the interposer and the Si chips are through microbumps. We powered one daisy chain of the microbumps under one Si chip; however, the un-powered microbumps in the neighboring chip are failed with big holes in the solder layer. We find that Joule heating from the powered microbumps is transferred horizontally to the bottom of the neighboring un-powered microbumps, and creates a large temperature gradient, in the order of 1000 °C/cm, through the un-powered microbumps in the neighboring chip, so the latter failed by thermomigration. In addition, we used synchrotron radiation tomography to compare three sets of microbumps in the test structure: microbumps under electromigration, microbumps under thermomigration, and microbumps under a constant temperature thermal annealing. The results show that the microbumps under thermomigration have the largest damage. Furthermore, simulation of temperature distribution in the test structure supports the finding of thermomigration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Structural and mechanical properties of monolayer amorphous carbon and boron nitride

Amorphous materials exhibit various characteristics that are not featured by crystals and can sometimes be tuned by their degree of disorder (DOD). Here, we report results on the mechanical properties of monolayer amorphous carbon (MAC) and monolayer amorphous boron nitride (maBN) with different DOD. The pertinent structures are obtained by kinetic-Monte-Carlo (kMC) simulations using machine-learning potentials (MLP) with density-functional-theory (DFT)-level accuracy. An intuitive order parameter, namely the areal fraction F x occupied by crystallites within the continuous random network, is proposed to describe the DOD. We find that F x captures the essence of the DOD: Samples with the same F x but different sizes and arrangements of crystallites, obtained using two distinct kMC procedures, have virtually identical radial distributions functions as well as bond-length and bond-angle distributions. Furthermore, by simulating the fracture process with molecular dynamics, we found that the mechanical responses of MAC and maBN before fracture are mainly determined by F x and are insensitive to the sizes and specific arrangements and to some extent the numbers and area distributions of the crystallites. The behavior of cracks in the two materials is analyzed and found to mainly propagate in meandering paths in the CRN region and to be influenced by crystallites in distinct ways that toughen the material. Furthermore, the present results reveal the relation between structure and mechanical properties in amorphous monolayers and may provide a universal toughening strategy for 2D materials.

2-dimensional systems↗
DOE OSTI2025

Deterministic High-Fidelity Neutronics Simulation of Pebble Bed Reactors Using Pebble Tracking Transport

The pebble tracking transport (PTT) algorithm offers a high-fidelity deterministic approach for neutron transport for pebble bed reactors (PBRs). This approach requires the mesh for the active-core region to consist exclusively of tetrahedral elements, where each node in the pebble-packing region represents a pebble centroid. This paper investigates the application of PTT for full-scale PBRs, considering both the isothermal and the temperature-dependent core conditions. Macroscopic cross sections are generated using Serpent 2 full-core eigenvalue simulations where pebbles are grouped into disjoint subsets using machine learning. To minimize the need for individual cross-section sets for each pebble in the core, K-means clustering is used to group pebbles by temperature and neutronic environment parameters. Here, we compare the multiplication factor and power rate distributions between PTT simulations using the Griffin reactor physics software and reference solutions from Serpent 2. Our analysis shows that a full-core, high-fidelity PTT calculation produces accurate results with minimal local (pebblewise) errors. Additionally, timing results indicate that PTT simulations converge rapidly on modern supercomputing platforms.

Griffin↗
DOE OSTI2017

Towards improved parameterization of a macroscale hydrologic model in a discontinuous permafrost boreal forest ecosystem

Modeling hydrological processes in the Alaskan sub-arctic is challenging because of the extreme spatial heterogeneity in soil properties and vegetation communities. Nevertheless, modeling and predicting hydrological processes is critical in this region due to its vulnerability to the effects of climate change. Coarse-spatial-resolution datasets used in land surface modeling pose a new challenge in simulating the spatially distributed and basin-integrated processes since these datasets do not adequately represent the small-scale hydrological, thermal, and ecological heterogeneity. The goal of this study is to improve the prediction capacity of mesoscale to large-scale hydrological models by introducing a small-scale parameterization scheme, which better represents the spatial heterogeneity of soil properties and vegetation cover in the Alaskan sub-arctic. The small-scale parameterization schemes are derived from observations and a sub-grid parameterization method in the two contrasting sub-basins of the Caribou Poker Creek Research Watershed (CPCRW) in Interior Alaska: one nearly permafrost-free (LowP) sub-basin and one permafrost-dominated (HighP) sub-basin. The sub-grid parameterization method used in the small-scale parameterization scheme is derived from the watershed topography. We found that observed soil thermal and hydraulic properties – including the distribution of permafrost and vegetation cover heterogeneity – are better represented in the sub-grid parameterization method than the coarse-resolution datasets. Parameters derived from the coarse-resolution datasets and from the sub-grid parameterization method are implemented into the variable infiltration capacity (VIC) mesoscale hydrological model to simulate runoff, evapotranspiration (ET), and soil moisture in the two sub-basins of the CPCRW. Simulated hydrographs based on the small-scale parameterization capture most of the peak and low flows, with similar accuracy in both sub-basins, compared to simulated hydrographs based on the coarse-resolution datasets. On average, the small-scale parameterization scheme improves the total runoff simulation by up to 50 % in the LowP sub-basin and by up to 10 % in the HighP sub-basin from the large-scale parameterization. This study shows that the proposed sub-grid parameterization method can be used to improve the performance of mesoscale hydrological models in the Alaskan sub-arctic watersheds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2010

Beam Beam Simulation Studies for NLC And SLC2000

We apply and modify the computer codes CAIN(P.Chen, G.Horton-Smith, T.Ohgaki, A.W.Weidemann, K.Yokoya, contributed to Workshop on Gamma-Gamma Colliders, Berkeley, CA, March 28-31, 1994; SLAC-PUB-6583.) and GUINEAPIG(D.Schulte (DESY), unpublished.) to make detailed studies of the beam-beam interaction in the SLAC design for a future linear collider (NLC), as well as to the proposed SLC2000 project. Using realistic beam distributions, we present simulations related to the tuning and optimization of luminosity in SLC2000 and NLC.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

A numerical simulation of CCN impacts on weather modification efficiency

Aerosols affect development of clouds and precipitation by serving as cloud condensation nuclei (CCN) and ice nuclei (IN). Considering the dramatically changing ambient aerosol concentration, it is important to examine the potential “side effect” of aerosol pollution on precipitation enhancement by weather modification. In this study, the cloud seeding was performed on a precipitation event in Beijing in the summer of 2008, which is simulated by the NSSL two-moment cloud scheme of the Weather Research and Forecasting (WRF) model. Sensitivity tests were conducted by modifying the ambient aerosol concentration and the ice crystal seeding amount to investigate the cloud seeding efficacy in different CCN concentration scenarios. There was a slight difference in the precipitation distribution between the simulations with two ambient CCN concentrations: the northern precipitation center in polluted scenario was weaker and the southern center was stronger. Compared with normal CCN scenario, the cloud liquid water mass and ice crystal mass in the severe pollution scenario is larger, and the total contents of snow and graupel were not sensitive to the CCN concentration. With the same amount of man-made ice crystals seeding, the precipitation enhancement was quite different under different CCN conditions. The higher the CCN concentration usually leads to stronger precipitation suppression. As CCN concentration increase, the deposition growth of snow, auto-conversion and accretion of ice crystals to snow were weakened, as well as the conversion of melting snow and graupel into rainwater.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2011

A 2D simulation study of Langmuir, whistler, and cyclotron maser instabilities induced by an electron ring-beam distribution

We carried out a series of 2D simulations to study the beam instability and cyclotron maser instability (CMI) with the initial condition that a population of tenuous energetic electrons with a ring-beam distribution is present in a magnetized background plasma. In this paper, weakly relativistic cases are discussed with the ring-beam kinetic energy ranging from 25 to 100 keV. The beam component leads to the two-stream or beam instability at an earlier stage, and the beam mode is coupled with Langmuir or whistler mode, leading to excitation of beam-Langmuir or beam-whistler waves. When the beam velocity is large with a strong beam instability, the initial ring-beam distribution is diffused in the parallel direction rapidly. The diffused distribution may still support CMI to amplify the X1 mode (the fundamental X mode). On the contrary, when the beam velocity is small and the beam instability is weak, CMI can amplify the Z1 (the fundamental Z mode) effectively while the O1 (the fundamental O mode) and X2 (the second harmonic X mode) modes are very weak and the X1 mode is not excited. In this report, different cases with various parameters are presented and discussed for a comprehensive understanding of ring-beam instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Optimizing mini-ridge filter thickness to reduce proton treatment times in a spot-scanning synchrotron system

Purpose: Study the contributors to treatment time as a function of Mini-Ridge Filter (MRF) thickness to determine the optimal choice for breath-hold treatment of lung tumors in a synchrotron-based spot-scanning proton machine. Methods: Five different spot-scanning nozzles were simulated in TOPAS: four with MRFs of varying maximal thicknesses (6.15–24.6 mm) and one with no MRF. The MRFs were designed with ridges aligned along orthogonal directions transverse to the beam, with the number of ridges (4–16) increasing with MRF thickness. The material thickness given by these ridges approximately followed a Gaussian distribution. Using these simulations, Monte Carlo data were generated for treatment planning commissioning. For each nozzle, standard and stereotactic (SR) lung phantom treatment plans were created and assessed for delivery time and plan quality. Results: Use of a MRF resulted in a reduction of the number of energy layers needed in treatment plans, decreasing the number of synchrotron spills needed and hence the treatment time. For standard plans, the treatment time per field without a MRF was 67.0 ± 0.1 s, whereas three of the four MRF plans had treatment times of less than 20 s per field; considered sufficiently low for a single breath-hold. For SR plans, the shortest treatment time achieved was 57.7 ± 1.9 s per field, compared to 95.5 ± 0.5 s without a MRF. There were diminishing gains in time reduction as the MRF thickness increased. Dose uniformity of the PTV was comparable across all plans; however, when the plans were normalized to have the same coverage, dose conformality decreased with MRF thickness, as measured by the lung V20%. Conclusions: Single breath-hold treatment times for plans with standard fractionation can be achieved through the use of a MRF, making this a viable option for motion mitigation in lung tumors. For stereotactic plans, while a MRF can reduce treatment times, multiple breath-holds would still be necessary due to the limit imposed by the proton extraction time. To balance treatment time and normal tissue dose, the ideal MRF choice was shown to be the thinnest option that is able to achieve the desired breath-hold timing.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Liquid li structure and dynamics: A comparison between OFDFT and second nearest-neighbor embedded-atom method

The structure and dynamics of liquid lithium are studied using two simulation methods: orbital-free (OF) first-principles molecular dynamics (MD), which employs OF density functional theory (DFT), and classical MD utilizing a second nearest-neighbor embedded-atom method potential. The properties we studied include the dynamic structure factor, the self-diffusion coefficient, the dispersion relation, the viscosity, and the bond angle distribution function. Our simulation results were compared to available experimental data when possible. Each method has distinct advantages and disadvantages. For example, OFDFT gives better agreement with experimental dynamic structure factors, yet is more computationally demanding than classical simulations. Classical simulations can access a broader temperature range and longer time scales. The combination of first-principles and classical simulations is a powerful tool for studying properties of liquid lithium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗