Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPECTATION VALUE”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

High-speed tunable generation of random number distributions using actuated perpendicular magnetic tunnel junctions

Perpendicular magnetic tunnel junctions (pMTJs) actuated by nanosecond pulses are emerging as promising devices for true random number generation (TRNG) due to their intrinsic stochastic behavior and high throughput. In this work, we demonstrate the tunability and quality of random number distributions generated by pMTJs operating at a frequency of 104 MHz. First, changing the pulse amplitude is used to systematically vary the probability bias. The variance of the resulting bitstreams closely matches the expected binomial distribution, demonstrating consistency with an underlying sequence of Bernoulli trials. Second, the quality of uniform distributions of 8-bit random numbers generated with a probability bias of 0.5 is considered. A reduced chi-square analysis of these data shows that only two XOR operations are sufficient to achieve this distribution with p-values greater than 0.05. Finally, we show that there is a correlation between long-term probability bias variations and pMTJ resistance. These findings suggest that variations in the characteristics of the pMTJ underlie the observed variation of probability bias. In conclusion, our results highlight the potential of stochastically actuated pMTJs for high-speed, tunable TRNG applications, showing the importance of the stability of pMTJ device characteristics in achieving reliable, long-term performance.

Magnetic tunnel junctions↗
DOE OSTI2024

Search for singly produced vectorlike top partners in multilepton final states with 139 fb −1 of 𝑝⁢𝑝 collision data at √𝑠 =13 TeV with the ATLAS detector

A search for the single production of a vectorlike top partner (𝑇) with mass greater than 1 TeV decaying into a 𝑍 boson and a top quark is presented, using the full Run 2 dataset corresponding to 139 fb −1 of 𝑝⁢𝑝 collisions at √𝑠 =13 TeV, collected in 2015–2018 with the ATLAS detector at the Large Hadron Collider. The targeted final state is characterized by the presence of a pair of electrons or muons with opposite-sign charges which form a 𝑍-boson candidate, as well as by the presence of 𝑏-tagged jets and forward jets. Events with exactly two or at least three leptons are categorized into two independently optimized analysis channels. No significant excess above the background expectation is observed and the results from the two channels are statistically combined to set exclusion limits at 95% confidence level on the masses and couplings of 𝑇. The results are interpreted in several benchmark scenarios to set limits on the mass and universal coupling strength (𝜅) of the vectorlike quark. For singlet 𝑇 quarks, 𝜅 values between 0.22 and 0.64 are excluded for masses between 1000 and 1975 GeV. For 𝑇 quarks in the doublet scenario, where the production cross section is much lower, 𝜅 values between 0.54 and 0.88 are excluded for masses between 1000 and 1425 GeV.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Collaborative Research: Properties and Dynamics of the Shallow Crust (Final Report)

Ground motions recorded at one location are often extrapolated to nearby regions within a given radius. Here, we explore the appropriateness of spatial extrapolation using data from seven small aperture seismic network deployments in southern California. Six of these deployments are linear arrays of 4-13 stations, and one is a 2D array of 13 stations at Pinyon Flats Observatory. The spatial footprint array diameters are 3 km or less, and each array was operational for a year or more. From our base catalog (M2.5+ earthquakes; 4038 events; September 2010 - June 2023), automated methods remove temporally overprinted waveforms from nearby events (< 5 km) in quick succession (< 5 min) and data with nonviable waveforms. These 200 samples per second data are filtered at 0.5-25 Hz and must have signal-tonoise ratios (SNRs) of 2.5+. Peak ground acceleration (PGA) and peak ground velocity (PGV) are derived individually from the maximum absolute values of each of the 3-component waveforms (vertical, northsouth, and east-west). Five of the seven arrays traverse the San Jacinto fault, and two do not. Ground motion observations are compared with theoretical estimates from Abrahamson et al., 2014. On average, arrays deployed within and across fault zones consistently record ground motions above theoretical expectations, whereas off-fault arrays record ground motions at or slightly below theoretical expectations. We attribute these differences to site conditions because these trends prevail for the full data suites. For each network and each individual channel, the coefficient of variation indicates that the standard deviations are ~30±6% of the mean. Exploring relative ground motion contributions from all three channels (ternary plots), as expected, most data show that vertical ground motions are attenuated compared to horizontal ones. However, this is not always the case for RA array data, where vertical motions can be ~2-3 times larger than horizontal motions for select events near Cahuilla, CA. These anomalously high vertical motions are focal mechanism-related. These results suggest that ternary plots created using only a small amount of data can be used as a data quality metric and a tool to find anomalous features in three-component data.

58 GEOSCIENCES↗
DOE OSTI2021

System Modeling of the HTTR and Economic Dispatch Model of the Secondary System

High Temperature Gas-cooled Reactors (HTGRs) can be used for the generation of electricity and their process heat can be used to improve the efficiency of chemical processes such as hydrogen production. The JAEA-operated High-Temperature engineering Test Reactor (HTTR-GT/H2) is exploring using the reactor for electricity and hydrogen production. A RELAP5-3D model of the HTTR-GT/H2 secondary system has been developed using design information. The various components and heat exchangers in the secondary system were modeled and results were compared to the design conditions. The results for the sole-power generation mode were shown to fit the design conditions very well. The largest temperature difference was on the order of 7 K, and the largest pressure difference was on the order of 0.05 MPa. The results for the hydrogen cogeneration mode did not match the design conditions nearly as well. The largest temperature difference was about 39 K and the largest pressure difference was about 0.27 MPa at the compressor outlet. The larger differences for the hydrogen cogeneration mode are attributed to the various complex components and the flow being split in the secondary loop. A transient reduction in heat removal capability of the secondary system was investigated. Reactor temperatures are anticipated to rise as a result. The core reactivity response due to this increase in temperature is investigated and is expected to add negative reactivity to the reactor. An economic dispatch model was developed for a nuclear-driven iodine-sulfur cycle system to determine hydrogen sale prices that would make such a system profitable. The study focuses on the development of the economic model and the role that input data plays on final calculated values. It was found that the input electricity prices, whether using historical data or a host of synthetic time histories, produce significantly different breakeven hydrogen sale prices. As such, great care should be used in these economic dispatch analyses to select reasonable input assumptions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Effect of Anoxic Iron Corrosion on WIPP Brine Geochemistry FY23 Final Report (U)

A 280-day study was completed to evaluate the effect of zero-valent iron (Fe 0 ) on the Waste Isolation Pilot Plant (WIPP) brine geochemistry under anticipated reducing conditions. Hydrogen (H 2 ) gas is expected to be present in the repository after closure due to the anoxic corrosion of a vast quantity of iron contained in the waste forms disposed at WIPP; therefore, a background argon atmosphere containing H 2 was chosen for this study. WIPP groundwater brine pH and E h will impact the mobility and fate of plutonium within the repository. Modeling and laboratory results for Castile WIPP brine indicate that equilibrium fa values relative to the standard hydrogen electrode (SHE) are 40 mV more reducing (i.e., more negative) than those for Salado WIPP brine (-480 mV vs. -440 mV, respectively) because of the higher pH of the Castile brine (pH 9 .3 for Castile vs. pH 8.8 for Salado). The E h and pH data were corrected for the effects of high ionic strength. The experimental results for both brines are consistent with thermodynamic predictions using OLI Systems' Mixed Solvent Electrolyte chemical equilibrium model. The measured and corrected pH and E h data from this study are provided in Table ES-I and Table ES-2, respectively. The experimental study, with four test conditions in triplicate, was performed in a dual glovebox with a nominally 3 vol.% H 2 in argon atmosphere (target H 2 range: 3 ± I vol.%). Simulants containing MgO only ( experimental control) and MgO+Fe 0 (WIPP base case) were prepared for both the Salado and Castile brines. MgO was included in all simulants to account for the use of bulk magnesium oxide in the WIPP repository. Fe 0 was included in some simulants to incorporate the effects of the anoxic corrosion of iron and in-situ hydrogen generation in the study. The brine compositions were developed by Sandia National Laboratory (SNL; Xiong, 2008) and have been used in previous WIPP evaluations. The test method (agitation, etc.) is partially based on ASTM D3987-12. Twelve rounds of periodic measurements of pH and E h were performed over the course of the study. Chemical analysis results for liquids and solids (ICP-MS, ICP-ES, IC Anion, TIC, SEM-EDX) are consistent with the pH, E h , and thermodynamic modeling results. This study included the following conditions that deviate from anticipated post-closure conditions following brine intrusion, but were selected to facilitate bench-scale testing to validate modeling of pH and E h for the post-closure WIP P repository: an anoxic glove box atmosphere containing ≤ 4 vol. % H 2 vs. substantially higher H 2 gas concentrations assumed in the WIPP Performance Assessment (PA); a significantly higher liquid-to-solid test ratio compared to the much lower phase ratio anticipated in the WIP P repository; agitation of the simulant bottles to maximize mass transfer; and finally the use of Fe 0 reagents having a much greater surface area than expected in the WIP P repository. Non-representative conditions were chosen for various reasons such as: to provide bounding conservative results, to provide a margin of safety for testing, or to facilitate simulant sub-sampling and analysis. In a parallel effort, aqueous electrolyte thermodynamic models were developed for the synthetic Salado and Castile brines to inform the experimental design, facilitate laboratory data interpretation, and allow extension of evaluations beyond the parameters tested. Thermodynamic modeling simulations including the MgO and Fe 0 additives that are directly relevant to the experimental measurements (e.g., pH calibration curve, ORP corrections) are included in this report. The measured fa of the simulants was close to the OLI model predictions for both brines and was largely controlled by the background H 2 partial pressure in the vapor phase as well as H 2 generated in situ in the aqueous phase by the Fe 0 corrosion. The H 2 gas-phase concentration tested and thermodynamically evaluated was much lower than is assumed in the WIPP PA; however, H 2 (g) concentrations significantly below this level are still predicted to result in very reducing conditions. In conclusion: • The experimental results are consistent with thermodynamic model predictions for fa, pH, and the effects of high ionic strength. • Evidence to date suggests that the H2 concentration in the glovebox atmosphere ultimately determined the final E h values of the simulants and resulted in highly reducing conditions. As a result, little difference was observed between the control simulants containing only MgO and the WIPP base-case simulants that contained MgO and Fe 0 . • This test methodology is recommended for future studies evaluating WIPP repository conditions. The methodology includes: (1) background H 2 in argon with agitation ( or could alternatively include in-situ-generated H 2 in sealed bottles); (2) carefully measured and corrected ORP data ( with much effort focused on allowing the probes to fully stabilize); and (3) ionic-strength-corrected pH data. Other best practices, such as simulant sparging/handling, ORP probe replacement, etc., should also be considered. • The coupling of experimental studies and thermodynamic modeling is also highly recommended because these methods inform and direct one another leading to greater confidence in and understanding of the results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Much of zero emissions commitment occurs before reaching net zero emissions

Abstract We explore the response of the Earth’s coupled climate and carbon system to an idealized sequential addition and removal of CO 2 to the atmosphere, following a symmetric and continuous emissions pathway, in contrast to the discontinuous emissions pathways that have largely informed our understanding of the climate response to net zero and net negative emissions to date. We find, using both an Earth system model and an ensemble of simple climate model realizations, that warming during the emissions reduction and negative emissions phases is defined by a combination of a proportionality of warming to cumulative emissions characterized by the transient climate response to emissions (TCRE), and a deviation from that proportionality that is governed by the zero emissions commitment (ZEC). About half of the ZEC is realized before reaching zero emissions, and the ZEC thus also controls the timing between peak cumulative CO 2 emissions and peak temperature, such that peak temperature may occur before peak cumulative emissions if ZEC is negative, underscoring the importance of ZEC in climate policies aimed to limit peak warming. Thus we argue that ZEC is better defined as the committed warming relative to the expected TCRE proportionality, rather than as the additional committed warming that will occur after reaching net zero CO 2 emissions. Once established, the combined TCRE and ZEC relationship holds almost to complete removal of prior cumulative CO 2 emissions. As cumulative CO 2 emissions approach zero through negative CO 2 emissions, CO 2 concentrations drop below preindustrial values, while residual long-term climate change continues, governed by multicentennial dynamical processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Shape, colour, and distance in weak gravitational flexion

ABSTRACT Canonically, elliptical galaxies might be expected to have a perfect rotational symmetry, making them ideal targets for flexion studies – however, this assumption has not been tested. We have undertaken an analysis of low- and high-redshift galaxy catalogues of known morphological type with a new gravitational lensing code, Lenser. Using colour measurements in the u − r bands and fit Sérsic index values, objects with characteristics consistent with early-type galaxies are found to have a lower intrinsic scatter in flexion signal than late-type galaxies. We find this measured flexion noise can be reduced by more than a factor of two at both low and high redshift.

Fabritius II, Joseph M.↗
DOE OSTI2024

Measurement of the W-boson mass and width with the ATLAS detector using proton–proton collisions at $\sqrt{s}=7$ TeV

Proton-proton collision data recorded by the ATLAS detector in 2011, at a centre-of-mass energy of 7 TeV, have been used for an improved determination of the W-boson mass and a first measurement of the W-boson width at the LHC. Recent fits to the proton parton distribution functions are incorporated in the measurement procedure and an improved statistical method is used to increase the measurement precision. The measurement of the W-boson mass yields a value of m w = 80,366.5 ± 9.8 (stat.) ± 12.5 (syst.) MeV = 80,366.5 ± 15.9 MeV, and the width is measured as Γ w = 2202 ± 32 (stat.) ± 34 (syst.) MeV = 2202 ± 47 MeV. The first uncertainty components are statistical and the second correspond to the experimental and physics-modelling systematic uncertainties. Both results are consistent with the expectation from fits to electroweak precision data. The present measurement of m w is compatible with and supersedes the previous measurement performed using the same data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A Technical and Economic Assessment of LWR Flexible Operation for Generation and Demand Balancing to Optimize Plant Revenue

With increased penetration of subsidized variable renewable energy (VRE) resources and competition from low natural gas prices, existing light water reactor (LWR) nuclear power plants (NPPs) are struggling to remain economically competitive. This work examines the potential economic competitiveness of various thermal energy storage (TES) technologies when coupled directly or indirectly with a NPP. To highlight their relative economic competitiveness, we contrast several energy storage solutions in stochastic dispatch optimization. We leverage data from recent work analyzing a range of TES technologies with varying capital costs, performance, and technology readiness level (TRL) to establish our case. We explore inserting these technologies into an electricity market with existing nuclear generation and large projected variable renewable energy (VRE) penetration. Although these technologies' projected capital costs may make them unlikely candidates in their current state, this analysis demonstrates a high-fidelity techno-economic analysis of energy storage. Furthermore, as the projected cost of energy storage technologies evolves, this analysis sets a precedent for similar future investigations. One region with projected trends that may be unfavorable for existing nuclear capacity is the New York Independent System Operator (NYISO) market. New York state’s baseload generation has been historically provided by fossil-fired and nuclear assets. However, amid economic pressures from subsidized VREs and low natural gas prices, the state has recently deactivated Indian Point nuclear power plant units 2 and 3. Furthermore, the state plans to meet its zero-emission generation target by 2040 by replacing fossil-fired capacity with significant investments in VRE resources like wind and solar photovoltaic (PV) and battery storage. Increased intermittent resource penetration lowers the baseload power requirement, adding further economic pressure to the state’s three remaining NPPs still in operation. With three NPPs still in operation in New York, this work analyzes potential economic benefits to NPPs on the New York grid when directly or indirectly coupled with various TES technologies. This work requires two modeling steps to analyze the potential economic benefits of various system configurations of the TES directly or indirectly coupled with nuclear. First, this analysis leverages capacity expansion modeling by experts at the Electric Power Research Institute (EPRI). Using their deterministic capacity expansion model, U.S. Regional Economy, Greenhouse Gas, and Energy (US-REGEN), EPRI analysts evaluated the capacity and generation evolution of the New York state energy market under four projection scenarios. These four projection scenarios were developed to represent the potential evolution of the capacity and generation in NYISO from 2015 to 2050 under various economic, technology, and policy constraints. The results from these capacity expansion models are then used as boundary conditions in the second modeling step. The second modeling step uses the Holistic Energy Resource Optimization Network (HERON) for a set of stochastic techno-economic analyses (STEAs) to investigate the potential increase in the economic viability of various configurations of the TES. With no current capacity expansion capabilities, HERON takes the data generated from US-REGEN for 2050 to generate synthetic load, solar, and wind data. Then HERON economically optimizes the capacity and dispatch of the various TES configurations. The potential economic benefit is the differential net present value (NPV) of the TES configurations from the no-TES baseline. As a stochastic techno-economic analysis package, HERON introduces uncertainty into the economic metrics, while US-REGEN trades resolution for reduced computational complexity. Using HERON also allows the modeling of direct thermal coupling, a feature not common in capacity and dispatch models. As expected, with high capital costs, the costs of introducing energy storage for all the technologies considered outweighed the potential economic benefit of this strategy for flexible plant operation. The benefit of this analysis is primarily in demonstrating a workflow that examines innovative solutions to increase NPP revenue via TES coupling. HERON’s stochastic capacity and dispatch optimization process used in this work has proven an effective tool in observing and evaluating the impact of introducing storage technologies in a grid energy system.

25 ENERGY STORAGE↗
DOE OSTI2021

Command and Control of Utilities Systems

The utilities systems that encompass the Y-12 National Security Complex are an integral part of the entire campus, and are essential to maintaining critical mission focus for the site. Although the utilities systems each differ in function, the central production and distribution model creates a common thread on which each system is run. Several incidents in FY 2018 and FY 2019 documented a lack of overall systems control and highlighted the need to create a better command system that could maintain system availability and decrease possible impacts to the surrounding environment. An initial system was started in September of 2019 to better control mechanical utility distribution systems. The measurement of the system can be modeled against a standard command and control related incident. The corresponding work hours from a typical control-related incident were used as key component to this system model. System effectiveness is modeled against utility’s incidents that demonstrate a loss of command and control. The increased simulations of the model also document control can only be effectively contained within the stated time (3-5 years) by yielding more engineers to complete alignments, or drastically reducing the amount of time for alignments. The total control points of the model document a clear correlation between alignments, control points, and the amount of man-hours per incident. Current data shows that actual incident man-hours (1.88 manhours/ incident since initial system implementation) are near expected incident man-hours (1.81 manhours/ incident) with command and control for October 2021, which documents a 33.3% reduction in man-hours/incident since initial system implementation. A conclusion can be drawn that focus on increasing the amount of controlled points in the system can reduce incidents and the decrease the amount of non-value added work to employee’s day-to-day lives, while also achieving a culture that enhances and sustains utilities reliability.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Modernizing the Nuclear Industry and New Ways of Working

Nuclear energy is recognized as the most feasible energy source towards achieving nation-wide net-zero goals (COP 208, MIT). Additionally, the demand for reliable, consistent energy supply is soaring with the construction of energy hungry data centers and AI tools. This means the construction of new nuclear and the continuation of current nuclear are a high national priority. However, this leaves the nuclear industry with a workforce challenge. The rising demand for nuclear power is threatened by an underpopulated workforce. There exist a few contributors to the decreased workforce such as skewed workforce demographics resulting in large-scale workforce retirement. Also, some plants have observed difficulty hiring and retaining skill sets that typically comprise the nuclear power workforce – specifically skilled craftsman. Those who are graduating with desirable skills are seeking employment in sectors that are more modern and culturally more aligned with younger generational values. Lastly, within the nuclear sector, skilled employees will be a competitive commodity as advanced plants come online offering work environments that incorporate modern technologies, skill sets, and opportunities for career advancement. These factors emphasize the need for legacy plants to modernize with a focus on the workplace environment and culture that meets younger generations’ skills, cultural expectations, and desire for advancement opportunities. Nuclear plants are actively engaged in deploying advancements that improve the management of systems, structures, and components as well as process improvements and technologies that can reduce the costs of operation and maintenance. However, the advancements must also be analyzed, reviewed and deployed in a manner that considers the impact on how people within the organization collaborate, communicate, make decisions, and solve problems. Neglecting to consider the cultural impact of modernization risks poor adoption, unrealized opportunities, or insignificant change toward helping attract new employees. The traditional way of working is not necessarily the way new generations want to engage with their employer. For instance, in a survey performed by North American Young Generation in Nuclear the top three reasons for younger nuclear employees to seek other job opportunities was seeking better work-life balance, a lack of advancement opportunities, and work culture and leadership style differences (Smyth et al. 2022). Integrated Operations for Nuclear (ION) is an approach for the nuclear industry to create long-term strategic modernization plans and analyze each advancement for the impact on people and processes and measure how those impacts flow up to support high-level, long-term plant goals and requirements. One principle of ION is to replace the current labor-centric operation style in legacy nuclear plants with data-centric collaborative operation styles. This paper will explain how transitioning operations to leverage centralized skills and responsibility, multi-skilled teams, and collaborative decision making can flatten the hierarchical structure in plants, enable greater individual efficacy while also offering more experiences and advancement opportunities to staff. Adopting the ION way of working shifts the industry mind-set towards more networked, collaborative, and modern way of working that will attract skilled, next-generation workers to operate the legacy nuclear fleet.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

Central vortex steady states and dynamics of Bose–Einstein condensates interacting with a microwave field

Here we study central vortex steady states and dynamics of a two-dimensional (2D), two-component, Gross–Pitaevskii equation (CGPE) system for two pseudo-spinor Bose Einstein condensates (BECs) interacting with an electromagnetic field (microwave) analytically and numerically. For the central vortex steady state at any given winding number S, we prove its existence in a reduced, single component detuning limit when contact-interaction strength β < β b , and nonexistence when β > β b , respectively, where β b is a threshold value, whose value is given in Theorem 3.1 in the paper. We extend the existence and nonexistence result to the general two pseudo-spinor case and prove that a central vortex steady state exists for any given S if β ≤ β b while it does not exist when β > 2β b . We then derive dynamical equations for some observables (expectations of the matter wave function) such as the center-of-mass, position of dispersion, the linear and angular momentum of the two pseudo-spinor CGPEs. Finally, numerical computations are brought in to validate and extend the existence of vortex steady state results to β b ≤ β < 2β b for the two pseudo-spinor case and to explore transient dynamics of the observables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Vessel Rupture Thresholds for Vessel–Bubble Interactions Using an Earthworm Vasculature Model

Intravenous microbubble oscillation in the presence of ultrasound has the potential to yield a wide range of therapeutic benefits. However, the likelihood of vessel damage caused by mechanical effects has not been quantified as a function of the numerous important parameters in therapeutic ultrasound procedures. In this study, we examined the effects of microbubbles injected into the vasculature of the earthworm. It was found that the elastic properties of earthworm blood vessels are similar to those of arteries in older humans, and that earthworms are well suited to the large number of experiments necessary to investigate safety of procedures involving microbubble oscillation in sonicated vessels. Microbubbles were infused into earthworm vessels, and the rupture time during sonication was recorded as a function of ultrasound frequency, pulse repetition frequency and acoustic pressure. A modified mechanical index (MMI) was defined that successfully captured the trends in rupture probability and rupture time for the different parameter values, creating a database of vessel rupture thresholds. In the absence of bubbles, the product of MMI squared and rupture time was approximately constant, indicating a possible radiation-force effect. Here, the MMI was an effective correlating parameter in the presence of bubbles, though the mathematical dependence is not yet apparent. The results of the study are expected to be valuable in designing more refined studies in vertebrate models, as well as informing computational models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Quantal diffusion approach for multinucleon transfer processes in the 58,64 Ni + 208 Pb reactions: Toward the production of unknown neutron-rich nuclei

In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Design and Analysis of Multifidelity Finite Element Simulations

Abstract The numerical accuracy of finite element analysis (FEA) depends on the number of finite elements used in the discretization of the space, which can be varied using the mesh size. The larger the number of elements, the more accurate the results are. However, the computational cost increases with the number of elements. In current practice, the experimenter chooses a mesh size that is expected to produce a reasonably accurate result, and for which the computer simulation can be completed in a reasonable amount of time. Improvements to this approach have been proposed using multifidelity modeling by choosing two or three mesh sizes. However, mesh size is a continuous parameter, and therefore, multifidelity simulations can be performed easily by choosing a different value for the mesh size for each of the simulations. In this article, we develop a method to optimally find the mesh sizes for each simulation and satisfy the same time constraints as a single or a double mesh size experiment. A range of different mesh sizes used in the proposed method allows one to fit multifidelity models more reliably and predict the outcome when meshes approach infinitesimally small, which is impossible to achieve in actual simulations. We illustrate our approach using an analytical function and a cantilever beam finite element analysis experiment.

Engineering↗
DOE OSTI2024

The PARADIGM Project: Case Study in Balancing Experiment Uncertainty with Design simplicity

Accurate nuclear data are required for simulations of many applications including nuclear criticality safety. Actinide nuclear data at intermediate energies (from 1 to 100s of keV) are imprecise and inaccurate, because of scarce differential data, and an insufficient theory approach to capture the structures expected in the data to yield evaluated nuclear data, and lack of integral data for proper validation. This is a known deficiency but has proved challenging to address. More specifically, only 5% of integral experiments in the International Criticality Safety Benchmark Evaluation Project (ICSBEP) benchmark suite address intermediate energies (Fig. 1). Associated calculated effective multiplication factor, k eff , values for these experiments are far outside the experimental uncertainties and are 25× further from experiment than for fast energies. These differences could either stem from systematic biases in nuclear data, experiments or both. The goal of the PARADIGM (PARallel Approach of Differential and InteGral Measurements) project is to significantly reduce (by more than tens of percent) the uncertainties of intermediate energy actinide nuclear data. The PARADIGM project designed and intends to execute LANSCE (Los Alamos Neutron Science CEnter) and NCERC (National Criticality Experiments Research Center) intermediate experiments in parallel. They will specifically address a high priority nuclear data need—reducing bias and uncertainty in intermediate plutonium nuclear data. The two experiment will achieve that by informing each other and nuclear theory. By doing all these steps in parallel, the timeline to deliver improved nuclear data to users will significantly be reduced. This work will focus on the integral experiment final design and the balance of design and modeling simplicity while minimizing experiment uncertainty.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Role of poloidal E × B drift in divertor heat transport in DIII‐D

Abstract Simulations for DIII‐D high confinement mode plasmas with the multifluid code UEDGE show a strong role of poloidal E × B drifts on divertor heat transport, challenging the paradigm of conduction‐limited scrape‐off layer (SOL) transport. While simulations with reduced drift magnitude are well aligned with the assumption that electron heat conduction dominates the SOL heat transport, simulations with drifts predict that the poloidal convective E × B heat transport dominates over electron heat conduction in both attached and detached conditions. As poloidal E × B flow propagates across magnetic field lines, poloidal transport with shallow magnetic pitch angles can reach values that are of the same order as would be provided by sonic flows parallel to the field lines. These flows can lead to strong convection‐dominated divertor heat transport, increasing the poloidal volume of radiative power front, consistent with previous measurements at DIII‐D. Due to these convective flows, the Lengyel integral approach, assuming zero convective fraction, is expected to provide a pessimistic estimate for the radiative capability of impurities in the divertor. For the DIII‐D simulations shown here, the Lengyel integral approach underestimates the radiated power by a factor of 6, indicating that, for reliable DIII‐D divertor power exhaust predictions, full two‐dimensional (2D) calculations, including drifts, would be necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

On the trapped magnetic moment in type-II superconductors

In this article, measurements of the trapped (remanent) magnetic moment, M trap (H) , when a small magnetic field H is turned off after cooling below the superconducting transition temperature, T c , or ramping a magnetic field up and down after cooling in a zero magnetic field, offer substantial advantages in difficult cases of small samples and large field-dependent backgrounds, which is relevant for hydrogen-based ultra-high-T c superconductors (UHTS). Until recently, there was no need for a separate paper on the trapped magnetic flux for well-known critical-state models due to the simplicity of the physics involved. However, recent publications showed the need for such an analysis. This note summarizes the expectations for the Bean model with constant critical current density and the Kim model with field-dependent critical currents. It is shown that if the trapped moment is fitted to the power law, M trap $\propto$ H a , the fixed exponent α = 2 is exact for the Bean model, while Kim models show a wide interval of possible values, 2$\leqslant$a$\leqslant$4. Furthermore, accounting for reversible magnetization expands the range of possible exponents to 1$\leqslant$a$\leqslant$4. In addition, demagnetizing factors are essential and make the trapped moment orientation dependent even in isotropic materials. As a concrete application, it is shown that flux trapping experiments performed on H 3 S UHTS compounds can be well described using this generalized approach, lending further support to the type-II superconducting nature of H 3 S under ultra-high pressure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗