SEARCH · Engineering Papers

Results for “merit function”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2020

ELM control optimization for various ITER scenarios based on linear and quasi-linear figures of merit

For the purpose of better understanding type-I edge localized mode (ELM) control in ITER with resonant magnetic perturbation (RMP) fields, the plasma response to RMP is computed by a resistive full magneto-hydrodynamic model in toroidal geometry. Five scenarios designed for ITER are considered, ranging from the pre-nuclear to nuclear phases. The plasma response to RMP is quantified by the plasma surface displacement near the X-point of the divertor plasma and at the outboard mid-plane. The optimal coil configurations between two high- Q deuterium-tritium (DT) scenarios (at the same plasma current of 15 MA and the same magnetic field of 5.3 T but different fusion gains, Q = 5 and 10) are predicted to be similar. For the other ITER scenarios with similar edge safety factor q 95 ~ 3 to that of the baseline scenario, the optimal coil phasing is also similar. The optimization results are different for a half-current full-field (7.5 MA/5.3 T) scenario, largely due to the difference in q 95 . The RMP coil currents are also optimized to tailor the core vs edge toroidal torques exerted by the 3D RMP fields on the plasma column. Torque optimization, with various objective functions proposed in the study, is useful for minimizing the side effects of RMP on the plasma core flow in ITER, while still maintaining the ELM control capability. Full utilization of three rows of ELM control coils in ITER is found to be essential to ensure both flexibility and robustness of ELM control, in terms of both linear and quasilinear plasma responses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Watching Polarons Move in the Energy and Frequency Domains Using Color Impedance Spectroscopy

The hybrid electronic–ionic transport property of π-conjugated polymers enables new (opto)electrochemical device constructs for energy conversion and storage and biosensing applications. One major challenge is separating the energy and frequency dependence of Faradaic events–those involving charge transfer and the redox processes of the conjugated backbone–from the non-Faradaic components, such as ionic motion. Herein, we combine optical spectroscopy with electrochemical impedance spectroscopy (EIS) to resolve the frequency response of ionic–electronic coupling as a function of electrochemical doping potential. First, using EIS, we identify two different frequency regimes resulting in potential-dependent capacitive elements on the order of ~10 μF/cm 2 in a high-frequency regime and ~50–150 μF/cm 2 in a low-frequency regime. Given the larger magnitude and greater potential dependence, we posit that polaronic motion is more likely to occur at low frequencies (<1 kHz) and overlaps with ionic motion. The use of color impedance spectroscopy (CIS) enables observation of polaronic motion with frequency modulation. We observe that higher doping potentials show a greater motion of polarons above the DC-bias baseline concentration for onset in electrochemical doping, but all potentials considered demonstrate a critical frequency at which the polaronic motion is “frozen” (~40 Hz). Furthermore, this critical information obtained from CIS in highly dielectric environments offers a unique figure of merit for future studies on electronic–ionic coupling by which to compare across polymer/electrolyte interfaces, including the role of a charge-supporting electrolyte, a solvent, and alternative Faradaic processes (e.g., electrocatalysis).

36 MATERIALS SCIENCE↗

DOE OSTI2024

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↗

DOE OSTI2023

A Comparative Evaluation and Selection of High-Temperature Heat Exchangers for Application to Integrated Energy Systems

The following report aims to create a refined and well-structured method for comparatively evaluating heat exchanger technologies for integrated energy systems that caters customers’ specific needs while meeting engineering requirements. For the evaluation, this study elevates previous evaluation metrics, enhances the knowledge base via literature and market surveys, and identifies the figures of merit with robust rationales to enhance the quality of decisions made throughout the proposed heat exchanger evaluation process. The heat exchanger designs evaluated as part of the case study are shell and tube heat exchangers, printed circuit heat exchangers, plate heat exchangers, spiral heat exchangers, and heat pipe heat exchangers. The information presented in this report is meant for industries interested in making a preliminary screening process to identify the most suitable heat exchanger design for their application of interest.

42 ENGINEERING↗

DOE OSTI2023

A Comparative Evaluation and Selection of High-Temperature Heat Exchangers for Application to Integrated Energy Systems

42 - ENGINEERING↗

DOE OSTI2021

Ultrawide Temperature Range Super-Invar Behavior of R 2 (Fe,Co) 17 Materials ( R = Rare Earth)

Super Invar (SIV), i.e., zero thermal expansion of metallic materials underpinned by magnetic ordering, is of great practical merit for a wide range of high precision engineering. However, the relatively narrow temperature window of SIV in most materials restricts its potential applications in many critical fields. Here, we demonstrate the controlled design of thermal expansion in a family of R 2 (Fe,Co) 17 materials (R=rare Earth). We find that adjusting the Fe-Co content tunes the thermal expansion behavior and its optimization leads to a record-wide SIV with good cyclic stability from 3–461 K, almost twice the range of currently known SIV. In situ neutron diffraction, Mössbauer spectra and first-principles calculations reveal the 3d bonding state transition of the Fe-sublattice favors extra lattice stress upon magnetic ordering. On the other hand, Co content induces a dramatic enhancement of the internal molecular field, which can be manipulated to achieve “ultrawide” SIV over broad temperature, composition and magnetic field windows. Furthermore, these findings pave the way for exploiting thermal-expansion-control engineering and related functional materials.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Harnessing Condorcet Methods to Improve Decision-making Based on Ranked Data

Decision-makers must often choose between multiple alternatives based on their relative merits across a variety of criteria. These multiple-criteria decision-making (MCDM) problems must be approached in an objective, measurable, and transparent fashion to obtain well-supported outcomes. The challenge of using sets of ranked data to identify the optimal choice is not unlike the challenge of using ranked ballots to identify the winner of an election in a preference voting system. As such, the methods of rank aggregation used in elections can be applied to help resolve MCDM problems. Previous research into the design of elections and ballots has yielded many algorithms with well-understood properties. A subset of these, called Condorcet methods, reliably identify the “Condorcet winner,” if it exists, giving the result that would defeat any other in a pairwise comparison. That, in addition to several other desirable qualities, makes Condorcet methods like the Schulze method, Ranked Pairs, and Copeland’s method, effective and scalable means for tackling MCDM problems. These rank aggregation methods have been implemented in JavaScript functions for incorporation into a web-based MCDM decision tool. These algorithms can be applied in a wide variety of contexts (for example, determining the “best” environmental remediation method or selecting a subcontractor) to facilitate effective decision-making. The performance of the JavaScript Condorcet prototype tool was compared across the implemented algorithms and with both a non-Condorcet rank aggregation method as well as an implementation of the Simple Multi-Attribute Rating Technique (SMART) algorithm.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2022

Dynamic calibration of differential equations using machine learning, with application to turbulence models

We present a methodology for calibration of parametric ordinary and partial differential equation models, using off-the-shelf software for back-propagation in Neural Networks (NN). As a prototypical example, we consider calibration of a Reynolds-averaged Navier-Stokes (RANS) turbulence closure model, against ground truth data from direct numerical simulations (DNS) of two different turbulent flows. Numerical time integration is represented as a custom NN, where only the RANS model parameters are trainable. A loss function is defined to quantify the mismatch between the NN prediction and the ground truth over a predefined, finite time integration window. This loss function is then minimized using a gradient descent method utilizing the back-propagation algorithm. Furthermore, this dynamic approach to training is to be contrasted with a static approach, wherein a least square regression estimate for parameters is obtained in the limit of an infinitesimal time integration window. In a first test of static and dynamic approaches against ground truth data generated by the model, the former proves to be significantly faster and more accurate than the latter at recovering the parameters. When both calibration approaches are tested against DNS data, for which it is known that the model cannot achieve a perfect fit, the static approach yields a good prediction only for short times, while the dynamic approach results in physical and stable predictions over the entire integration window. After optimization of the dynamic approach for time step, spatial resolution, stability, and physics-based constraints, we obtain a 50% improvement of outcomes over those obtained from the existing, manually calibrated set of parameters, demonstrating the merits of this systematic and automated procedure.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Deterministic modeling of hybrid nonlinear effects in epsilon-near-zero thin films

In nonlinear optics, significant effort is concentrated on improving the strength and efficiency of interactions; however, experimentally investigating nonlinear materials is a complex, time-consuming, and costly investment. Moreover, it is often challenging to isolate, study, and optimize material parameters in an experiment due to complexities in the growth process. Recently, epsilon-near-zero materials have received a great deal of attention as promising nonlinear optical materials, but like many up-and-coming materials, the ability to explore and optimize their properties has been challenging. Here, we establish a framework to rapidly evaluate the performance of nonlinear epsilon-near-zero materials for both inter- and intraband effects in silico, requiring only an energy-momentum (E-k) diagram, linear optical properties, and experimental conditions. Measured nonlinear reflection and transmission in gallium-doped zinc oxide films are compared to the numerical framework for both intra- and interband excitation to verify accuracy across wavelength and irradiance while two figures of merit (FoMs) are introduced to quickly evaluate the performance of films without a full numerical framework. This capability is used to predict the performance of highly doped gallium nitride, cadmium oxide, zinc oxide, and indium tin oxide films, and efficient intra- and interband operation conditions are identified. Through this numerical framework and the FoMs, the exploration of unstudied epsilon-near-zero materials is enabled without the need for a nonlinear experiment, thereby accelerating the search for more efficient nonlinear materials and excitation conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Improving the NOvA 3-Flavour Neutrino Oscillation Analysis

NOvA is a long-baseline neutrino experiment studying neutrino oscillations, a quantum mechanical interference phenomena where the observed neutrino flavour differs from that measured earlier, stemming from neutrino mass and flavour states mixing. NOvA consists of two functionally identical tracking calorimeter detectors deployed in the Fermilab NuMI beam. Both detectors are placed 14.6 mrad off the beam axis to achieve a narrow energy peak at 1.8 GeV at an oscillation maximum. The NOvA 3-flavour oscillation analysis measures the neutrino oscillation parameters sin2θ23 and ∆m232 as well as sets limit to δCP , the octant of θ23 and the sign of ∆m32. The event selection for 3-flavour neutrino oscillation analysis ensures the maximum quantity of signal is made available and a minimal amount of background is present. The current disappearance analysis selection has an efficiency of 80% for selecting νμ CC events, and with some improvements additional events could be recovered into the analysis to improve the sensitivity to the aforementioned oscillation parameters. In an effort to recover the currently rejected signal event to the analysis, these events were trained in a classification neural network. The aim of the network was to divide the data into signal (νμ CC events) and background (NC and νe CC events). The highest performing network gave an additional figure of merit gain of 2.34 increasing the sensitivity by 4.3% in effective POT equal to 33 days of additional data taking. This was compared to changing the current particle identification event selection cuts, the best result out of the tested cut combinations gave an additional FOM of 3.36 equal to 3.7% increase in effective POT equivalent to 28 days of data taking.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Demonstration of Energy-Resolved γ-Ray Detection at Room Temperature by the CsPbCl 3 Perovskite Semiconductor

The detection of gamma-rays at room temperature with high-energy resolution using semiconductors is one of the most challenging applications. The presence of even the smallest amount of defects is sufficient to kill the signal generated from gamma-rays which makes the availability of semiconductors detectors a rarity. Lead halide perovskite semiconductors exhibit unusually high defect tolerance leading to outstanding and unique optoelectronic properties and are poised to strongly impact applications in photoelectric conversion/detection. Here we demonstrate for the first time that large size single crystals of the all-inorganic perovskite CsPbCl 3 semiconductor can function as a high-performance detector for gamma-ray nuclear radiation at room temperature. CsPbCl 3 is a wide-gap semiconductor with a bandgap of 3.03 eV and possesses a high effective atomic number of 69.8. We identified the two distinct phase transitions in CsPbCl 3 , from cubic (Pm-3m) to tetragonal (P4/mbm) at 325 K and finally to orthorhombic (Pbnm) at 316 K. Despite crystal twinning induced by phase transitions, CsPbCl 3 crystals in detector grade can be obtained with high electrical resistivity of similar to 1.7 X 10 9 Ω∙cm. The crystals were grown from the melt with volume over several cubic centimeters and have a low thermal conductivity of 0.6 W m -1 K -1 . The mobilities for electron and hole carriers were determined to similar to 30 cm 2 /(V s). Using photoemission yield spectroscopy in air (PYSA), we determined the valence band maximum at 5.66 +/- 0.05 eV. Under gamma-ray exposure, our Schottky-type planar CsPbCl 3 detector achieved an excellent energy resolution (similar to 16% at 122 keV) accompanied by a high figure-of-merit hole mobility-lifetime product (3.2 x 10 -4 cm 2 /V) and a long hole lifetime (16 mu s). The results demonstrate considerable defect tolerance of CsPbCl 3 and suggest its strong potential for gamma-radiation and X-ray detection at room temperature and above.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Automated quantum error mitigation based on probabilistic error reduction

Current quantum computers suffer from a level of noise that prohibits extracting useful results directly from longer computations. The figure of merit in many near-term quantum algorithms is an expectation value measured at the end of the computation, which experiences a bias in the presence of hardware noise. A systematic way to remove such bias is probabilistic error cancellation (PEC). PEC requires a full characterization of the noise and introduces a sampling overhead that increases exponentially with circuit depth, prohibiting high-depth circuits at realistic noise levels. Probabilistic error reduction (PER) is a related quantum error mitigation method that systematically reduces the sampling overhead at the cost of reintroducing bias. In combination with zero-noise extrapolation, PER can yield expectation values with an accuracy comparable to PEC.Noise reduction through PER is broadly applicable to near-term algorithms, and the automated implementation of PER is thus desirable for facilitating its widespread use. To this end, we present an automated quantum error mitigation software framework that includes noise tomography and application of PER to user-specified circuits. We provide a multi-platform Python package that implements a recently developed Pauli noise tomography (PNT) technique for learning a sparse Pauli noise model and exploits a Pauli noise scaling method to carry out PER.We also provide software tools that leverage a previously developed toolchain, employing PyGSTi for gate set tomography and providing a functionality to use the software Mitiq for PER and zero-noise extrapolation to obtain error-mitigated expectation values on a user-defined circuit.

McDonough, Benjamin↗

DOE OSTI2024

Superconductivity in CH 4 and BH – 4 containing compounds derived from the high-pressure superhydrides

Inspired by the synthesis of the high-pressure Fm3m LaH 10 superconducting superhydride, systematic density functional theory (DFT) calculations are performed to study ternaries that could be derived from it by replacing two of the hydrogen atoms with boron or carbon and varying the identity of the electropositive element. Though many of the resulting alkali-metal and alkaline-earth MC 2 H 8 phases are predicted to be dynamically stable at mild pressures, their superconducting critical temperatures (T c s) are low because their metallicity results from the filling of an electride-like band. Substitution with a trivalent element leads to phases with substantial metal d- character at the Fermi level whose T c s are typically above 40 K. Here, among the MB 2 H 8 phases examined, KB 2 H 8 , RbB 2 H 8 and CsB 2 H 8 are predicted to be dynamically stable at very mild pressures, and their stability is rationalized by a DFT-Chemical Pressure analysis that elucidates the role of the M atom size. Quantum anharmonic effects strongly affect the properties of KB 2 H 8 , the highest predicted T c compound, near 10 GPa, but molecular dynamics simulations reveal it would decompose below its T c at this pressure. Nonetheless, at ca. 50 GPa KB 2 H 8 is predicted to be thermally stable with a superconducting figure of merit surpassing that of the recently synthesized LaBeH 8 .

36 MATERIALS SCIENCE↗

DOE OSTI2024

Achieving environmental stability in an atomically thin quantum spin Hall insulator via graphene intercalation

Atomic monolayers on semiconductor surfaces represent an emerging class of functional quantum materials in the two-dimensional limit — ranging from superconductors and Mott insulators to ferroelectrics and quantum spin Hall insulators. Indenene, a triangular monolayer of indium with a gap of ~ 120 meV is a quantum spin Hall insulator whose micron-scale epitaxial growth on SiC(0001) makes it technologically relevant. However, its suitability for room-temperature spintronics is challenged by the instability of its topological character in air. It is imperative to develop a strategy to protect the topological nature of indenene during ex situ processing and device fabrication. Here we show that intercalation of indenene into epitaxial graphene provides effective protection from the oxidising environment, while preserving an intact topological character. Our approach opens a rich realm of ex situ experimental opportunities, priming monolayer quantum spin Hall insulators for realistic device fabrication and access to topologically protected edge channels.

36 MATERIALS SCIENCE↗

DOE OSTI2020

STRUCTURE AND PROPERTIES STUDY ON ENERGY MATERIALS: THERMOELECTRIC MATERIAL TETRAHEDRITE AND LITHIUM ION CONDUCTOR

Development of efficient energy materials is critical in order to ease the energy demand and reduce our dependence of fossil fuel. Thermoelectric materials are promising due to their capability of generating electrical power by recovering waste heat. The performance of thermoelectric materials is quantified by a dimensionless figure of merit zT, which depends on their properties such as electrical conductivity, Seebeck coefficient and thermal conductivity. Tetrahedrites, a copper antimony sulfosalt mineral, typified by Cu 12-x MxSb 4 S 13 , where M is a transition metal element such as Ni, Zn, Fe or Mn, have great potential for thermoelectric application due to their relatively high zT (close to 1 at 700 K), earth-abundance, environmental friendliness, favorable electrical properties, and most importantly intrinsic low lattice thermal conductivity (less than 1 W m -1 K -1 ) in wide temperature. In addition to energy recovery, reliable energy storage devices are also emerging to relieve the energy demand and improve the efficiency of consuming energy resources. Lithium-ion batteries are known to be reliable and successful electrochemical energy storage devices and appliable in various aspects, including laptops, smartphones and electrical vehicles. Lithium phosphorous oxynitride (LiPON) are widely used as thin-film solid-state electrolytes in Li-ion battery, which is the only demonstrated solid-state electrolyte that is quite stable in direct contact with Li metal at potentials from 0-5 V. However, the structure of LiPON, the effects of N doping, and the origin of its good electrochemical stability remains inconclusive. In this thesis, reliable modeling techniques accompanied with experimental tools, are applied to study the thermoelectric material tetrahedrite and the ionic conductor LiPON, in order to study their structural and dynamical properties. Accurate and efficient density-functional theory (DFT) and density-functional tight-binding (DFTB) methods, combined with molecular dynamics (MD) simulations are utilized in order to investigate the structures and properties of these energy materials. The incoherent and coherent atomic dynamics study of tetrahedrite Cu 10.5 NiZn 0.5 Sb 4 S 13 provides the origin of softening upon cooling by investigate the motion of Cu12e at different temperatures. The dynamic structure factors in the longitudinal and transverse direction will also be discussed. The Cu movement of Cu-rich tetrahedrite Cu 14 Sb 4 S 13 is revealed by Cu self-diffusivity, nuclear density map and “nudged elastic band” (NEB). Moreover, we investigate the effect of simulation cell size and basis sets on the DFT-based MD simulation results using tetrahedrite Cu1 0 Zn 2 Sb 4 S 13 thermoelectric as a model material, showing the advantage of larger cell by accessing smaller Q range. In addition, the low-temperature structural properties of Cu 12 Sb 4 S 13 is measured by neutron diffraction, which indicates that no cubic to tetragonal transition occurs at metal-semiconductor transition (MST) temperature. Thermoelectric properties such as Seebeck coefficient, electrical resistivity and electrical thermal conductivity will also be investigated. DFTB method is implemented to study the structure and transport properties of Li 3 PO 4 and LiPON, while the exploration of N doping effect is included. Lastly, the LiPON/Li interphase will be revealed in order to study the origin of electrochemical stability.

Li, Junchao↗

DOE OSTI2024

High-Fidelity Arc-Discharge Model for Hydrogen-Plasma-Smelting-Reduction of Iron Ore

Electrification and use of renewable hydrogen is currently a necessity for decarbonizing the iron-and-steel industry. In this regard, hydrogen plasma smelting reduction (HPSR) is a novel pathway that is being explored for reduction of iron ore. HPSR provides several decarbonization merits compared to conventional blast furnaces. Firstly, the use of renewable hydrogen drastically reduces the CO2 emissions compared to the use of coke. Secondly, renewable electricity in the form of a thermal plasma for making reactive hydrogen species (radicals, ions) are more efficient at reducing iron ore compared to neutral H2. Thirdly, a molten product compatible with downstream processes is obtained from the intense heat transfer from the plasma. However, the scale-up of this technology requires fundamental exploration of hydrogen plasma dynamics and its interaction with complex solid material that include phase changing iron-ore and slag. In this work, we present a first principles continuum scale model for thermal plasmas in Ar/H2 gas mixtures typically used for HPSR. The thermal plasma governing equations for mass, momentum and energy with Lorentz force and Joule heating source terms are solved along with electromagnetic equations for electrostatic and magnetic vector potential. Our solver will be based on Pele, a suite of reacting flow solvers designed for advanced scientific computing architectures (Henry De Frahan et al., Proceedings of SIAM Parallel Processing, 13-25, 2024), and will utilize adaptive mesh generation for enhanced resolutions at locations of intense physicochemical interactions. This study will present the impact of Ar to H2 ratios on excited/dissociated hydrogen species concentrations, plasma temperature and conductivity along with the impact of outgassed species (water, metal vapor, O, OH radicals) from ore surface on gas phase chemistry. Furthermore, the heat and species flux to the surface will be quantified as a function of applied voltages in a transferred arc configuration.

hydrogen plasma↗

DOE OSTI2025

Implementation of the D1S Methodology for Shutdown Dose Rate Calculations in the OpenMC Monte Carlo Particle Transport Code

We present an implementation of the direct one-step (D1S) methodology for shutdown dose rate (SDR) calculations in the OpenMC Monte Carlo particle transport code. In addition to being the first fully open-source D1S implementation, it is also the first to require no ad hoc source code or nuclear data library modifications. The code can seamlessly switch between production of prompt and decay photons based on a user input parameter, and the decay data needed for decay photon generation are made available through a depletion chain file, which is already used for OpenMC’s built-in depletion/activation solver. A set of Python functions significantly eases the burden of computing and applying time correction factors needed to properly account for the time dependence of radionuclide activity. To assess the accuracy of the D1S implementation, SDR calculations have been carried out for three problems: a prism of iron irradiated by 14-MeV neutrons, the ITER port plug computational benchmark, and the Frascati Neutron Generator (FNG) ITER dose rate benchmark problem from the Shielding INtegral Benchmark Archive and Database (SINBAD). For each of these problems, comparisons were made to calculations using the rigorous two-step (R2S) method. The results on the iron prism problem illustrate how the D1S method achieves superior spatial resolution compared to the R2S method without the need for spatial discretization of the activation regions. The D1S and R2S results for the ITER port plug benchmark agree well with previously reported results in the literature. While the D1S results are 10% to 15% lower than the R2S results, this may be due to stochastic uncertainty and/or spatial discretization in the R2S calculations. On the FNG dose rate benchmark problem, the D1S method produces dose rate estimates that are within 4% of the dose rates predicted using a cell-based R2S workflow. The D1S estimates of the SDR are also in reasonable agreement with the experimental measurements and show the same basic trends that have been observed in previous works. A qualitative analysis of the execution time and uncertainty for the R2S and D1S workflows suggests that the D1S method would attain a higher figure of merit.

D1S method↗

DOE OSTI2021

Shake loss intensities in x-ray photoelectron spectroscopy: Theory, experiment, and atomic composition accuracy for MgO and related compounds

The relative intensities of XPS core levels, scaled by their photoionization cross sections, are regularly used to determine sample atomic composition. Cross sections, however, give the intensity to all possible final states for the core ionizations, not just to the main peak. This includes all intrinsic satellite structure (shake states and, for open shell systems, the different ionic multiplets). In practice, for solids, this is usually experimentally impossible to determine accurately because such a satellite structure sits on the inelastically scattered electron background and cannot be easily separated. Therefore, usually, only the intensity of the main peak is used. This limits the ultimate possible accuracy of XPS composition determination. The purpose of the present paper is to examine the contributions that a theoretical analysis of losses of intensity can make to improve quantitation. For an MgO single crystal, we show that the correct stoichiometry of 1:1 can be recovered using the theoretical analysis of the experimental MgO peak ratio intensities. For materials with a sufficient bandgap for the XPS main peaks to be separated from the scattered background, the intensity of main peaks can often be accurately determined. Thus, if one uses theory to calculate that fraction of the total intensity lost from a main peak into all its satellite structure, the intensity of just main peaks could then be used to more accurately determine relative atom % composition. This work tests this approach using a single crystal MgO (50% Mg, 50% O) standard. Ab initio electronic structure theory of representative MgO clusters is used to determine Hartree–Fock wave functions for the ground state and final ionized states corresponding to the main Mg 2p and O1s XPS peaks of the oxide. The sudden approximation, SA, is used to determine the fractional losses from these main peaks to shake satellites, which is found to be greater for O1s than Mg2p. This results in predicted “apparent composition” for stoichiometric MgO of 55.2% Mg, 44.8% O instead of the true 50% Mg, 50% O. Equivalent theory for CaO results in a predicted apparent Ca value of 53.4%. Experimentally, using Mg2s or 2p intensity ratio to O1s, we find values between 52.2% and 56.0% Mg using two crystals and four different instrument electron pass energies. The average value of the measurements is 54.5% Mg when corrected for the presence of an adventitious carbon overlayer and slight surface hydroxide. Though this agreement with theory may be somewhat fortuitous, given the potential experimental errors, which are fully discussed, it is similar to that in our earlier study on LiF. We also present preliminary experimental data on Mg(OH) 2 and MgSO 4 , which show a similar trend of apparently higher than 50% Mg, but we have no theory values. We are not yet able to experimentally test for validation of the difference between apparent composition for MgO (55.2% Mg) and CaO (53.4% Ca), owing to significant carbonate formation at the surface of the single crystal CaO. Here, an important conclusion is that the theoretical determination of shake losses, obtained with ab initio wavefunctions and the SA, is likely to be a useful way to calibrate the accuracy and reliability of compositions obtained from XPS intensities and merits further study.

47 OTHER INSTRUMENTATION↗