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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Uncertainty quantification for nuclear forensics with population analyses

Although neural networks offer cutting-edge predictive power, their deployment in high-consequence nuclear forensic applications is limited, partly because of their black-box nature. Incorporating robust uncertainty quantification methods into the predictive frameworks of neural networks is progress towards their future deployment in such scenarios. This work integrates uncertainty quantification into neural networks for nuclear reactor core-average burnup estimation from simulated environmental samples. We test two regimes (homogeneous and heterogeneous events) on DeepSets and Set Transformer architectures, we find both quantify predictive uncertainty effectively, but Set Transformer excels in partitioning latent events, offering superior predictive power and more informative uncertainty estimates.

Hatton, Conner [ORNL] (ORCID:0009000804970959)↗
DOE OSTI2019

Experimental validation of a nuclear forensics methodology for source reactor-type discrimination of chemically separated plutonium

An experimental validation of a nuclear forensics methodology for the source reactor-type discrimination of separated weapons-useable plutonium is presented. The methodology uses measured values of intra-element isotope ratios of plutonium and fission product contaminants. MCNP radiation transport codes were used for various reactor core modeling and fuel burnup simulations. A reactor-dependent library of intra-element isotope ratio values as a function of burnup and time since irradiation was created from the simulation results. The experimental validation of the methodology was achieved by performing two low-burnup experimental irradiations, resulting in distinct fuel samples containing sub-milligram quantities of weapons-useable plutonium. The irradiated samples were subjected to gamma and mass spectrometry to measure several intra-element isotope ratios. For each reactor in the library, a maximum likelihood calculation was utilized to compare the measured and simulated intra-element isotope ratio values, producing a likelihood value which is proportional to the probability of observing the measured ratio values, given a particular reactor in the library. The measured intra-element isotope ratio values of both irradiated samples and its comparison with the simulation predictions using maximum likelihood analyses are presented. Finally, the analyses validate the nuclear forensics methodology developed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2013

Computational and experimental prediction of dust production in pebble bed reactors, Part II

This paper is the continuation of Part I, which describes the high temperature and high pressure helium environment wear tests of graphite–graphite in frictional contact. In the present work, it has been attempted to simulate a Pebble Bed Reactor core environment as compared to Part I. The experimental apparatus, which is a custom-designed tribometer, is capable of performing wear tests at PBR relevant higher temperatures and pressures under a helium environment. This environment facilitates prediction of wear mass loss of graphite as dust particulates from the pebble bed. The experimental results of high temperature helium environment are used to anticipate the amount of wear mass produced in a pebble bed nuclear reactor.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2018

Assessment of the analysis capability for core-wide PWR pellet-clad interaction screening of Watts Bar Unit 1

This report illustrates the application of the Virtual Environment for Reactor Applications simulation environment to pressurized water reactor core-wide pellet-clad interaction screening, with emphasis on the multi-dimensional analysis capabilities of the BISON fuel performance code. The Virtual Environment for Reactor Applications codes utilized in this effort include the BISON fuel performance code and the VERA-CS sub-element, which combines the MPACT neutron transport code and the COBRA-TF thermal-hydraulics code. The target reactor for this work is Watts Bar Unit 1, cycle 6 and cycle 7, selected on the basis that pellet-clad interaction failures had occurred during power start up in both cycles, despite the moderate power levels. These instances of failure at Watts Bar Unit 1 offered an interesting test of the code’s veracity and a particularly challenging condition for this assessment.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Numerical solution of large scale Hartree–Fock–Bogoliubov equations

The Hartree–Fock–Bogoliubov (HFB) theory is the starting point for treating superconducting systems. However, the computational cost for solving large scale HFB equations can be much larger than that of the Hartree–Fock equations, particularly when the Hamiltonian matrix is sparse, and the number of electrons N is relatively small compared to the matrix size N b . We first provide a concise and relatively self-contained review of the HFB theory for general finite sized quantum systems, with special focus on the treatment of spin symmetries from a linear algebra perspective. We then demonstrate that the pole expansion and selected inversion (PEXSI) method can be particularly well suited for solving large scale HFB equations. For a Hubbard-type Hamiltonian, the cost of PEXSI is at most $\mathcal{O}$(N b 2 ) for both gapped and gapless systems, which can be significantly faster than the standard cubic scaling diagonalization methods. We show that PEXSI can solve a two-dimensional Hubbard-Hofstadter model with N b up to 2.88 × 10 6 , and the wall clock time is less than 100 s using 17 280 CPU cores. Finally, this enables the simulation of physical systems under experimentally realizable magnetic fields, which cannot be otherwise simulated with smaller systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Photoinduced molecular chirality probed by ultrafast resonant X-ray spectroscopy

Recently developed circularly polarized X-ray light sources can probe the ultrafast chiral electronic and nuclear dynamics through spatially localized resonant core transitions. Here, we present simulations of time-resolved circular dichroism signals given by the difference of left and right circularly polarized X-ray probe transmission following an excitation by a circularly polarized optical pump with the variable time delay. Application is made to formamide which is achiral in the ground state and assumes two chiral geometries upon optical excitation to the first valence excited state. Probes resonant with various K-edges (C, N, and O) provide different local windows onto the parity breaking geometry change thus revealing the enantiomer asymmetry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Two extreme young objects in Barnard 1-b

Two submillimeter/millimeter sources in the Barnard 1b (B1-b) core, B1-bN and B1-bS, have been studied in dust continuum, H{sup 13}CO{sup +} J = 1-0, CO J = 2-1, {sup 13}CO J = 2-1, and C{sup 18}O J = 2-1. The spectral energy distributions of these sources from the mid-IR to 7 mm are characterized by very cold temperatures of T {sub dust} < 20 K and low bolometric luminosities of 0.15-0.31 L {sub ☉}. The internal luminosities of B1-bN and B1-bS are estimated to be <0.01-0.03 L {sub ☉} and ∼0.1-0.2 L {sub ☉}, respectively. Millimeter interferometric observations have shown that these sources have already formed central compact objects of ∼100 AU sizes. Both B1-bN and B1-bS are driving the CO outflows with low characteristic velocities of ∼2-4 km s{sup –1}. The fractional abundance of H{sup 13}CO{sup +} at the positions of B1-bN and B1-bS is lower than the canonical value by a factor of four to eight. This implies that a significant fraction of CO is depleted onto dust grains in the dense gas surrounding these sources. The observed physical and chemical properties suggest that B1-bN and B1-bS are in an earlier evolutionary stage than most of the known class 0 protostars. In particular, the properties of B1-bN agree with those of the first hydrostatic core predicted by the MHD simulations. The CO outflow was also detected in the mid-IR source located at ∼15'' from B1-bS. Since the dust continuum emission was not detected in this source, the circumstellar material surrounding this source is less than 0.01 M {sub ☉}. It is likely that the envelope of this source was dissipated by the outflow from the protostar that is located to the southwest of B1-b.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

A comparison of different continuum approaches in modeling mixed-type dislocations in Al

Mixed-type dislocations are common in metals and play a vital role in their plastic deformation. Key characteristics of mixed-type dislocations cannot simply be extrapolated from those of dislocations with pure edge or pure screw characters. Yet, mixed-type dislocations traditionally received disproportionately less attention in the modeling and simulation community. Here, we explore core structures of mixed-type dislocations in Al using three continuum approaches, namely, the phase-field dislocation dynamics (PFDD) method, the atomistic phase-field microelasticity (APFM) method, and the concurrent atomistic-continuum (CAC) method. Results are benchmarked against molecular statics. We advance the PFDD and APFM methods in several aspects such that they can better describe the dislocation core structure. In particular, in these two approaches, the gradient energy coefficients for mixed-type dislocations are determined based on those for pure-type ones using a trigonometric interpolation scheme, which is shown to provide better prediction than a linear interpolation scheme. The dependence of the in-slip-plane spatial numerical resolution in PFDD and CAC is also quantified.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Probabilistic structure of the geodynamo

One of the most intriguing features of Earth's axial magnetic dipole field, well known from the geological record, is its occasional and unpredictable reversal of polarity. Understanding the phenomenon is rendered very difficult by the highly nonlinear nature of the underlying magnetohydrodynamic problem. Numerical simulations of the liquid outer core, where regeneration occurs, are only able to model conditions that are far from Earth-like. On the analytical front, the situation is not much better; basic calculations, such as relating the average rate of reversals to various core parameters, have apparently been intractable. Here, we present a framework for solving such problems. Starting with the magnetic induction equation, we show that by considering its sources to be stochastic processes with fairly general properties, we can derive a differential equation for the joint probability distribution of the dominant toroidal and poloidal modes. This can be simplified to a Fokker-Planck equation and, with the help of an adiabatic approximation, reduced even further to an equation for the dipole amplitude alone. From these equations various quantities related to the magnetic field, including the average reversal rate, field strength, and time to complete a reversal, can be computed as functions of a small number of numerical parameters. As a result, these parameters in turn can be computed from physical considerations or constrained by paleomagnetic, numerical, and experimental data.

58 GEOSCIENCES↗
DOE OSTI2015

A voxel-based multiscale model to simulate the radiation response of hypoxic tumors

Purpose: In radiotherapy, it is important to predict the response of tumors to irradiation prior to the treatment. This is especially important for hypoxic tumors, which are known to be highly radioresistant. Mathematical modeling based on the dose distribution, biological parameters, and medical images may help to improve this prediction and to optimize the treatment plan. Methods: A voxel-based multiscale tumor response model for simulating the radiation response of hypoxic tumors was developed. It considers viable and dead tumor cells, capillary and normal cells, as well as the most relevant biological processes such as (i) proliferation of tumor cells, (ii) hypoxia-induced angiogenesis, (iii) spatial exchange of cells leading to tumor growth, (iv) oxygen-dependent cell survival after irradiation, (v) resorption of dead cells, and (vi) spatial exchange of cells leading to tumor shrinkage. Oxygenation is described on a microscopic scale using a previously published tumor oxygenation model, which calculates the oxygen distribution for each voxel using the vascular fraction as the most important input parameter. To demonstrate the capabilities of the model, the dependence of the oxygen distribution on tumor growth and radiation-induced shrinkage is investigated. In addition, the impact of three different reoxygenation processes is compared and tumor control probability (TCP) curves for a squamous cells carcinoma of the head and neck (HNSSC) are simulated under normoxic and hypoxic conditions. Results: The model describes the spatiotemporal behavior of the tumor on three different scales: (i) on the macroscopic scale, it describes tumor growth and shrinkage during radiation treatment, (ii) on a mesoscopic scale, it provides the cell density and vascular fraction for each voxel, and (iii) on the microscopic scale, the oxygen distribution may be obtained in terms of oxygen histograms. With increasing tumor size, the simulated tumors develop a hypoxic core. Within the model, tumor shrinkage was found to be significantly more important for reoxygenation than angiogenesis or decreased oxygen consumption due to an increased fraction of dead cells. In the studied HNSSC-case, the TCD{sub 50} values (dose at 50% TCP) decreased from 71.0 Gy under hypoxic to 53.6 Gy under the oxic condition. Conclusions: The results obtained with the developed multiscale model are in accordance with expectations based on radiobiological principles and clinical experience. As the model is voxel-based, radiological imaging methods may help to provide the required 3D-characterization of the tumor prior to irradiation. For clinical application, the model has to be further validated with experimental and clinical data. If this is achieved, the model may be used to optimize fractionation schedules and dose distributions for the treatment of hypoxic tumors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Attention to quantum complexity

The imminent era of error-corrected quantum computing demands robust methods to characterize quantum state complexity from limited, noisy measurements. We introduce the Quantum Attention Network (QuAN), a classical artificial intelligence (AI) framework leveraging attention mechanisms tailored for learning quantum complexity. Inspired by large language models, QuAN treats measurement snapshots as tokens while respecting permutation invariance. Combined with our parameter-efficient miniset self-attention block, this enables QuAN to access high-order moments of bit-string distributions and preferentially attend to less noisy snapshots. We test QuAN across three quantum simulation settings: driven hard-core Bose-Hubbard model, random quantum circuits, and toric code under coherent and incoherent noise. QuAN directly learns entanglement and state complexity growth from experimental computational basis measurements, including complexity growth in random circuits from noisy data. In regimes inaccessible to existing theory, QuAN unveils the complete phase diagram for noisy toric code data as a function of both noise types, highlighting AI’s transformative potential for assisting quantum hardware.

Kim, Hyejin [Cornell Univ., Ithaca, NY (United Sta↗
DOE OSTI2018

Quasi-Stable Configurations of Torus Vortex Knots and Links

The dynamics of torus vortex configurations V{sub n, p, q} in a superfluid liquid at zero temperature (n is the number of quantum vortices, p is the number of turns of each filament around the symmetry axis of the torus, and q is the number of turns of the filament around its central circle; radii R{sub 0} and r{sub 0} of the torus at the initial instant are much larger than vortex core width ξ) has been simulated numerically based on the regularized Biot–Savart law. The lifetime of vortex systems till the instant of their substantial deformation has been calculated with a small step in parameter B{sub 0} = r{sub 0}/R{sub 0} for various values of parameter Λ = ln(R{sub 0}/ξ). It turns out that for certain values of n, p, and q, there exist quasi-stability regions in the plane of parameters (B{sub 0}, Λ), in which the vortices remain almost invariable during dozens and even hundreds of characteristic lengths.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

A method for development of efficient 3D models for neutronic calculations of ASTRA critical facility using experimental information

The application of experimental information on measured axial distributions of fission reaction rates for development of 3D numerical models of the ASTRA critical facility taking into account azimuthal asymmetry of the assembly simulating a HTGR with annular core is substantiated. Owing to the presence of the bottom reflector and the absence of the top reflector, the application of 2D models based on experimentally determined buckling is impossible for calculation of critical assemblies of the ASTRA facility; therefore, an alternative approach based on the application of the extrapolated assembly height is proposed. This approach is exemplified by the numerical analysis of experiments on measurement of efficiency of control rods mockups and protection system (CPS).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Communication-Efficient Property Preservation in Tracer Transport

Atmospheric tracer transport is a computationally demanding component of the atmospheric dynamical core of weather and climate simulations. Simulations typically have tens to hundreds of tracers. A tracer field is needed to preserve several properties, including mass, shape, and tracer consistency. To improve computational efficiency, it is common to apply different spatial and temporal discretizations to the tracer transport equations than to the dynamical equations. Using different discretizations increases the difficulty of preserving properties. This paper provides a unified framework to analyze the property preservation problem and classes of algorithms to solve it. We observe the primary problem and a safety problem; describe three classes of algorithms to solve these; introduce new algorithms in two of these classes; make connections among the algorithms; analyze each algorithm in terms of correctness, bound on its solution magnitude, and its communication efficiency; and study numerical results. A new algorithm, QLT, has the smallest communication volume, and in an important case it redistributes mass approximately locally. These algorithms are only very loosely coupled to the underlying discretizations of the dynamical and tracer transport equations and thus are broadly and efficiently applicable. In addition, they may be applied to remap problems in applications other than tracer transport.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Consist v0.1.0

A Python library for provenance tracking, intelligent caching, and data virtualization in scientific simulation workflows. It automatically records code, configuration, and input data to skip redundant computations and enables querying results across many runs without manual bookkeeping. Designed to support multi-model simulation workflows like the BEAM CORE toolset at LBL, but designed to be extensible to a wide range of research workflows. Combines lineage tracking features as provided by OpenLineage with deterministic hashing like SnakeMake, and adds powerful analysis tools on model outputs.

Needell, Zachary [Lawrence Berkeley National Labor↗
DOE OSTI2019

Interaction of the Westerlies with the Tibetan Plateau in Determining the Mei-Yu Termination

This article explores how the termination of the mei-yu is dynamically linked to the westerlies impinging on the Tibetan Plateau. It is found that the mei-yu stage terminates when the maximum upper-tropospheric westerlies shift beyond the northern edge of the plateau, around 40°N. This termination is accompanied by the disappearance of tropospheric northerlies over northeastern China. The link between the transit of the jet axis across the northern edge of the plateau, the disappearance of northerlies, and termination of the mei-yu holds on a range of time scales from interannual through seasonal and pentad. Diagnostic analysis indicates that the weakening of the meridional moisture contrast and meridional wind convergence, mainly resulting from the disappearance of northerlies, causes the demise of the mei-yu front. The authors propose that the westerlies migrating north of the plateau and consequent weakening of the extratropical northerlies triggers the mei-yu termination. Model simulations are employed to test the causality between the jet and the orographic downstream northerlies by repositioning the northern edge of the plateau. As the plateau edge extends northward, orographic forcing on the westerlies strengthens, leading to persistent strong downstream northerlies and a prolonged mei-yu. Idealized simulations with a dry dynamical core further demonstrate the dynamical link between the weakening of orographically forced downstream northerlies with the positioning of the jet from south to north of the plateau. Changes in the magnitude of orographically forced stationary waves are proposed to explain why the downstream northerlies disappear when the jet axis migrates beyond the northern edge of the plateau.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

A Molecular Dynamics Study of Subsurface Hydrogen-Helium Bubbles in Tungsten

Fusion reactor materials experience high ion fluxes and operating temperatures, which will ultimately produce subsurface helium and hydrogen bubbles in the tungsten divertor that can cause surface degradation and impact core plasma performance. Molecular dynamics simulations have been used to evaluate the behavior of hydrogen and helium near a 2-nm bubble or void below a tungsten surface as a function of surface orientation, temperature, gas atom concentration, initial hydrogen distribution, and depth below the surface. A clear tendency for hydrogen to segregate to the bubble-matrix interface is observed in these simulations, regardless of the initial spatial distribution of the hydrogen or simulation parameters. This segregation is due in part to a local minimum in the hydrogen energy at the periphery of the bubble. Further work is required to fully characterize the mechanism of this behavior and to assess the quantities of hydrogen in the bubble and at the bubble periphery.

Bergstrom, Z. J.↗
DOE OSTI2013

Full Core Reactor Analysis: Running Denovo on Jaguar

Fully-consistent, full-core, 3D, deterministic neutron transport simulations using the orthogonal mesh code Denovo were run on the massively parallel computing architecture Jaguar XT5. Using energy and spatial parallelization schemes, Denovo was able to efficiently scale to more than 160k processors. Cell-homogenized cross sections were used with step-characteristics, linear-discontinuous finite element, and trilinear-discontinuous finite element spatial methods. It was determined that using the finite element methods gave considerably more accurate eigenvalue solutions for large-aspect ratio meshes than using step-characteristics.

Denovo↗