Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulated flux procedure”

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.

61 records · Page 4

Double-differential Measurements of Mesonless Charged-current Muon Neutrino Interactions on Argon with Final-State Protons using the MicroBooNE Detector

This note reports the first double-differential measurements of charged-current $ν_µ$ scattering on argon leading to final states containing zero mesons and one or more protons. This event topology (hereafter abbreviated as CC0$πNp$) is the most common at the neutrino energies typically produced by the Fermilab Booster Neutrino Beam. A detailed understanding of neutrino-argon scattering in the CC0$πNp$ channel is therefore crucial for the success of the precision neutrino oscillation analyses planned for the Short-Baseline Neutrino (SBN) program. This remains true for the upcoming Deep Underground Neutrino Experiment (DUNE), but the higher mean neutrino energy used there will ensure that more inelastic reaction modes, such as single pion production, will also play a major role. The analysis described in this note builds on a previous MicroBooNE study of the CC0$πNp$ channel which obtained the first single-differential cross-section measurements on an argon target. Since that foundational work, significant improvements have been made to MicroBooNE’s simulation software, event reconstruction algorithms, and procedure for calculating systematic uncertainties. When combined with a larger dataset (corresponding to a beam exposure of 6.79 × 10 20 protons-on-target versus 1.60 × 10 20 in Ref. [3]), these enhancements allow the important CC0$πNp$ channel to be studied in more detail. This note begins with a description of the data and simulation samples used as input to the analysis. Section 3 then defines the CC0$πNp$ signal event topology, and Section 4 describes a set of selection criteria designed to identify these events in MicroBooNE data. Binning schemes are then defined in Section 5 for two double-differential measurements of event rates. The first of these considers the momentum and scattering cosine of the outgoing muon, while the second reports the same observables for the leading proton, i.e., the final-state proton with the largest momentum. After a discussion of systematic uncertainties in Section 6, the note concludes by comparing the predictions of MicroBooNE Monte Carlo (MC) simulations to the measured double-differential distributions. These results will form the basis for a future extraction of flux-averaged double-differential CC0$πNp$ cross sections that will be immediately comparable to the theoretical predictions of multiple neutrino event generators. The selection described herein may also be used to study various other observables in CC0$πNp$ events, including those which are sensitive to correlations between leptonic and hadronic kinematics in the final state.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Low-Cost Heliostat for High-Flux Small-Area Receivers (Final Technical Report)

This project analyzed a two-stage heliostat concept consisting of a tracking stage and a concentrating stage. The tracking stage uses mirrors mounted on a common drive that move to track the sun. The concentrating stage consists of stationary mirrors that each have a unique angle to direct rays towards a small-area, high-flux, point-focused receiver. By splitting the collection and concentrating process into two stages, multiple small, inexpensive mirrors can share a structure and be controlled by a single drive in the tracking stage. The project effort developed modeling techniques that were specifically relevant to this two-stage heliostat concept. Both field-level and unit-level models were developed. The field-level model does not explicitly consider unit-level losses which are predicted by the unit-level model and then integrated into the field-level model through a correlation referred to as an efficiency modifier. This approach is referred to as the two-model approach; the development and demonstration of this two-model approach for a multi-stage heliostat technology is a key outcome of this work. The field-level model is used to design a field that hits a specific design day power given a set of heliostat design parameters. An oversized field is simulated and then heliostat units are removed based on their annual energy production in order to generate the highest performing field. The field reduction procedure fits a smooth curve fit to annual energy production as a function of position in the field which has the effect of reducing the noise that is otherwise caused by the Monte Carlo ray tracing technique. This approach is referred to as the annual energy fit method and substantially reduces computational run time for a given field level modeling accuracy. The annual energy fit approach enables the selection of a properly sized, high-performing field using orders of magnitude fewer rays than would otherwise be possible and the development of this approach is a second key outcome of this work. These models are used within a genetic optimization algorithm in order to optimize the geometric parameters associated with a heliostat in order to achieve the lowest cost per unit of collected design day power. The cost modeling that underlies the optimization is a simple, scaling type analysis backed up by a much more detailed Design for Manufacture and Assembly (DFMA) analysis. Although the figure of merit used for optimization was not cost per mirror area, this metric is reasonable to use as a means of comparison. The optimally designed 500 kW design has a tracking mirror specific cost of $181.85/m 2 , which is significantly larger than the target value and also larger than the current state of the art. The cost of the torque-tube type linkages contributed substantially to the overall cost. Based on this observation, potentially attractive alternative design configuration utilizing a capstan type actuation system should be investigated. Finally, NREL compared the performance of the two-stage heliostat to the performance of a focused and different sized flat conventional heliostats and showed that, as expected, additional losses versus the convention heliostat caused by a worse cosine efficiency, two stages of reflection, and interstage interactions. The two-stage heliostat requires around 75% more reflective area than a flat 1x1 meter conventional heliostat (similar to a focused heliostat) and 40% more than a flat 2x2 meter conventional heliostat.

14 SOLAR ENERGY↗

The gallium solar neutrino capture cross section revisited

Solar neutrino flux constraints from the legacy GALLEX/GNO and SAGE experiments continue to influence contemporary global analyses of neutrino properties. The constraints depend on the neutrino absorption cross sections for various solar sources. Following recent work updating the 51 Cr and 37 Ar neutrino source cross sections, we reevaluate the 71 Ga solar neutrino cross sections, focusing on contributions from transitions to 71 Ge excited states, but also revising the ground-state transition to take into account new 71 Ge electron-capture lifetime measurements and various theory corrections. The excited-state contributions have been traditionally taken from forward-angle (𝑝, 𝑛) cross sections. Here we correct this procedure for the ≈ 10%–20% tensor operator contribution that alters the relationship between Gamow-Teller and (𝑝, 𝑛) transition strengths. Using state-of-the-art nuclear shell-model calculations to evaluate this correction, we find that it lowers the 8 B and hep neutrino cross sections. However, the addition of other corrections, including contributions from near-threshold continuum states that radiatively decay, leads to an overall increase in the 8 B and hep cross sections of ≈ 10% relative to the values recommended by Bahcall. Uncertainties are propagated using Monte Carlo simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Estimating List-Mode Data Sensitivities to Nuclear Data with MCNP6

Nuclear data are a vital component of predictive simulations used in applications like experiment design, stockpile stewardship, nuclear nonproliferation/safeguards, health physics, and criticality safety. A singular simulation requires the coalescence of different areas of nuclear data such as cross sections, angular distributions, and energy distributions of emitted neutrons for different materials and energy ranges. Improving nuclear data and thus reducing the uncertainty in simulated parameters could enable smaller, better-informed safety factors and ultimately reduce operational and procedural costs. There is a constant effort to garner a better understanding of the physical quantities represented by nuclear data through experiments. Integral experiment benchmarks use simulated and measured results to validate current nuclear data values. In the past, benchmarks primarily focused on the effective multiplication factor (k eff ); however, this limited scope has caused compensating errors and areas of nuclear data that lack validation. Compensating errors are inaccuracies in nuclear data that are obfuscated by cancellation when observing integrated values such as k eff . Diverse integral benchmark experiments that look for quantities of interest other than k eff and include multiple responses minimize the possibility of compensating errors and provides validation to areas of nuclear data previously lacking experimental validation. Benchmark experiments can be optimized during the design process to be highly dependent on specific areas of nuclear data. The dependence of a response in an experiment to a specific area/type of nuclear data is defined as sensitivity. A larger sensitivity means that nuclear data uncertainties will play a larger role in the response(s) resulting in larger bias. Currently, the sensitivity capabilities of the Monte Carlo N-Particle (MCNP ®1 ) transport code are limited to responses of k eff and tallied values (e.g., flux, surface current). As a part of the EUCLID project, this work explores estimating list-mode nuclear data sensitivities that can be used to design experiments aimed to constrain and reduce compensating errors in nuclear data by focusing on responses other than k eff . Tallied values are ideal quantities that are estimated with detectors during experiments. List-mode data (a list of neutron collection times) are the direct output of detector systems in subcritical neutron noise experiments. Expanding MCNP sensitivity capabilities to include the sensitivity of responses estimated from list-mode data, such as the prompt neutron decay constant (α) and multiplicity estimates (S and D), enables more direct comparison of simulated and measured experimental quantities. Additionally, deterministic tools such as SENSMG are capable of obtaining sensitivities to a wide variety of responses; however, these tools cannot handle complex geometries due to the assumptions made in discretizing the phase-space variables of the Boltzman transport equation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Deep potential molecular dynamics simulations of ion-enhanced etching of silicon by atomic chlorine

The continued development of plasma-assisted processing techniques requires a fundamental understanding of plasma-surface interactions. Molecular dynamics (MD) simulations have been employed to complement experimental studies and better understand the properties of such systems. Recently, machine learning (ML) methods have enabled the development of ab initio-based interatomic potentials, which can be generalized to complex combinations of multiple atom types. In this work, we use ML potentials developed using the Deep Potential Molecular Dynamics (DeepMD) framework to provide a model of ion-enhanced etching of Si by Cl atoms. We demonstrate the importance of proper selection of the training data set to the accuracy of the DeepMD model and compare our results to MD results using empirical potentials, as well as to experimental measurements. Exposure of undoped Si at 300 K to thermal Cl atoms yields a steady-state Cl coverage of 1.25 monolayers, which is slightly lower than the value obtained in previous experimental studies. Predictions of Si etch yields by simultaneous Cl atom and Ar + ion impacts as a function of ion energy, neutral to ion flux ratio, and angle of incidence of the ions are in reasonably good agreement with classical MD results and experimental measurements. Finally, etch yields and SiCl x mixed layer thicknesses during simultaneous bombardment of the Si(100) surface by Cl atoms and Cl + ions are in good agreement with experimental data. In conclusion, the present work is a necessary condition for the extension of the DeepMD procedure to more complex systems of interest in plasma-surface interactions.

Artificial neural networks↗

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Revisions of ORNL 188 W Process Based on Nonradiological Experiments

Tungsten-188 is in widespread use in the 188 W(t 1/2 = 69 d)/ 188 Re(t 1/2 =16.9 h) biomedical generator. Oak Ridge National Laboratory (ORNL) has been providing this product to the world since 1999. At ORNL, 188 W is produced via irradiation in ORNL’s High Flux Isotope Reactor (HFIR). Enriched 186 W targets in the form of sintered metallic pellets or rings achieve a compact loading in the irradiation vessel, providing a high yield per unit target. The enrichment of the target is >90% 186 W and this isotope undergoes double neutron capture to produce the desired 188W product. While 188 W is produced by neutron bombardment, 191 Os(t 1/2 = 15.4 d) is simultaneously produced as a byproduct requiring separation from 188 W by post irradiation treatment. In the current processing pathway, the irradiated W metal pellets or rings are first converted into an oxide form of WO 3 by heating the irradiated W metal target at 750°C in a quartz reaction vessel inside a vertical furnace under a constant flow of air. During heating, W metal reacts with oxygen in the air to produce WO 3 , which is soluble in 6 M NaOH for later purification process. This oxidation process also converts 188 Os (the decay daughter of 188 W) and 191 Os (the irradiation produced byproduct) into OsO 4 , a highly volatile and toxic gas. The gaseous effluents driven from the quartz reaction vessel are passed through a scrubbing array to remove OsO 4 before the air is discharged from the process. This heterogeneous oxidation method simultaneously achieves two goals: 1). converting metal target to a soluble oxide form and 2). separating volatile OsO 4 away from the solid WO 3 product, although harmful 191 OsO 4 is unfavorable but being taken care of by the down road scrubbing array. In the past twenty years the existing OsO 4 scrubbing array served well in preventing OsO 4 from being released into the environment, until July of 2020 when a minor amount of 191 Os was found to have been released into the environment which resulted in a standby of building 4501 from July 2020 to February 2021.Since October 20 of 2020 a team of researchers from groups of RSTD and NEFD were organized to work for the following tasks: (1) Fully understand the root causes of the July incident of 191 Os release and create a redesigned 191 Os scrubbing array to eliminate the possibility of 191 Os release in future 188 W processes; (2) In Phase-1 tests, select a correct air (O 2 ) flow rate in 188 W process, that provides sufficient oxygen for the oxidation of tungsten metal, while still allowing sufficient resonance time for 191 Os removal by the scrubbers; (3) In Phase-2 tests, confirm the efficiency of the redesigned scrubbing array to absorb excessive amount of non-rad OsO 4 (generated from Os metal powder) at the selected air flow; (4) In Phase-3 tests, simulate the hot cell 188 W process by heating pressed non-rad W-Os pellets in a quartz vessel of new design and the new array under operation parameters selected in Phase-1 and -2; (5) Complete a new operation procedure for 188 W process in hot cell and a TM report as a summary of the Os Mitigation Project. This TM report summarizes the investigation on root causes from technical aspects in July 191 Os release and the correspondent improvement towards the redesigned 191 Os scrubbing system. The report presents the efforts made on selection of optimal operational parameter with the new scrubbing array and the quartz vessel of new design, based on experiment data obtained in Phase-1, -2 and -3 tests. Further renovations for 188 W process and additional improvement of W target treatment will also be discussed.

07 ISOTOPE AND RADIATION SOURCES↗