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

Modeling of standing wave RF cavities for tracking through multi-pass energy recovery linac

Short bunches, high current and multiple linac pass are all characteristics of Energy Recovery Linacs (ERLs), which may result in collective effects. They in turn, may affect the beam, degrading its quality, or even yield to instabilities causing a beam loss. To study and mitigate these effects one needs a numerical simulation code, that can take into account both the collective effects, as well as, particular ERL features, such as a multi-turn design that does not reach a steady state or the multiple passages of the beam through Radio-Frequency (RF) cavities at different energies. CODAL [1], a code developed by SOLEIL in collaboration with IJCLab, enables such studies. It is a 6 dimensional (6D) tracking code applying 'kicks' based on the integration of the local Hamiltonian for each element of the lattice. It is also capable of simulating space charge, wakefields and coherent synchrotron radiation. However, to correctly take into account the ERL dynamics, an upgrade had to be made to include the effect of a standing wave RF cavity in 6D. In this paper, we will concentrate on the implementation and benchmarking (with DESY’s tracking code ASTRA [2]) of both the longitudinal and the transverse models (by J.B. Rosenzweig and L. Serafini [3]), which we use to carry out tracking of fully analytical 6D RF cavity.

Bogacz, A.↗

[Presentation Slides] Code-to-Code Benchmark Study for Thermal Stress Modeling and Preliminary Analysis of the High-temperature Single Heat-Pipe Experiment

In the heat-pipe-cooled microreactor, heat pipes remove heat from the reactor core as a passive heat-transfer device, so the fluid circulation is not required for cooling, which can substantially simplify the overall reactor design. However, given the extremely high temperatures in the core region and potentially large temperature gradients across the structure materials, thermal stresses need to be well-analyzed to ensure structural integrity during normal operations and accident scenarios. This presentation slides discuss finite element method-based thermal-stress analysis for the high-temperature single heat-pipe test article in the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory (INL), using two commercial software packages, Abaqus and Star-CCM+. A code-to-code benchmark study was performed to crosscheck the model setup and capability of each code and to gain preliminary insights into the potential thermal stress concerns from the current experimental setup. It is observed the significant thermal stresses happen at the inner surface of the heat pipe hole surrounded by electric heaters where the largest temperature gradients appear. The temperature fields have good agreements between Abaqus and Star-CCM+, while the induced thermal stresses show modest deviations probably due to differences of meshing engines used in these two codes. It is found that the local maximum thermal stresses may reach close to the ultimate tensile strength and yield strength of structural material depending on the heater power. Ultimately, the coupled thermal-structural analysis will help guide the current experimental plan and ensure the facility safety for future experimental study.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Permutation codes for the state assignment of fault tolerant sequential machines

A new fault-tolerant state assignment method is suggested for synchronous sequential machines. It is assumed that the inputs are fault free and that for no input it is possible to reach all or most of the states, whose number may be fairly large. Error correcting codes for the state assignment are generated by permutations of a chosen linear code. A state assignment algorithm is developed and its computational complexity is estimated. Examples are given.

Chen, M.↗

Electron-scale current sheets and energy dissipation in 3D kinetic-scale plasma turbulence with low electron beta

ABSTRACT Three-dimensional kinetic-scale turbulence is studied numerically in the regime where electrons are strongly magnetized (the ratio of plasma species pressure to magnetic pressure is βe = 0.1 for electrons and βi = 1 for ions). Such a regime is relevant in the vicinity of the solar corona, the Earth’s magnetosheath, and other astrophysical systems. The simulations, performed using the fluid-kinetic spectral plasma solver (sps) code, demonstrate that the turbulent cascade in such regimes can reach scales smaller than the electron inertial scale, and results in the formation of electron-scale current sheets (ESCS). Statistical analysis of the geometrical properties of the detected ESCS is performed using an algorithm based on the medial axis transform. A typical half-thickness of the current sheets is found to be on the order of electron inertial length or below, while their half-length falls between the electron and ion inertial length. The pressure–strain interaction, used as a measure of energy dissipation, exhibits high intermittency, with the majority of the total energy exchange occurring in current structures occupying approximately 20 per cent of the total volume. Some of the current sheets corresponding to the largest pressure–strain interaction are found to be associated with Alfvénic electron jets and magnetic configurations typical of reconnection. These reconnection candidates represent about 1 per cent of all the current sheets identified.

79 ASTRONOMY AND ASTROPHYSICS↗

Coronal Magnetic Structures Observing Campaign. II - Magnetic and plasma properties of a solar active region

Simultaneous soft X-ray, microwave, and photospheric magnetic field observations were taken during the Coronal Magnetic Structures Observing Campaign. The plasma electron temperatures and emission measures determined from the X-ray data are used to predict the intensity and structure of the thermal bremsstrahlung emission at 20 and 6 cm. Comparing these predictions with the microwave observations, it is found that the 20 cm structure is very similar to that expected from the X-rays, but a substantial amount of the 6 cm emission was resolved out. The predicted 20 cm brightness temperatures are higher than the observed, requiring cool absorbing material (not greater than 500,000 K) between the hot X-ray plasma and the observer. The absorption mechanism in the cool plasma at 20 cm is most likely thermal bremsstrahlung, requiring coronal magnetic fields as high as 150 G. 'Coronal Magnetograms', made by extrapolating the photospheric longitudinal field using the Sakurai code, show that appropriate values of the total field are reached at heights of 6000-10,000 km above the photosphere (at many but not all locations).

Schmelz, J. T.↗

A Simplified Building Modeling Approach for Identifying Whole Building Retrofits

Building Energy Modeling (BEM) is an effective strategy for optimizing new and retrofit building design, evaluating its energy savings potential and rating its energy performance. However, because of the complexity and expense of modeling, only a small percentage of buildings are simulated. Over 80% of buildings are <25,000 ft2, and the cost associated with BEM can be a big deterrent for reaching this subsector, which usually ends up complying with energy codes and qualifying for incentives using prescriptive measures or generalized design guides. Similarly, efficiency considerations for small building retrofits are often identified through deemed measures or Technical Reference Manuals, which though a scalable and convenient approach, can limit innovation in design and optimization of measures. Simplification of the modeling process, where appropriate, can lead to increased energy savings and more informed decision making. This paper will discuss the development of a ruleset for a simplified approach based on the Performance Rating Method (PRM) contained in ASHRAE Standard 90.1. The approach aims to simplify and lower the cost of the energy modeling process in a controlled and documented manner. An anticipated outcome is an increased use of BEM and its effectiveness in the design, retrofit and operation of commercial buildings. In addition, simplification of the modeling process will result in fewer errors from misinterpretations of both program requirements and program intent in simulation software. The paper will discuss the need and advantages of a simplified modeling approach for whole building retrofits and the technical process for achieving the same.

Goel, Supriya↗

High compression image and image sequence coding

The digital representation of an image requires a very large number of bits. This number is even larger for an image sequence. The goal of image coding is to reduce this number, as much as possible, and reconstruct a faithful duplicate of the original picture or image sequence. Early efforts in image coding, solely guided by information theory, led to a plethora of methods. The compression ratio reached a plateau around 10:1 a couple of years ago. Recent progress in the study of the brain mechanism of vision and scene analysis has opened new vistas in picture coding. Directional sensitivity of the neurones in the visual pathway combined with the separate processing of contours and textures has led to a new class of coding methods capable of achieving compression ratios as high as 100:1 for images and around 300:1 for image sequences. Recent progress on some of the main avenues of object-based methods is presented. These second generation techniques make use of contour-texture modeling, new results in neurophysiology and psychophysics and scene analysis.

Kunt, Murat↗

DFT-FE – A massively parallel adaptive finite-element code for large-scale density functional theory calculations

We present an accurate, efficient and massively parallel finite-element code, DFT-FE, for large-scale ab-initio calculations (reaching ~100,000 electrons) using Kohn-Sham density functional theory (DFT). DFT-FE is based on a local real-space variational formulation of the Kohn-Sham DFT energy functional that is discretized using a higher-order adaptive spectral finite-element (FE) basis, and treats pseudopotential and all-electron calculations in the same framework, while accommodating non-periodic, semi-periodic and periodic boundary conditions. We discuss the main aspects of the code, which include, the various strategies of adaptive FE basis generation, and the different approaches employed in the numerical implementation of the solution of the discrete Kohn-Sham problem that are focused on significantly reducing the floating point operations, communication costs and latency. We demonstrate the accuracy of DFT-FE by comparing the energies, ionic forces and periodic cell stresses on a wide range of problems with popularly used DFT codes. Further, we demonstrate that DFT-FE significantly outperforms widely used plane-wave codes---both in CPU-times and wall-times, and on both non-periodic and periodic systems---at systems sizes beyond a few thousand electrons, with over 5-10 fold speedups in systems with more than 10,000 electrons. In conclusion, the benchmark studies also highlight the excellent parallel scalability of DFT-FE, with strong scaling demonstrated on up to 192,000 MPI tasks.

97 MATHEMATICS AND COMPUTING↗

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↗

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↗

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↗

Evaluation of the Land Surface Reflectance Fundamental Climate Data Record

The land surface reflectance is a fundamental climate data record at the basis of the derivation of other climate data records (Albedo, LAI/Fpar, Vegetation indices) and has been recognized as a key parameter in the understanding of the land-surface-climate processes. In this presentation, we present the validation of the Land surface reflectance used for MODIS, VIIRS, Landsat 8 and Sentinel 2 data. This methodology uses the 6SV Code and data from the AERONET network. The overall accuracy clearly reaches the satellite specifications. To understand how to improve the validation, we developed an exhaustive error budget. Results show an impact of the absorption of aerosol and of the fine mode volume concentration.

Atmospheric correction↗

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↗

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↗

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↗

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↗

Accuracy versus convergence rates for a three dimensional multistage Euler code

Using a central difference scheme, it is necessary to add an artificial viscosity in order to reach a steady state. This viscosity usually consists of a linear fourth difference to eliminate odd-even oscillations and a nonlinear second difference to suppress oscillations in the neighborhood of steep gradients. There are free constants in these differences. As one increases the artificial viscosity, the high modes are dissipated more and the scheme converges more rapidly. However, this higher level of viscosity smooths the shocks and eliminates other features of the flow. Thus, there is a conflict between the requirements of accuracy and efficiency. Examples are presented for a variety of three-dimensional inviscid solutions over isolated wings.

Turkel, Eli↗

Accuracy versus convergence rates for a three dimensional multistage Euler code

Using a central difference scheme, it is necessary to add an artificial viscosity in order to reach a steady state. This viscosity usually consists of a linear fourth difference to eliminate odd-even oscillations and a nonlinear second difference to suppress oscillations in the neighborhood of steep gradients. There are free constants in these differences. As one increases the artificial viscosity, the high modes are dissipated more and the scheme converges more rapidly. However, this higher level of viscosity smooths the shocks and eliminates other features of the flow. Thus, there is a conflict between the requirements of accuracy and efficiency. Examples are presented for a variety of three-dimensional inviscid solutions over isolated wings.

Turkel, Eli↗