Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reach code”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Fault-Tolerant Operation of Bosonic Qubits with Discrete-Variable Ancillae

Fault-tolerant quantum computation with bosonic qubits often necessitates the use of noisy discrete-variable ancillae. In this work, we establish a comprehensive and practical fault-tolerance framework for such a hybrid system and synthesize it with fault-tolerant protocols by combining bosonic quantum error correction (QEC) and advanced quantum control techniques. We introduce essential building blocks of error-corrected gadgets by leveraging ancilla-assisted bosonic operations using a generalized variant of path-independent quantum control. Using these building blocks, we construct a universal set of error-corrected gadgets that tolerate a single-photon loss and an arbitrary ancilla fault for four-legged cat qubits. Notably, our construction requires only dispersive coupling between bosonic modes and ancillae, as well as beam-splitter coupling between bosonic modes, both of which have been experimentally demonstrated with strong strengths and high accuracy. Moreover, each error-corrected bosonic qubit is comprised of only a single bosonic mode and a three-level ancilla, featuring the hardware efficiency of bosonic QEC in the full fault-tolerant setting. We numerically demonstrate the feasibility of our schemes using current experimental parameters in the circuit-QED platform. Finally, we present a hardware-efficient architecture for fault-tolerant quantum computing by concatenating the four-legged cat qubits with an outer qubit code utilizing only beam-splitter couplings. Our estimates suggest that the overall noise threshold can be reached using existing hardware. These developed fault-tolerant schemes extend beyond their applicability to four-legged cat qubits and can be adapted for other rotation-symmetrical codes, offering a promising avenue toward scalable and robust quantum computation with bosonic qubits. Published by the American Physical Society 2024

Physics↗
DOE OSTI2014

Modeling and simulation of hydrogen behavior in Zircaloy-4 fuel cladding

As a result of corrosion during normal operation in nuclear reactors, hydrogen can enter the zirconium-alloy fuel cladding and precipitate as brittle hydride platelets, which can severely degrade the cladding ductility. Under a heterogeneous temperature distribution, hydrides tend to accumulate in the colder areas, creating local spots of degraded cladding that can favor crack initiation. Therefore, an estimation of the local hydride distribution is necessary to help predict the risk of cladding failure. The hydride distribution is governed by three competing phenomena. Hydrogen in solid solution diffuses under a concentration gradient due to Fick’s law and under a temperature gradient due to the Soret effect. Precipitation of the hydride platelets occurs once the hydrogen solubility limit is reached. A model of these phenomena was implemented in the 3D fuel performance code BISON in order to calculate the hydrogen distribution for arbitrary geometries, such as a nuclear fuel rod, and is now available for BISON users. Simulations have been performed on simple geometries to validate the model and its implementation. The simulations predict that before precipitation occurs, hydrogen tends to accumulate in the colder spots due to the Soret effect. Once the solubility limit is reached, hydrogen precipitates and forms a rim close to the outer edge of the cladding. The simulations also predict that the reactor shut down has little effect on already precipitated hydrides but causes the remaining hydrogen to precipitate homogeneously into hydrides.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Comparing Measured (a,n) Thick Target Yields to MCNP6 using TENDL 2012 Libraries - Paper 31

MCNP6, due to the heritage of MCNPX, has always possessed the capability to transport the light ions, such as alpha particles. Two types of competing physics take place during the particle transport: atomic and nuclear physics. Charged particles lose energy to media due to Coulombic interactions. MCNP6 models Coulombic interactions through the use of the continuous slowing down approximation. Stopping powers are calculated at a variety of energies using the Bethe-Bloch formula with Sternheimer and Peierls density effect correction. For energies less than 1.31 MeV/AMU, a low-energy SPAR model is used. When the energy of the particle is high enough, there is a transition from nuclear interactions to nucleon interactions. MCNP6 possess various physics models to account for these types of intra-nuclear cascades. These models are less accurate at low energies. Even though MCNP6 will transport these light ions to low energies (i.e. less than 20 MeV), most of the available model physics in MCNP6, for nuclear interactions, has a suggested minimum energy of ∼20-200 MeV. The result of the poorer prediction at low energies causes important light ion events, such as F-19(a,n)Na-22 from U-234 decay in UF{sub 6}, to produce orders of magnitude fewer neutrons in MCNP6 than should be emitted. As with neutrons, ACE (A Compact ENDF) formatted light ion libraries have been developed from several different evaluations. These libraries possess better physics at low energies where the physics models break down. The TENDL 2012 evaluations possess alpha transport libraries, in ACE format, such that they can easily be used with codes like MCNP6 that read ACE formatted cross sections. As a result, a capability was developed for MCNP6 to read and implement these library cross sections when the light ion reached energies below the library upper energy cutoff. Perfectly modelling an experiment for benchmark comparisons to a code calculation can be tricky. For complex setups, systematic error is always possible; and usually a hefty uncertainty analysis is required to quantify just the geometry setup and material combination choices. For an (a,n) interaction, a thick target yield is defined as the amount of neutrons emitted given that the alpha ranges out in the material. For alpha transport at low energy <10 MeV, in reasonably dense materials (>1-2 g/cc), with a material dimension length larger than a few hundred microns, alpha particles will tend to completely stop and range out in the material. For Li, Be, B, C, N, O, F, Na, Mg and Al thick target yields exist in 200 keV intervals from 1 to 9.8 MeV (most of the elements do have (a,n) thresholds above 1 MeV). Because modelling a thick target yield in MCNP6 is much easier than attempting to account for the complexity of systematic error from a full out experiment, we present here comparisons of (a,n) thick target yields to those calculated in MCNP6 using the TENDL 2012 alpha library evaluations.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

BEYONDPLANCK VIII. Efficient sidelobe convolution and corrections through spin harmonics

We introduce a new formulation of the Conviqt convolution algorithm in terms of spin harmonics, and apply this to the problem of sidelobe correction for BEYONDPLANCK, the first end-to-end Bayesian Gibbs sampling framework for CMB analysis. We compare our implementation to the previous Planck LevelS implementation, and find good agreement between the two codes in terms of accuracy, but with a speed-up reaching a factor of 3–10, depending on the frequency bandlimits, l max and m max . The new algorithm is significantly simpler to implement and maintain, since all low-level calculations are handled through an external spherical harmonic transform library. We find that our mean sidelobe estimates for Planck LFI are in good agreement with previous efforts. Additionally, we present novel sidelobe rms maps that quantify the uncertainty in the sidelobe corrections due to variations in the sky model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Nonlocal effects on thermal transport in hydrodynamic simulations of unmagnetized MagLIF-relevant gaspipes on NIF

We present simulations of heat flow relevant to gaspipe experiments on the National Ignition Facility to investigate kinetic effects on transport phenomena. D 2 and neopentane (C 5 H 12 ) filled targets are used to study the laser preheat stage of a MagLIF scheme where an axial magnetic field is sometimes applied to the target. Simulations were done with the radiation-MHD code HYDRA with a collision-dominated fluid model and the SNB nonlocal electron thermal conduction model. Using the SNB model to evolve the electron temperature increased the heat front propagation of neopentane gas targets compared to a local model by limiting radial heat flow. This increases electron temperature near the axis, which decreases laser absorption. We find that the effect of heat flow models on temperature profiles and laser propagation is modest. Beyond the SNB model, we utilize HYDRA to initialize plasma conditions for the Vlasov–Fokker–Planck K2 code. We run K2 until a quasi-steady state is reached and examine the impact of kinetic effects on heat transport. Although axial heat flow is well predicted by fluid models, the fluid model consistently overpredicts radial heat flow up to 150% in regions with the largest temperature gradient of D 2 filled gaspipes. On the other hand, the SNB nonlocal electron conduction model is found to be adequate for capturing kinetic heat flow in gaspipes.

Lau, Ryan Y. [Univ. of Colorado, Boulder, CO (Unit↗
DOE OSTI2020

Counterweight for Long Reach Robotic Arm

Maintenance for vital systems proves challenging due to their location: tucked up and behind components. Pinpointing the precise location of errors presents difficulty when relying solely on the human eye. The solution to this issue arises in the form of a Long Reach Robotic Arm. An adequately balanced robot with several individually controlled links contains the capability to maneuver around obstructing components to locate the points of failure with accuracy. The design process for such a robot begins with the writing of a program. The program intends to visualize how the Long Reach Robotic Arm will move in 2D space by printing a select number of lines representing each arm link to the screen and controlling each link through the keyboard. This code substantially simplifies the design process for the multi-linked arm. The next step in the Long Reach Robotic Arm design process is the development of the arm and its counterweight system. The clear division of this work allows two interns to work on this project: one assumes responsibility for the arm, and the other for the counterweight. I undertook the design for the counterbalance. My half of the project is essential for maintaining safety as an extended multi-linked robot arm presents the challenge of preserving stability. The best method to supply a balance to the robot, even in extension, is a mobile counterweight with a compact design.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

2D radial-azimuthal particle-in-cell benchmark for E × B discharges

In this paper we propose a representative simulation test-case of E × B discharges accounting for plasma wall interactions with the presence of both the electron cyclotron drift instability and the modified-two-stream-instability. Seven independently developed particle-in-cell (PIC) codes have simulated this benchmark case, with the same specified conditions. The characteristics of the different codes and computing times are given. Here, results show that both instabilities were captured in a similar fashion and good agreement between the different PIC codes is reported as main plasma parameters were closely related within a 5% interval. The number of macroparticles per cell was also varied and statistical convergence was reached. Detailed outputs are given in the supplementary data, to be used by other similar groups in the perspective of code verification.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

The staged z-pinch as a potential high gain fusion energy source: An independent review, a negative conclusion

Previously published calculations predict that the “staged z-pinch” (SZP) can achieve 400 MJ of fusion yield on a Z-class machine. The SZP is touted to need no external preheat mechanism and no external pre-magnetization method. Instead, it is claimed that the imploding liner can adequately “shock preheat” the fuel and magnetic field diffusion through the liner can adequately magnetize the fuel. Here in this paper, we analyze a number of published SZP calculations and demonstrate that the calculations have major errors—the computer code used to do the calculations does not appear to be accurately solving the physical model it is intended to solve. A variety of independent analyses lead to this conclusion. This conclusion is confirmed by detailed one-dimensional magnetohydrodynamic (MHD) calculations conducted on different computer codes using a variety of proposed SZP operating parameters. Although using parameters similar or identical to the published calculations, our MHD calculations do not reach fusion conditions; there is no conceivable modification of the parameters that would lead to high-gain fusion conditions using these other codes. Our analyses and a review of the magnetized target parameter space leads to further conclusion that the SZP should not be considered to be a potential high-gain fusion source.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Long Reach Robotic Arm

The Accelerator Division at Fermilab began a robotics program, known as the AD Robotics Initiative. A multi-axis robotic arm was proposed to be included in the program. This summer internship, our team began the planning and design for this project, called the Long Reach Robotic Arm. The problem is that critical systems around the magnets require analysis when water leaks cause them to fail. However, the tightly compacted parts and the quadrupole magnet, weighing approximately 1,000 pounds, cause access viewing problems. The process of having to dismantle and remove the magnet is time-consuming; it is not an efficient solution. The goal for the Long Reach Robotic Arm was to extend and reach behind the component of the Main Injector while "snaking" around and avoiding obstructions. At the start of the project, my teammate, Amanda Hoeksema, and I coded a program in the C language that will determine the number of links and their respective lengths for the robot arm. After, we began the 3D CAD model of the robot arm, base, and counterweight system. The design of the Long Reach Robotic Arm provided a learning experience in learning how to code in the C programming language, using SDL libraries and understanding mechanical engineering.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

Addressing Over-Solving in Multiphysics Codes with Nested Solvers

Examples of multiphysics problems in the field of nuclear engineering are numerous. For light water reactors the coupling of neutronics with thermal-hydraulics is the most predominant. Often this coupling is modeled by combining existing codes in a Picard Iteration scheme, that is, by solving one physics then using that answer to solve the other and iterating until both answers stabilize. However, simply solving to a high precision on every Picard (multiphysics) iteration is much less efficient than progressively solving to higher precisions. This unneeded work is known as over-solving. Over-solving in single-physics solvers has been studied previously for the Newton method and Chebyshev iteration. Both of these methods solve an inner linear problem that determines the next update step for the outer loop. If the inner problem is solved to a high precision in every outer iteration, the result will be substantially more computational expense (over-solving) and a higher probability of reaching the 'wrong' solution. Methods for reducing over-solving in partitioned multiphysics codes have only been developed recently. This work introduces further improvements, including an investigation of over-solving methods on both the multiphysics and single-physics levels.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

The thermal state of KS 1731−260 after 14.5 years in quiescence

Crustal cooling of accretion-heated neutron stars provides insight into the stellar interior of neutron stars. The neutron star X-ray transient, KS 1731−260, was in outburst for 12.5 years before returning to quiescence in 2001. We have monitored the cooling of this source since then through Chandra and XMM-Newton observations. Here we present a 150 ks Chandra observation of KS 1731−260 taken in 2015 August, about 14.5 years into quiescence and 6 years after the previous observation. We find that the neutron star surface temperature is consistent with the previous observation, suggesting that crustal cooling has likely stopped and the crust has reached thermal equilibrium with the core. Using a theoretical crust thermal evolution code, we fit the observed cooling curves and constrain the core temperature (T {sub c} = 9.35 ± 0.25 × 10{sup 7} K), composition (Q {sub imp}=4.4{sub −0.5}{sup +2.2}), and level of extra shallow heating required (Q {sub sh} = 1.36 ± 0.18 MeV/nucleon). We find that the presence of a low thermal conductivity layer, as expected from nuclear pasta, is not required to fit the cooling curve well, but cannot be excluded either.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

FleCSPH notes

FleCSI is a compile-time con gurable framework designed to support multi-physics application development. As such, FleCSI provides a very general set of infrastructure design patterns that can be specialized and extended to suit the needs of a broad variety of solver and data requirements. FleCSI currently supports multi-dimensional mesh topology, geometry, and adjacency information, as well as n-dimensional hashed-tree data structures, graph partitioning interfaces, and dependency closures. FleCSI introduces a functional programming model with control, execution, and data abstractions that are consistent both with MPI and with state-of-the-art, task-based runtimes such as Legion and Charm++. The abstraction layer insulates developers from the underlying runtime, while allowing support for multiple runtime systems including conventional models like asynchronous MPI. The intent is to provide developers with a concrete set of user-friendly programming tools that can be used now, while allowing exibility in choosing runtime implementations and optimization that can be applied to future architectures and runtimes. FleCSI's control and execution models provide formal nomenclature for describing poorly understood concepts such as kernels and tasks. FleCSI's data model provides a low-buy-in approach that makes it an attractive option for many application projects, as developers are not locked into particular layouts or data structure representations. FleCSI currently provides a parallel but not distributed implementation of Binary, Quad and Oct-tree topology. This implementation is base on space lling curves domain decomposition, the Morton order. The current FleCSI version requires the implementation of a driver and a specialization driver. The role of the specialization driver is to provide the data distribution. This feature is not complete in FleCSI code and we provide it. The next step will be to incorporate it directly from FleCSPH to FleCSI as we reach a good level of performance. Then the driver represent the general execution of the resolution without worrying of the data locality and communications. As FleCSI is an On-Development code the structure may change in the future and we keep track of these changes in FleCSPH.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

CONVECTIVE PROPERTIES OF ROTATING TWO-DIMENSIONAL CORE-COLLAPSE SUPERNOVA PROGENITORS

We explore the effects of rotation on convective carbon, oxygen, and silicon shell burning during the late stages of evolution in a 20 M {sub ⊙} star. Using the Modules for Experiments in Stellar Astrophysics we construct one-dimensional (1D) stellar models both with no rotation and with an initial rigid rotation of 50% of critical. At different points during the evolution, we map the 1D models into 2D and follow the multidimensional evolution using the FLASH compressible hydrodynamics code for many convective turnover times until a quasi-steady state is reached. We characterize the strength and scale of convective motions via decomposition of the momentum density into vector spherical harmonics. We find that rotation influences the total power in solenoidal modes, with a slightly larger impact for carbon and oxygen shell burning than for silicon shell burning. Including rotation in 1D stellar evolution models alters the structure of the star in a manner that has a significant impact on the character of multidimensional convection. Adding modest amounts of rotation to a stellar model that ignores rotation during the evolutionary stage, however, has little impact on the character of the resulting convection. Since the spatial scale and strength of convection present at the point of core collapse directly influence the supernova mechanism, our results suggest that rotation could play an important role in setting the stage for massive stellar explosions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Worst Case Sabotage Consequence Modeling for Advanced Reactors

Under proposed rulemaking by the United States Nuclear Regulatory Commission, advanced reactor licensees may qualify for a reduction in onsite responders if the consequences from a sabotage attack do not exceed 25 rem at the reactor site boundary 2 hours after a release occurs. This work presents a preliminary assessment of the feasibility of meeting the 25 rem at the site boundary criteria. Using advanced reactor radionuclide inventories radiological release estimates created by collaborators at Oak Ridge National Laboratory, the MACCS consequence analysis code was used to investigate the distance at the dose from hypothetical sabotage events reaches 25 rem for different advanced reactor technologies. It was found that the distance to 25 rem depended heavily on sampled weather conditions and sabotage release fractions, and further refinement of release fractions and MACCS modeling parameters is recommended.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

MULTI-PARAMETRIC STUDY OF RISING 3D BUOYANT FLUX TUBES IN AN ADIABATIC STRATIFICATION USING AMR

We study the buoyant rise of magnetic flux tubes embedded in an adiabatic stratification using two-and three-dimensional, magnetohydrodynamic simulations. We analyze the dependence of the tube evolution on the field line twist and on the curvature of the tube axis in different diffusion regimes. To be able to achieve a comparatively high spatial resolution we use the FLASH code, which has a built-in Adaptive Mesh Refinement (AMR) capability. Our 3D experiments reach Reynolds numbers that permit a reasonable comparison of the results with those of previous 2D simulations. When the experiments are run without AMR, hence with a comparatively large diffusivity, the amount of longitudinal magnetic flux retained inside the tube increases with the curvature of the tube axis. However, when a low-diffusion regime is reached by using the AMR algorithms, the magnetic twist is able to prevent the splitting of the magnetic loop into vortex tubes and the loop curvature does not play any significant role. We detect the generation of vorticity in the main body of the tube of opposite sign on the opposite sides of the apex. This is a consequence of the inhomogeneity of the azimuthal component of the field on the flux surfaces. The lift force associated with this global vorticity makes the flanks of the tube move away from their initial vertical plane in an antisymmetric fashion. The trajectories have an oscillatory motion superimposed, due to the shedding of vortex rolls to the wake, which creates a Von Karman street.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Integration of Solar Training into Allied Professional Development Platforms - Final Technical Report

The Interstate Renewable Energy Council (IREC) designed the “Integration of Solar Training into Allied Professional Development Platforms” project to surmount barriers to the expansion of solar technology caused by gaps in knowledge impacting on-the-job performance of allied professionals. At project start, questionnaires deployed nationally indicated that only 9% of firefighters felt certain they could shut down a PV system to the extent possible, and only a third of code officials thought they could even discern the difference between a PV and solar thermal system. Therefore, the project objective was to measurably and significantly increase the number of code officials and fire service professionals with solar knowledge directly applicable to their job function. IREC identified specific gaps in current knowledge, and integrated research-proven adult learning principles and technology solutions to enable a cost-effect but impactful reach. The project design ensured that solar training is a continuing component of the targeted sector professional development toolkit. To this end, the key partners recruited and cultivated were existing, trusted training partners in the targeted sectors: International Code Council Solar Rating and Certification Corporation (ICC-SRCC), International Association of Electrical Inspectors (IAEI), International Association of Fire Fighters (IAFF) and National Association of State Fire Marshals (NASFM). During the project, IREC held 79 courses in 25 states that trained 4,156 people, offered a PV in the 2017 National Electric Code webinar attended by 787 people, developed on-line training for code officials with 2,298 on-line learners, and developed an on-line training for firefighters with 4,930 on-line learners. In addition, the project participated in 28 conferences and reached over 537,000 professionals with solar information.

14 SOLAR ENERGY↗
DOE OSTI2023

Design and construction of the MUSE permanent magnet stellarator

This paper documents the design and construction of MUSE, the world's first permanent magnet (PM) stellarator and the first quasi-axisymmetric experiment. The purpose of MUSE is to develop and assess a new way of building optimised stellarators that uses simple planar coils PMs. Our PM optimisation algorithm consists of initialising a geometry to pack dipoles densely, running the FAMUS code to minimise surface field error subject to PM constraints and applying discrete jumps to reach a physically realisable solution. FAMUS treats the PM system as a set of ideal point dipoles. From there we construct finite-volume magnet towers to be housed in 3D-printed PM holders. We describe the design of the PM holders, which were validated by laser metrology. We analyse the effects of finite permeability, sensitivity to perturbations and magnetostatic forces. An exact analytic formula for the magnetic field from a finite-volume PM tower is presented to compute PM–PM forces and stress on the PM holder. Stellarator construction is complete and experiments are underway.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Computational design of short pulse laser driven iron opacity experiments

Here, the resolution of current disagreements between solar parameters calculated from models and observations would benefit from the experimental validation of theoretical opacity models. Iron's complex ionic structure and large contribution to the opacity in the radiative zone of the sun make iron a good candidate for validation. Short pulse lasers can be used to heat buried layer targets to plasma conditions comparable to the radiative zone of the sun, and the frequency dependent opacity can be inferred from the target's measured x-ray emission. Target and laser parameters must be optimized to reach specific plasma conditions and meet x-ray emission requirements. The HYDRA radiation hydrodynamics code is used to investigate the effects of modifying laser irradiance and target dimensions on the plasma conditions, x-ray emission, and inferred opacity of iron and iron-magnesium buried layer targets. It was determined that plasma conditions are dominantly controlled by the laser energy and the tamper thickness. The accuracy of the inferred opacity is sensitive to tamper emission and optical depth effects. Experiments at conditions relevant to the radiative zone of the sun would investigate the validity of opacity theories important to resolving disagreements between solar parameters calculated from models and observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗