Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed simulation”

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,117 records · Page 62

Demand Response Availability Profiles for California in the Year 2020

Demand response (DR) is being considered as a valuable resource for keeping the electrical grid stable and efficient, and deferring upgrades to generation, transmission, and distribution systems. However, simulations to determine how much infrastructure upgrades can be deferred are necessary in order to plan optimally. Production cost modeling is a technique, which simulates the dispatch of generators to meet demand and reserves in each hour of the year, at minimal cost. By integrating demand response resources into a production cost model (PCM), their value to the grid can be estimated and used to inform operations and infrastructure planning. DR availability profiles and constraints for 13 end-uses in California for the year 2020 were developed by Lawrence Berkeley National Laboratory (LBNL), and integrated into a production cost model by Lawrence Livermore National Laboratory (LLNL), for the California Energy Commission’s Value of Energy Storage and Demand Response for Renewable Integration in California Study. This report summarizes the process for developing the DR availability profiles for California, and their aggregate capabilities. While LBNL provided potential DR hourly profiles for regulation product in the ancillary services market and five-minute load following product in the energy market for LLNL’s study, additional results in contingency reserves and an assumed flexible product are also defined. These additional products are included in the analysis for managing high ramps associated with renewable generation and capacity products and they are also presented in this report.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Phase Transformation of Thermally Aged and Neutron Irradiated Duplex Stainless Steel Used in LWRs. Final report

The lifetime of reactor components made of duplex stainless steels can be limited by the embrittlement from thermal aging, neutron irradiation or a synergistic effect between thermal aging and neutron irradiation. There is still a large scientific knowledge gap in understanding the phase evolution in duplex stainless steels upon thermal aging with or without neutron irradiation, the synergistic effect between thermal aging and neutron irradiation, and the impacts of structural evolution on material mechanical responses. In this project, totally, five major tasks were completed through this project: a) the in-situ WAXS tensile tests of thermally aged cast duplex stainless steels, b) the EXAFS investigation of thermally aged cast austenitic stainless steel and welds of austenite sainless steel, c) the characterization of irradiation effect on neutron irradiated CASS using APT, d) characterize the ferrite hardening of neutron irradiated CASS using nanoindentation, and e) the FEM modelling to validate the in-situ WAXS tensile tests and to simulate the stress/strain distribution during deformation for thermally aged duplex stainless steel. In summary, the thermal aging without neutron irradiation can significantly increase the strength of ferrite phase through the spinodal decomposition and G-phase precipitations. The in-situ X-ray experiments show that the hardening effect in ferrite phase can significantly alter the load partition between the ferrite and austenite phases. Atom probe tomography characterizations on samples irradiated to 0.08, 3, 5, 10, 20, and 40 dpa shows that neutron irradiation has dominated the microstructural evolution after 3 dpa, and the spinodal decomposition has saturated yet even at 40 dpa based on the quantified wavelength and amplitude.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

WHOLESCALE - Water & Hole Observations Leverage Effective Stress Calculations And Lessen Expenses (Final Technical Report 2020 - 2024)

The WHOLESCALE acronym stands for Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses. The goal of the WHOLESCALE project is to simulate the spatial distribution and temporal evolution of stress in the geothermal system at San Emidio in Nevada, United States. To reach this goal, the WHOLESCALE team has developed a methodology to incorporate and interpret data from four methods of measurement into a multi-physics model that couples thermal, hydrological, and mechanical (T H-M) processes. The WHOLESCALE team has applied this methodology at the San Emidio geothermal field, located ~100 km north of Reno, Nevada in the northwestern Basin and Range province. The WHOLESCALE team includes 30 individuals working at two universities, two national laboratories, and one industry partner. Two master-degree students and five post-doctoral researchers have gained professional experience and earned partial financial support via the WHOLESCALE project. The WHOLESCALE team has taken advantage of the perturbations created by changes in pumping operations during planned shutdowns in 2016, 2021, and 2022 to infer temporal changes in the state of stress in the geothermal system at San Emidio, Nevada, U.S. The WHOLESCALE results support the working hypothesis that increasing pore-fluid pressure reduces the effective normal stress acting across fault zones. During normal operations, pumping in deep production wells decreases fluid pressures and thus increases the effective normal stresses on faults, reducing microseismicity. During planned shutdowns, the cessation of production increases pore-fluid pressure and reduces effective normal stress. The WHOLESCALE products generated during the 4-year period between 2020 and 2024 include: three articles published in the open-access, peer-reviewed scientific literature, two master’s theses, 20 presentations or papers at scientific conferences, and 17 data sets available on public repositories. The WHOLESCALE project has been completed in two phases that included three performance periods separated by two Go/No-go Stage Gate Reviews. Tasks were classified by data type (i.e., Geologic Structure, Borehole, Geodesy, Hydrology, Seismology, and Modeling). The first phase of the project started July 31, 2020 and included ongoing project coordination (Task 1), a project kickoff (Task 2), analysis of existing data (Task 3), development of the initial stress model & deployment design (Task 4), and Go/No-go Decision Point #1 (Task 5). Phase II began with implementing the 2022 deployment (Task 6), followed by Go/No-go Decision Point #2 (Task 7) The remainder of Phase II consisted of analyzing data collected during deployment (Task 8), calibration of the stress model on all observations (Task 9), and the Final Review (August 23, 2024) & Reporting (Task 10).

15 GEOTHERMAL ENERGY↗

Optimization-Based Resiliency Verification in Microgrids via Maximal Adversarial Set Characterization

Critical energy infrastructures are increasingly relying on advanced sensing and control technologies for efficient and optimal utilization of flexible energy resources. Algorithmic procedures are needed to ensure that such systems are designed to be resilient to a wide range of cyber-physical adversarial events. This paper provides a robust optimization framework to quantify the largest adversarial perturbation that a system can accommodate without violating pre-specified resiliency metrics. We formulate the maximal adversarial set characterization as a bi-level optimization problem which is solved via Lagrangian relaxations.We illustrate the proposed algorithm on an islanded microgrid example: a modified IEEE 123-node feeder with distributed energy resources. Simulations are carried out to characterize the tolerable adversarial perturbations for varying levels of available flexibility (energy reserves).

Nazir, Mohammad Nawaf↗

Application of Phase-field Method in Predicting Gas Bubble Microstructure Evolution in Nuclear Fuels

Fission product accumulation and gas bubble microstructure evolution in nuclear fuels strongly affect thermo-mechanical properties such as thermal conductivity, gas release, volumetric swelling and cracking, and hence the fuel performance. In this paper, a general phase-field model is developed to predict gas bubble formation and evolution. Important materials processes and thermodynamic properties including the generation of gas atoms and vacancies, sinks for vacancies and gas atoms, the elastic interaction among defects, gas re-solution, and inhomogeneity of elasticity and diffusivity are accounted for in the model. The simulations demonstrate the potential application of the phase-field method in investigating 1) heterogeneous nucleation of gas bubbles at defects; 2) effect of elastic interaction, inhomogeneity of material properties, and gas re-solution on gas bubble microstructures; and 3) effective properties from the output of phase-field simulations such as distribution of defects, gas bubbles, and stress fields.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of the Next-Generation Caustic-Side Solvent Extraction (NG-CSSX) Process for Cesium Removal from High-Level Tank Waste

This paper describes the chemical performance of the Next-Generation Caustic-Side Solvent Extraction (NG-CSSX) process in its current state of development for removal of cesium from the alkaline high-level tank wastes at the Savannah River Site (SRS) in the US Department of Energy (USDOE) complex. Overall, motivation for seeking a major enhancement in performance for the currently deployed CSSX process stems from needs for accelerating the cleanup schedule and reducing the cost of salt-waste disposition. The primary target of the NG-CSSX development campaign in the past year has been to formulate a solvent system and to design a corresponding flowsheet that boosts the performance of the SRS Modular CSSX Unit (MCU) from a current minimum decontamination factor of 12 to 40,000. The chemical approach entails use of a more soluble calixarene-crown ether, called MaxCalix, allowing the attainment of much higher cesium distribution ratios (DCs) on extraction. Concurrently decreasing the Cs-7SB modifier concentration is anticipated to promote better hydraulics. A new stripping chemistry has been devised using a vitrification-friendly aqueous boric acid strip solution and a guanidine suppressor in the solvent, resulting in sharply decreased DCs on stripping. Results are reported herein on solvent phase behavior and batch Cs distribution for waste simulants and real waste together with a preliminary flowsheet applicable for implementation in the MCU. The new solvent will enable MCU to process a much wider range of salt feeds and thereby extend its service lifetime beyond its design life of three years. Other potential benefits of NG-CSSX include increased throughput of the SRS Salt Waste Processing Facility (SWPF), currently under construction, and an alternative modular near-tank application at Hanford.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Predicting Reactor Antineutrino Emissions Using New Precision Beta Spectroscopy

Neutrino experiments at nuclear reactors are currently vital to the study of neutrino oscillations. The observed antineutrino rates at reactors are typically lower than model expectations. This observed deficit is called the “reactor neutrino anomaly”. A new understanding of neutrino physics may be required to explain this deficit, though model estimation uncertainties may also play a role in the apparent discrepancy. PNNL is currently investigating an experimental technique that promises reduced uncertainties for measured data to support these hypotheses and interpret reactor antineutrino measurements. The experimental approach is to 1) direct a proton accelerator beam on a metal target to produce a source of neutrons, 2) use spectral tailoring to modify the neutron spectrum to closely simulate the energy distribution of a power reactor neutron spectrum, 3) irradiate isotopic fission foils (235U, 238U, 239Pu, 241Pu) in this neutron spectrum so that fissions occur at energies representative of a reactor, 4) transport the beta particles released by the fission products in the foils to a beta spectrometer, 5) measure the beta energy spectrum, and 6) invert the measured beta energy spectrum to an antineutrino energy spectrum. A similar technique using a beta spectrometer and isotopic fission foils was pioneered in the 1980’s at the ILL thermal reactor. Those measurements have been the basis for interpreting all subsequent antineutrino measurements at reactors. A basic constraint in efforts to reduce uncertainties in predicting the antineutrino emission from reactor cores is any underlying limitation of the original measurements. This may include beta spectrum energy resolution, the absolute normalization of beta emission to number of fission, statistical counting uncertainties, lack of 238U data, the purely thermal nature of the IIL reactor neutrons used, etc. An accelerator-based neutron source that can be tailored to match various reactor neutron spectra provides an advantage for control in studying how changes in the neutron spectra (i.e. "in the reactor core") affects the resulting fission product beta spectrum. Furthermore, the 238U antineutrino spectrum, which has not been measured, can be studied directly because of the enhanced 1 MeV fast neutron flux available at the accelerator source. A facility such as the Project X Injector Experiment (PXIE) 30 MeV proton linear accelerator at Fermilab is being considered for this experiment. The hypothesis is that a new approach utilizing the flexibility of an accelerator neutron source with spectral tailoring coupled with a careful design of an isotopic fission target and beta spectrometer and the inversion of the beta spectrum to the neutrino spectrum will allow further reduction in the uncertainties associated with prediction of the reactor antineutrino spectrum.

Asner, David M.↗

SimulatorToFMU v0.1

SimulatorToFMU is a software package written in Python which allows users to export a memoryless Python-driven simulation program or script as a Functional Mock-up Unit (FMU) for model exchange or co-simulation.In CyDER (Cyber Physical Co-simulation Platform for Distributed Energy Resources in Smart Grids), SimulatorToFMU will allow exporting OPAL-RT as an FMU. This will enable OPAL-RT to be linked to CYMDIST and GridDyn FMUs through a standardized open source interface.

Nouidui, Thierry↗

The 3D clustering of Lyman Alpha Emitters measured with DESI

We present a clustering analysis of Lyman-$α$ emitters (LAEs) using spectroscopic observations from the Dark Energy Spectroscopic Instrument (DESI) of candidates selected from the Blanco/DECam Intermediate-Band Imaging Survey (IBIS). We measure the two-point correlation function and the power spectrum, including cross-correlations with DESI quasars. Using both analytical and halo occupation distribution (HOD) simulation-based modeling, we find a linear bias of $b \sim 2.31$--$2.62$ for LAEs over the redshift range $2.26 < z < 3.41$. The analytical modeling also provides constraints on the strength of radiative transfer effects, while the HOD analysis characterizes the LAE-halo connection across multiple models. Finally, we quantify the magnitude of non-perturbative clustering effects such as Fingers of God in the LAE population, providing essential input for the accurate modeling of LAE-based cosmological analyses in forthcoming high-redshift surveys such as DESI-II.

Ebina, H. [UC, Berkeley; LBL, Berkeley] (ORCID:000↗

Asymmetry of the Balmer-alpha line shape and recovery of the effective hydrogen temperature in the tokamak scrape-off layer

An algorithm for recovering the effective temperature of atoms of hydrogen (and its isotopes) in the tokamak scrape-off layer from the asymmetry of the Balmer-alpha line shape is proposed. The algorithm is based on the parametrization of the asymmetry of the line shape caused by the nonlocal character of neutral hydrogen flux from the wall into the tokamak plasma. The accuracy of the algorithm is tested against the results of simulations of the velocity distribution function of deuterium neutrals in the scrape-off layer by the EIRENE code with the use of the source data on the main plasma component in the quasi-stationary stage of the inductive mode of ITER operation calculated by the SOLPS4.3 (B2-EIRENE) code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of realistic PET simulations incorporating tumor patient's specificity using anthropomorphic models: Creation of an oncology database

Purpose: The GATE Monte Carlo simulation toolkit is used for the implementation of realistic PET simulations incorporating tumor heterogeneous activity distributions. The reconstructed patient images include noise from the acquisition process, imaging system's performance restrictions and have limited spatial resolution. For those reasons, the measured intensity cannot be simply introduced in GATE simulations, to reproduce clinical data. Investigation of the heterogeneity distribution within tumors applying partial volume correction (PVC) algorithms was assessed. The purpose of the present study was to create a simulated oncology database based on clinical data with realistic intratumor uptake heterogeneity properties.Methods: PET/CT data of seven oncology patients were used in order to create a realistic tumor database investigating the heterogeneity activity distribution of the simulated tumors. The anthropomorphic models (NURBS based cardiac torso and Zubal phantoms) were adapted to the CT data of each patient, and the activity distribution was extracted from the respective PET data. The patient-specific models were simulated with the Monte Carlo Geant4 application for tomography emission (GATE) in three different levels for each case: (a) using homogeneous activity within the tumor, (b) using heterogeneous activity distribution in every voxel within the tumor as it was extracted from the PET image, and (c) using heterogeneous activity distribution corresponding to the clinical image following PVC. The three different types of simulated data in each case were reconstructed with two iterations and filtered with a 3D Gaussian postfilter, in order to simulate the intratumor heterogeneous uptake. Heterogeneity in all generated images was quantified using textural feature derived parameters in 3D according to the ground truth of the simulation, and compared to clinical measurements. Finally, profiles were plotted in central slices of the tumors, across lines with heterogeneous activity distribution for visual assessment.Results: The accuracy of the simulated database was assessed against the original clinical images. The PVC simulated images matched the clinical ones best. Local, regional, and global features extracted from the PVC simulated images were closest to the clinical measurements, with the exception of the size zone variability and the mean intensity values, where heterogeneous tumors showed better reproducibility. The profiles on PVC simulated tumors after postfiltering seemed to represent the more realistic heterogeneous regions with respect to the clinical reference.Conclusions: In this study, the authors investigated the input activity map heterogeneity in the GATE simulations of tumors with heterogeneous activity distribution. The most realistic heterogeneous tumors were obtained by inserting PVC activity distributions from the clinical image into the activity map of the simulation. Partial volume effect (PVE) can play a crucial role in the quantification of heterogeneity within tumors and have an important impact on applications such as patient follow-up during treatment and assessment of tumor response to therapy. The development of such a database incorporating patient anatomical and functional variability can be used to evaluate new image processing or analysis algorithms, while providing control of the ground truth, which is not available when dealing with clinical datasets. The database includes all images used and generated in this study, as well as the sinograms and the attenuation phantoms for further investigation. It is freely available to the interested reader of the journal at http://www.med.upatras.gr/oncobase/.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Numerical experiments on the detailed energy conversion and spectrum studies in a corona current sheet

In this paper, we study the energy conversion and spectra in a corona current sheet (CS) by 2.5 dimensional MHD numerical simulations. Numerical results show that many Petschek-like fine structures with slow-mode shocks mediated by plasmoid instabilities develop during the magnetic reconnection process. The termination shocks can also be formed above the primary magnetic island and at the head of secondary islands. These shocks play important roles in generating thermal energy in a corona CS. For a numerical simulation with initial conditions close to the solar corona environment, the ratio of the generated thermal energy to the total dissipated magnetic energy is around 1/5 before secondary islands appear. After secondary islands appear, the generated thermal energy starts to increase sharply and this ratio can reach a value of about 3/5. In an environment with a relatively lower plasma density and plasma β, the plasma can be heated to a much higher temperature. After secondary islands appear, the one-dimensional energy spectra along the CS do not behave as a simple power law and the spectrum index increases with the wave number. The average spectrum index for the magnetic energy spectrum along the CS is about 1.8. The two-dimensional spectra intuitively show that part of the high energy is cascaded to large kx and ky space after secondary islands appear. The plasmoid distribution function calculated from numerical simulations behaves as a power law closer to f(ψ)∼ψ{sup −1} in the intermediate ψ regime. By using η{sub eff}=v{sub inflow}⋅L, the effective magnetic diffusivity is estimated to be about 10{sup 11} ∼ 10{sup 12} m{sup 2} s{sup −1}.

79 ASTRONOMY AND ASTROPHYSICS↗

Attack-resilient algorithms and testbed federation for wide-area protection and control in smart grid

Today’s energy infrastructure is undergoing a massive transformation across all generation, transmission, and distribution systems to provide reliability, efficiency, and sustainability to the power system network. During the past 10 years, there has been a remarkable enhancement in developing close-loop wide-area protection and control (WAPAC) applications, such as automatic generation control (AGC), wide-area protection scheme (WAPS), synchrophasor-based wide-area voltage control system (WAVCS), etc., in the energy management system (EMS) to provide an adequate situational awareness to control center operators during physical disturbances and natural events. The rapid digitalization, unencrypted communication, and numerous data-sharing devices have increased digital access points that have exposed the grid network to countless cyberattacks pertinent to the WAPAC cybersecurity. Given the amount of complex grid network coupled with legacy infrastructures, it is untrustworthy to completely rely on information technology (IT)-based conventional security measures to defend against evolving cyberattacks. Hence, there is a compelling need to go beyond the traditional paradigm of ‘security by obscurity’ and ‘bolt-on’ security measures, and develop a suite of innovative security solutions at multiple layers to transform the current “fault-resilient” grid into “fault and attack-resilient” grid of the future. This dissertation specifically focuses on developing application-specific attack-resilient solutions for critical WAPAC applications in the smart grid. The first part of the study, Attack-Resilient Wide-Area Protection Scheme, presents a cyber-physical attack-resilient system (CPARS) for WAPS, where two anomaly detection methods are proposed to detect multi-level cyberattacks followed by the rules-based intrusion mitigation system (RIMS) to initiate effective mitigation strategies, tailored to these attacks, to restore the normal grid operation after disturbances. The first anomaly detection component, Machine learning-based Anomaly Detection for Centralized WAPS, discusses a novel architecture and methodology for developing a machine-learning-based anomaly detection system (ADS) for the centralized wide-area protection system (CWAPS) by applying the variational mode decomposition (VMD) technique and decision tree (DT) algorithms using cyber logs and synchrophasor measurements. The second component, Multi-Agent-based Anomaly Detection for Decentralized WAPS, presents a two-level hierarchical multi-agent-based decentralized WAPS against the system-aware stealthy cyber-attacks. Finally, for the RIMS, the state transition diagram-based mitigation measures are defined based on the grid operation states and experimental testing and evaluation are presented in a cyber-physical testbed environment. The second part of the study, Attack-Resilient Wide-Area Voltage Control System, addresses the existing vulnerabilities in the wide-area voltage control system (WAVCS), presents the methodology for performing an impact analysis during data integrity attacks, and proposes a data-driven attackresilient system (DARS) by incorporating machine learning-based ADS and rules-based anomaly mitigation system (RAMS) to detect these attacks and provide necessary mitigation actions to maintain the transient voltage stability after disturbances. Finally, the third part, Networked Federation Testbed (NEFTSec) for Smart Grid Cybersecurity, focuses on developing an interconnected federated testbed through a common network like the internet, describes its conceptual architectures, and also presents its design components for experimental testing, validation, and evaluation. Further, it presents the co-simulation interface algorithm (CIA) to facilitate the geographically-dispersed real-time distributed simulation (GD-RTDS) during the cyber-physical federation. In general, this dissertation presents innovative attack-resilient solutions for WAPAC applications that are experimentally tested and evaluated using the PowerCyber testbed and NEFTSec platform to accelerate the transition of state-of-the-art research works to the real-field deployments.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The Quantum Liquid of Alpha Clusters - A Variational Approach

Within the variational approach of Bose liquids we analyze the g.s. energy of charge neutral alpha matter at T = 0. As a prerequisite for such calculation we take from the literature or propose new {alpha}-{alpha} potentials that are particularly suitable for this task, i.e. posses a repulsive core and/or reproduce the low energy scattering data and the resonance properties of the {alpha}-{alpha} system. The alpha matter EOS is then obtained with the HNC method using Pandharipande-Bethe correlation derived variationally in the lowest order expansion of the energy functional or a simple gaussian function with a healing range determined by the normalization of the radial distribution function in the lowest order. We show that saturation is achieved only via repulsive and shallow potentials that are not consistent with the scattering and resonance constraints.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hybrid Simulation of Proton Cyclotron Waves Upstream of Mars Generated by Pickup Ion Beam Distribution

Linear instability analysis as well as a corresponding two‐dimensional hybrid simulation are performed to examine the excitation of the proton cyclotron waves observed upstream of Mars. The waves are believed to be excited by the pickup ions produced from the ionization of the Martian hydrogen exosphere. And previous statistical analysis of wave observations suggested that the waves are mostly related to pickup ion beam velocity distributions. While earlier linear instability analysis of pickup ion beam distributions has mainly been focused on the parallel unstable modes, our analysis reveals that the maximum growth rate occurs at very oblique propagation. The corresponding hybrid simulation confirms the linear analysis results and further demonstrates that the pickup ions are scattered toward an isotropic shell velocity distribution by the waves excited. Interestingly, the waves at oblique propagation gradually damp out and the system is eventually dominated by waves of quasi‐parallel propagation.

79 ASTRONOMY AND ASTROPHYSICS↗

Monte Carlo simulation of simultaneous radiation detection in the hybrid tomography system ClearPET-XPAD3/CT

Using the Geant4 based simulation framework SciFW1, a detailed simulation was performed for a detector array in the hybrid tomography prototype for small animals called ClearPET / XPAD, which was built in the Centre de Physique des Particules de Marseille. The detector system consists of an array of phoswich scintillation detectors: LSO (Lutetium Oxy-ortosilicate doped with cerium Lu{sub 2}SiO{sub 5}:Ce) and LuYAP (Lutetium Ortoaluminate of Yttrium doped with cerium Lu{sub 0.7}Y{sub 0.3}AlO{sub 3}:Ce) for Positron Emission Tomography (PET) and hybrid pixel detector XPAD for Computed Tomography (CT). Simultaneous acquisition of deposited energy and the corresponding time - position for each recorded event were analyzed, independently, for both detectors. interference between detection modules for PET and CT. Information about amount of radiation reaching each phoswich crystal and XPAD detector using a phantom in order to study the effectiveness by radiation attenuation and influence the positioning of the radioactive source {sup 22}Na was obtained. The simulation proposed will improve distribution of detectors rings and interference values will be taken into account in the new versions of detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Simple, Secure, Internet Delivery of MOOSE-based Applications

Application packaging and distribution are the final steps for delivering software to end-users; both are frequently neglected when creating scientific software. Commercial businesses rely on electronic distribution systems that have rendered disk drives obsolete. Still, national laboratories continue to rely heavily on removable media to distribute and limit access to controlled applications. With increasing concerns of unauthorized copying of sensitive applications, a modern distribution system that utilizes cryptographically secure communication and authentication protocols has been developed. This new distribution system will secure the chain of custody for nuclear software while simultaneously simplifying access to these tools. This report summarizes four primary advancements made toward the secure distribution of Nuclear Energy Advanced Modeling and Simulation (NEAMS)-developed, Multiphysics Object Oriented Simulation Environment (MOOSE)-based applications: application installation, package distribution, automated package building, and distribution of documentation. NEAMS is currently developing more than ten separate applications based on the open-source MOOSE Framework. Distribution of these applications has primarily been accomplished by distributing source code, with end-users compiling the applications themselves. This work created a mechanism where MOOSE applications can be installed in a similar way to any other software. This allows both administrators and end-users simplified access to runnable executables. With this new installation capability, it was then possible to rethink distribution. A new, secure capability for delivering MOOSE-based applications over the internet has been created. This system requires unique cryptographic tokens for authentication, greatly securing the custody chain for software. Once granted access, installation of any NEAMS code can be accomplished with these terminal commands: "conda install ncrc" "ncrc install ncrc-bison." After these two commands (and authenticating) the BISON application will be securely down- loaded from Idaho National Laboratory (INL)’s servers, installed, and ready to use. To enable this new distribution capability to be successful, the open-source Continuous Integration, Verification, Enhancement, and Testing (CIVET) Continuous Integration (CI) capability was augmented to add Continuous Delivery (CD). CD enables the automated building and packaging of MOOSE-based applications as they are modified by development teams, ensuring that our customers can obtain up-to-date versions of the software at any time. The need for instruction on how to use these applications was addressed through modifications to the MOOSE documentation system. The MooseDocs capability, which enables robust documentation of MOOSE-based applications, has been extended to allow both for the installation of documentation and the packaging of documentation with installed applications. Together, these enhancements form the core of a new, secure distribution mechanism for nuclear simulation tools. In concert with the Nuclear Computational Resource Center (NCRC), NEAMS- developed applications will now be straightforward to obtain securely.

97 MATHEMATICS AND COMPUTING↗

Emittance study of a 28 GHz electron cyclotron resonance ion source for the Rare Isotope Science Project superconducting linear accelerator

A 28 GHz electron cyclotron resonance (ECR) ion source is being developed for use as an injector for the superconducting linear accelerator of the Rare Isotope Science Project. Beam extraction from the ECR ion source has been simulated using the KOBRA3-INP software. The simulation software can calculate charged particle trajectories in three dimensional complex magnetic field structures, which in this case are formed by the arrangement of five superconducting magnets. In this study, the beam emittance is simulated to understand the effects of plasma potential, mass-to-charge ratio, and spatial distribution. The results of these simulations and their comparison to experimental results are presented in this paper.

43 PARTICLE ACCELERATORS↗