Results for “script” — Page 64

DOE OSTI2016

FM-BEM simulation of 3-D microvoid and microcrack graphite models with Python and ABAQUS/CAE pre-processing - 18603

The mechanical property and structural integrity of graphite structure have an important significance on the safe operation of high temperature gas cooled reactor. The graphite is porous media, the structural deficiencies such as microcracks and microvoids in it have great influence on its material mesoscopic structure and mechanical parameter. The shape and size distribution of microvoids and microcracks at different oxidation degree have a great influence on the overall performance. In order to simulate the geometric structure of the oxidized graphite, we used Python script to pre-process in Abaqus/CAE and this leads to different graphite models with different oxidation degrees. Then, these models were set with our pre-designed force boundary conditions and displacement boundary conditions. The models information was output to a FMBEM solver written by C++ code and the graphite s equivalent young s modulus and other mechanic parameters were calculated. In this paper, based on the 3- D microvoid and microcrack graphite models simulation, we studied the influence of the change of mesoscopic structure on graphite's mechanical property. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2016

Solution of the C5G7 Benchmark using the MOCUM Transport Code with ANSYS Unstructured Mesh

The objective of this research is to use finite element analysis (FEA) software to perform the unstructured meshing for reactor geometries and port the mesh information to MOCUM code. This workflow takes the advantage of the mature meshing module of the mainstream FEA; reduces the development effort for nuclear codes; and increases the nuclear code accuracy by using high quality meshes. In this work, ANSYS is used to build the reactor model and perform the unstructured meshing. The mesh information is then imported to the MOCUM code, a method of characteristic (MOC) transport program supporting unstructured meshes, using a developed MATLAB script. This procedure is verified by modeling the C5G7 benchmark problem. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2016

Dakota Uncertainty Quantification Methods Applied to the CFD code Nek5000

Predictive modeling and simulation of the performance of a nuclear reactor and its fuel is a challenging task because of the large number of coupled physical phenomena that must be addressed. In addition to this intrinsic complexity, model uncertainty must be accounted for in any analysis if the model will be used to facilitate design or operational decisions which may impact safety and performance. Rigorous, structured uncertainty analyses are performed by first characterizing the model's input uncertainties and, by propagating the uncertainty through the model in order to estimate the output's uncertainty. This paper presents the state of advancement of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) project to characterize the uncertainty of the CFD code Nek5000 using the Dakota package. Nek5000 is a high order spectral element computational fluid dynamics (CFD) code developed at Argonne National Laboratory for high resolution spectral-filtered large eddy simulations (LESs) and unsteady Reynolds averaged Navier-Stokes (RANS) simulations. The Dakota package developed at Sandia National Laboratory can be integrated with many scientific and engineering codes to facilitate efficient, effective uncertainty quantification (UQ) and sensitivity analyses. Dakota is designed to interface with external packages/codes such as Nek5000 through a script-driven interface called fork. Such interface provides maximum flexibility to the users for pre- and post-processing tasks preceding and following each function evaluation. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2015

NJOY21: Making NJOY Suitable for the 21. Century - 14409

NJOY has a long history of nuclear data processing. NJOY began as a successor to MINX (Multigroup Interpretation of Nuclear X-Sections), a Los Alamos code and was first publicly released in the summer of 1977. The initial version included modules such as RECONR for resonance reconstruction on a union energy grid, BROADR for Doppler broadening plus a new Multigroup averaging program, GROUPR and an associated output module, DTR. Soon thereafter the HEATR, ERRORR, and THERMR modules were added to NJOY. By the mid 1980's, NJOY had a worldwide user community. The current version of NJOY2012 was released in December 2012. NJOY2012 builds on the success of more than 45 years of nuclear data processing at Los Alamos National Laboratory. NJOY2012 was rewritten using Fortran 90/95, utilizing among other things the concepts of public and private variables, eliminating common blocks in lieu of Fortran modules and taking advantage of allocatable arrays. A new version of NJOY is being developed at Los Alamos National Laboratory; NJOY21: NJOY for the 21. Century. NJOY21 uses modern programming standards of quality, reliability, and maintainability. It is being designed from the ground up to utilize modern methods of interacting with computer codes by providing an interactive and scripting (e.g., Python) interface. This will allow scientists the ability to more easily work with the results of the various steps in the processing of nuclear data. NJOY21 will continue the legacy of NJOY by providing a trusted and verified processing code that is fast, flexible, easy to use and customize and capable of handling both modern and established nuclear data formats. This paper introduces NJOY21 by first relating the history of NJOY; the majority of the paper will discuss NJOY21. NJOY21 will take advantage of the many decades of nuclear data processing experience that is contained in NJOY2012 and previous versions of NJOY and extend the capabilities and convenience of processing nuclear data. We will discuss how NJOY21 will remain backwards compatible with NJOY2012 - i.e., how current NJOY2012 input decks will continue to work with NJOY21. We will also discuss the verification and validation of the NJOY21 code base. Finally we will discuss new features in NJOY21 that do not exist in NJOY2012. These features include changes in how a scientist can interact with NJOY to process and verify and validate nuclear data; and the ability to read and write new nuclear data formats. (authors)

42 ENGINEERING↗

DOE OSTI2022

Native Natural circulation in VERA

This work seeks to introduce a native natural circulation solver to VERA's thermohydraulics capabilities. This work begins with initial investigations into the viability of such a solver by determining the fundamental equations as well as necessary inputs. The method was then implemented in the form of a Python script coupled to CTF which was compared to TRACE. Once viability of the natural circulation solver was established, it was then implemented natively in MPACT for use with previously created simplified thermohydraulics solvers. This new implementation of the native natural circulation solver in MPACT is then tested for single physics problems for both a BWR and PWR single channel model and compared to TRACE where it shows an error of about 30% for BWRs and about 15% for PWRs. The solver is additionally tested for NuScale based models where it demonstrates an error to the reported flow rate of about 10%, resulting in an eigenvalue error of only 138 pcm when comparing fixed flow rate multiphysics runs. Finally, the solver is tested in multiphysics systems for a small PWR assembly, a small BWR assembly, and a full core approximate model of the NuScale design. In these multiphysics runs, it demonstrates no impact on multiphysics convergence stability, as well as negligible impact on runtimes. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Implementation of a Time-domain Cosmic-ray-muon Tagger for the NEXUS Low-background Cryogenic Facility

Quantifying the effects of radiation on the operation of qubits both as quantum information systems as well as particle detectors has emerged as a pressing issue in quantum science in recent years. We present an overview of the design, operation, and deployment of a 90-$\mathrm{in}^2$ three-panel muon detector for use in the NEXUS experimental facility at Fermilab to temporally isolate correlated errors in qubits and determine if they possess an astrophysical origin. Constructed with three scintillator-attached PMTs read out with NIM modules in a triply-coincident logic scheme, we measure a surface-level muon luminosity of 9.7425 muons per second---consistent with an average surface-level cosmic-ray flux of approximately one muon per square centimeter per minute. Integrating the NIM modules with an MCC 128 DAQ HAT and Raspberry Pi, a Python script records exactly when a muon struck the detector and writes a timestamp to a log file for follow-up cross referencing. This experiment will broadly contribute to further studies aimed at understanding the source and mitigation of information loss in qubits.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

MLOps for Beam Controls

Machine learning operations (MLOps) is the standardization and streamlining of the ML development lifecycle to address the challenges associated with large-scale machine learning applications. The full MLOps pipeline consists of open-source tools: DataHub, MinIO and MLflow. It is being used for dataset management and model development to handle changing data dependencies, varying business needs, reproducibility, and diverse teams working with differing tools and skills. To demonstrate the completion of an MLOps pipeline for particle accelerator operations, we are deploying a simple script that computes settings for the Booster’s gradient magnet power supply. Once the demonstration is complete, we will develop and deploy ML-based optimization algorithms to improve Booster’s overall efficiency. This MLOps pipeline opens the gate to systematically develop and deploy ML applications for accelerator controls and diagnostics.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Alfalfa Virtual Building Service: Software Engineering Best Practices Applied to Runtime Interaction with Building Energy Models

Buildings are active participants in increasingly complex energy systems. Building Energy Modeling (BEM) has a key role to play in planning and de-risking an equitable energy transition, with BEM-backed "virtual buildings" critical path for diverse applications that include workforce training tools, Hardware-in-the-Loop (HIL) experimentation to study equipment performance under a range of conditions, Control-Hardware-in-the-Loop (CHIL) experimentation to de-risk commercial control implementations at equipment through grid orchestration levels, and integration of dynamic load profiles into grid modeling tools for energy system experimentation at the urban scale. Modeling requirements vary across these applications, but many software engineering tasks do not. The Alfalfa Virtual Building Service (AVBS, see https://github.com/NREL/alfalfa/wiki) is an open-source web service that solves these common tasks robustly in one place, providing a foundational platform for power users to bootstrap their own applications. AVBS abstracts the specifics of runtime interaction with OpenStudio, Modelica, and Spawn of EnergyPlus models behind a unified REST API. Additionally, AVBS provides resources for cloud deployment and scaling to 100s of parallel simulations, a growing library of modular Operational Technology (OT) integrations for emulation of real-world interfaces, and scripts to automate the population of communities of virtual buildings from URBANopt, ResStock and ComStock.

building automation↗

DOE OSTI2024

MLOps for Beam Controls

Machine learning operations (MLOps) is the standardization and streamlining of the ML development lifecycle to address the challenges associated with large-scale machine learning applications. The full MLOps pipeline consists of open-source tools: DataHub, MinIO and MLflow. It is being used for dataset management and model development to handle changing data dependencies, varying business needs, reproducibility, and diverse teams working with differing tools and skills. To demonstrate the completion of an MLOps pipeline for particle accelerator operations, we are deploying a simple script that computes settings for the Booster’s gradient magnet power supply. Once the demonstration is complete, we will develop and deploy ML-based optimization algorithms to improve Booster’s overall efficiency. This MLOps pipeline opens the gate to systematically develop and deploy ML applications for accelerator controls and diagnostics.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Characterizing Supernovae Neutrinos in MARLEY

This summer I worked on characterizing the uncertainties in neutrino-argon interactions due to the optical model potential. Functionally, this was a great deal of new experience coding in C++. I wrote a script using MARLEY(Gardiner 2019) to compute both the branching ratios and exclusive cross sections for a whole host of different decay processes from the unstable Potassium nucleus created by a v+Ar charged current interaction. These decays were computed through MARLEY, governed by a global parametrization of the Optical Model potential(Koning-Delaroche 2003). Helpfully, Pruitt+ 2023 provided public supplemental material outlining 400 separate "universes" for the optical model potential parameters. These "universes" represent uncertainties on the optical potential itself, and I wrote code to incorporate them into the MARLEY framework and quantify their impact on the predictions of neutrino-argon cross sections. This stands as the first computation of nuclear de-excitation systematic errors on cross sections in this energy band for v+Ar interactions, a critical step to understanding the science we can do with DUNE in the tens-of-MeV band.

Abu El-Haj, Luca↗

DOE OSTI2024

Axion Dark Matter EXperiment 2A Cavity Characterization

We cooled the ADMX 2A cavities down to 5K and found the Q-factors for various tunings, which were about a factor of five larger than the Q-factors from the previous cooldown. We also found that the tuning mechanism with the piezo motor and gearbox is not currently equipped to handle cryogenic temperatures. We created an experimental warm bead pull setup using piezo motors, an alumina bead, a pulley system, and a distance calibrating script. We found that the bead pull method can successfully show us the reflection-induced electric field mode structure. This technique allowed us to differentiate between different TE and TM modes near mode crossings. However, some bead pull sweeps produced noisy or asymmetrical results, possibly due to asymmetries in the tuning rod orientation or antenna. Future iterations of the bead pull method can help us estimate the cavity form factor and maximize our scan range near mode crossings.

Maglione, D.↗

DOE OSTI2024

Simulations, Modeling and Data Analysis of Parity Violating Electron Scattering Experiments

In the Standard Model (SM) of nuclear and particle physics, parity violation is incorporated through the representation of the weak interaction as a chiral gauge interaction. Only the left-handed components of particles and right-handed components of antiparticles participate in weak interactions in the Standard Model. This implies that parity is asymmetric for the weak interaction. Parity violating electron scattering (PVES) experiments are designed to probe the physics parameters related to the SM, with the possibility to discover physics beyond the SM (BSM) by measuring the parity violating asymmetry ¿¿¿ of longitudinally polarized electrons scattered off unpolarized targets with high precision. This dissertation will be focused on two PVES experiments, the next 208Pb Lead Radius Experiment (PREX-II), and the Measurement of a Lepton-Lepton Electroweak Reaction (MOLLER) experiment, as well as in some small sections, the Calcium Radius Experiment (CREX) and P2 experiment which are also PVES experiments). PREX-II and CREX experiments, performed in Hall A at the Thomas Jefferson National Accelerator Facility (Jefferson Lab), measured ¿¿¿ in the elastic scattering of longitudinally polarized electrons from 208Pb and 48Ca targets to provide a precise model independent determination of the neutron skin thickness of 208Pb and 48Ca nuclei, respectively. The MOLLER experiment, proposed to start in 2027 and also to be performed in Hall A at Jefferson Lab, is to measure ¿¿¿ of longitudinally polarized electrons scattered off unpolarized electrons (Møller scattering) to determine the weak charge of electrons ¿¿¿ and the weak mixing angle ¿¿ with high precision. As for the P2 experiment, which will be performed at the upcoming MESA accelerator in Mainz Germany, it is to measure the weak charge of proton ¿¿¿ using ¿¿¿ in the elastic electron-proton scattering of polarized electrons off unpolarized protons. The final results from the PREX-II experiment are presented as ¿¿¿=550±16 (¿¿¿¿)±8 (¿¿¿¿) parts-per-billion (ppb). Combining the PREX-I and PREX-II results, the neutron skin thickness from PREX experiments is determined as ¿¿-¿¿=0.283±0.071 ¿¿ in 208Pb. This thesis lists the software and computational contribution of the author to these PVES experiments, including writing scripts and software to help with the PREX-II/CREX experiments, analyzing data to provide useful information and systematic uncertainty for the PREX-II experiment, modeling and simulations for the MOLLER, and providing an alternative design of an electronic equipment for the P2 experiment.

Chen, Yufan↗

DOE OSTI2024

Simulations, Modeling and Data Analysis of Parity Violating Electron Scattering Experiments

In the Standard Model (SM) of nuclear and particle physics, parity violation is incorporated through the representation of the weak interaction as a chiral gauge interaction. Only the left-handed components of particles and right-handed components of antiparticles participate in weak interactions in the Standard Model. This implies that parity is asymmetric for the weak interaction. Parity violating electron scattering (PVES) experiments are designed to probe the physics parameters related to the SM, with the possibility to discover physics beyond the SM (BSM) by measuring the parity violating asymmetry ¿¿¿ of longitudinally polarized electrons scattered off unpolarized targets with high precision. This dissertation will be focused on two PVES experiments, the next 208Pb Lead Radius Experiment (PREX-II), and the Measurement of a Lepton-Lepton Electroweak Reaction (MOLLER) experiment, as well as in some small sections, the Calcium Radius Experiment (CREX) and P2 experiment which are also PVES experiments). PREX-II and CREX experiments, performed in Hall A at the Thomas Jefferson National Accelerator Facility (Jefferson Lab), measured ¿¿¿ in the elastic scattering of longitudinally polarized electrons from 208Pb and 48Ca targets to provide a precise model independent determination of the neutron skin thickness of 208Pb and 48Ca nuclei, respectively. The MOLLER experiment, proposed to start in 2027 and also to be performed in Hall A at Jefferson Lab, is to measure ¿¿¿ of longitudinally polarized electrons scattered off unpolarized electrons (Møller scattering) to determine the weak charge of electrons ¿¿¿ and the weak mixing angle ¿¿ with high precision. As for the P2 experiment, which will be performed at the upcoming MESA accelerator in Mainz Germany, it is to measure the weak charge of proton ¿¿¿ using ¿¿¿ in the elastic electron-proton scattering of polarized electrons off unpolarized protons. The final results from the PREX-II experiment are presented as ¿¿¿=550±16 (¿¿¿¿)±8 (¿¿¿¿) parts-per-billion (ppb). Combining the PREX-I and PREX-II results, the neutron skin thickness from PREX experiments is determined as ¿¿-¿¿=0.283±0.071 ¿¿ in 208Pb. This thesis lists the software and computational contribution of the author to these PVES experiments, including writing scripts and software to help with the PREX-II/CREX experiments, analyzing data to provide useful information and systematic uncertainty for the PREX-II experiment, modeling and simulations for the MOLLER, and providing an alternative design of an electronic equipment for the P2 experiment.

Chen, Yufan↗

DOE OSTI2025

Reassessing the Origins and Contemporary Relevance of ck Acceptability Parameters: Evolving Perspectives on Similarity

“Sensitivity and Uncertainty Analyses Applied to Criticality Safety Validation,” introduces sensitivity and uncertainty methods to address challenges in defining and extending areas of applicability for criticality safety validation. These areas are traditionally defined by the bounds or limits on key parameters, but establishing valid ranges and managing complex parameter variations remain challenging. NUREG/CR-6655 introduces ck and other integral indices, as well as concepts such as the completeness of benchmark coverage, to better quantify system similarities. The work proposed herein seeks to evaluate these foundational concepts to ensure that the bounds remain effective in guiding the assessment of similarity and applicability in modern applications. The concept of completeness, along with other parameters envisioned within the framework, serves as an example of the foundational ideas that have been established, though their effectiveness in practice may not be fully understood. Advancements in scripting tools, coupled with the speed and efficiency of modern computing and statistical models, now allow for faster and more thorough assessments than previously possible. These advancements also enable the identification of trends within the data, which could provide additional insight into system behavior and further broaden the scope of previously performed benchmarks. By leveraging these capabilities, we will revisit and expand the scope of these foundational methods to determine whether the necessary elements for robust similarity evaluation are already embedded, partially realized, or remain untapped.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Characterization of the X-ARAPUCA systems' response for SBND's light detection

X-ARAPUCA waveform data, collected at a test stand at the Noble Liquid Test Facility (NLTF) and the Short Baseline Near Detector (SBND), were analyzed via python scripts. For the ARARA system, various distributions and correlations were studied, including the area and amplitude distributions of the waveform peaks for identifying single photoelectron peaks, which was able to be accomplished. The behavior of the baseline artifact was also characterized, which appeared only at high voltages/bright events, and separates further from the prompt peak as the voltage/brightness increases. The ARARA system is being studied in order to be implemented into SBND, just as the APSAIA X-ARAPUCAs have been.

Novotny, Andrew [Fermilab]↗

DOE OSTI2025

Treyson Ricks - Intern Showcase Poster

Quinone-based sorbents offer a tunable, energy-efficient route to electrochemical CO2 capture, but systematic guidance for molecular design is lacking. Here, we report a high-throughput computational workflow that combines density functional theory (DFT) screening with machine-learning (ML) modeling to evaluate CO2 binding thermodynamics across several quinone derivatives, spanning benzoquinones, naphthoquinones, and anthraquinones. In addition to using solvents to stabilize the quinone anion and dianion, we studied the effect of ion-pairing on the reduction potentials and the CO2 binding energy. Automated Python scripts handled geometry optimizations and adduct-formation energies on an HPC cluster, reducing manual effort significantly. This integrated platform can uncover structure–property relationships and enables rapid in silico evaluation of untested candidates. We present one example from our workflow to showcase the capability of using quinones with ion-pairing to effectively capture CO2. Our approach paves the way for the rational selection of optimal quinone sorbents and can be extended with experimental thermochemical and kinetic data, alternative redox cycles, and stability assessments to accelerate development of next-generation electrochemical CO2 capture materials.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗

DOE OSTI2025

Immersive Scientific Visualization of Molten-Salt Reactor Waste Characteristics Using Virtual Reality

Immersive visualization is changing how we explore, communicate, and understand complex scientific systems. In nuclear energy, an area in which data are often multidimensional, time-dependent, and difficult to interpret, virtual reality (VR) represents a powerful and intuitive informational medium. This work introduces a VR-based platform that visualizes the post-shutdown behavior and waste management lifecycle of molten-salt reactors (MSRs), a next-generation reactor type with unique operational and safety characteristics. The platform, built in Unity, is streamed on the Meta Quest 3 headset. It transforms high-fidelity simulation data into an interactive, immersive experience. Users can explore time-dependent reactor characteristics such as nuclide decay, which is a key factor for evaluating reactor waste strategies. The datasets were generated using the MOOSE (Multiphysics Object-Oriented Simulation Environment) framework and then processed through ParaView scripting for smooth integration into Unity. From a visualization standpoint, the platform emphasizes spatial storytelling, temporal exploration, and user-centered interaction. Users can navigate 3D reactor geometries, slice through volumetric data, and manipulate time to observe how physical phenomena evolve. Real-scale rendering and embodied interaction make the experience feel tangible. The interface is designed to be accessible, even to those without nuclear or simulation expertise. This lowers the barrier for stakeholders, policymakers, and the general public, while still supporting expert analysis and collaborative decision-making. This work shows how immersive visualization can function as both a scientific tool and a communication interface. By integrating simulation, processing, and visualization into a cohesive workflow, we offer a scalable framework for immersive scientific storytelling. The modular design supports future extensions to other reactor types and lifecycle stages, from shutdown to long-term storage, making the platform adaptable for both research and outreach.

99 - GENERAL AND MISCELLANEOUS↗

DOE OSTI2026

Validation and Verification of Python based Neutron Spectrum Unfolding Software

To validate and verify the python-based code (PySL), designed to replicate the programs used by STAYSL for Beam Correction Factor (BCF) and Self-Shielding Factor (SHIELD), a series of tests were performed. To test BCF a python script was written to generate a random flux history file and both versions of the code processed the data. The test verified matching values up to at least one decimal place, approximately 10,000 tests where run and each one passed. Isotopes began to fail the tests once neutron saturation was reached. To verify this the total time of exposure was varied the isotopes that failed were compared to a list of their half-lives. The test process for SHIELD was very similar but, in this case, the code began by producing an input file with varying thickness and device type/environment for the SHIELD input. The failure condition for this test was if any of the data points for an isotope had a difference above 3%. Approximately 40 of these tests were run and there were only 3 isotopes that had reoccurring failures but only 2% of their points were above the 3% difference. A visual comparison was conducted by plotting the results from both programs. Although the test failed, the differences between their values were minuscule, and the self-shielding factor’s shape was preserved when plotted. Next steps for this project will be validating and verifying the python-based SigPhi code and then reproducing and testing the least squares unfolding performed by STAYSL.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗