Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer Codes Computer Programs”

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,711 records · Page 95

Numerical simulation projects in micromagnetics with Jupyter

We report a case study where an existing materials science course was modified to include numerical simulation projects on the micromagnetic behavior of materials. The Ubermag micromagnetic simulation software package is used in order to solve problems computationally. The simulation software is controlled through the Python code in Jupyter notebooks. Our experience is that the self-paced problem-solving nature of the project work can facilitate a better in-depth exploration of the course contents. We discuss which aspects of the Ubermag and the project Jupyter ecosystem have been beneficial for the students' learning experience and which could be transferred to similar teaching activities in other subject areas.

97 MATHEMATICS AND COMPUTING↗

COFFEA - Columnar Object Framework For Effective Analysis [Slides]

The COFFEA Framework provides a new approach to HEP analysis, via columnar operations, that improves time-to-insight, scalability, portability, and reproducibility of analysis. It is implemented with the Python programming language and commodity big data technologies such as Apache Spark and NoSQL databases. To achieve this suite of improvements across many use cases, COFFEA takes a factorized approach, separating the analysis implementation and data delivery scheme. All analysis operations are implemented using the NumPy or awkward-array packages which are wrapped to yield user code whose purpose is quickly intuited. Various data delivery schemes are wrapped into a common front-end which accepts user inputs and code, and returns user defined outputs. We will present published results from analysis of CMS data using the COFFEA framework along with a discussion of metrics and the user experience of arriving at those results with columnar analysis.

97 MATHEMATICS AND COMPUTING↗

Modeling Enhancements, Cross-Section Generation Updates, and Benchmarking with Shift

This technical report documents the modeling enhancements, cross-section generation updates, and bench marking with the Shift Monte Carlo code performed under the US Department of Energy Nuclear Energy Advanced Modeling and Simulation Program in FY 2024. The work performed included several modeling enhancements, such as integration of cross-section generation in Titan and the ability to produce microscopic multigroup cross sections with Shift. Benchmarking of the cross sections produced by Shift and the two-step workflow with Griffin was performed for three problems: the Advanced Breeder Test Reactor, a generic pebble bed reactor, and a TRISO heat pipe microreactor. Comparisons of results from these benchmark problems were done with Serpent, OpenMC, and Griffin. These enhancements provide a robust foundation for applying Shift for both reference and two-step neutronics analysis for advanced reactor simulation.

97 MATHEMATICS AND COMPUTING↗

Data Analyst - A Device for Continuous Repeating Real-Time Assay of HPGe, CZT, or Scintillation Gamma Spectra - 18574

There are many situations where a continuous repeating series of quantitative gamma spectral assays is desirable. These include measurements of fluid going through a pipe, of items on a conveyor passing under a detector, a detector aimed at the ground on a moving platform, or a detector in free air as a plume of radioactivity drifts by. While Canberra and many others have built and deployed many of these devices, they have all been custom built for that particular application, and rather expensive. The Data Analyst [DA] is a general purpose device that is designed to fulfill many of these applications without modification, but can be fairly-easily modified if needed. The DA is a small [12 x 18 x 5 cm], low-power [3 W] box that accepts the input from external MCAs, controls the acquisition cycle, analyzes the results, and transmits them to the outside world when needed. The detector-MCAs supported are HPGe with Lynx MCA, NaI or CeBr or LaBr with Osprey MCA, and CZT or CsI with internal MCAs. The DA internal software used for analysis is the full standard Genie analysis software, which is very widely used and well accepted in the industry. This allows exported spectral files to be reviewed and confirmed by Experts, and reanalyzed if needed. The DA hardware is totally autonomous. As soon as power is applied, it starts the pre-programmed sequences of acquisition and analysis. Multiple analysis sequences [called Workflows] can be in operation simultaneously, e.g. with different count times and different nuclide libraries. Only the detector with MCA, and the DA are required in the field. A PC can be connected to the DA locally or remotely via WiFi, USB, or Ethernet. The PC is needed to setup the DA, and to readout the results. The DA also supports [with application-specific programming] external Inputs/Outputs either for recording with the analysis results or for use in the computation of the result, e.g. temperature, flow rate, pressure, ... ]. The DA also has an internal GPS receiver with external antenna. Applications of the DA that are discussed include stack gas effluent activity where a portion of the effluent stream was extracted to flow through a shield with a HPGe detector and Lynx MCA, and measurement of reactor coolant activity where a shielded CZT detector and MCA were aimed at various sections of reactor primary coolant during a maintenance outage. Potential other applications are also described. (author)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

GNDS v2.0 release and future developments [Slides]

This presentation begins by discussing the major changes in v2.0 of GNDS and that the GNDS-2.0 formatted libraries are now available. It also examines how code support for GNDS is growing with processing codes, APIs and transport codes, and evaluation/testing codes. The presentation states that while migrating to GNDS can be a major effort it worthwhile, especially for anyone who is just getting started as nuclear data creator or user. Additionally, this presentation outlines the established milestones that the EG-GNDS (the group in charge of GNDS specs) has for the next GNDS release and elaborates on the first goal, which is to migrate GNDS specifications to the new JSON schema language.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Streamlining resummed QCD calculations using Monte Carlo integration

Some of the most arduous and error-prone aspects of precision resummed calculations are related to the partonic hard process, having nothing to do with the resummation. In particular, interfacing to parton-distribution functions, combining various channels, and performing the phase space integration can be limiting factors in completing calculations. Conveniently, however, most of these tasks are already automated in many Monte Carlo programs, such as MadGraph [1], Alpgen [2] or Sherpa [3]. In this paper, we show how such programs can be used to produce distributions of partonic kinematics with associated color structures representing the hard factor in a resummed distribution. These distributions can then be used to weight convolutions of jet, soft and beam functions producing a complete resummed calculation. In fact, only around 1000 unweighted events are necessary to produce precise distributions. A number of examples and checks are provided, including e + e – two- and four-jet event shapes, n -jettiness and jet-mass related observables at hadron colliders at next-to-leading-log (NLL) matched to leading order (LO). Furthermore, the attached code can be used to modify MadGraph to export the relevant LO hard functions and color structures for arbitrary processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Magnetism and transport properties of layered rare-earth cobaltates Ln{sub 0.3}CoO{sub 2}

The ab-initio (GGA+U) electronic structure calculations of layered cobaltates Ln{sub 0.3}CoO{sub 2} (Ln = La, Pr, Nd) prepared by ionic exchange from Na{sub 0.90}CoO{sub 2} precursors have been performed. The data are used for numerical modeling of Seebeck coefficient within Boltzmann transport theory using BoltzTraP program [G. K. H. Madsen and D. J. Singh, Comput. Phys. Commun. 175, 67 (2006)], as well as for determination of the crystal field split levels of rare-earth ions using a method based on a transformation of Bloch states into the basis of Wannier functions [P. Novák et al., Phys. Rev. B 87, 205139 (2013)]. An overall agreement with observed magnetism and transport properties is obtained. In particular, the high p-type thermopower is well reproduced in a broad temperature range, but instead of theoretical linear decrease down to the lowest temperatures, the real systems exhibit an anomalous change of Seebeck sign, which might be related to the change of bare metallic carriers into the polaronic ones.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Argonne National Laboratory Site Environmental Report for Calendar Year 2023

This report discusses the status and the accomplishments of the environmental protection program at Argonne National Laboratory for calendar year 2023. The status of Argonne environmental protection activities with respect to compliance with the various laws and regulations is discussed, along with environmental management, sustainability efforts, environmental corrective actions, and habitat restoration. To evaluate the effects of Argonne operations on the environment, samples of environmental media collected on the site, at the site boundary, and off the Argonne site were analyzed and compared with applicable guidelines and standards. A variety of radionuclides were measured in air, surface water, groundwater, and bottom sediment samples. In addition, chemical constituents in surface water, groundwater, and wastewater were analyzed. External penetrating radiation doses were measured, and the potential for radiation exposure to off site population groups was estimated. Results are interpreted with respect to the origin of the radioactive and chemical substances (i.e., natural, Argonne, and other) and are compared with applicable standards intended to protect human health and the environment. A U.S. Department of Energy (DOE) dose calculation methodology, based on International Commission on Radiological Protection (ICRP) recommendations and the U.S. Environmental Protection Agency’s (EPA) CAP 88 computer code, was used in preparing this report.

54 ENVIRONMENTAL SCIENCES↗

Implementation of Perturbation Theory and Sensitivity Capabilities in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor Multiphysics analysis application, jointly developed by Argonne and Idaho National Laboratories under the DOE-NE NEAMS program. This fiscal year, capabilities for reactivity and sensitivity evaluation using perturbation methods were implemented and verified. The First Order Perturbation Method (FOPT) was employed to compute reactivity worth resulting from small perturbations in input parameters, while the Generalized Perturbation Theory (GPT) was used to evaluate sensitivities of a range of response types, including reaction rate ratio, k-eigenvalue, neutron generation time, and effective delayed neutron fraction. These perturbation methods enable users to quantify how response quantities change due to a perturbation in a input parameter without explicitly performing an additional transport simulation for each perturbed state. In particular, the GPT formulation accounts for indirect effects arising from flux changes by solving generalized inhomogeneous equations, for which a Neumann series-based iterative solution method was developed and implemented in Griffin. The implemented reactivity and sensitivity evaluation capabilities were verified using two test problems: an infinite homogeneous system and a two-dimensional hexagonal core. The results showed excellent agreement with reference solutions obtained by a direct method based on finite difference approximation as well as GPT-based results from the PERSENT code, confirming the accuracy of both reactivity and sensitivity evaluations. Additionally, preliminary uncertainty quantification (UQ) results were obtained by combining the sensitivity values computed using GPT and external covariance data, demonstrating that the implemented sensitivity results can be reliably used for uncertainty calculations. To further demonstrate the generality and practical strength of the implementation, the sensitivity evaluation capability was successfully applied to the Empire microreactor with a geometrically complex design that poses significant modeling challenges. The results confirm that Griffin enables sensitivity evaluations even for irregular and highly heterogeneous reactor configurations, thereby establishing a foundation for UQ applications in advanced reactor designs and analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Generic Biogeochemical Module for Earth System Models: Next Generation BioGeoChemical Module (NGBGC), Version 1.0

Physical and biogeochemical processes regulate soil carbon dynamics and CO2 flux to and from atmosphere, influencing global climate changes. Integration of these processes into earth system models (e.g., community land models (CLM)), however, currently faces three major challenges: 1) extensive efforts are required to modify modeling structures and to rewrite computer programs to incorporate new or updated processes as new knowledge is being generated, 2) computational cost is prohibitively expensive to simulate biogeochemical processes in land models due to large variations in the rates of biogeochemical processes, and 3) various mathematical representations of biogeochemical processes exist to incorporate different aspects of fundamental mechanisms, but systematic evaluation of the different mathematical representations is difficult, if not possible. To address these challenges, we propose a new computational framework to easily incorporate physical and biogeochemical processes into land models. The new framework consists of a new biogeochemical module with a generic algorithm and reaction database so that new and updated processes can be incorporated into land models without the need to manually set up the ordinary differential equations to be solved numerically. The reaction database consists of processes of nutrient flow through the terrestrial ecosystems in plants, litter and soil. This framework facilitates effective comparison studies of biogeochemical cycles in an ecosystem using different conceptual models under the same land modeling framework. The approach was first implemented in CLM and benchmarked against simulations from the original CLM-CN code. A case study was then provided to demonstrate the advantages of using the new approach to incorporate a phosphorus cycle into the CLM model. To our knowledge, the phosphorus-incorporated CLM is a new model that can be used to simulate phosphorus limitation on the productivity of terrestrial ecosystems.

Fang, Yilin↗

3D simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C1, and NIMROD developed the capability of modeling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state-of-the-art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. There are differences between the physics models and numerical methods, and the VDE evolution leads to sensitivities on the initial conditions that cannot be avoided as can be done in edge localized modes (ELM) and sawtooth simulations (due to the non-cyclical nature of VDEs). Nonetheless, the comparison serves to quantify the level of agreement in the relevant quantities used to characterize disruptions, such as the 3D wall forces and energy decay. The results bring confidence regarding the use of the mentioned codes for disruption studies, and they distinguish aspects that are specific to the models used (e.g., reduced vs full MHD models). Here, the simulations show important 3D fea- tures for a NSTX plasma, such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz-based JOREK reduced MHD model allows capturing many aspects of the 3D dynamics even in the spherical tokamak limit considered here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SU-E-T-99: Design and Development of Isocenter Parameter System for CT Simulation Laser Based On DICOM RT

Purpose: In order to receive DICOM files from treatment planning system and generate patient isocenter positioning parameter file for CT laser system automatically, this paper presents a method for communication with treatment planning system and calculation of isocenter parameter for each radiation field. Methods: Coordinate transformation and laser positioning file formats were analyzed, isocenter parameter was calculated via data from DICOM CT Data and DICOM RTPLAN file. An in-house software-DicomGenie was developed based on the object-oriented program platform-Qt with DCMTK SDK (Germany OFFIS company DICOM SDK) . DicomGenie was tested for accuracy using Philips CT simulation plan system (Tumor LOC, Philips) and A2J CT positioning laser system (Thorigny Sur Marne, France). Results: DicomGenie successfully established DICOM communication between treatment planning system, DICOM files were received by DicomGenie and patient laser isocenter information was generated accurately. Patient laser parameter data files can be used for for CT laser system directly. Conclusion: In-house software DicomGenie received and extracted DICOM data, isocenter laser positioning data files were created by DicomGenie and can be use for A2J laser positioning system.

60 APPLIED LIFE SCIENCES↗

Design of 95 GHz gyrotron based on continuous operation copper solenoid with water cooling

The design work for 2nd harmonic 95 GHz, 50 kW gyrotron based on continuous operation copper solenoid is presented. Thermionic magnetron injection gun specifications were calculated according to the linear trade off equation, and simulated with CST program. Numerical code is used for cavity design using the non-uniform string equation as well as particle motion in the “cold” cavity field. The mode TE02 with low Ohmic losses in the cavity walls was chosen as the operating mode. The Solenoid is designed to induce magnetic field of 1.8 T over a length of 40 mm in the interaction region with homogeneity of ±0.34%. The solenoid has six concentric cylindrical segments (and two correction segments) of copper foil windings separated by water channels for cooling. The predicted temperature in continuous operation is below 93 °C. The parameters of the design together with simulation results of the electromagnetic cavity field, magnetic field, electron trajectories, and thermal analyses are presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

ACDT: Architected Composite Data Types Trading-in Unfettered Data Access for Improved Execution

With Exascale performance and its challenges in mind, one ubiquitous concern among architects is energy efficiency. Petascale systems projected to Exascale systems are unsustainable at current power consumption rates. One major contributor to system-wide power consumption is the number of memory operations leading to data movement and management techniques applied by the runtime system. To address this problem, we present the concept of the Architected Composite Data Types (ACDT) framework. The framework is made aware of data composites, assigning them a specific layout, transformations and operators. Data manipulation overhead is amortized over a larger number of elements and program performance and power efficiency can be significantly improved. We developed the fundamentals of an ACDT framework on a massively multithreaded adaptive runtime system geared towards Exascale clusters. Showcasing the capability of ACDT, we exercised the framework with two representative processing kernels - Matrix Vector Multiply and the Cholesky Decomposition – applied to sparse matrices. As transformation modules, we applied optimized compress/decompress engines and configured invariant operators for maximum energy/performance efficiency. Additionally, we explored two different approaches based on transformation opaqueness in relation to the application. Under the first approach, the application is agnostic to compression and decompression activity. Such approach entails minimal changes to the original application code, but leaves out potential applicationspecific optimizations. The second approach exposes the decompression process to the application, hereby exposing optimization opportunities that can only be exploited with application knowledge. The experimental results show that the two approaches have their strengths in HW and SW respectively, where the SW approach can yield performance and power improvements that are an order of magnitude better than ACDT-oblivious, hand-optimized implementations.We consider the ACDT runtime framework an important component of compute nodes that will lead towards power efficient Exascale clusters.

exascale↗

Advanced Computational Framework for Environmental Management ZEM, Version 1.x

Typically environmental management problems require analysis of large and complex data sets originating from concurrent data streams with different data collection frequencies and pedigree. These big data sets require on-the-fly integration into a series of models with different complexity for various types of model analyses where the data are applied as soft and hard model constraints. This is needed to provide fast iterative model analyses based on the latest available data to guide decision-making. Furthermore, the data and model are associated with uncertainties. The uncertainties are probabilistic (e.g. measurement errors) and non-probabilistic (unknowns, e.g. alternative conceptual models characterizing site conditions). To address all of these issues, we have developed an integrated framework for real-time data and model analyses for environmental decision-making called ZEM. The framework allows for seamless and on-the-fly integration of data and modeling results for robust and scientifically-defensible decision-making applying advanced decision analyses tools such as Bayesian- Information-Gap Decision Theory (BIG-DT). The framework also includes advanced methods for optimization that are capable of dealing with a large number of unknown model parameters, and surrogate (reduced order) modeling capabilities based on support vector regression techniques. The framework is coded in Julia, a state-of-the-art high-performance programing language (http://julialang.org). The ZEM framework is open-source and can be applied to any environmental management site. The framework will be open-source and released under GPL V3 license.

Vesselinov, VelimirV.↗

Open-source release of CGMF 1.1 and Integration into the MCNP6.3 ® Code [Slides]

As a result of a multi-year NA-22 project, CGMF was integrated into MCNP6.2 and publicly released. CGMF was open-sourced and publicly released and MCNP6.3 was updated to include the latest version and is in the process of being publicly released. Current and future plans include global optimization and uncertainty quantification within CGMF, model parameter fitting such that CGMF may be used in ENDF/B evaluations, and improving both standalone and MCNP-integrated CGM (non-fission) simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Methods, media, and systems for detecting attack on a digital processing device

Methods, media, and systems for detecting attack are provided. In some embodiments, the methods include: comparing at least part of a document to a static detection model; determining whether attacking code is included in the document based on the comparison of the document to the static detection model; executing at least part of the document; determining whether attacking code is included in the document based on the execution of the at least part of the document; and if attacking code is determined to be included in the document based on at least one of the comparison of the document to the static detection model and the execution of the at least part of the document, reporting the presence of an attack. In some embodiments, the methods include: selecting a data segment in at least one portion of an electronic document; determining whether the arbitrarily selected data segment can be altered without causing the electronic document to result in an error when processed by a corresponding program; in response to determining that the arbitrarily selected data segment can be altered, arbitrarily altering the data segment in the at least one portion of the electronic document to produce an altered electronic document; and determining whether the corresponding program produces an error state when the altered electronic document is processed by the corresponding program.

97 MATHEMATICS AND COMPUTING↗

Methods, media, and systems for detecting attack on a digital processing device

Methods, media, and systems for detecting attack are provided. In some embodiments, the methods include: comparing at least part of a document to a static detection model; determining whether attacking code is included in the document based on the comparison of the document to the static detection model; executing at least part of the document; determining whether attacking code is included in the document based on the execution of the at least part of the document; and if attacking code is determined to be included in the document based on at least one of the comparison of the document to the static detection model and the execution of the at least part of the document, reporting the presence of an attack. In some embodiments, the methods include: selecting a data segment in at least one portion of an electronic document; determining whether the arbitrarily selected data segment can be altered without causing the electronic document to result in an error when processed by a corresponding program; in response to determining that the arbitrarily selected data segment can be altered, arbitrarily altering the data segment in the at least one portion of the electronic document to produce an altered electronic document; and determining whether the corresponding program produces an error state when the altered electronic document is processed by the corresponding program.

97 MATHEMATICS AND COMPUTING↗