Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “production simulation codes”

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 199 records · Page 11

Improving resbos for the precision needs of the LHC

The resummation calculation (esos) is a widely used tool for the simulation of single vector boson production at colliders. In this work, we develop a significant improvement over the esos code by increasing the accuracy from NNLL + NLO to N 3 LL + NNLO and release the esos v2.0 code. Furthermore, we propose a new nonperturbative function that includes information about the rapidity of the system (IFY). The IFY functional form was fitted to data from fixed target experiments, the Tevatron, and the LHC. We find that the nonperturbative function has mild rapidity dependence based on the results of the fit. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Determining mode excitations of vacuum electronics devices via three-dimensional simulations using the SOS code

The SOS code is used to compute the resonance modes (frequency-domain information) of sample devices and separately to compute the transient behavior of the same devices. A code, DOT, is created to compute appropriate dot products of the time-domain and frequency-domain results. The transient behavior of individual modes in the device is then plotted. Modes in a coupled-cavity traveling-wave tube (CCTWT) section excited beam in separate simulations are analyzed. Mode energy vs. time and mode phase vs. time are computed and it is determined whether the transient waves are forward or backward waves for each case. Finally, the hot-test mode frequencies of the CCTWT section are computed.

Warren, Gary↗

Geothermal Reservoir Simulation Results in support of Feasibility Study of Direct District Heating for the Cornell Campus Utilizing Deep Geothermal Energy

This dataset contains input data, code, ReadMe files, output data, and figures that summarize the results of a stochastic analysis of geothermal reservoir production from two potential geothermal reservoirs that were evaluated for the Cornell University Deep Direct-Use project. These potential reservoirs are the Trenton-Black River (TBR) from 2.27-2.3 km depth, and basement rocks from 3.0-3.5 km depth and 3.5-4.0 km depth. Several utilization scenarios consisting of different injection fluid temperatures and flow rates were evaluated for each reservoir. Uncertainty in geologic properties, thermal properties, economic costs, and utilization efficiencies were evaluated using a Monte Carlo analysis of the reservoir simulations. Some reservoir simulations of the TBR were completed using the TOUGH2 software, as implemented in PetraSIM. The PetraSIM run files and associated data are provided with this submission. All other reservoir simulations were completed using the GEOPHIRES software, with some modifications to complete the uncertainty analyses. ReadMe files that describe additions to GEOPHIRES, the GEOPHIRES input data, and the output data are all provided, and references are provided to the code repository. Figures that summarize the reservoir heat production, temperature drawdown, and the probability of meeting targeted building heating demands with the produced heat and fluid temperatures are provided.

15 GEOTHERMAL ENERGY↗

Final Report - PyVERAIn: An Open Source, Python-based CASL Preprocessing Platform for Reactor Geometry, Code Discretizations, and Parallel Domain Decompositions

The introduction of advanced modeling and simulation practices over the last several decades has proven to produce better products, more quickly, and at the reduced cost. This computational experimentation has enabled the design and development of a variety of innovative products by reducing the need to build expensive prototypes, and by allowing rapid, comprehensive design studies across a wide range of products. Despite these successes, there are significant barriers to widespread adoption of these tools. This project aims to address several of these barriers.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

LACEwING: A New Moving Group Analysis Code

We present a new nearby young moving group (NYMG) kinematic membership analysis code, LocAting Constituent mEmbers In Nearby Groups (LACEwING), a new Catalog of Suspected Nearby Young Stars, a new list of bona fide members of moving groups, and a kinematic traceback code. LACEwING is a convergence-style algorithm with carefully vetted membership statistics based on a large numerical simulation of the Solar Neighborhood. Given spatial and kinematic information on stars, LACEwING calculates membership probabilities in 13 NYMGs and three open clusters within 100 pc. In addition to describing the inputs, methods, and products of the code, we provide comparisons of LACEwING to other popular kinematic moving group membership identification codes. As a proof of concept, we use LACEwING to reconsider the membership of 930 stellar systems in the Solar Neighborhood (within 100 pc) that have reported measurable lithium equivalent widths. We quantify the evidence in support of a population of young stars not attached to any NYMGs, which is a possible sign of new as-yet-undiscovered groups or of a field population of young stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Integration of the Nek5000 Computational Fluid Dynamics Code with the NEAMS Workbench

The mission of the US Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program is to develop, apply, deploy, and support state-of-the-art predictive modeling and simulation tools for the design and analysis of current and future nuclear energy systems using computing architectures ranging from laptops to leadership-class facilities. NEAMS is organized into three product lines: the Fuels Product Line (FPL),the Reactors Product Line (RPL), and the Integration Product Line (IPL). The FPL and RPL focus on the development of high-fidelity tools such as the MOOSE-based applications BISON and MARMOT [3, 9] at Idaho National Laboratory, and the PROTEUS neutronics code and the NEK5000 computational fluid dynamics code from the Simulation-based High-efficiency Advanced Reactor Prototyping (SHARP) framework [10]. The NEAMS IPL will enable users to apply high-fidelity simulations to inform lower order models for the design, analysis, and licensing of advanced nuclear systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Optimization of the beam extraction systems for the Linac4 H{sup −} ion source

The development of the Linac 4 and its integration into CERN’s acceleration complex is part of the foreseen luminosity upgrade of the Large Hadron Collider (LHC). The goal is to inject a 160 MeV H{sup −} beam into the CERN PS Booster (PSB) in order to increase the beam brightness by a factor of 2 compared to the 50 MeV proton linac, Linac 2, that is currently in operation. The requirements for the ion source are a 45 keV H{sup −} beam of 80 mA intensity, 2 Hz repetition rate and 0.5 ms pulse length within a normalized rms-emittance of 0.25 mm· mrad. The previously installed beam extraction system has been designed for an H{sup −} ion beam intensity of 20 mA produced by an RF-volume source with an electron to H{sup −} ratio of up to 50. For the required intensity upgrades of the Linac4 ion source, a new beam extraction system is being produced and tested; it is optimized for a cesiated surface RF-source with a nominal beam current of 40 mA and an electron to H{sup −} ratio of 4. The simulations, based on the IBSIMU code, are presented. At the Brookhaven National Laboratory (BNL), a peak beam current of more than 100 mA was demonstrated with a magnetron H{sup −} source at an energy of 35 keV and a repetition rate of 2 Hz. A new extraction system is required to operate at an energy of 45 keV; simulation of a two stage extraction system dedicated to the magnetron is presented.

43 PARTICLE ACCELERATORS↗

Nonlinear light-output calibration of the oxygenated xylene scintillators used in OMEGA neutron time-of-flight spectrometers

Neutron time-of-flight (nTOF) spectrometers are essential instruments for measuring and evaluating the performance of inertial confinement fusion implosions. The neutron spectrometers utilized for the OMEGA laser include two liquid-based scintillators, each consisting of a large volume filled with xylene that is coupled to four photomultiplier tubes. Analysis of the signal from these detectors requires detailed knowledge of the scintillator’s light output, which is needed to fit the nTOF spectrum, from which the neutron energy spectrum is informed. The light output is nonlinearly proportional to the neutron energy, which, in turn, affects the interpretation of the neutron energy spectrum from a TOF signal. A recent campaign on OMEGA was performed to calibrate the xylene detectors and infer the shape of the light-output curve. The campaign utilized materials with increasing Z placed in the OMEGA target chamber to initiate scattering events with the 14 MeV fusion neutrons. This process leads to the production of backscatter neutrons of varying energies that appear as peaks in the nTOF data. Simulations using a neutron transport code were combined with the measured deuterium–tritium neutron yields to calculate the expected backscattered neutron yields from the well-known scattering cross sections of each material. Furthermore, the neutron-energy dependent light output of the scintillator inferred from the experiment is compared to the light-output curve simulated with a neutron transport code for the following neutron energies: 1.5, 2.5, 6, and 14 MeV.

47 OTHER INSTRUMENTATION↗

Code Generators for Floating-Point Unit Design in Integrated Circuits (OpenFloat) v1.0

This IP provides a comprehensive set of code generators for various floating-point units (FPUs) essential for integrated circuit design and integration, targeting a broad spectrum of applications, including machine learning and scientific computing. The suite includes FP adders, multipliers, subtractors, dividers, reciprocals, exponentials, square roots, trigonometric functions (sine, cosine, arctangent), and more. It supports customizable hardware design parameters, such as precision (16, 32, 64, and 128 bits) and pipeline depths, offering users enhanced flexibility and productivity. The generated code is in an industry-standard hardware description language, ensuring compatibility with standard design flows, including simulation, verification, synthesis, and implementation on both field-programmable gate arrays (FPGAs) and application-specific integrated circuits (ASICs).

Shalf, JohnM. [Lawrence Berkeley National Laborato↗

MGGPOD: a Monte Carlo Suite for Modeling Instrumental Line and Continuum Backgrounds in Gamma-Ray Astronomy

Intense and complex instrumental backgrounds, against which the much smaller signals from celestial sources have to be discerned, are a notorious problem for low and intermediate energy gamma-ray astronomy (approximately 50 keV - 10 MeV). Therefore a detailed qualitative and quantitative understanding of instrumental line and continuum backgrounds is crucial for most stages of gamma-ray astronomy missions, ranging from the design and development of new instrumentation through performance prediction to data reduction. We have developed MGGPOD, a user-friendly suite of Monte Carlo codes built around the widely used GEANT (Version 3.21) package, to simulate ab initio the physical processes relevant for the production of instrumental backgrounds. These include the build-up and delayed decay of radioactive isotopes as well as the prompt de-excitation of excited nuclei, both of which give rise to a plethora of instrumental gamma-ray background lines in addition t o continuum backgrounds. The MGGPOD package and documentation are publicly available for download. We demonstrate the capabilities of the MGGPOD suite by modeling high resolution gamma-ray spectra recorded by the Transient Gamma-Ray Spectrometer (TGRS) on board Wind during 1995. The TGRS is a Ge spectrometer operating in the 40 keV to 8 MeV range. Due to its fine energy resolution, these spectra reveal the complex instrumental background in formidable detail, particularly the many prompt and delayed gamma-ray lines. We evaluate the successes and failures of the MGGPOD package in reproducing TGRS data, and provide identifications for the numerous instrumental lines.

Weidenspointner, G.↗

Multi-scale fission product release model with comparison to AGR data

TRistructural ISOtropic (TRISO) particle fuel is central to several advanced, high-temperature reactor designs. Each particle consists of a fuel kernel encapsulated by three layers of carbon and ceramics that prevent the release of fission products and ensure physical integrity. Despite outstanding retention properties, fission product release has been observed from intact particles. To better understand and quantify fission product release from TRISO particles, a multiscale, mechanistic model of fission product transport is being developed by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. Previous work focused on silver (Ag) transport and improved Ag release predictions. The work described in this report builds on this experience to better understand cesium (Cs) transport in silicon carbide (SiC), the main barrier to the release of fission products. Atomistic simulations provide bulk and grain boundary (GB) Cs diffusivities in SiC, which are used by phase field simulations in the mesoscale code Marmot to determine the temperature, microstructure, and irradiation-dependent Cs diffusivity at the mesoscale in SiC. This approach attributes the different temperature regimes experimentally observed for Cs diffusivities in SiC to a transition from bulk-dominated diffusivity at high temperatures to a GB-dominated regime at low temperatures, providing new insight. The multiscale, mechanistic effective diffusivity is then implemented in the fuel performance code BISON and further validated by comparing Cs release predictions from Advanced Gas Reactor (AGR)-1 and AGR-2 post-irradiation measurements. The new model improves BISON’s predictability. This document also reports improvements made on Ag transport modeling by accounting for different GB types having different diffusivities. Moreover, this report details preliminary efforts to model palladium (Pd) attack of the SiC at the mesoscale using a phase field approach. Pd attack and its impact on accelerated Ag transport remains a misunderstood phenomenon, and we use the model to demonstrate that the formation of lamellae that has been observed in experiments can be explained by the reaction of Pd with SiC to form alternating layers of graphite and Pd 2 Si. This effort aims to improve our understanding of the reaction and eventually provide a model for BISON to account for Pd penetration and its effects on fission product release.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Non-adiabatic Excited State Molecule Dynamics Modeling of Photochemistry and Photophysics of Materials

Understanding and controlling excited state dynamics lies at the heart of all our efforts to design photoactive materials with desired functionality. This tailor-design approach has become the standard for many technological applications (e.g., solar energy harvesting) including the design of organic conjugated electronic materials with applications in photovoltaic and light-emitting devices. Over the years, our team has developed efficient LANL-based codes to model the relevant photophysical processes following photoexcitation (spatial energy transfer, excitation localization/delocalization, and/or charge separation). The developed approach allows the non-radiative relaxation to be followed on up to ~10 ps timescales for large realistic molecules (hundreds of atoms in size) in the realistic solvent dielectric environment. The Collective Electronic Oscillator (CEO) code is used to compute electronic excited states, and the Non-adiabatic Excited State Molecular Dynamics (NA-ESMD) code is used to follow the non-adiabatic dynamics on multiple coupled Born-Oppenheimer potential energy surfaces. Our preliminary NA-ESMD simulations have revealed key photoinduced mechanisms controlling competing interactions and relaxation pathways in complex materials, including organic conjugated polymer materials, and have provided a detailed understanding of photochemical products and intermediates and the internal conversion process during the initiation of energetic materials. This project will be using LANL-based CEO and NA-ESMD codes to model nonradiative relaxation in organic and energetic materials. The NA-ESMD and CEO codes belong to a class of electronic structure/quantum chemistry codes that require large memory, “long-queue-few-core” distribution of resources in order to make useful progress. The NA-ESMD simulations are trivially parallelizable requiring ~300 processors for up to one week runtime to reach a meaningful restart point.

36 MATERIALS SCIENCE↗

Modelling gamma-ray photon emission and pair production in high-intensity laser–matter interactions

In high-intensity (>10{sup 21} Wcm{sup −2}) laser–matter interactions gamma-ray photon emission by the electrons can strongly affect the electron's dynamics and copious numbers of electron–positron pairs can be produced by the emitted photons. We show how these processes can be included in simulations by coupling a Monte Carlo algorithm describing the emission to a particle-in-cell code. The Monte Carlo algorithm includes quantum corrections to the photon emission, which we show must be included if the pair production rate is to be correctly determined. The accuracy, convergence and energy conservation properties of the Monte Carlo algorithm are analysed in simple test problems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Using Polarized PCRTM-Solar Model to Validate EMIT Surface Reflectance L2A Product

Using the vector linearized discrete ordinate radiative transfer model (VLIDORT), a multiple-scattering multi-layer scattering code, coupled with a fast principal component (PC) radiative transfer model (PCRTM), we can simulate hyperspectral radiation fields from visible to shortwave infrared (VSWIR), i.e., 250 nm to 2500 nm. Since only limited accurate line-by-line calculation using VLIDORT is needed, this coupled PCRTM-VLIDORT approach is very efficient for simulating hyperspectral VSWIR spectra. These simulations incorporate accurately defined optical properties (gas absorptions, molecular and aerosol/cloud particle scattering) and detailed surface reflectance inputs. In this study we used the current L2A products of the Earth Surface Mineral Dust Source Investigation (EMIT) imaging spectrometer as the inputs to PCRTM-VLIDORT to simulate EMIT spectra. Comparing the simulated spectra with EMIT measurements provides an efficient method to validate the current EMIT product.

Xiaozhen Xiong↗

Validation of nuclear data based on the ASPIS experiments from the SINBAD database - 25428

SINBAD, the Shielding Integral Benchmark Archive and Database, includes the experimental data relative to integral benchmark experiments relevant for radiation shielding applications. The ASPIS experiments, included in the Reactor Shielding section of SINBAD, provide neutron activation foils and spectra measurements performed before/between shielding slabs of Graphite, Sodium, Iron, Steels, Water and combinations of them. Benchmark models for the measurement simulations with the MCNP code are usually provided in SINBAD. Nonetheless, a user-friendly tool is presented, which would be useful for an independent production of MCNP benchmark models. It is based on GPL C++ code called CombLayer. The reports issued for SINBAD on the quality assessment of the ASPIS experiments are also described. The discussion on the quality of the ASPIS experiments for their use in nuclear data evaluations had a positive outcome, except for two experiments out of ten. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) - NIMROD DIII-D Shattered Pellet Injection Disruption Mitigation and RE Modeling. Final report

Final technical report for UCSD and subcontractor contributions to DE-SC0016420, Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM). Project description. MHD modeling includeing runaway electron drift orbit calculations will be carried out for DII-D disruption mitigation scenarios. The effects of pitch angle scattering on the energy distribution of runaway electron population will be evaluated using one or more collision models integrated / coupled with the runaway orbit following model in the nimrod code. An extension of the test particle model will be developed to predict runaway electron generation considering both primary and secondary sources of runaway electron production. Simulations simultaneously modelling the effects of generation, confinement, and moderation-space losses during a DII-D rapid shutdown (mitigated disruption) will be completed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

General relativistic simulations of collapsing binary neutron star mergers with Monte Carlo neutrino transport

Recent gravitational wave observations of neutron-star-neutron-star and neutron-star-black-hole binaries appear to indicate that massive neutron stars may not be too uncommon in merging systems. These discoveries have led to an increased interest in the simulation of merging compact binaries involving massive stars. In this paper, we present a first set of evolution of massive neutron star binaries using Monte Carlo radiation transport for the evolution of neutrinos. Here, we study a range of systems, from nearly symmetric binaries that collapse to a black hole before forming a disk or ejecting material, to more asymmetric binaries in which tidal disruption of the lower mass star leads to the production of more interesting postmerger remnants. For the latter type of systems, we additionally study the impact of viscosity on the properties of the outflows, and compare our results to two recent simulations of identical binaries performed with the whiskyTHC code. We find agreement on the black hole properties, disk mass, and mass and velocity of the outflows within expected numerical uncertainties, and some minor but noticeable differences in the evolution of the electron fraction when using a subgrid viscosity model, with viscosity playing a more minor role in our simulations. The method used to account for r-process heating in the determination of the outflow properties appears to have a larger impact on our result than those differences between numerical codes. We also use the simulation with the most ejected material to verify that our newly implemented Lagrangian tracers provide a reasonable sampling of the matter outflows as they leave the computational grid. We note that, given the lack of production of hot outflows in these mergers, the main role of neutrinos in these systems is to set the composition of the postmerger remnant. One of the main potential uses of our simulations is, thus, as improved initial conditions for longer evolutions of such remnants.

79 ASTRONOMY AND ASTROPHYSICS↗

Then and Now: Improving Software Portability, Productivity, and 100× Performance

The US Exascale Computing Project (ECP) has succeeded in preparing applications to run efficiently on the first reported Exascale supercomputers in the world. To achieve this, it modernized the whole leadership software stack, from libraries to simulation codes. In this article, we contrast selected leadership software before and after ECP. We discuss how sustainable research software development for leadership computing can embrace the conversation with the hardware vendors, the leadership computing facilities, the software community, and the domain scientists who are the application developers and integrators of software products. We elaborate on how software needs to take portability as a central design principle and to benefit from interdependent teams; we also demonstrate how moving to programming languages with high momentum, like modern C++, can help improve the sustainability, interoperability, and performance of research software. Finally, we showcase how cross-institutional efforts can enable algorithm advances that are beyond incremental performance optimization.

97 MATHEMATICS AND COMPUTING↗