Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulator”

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 685 records · Page 38

Secrets of the Dark Universe: Simulating the Sky on the Blue Gene/Q, The Outer Rim Simulation

An astonishing 99.6% of our Universe is dark. Observations indicate that the Universe consists of 70% of a mysterious dark energy and 25% of a yet-unidentified dark matter component, and only 0.4% of the remaining ordinary matter is visible. Understanding the physics of this dark sector is the foremost challenge in cosmology today. Sophisticated simulations of the evolution of the Universe play a crucial task in this endeavor. This movie shows an intermediate stage in a large simulation of the distribution of matter in the Universe, the so-called cosmic web, accounting for the influence of dark energy. The simulation is evolving 1.1 trillion particles. The movie shows a snapshot of the Universe when it was 1.6 billion years old. Credits: Science: Hal Finkel, Salman Habib, Katrin Heitmann, Kalyan Kumaran, Vitali Morozov, Tom Peterka, Adrian Pope, Tim Williams, David Daniel, Patricia Fasel, Nicholas Frontiere and Zarija Lukic. Visualization: Mark Herald, Joseph A. Insley, Michael E. Papka and Venkatram Vishwanath

79 ASTRONOMY AND ASTROPHYSICS↗

Adapting the serial Alpgen parton-interaction generator to simulate LHC collisions on millions of parallel threads

As the LHC moves to higher energies and luminosity, the demand for computing resources increases accordingly and will soon outpace the growth of the World- wide LHC Computing Grid. To meet this greater demand, event generation Monte Carlo was targeted for adaptation to run on Mira, the supercomputer at the Argonne Leadership Computing Facility. Alpgen is a Monte Carlo event generation application that is used by LHC experiments in the simulation of collisions that take place in the Large Hadron Collider. This paper details the process by which Alpgen was adapted from a single-processor serial-application to a large-scale parallel-application and the performance that was achieved.

HEP↗

Validation of Phasor-Domain Transmission and Distribution Co-simulation Against Electromagnetic Transient Simulation

The rapid deployment of renewable energy resources has led to the widespread use of power electronics in modern power systems. As these systems transition from being dominated by large synchronous machines to increasingly incorporating inverter-based resources (IBRs), traditional methods are becoming inadequate. Addressing this challenge, this paper introduces a scalable phasor-domain T\&D co-simulation framework based on open-source software. It focuses on the framework's validation against the PSCAD Electromagnetic Transient (EMT) analysis tool. The validation results demonstrate the framework's high-fidelity and a computational time speed-up of 60 to 100 times, marking a pioneering validation effort in T\&D co-simulation research.

Inverter-based resources, co-simulation, Electroma↗

Cascade annealing simulations of bcc iron using object kinetic Monte Carlo

Simulations of displacement cascade annealing were carried out using object kinetic Monte Carlo based on an extensive MD database including various primary knock-on atom energies and directions. The sensitivity of the results to a broad range of material and model parameters was examined. The diffusion mechanism of interstitial clusters has been identified to have the most significant impact on the fraction of stable interstitials that escape the cascade region. The maximum level of recombination was observed for the limiting case in which all interstitial clusters exhibit 3D random walk diffusion. The OKMC model was parameterized using two alternative sets of defect migration and binding energies, one from ab initio calculations and the second from an empirical potential. The two sets of data predict essentially the same fraction of surviving defects but different times associated with the defect escape processes. This study provides a comprehensive picture of the first phase of long-term defect evolution in bcc iron and generates information that can be used as input data for mean field rate theory (MFRT) to predict the microstructure evolution of materials under irradiation. In addition, the limitations of the current OKMC model are discussed and a potential way to overcome these limitations is outlined.

36 MATERIALS SCIENCE↗

Monte Carlo simulations of channeling spectra recorded for samples containing complex defects

The main aim of the present paper is to describe the current status of the development of McChasy, a Monte Carlo simulation code, to make it suitable for the analysis of dislocations and dislocation loops in crystals. Such factors like the shape of the bent channel and geometrical distortions of the crystalline structure in the vicinity of dislocation has been discussed. Several examples of the analysis performed at different energies of analyzing ions are presented. The results obtained demonstrate that the new procedure applied to the spectra recorded on crystals containing dislocation yields damage profiles which are independent of the energy of the analyzing beam.

Carlo↗

A Stereo Camera Simulator for Large-Eddy Simulations of Continental Shallow Cumulus Clouds Based on Three-Dimensional Path-Tracing

The complex spatial and temporal structure of cumulus clouds complicates their representation in weather and climate models. Classic meteorological instrumentation struggles to fully capture these features. Networks of multiple high-resolution hemispheric cameras are increasingly used to fill this data gap, and provide information on this missing multi-dimensional spatial information. In this study, a path-tracing algorithm is used to generate virtual camera images of resolved clouds in large-eddy simulations (LES). These images are then used as a camera network simulator, allowing reconstructions of three-dimensional cloud edges from the model output. Because the actual LES cloud field is fully known, the combined path-tracing and reconstruction method can be statistically analyzed. The method is applied to LES realizations of summertime shallow cumulus at the Jülich Observatory for Cloud Evolution (JOYCE), Germany, which also routinely operates a camera network. We find that the path-tracing method allows accurate reconstruction of up to 70% of the visible cloud edges. Additional sensitivity tests show that the method is robust for changes in its hyperparameters. The sensitivity to cloud optical thickness is also investigated, finding a cloud boundary placement error of approximately 182 m. This error can be considered typical for cloud boundary reconstruction using real stereo camera imagery. The results provide proof of principle for future use of the method for evaluating LES clouds against camera network imagery, and for further optimizing the configuration of such camera networks.

54 ENVIRONMENTAL SCIENCES↗

Modeling of 10 GeV-1 TeV laser-plasma accelerators using Lorentz boosted simulations

Modeling of laser-plasma wakefield accelerators in an optimal frame of reference [J.-L. Vay, Phys. Rev. Lett. 98, 130405 (2007)] allows direct and efficient full-scale modeling of deeply depleted and beam loaded laser-plasma stages of 10 GeV-1 TeV (parameters not computationally accessible otherwise). This verifies the scaling of plasma accelerators to very high energies and accurately models the laser evolution and the accelerated electron beam transverse dynamics and energy spread. Over 4, 5, and 6 orders of magnitude speedup is achieved for the modeling of 10 GeV, 100 GeV, and 1 TeV class stages, respectively. Agreement at the percentage level is demonstrated between simulations using different frames of reference for a 0.1 GeV class stage. Obtaining these speedups and levels of accuracy was permitted by solutions for handling data input (in particular, particle and laser beams injection) and output in a relativistically boosted frame of reference, as well as mitigation of a high-frequency instability that otherwise limits effectiveness.

43 PARTICLE ACCELERATORS↗

Monte Carlo simulations of electron thermalization in alkali iodide and alkaline-earth fluoride scintillators

A Monte Carlo model of electron thermalization in inorganic scintillators, which was developed and applied to CsI in a previous publication [Wang et al., J. Appl. Phys. 110, 064903 (2011)], is extended to another material of the alkali halide class, NaI, and to two materials from the alkaline-earth halide class, CaF{sub 2} and BaF{sub 2}. This model includes electron scattering with both longitudinal optical (LO) and acoustic phonons as well as the effects of internal electric fields. For the four pure materials, a significant fraction of the electrons recombine with self-trapped holes and the thermalization distance distributions of the electrons that do not recombine peak between approximately 25 and 50 nm and extend up to a few hundreds of nanometers. The thermalization time distributions of CaF{sub 2}, BaF{sub 2}, NaI, and CsI extend to approximately 0.5, 1, 2, and 7 ps, respectively. The simulations show that the LO phonon energy is a key factor that affects the electron thermalization process. Indeed, the higher the LO phonon energy is, the shorter the thermalization time and distance are. The thermalization time and distance distributions show no dependence on the incident {gamma}-ray energy. The four materials also show different extents of electron-hole pair recombination due mostly to differences in their electron mean free paths (MFPs), LO phonon energies, initial densities of electron-hole pairs, and static dielectric constants. The effect of thallium doping is also investigated for CsI and NaI as these materials are often doped with activators. Comparison between CsI and NaI shows that both the larger size of Cs{sup +} relative to Na{sup +}, i.e., the greater atomic density of NaI, and the longer electron mean free path in NaI compared to CsI contribute to an increased probability for electron trapping at Tl sites in NaI versus CsI.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Numerical simulation of cathode plasma dynamics in magnetically insulated vacuum transmission lines

A novel algorithm for the simulation of cathode plasmas in particle-in-cell codes is described and applied to investigate cathode plasma evolution in magnetically insulated transmission lines (MITLs). The MITL electron sheath is modeled by a fully kinetic electron species. Electron and ion macroparticles, both modeled as fluid species, form a dense plasma which is initially localized at the cathode surface. Energetic plasma electron particles can be converted to kinetic electrons to resupply the electron flux at the plasma edge (the “effective” cathode). Using this model, we compare results for the time evolution of the cathode plasma and MITL electron flow with a simplified (isothermal) diffusion model. Simulations in 1D show a slow diffusive expansion of the plasma from the cathode surface. But in multiple dimensions, the plasma can expand much more rapidly due to anomalous diffusion caused by an instability due to the strong coupling of a transverse magnetic mode in the electron sheath with the expanding resistive plasma layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mechanism of Void Nucleation and Growth in bcc Fe: Atomistic Simulations at Experimental Time Scales

Evolution of small-vacancy clusters in bcc Fe is simulated using a multiscale approach coupling an atomistic activation-relaxation method for sampling transition-state pathways with environment-dependent reaction coordinate calculations and a kinetic Monte Carlo simulation to reach time scales on the order of {approx}10{sup 4} s. Under vacancy-supersaturated condition, di- and trivacancy clusters form and grow by coalescence (Ostwald ripening). For cluster size greater than four we find a transition temperature of 150 deg. C for accelerated cluster growth, as observed in positron annihilation spectroscopy experiments. Implications for the mechanism of stage-IV radiation-damage-recovery kinetics are discussed.

36 MATERIALS SCIENCE↗

Monte Carlo simulations of photoneutron activation in a medical LINAC

The first step of the present study was to simulate the γ-radioactivity produced by photoneutrons in the treatment head of a decommissioned medical LINAC. These results, obtained with MCNPX and CINDER'90, are then compared with γ-ray spectrometric measurements of parts of the treatment head. Our ultimate objective is to simulate the associated production by photonuclear reactions of pure β-emitters, which are more difficult to access experimentally than the γ-ray emitters. (authors)

43 PARTICLE ACCELERATORS↗

Novel Approach to PV Inverter Modeling and Simulation Leveraging Experiments, Learning Based Modeling and Co-Simulation

Photovoltaic (PV) inverter manufacturers use custom, proprietary control approaches and topologies in their inverter design. The proprietary nature of these approaches makes it challenging to share electromagnetic transients (EMT) domain models for system studies. This research work presents an approach to develop EMT models from experimental data. We use novel approach in experimental design, high fidelity data collection, use of learning-based modeling, and co-simulation to reduce the time taken to develop an EMT model for an inverter under test (IUT). We used a 20 kW off-the-shelf grid following PV inverter and subjected the inverter to controlled tests. The tests include voltage and frequency step changes, as well as solar irradiance variations. The recorded high frequency data were used to train a neural network model representing the dynamic behavior of the IUT. The model was subsequently imported into an EMT tool using co-simulation techniques, and thus completing the modeling effort.

black box inverter modeling↗

Impact of Detailed Parameter Modeling of Open-Cycle Gas Turbines on Production Cost Simulation: Preprint

Flexible resources are increasingly important as variable renewable energy deployment in the power system increases. Although many systems are transitioning away from fossil fuels, open-cycle gas turbines are likely to play an important balancing role for some time, thus requiring accurate modeling of their operational parameters. This paper explores the impact of detailed representation of three operational parameters - start- up costs, run-up rates, and forced outage rates - in the production cost model of a system as it adopts higher levels of wind and solar. Using PLEXOS simulations of the NREL-118 bus test system, the study examines how more detailed parameter modeling affects outcomes such as the number of start-ups and shutdowns, ramping and total generation costs for open-cycle gas turbines, as renewable energy levels increase. The results suggest the value of detailed parameter modeling and continued research on combustion turbines' ability to provide flexibility.

economic dispatch↗

A Blockchain Based Co-Simulation Framework for Integrating DERs to Wholesale Electricity Markets

As the number of distributed energy resources (DERs) continue to increase across energy-delivery systems, there remains a need for integrating their capabilities into traditional grid operations. In this paper, a blockchain-based solution is proposed to facilitate FERC's Order No. 2222 implementations. The presented use-case enables small-scale DERs to participate in wholesale market operations through DER aggregators, while also enabling local distribution system operators to enforce distribution system constraints in a secure and traceable manner. The presented use case is built around the Unified Testing Platform (UTP) being developed as a part of the Blockchain for Optimized Security and Energy Management (BLOSEM) project. This is a multi-lab effort intended to simplify the deployment of blockchain-powered grid solutions by enabling the integration of simulation tools, and blockchain technologies through the use of system-agnostic interfaces that provide a modular, interoperable, and reusable connectivity layer.

blockchain interoperability↗