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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

System-of-Systems Approach for Integrated Energy Systems Modeling and Simulation: Preprint

Today’s electricity grid is the most complex system ever built—and the future grid is likely to be even more complex because it will incorporate distributed energy resources (DERs) such as wind, solar, and various other sources of generation and energy storage. The complexity is further augmented by the possible evolution to new retail market structures that provide incentives to owners of DERs to support the grid. To understand and test new retail market structures and technologies such as DERs, demand-response equipment, and energy management systems while providing reliable electricity to all customers, an Integrated Energy System Model (IESM) is being developed at NREL. The IESM is composed of a power flow simulator (GridLAB-D), home energy management systems implemented using GAMS/Pyomo, a market layer, and hardware-in-the-loop simulation (testing appliances such as HVAC, dishwasher, etc.). The IESM is a system-of-systems (SoS) simulator wherein the constituent systems are brought together in a virtual testbed. We will describe an SoS approach for developing a distributed simulation environment. We will elaborate on the methodology and the control mechanisms used in the co-simulation illustrated by a case study.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

THE EVOLUTION OF ASTEROIDS IN THE JUMPING-JUPITER MIGRATION MODEL

In this work, we investigate the evolution of a primordial belt of asteroids, represented by a large number of massless test particles, under the gravitational effect of migrating Jovian planets in the framework of the jumping-Jupiter model. We perform several simulations considering test particles distributed in the Main Belt, as well as in the Hilda and Trojan groups. The simulations start with Jupiter and Saturn locked in the mutual 3:2 mean motion resonance plus three Neptune-mass planets in a compact orbital configuration. Mutual planetary interactions during migration led one of the Neptunes to be ejected in less than 10 Myr of evolution, causing Jupiter to jump by about 0.3 AU in semimajor axis. This introduces a large-scale instability in the studied populations of small bodies. After the migration phase, the simulations are extended over 4 Gyr, and we compare the final orbital structure of the simulated test particles to the current Main Belt of asteroids with absolute magnitude H < 9.7. The results indicate that, in order to reproduce the present Main Belt, the primordial belt should have had a distribution peaked at ∼10° in inclination and at ∼0.1 in eccentricity. We discuss the implications of this for the Grand Tack model. The results also indicate that neither primordial Hildas, nor Trojans, survive the instability, confirming the idea that such populations must have been implanted from other sources. In particular, we address the possibility of implantation of Hildas and Trojans from the Main Belt population, but find that this contribution should be minor.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Enrichment Determination of Uranium in Shielded Configurations

The determination of the enrichment of uranium is required in many safeguards and security applications. Typical methods of determining the enrichment rely on detecting the 186 keV gamma ray emitted by 235 U. In some applications, the uranium is surrounded by external shields, and removal of the shields is undesirable. In these situations, methods relying on the detection of the 186 keV gamma fail because the gamma ray is shielded easily. Oak Ridge National Laboratory (ORNL) has previously measured the enrichment of shielded uranium metal using active neutron interrogation. The method consists of measuring the time distribution of fast neutrons from induced fissions with large plastic scintillator detectors. To determine the enrichment, the measurements are compared to a calibration surface that is created from Monte Carlo simulations where the enrichment in the models is varied. In previous measurements, the geometry was always known. ORNL is extending this method to situations where the geometry and materials present are not known in advance. In the new method, the interrogating neutrons are both time and directionally tagged, and an array of small plastic scintillators measures the uncollided interrogating neutrons. Therefore, the attenuation through the item along many different paths is known. By applying image reconstruction techniques, an image of the item is created which shows the position-dependent attenuation. The image permits estimating the geometry and materials present, and these estimates are used as input for the Monte Carlo simulations. As before, simulations predict the time distribution of induced fission neutrons for different enrichments. Matching the measured time distribution to the closest prediction from the simulations provides an estimate of the enrichment. This presentation discusses the method and provides results from recent simulations that show the importance of knowing the geometry and materials from the imaging system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2011

Towards High Performance Discrete-Event Simulations of Smart Electric Grids

Future electric grid technology is envisioned on the notion of a smart grid in which responsive end-user devices play an integral part of the transmission and distribution control systems. Detailed simulation is often the primary choice in analyzing small network designs, and the only choice in analyzing large-scale electric network designs. Here, we identify and articulate the high-performance computing needs underlying high-resolution discrete event simulation of smart electric grid operation large network scenarios such as the entire Eastern Interconnect. We focus on the simulator's most computationally intensive operation, namely, the dynamic numerical solution for the electric grid state, for both time-integration as well as event-detection. We explore solution approaches using general-purpose dense and sparse solvers, and propose a scalable solver specialized for the sparse structures of actual electric networks. Based on experiments with an implementation in the THYME simulator, we identify performance issues and possible solution approaches for smart grid experimentation in the large.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2016

Evaluation of the aerosol vertical distribution in global aerosol models through comparison against CALIOP measurements: AeroCom phase II results

The ability of 11 models in simulating the aerosol vertical distribution from regional to global scales, as part of the second phase of the AeroCom model intercomparison initiative (AeroCom II), is assessed herein and compared to results of the first phase. The evaluation is performed using a global monthly gridded data set of aerosol extinction profiles built for this purpose from the CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) Layer Product 3.01. Results over 12 subcontinental regions show that five models improved, whereas three degraded in reproducing the interregional variability in Z α 0–6 km , the mean extinction height diagnostic, as computed from the CALIOP aerosol profiles over the 0–6 km altitude range for each studied region and season. While the models' performance remains highly variable, the simulation of the timing of the Z α 0–6 km peak season has also improved for all but two models from AeroCom Phase I to Phase II. The biases in Z α 0–6 km are smaller in all regions except Central Atlantic, East Asia, and North and South Africa. Most of the models now underestimate Z α 0–6 km over land, notably in the dust and biomass burning regions in Asia and Africa. At global scale, the AeroCom II models better reproduce the Z α 0–6 km latitudinal variability over ocean than over land. Hypotheses for the performance and evolution of the individual models and for the intermodel diversity are discussed. We also provide an analysis of the CALIOP limitations and uncertainties contributing to the differences between the simulations and observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Coupled TD Simulation project (Final Report)

This project implements coupled simulation of electric power grid transmission and distribution (T&D) by instantiating the co-simulation platform HELICS using a commercial distribution planning software CYMDIST1 by Eaton Inc., a power management company, and by leveraging the open-source transmission grid simulator GridDyn and LLNL’s High Performance Computing (HPC). The project is funded by DOE’s Technology Commercialization Funding (TCF) and supported by DOE Office of Technology Transitions (OTT).

14 SOLAR ENERGY↗
DOE OSTI2016

The charge imbalance in ultracold plasmas

Ultracold plasmas are regarded as quasineutral but not strictly neutral. The results of charge imbalance in the expansion of ultracold plasmas are reported. The calculations are performed by a full molecular-dynamics simulation. The details of the electron velocity distributions are calculated without the assumption of electron global thermal equilibrium and Boltzmann distribution. Spontaneous evolutions of the charge imbalance from the initial states with perfect neutrality are given in the simulations. The expansion of outer plasma slows down with the charge imbalance. The influences of plasma size and parameters on the charge imbalance are discussed. The radial profiles of electron temperature are given for the first time, and the self-similar expansion can still occur even if there is no global thermal equilibrium. The electron disorder induced heating is also found in the simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Monte Carlo calculation of the maximum therapeutic gain of tumor antivascular alpha therapy

Purpose: Metastatic melanoma lesions experienced marked regression after systemic targeted alpha therapy in a phase 1 clinical trial. This unexpected response was ascribed to tumor antivascular alpha therapy (TAVAT), in which effective tumor regression is achieved by killing endothelial cells (ECs) in tumor capillaries and, thus, depriving cancer cells of nutrition and oxygen. The purpose of this paper is to quantitatively analyze the therapeutic efficacy and safety of TAVAT by building up the testing Monte Carlo microdosimetric models. Methods: Geant4 was adapted to simulate the spatial nonuniform distribution of the alpha emitter {sup 213}Bi. The intraluminal model was designed to simulate the background dose to normal tissue capillary ECs from the nontargeted activity in the blood. The perivascular model calculates the EC dose from the activity bound to the perivascular cancer cells. The key parameters are the probability of an alpha particle traversing an EC nucleus, the energy deposition, the lineal energy transfer, and the specific energy. These results were then applied to interpret the clinical trial. Cell survival rate and therapeutic gain were determined. Results: The specific energy for an alpha particle hitting an EC nucleus in the intraluminal and perivascular models is 0.35 and 0.37 Gy, respectively. As the average probability of traversal in these models is 2.7% and 1.1%, the mean specific energy per decay drops to 1.0 cGy and 0.4 cGy, which demonstrates that the source distribution has a significant impact on the dose. Using the melanoma clinical trial activity of 25 mCi, the dose to tumor EC nucleus is found to be 3.2 Gy and to a normal capillary EC nucleus to be 1.8 cGy. These data give a maximum therapeutic gain of about 180 and validate the TAVAT concept. Conclusions: TAVAT can deliver a cytotoxic dose to tumor capillaries without being toxic to normal tissue capillaries.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2026

Use of Digital Real-Time Simulation and Optimization to Identify Maximum Real Power Injection on Banshee Distribution Network

This study investigates the hosting capacity of the Banshee Distribution Network by optimizing the real power injection at carefully selected Distributed Energy Resource (DER) locations. The analysis is conducted within the framework of power system operational constraints, including bus voltage ranges, thermal line ratings, and transformer loading limits. A Python-based Genetic Algorithm (GA), implemented using the PyGAD library, is employed to iteratively identify the optimal power injection configuration that maximizes network utilization while preserving system reliability. The methodology integrates a real-time simulation environment using the Real-Time Digital Simulator (RTDS), allowing high-fidelity evaluation of power flow and voltage behavior under each proposed injection scenario. By coupling the optimization algorithm with real-time simulation feedback, this approach ensures that both static and dynamic constraints are enforced during the evaluation process. The GA leverages evolutionary operators such as selection, crossover, and mutation to navigate the nonlinear search space efficiently. The results of the study delineate the feasible hosting capacity at three targeted buses, reflecting maximum real power levels that can be injected without causing voltage violations, transformer overloading, or line congestion. These findings provide a decision- support tool for distribution planners and utilities aiming to integrate higher penetrations of DERs in existing infrastructure. Additionally, the work lays the foundation for extending such optimization techniques to multi-objective formulations, including economic dispatch and reactive power coordination, in future studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Improving the Performance of Integrated Power-Hardware-in-the-Loop and Quasi-Static Time-Series Simulations

This paper introduces a novel multirate cosimulation architecture that overcomes previous challenges integrating faster real-time, microseconds-scale, power-hardware-in-the-loop (PHIL) with larger scale but slower, real-time, seconds-scale, quasi-static time-series (QSTS) simulation. Specifically, an intermediate reduced-equivalent electromagnetic transient (EMT) model duplicates the key power system topology to capture high-speed dynamics and converge hardware-power system interactions between QSTS updates. Exchanging multidimensional parameter vectors (loads, control status, etc.) between the QSTS and EMT models enables capturing the interactions of the rich, high-node-count (thousands of electrical nodes) QSTS model with the hardware. This architecture offers full spatial resolution from QSTS, including individual load dynamics, actual distribution management system controls, and nodal voltages. Simultaneously, the intermediate EMT model provides more detailed high-speed transients of key distributed energy resource hardware interactions. In addition, we use careful PHIL interface design and the exchange of complex power data, rather than current, to further improve performance. This cosimulation architecture is demonstrated by testing a simulated distribution system with an interconnected 500-kVA advanced photovoltaic (PV) inverter in PHIL. In conclusion, the architecture successfully captured PV local volt-volt ampere reactive (volt-VAR) control interactions with the larger network under time-varying electrical and weather conditions.

42 ENGINEERING↗
DOE OSTI2011

Distributed non-equilibrium Green’s function algorithms for the simulation of nanoelectronic devices with scattering

Through the non-equilibrium Green’s function (NEGF) formalism, quantum-scale device simulation can be performed with the inclusion of electron-phonon scattering. Yet, the simulation of realistically sized devices under the NEGF formalism typically requires prohibitive amounts of memory and computation time. Two of the most demanding computational problems for NEGF simulation involve mathematical operations with structured matrices called semiseparable matrices. In this work, we present parallel approaches for these computational problems which allow for efficient distribution of both memory and computation based upon the underlying device structure. This is vital when simulating realistically sized devices due to the aforementioned computational burdens. First, we consider determining a distributed compact representation for the retarded Green’s function matrix G R . This compact representation is exact and allows for any entry in the matrix to be generated through the inherent semiseparable structure. The second parallel operation allows for the computation of electron density and current characteristics for the device. Specifically, matrix products between the distributed representation for the semiseparable matrix G R and the self-energy scattering terms in Σ < produce the less-than Green’s function G < . As an illustration of the computational efficiency of our approach, we stably generate the mobility for nanowires with cross-sectional sizes as large as 4.5 nm, assuming an atomistic model with scattering.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Wave-optical assessment of alignment tolerances in nano-focusing with ellipsoidal mirror

High-precision ellipsoidal mirrors, which can efficiently focus X-rays to the nanometer dimension with a mirror, have not been realized because of the difficulties in the fabrication process. The purpose of our study was to develop nano-focusing ellipsoidal mirrors in the hard X-ray region. We developed a wave-optical focusing simulator for investigating alignment tolerances in nano-focusing with a designed ellipsoidal mirror, which produce a diffraction-limited focus size of 30 × 35 nm{sup 2} in full width at half maximum at an X-ray energy of 7 keV. The simulator can calculate focusing intensity distributions around the focal point under conditions of misalignment. The wave-optical simulator enabled the calculation of interference intensity distributions, which cannot be predicted by the conventional ray-trace method. The alignment conditions with a focal length error of ≲ ±10 µm, incident angle error of ≲ ±0.5 µrad, and in-plane rotation angle error of ≲ ±0.25 µrad must be satisfied for nano-focusing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2010

Convective stability of a plasma in a system of coupled adiabatic open cells in the Kruskal-Oberman model

The Kruskal-Oberman kinetic model is used to determine the conditions for the convective stability of a plasma in a system of coupled axisymmetric adiabatic open cells in which the magnetic field curvature has opposite signs. For a combination of a nonparaxial simple mirror cell and a semicusp, the boundaries of the interval of values of the flux coordinate where the plasma can be stable are determined, as well as the range in which the ratio of the pressures in the component cells should lie. Numerical simulations were carried out for different particle distributions over the pitch angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

TYPE III RADIO BURSTS IN CORONAL PLASMAS WITH KAPPA PARTICLE DISTRIBUTIONS

We present the first simulations of type III bursts produced in the corona with suprathermal non-Maxwellian background particles, as inferred from solar wind data and proposed by theories for the corona and solar wind. The coronal background particles are assumed to follow kappa ({kappa}) distributions. The predicted f{sub p} emission of type III bursts is sensitive via the {kappa} index to the presence of suprathermal background particles, where f{sub p} is the local plasma frequency. The simulations show that (1) the speeds v{sub b} of type III beams are much larger (e.g., v{sub b} Almost-Equal-To 0.58c for {kappa} = 5) and so type III bursts drift much faster for low {kappa} ({<=}5) background plasmas than for Maxwellian backgrounds (producing v{sub b} < 0.3c), and (2) f{sub p} emission generated in a {kappa}-distributed background corona has a larger total bandwidth than in a Maxwellian background, for similar onset frequencies. Type III beams are thus more persistent, i.e., extending over larger distances, in {kappa}-distributed corona. Consequently, observations of fast-drifting coronal type III bursts and associated fast electron beams suggest that the ambient electrons in the corona are {kappa}-distributed, at least when such bursts are observed. These results support, from the new viewpoint of nonthermal radio emission, the occasional presence of suprathermal background electrons in the corona and the associated mechanisms (e.g., 'velocity filtration') for coronal heating and solar wind acceleration. The new results also help resolve longstanding issues regarding the speeds and persistence of type III beams, and the production of remotely observable levels of f{sub p} emission despite severe losses during propagation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Large-scale atomistic model construction of subbituminous and bituminous coals for solvent extraction simulations with reactive molecular dynamics

Large-scale atomistic models for complex polycyclic aromatic hydrocarbon systems help understand the chemical properties and behaviors of complex feedstocks such as coal or petroleum. However, the development and utilization of large-scale models remain limited due to the difficulty in achieving the varied structural characteristics necessary to capture stochastic nature of these feedstocks. Here we demonstrate a systematic workflow to construct stochastic molecular systems from a broad analytical suite: high-resolution transmission electron microscopy (HRTEM), carbon-13 nuclear magnetic resonance spectroscopy ( 13 C NMR), laser desorption ionization mass spectroscopy (LDI-MS), and elemental analysis. We present a model construction and analysis utility of a new Python-based module. We selected one subbituminous and three high-volatile bituminous coals to construct large-scale models (~40,000 atoms). The constructed models were utilized to examine the affinity for solvent extraction (naphthalene or tetralin) and the effect of structural properties (e.g., aromatic cluster size, functional groups, and cross-linking) in reactive molecular dynamics simulations. Complex chemical reactions were monitored with bond order transitions, intermediates formation, and mass distributions. Reactive molecular dynamics simulations suggest a plausible chemical extraction process and products for the complex fossil feedstocks. The results indicated that radical formations with bond breaking of bridging oxygens and carbons were required at high temperatures to facilitate hydrogeneration and extraction of gas molecules from radical-free molecules. We observed that aliphatic chains of tetralin were easily decomposed and combined with radicals to form small size of molecules with aryl bonding, mainly increasing molecules in the 500–1000 Da, while naphthalene had little impact on chemical extraction process.

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

Slow and static spin correlations in Dy(2+x)Ti(2-x)O(7-d)

The static and dynamic spin correlations in the spin ices Dy{sub 2.3}Ti{sub 1.7}O{sub 6.85} and Dy{sub 2}Ti{sub 2}O{sub 7} have been studied in polarized neutron diffraction and neutron spin echo experiments. The measurements reveal that, below 100 mK, the magnetic scattering braodens and shifts to higher |Q| upon stuffing the pyrochlore lattice with additional Dy{sup 3+} ions. These observations can be related, by means of reverse Monte Carlo simulation, to the modified distribution of near-neighbour distances and an overall more antiferromagnetic character of the near-neighbour couplings. The dynamic measurements show that the spin correlations are slower in the stuffed system. These results will be discussed and compared to the holmium analogues.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Approximating Density Probability Distribution Functions Across Cosmologies

Using a suite of self-similar cosmological simulations, we measure the probability distribution functions (PDFs) of real-space density, redshift-space density, and their geometric mean. We find that the real-space density PDF is well-described by a function of two parameters: ns, the spectral slope, and σL, the linear rms density fluctuation. For redshift-space density and the geometric mean of real- and redshift-space densities, we introduce a third parameter, ${s}_{L}=\sqrt{\langle {\left({{dv}}_{\mathrm{pec}}^{L}/{dr}\right)}^{2}\rangle }/H$. We find that density PDFs for the LCDM cosmology is also well-parameterized by these three parameters. As a result, we are able to use a suite of self-similar cosmological simulations to approximate density PDFs for a range of cosmologies. We make the density PDFs publicly available and provide an analytical fitting formula for them.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

An assessment of tropical cyclones in North American $\mathrm{CORDEX}$ $\mathrm{WRF}$ simulations

This work presents an assessment of tropical cyclones (TCs) in the 25 km and 50 km resolution reanalysis-forced and baseline and future (RCP8.5) global climate model (GCM) forced simulations produced for the North American branch of the international Coordinated Regional climate Downscaling Experiment (NA-CORDEX) using the Weather Research and Forecasting (WRF) model. A set of complementary 12 km resolution simulations produced as a part of a different project is also included in this assessment. Before examining the projections from the GCM-driven simulations, the ability of the simulations to minimally produce a realistic spatial distribution of historical TC occurrence was assessed in simulations forced by reanalysis and the three different GCMs used herein. Then, projections for occurrence, TC related mean precipitation and precipitation intensity, storm duration, the intensity measured by minimum pressure and maximum wind speed, storm size, and translation speed were examined. Several of these characteristics show little to no change in the future in trend or in distribution across the ensemble. However, many simulations suggest a westward shift or increase in TC occurrence over the East Pacific basin and a decrease in occurrence over the Caribbean and Gulf of Mexico. Increases (decreases) in total storm-related precipitation are projected where TC occurrence increases (decreases). TC precipitation intensity is found to increase in all simulations over the East Pacific, but projections are mixed over the North Atlantic. Finally, the ensemble projects a distribution shift towards more intense TC over the East Pacific, and a shift toward faster translation speeds over the North Atlantic.

54 ENVIRONMENTAL SCIENCES↗