Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulators”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Hanford Composite Analysis Special Studies FY2025 - Groundwater Flow Simulation for use in a Comparison of Simulated Concentrations in the Saturated Zone Estimated by the Plateau-to-River Model Versions 8.3 and 9.1

In order to meet the requirements of U.S. Department of Energy (DOE) in DOE O 435.1, Radioactive Waste Management, a composite analysis (CA) must be completed for the Central Plateau at the Hanford Site. The CA requires estimates of fate and transport of radionuclides in the groundwater from multiple sources within the modeling domain. This Environmental Calculation File (ECF) details the application of the Plateau-to-River (P2R) Model (CP-57037, Rev. 3, Model Package Report for the P2R Model: Version 9.1) to predict the flow of groundwater on the Central Plateau for the 10,000-year simulation to support a special analysis associated with the Hanford Site CA. The simulated flow field will support the simulation of fate and transport of contaminants for comparison to the results obtained as part of the CA to evaluate impact of the updated version of the P2R Model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

ILAMBv2.7 benchmarking results comparing E3SMv2.1 land-atmosphere coupled (BGCv2LNDATM) and stand alone land (ELM) simulations with CMIP6 emission driven historical simulations

This dataset contains land model benchmarking results for the Energy Exascale Earth System Model version 2.1 (E3SMv2.1), including outputs from both coupled biogeochemistry simulations and stand-alone land model simulations. These results are compared against several emission-driven historical simulations from the Coupled Model Intercomparison Project Phase 6 (CMIP6). Benchmarking was conducted using the International Land Model Benchmarking (ILAMB) package, version 2.7 (ILAMBv2.7). CMIP6 model outputs were sourced from the Earth System Grid Federation (ESGF), while the E3SMv2.1 results were derived from raw model outputs. These outputs underwent processing steps such as time serialization, conservative regridding, and data standardization to ensure comparability. For spatial interpolation, the Earth System Modeling Framework (ESMF) tool, ESMF_RegridWeightGen, was employed to generate regridding weights, enabling the transformation of E3SM’s native cubed-sphere grid to a regular latitude-longitude grid.

Feng, Sha [PNNL]↗
DOE OSTI2025

The AGORA High-resolution Galaxy Simulations Comparison Project. VIII. Disk Formation and Evolution of Simulated Milky Way Mass Galaxy Progenitors at 1 < z < 5

We investigate how differences in the stellar feedback produce disks with different morphologies in Milky Way–like progenitors over 1 ≤ z ≤ 5, using eight state-of-the-art cosmological hydrodynamics simulation codes in the AGORA project. In three of the participating codes, a distinct, rotation-dominated inner core emerges with a formation timescale of ≲300 Myr, largely driven by a major merger event, while two other codes exhibit similar signs of wet compaction—gaseous shrinkage into a compact starburst phase—at earlier epochs. The remaining three codes show only weak evidence of wet compaction. Consequently, we divide the simulated galaxies into two groups: those with strong compaction signatures and those with weaker ones. Galaxies in these two groups differ in size, stellar age gradients, and disk-to-total mass ratios. Specifically, codes with strong wet compaction build their outer disks in an inside-out fashion, leading to negative age gradients, whereas codes with weaker compaction feature flat or positive age gradients caused primarily by outward stellar migration. Although the stellar half-mass radii of these two groups diverge at z ∼ 3, the inclusion of dust extinction brings their sizes and shapes in mock observations closer to each other and to observed galaxies. We attribute the observed morphological differences primarily to variations in the stellar feedback implementations—such as delayed cooling timescales, and feedback strengths—that regulate both the onset and duration of compaction. Overall, our results suggest that disk assembly at high redshifts is highly sensitive to the details of the stellar feedback prescriptions in simulations.

Jung, Minyong [Seoul National Univ. (Korea, Republ↗
DOE OSTI2013

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↗
DOE OSTI2017

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↗
DOE OSTI2024

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↗
DOE OSTI2012

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↗
DOE OSTI2012

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↗
DOE OSTI2024

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↗
DOE OSTI2011

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↗
DOE OSTI2012

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↗
DOE OSTI2015

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↗
DOE OSTI2015

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↗
DOE OSTI2025

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↗
DOE OSTI2023

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↗
DOE OSTI2022

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↗