Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Singularity”

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 829 records · Page 46

NA-213 Presentation.pdf

The Savannah River National Laboratory (SRNL) is one of 17 United States Department of Energy (DOE) National Laboratories. SRNL, a multi-program national laboratory, is a leading research and development institution for the Offices of Environmental Management and Legacy Management at the U.S. Department of Energy and the Weapons and Nonproliferation programs for the National Nuclear Security Administration. Battelle Savannah River Alliance, LLC (BSRA), a not-for-profit limited liability company, manages and operates SRNL for the DOE. BSRA board leadership includes Battelle Memorial Institute, Clemson University, University of South Carolina, South Carolina State University, University of Georgia, and Georgia Institute of Technology. Battelle Memorial Institute and the five universities are joined in partnership with preferred subcontractors TechSource and Longenecker & Associates with the singular purpose of maintaining SRNL as a best-in-class national laboratory. Scientists and Engineers at SRNL use leading edge science and technology to advance the Department of Energy’s critical mission outcomes.

Hasty, Timothy C.↗

Managing Hanford's Direct Feed Low Activity Waste Program - 20002

The Office of River Protection (ORP) is responsible for the management and completion of the River Protection Project (RPP) mission, which comprises both the Hanford Site tank farms operations and the completion and operation of the Waste Treatment and Immobilization Plant Project. The RPP mission is to safely retrieve, treat, and immobilize Hanford's high-level and low-activity tank waste and close the tank farms to protect the Columbia River. The Direct-Feed Low-Activity Waste (DFLAW) Program is a major subset of the overall RPP mission, providing disposition of low-activity waste. The scope is to retrieve, treat, and immobilize the low-activity waste at the Low-Activity Waste Facility. DFLAW will begin operations no later than December 2023. The DFLAW Program overarches a suite of individual projects, activities, and infrastructure upgrades. The major projects and activities include multiple contractors under both ORP and the Richland Operations Office, underscoring the need for effective integration, coordination, and collaboration. The objectives of the DFLAW Program leadership team are to: - Effectively coordinate and integrate the projects that comprise the DFLAW Program - Manage the interfaces between the projects so that the integrated DFLAW Program is completed successfully - Ensure the DFLAW portfolio of projects operate as required without gaps or conflicts at the interfaces. DOE is accomplishing these objectives through a newly established leadership model that integrates the DOE and contractor work streams and promotes a teamwork model to enable collaborative success. Long considered one of the most formidable cleanup challenges at Hanford, the Department of Energy and its contractors are on the verge of achieving a cleanup commitment that has been decades in the making. DOE is preparing to vitrify (turn to glass) Hanford's chemical and radioactive tank waste using the DFLAW process. The Modified Consent Decree milestone for completing hot commissioning of the Hanford Site's Low- Activity Waste Facility is Dec. 31, 2023. While this milestone may seem distant under normal project management circumstances, it is a relatively short period of time for the operational and cultural transformation necessary to successfully begin treating tank waste at Hanford. Direct feed means separating the waste at a tank farm to remove the more radioactive portion (i.e., solids and cesium) so that the resulting low-activity (less radioactive) waste can be fed directly to the Waste Treatment and Immobilization Plant's Low-Activity Waste Facility. Starting DFLAW will require a singular, intensive leadership focus, along with a sustained, collective commitment to excellence and teamwork by the Department and its contractors. Major upgrades to Hanford's infrastructure will occur prior to startup. Supporting DFLAW operations requires significant integration between the DFLAW facilities - the Waste Treatment Plant, Hanford's tank farms, effluent treatment facilities, disposal facilities - and between all of the site contractors to achieve success. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Oxidation of Glycolate in the Defense Waste Processing Facility (DWPF) Recycle Collection Tank - 20305

The Savannah River Site's Defense Waste Processing Facility (DWPF) operations are being upgraded with the introduction of the Nitric-Glycolic Flowsheet. Glycolic acid has been shown superior to formic acid as the reducing acid used during chemical processing. The new flowsheet improves or maintains necessary parameters such as 1) reduction of mercury, 2) adjustment of feed rheology and 3) adjustment of melter oxidation/reduction potential. Further, the potential for catalytic hydrogen generation in DWPF processing is virtually eliminated. DWPF process condensates are collected and returned to the SRS Concentration, Storage and Transfer Facilities (CSTF). The Recycle Collection Tank (RCT) collects off-gas condensate during chemical processing, vitrification, and other unit operations performed in DWPF and is the singular return vessel delivering recycle effluent back to CSTF. Each batch of recycle may contain a small amount of glycolate from chemical processing and melter off-gas condensates. To avoid potential flammability issues due to thermolysis of glycolate in the CSTF, chemical oxidation within the RCT has been investigated as an option for mitigating the transfer of glycolate. Sodium permanganate has been down-selected as the best option for oxidation of glycolate. Testing was performed using both 2-L and 22-L reactors (16,800:1 and 1,530:1 scale by volume) with non-radioactive waste simulants to approximate the expected RCT compositions. RCT simulants were evaluated at various process pH and temperature conditions. Also, RCT operations, namely the sequence of addition of corrosion inhibitors (NaOH and NaNO{sub 2}) versus a permanganate strike, were evaluated. Glycolate was introduced via a sludge simulant to mimic both expected entrainment and abnormal process foam-over conditions - the range being between 68 and 5100 mg/kg glycolate. Glycolate destruction was monitored by ion chromatography (IC). The corresponding manganese behavior was monitored in real-time using in situ ultraviolet-visible (UV-Vis) spectroscopy. RCT glycolate content can be reduced to below the IC detection limit within 90 minutes for all concentrations investigated. Ion Chromatography analysis revealed that under alkaline conditions, glycolate is primarily oxidized to oxalate with no significant formation of CO{sub 2} or carbonate, and nitrite is not oxidized to nitrate. Initially, complete oxidation of organics species and nitrite was assumed. Determination of the mechanistic chemical reaction has allowed the required amount of permanganate to be more accurately predicted and the total addition to be significantly reduced. UV-Vis measurements reveal that permanganate (Mn{sup 7+}) is reduced to manganate (Mn{sup 6+}) in the RCT. The oxidant stoichiometry is defined by using the initial permanganate to glycolate (P/G) molar ratio. At low initial glycolate concentration (68 and 140 mg/kg), the minimum required initial permanganate to glycolate (P/G) molar ratio was found to be 5-6. With high initial glycolate concentrations (5100 mg/kg) a lower (P/G) molar ratio of ∼2.5 was needed. The final portion of this effort supporting the nitric/glycolic flowsheet will be to test actual (fully radioactive) RCT samples as per the above simulant tests. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Transporting Non-Compliant TRU Wastes in the OPTIMUS{sup TM} Packages - 20493

The OPTIMUS{sup TM} product line of packages was developed with the purpose of providing a versatile and modular packaging option for shipping problematic wastes and fissile material contents. The product line includes a high-activity design in the OPTIMUS-H and a low activity design in the OPTIMUS-L. The two variants utilize a singular containment vessel design, that is enclosed in different protective components that offer different levels of shielding along with thermal and impact protection. The thick shielding provided by the OPTIMUS-H packaging allows for higher activity contents than the OPTIMUS-L to be transported. But the lightweight design of the OPTIMUS-L packaging allows for more packages to be transported in a single shipment. The containment vessel design shared between the two packages provides leak-tight containment of all radioactive contents along with the capabilities to easily leak test and the option to inert and backfill the contents, as necessary, prior to each shipment of the package. Together the OPTIMUS package designs offer significant flexibility and are capable of handling a wide variety of waste materials and other radioactive contents. One of the primary contents covered in the initial design effort for these packages is TRU waste contents beyond the standard waste materials that are compliant with the WIPP Waste Acceptance Criteria. More specifically, the initial content of interest is TRU Waste drums containing sealed containers with potentially flammable gases. While there are multiple other packaging options available for transporting standard WIPP compliant TRU wastes, the goal of the OPTIMUS packages is to offer a superior option for transporting both standard TRU wastes as well as the more problematic non-compliant wastes, among other contents. The non-compliant TRU waste contents explicitly included for the OPTIMUS packages are standard aerosol cans and DOT 3E lecture bottles. However, the methods applied for the package containment can be easily adjusted to cover other sealed container types with potentially flammable gases. The primary challenges with including these items are both the obvious issue of potential for flammable gases (e.g. aerosol propellants) in the package, but also the uncertainty in the exact state of the contents. Because these desired contents are waste materials, the characterization of the materials present in the waste may include some uncertainties. For example, though it may be known that there are one or more aerosol cans in a TRU waste drum, it may not be known if this can is full, spent, or anywhere in between. Also there may be equal uncertainty in the contents of the aerosol can, specifically in the potentially flammable propellant remaining in the can. To include these types of non-compliant items as acceptable contents for the package, the methods utilized in the safety analyses of the OPTIMUS packages must consider all of the uncertainties in the characterization of the waste. The methods utilized cover a range of concerns with transporting contents of this nature including pressure buildup, gas generation, and gaseous combustion for demonstrating the containment of the package. The ability to transport these problematic contents in the OPTIMUS packages gives TRU waste generating sites the option to relocate these drums to an offsite location where the non-compliant items can be properly managed through removal or destruction. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Neutron diffusion calculation in heterogeneous geometry based on local/global iteration using proper orthogonal decomposition

This study newly proposes a heterogeneous core calculation method based on local/global iteration using proper orthogonal decomposition (POD). By using the singular value decomposition (SVD) and the low-rank approximation, appropriate POD bases for expanding the neutron flux can be obtained from snapshot data of the neutron flux obtained by fine mesh calculations. By projection using the POD bases, the dimension of the target equation (e.g., discretized neutron diffusion equation) can be dramatically reduced. In the proposed method, POD is effectively applied to each single assembly calculation (local calculation). Furthermore, using the local/global iteration, the effective neutron multiplication factor and the neutron flux distribution in the whole core geometry can be obtained by combining the numerical results of the local calculation for each fuel assembly and the global calculation for the whole core. As a feasibility study, the proposed method is applied to a one-dimensional heterogeneous core analysis, and the accuracy is investigated by changing the total number of POD bases. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Quantum Simulation of an Open System: A Dissipative 1+1D Ising Model

The 1+1D Ising model is an ideal benchmark for quantum algorithms, as it is very well understood theoretically. This is true even when expanding the model to include complex coupling constants. In this work, we implement quantum algorithms designed for the simulation of open or complex coupling quantum field theories on IBM devices with a focus on the measurement of the Lee-Yang edge singularity. This feature corresponds (at large volumes) to a phase transition, and our successful reproduction of the transition represents a non-trivial test for current hardware and its ability to distinguish features of interest in quantum field theories.

Gustafson, Erik↗

Powering Data Centers with Clean Energy: A Techno-Economic Case Study of Nuclear and Renewable Energy Dependability

Rising data demands from artificial intelligence (AI) and large language models (LLMs) generating images, videos, and text have prompted increased need for larger and more robust data centers in the United States. Major companies interested in these larger data centers face the choice of linking them to existing regional grids, building stand-alone power supplies onsite, or a combination of both. The request, review, and approval process for new transmission lines to grids in the United States, however, has grown in recent years to times spans rivaling those of new construction for nuclear power plants. Building an islanded power supply for each data center is therefore becoming a prominent option. In this case study, several technologies are modeled in techno-economic simulations for long-term system costs subject to fixed electricity demand from a singular data center. A 250 MWe data center is assumed with additional 50 MWe for resiliency. Techno-economic simulations are conducted using the Holistic Energy Resource Optimization Network (HERON) software, which is a part of the Framework for Optimization of Resources and Economics (FORCE) tool suite. Technologies considered include solar, wind, lithium-ion batteries, and several types of nuclear reactors: large-scale reactors, small modular reactors, and microreactors. A low- and high-cost estimate for each technology is assumed to develop a range of expected economic performance. Low-cost estimates included several clean energy production tax credits. Different combinations of renewable energy generators with nuclear reactors are considered, ranging from a fully renewable-powered data center to a fully nuclear-powered data center. Historic time series of wind and solar availability from the Texas grid are used to train a reduced order model; this model then generates unique time series with similar characteristics of the training dataset. Multiple scenarios of weather and subsequent operations are simulated for each renewable-nuclear combination to determine total costs throughout the project lifetime. Fully renewable-powered configurations required large amounts of installed capacity (GW scale) in the simulations to meet the fixed demand of the data center. This is due to some scenarios in the historical dataset which captured low-wind and low-solar days, requiring over-building of these technologies as well as batteries to compensate for the low amounts of electricity generation. Fully nuclear-powered configurations outperformed the fully renewable and mixed renewable-nuclear configurations in terms of cost, with ranges between $1B and $10B in 2023 USDs compared to $40B+ for fully renewable configurations. Of the nuclear technologies, small modular reactors performed better economically than large-scale nuclear models due to lower projected capital costs, and both performed better than the microreactor models. These results demonstrate the applicability of firm, dispatchable electricity resources from baseload generators like nuclear power plants for operating facilities that run at constant power without daily variability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Regularized Differentiation for Bioburden Density Estimation in Planetary Protection

In this paper, we propose and investigate the performance of two novel shrinkage estimators for bioburden density estimation in planetary protection. The estimators are based on the regularized differentiation of a cumulative count of colony forming units collected throughout the data collecting session or the life cycle of the entire mission. The regularized differentiation recasts the problem of bioburden density estimation as a linear least squares problem. The least squares problem is then solved through regularization techniques, such as truncated singular value decomposition and penalized least squares. The regularization is necessary to avoid noise amplification during the differentiation of noisy data. The two regularization estimators are compared with four other commonly used estimators to simultaneously evaluate the means of multivariable independent Poisson distributions: the maximum likelihood, noninformative Bayes estimator with Jeffreys prior, Empirical Bayes using conjugate gamma-Poisson model with gamma parameters selected by method of moments, and the Clevenson-Zidek estimator. It is shown through computer-simulated data that the regularized differentiation based on ridge regression has the smallest mean-squared error among all estimators. The analysis of shrinkage mechanism implemented by regularized differentiation is performed, and it is shown that the regularized differentiation amounts to performing a weighted averaging of all the samples. The weights are determined by the regularization parameter automatically selected by the L-curve technique. Since the method of least squares makes no distributional assumptions about the data, it presents an attractive technique for bioburden density estimation when there are concerns about the misspecification of the distributional model. The paper concludes with the analysis of the bioburden data collected during InSight mission and directions for future work.

97 - MATHEMATICS AND COMPUTING↗

Generative models on phase space

Deep generative models such as diffusion and flow matching are powerful machine learning tools capable of learning and sampling from high-dimensional distributions. They are particularly useful when the training data appears to be concentrated on a submanifold of the data embedding space. For high-energy physics data, consisting of collections of relativistic energy-momentum 4-vectors, this submanifold can enforce extremely strong physically-motivated priors, such as energy and momentum conservation. If these constraints are learned only approximately, rather than exactly, this can inhibit the interpretability and reliability of such generative models. To remedy this deficiency, we introduce generative models which are, by construction, confined at every step of their sampling trajectory to the manifold of massless N-particle Lorentz-invariant phase space in the center-of-momentum frame. In the case of diffusion models, the "pure noise" forward process endpoint corresponds to the uniform distribution on phase space, which provides a clear starting point from which to identify how correlations among the particles emerge during the reverse (de-noising) process. We demonstrate that our models are able to learn both few-particle and many-particle distributions with various singularity structures, paving the way for future interpretability studies using generative models trained on simulated jet data.

Bogorad, Zachary [Fermilab]↗

The Verification and Validation of a Magnetic Plasma Fluid Model Utilizing the MOOSE (Multiphysics Object Oriented Simulation Environment) Framework

As the goal of achieving fusion power on the grid comes closer to fruition, fully coupled multiphysics models of fusion devices will be crucial. Currently, there are two main approaches to developing these platforms: (1) loosely coupled, where one couples existing codes and solvers together through input and output parameters and data, and (2) tightly coupled, where one develops the necessary models within a singular, integrated framework. This work focuses on the latter approach for magnetically confined fusion devices by developing a fluid-based plasma-edge model within the Multiphysics Object Oriented Simulation Environment (MOOSE) Framework. This effort is coordinated with other efforts to develop, test, demonstrate, and deploy fusion relevant multiphysics capabilities including electromagnetics, particle-in-cell plasma, tritium transport, and fusion blanket design. This new model is an expansion of the MOOSE-based plasma application, Zapdos, which was originally formulated to model low-temperature, non-magnetized plasma processes. Verification, benchmarking, and validation studies have been conducted. Verification studies involved utilizing the method of manufactured solutions and comparing the convergence slope of a known solution to the theoretical slope. Benchmarking consists of comparisons to existing edge codes, namely BOUT++ and SOLEDGE3X. Validation efforts focused on comparisons against open-source data from the TCV tokamak.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Double-real corrections to color singlet decay in a parton-shower inspired scheme

We introduce a local infrared subtraction method for next-to-next-to-leading order QCD calculations in color singlet decays, with counterterms based on scalar radiators and pure splitting functions. Overlapping singularities in the multipole radiation pattern are disentangled by partial fractioning, and the kinematics mapping corresponds to iterated next-to-leading order kinematics. We verify that the double-real remainder to $e^+e^-\to\;q\bar{q}$ is rendered finite in the single and double unresolved limits and investigate the numerical convergence of the Monte-Carlo integral. We compute the phase-space integrals of the scalar counterterms in the back-to-back configuration, both analytically and with the help of numerical techniques based on sector decomposition.

Campbell, John M. [Fermilab]↗

Precise 2D electric field density simulations for superconducting quantum devices

Dielectric loss due to two-level systems is a limiting factor for superconducting qubit relaxation times. These losses arise mostly from nanometer-scale interfacial defect regions in superconducting devices with planar dimensions of microns to millimeters, thus making it resource intensive to accurately simulate the electric field density in these regions with traditional electromagnetic solvers. In this work, we demonstrate a fast boundary integral equation solver that allows precise simulation of electric field density in these thin regions, showing a speedup of around two orders of magnitude over traditional solvers, with relative errors around $10^{-7}$ for a ten-minute solution runtime. By computing participation ratios through Green's first identity without squaring the electric field, our approach is less susceptible to the field singularities near conductor corners. We apply this solver to a basic untrenched coplanar waveguide cross-section, showing that the common assumption of participation ratio linearity with dielectric constant holds well for some interfaces and not others; in particular, while the metal-air (MA) top and corner follow this linear relationship strongly, the MA sidewall does not. We then compare isotropic and anisotropic etching, showing that the MA sidewall and the metal-air-substrate triple junction are the most strongly affected. We are currently leveraging this solver to explore geometries that will uniquely isolate the participation ratios of the different dielectrics. Finally, we are working to combine this solver framework with a full 3D microwave solver to accurately calculate participation ratios for the thin dielectrics that are known sources of loss in superconducting qubits.

Gimbutas, Z. [NIST, Boulder] (ORCID:00000003320982↗

Encoding beyond cosmological horizons in de Sitter JT gravity

Black hole event horizons and cosmological event horizons share many properties, making it natural to ask whether our recent advances in understanding black holes generalize to cosmology. To this end, we discuss a paradox that occurs if observers can access what lies beyond their cosmological horizon in the same way that they can access what lies beyond a black hole horizon. In particular, distinct observers with distinct horizons may encode the same portion of spacetime, violating the no-cloning theorem of quantum mechanics. This paradox is due precisely to the observer-dependence of the cosmological horizon — the sharpest difference from a black hole horizon — although we will argue that the gravity path integral avoids the paradox in controlled examples.

2D gravity↗

The radial action from probe amplitudes to all orders

We extract the relativistic classical radial action from scattering amplitudes, to all orders in perturbation theory, in the probe limit. Our sources include point charges and monopoles, as well as the Schwarzschild and pure-NUT gravitational backgrounds. A characteristic relativistic effect, that scattering trajectories may wind around these sources any number of times, can be recovered when all-order amplitudes are available. We show that the amplitude for scattering a probe off a pure NUT is given by the solution of a transcendental equation involving continued fractions, and explain how to solve this equation to any desired loop order.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Heavy states in 3d gravity and 2d CFT

We discuss correlators of light fields in heavy states in AdS 3 gravity and holographic 2d CFTs. In the bulk, the propagator of free fields in AdS backgrounds containing a conical defect or a BTZ black hole can be obtained by solving a wave equation, as well as by the method of images. On the boundary, these geometries are sourced by heavy operator insertions, and the propagator is dual to a heavy-light (HHLL) correlator. By matching its expansion in Virasoro blocks to our bulk results, we determine the OPE coefficients of all contributing states in both the s and t channels. In the s channel, these states are excitations of the light field on top of the heavy state, and their OPE coefficients are the amplitudes to create them. The t-channel OPE is dominated by the Virasoro vacuum block, but there is also an infinite family of light two-particle states that contribute to the correlator. The OPE coefficients that couple these states to heavy operators represent their expectation values in heavy backgrounds. We determine them exactly, derive their asymptotic form at large twist, and discuss their behavior near and above the BTZ threshold, where they become thermal one-point functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Collisions of localized shocks and quantum circuits

We study collisions between localized shockwaves inside a black hole interior. We give a holographic boundary description of this process in terms of the overlap of two growing perturbations in a shared quantum circuit. The perturbations grow both exponentially as well as ballistically. Due to a competition between different physical effects, the circuit analysis shows dependence on the transverse locations and exhibits four regimes of qualitatively different behaviors. On the gravity side we study properties of the post-collision geometry, using exact calculations in simple setups and estimations in more general circumstances. We show that the circuit analysis offers intuitive and surprisingly accurate predictions about gravity computations involving non-linear features of general relativity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Perturbative unitarity and NEC violation in genesis cosmology

Explorations of the violation of null energy condition (NEC) in cosmology could enrich our understanding of the very early universe and the related gravity theories. Although a fully stable NEC violation can be realized in the “beyond Horndeski” theory, it remains an open question whether a violation of the NEC is allowed by some fundamental properties of UV-complete theories or the consistency requirements of effective field theory (EFT). We investigate the tree-level perturbative unitarity for stable NEC violations in the contexts of both Galileon and “beyond Horndeski” genesis cosmology, in which the universe is asymptotically Minkowskian in the past. We find that the constraints of perturbative unitarity imply that we may need some unknown new physics below the cut-off scale of the EFT other than that represented by the “beyond Horndeski” operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Generalized symmetry in dynamical gravity

We explore generalized symmetry in the context of nonlinear dynamical gravity. Our basic strategy is to transcribe known results from Yang-Mills theory directly to gravity via the tetrad formalism, which recasts general relativity as a gauge theory of the local Lorentz group. By analogy, we deduce that gravity exhibits a one-form symmetry implemented by an operator U α labeled by a center element α of the Lorentz group and associated with a certain area measured in Planck units. The corresponding charged line operator W ρ is the holonomy in a spin representation ρ, which is the gravitational analog of a Wilson loop. The topological linking of U α and W ρ has an elegant physical interpretation from classical gravitation: the former materializes an exotic chiral cosmic string defect whose quantized conical deficit angle is measured by the latter. We verify this claim explicitly in an AdS-Schwarzschild black hole background. Notably, our conclusions imply that the standard model exhibits a new symmetry of nature at scales below the lightest neutrino mass. More generally, the absence of global symmetries in quantum gravity suggests that the gravitational one-form symmetry is either gauged or explicitly broken. The latter mandates the existence of fermions. Finally, we comment on generalizations to magnetic higher-form or higher-group gravitational symmetries.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗