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,261 records · Page 70

Synergizing human expertise and AI efficiency with language model for microscopy operation and automated experiment design

With the advent of large language models (LLMs), in both the open source and proprietary domains, attention is turning to how to exploit such artificial intelligence (AI) systems in assisting complex scientific tasks, such as material synthesis, characterization, analysis and discovery. Here, we explore the utility of LLMs, particularly ChatGPT4, in combination with application program interfaces (APIs) in tasks of experimental design, programming workflows, and data analysis in scanning probe microscopy, using both in-house developed APIs and APIs given by a commercial vendor for instrument control. We find that the LLM can be especially useful in converting ideations of experimental workflows to executable code on microscope APIs. Beyond code generation, we find that the GPT4 is capable of analyzing microscopy images in a generic sense. At the same time, we find that GPT4 suffers from an inability to extend beyond basic analyses for more in-depth technical experimental design. We argue that an LLM specifically fine-tuned for individual scientific domains can potentially be a better language interface for converting scientific ideations from human experts to executable workflows. Such a synergy between human expertise and LLM efficiency in experimentation can open new doors for accelerating scientific research, enabling effective experimental protocols sharing in the scientific community.

97 MATHEMATICS AND COMPUTING↗

Generalized Bayesian Framework for Evaluation of Integral Benchmark Experiments [Slides]

This presentation is on generalized Bayesian framework for evaluation of integral benchmark experiments. This presentation starts off with assumptions and approximations used with Bayes Theorem. then an overview of approximations used by ORNL codes, and Generalized Bayesian Monto Carlo (GBMC). The presentation then details out a precise framework. This presentation then concludes with considerations.

97 MATHEMATICS AND COMPUTING↗

INVESTIGATION OF BREAKDOWN INDUCED SURFACE DAMAGE ON 805 MHZ PILLBOX CAVITY INTERIOR SURFACES

The MuCool Test Area (MTA) at Fermilab is a facility to develop the technology required for ionization cooling for a future Muon Collider and/or Neutrino Factory. As part of this research program, we have tested two 805 MHz vacuum RF cavities in a multi-Tesla magnetic field to study the effects of the static magnetic field on the cavity operation. This study gives useful information on field emitters in the cavity, dark current, surface conditioning, breakdown mechanisms and material properties of the cavity. All these factors determine the maximum accelerating gradient in the cavity. This paper discusses the image processing technique for quantitative estimation of spark damage spot distribution on cavity interior surfaces. The distribution is compared with the electric field distribution predicted by a computer code calculation. The local spark density is proportional to probability of surface breakdown and shows a power law dependence on the maximum electric field (E). This E dependence is consistent with the dark current calculated from the Fowler-Nordheim equation.

43 PARTICLE ACCELERATORS↗

Simulation of the RPI Reactor Critical Facility Core Using High-Fidelity Neutronics Code PROTEUS

The PROTEUS code is a high-fidelity three-dimensional (3D) deterministic neutron transport code developed by ANL under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The code contains the SN2ND (second order discrete ordinates method) and two MOC transport solvers, MOCFE and MOCEX (first order methods of characteristics), based on unstructured finite element meshes which allow users to model complex or unconventional geometry reactor problems. Among the two MOC solvers which normally perform better than SN for heterogeneous geometry problems, MOCEX is based on a rigorous formulation using the two-dimensional (2D) MOC radially and the discontinuous Galerkin finite element method axially to overcome the drawbacks of MOCFE with the full 3D MOC representation which requires enormous memory and computational efforts. The Reactor Critical Facility (RCF), an open tank research reactor at the Rensselaer Polytechnic Institute (RPI), became operational in 1956 under the American Locomotive company and has been operating using low-enriched uranium fuel pins, which makes it the only university research facility to use fuel rods similar to operating commercial light water reactors (LWRs). This paper presents modeling and simulation results of the RCF using PROTEUS for code verification and validation which also supports an ongoing NEUP project in which a set of benchmark validation experiments for multiphysics coupling is being developed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Performance of an Astrophysical Radiation Hydrodynamics Code under Scalable Vector Extension Optimization

We present results of a performance study of an astrophysical radiation hydrodynamics code, V2D, on the Arm-based A64FX processor developed by Fujitsu. The code solves sparse linear systems, a task for which the A64FX architecture should be well suited. Here, we performed the performance analysis study on Ookami, an Apollo 80 platform utilizing the A64FX processor. We explored several compilers and performance anal-ysis packages and found the code did not perform as expected under scalable vector extension optimization, suggesting that a “deeper dive” into analyzing the code is worthwhile. However, a simple driver program that exercised basic sparse linear algebra routines used by V2D did show significant speedup with the use of the scalable vector extension optimization. We present the initial results from the study which used V2D on a relatively simple test problem that emphasized the repeated solution of sparse linear systems.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-T-254: Optimization of GATE and PHITS Monte Carlo Code Parameters for Uniform Scanning Proton Beam Based On Simulation with FLUKA General-Purpose Code

Purpose: Monte Carlo codes are becoming important tools for proton beam dosimetry. However, the relationships between the customizing parameters and percentage depth dose (PDD) of GATE and PHITS codes have not been reported which are studied for PDD and proton range compared to the FLUKA code and the experimental data. Methods: The beam delivery system of the Indiana University Health Proton Therapy Center was modeled for the uniform scanning beam in FLUKA and transferred identically into GATE and PHITS. This computational model was built from the blue print and validated with the commissioning data. Three parameters evaluated are the maximum step size, cut off energy and physical and transport model. The dependence of the PDDs on the customizing parameters was compared with the published results of previous studies. Results: The optimal parameters for the simulation of the whole beam delivery system were defined by referring to the calculation results obtained with each parameter. Although the PDDs from FLUKA and the experimental data show a good agreement, those of GATE and PHITS obtained with our optimal parameters show a minor discrepancy. The measured proton range R90 was 269.37 mm, compared to the calculated range of 269.63 mm, 268.96 mm, and 270.85 mm with FLUKA, GATE and PHITS, respectively. Conclusion: We evaluated the dependence of the results for PDDs obtained with GATE and PHITS Monte Carlo generalpurpose codes on the customizing parameters by using the whole computational model of the treatment nozzle. The optimal parameters for the simulation were then defined by referring to the calculation results. The physical model, particle transport mechanics and the different geometrybased descriptions need accurate customization in three simulation codes to agree with experimental data for artifact-free Monte Carlo simulation. This study was supported by Grants-in Aid for Cancer Research (H22-3rd Term Cancer Control-General-043) from the Ministry of Health, Labor and Welfare of Japan, Grants-in-Aid for Scientific Research (No. 23791419), and JSPS Core-to-Core program (No. 23003). The authors have no conflict of interest.

60 APPLIED LIFE SCIENCES↗

A Deterministic and Probabilistic Framework for Forecasting of Time-Series Damage States and Associated End-of-lIfe of a Pressurizer Water Reactor Surge Line under Design-Basis and Grid-Load-Following Loading Conditions

Argonne National Laboratory in support of the Department of Energy’s Light Water Reactor Sustainability (LWRS) program. The present semi-annual report covers the period between May 2018 and September 2018. In this report, we present the work performed to further improve the capability for structural integrity prediction. Part of this work focused on validating a computational fluid dynamics model for modeling thermal stratification. We also present finite element (FE) model based thermal-mechanical stress analysis of a pressurized water reactor (PWR) surge line under grid-load-following loading condition. In addition, we present work related to end-of-life prediction of a PWR surge line based on ASME code and NUREG-6909 deterministic approaches under design-basis and grid-load-following conditions. Our research in this report period also majorly involved with the development of a strategy for probabilistic fatigue life estimation. First, we applied a Weibull probabilistic modeling approach based on end-of-life fatigue data (traditional S~N data under in-air and PWR water conditions) to estimate the probabilistic life for a given strain amplitude. Second, we developed a data analytics based Markov-Chain-Monte-Carlo (MCMC) model for probabilistic modeling of time-series fatigue damage states and the associated probabilistic end-of-life of a component. Unlike the Weibull type approach, the MCMC type may not require a large number of fatigue tests for estimating the end-of-life for given loading and environmental conditions. We demonstrated the use of this framework through symmetric fatigue tests (stress- or strain-ontrolled tests either under in-air or PWR-water conditions) obtained from our earlier experiments. The MCMC model can use base time-series data either obtained from a model or fatigue experiment. Since experimental data are typically more accurate than model-calculated data, we used experiment-based time-series data to demonstrate the MCMC framework. However, when performing a complex experiment is not possible (due to time and costs involved), a model-based approach can be used to generate the time-series data. Similar model-based data are presented in our earlier report (e.g. for cyclic hardening damage states of a PWR surge line under constant amplitude and design-basis loading [1]). However, for the MCMC work in this report, we used experiment-based time-series damage states obtained from fatigue experiments under design-basis and grid-load-following loading conditions. The experimental loading inputs were based on thermal-mechanical FE model results of a typical PWR surge line. Then, we used the experimental time-series damage state data for forecasting the probabilistic time-series damage states and end-of-life of the PWR surge line. Based on this MCMC modeling strategy, and with the assumed loading, environment and thermal-mechanical boundary conditions, we have estimated that a 316SS-pure base metal PWR surge line would survive at least 159 years with zero failure probability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Enhanced Visualization Capabilities in OpenMC: Report of Activities under a DOE GAIN Voucher

This work was performed under a Cooperative Research and Development Agreement (CRADA) through the U.S. Department of Energy’s Gateway for Accelerated Innovations in Nuclear (GAIN) program in partnership with Oklo Inc., the industry participant. The broad goal of this project was to improve the visualization capabilities of OpenMC related to model development and simulation results. OpenMC is a Monte Carlo particle transport code that is developed and maintained by a broad community of contributors, of which Argonne National Laboratory plays a leading role.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Validation of Numerical Tools for Calculating Reactivity Feedback in Sodium Fast Reactors Using SEFOR Experimental Data

The Southwest Experimental Fast Oxide Reactor (SEFOR) was an experimental sodium-cooled fast breeder reactor operated from 1969 to 1972 with experiments designed to measure Doppler reactivity feedback in a wide temperature range from around 350 °F to temperatures approaching the melting point of mixed oxide fuel of around 5000 °F, providing valuable data for code validations. Co-supported by the Department of Energy (DOE) Fast Reactor Program (FRP) and the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the SEFOR benchmark project focused on using the experimental data to validate numerical tools that are used in industry and academia to design and license sodium-cooled fast reactors (SFRs). By the end of FY-25, substantial progress was achieved in the SEFOR benchmark study. A variety of numerical tools commonly used for modeling SFRs were applied to develop models for SEFOR core configurations I-D, I-E, I-I, and I-J. These included Monte Carlo codes such as MCNP, Serpent, and Shift; deterministic codes such as the legacy Argonne Reactor Computation (ARC) suite and the high-fidelity NEAMS code Griffin; and the system analysis code SAS4A/SASSYS-1 (SAS). Using these models, both SEFOR zero-power experiments and power-ascending tests were successfully simulated. Comparisons were performed against experimental measurements of core criticalities, reflector worth, kinetics parameters (Λ/βeff), isothermal reactivity feedback (from 350 °F to 760 °F at zero power), and power-ascending reactivity feedback (as power increased from 0.4 MW to 17 MW). In general, these comparisons demonstrated very good agreement between numerical results and experimental data. In Fiscal Year 26 (FY-26), the SEFOR benchmark project will continue to address the modeling issues identified in FY-25. Effort will focus on the simulation of reactivity insertion transients in SEFOR core II using the ARC/SAS model. Future work will also focus on incorporating BISON into the SEFOR core modeling process to enable the first Multiphysics simulations of the isothermal tests based on the MOOSE framework.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

SU-F-R-48: Early Prediction of Pathological Response of Locally Advanced Rectal Cancer Using Perfusion CT:A Prospective Clinical Study

Purpose: To prospectively evaluate the tumor vascularity assessed by perfusion CT for prediction of chemo-radiation treatment (CRT) response in locally advanced rectal cancer (LARC). Methods: Eighteen consecutive patients (61.9±8.8 years, from March–June 2015) diagnosed with LARC who underwent 6–8 weeks CRT followed by surgery were included. The pre-treatment perfusion CT was acquired after a 5s delay of contrast agent injection for 45s with 1s interval. A total of 7-cm craniocaudal range covered the tumor region with 3-mm slice thickness. The effective radiation dose is around 15mSv, which is about 1.5 the conventional abdomen/pelvis CT dose. The parametric map of blood flow (BF), blood volume (BV), mean transit time (MTT), permeability (PMB), and maximum intensity map (MIP) were obtained from commercial software (Syngo-CT 2011A, Siemens). An experienced radiation oncologist outlined the tumor based on the pre-operative MR and pathologic residual region, but was blinded with regards to pathological tumor stage. The perfusion parameters were compared to histopathological response quantified by tumor regression grade as good-responder (GR, TRG 0-1) vs. non-good responder (non-GR). Furthermore, the predictive value for pathological complete response (pCR) was also investigated. Results: Both BV (p=0.02) and MTT (P=0.02) was significantly higher and permeambility was lower (p=0.004) in the good responders. The BF was higher in GR group but not statistically significant. Regarding the discrimination of pCR vs non-pCR, the BF was higher in the pCR group (p=0.08) but none of those parameters showed statistically significant differences. Conclusion: BV and MTT can discriminate patients with a favorable response from those that fail to respond well, potentially selecting high-risk patients with resistant tumors that may benefit from an aggressive preoperative treatment approach. However, future studies with more patient data are needed to verify the prognostic value of perfusion CT especially for pCR prediction. This work is supported by the National High-tech R&D program for Young Scientists by the Ministry of Science and Technology of China (Grant No. 2015AA020917), Natural Science Foundation of China (NSFC Grant No. 81201091).

60 APPLIED LIFE SCIENCES↗

FLASSH 1.0: Full Law Analysis Scattering System Hub

Nuclear technologies, such as thermal reactors and critical systems, are driven by reaction rates which are directly dependent on thermal neutrons. The ability to calculate reaction rates is directly linked to the capability to predict and understand the cross sections. For atoms bound within a material structure, the available modes of energy exchange are captured in the Thermal Scattering Law (TSL, i.e. S(α,β)). The TSL is fundamentally a material property and can be defined explicitly to account for the excitation and deexcitation of the lattice energy states. This work describes the development of the Full Law Analysis Scattering System Hub (FLASSH) which calculates the TSL and cross sections. FLASSH is a modern code which has been designed to produce accurate TSL evaluations, remove historical approximations, and improve formulations while making advance physics accessible and simple to use.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling Approach to Critical in the Upcoming CERBERUS Experiment

The planned Critical Experiment Reflected By coppEr to betteR Understand Scattering (CERBERUS) seeks to maximize sensitivity to elastic neutron scattering in the intermediate energy region (0.625 eV to 100 keV). It will be performed near the end of FY 2023 at the National Criticality Experiments Research Center (NCERC). Very few International Criticality Safety Benchmark Evaluation Project (ICSBEP) benchmarks are sensitive to copper in this region. Creating a benchmark sensitive in this area will make the Zeus benchmark series, an intermediate benchmark evaluation that uses a copper reflector, more useful for code and nuclear data validation. Approaching criticality in a safe manner is of utmost importance to avoid a criticality accident, which would present a safety concern and could cause damage to equipment. The two rules that are followed closely to ensure that a criticality accident does not occur are the 3/4 rule and the 1/2 rule. The 3/4 rule states that no more than 3/4 of a critical mass can be assembled by hand, and the 1/2 rule states that no more than 1/2 of the material expected to reach criticality, or 1/2 the distance needed to reach criticality, can be added before another measurement of the count rate has been taken. This work will discuss the approach to criticality as modeled in MCNP6® 1 particle transport code with the ENDF/B- VIII.0 cross-section library and the .00c data library.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

IRIS Reimagined: Advancements in Intelligent Runtime System for Task-Based Programming

Task-based programming models are gaining traction in scientific computing. IRIS is a portable runtime system that exploits multiple heterogeneous programming systems and can discover available resources and manage multiple diverse programming systems (e.g., CUDA, Hexagon, HIP, Level Zero, OpenCL, and OpenMP) simultaneously. It accounts for the constraints of task dependencies and provides customizable scheduling policies to map those tasks to heterogeneous devices. In this paper, we present new capabilities added to IRIS to improve its portability for heterogeneous programming, build-friendliness, and performance efficiency. The new additions include vendor-specific kernel support, a runtime system with a foreign function interface to eliminate writing wrapper or boilerplate code for heterogeneous kernels, an easy-to-use and configurable CMake-based build environment, automatic and efficient data transfers and orchestration, and the Hunter and DAGGER toolchains to evaluate IRIS’s task scheduling algorithms.

Miniskar, Narasinga Rao↗

tapir: A tool for topologies, amplitudes, partial fraction decomposition and input for reductions

The demand for precision predictions in the field of high energy physics has dramatically increased over recent years. Experiments conducted at the LHC, as well as precision measurements at the intensity frontier such as Belle II require equally precise theoretical predictions to make full use of the acquired data. To match the experimental precision, second-, third- and, for certain quantities, even higher-order calculations in perturbative quantum field theory are required. To facilitate such calculations, computer software automating as many steps as possible is required. Yet, each calculation poses different challenges and thus, a high level of configurability is required. In this context we present tapir: a tool for identification, manipulation and minimization of Feynman integral families. It is designed to integrate in toolchains based on the computer algebra system FORM, the use of which is common practice in the field. tapir can be used to reduce the complexity of multi-loop problems with cut-filters, topology mapping, partial fraction decomposition and alike. Program Title:tapir CPC Library link to program files:https://doi.org/10.17632/ptc9t46xyn.1 Developer's repository link:https://gitlab.com/tapir-devs/tapir Licensing provisions: GPLv3 Programming language:python 3, C++ Nature of problem: Multi-loop computations require the automatization of a large number of different tasks related to Feynman integral topologies. Among them are the identification and minimization of integral topologies, partial fraction decomposition of topologies in the case of linearly dependent propagators as well as mapping scalar products of loop momenta to scalar functions. Solution method: The minimization of topologies is performed by comparison of their respective Nickel indices [1], even further minimization utilizes Pak's algorithm [2]. To efficiently map scalar products of loop momenta to scalar functions FORM [3] code is generated. Additional comments including restrictions and unusual features: Minimization based on Pak's algorithm slows down for many lines and scales. A coarser minimization using the Nickel indices, however, is still possible. [1]B. Nickel, D. Meiron, G.A.J. Baker, Compilation of 2-pt and 4-pt graphs for continuous spin model, Report, University of Guelph, 1977.[2]A. Pak, J. Phys. Conf. Ser. 368 (2012) 012049, https://doi.org/10.1088/1742-6596/368/1/012049, arXiv:1111.0868.[3]B. Ruijl, T. Ueda, J. Vermaseren, FORM version 4.2, arXiv:1707.06453, 7 2017.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Wootz: a compiler-based framework for fast CNN pruning via composability

Convolutional Neural Networks (CNN) are widely used for Deep Learning tasks. CNN pruning is an important method to adapt a large CNN model trained on general datasets to fit a more specialized task or a smaller device. The key challenge is on deciding which filters to remove in order to maximize the quality of the pruned networks while satisfying the constraints. It is time-consuming due to the enormous configuration space and the slowness of CNN training.The problem has drawn many efforts from the machine learning field, which try to reduce the set of network configurations to explore. This work tackles the problem distinctively from a programming systems perspective, trying to speed up the evaluations of the remaining configurations through computation reuse via a compiler-based framework. We empirically uncover the existence of composability in the training of a collection of pruned CNN models, and point out the opportunities for computation reuse. We then propose composability-based CNN pruning, and design a compression-based algorithm to efficiently identify the set of CNN layers to pre-train for maximizing their reuse benefits in CNN pruning. We further develop a compiler-based framework named Wootz, which, for an arbitrary CNN, automatically generates code that builds a Teacher-Student scheme to materialize composability-based pruning. Experiments show that network pruning enabled by Wootz shortens the state-of-art pruning process by up to 186X while producing significantly better pruning results.

Guan, Hui↗

Seismic Soil-Structure Interaction Analyses of a Deeply Embedded Model Reactor – SASSI Analyses

This report summarizes the SASSI analyses of a deeply embedded reactor model performed by BNL and CJC and Associates, as part of the seismic soil-structure interaction (SSI) simulation capability project for the NEAMS (Nuclear Energy Advanced Modeling and Simulation) Program of the Department of Energy. The SASSI analyses included three cases: 0.2 g, 0.5 g, and 0.9g, all of which refer to nominal peak accelerations at the top of the bedrock. The analyses utilized the modified subtraction method (MSM) for performing the seismic SSI evaluations. Each case consisted of two analyses: input motion in one horizontal direction (X) and input motion in the vertical direction (Z), both of which utilized the same in-column input motion. Besides providing SASSI results for use in comparison with the time domain SSI results obtained using the DIABLO computer code, this study also leads to the recognition that the frequency-domain method should be modernized so that it can better serve its mission-critical role for analysis and design of nuclear power plants.

42 ENGINEERING↗

Performance Improvement of the Linear System Solver for CMFD Acceleration in PROTEUS-MOC

The PROTEUS code developed under the DOE NEAMS program is a high-fidelity capable whole core transports solver based on unstructured finite element meshes [1]. PROTEUS includes three distinct transport solvers: SN, MOC and NODAL. Among them, the PROTEUS-MOC is the 3D method of characteristics (MOC) solver for extruded geometries. In order to accelerate the convergence of 3D MOC calculation, the CMFD acceleration method was implemented in PROTEUS-MOC [2]. For solving the CMFD linear system, the unpreconditioned GMRES algorithm was used in the previous work [3, 4] for implementing the CMFD scheme because of low computation cost of coarse mesh calculations. In the performance tests, we found that the overhead of the CMFD acceleration became non-negligible as the condition number of the CMFD linear system and the size of problem increase. This indicated that the CMFD overhead should be reduced in the preconditioning technique in order to improve the computational performance. In this work, the performance of CMFD acceleration was notably improved by developing and implementing the new formulation of CMFD linear system and preconditioning technique.

Jung, Yeon Sang↗

Advancement of Functional Expansion Tallies Capabilities in Serpent 2

The need for accident tolerant fuels (ATFs) became self-evident after the series of unfortunate nuclear events at Fukushima Daiichi in 2011. Restarting the Transient Reactor Test Facility (TREAT) at INL is one component in this concerted effort to establish the ATF development and testing program under the United States Department of Energy. Prior to the closure of TREAT in 1994, each experiment required multiple low-power calibration tests before and after each run. Since then, computing advancements have set computer-based modeling and simulation (M and S) as the preferred calibration method. The anticipated savings are both temporal and monetary. MAMMOTH, a collection of M and S tools built from the finite element-based MOOSE Framework, has been established to support the designed experiments intended for TREAT. This collection of tools includes: 1) Rattlesnake, a striking diffusion-based neutronics code, 2) BISON, an advanced materials and fuels performance code, and 3) RELAP-7, a powerful thermal-hydraulics solver. Each is capable of independently performing its function accurately and effectively; unified, they create a powerful multiphysics reactor core M and S application. One unfortunate weakness in this configuration lies in the unique geometry and design of TREAT's core. Although simple to model, this poses a challenge to the diffusion solver in Rattlesnake. One alternative is to replace Rattlesnake with Serpent, a fully-fledged Monte Carlo (MC) neutronics code. Coincidentally, Rattlesnake calculations are frequently validated against Serpent simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗