Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “core simulation”

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

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

At least 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A new multi-energy neutrino radiation-hydrodynamics code in full general relativity and its application to the gravitational collapse of massive stars

We present a new multi-dimensional radiation-hydrodynamics code for massive stellar core-collapse in full general relativity (GR). Employing an M1 analytical closure scheme, we solve spectral neutrino transport of the radiation energy and momentum based on a truncated moment formalism. Regarding neutrino opacities, we take into account a baseline set in state-of-the-art simulations, in which inelastic neutrino–electron scattering, thermal neutrino production via pair annihilation, and nucleon–nucleon bremsstrahlung are included. While the Einstein field equations and the spatial advection terms in the radiation-hydrodynamics equations are evolved explicitly, the source terms due to neutrino–matter interactions and energy shift in the radiation moment equations are integrated implicitly by an iteration method. To verify our code, we first perform a series of standard radiation tests with analytical solutions that include the check of gravitational redshift and Doppler shift. A good agreement in these tests supports the reliability of the GR multi-energy neutrino transport scheme. We then conduct several test simulations of core-collapse, bounce, and shock stall of a 15M{sub ⊙} star in the Cartesian coordinates and make a detailed comparison with published results. Our code performs quite well to reproduce the results of full Boltzmann neutrino transport especially before bounce. In the postbounce phase, our code basically performs well, however, there are several differences that are most likely to come from the insufficient spatial resolution in our current 3D-GR models. For clarifying the resolution dependence and extending the code comparison in the late postbounce phase, we discuss that next-generation Exaflops class supercomputers are needed at least.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Second release of the CoRe database of binary neutron star merger waveforms

Abstract We present the second data release of gravitational waveforms from binary neutron star (BNS) merger simulations performed by the Computational Relativity ( CoRe ) collaboration. The current database consists of 254 different BNS configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to ℓ = m = 4 . They span a significant portion of the mass, mass-ratio, spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Assessment of the Analysis Capability for Core-wide PWR Pellet-Clad Interaction (PCI) Screening

The Level 1 Milestone, L1:CASL.P13.03, “Assess the analysis capability for core-wide PWR Pellet-Clad Interaction (PCI) screening and demonstrate detailed 3-D analysis on selected sub-region.”, has been completed by the Fuels, Materials and Chemistry and Physics Integration focus areas of CASL. This report illustrates the application of VERA simulation to PWR core-wide PCI screening, with emphasis on the multi-dimensional analysis capabilities of BISON fuel performance code. The VERA codes utilized in this effort are: the fuel performance code BISON and the VERA-CS sub-element, which combines the neutron transport code MPACT and the thermal-hydraulics code COBRA-TF (CTF). The target reactor for this work is Watts Bar Unit 1, Cycles 6 and 7, selected on the basis that PCI failures had occurred during power start-up in both cycles, despite the moderate power levels, which offered an interesting test of the code’s veracity and a particularly challenging condition for this Milestone’s assessment. The analysis began by performing BISON-IN-VERA, i.e. VERA-CS + BISON, quarter-core depletion calculations, from which limiting assemblies and limiting fuel rods were selected using a four-parameter selection criterion: high power, high temperature, small gap thickness and high cladding hoop stress. The cladding material being Westinghouse Zirlo, the analysis utilized a new material model developed for this cladding by CASL’s FMC team using open literature sources. The VERA-CS + BISON analysis resulted in identifying 51 limiting rods as candidates for more detailed PCI analysis using stand-alone BISON, employing local 2D and 3D geometric models. This local-effects BISON analysis provided detailed PCI failure assessment, considering the impact of power histories on the two potential sources of PCI failure: the classical SCC-PCI and the missing pellet surface geometry MPS-PCI. The 51-rod analyses of Watts Bar Unit 1 cycle 6 and 7 revealed that the classical SCC-PCI stress levels were well below the failure threshold determined in a previous CASL study, which implicated the other mode of failure, namely MPS-PCI. Considering MPS defects to be the root cause for the failure, the local-effects calculations was repeated using MPS defect sizes equal to and larger than Westinghouse quality assurance standards. Initial results of this new analysis showed a significant increase in the cladding stress, indicating that an MPS may have been present in cycle 6 and 7 for the failure to occur. In order to determine the size of MPS needed for cycle 6 and 7 rods to fail, a sensitivity study using 2-D and 3-D models was performed. The results showed cladding hoop stresses, while remaining below the failure level for an MPS of equal size to the acceptance standard, continued to increase as the size of the MPS increased. This report has demonstrated the use of BISON, in consort with VERA’s simulation methodology, for core-wide PCI failure screening capability. Utilizing this demonstration, and emulating the PCI margin analysis methodology developed in the EPRI PCI guidelines, BISON-IN-VERA simulation can be used very effectively to develop core power histories to minimize PCI failure potential.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Multiscale simulations of electron and ion dynamics in self-irradiated silicon

The interaction of energetic ions with the electronic and ionic system of target materials is an interesting but challenging multiscale problem, and understanding of the early stages after impact of heavy, initially charged ions is particularly poor. At the same time, energy deposition during these early stages determines later formation of damage cascades. In this work, we address the multiscale character by combining real-time time-dependent density functional theory for electron dynamics with molecular dynamics simulations of damage cascades. Our first-principles simulations prove that core electrons affect electronic stopping and have an unexpected influence on the charge state of the projectile. We show that this effect is absent for light projectiles, but dominates the stopping physics for heavy projectiles. By parametrizing an inelastic energy loss friction term in the molecular dynamics simulations using our first-principles results, we also show a qualitative influence of electronic stopping physics on radiation-damage cascades.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2019

Building Booster Separation Aerodynamic Databases for Artemis II

NASA’s Artemis II mission will mark the return of humans to near-lunar space for the first time since Apollo. Shortly after launch on the Space Launch System (SLS), a critical phase of ascent occurs when 16 small rockets fire to push the boosters away from the core. Minimizing the risk of failure during separation requires the construction of multiple 13-dimensional databases, including perturbations in position, flight conditions, and engine thrust. The SLS Computational Fluid Dynamics team used NASA’s FUN3D flow solver on the Pleiades and Electra supercomputers to run 5,780 simulations at nominal conditions and over 8,000 simulations with a core stage engine failure to generate the databases needed to verify the booster separation system for Artemis II.

Meeroff, Jamie G.↗
DOE OSTI2022

Computer Science Research Needs for Parallel Discrete Event Simulation (PDES)

Historically, scientific computing efforts have demonstrated the clear need for, and effective use of, supercomputing with traditional time-stepped simulations. Nevertheless, there are several areas in the mission spaces of the U.S. Department of Energy and other agencies waiting to tap advanced computing research using a different, discrete event style of modeling, simulation, and analysis. These span a wide spectrum of applications including energy grid resilience, urban planning and policy, transportation science, building technologies, emergency response and planning, environmental impact analysis, computational epidemiology, Internet communications, cyber security, and cyber-physical systems, to name only a few. Even within traditional scientific applications, the role of discrete event modes of execution is increasing in the form of new event-based mathematical solvers such as quantized state integration methods and discrete-continuous hybrid system solvers. Co-design of advanced supercomputing hardware systems is another area that exploits discrete event simulation at its core for effective analyses. Complex systems, entity behaviors and interconnections play a significant role in all these applications, which are mapped to large-scale models with discrete event formulations. To make advancements in all the aforementioned scientific areas, many technical aspects need to be more thoroughly studied and deeply understood in parallel discrete event simulation (PDES). The unique dynamics inherent in a discrete event modeling approach, by their very nature, intersect and influence the entire stack of the computing system, including (a) the unique nature of the instruction sets exercised in PDES workloads without a predominance of high-precision floating point operations, (b) virtual time-constrained multi-threaded execution of many logical processes per processor, (c) extremely variable and difficult to predict network traffic characteristics, (d) interfaces and inter-dependencies with machine learning and artificial intelligence codes at higher software layers, and (e) highly challenging load balancing needs, especially in effectively accounting for accelerated/extremely heterogeneous computing in current and future high-performance computing systems. Efficient and accurate parallel execution of PDES workloads is also dominated by challenges in dealing with their asynchronous concurrency fundamentally present at the model level. Conservative synchronization, optimistic/speculative synchronization, and their hybrid schemes open new questions in fundamental computer science with respect to reversibility of computation and prediction (lookahead) of behaviors inherent within model codes. On the implementation front, there are relatively few scalable, general-purpose parallel discrete event simulators in the world, and even fewer have been studied on emerging hardware platforms. To enable scientific advances using PDES, the research needs in computer science must also be pursued and met in the intersection of the algorithmic and hardware-aware aspects of scalable PDES engines. This report is aimed at capturing a computer science-oriented view of this important area of research in PDES, presenting a sample of important applications with their inherent discrete event technology elements. Needs are outlined in core areas of parallel discrete event research as well as cross-cutting directions in computer science research that positively impact scientific advancements across several important application areas. A selection of priority research opportunities in advanced computing for PDES is identified to serve as reference for key research topics and their order of importance for scientific advancements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

GenASiS: General Astrophysical Simulation System. II. Self-gravitating Baryonic Matter*

GenASiS (General Astrophysical Simulation System) is a code being developed initially and primarily, though not exclusively, for the simulation of core-collapse supernovae on the world's leading capability supercomputers. This paper---the second in a series---documents capabilities for Newtonian self-gravitating fluid dynamics, including tabulated microphysical equations of state treating nuclei and nuclear matter (`baryonic matter'). Computation of the gravitational potential of a spheroid, and simulation of the gravitational collapse of dust and of an ideal fluid, provide tests of self-gravitation against known solutions. In multidimensional computations of the adiabatic collapse, bounce, and explosion of spherically symmetric pre-supernova progenitors---which we propose become a standard benchmark for code comparisons---we find that the explosions are prompt and remain spherically symmetric (as expected), with an average shock expansion speed and total kinetic energy that are inversely correlated with the progenitor mass at the onset of collapse and the compactness parameter.

Cardall, Christian [ORNL] (ORCID:000000020086105X)↗
DOE OSTI2023

Three approaches for the classification of protoneutron star oscillation modes

ABSTRACT The future detection of gravitational waves (GWs) from a Galactic core-collapse supernova will provide information on the physics inside protoneutron stars (PNS). In this work, we apply three different classification methods for the PNS non-radial oscillation modes: Cowling classification, Generalized Cowling Nomenclature (GCN), and a classification based on modal properties (CBMP). Using PNS models from 3D simulations of core-collapse supernovae, we find that in the early stages of the PNS evolution, typically 0.4 s after the bounce, the Cowling classification is inconsistent, but the GCN and the CBMP provide complementary information that helps to understand the evolution of the modes. In the GCN, we note several avoided crossings as the mode frequencies evolve at early times, while the CBMP tracks the modes across the avoided crossings. We verify that the strongest emission of GWs by the PNS corresponds to the f mode in the GCN, indicating that the mode trapping region alternates between the core and the envelope at each avoided crossing. At later times, approximately 0.4 s after the bounce, the three classification methods present a similar description of the mode spectrum. We use our results to test universal relations for the PNS modes according to their classification and find that the behaviour of the universal relations for f and p modes is remarkably simple in the CBMP.

Astronomy & Astrophysics↗
NASA NTRS2006

Comprehensive Report of Fan Performance From Duct Rake Instrumentation on 1.294 Pressure Ratio, 806 ft/sec Tip Speed Turbofan Simulator Models

A large scale model representative of an advanced ducted propulsor-type, low-noise, very high bypass ratio turbofan engine was tested for acoustics, aerodynamic performance, and off-design operability in the NASA Glenn 9- by 15-Foot Low-Speed Wind Tunnel. The test was part of NASA s Advanced Subsonic Technology Noise Reduction Program. The low tip speed fan, nacelle, and un-powered core passage were simulated. As might be expected, the effect of stall management casing treatment was a performance penalty. Reducing the recirculating flow at the fan tip reduced the penalty while still providing sufficient stall margin. Two fans were tested with the same aerodynamic design; one with graphite composite material, and the other with solid titanium. There were surprising performance differences between the two fans, though both blades showed some indication of transitional flow near the tips. Though the pressure and temperature ratios were low for this fan design, the techniques used to improve thermocouple measurement accuracy gave repeatable data with adiabatic efficiencies agreeing within 1 percent. The measured fan adiabatic efficiency at simulated takeoff conditions was 93.7 percent and matched the design intent.

Jeracki, Robert J.↗
DOE OSTI2014

MO-E-18C-04: Advanced Computer Simulation and Visualization Tools for Enhanced Understanding of Core Medical Physics Concepts

Purpose: Most medical physics programs emphasize proficiency in routine clinical calculations and QA. The formulaic aspect of these calculations and prescriptive nature of measurement protocols obviate the need to frequently apply basic physical principles, which, therefore, gradually decay away from memory. E.g. few students appreciate the role of electron transport in photon dose, making it difficult to understand key concepts such as dose buildup, electronic disequilibrium effects and Bragg-Gray theory. These conceptual deficiencies manifest when the physicist encounters a new system, requiring knowledge beyond routine activities. Methods: Two interactive computer simulation tools are developed to facilitate deeper learning of physical principles. One is a Monte Carlo code written with a strong educational aspect. The code can “label” regions and interactions to highlight specific aspects of the physics, e.g., certain regions can be designated as “starters” or “crossers,” and any interaction type can be turned on and off. Full 3D tracks with specific portions highlighted further enhance the visualization of radiation transport problems. The second code calculates and displays trajectories of a collection electrons under arbitrary space/time dependent Lorentz force using relativistic kinematics. Results: Using the Monte Carlo code, the student can interactively study photon and electron transport through visualization of dose components, particle tracks, and interaction types. The code can, for instance, be used to study kerma-dose relationship, explore electronic disequilibrium near interfaces, or visualize kernels by using interaction forcing. The electromagnetic simulator enables the student to explore accelerating mechanisms and particle optics in devices such as cyclotrons and linacs. Conclusion: The proposed tools are designed to enhance understanding of abstract concepts by highlighting various aspects of the physics. The simulations serve as virtual experiments that give deeper and long lasting understanding of core principles. The student can then make sound judgements in novel situations encountered beyond routine clinical activities.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1995

Microwave and infrared simulations of an intense convective system and comparison with aircraft observations

A three-dimensional cloud model, radiative transfer model-based simulation system is tested and validated against the aircraft-based radiance observations of an intense convective system in southeastern Virginia on 29 June 1986 during the Cooperative Huntsville Meteorological Experiment. NASA's ER-2, a high-altitude research aircraft with a complement of radiometers operating at 11-micrometer infrared channel and 18-, 37-, 92-, and 183-GHz microwave channels provided data for this study. The cloud model successfully simulated the cloud system with regard to aircraft- and radar-observed cloud-top heights and diameters and with regard to radar-observed reflectivity structure. For the simulation time found to correspond best with the aircraft- and radar-observed structure, brightness temperatures T(sub b) are simulated and compared with observations for all the microwave frequencies along with the 11-micrometer infrared channel. Radiance calculations at the various frequencies correspond well with the aircraft observations in the areas of deep convection. The clustering of 37-147-GHz T(sub b) observations and the isolation of the 18-GHz values over the convective cores are well simulated by the model. The radiative transfer model, in general, is able to simulate the observations reasonably well from 18 GHz through 174 GHz within all convective areas of the cloud system. When the aircraft-observed 18- and 37-GHz, and 90- and 174-GHz T(sub b) are plotted against each other, the relationships have a gradual difference in the slope due to the differences in the ice particle size in the convective and more stratiform areas of the cloud. The model is able to capture these differences observed by the aircraft. Brightness temperature-rain rate relationships compare reasonably well with the aircraft observations in terms of the slope of the relationship. The model calculations are also extended to select high-frequency channels at 220, 340, and 400 GHz to simulate the Millimeter-wave Imaging Radiometer aircraft instrument to be flown in the near future. All three of these frequencies are able to discriminate the convective and anvil portions of the system, providing useful information similar to that from the frequencies below 183 GHz but with potentially enhanced spatial resolution from a satellite platform. In thin clouds, the dominant effect of water vapor is seen at 174, 340, and 400 GHz. In thick cloudy areas, the scattering effect is dominant at 90 and 220 GHz, while the overlaying water vapor can attenuate at 174, 340, and 400 GHz. All frequencies (90-400 GHz) show strong signatures in the core.

Prasad, N.↗
DOE OSTI2018

Equation of State Dependent Dynamics and Multi-messenger Signals from Stellar-mass Black Hole Formation

We investigate axisymmetric black hole (BH) formation and its gravitational wave (GW) and neutrino signals with self-consistent core-collapse supernova simulations of a non-rotating 40 M ⊙ progenitor star using the isotropic diffusion source approximation for the neutrino transport and a modified gravitational potential for general relativistic effects. We consider four different neutron star (NS) equations of state (EoS): LS220, SFHo, BHBΛ$\phi$, and DD2, and study the impact of the EoS on BH formation dynamics and GW emission. We find that the BH formation time is sensitive to the EoS from 460 to >1300 ms and is delayed in multiple dimensions for ~100–250 ms due to the finite entropy effects. Depending on the EoS, our simulations show the possibility that shock revival can occur along with the collapse of the proto-neutron star (PNS) to a BH. The gravitational waveforms contain four major features that are similar to previous studies but show extreme values: (1) a low-frequency signal (~300–500 Hz) from core-bounce and prompt convection, (2) a strong signal from the PNS g-mode oscillation among other features, (3) a high-frequency signal from the PNS inner-core convection, and (4) signals from the standing accretion shock instability and convection. The peak frequency at the onset of BH formation reaches to ~2.3 kHz. Thus, the characteristic amplitude of a 10 kpc object at peak frequency is detectable but close to the noise threshold of the Advanced LIGO and KAGRA, suggesting that the next-generation GW detector will need to improve the sensitivity at the kHz domain to better observe stellar-mass BH formation from core-collapse supernovae or failed supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Modeling Occupant Core Temperatures Across the Boston Building Stock to Advance Public Health

In recent history, extreme heat has been the cause of most deaths from a natural disaster. Exposure to extreme heat can aggravate preexisting conditions, increase hospitalization, and even cause death. Thus, modeling the thermal resilience of households across the United States will allow for a quantitative assessment of the health and safety risks posed by extreme heat. For the first time, we simulate occupant comfort in representative households across Boston by combining the granular results of the ResStock(TM) model with a two-node heat strain model. This model calculates occupants' core body temperature in each simulated household over a year. We compare the simulated core temperatures to two public health metrics: hyperthermia and heat stroke. Our results show for households in Boston that do not have or use aid conditioning, thousands potentially experience many dangerous heat events each summer and these events can last for more than a day at a time. These heat events peak during the late afternoon and evening, just as residents are coming home, cooking meals, and trying to go to sleep. We have found that multi-family buildings, renters, and low-income homes in more airtight and insulated homes are the most at risk for these events. These results demonstrate the scale and urgency of exposure to extreme heat.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2016

Quantifying Transient Uncertainty in the BEAVRS Benchmark Using Time Series Analysis Methods

Advances in computation have brought about significant improvements in creating fast-running high-fidelity simulations of nuclear cores. The BEAVRS benchmark is a highly-detailed PWR specification with two cycles of measured operational data used to validate high-fidelity core analysis methods. This PWR depletion benchmark captures the fine details of the LWR fuel assemblies, burnable absorbers, in-core fission detectors, core loading and shuffling patterns. Specifically, 58 of the 193 assemblies contain in-core detectors with measurements taken over 61 axial positions every month. These detectors are U-235 fission chambers with slightly varying mass of U-235. The collected signals are normalized on a given assembly permitting full core comparisons. The fuel layout for cycle 1 and instrument tube locations for the reactor are given in figures 1 and 2 respectively. Through a series of data processing and comparisons, it was shown that axially integrated radial maps of reaction rates were in close agreement between provided detector data and calculated data. More recently, the focus of the BEAVRS project has been on quantifying uncertainty to further assess the validity of aforementioned results. A close investigation of sources of error shows that uncertainties from operational nuclear data arise primarily from data measurements and processing. While the errors in axial realignment and interpolation seem significant, these sources affect only a small subset of data points near grid spacers and do not persist throughout the entire core. Instead, measurement uncertainty from detector count rates and radial integration dominate the uncertainty in measured data are determined as the most significant sources of error. This uncertainty analysis treated each given burnup step as independent of comparable data at neighboring burnup steps. More recent work has been targeted at understanding how reaction rates vary over time, in order to determine whether calculated reaction rates follow any observable trends. This abstract hones in on transient uncertainty quantification, in order to compare observed data against models for transient behavior. The first section of the paper uses linear regression tools to fit operational data, while the latter portion of the paper explores more complex simulation tools to fit operational data. Ultimately, the BEAVRS benchmark aims to serve as a true non-proprietary international benchmark for the validation of high-fidelity tools. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Fuel Cycle Performance of Fast Spectrum Molten Salt Reactor Designs

In the search for new ways to generate low-carbon, reliable base-load power, a resurgence of interest in advanced nuclear energy technologies, including Molten Salt Reactors (MSRs), has produced multiple new conceptual MSRs including fast neutron spectrum designs. The fuel cycle performance of four historical fast MSR designs is analyzed using a recently developed SCALE/TRITON 6.2.4 Alpha with a continuous online reprocessing functionality. The fast spectrum and continuous feed and removal of material enable these concepts to have remarkable fuel cycle metrics: (1) resource utilization is approximately 18 times better than for a typical low-enriched thermal spectrum once-through fuel cycle (i.e., from approximately 180 t/GWe-year to 1 t/GWe-year); (2) fast MSRs generate approximately 25 times less nuclear waste than the current once-through fuel cycle. These metrics are consistent with the Evaluation and Screening Study, which produced a technology-agnostic quantification of the characteristic performance of alternate fuel cycles. Additionally, full-core and unit cell transport models were created and compared to verify the viability of using simplified unit cell geometries for long-term depletion simulation. The unit cell approximation provided a speedup of 20 times relative to the full-core simulation, with depleted mass relative error for major isotopes of less than 2%. Additional fast MSRs design and analysis challenges associated with different fuel cycles and the use of MSR technology are addressed and discussed.

Rykhlevskii, Andrei↗
DOE OSTI2017

Thermal Hydraulic Computational Fluid Dynamics Simulations and Experimental Investigation of Deformed Fuel Assemblies

The project, Toward a Longer Life Core: Thermal Hydraulic CFD Simulations and Experimental Investigation of Deformed Fuel Assemblies, DOE Project code DE-NE0008321, was a verification and validation project for flow and heat transfer through wire wrapped simulated liquid metal fuel assemblies that included both experiments and computational fluid dynamics simulations of those experiments. This project was a two year collaboration between AREVA, TerraPower, Argonne National Laboratory and Texas A&M University. Experiments were performed by AREVA and Texas A&M University. Numerical simulations of these experiments were performed by TerraPower and Argonne National Lab. Project management was performed by AREVA Federal Services. The first of a kind project resulted in the production of both local point temperature measurements and local flow mixing experiment data paired with numerical simulation benchmarking of the experiments. The project experiments included the largest wire-wrapped pin assembly Mass Index of Refraction (MIR) experiment in the world, the first known wire-wrapped assembly experiment with deformed duct geometries and the largest numerical simulations ever produced for wire-wrapped bundles.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2012

Transmutation Analysis of Enriched Uranium and Deep Burn High Temperature Reactors

High temperature reactors (HTRs) have been under consideration for production of electricity, process heat, and for destruction of transuranics for decades. As part of the transmutation analysis efforts within the Fuel Cycle Research and Development (FCR&D) campaign, a need was identified for detailed discharge isotopics from HTRs for use in the VISION code. A conventional HTR using enriched uranium in UCO fuel was modeled having discharge burnup of 120 GWd/MTiHM. Also, a deep burn HTR (DB-HTR) was modeled burning transuranic (TRU)-only TRU-O2 fuel to a discharge burnup of 648 GWd/MTiHM. For each of these cases, unit cell depletion calculations were performed with SCALE/TRITON. Unit cells were used to perform this analysis using SCALE 6.1. Because of the long mean free paths (and migration lengths) of neutrons in HTRs, using a unit cell to represent a whole core can be non-trivial. The sizes of these cells were first set by using Serpent calculations to match a spectral index between unit cell and whole core domains. In the case of the DB-HTR, the unit cell which was arrived at in this way conserved the ratio of fuel to moderator found in a single block of fuel. In the conventional HTR case, a larger moderator-to-fuel ratio than that of a single block was needed to simulate the whole core spectrum. Discharge isotopics (for 500 nuclides) and one-group cross-sections (for 1022 nuclides) were delivered to the transmutation analysis team. This report provides documentation for these calculations. In addition to the discharge isotopics, one-group cross-sections were provided for the full list of 1022 nuclides tracked in the transmutation library.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Can we observe the ion-neutral drift velocity in prestellar cores?

Given the low-ionization fraction of molecular clouds, ambipolar diffusion is thought to be an integral process in star formation. However, chemical and radiative-transfer effects, observational challenges, and the fact that the ion-neutral drift velocity is inherently very small render a definite detection of ambipolar diffusion extremely non-trivial. Here, we study the ion-neutral drift velocity in a suite of chemodynamical, non-ideal magnetohydrodynamic (MHD), two-dimensional axisymmetric simulations of prestellar cores where we alter the temperature, cosmic-ray ionization rate, visual extinction, mass-to-flux ratio, and chemical evolution. Subsequently, we perform a number of non-local thermodynamic equilibrium (non-LTE) radiative-transfer calculations considering various idealized and non-idealized scenarios in order to assess which factor (chemistry, radiative transfer, and/or observational difficulties) is the most challenging to overcome in our efforts to detect the ion-neutral drift velocity. We find that temperature has a significant effect in the amplitude of the drift velocity with the coldest modelled cores (T = 6 K) exhibiting drift velocities comparable to the sound speed. Against expectations, we find that in idealized scenarios (where two species are perfectly chemically co-evolving) the drift velocity ‘survives’ radiative-transfer effects and can in principle be observed. However, we find that observational challenges and chemical effects can significantly hinder our view of the ion-neutral drift velocity. Finally, we propose that HCN and HCNH + ⁠, being chemically co-evolving, could be used in future observational studies aiming to measure the ion-neutral drift velocity.

79 ASTRONOMY AND ASTROPHYSICS↗