SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,765 records · Page 98

All sources NASA NTRS DOE OSTI

DOE OSTI2018

SPALLOOP

Boiler tubes in steam power plants experience tube blockages due to exfoliation of oxide grown on the inner side of the tubes. In extreme cases, significant tube blockages can lead to forced power plant outages. It is thus desired to predict through modeling the amount of tube blockage in order to inform power plant operators of possible forced outages. SpalLoop solves for the stress-strain equations in an axisymmetric geometry, tracking the stress/strain evolution during boiler operation including outages for the entire boiler tube length. At each operational outage, i.e., temperature excursions down to room temperature, the amount of exfoliated area for the entire tube loop is estimated the amount of tube blockage is predicted based assumed blockage geometry and site. The SpaLLoop code contains modules developed for oxide growth, stress analysis, tube loop geometry, blockage area by taking into account the following phenomena and features, (a) Plant operation schedule with periodic alternate full-load and partial-load regimes and shut-downs, i.e., temperature excursions from high-load to room temperature, (b) axisymmetric formulation for cylindrical tubes, (c) oxide growth in a temperature gradient with multiple oxide layers, (d) geometry of a boiler tube with a single tube loop or two tube loops, (e) temperature variation along the tube length based on hot gas temperature distribution outside the tube and inlet steam temperature, (f) non-uniform oxide growth along the tube length according to the local steam tube temperature, (g) exfoliated area module: at each operational outage considered, the amount of exfoliated area and exfoliated volume along the tube is estimated, (h) blockage module: at each operational outage considered, the exfoliated volume/mass for each tube loop is estimated from which the amount of tube blockage is predicted based on given blockage geometry (length, location, and geometry). The computer program is written in FORTRAN90. Its modular structure was sought for allowing the best flexibility in updating the program by implementing new constitutive equations due to availability of new material property data and/or new physical phenomena.

Sabau, Adrian↗

DOE OSTI2018

Hybrid simulations of fishbone instabilities and Alfvén eigenmodes in DIII-D tokamak

In DIII-D hybrid discharges, the intense Alfvén eigenmode (AE) activity driven by Neutral Beam Injection (NBI) that is typically observed can be suppressed and replaced by fishbone modes when Electron Cyclotron Current Drive (ECCD) is centrally applied. Simulations have been carried out with the kinetic/magnetohydrodynamic (MHD) hybrid code M3D-K based on DIII-D discharges #161401 and #161403 without and with ECCD, respectively. In both cases unstable modes are found – the mode frequency and mode structure indicate the instability excited in #161403 is of fishbone type, while that in #161401 is identified as the beta-induced Alfvén eigenmode (BAE)-like mode. Moreover, we find that the calculated mode frequencies of these two shots are consistent with experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2011

A PENELOPE-based system for the automated Monte Carlo simulation of clinacs and voxelized geometries - application to far-from-axis fields

Purpose: Two new codes, PENEASY and PENEASYLINAC, which automate the Monte Carlo simulation of Varian Clinacs of the 600, 1800, 2100, and 2300 series, together with their electron applicators and multileaf collimators, are introduced. The challenging case of a relatively small and far-from-axis field has been studied with these tools. Methods: PENEASY is a modular, general-purpose main program for the PENELOPE Monte Carlo system that includes various source models, tallies and variance-reduction techniques (VRT). The code includes a new geometry model that allows the superposition of voxels and objects limited by quadric surfaces. A variant of the VRT known as particle splitting, called fan splitting, is also introduced. PENEASYLINAC, in turn, automatically generates detailed geometry and configuration files to simulate linacs with PENEASY. These tools are applied to the generation of phase-space files, and of the corresponding absorbed dose distributions in water, for two 6 MV photon beams from a Varian Clinac 2100 C/D: a 40 x 40 cm{sup 2} centered field; and a 3 x 5 cm{sup 2} field centered at (4.5, -11.5) cm from the beam central axis. This latter configuration implies the largest possible over-traveling values of two of the jaws. Simulation results for the depth dose and lateral profiles at various depths are compared, by using the gamma index, with experimental values obtained with a PTW 31002 ionization chamber. The contribution of several VRTs to the computing speed of the more demanding off-axis case is analyzed. Results: For the 40 x 40 cm{sup 2} field, the percentages {gamma}{sub 1} and {gamma}{sub 1.2} of voxels with gamma indices (using 0.2 cm and 2% criteria) larger than unity and larger than 1.2 are 0.2% and 0%, respectively. For the 3 x 5 cm{sup 2} field, {gamma}{sub 1} = 0%. These figures indicate an excellent agreement between simulation and experiment. The dose distribution for the off-axis case with voxels of 2.5 x 2.5 x 2.5 mm{sup 3} and an average standard statistical uncertainty of 2% (1{sigma}) is computed in 3.1 h on a single core of a 2.8 GHz Intel Core 2 Duo processor. This result is obtained with the optimal combination of the tested VRTs. In particular, fan splitting for the off-axis case accelerates execution by a factor of 240 with respect to standard particle splitting. Conclusions: PENEASY and PENEASYLINAC can simulate the considered Varian Clinacs both in an accurate and efficient manner. Fan splitting is crucial to achieve simulation results for the off-axis field in an affordable amount of CPU time. Work to include Elekta linacs and to develop a graphical interface that will facilitate user input is underway.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2024

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↗

DOE OSTI2019

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↗

DOE OSTI2024

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↗

DOE OSTI2018

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↗

DOE OSTI2022

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↗

DOE OSTI2016

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↗

DOE OSTI2015

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↗

DOE OSTI2024

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↗

DOE OSTI2025

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↗

DOE OSTI2013

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↗

DOE OSTI2021

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↗

DOE OSTI2014

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↗

DOE OSTI2014

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 ↗

DOE OSTI2014

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↗

DOE OSTI2016

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.↗