Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “D CODES”

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 37 records · Page 2

Computational Modeling and Simulation for Nonproliferation: The History (and Present) of Monte Carlo and the MCNP(R) Code at Los Alamos [Slides]

The emergence of the Monte Carlo method as a research tool springs from work done at Los Alamos in the 1940s.Monte Carlo and the MCNP Code have been and continue to be developed at Los Alamos for many decades. From basic science in support of understanding nuclear interaction physics to global security applications, application uses of the code are extensive. Recent R&D projects and code modernization efforts make the MCNP code a great tool for nuclear nonproliferation applications. In collaboration with nuclear safeguards experts, new training has just recently been developed to help new practitioners learn how to use the code for nuclear safeguards applications.

97 MATHEMATICS AND COMPUTING↗

Error-correcting codes for fermionic quantum simulation

Utilizing the framework of \mathbb{Z}_2 ℤ 2 lattice gauge theories in the context of Pauli stabilizer codes, we present methodologies for simulating fermions via qubit systems on a two-dimensional square lattice. We investigate the symplectic automorphisms of the Pauli module over the Laurent polynomial ring. This enables us to systematically increase the code distances of stabilizer codes while fixing the rate between encoded logical fermions and physical qubits. We identify a family of stabilizer codes suitable for fermion simulation, achieving code distances of d=2,3,4,5,6,7, allowing correction of any \lfloor \frac{d-1}{2} \rfloor ⌊ d − 1 2 ⌋ -qubit error. In contrast to the traditional code concatenation approach, our method can increase the code distances without decreasing the (fermionic) code rate. In particular, we explicitly show all stabilizers and logical operators for codes with code distances of d=3,4,5. We provide syndromes for all Pauli errors and invent a syndrome-matching algorithm to compute code distances numerically.

Physics↗

Prediction of PbLi fluid flow and temperature field in a thermal convection loop for qualification of fusion materials

Two modeling tools, a thermohydraulics code and a computational model in COMSOL Multiphysics, have been developed, tested and then applied to the analysis of fluid flow and heat transfer in a thermal convection loop (TCL). Such a device has recently been used to experimentally evaluate corrosion compatibility of APMT (Fe–21Cr–5Al–3Mo) steel with high-temperature molten eutectic lead-lithium (PbLi) alloy at Oak Ridge National Lab, TN, USA. Here, the 1-D thermohydraulics code allows for rapid calculations of the loop parameters as a function of the applied heat flux. The supplemental COMSOL finite-element computations provide detailed 3-D velocity and temperature field data but are much more time consuming. Both modelling tools demonstrate an excellent agreement in the computed circulation velocity as well as maximum and minimum temperatures. The performed TCL analysis focuses on flow development in the “hot” and “cold” legs, formation of Dean vortices in corner regions, blocking effect of the immersed samples, and radiative and convective cooling effects under the experiment-relevant conditions for Prandtl number Pr = 0.015, Grashof numbers Gr ~10 7 , and Reynolds numbers Re ~10 4 .

36 MATERIALS SCIENCE↗

Low-overhead transversal fault tolerance for universal quantum computation

Fast, reliable logical operations are essential for realizing useful quantum computers. By redundantly encoding logical qubits into many physical qubits and using syndrome measurements to detect and correct errors, we can achieve low logical error rates. However, for many practical quantum error correction codes such as the surface code, owing to syndrome measurement errors, standard constructions require multiple extraction rounds—of the order of the code distance d—for fault-tolerant computation, particularly considering fault-tolerant state preparation. Here we show that logical operations can be performed fault-tolerantly with only a constant number of extraction rounds for a broad class of quantum error correction codes, including the surface code with magic state inputs and feedforward, to achieve ‘transversal algorithmic fault tolerance’. Through the combination of transversal operations7 and new strategies for correlated decoding, despite only having access to partial syndrome information, we prove that the deviation from the ideal logical measurement distribution can be made exponentially small in the distance, even if the instantaneous quantum state cannot be made close to a logical codeword because of measurement errors. We supplement this proof with circuit-level simulations in a range of relevant settings, demonstrating the fault tolerance and competitive performance of our approach. Furthermore, our work sheds new light on the theory of quantum fault tolerance and has the potential to reduce the space–time cost of practical fault-tolerant quantum computation by over an order of magnitude.

Zhou, Hengyun [QuEra Computing, Boston, MA (United↗

Constraining Physical Models at Gigabar Pressures

High-energy-density (HED) experiments in convergent geometry are able to test physical models at pressures beyond hundreds of millions of atmospheres. The measurements from these experiments are generally highly integrated and require unique analysis techniques to procure quantitative information. This work describes a methodology to constrain the physics in convergent HED experiments by adapting the methods common to many other fields of physics. As an example, a mechanical model of an imploding shell is constrained by data from a thin-shelled direct-drive exploding-pusher experiment on the OMEGA Laser System using Bayesian inference, resulting in the reconstruction of the shell dynamics and energy transfer during the implosion. The model is tested by analyzing synthetic data from a 1-D hydrodynamics code and is sampled using a Markov chain Monte Carlo to generate the posterior distributions of the model parameters. The goal of this work is to demonstrate a general methodology that can be used to draw conclusions from a wide variety of HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Constraining physical models at gigabar pressures

High-energy-density (HED) experiments in convergent geometry are able to test physical models at pressures beyond hundreds of millions of atmospheres. The measurements from these experiments are generally highly integrated and require unique analysis techniques to procure quantitative infor mation. This work describes a methodology to constrain the physics in convergent HED experiments by adapting the methods common to many other fields of physics. As an example, a mechanical model of an imploding shell is constrained by data from a thin-shelled direct-drive exploding-pusher experiment on the OMEGA Laser System using Bayesian inference, resulting in the reconstruction of the shell dynamics and energy transfer during the implosion. The model is tested by analyzing synthetic data from a 1-D hydrodynamics code and is sampled using a Markov chain Monte Carlo to generate the posterior distributions of the model parameters. As a result, the goal of this work is to demonstrate a general methodology that can be used to draw conclusions from a wide variety of HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parallel transport sweeps on two-dimensional cartesian and hexagonal grids

This paper aims to provide a proof of concept for parallel transport sweeps on two-dimensional hexagonal grids for the discrete ordinates transport equation. While the method is an extension of the popular and well-established Koch-Baker-Alcoulffe (KBA) algorithm, there are significant differences between the cartesian and hexagonal grid and thereafter sweep. The most important is the three-way connectivity of hexagons within the grid which creates greater dependencies between the elements. The KBA method in structured orthogonal grids was first implemented in the DRAGON5 code and the method is first described here. The differences in implementation for the hexagonal grid are also described. Benchmark results are also presented, showing roughly 10 times speedup in computational times with roughly 100 processors, in both cases. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Scattering of radiofrequency waves by randomly modulated density interfaces in the edge of fusion plasmas

In the scrape-off layer and the edge region of a tokamak, the plasma is strongly turbulent and scatters the radiofrequency (RF) electromagnetic waves that propagate through this region. It is important to know the spectral properties of these scattered RF waves, whether used for diagnostics or for heating and current drive. The spectral changes influence the interpretation of the obtained diagnostic data, and the current and heating profiles. A full-wave, three-dimensional (3-D) electromagnetic code ScaRF has been developed for studying the RF wave propagation through turbulent plasma. ScaRF is a finite-difference frequency-domain (FDFD) method used for solving Maxwell's equations. The magnetized plasma is defined through the cold plasma by the anisotropic permittivity tensor. As a result, ScaRF can be used to study the scattering of any cold plasma RF wave. It can also be used for the study of the scattering of electron cyclotron waves in ITER-type and medium-sized tokamaks such as TCV, ASDEX-U and DIII-D. For the case of medium-sized tokamaks, there is experimental evidence that drift waves and rippling modes are present in the edge region. Hence, we have studied the scattering of RF waves by periodic density interfaces (plasma gratings) in the form of a superposition of spatial modes with varying periodicity and random amplitudes. The power reflection coefficient (a random variable) is calculated for different realizations of the density interface. In this work, the uncertainty of the power reflection coefficient is rigorously quantified by use of the Polynomial Chaos Expansion method in conjunction with the Smolyak sparse-grid integration, which is known as the PCE-SG method. The PCE-SG method is proven to be accurate and more efficient compared with alternative methods such as the Monte Carlo (MC) approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Unabsorbed light beamlets for diagnosing coronal density profiles and absorption nonuniformity in direct-drive implosions on OMEGA

Laser light scattered from a target is the most-direct measurement for diagnosing laser absorption in a direct-drive implosion. The 3ωGOI (gated optical imager) beamlets diagnostic images unabsorbed light from all 60 OMEGA beams as distinct “beamlet” spots for each beam. The implosion can be diagnosed using the position and intensity of these beamlets. The position of each beamlet in the image is determined by refraction and can be used to fit the coronal plasma density profile of the implosion. The inferred plasma density profiles are comparable to profiles predicted by the 1-D hydrodynamics code LILAC but suggest the predictions underestimate the density farther out in the corona. The intensity of light in each spot depends on the cumulative effects of absorption and cross-beam energy transfer along the beamlet’s path through the corona. Furthermore, the measured variation in intensity and polarization between similar spots indicates that absorption during OMEGA implosions is less uniform than previously known.

47 OTHER INSTRUMENTATION↗

Transient analysis of a micro-reactor using the DireWolf code suite

Transient analyses of heat pipe micro-reactors are necessary to ensure that hypothetical accident scenarios do not comprise reactor safety. Due to its small size and reliance on heat-pipes for cooling, the micro-reactor design introduced in this paper is a tightly coupled system which requires multi-physics tools to accurately model transient events. Idaho National Laboratory's DireWolf code suite based on the MOOSE framework is tailor-built to model heat-pipe reactors. This paper demonstrates DireWolf's ability to simulate the coupled thermal-neutronics transient behavior of a heat-pipe micro-reactor. The transient events presented here include an inadvertent rotation of all control drums simultaneously and a sudden complete rotation of a single control drum. A detailed description of each event is provided along with simulation results, including time dependent power and temperature distributions, and discussion. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Logical quantum processor based on reconfigurable atom arrays

Suppressing errors is the central challenge for useful quantum computing, requiring quantum error correction (QEC) for large-scale processing. However, the overhead in the realization of error-corrected ‘logical’ qubits, in which information is encoded across many physical qubits for redundancy, poses substantial challenges to large-scale logical quantum computing. Here we report the realization of a programmable quantum processor based on encoded logical qubits operating with up to 280 physical qubits. Using logical-level control and a zoned architecture in reconfigurable neutral-atom arrays, our system combines high two-qubit gate fidelities, arbitrary connectivity, as well as fully programmable single-qubit rotations and mid-circuit readout. Operating this logical processor with various types of encoding, we demonstrate improvement of a two-qubit logic gate by scaling surface-code distance from d = 3 to d = 7, preparation of colour-code qubits with break-even fidelities, fault-tolerant creation of logical Greenberger–Horne–Zeilinger (GHZ) states and feedforward entanglement teleportation, as well as operation of 40 colour-code qubits. Finally, using 3D [[8,3,2]] code blocks, we realize computationally complex sampling circuits with up to 48 logical qubits entangled with hypercube connectivity with 228 logical two-qubit gates and 48 logical CCZ gates. We find that this logical encoding substantially improves algorithmic performance with error detection, outperforming physical-qubit fidelities at both cross-entropy benchmarking and quantum simulations of fast scrambling. These results herald the advent of early error-corrected quantum computation and chart a path towards large-scale logical processors.

74 ATOMIC AND MOLECULAR PHYSICS↗

ITER ECH Transmission Line System Design and Status

The electron cyclotron (EC) heating & current drive (H&CD) system on ITER provides plasma heating by generating, transmitting, and launching high-intensity, high-frequency (170 GHz) electromagnetic wave energy steerable across the plasma cross-section. The transmission line (TL) subsystem connects the Matching Optics Unit (MOU) on each of the 24 gyrotrons to the 32 feed points in the four upper launchers and the 24 feed points in the equatorial launcher. Each TL must be able to operate at up to 1.2 MW of input power for up to 1 hour pulse lengths. The TL system contains 50 mm water-cooled corrugated waveguide, 90o miter bends, 140o miter bends, polarizer miter bend pairs, switches, expansion units, pumpouts, DC breaks, MOU-TL adapters, Radio Frequency (RF) loads, and isolation shutter valves. A detailed finite element analysis has been used to verify the thermo-mechanical performance of each component. The microwave performance has been analyzed using a 2-D electromagnetic code combined with a Monte Carlo code. This approach allows the impact of manufacturing and installation tolerances to be assessed and optimized to provide a high probability of achieving the system performance requirements. Prototypes of the waveguide and TL components have been fabricated and tested at high-power. Production contracts are now being issued for fabrication and delivery of the waveguide and components to ITER.

Hanson, Gregory↗

Higher-Form Anomalies on Lattices

We show that generic gapped quantum many-body states which respect an anomalous finite higher-form symmetry have an exponentially small overlap with any short-range entangled (SRE) state. Hence, anomalies of higher-form symmetries enforce $intrinsic$ long-range entanglement, which is in contrast with anomalies of ordinary (0-form) symmetries which are compatible with symmetric SRE states (specifically, symmetric cat states). As an application, we show that the anomalies of strong higher-form symmetries provide a diagnostic for mixed-state topological order in $d \geq 2$ spatial dimensions. We also identify a new (3+1)D intrinsic mixed-state topological order that does not obey remote-detectability by local decoherence of the (3+1)D Toric Code with fermionic loop excitations. This breakdown of remote detectability, as encoded in anomalies of strong higher-form symmetries, provides a partial characterization of intrinsically mixed-state topological order.

Feng, Yitao [Peking University, Beijing (China)] (↗

cheny-lanl/PreMevE-MEO-2024

This package contains R&D python codes that can nowcast >2 MeV electron distributions in the Earth's outer radiation belt, using LANL's GPS electron data (publicly available) as the primary driver. Functions of this PreMevE-MEO model has been described with details in a manuscript titled "PreMevE-MEO: Predicting Ultra-relativistic Electrons Using Observations from GPS Satellites" (LA-UR-24-23843), which has been submitted to Space Weather Journal.

Chen, Yue↗

LANL contributions to Rayleigh

Rayleigh is a 3-D convection code designed for the study of dynamo behavior in spherical geometry. It evolves the incompressible and anelastic MHD equations in spherical geometry using a pseudo-spectral approach.

Edelmann, Philipp↗

Noisy quantum trees: infinite protection without correction

We study quantum networks with tree structures, in which information propagates from a root to leaves. At each node in the network, the received qubit unitarily interacts with fresh ancilla qubits, after which each qubit is sent through a noisy channel to a different node in the next level. Therefore, as the tree depth grows, there is a competition between the irreversible effect of noise and the protection against such noise achieved by the delocalization of information. In the classical setting, where each node simply copies the input bit into multiple output bits, this model has been studied as the broadcasting or reconstruction problem on trees, which has broad applications. In this work, we study the quantum version of this problem. We consider a Clifford encoder at each node that encodes the input qubit in a stabilizer code, along with a single qubit Pauli noise channel at each edge. Such noisy quantum trees describe a scenario in which one has access to a stream of fresh (low-entropy) ancilla qubits, but cannot perform error correction. Therefore, they provide a different perspective on quantum fault tolerance. Furthermore, they provide a useful model for describing the effect of noise within the encoders of concatenated codes. We prove that above certain noise thresholds, which depend on the properties of the code such as its distance, as well as the properties of the encoder, information decays exponentially with the depth of the tree. On the other hand, by studying certain efficient decoders, we prove that for codes with distance d ≥ 2 and for sufficiently small (but non-zero) noise, classical information and entanglement propagate over a noisy tree with infinite depth. Indeed, we find that this remains true even for binary trees with certain 2-qubit encoders at each node, which encodes the received qubit in the binary repetition code with distance d = 1.

Quantum information↗

Quantitative proton radiography and shadowgraphy for arbitrary intensities

Charged-particle radiography and shadowgraphy data can be directly inverted to obtain a line-integrated transverse Lorentz force or a line-integrated transverse refractive index gradient if intensity modulations due to scattering and absorption are negligible, and angular deflections are small. We develop a new direct-inversion algorithm based on plasma physics and compare it to a new Monge–Ampère code and an existing power diagram code. The measured or source intensity is represented by electrons subject to drag, and the other intensity by fixed ions. The decrease in kinetic plus electrostatic energy determines convergence. The displacement of the electrons from their initial to their equilibrium positions determines the line-integrated force or refractive index gradient. We have implemented two approaches: PIC (particle in cell) and Lagrangian fluid, in 1-D and 2-D. The PIC code works for arbitrary intensities, can work efficiently in parallel, and can make use of existing codes. The Lagrangian code requires less memory and is faster than the PIC code without massively parallel processing, but fails in 2-D for large intensity modulations. The Monge–Ampère code is by far the fastest in 2-D, without massively parallel processing, but fails for intensities with large voids, high contrast ratios and large deflections across the boundaries, and could not obtain the degree of convergence possible with the PIC code. As a result, the power diagram code was by far the slowest and most memory intensive, and failed for large peaks in the measured intensity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE ART Graphite Program Intro Status 2023

General description of the ART Graphite Program with Status for 2023. This includes list of staff members, and brief description of behavior models, licensing & code, graphite R&D, As-Fab'd Properties, Machanism and Analysis, Irradiation. Also detail FY23 Activites, with AGC Update, Oxidation Activiteis, Baseline STatus, ASME Code Development, VIC Project Update, Molten salt studies. Slide included to highlight each employee contributions for work and papers and NEUPs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗