SEARCH · Search NASA
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.
It is a serial/parallel code for fast electron beam longitudinal forward or backward tracking through an electron linear accelerator. This codes includes longitudinal space-charge effect, structure and resistive wall wakefields, and longitudinal coherent synchrotron radiation effect. To use it as a parallel code, please comment out the "use mpistub" in BeamBunch.f90, rename mpif.h as mpif.hh, and use the Makefile_parallel as Makefile.
LW3D is code for modeling beam dynamics in particle accelerators including Lienard-Wiechert self-fields. It is similar to codes like IMPACT-T except that the Poisson solver is replaced by a Lienard-Wiechert solver. Along with space-charge effects, the code can model effects like coherent synchrotron radiation.
As higher intensity accelerators are designed, built and operated, the resonances and instabilities caused by space charge can become dominant effects in limiting the intensity of the beam that can be delivered. Particle loss has many causes but space charge can be a major factor in reducing beam intensity and generating emittance growth, particularly at low energies and in the lower energy stages of accelerators such as the AGS Booster. Successful operation relies on a good understanding of basic particle dynamics and it becomes an important challenge to ascertain how the beam evolves when charged particles experience both external forces through the fields of the magnets combined with internal forces through Coulomb interactions. A possible approach to such a study is to include calculations of space charge in computer simulation codes that track the motions of charged particles in 6D phase space. Unfortunately, repeated space charge calculations in 6D phase space are computationally intensive even with the most modern computing facilities. Study with traditional models of space charge is then either limited to a short real-time scale, or sacrifices accuracy and often omits some of the physics. The study described in this report attempts to explore all aspects of the Booster from injection, through RF capture and acceleration to the extraction energy. The simulations are allied to real machine settings. They do not attempt to model effects that are well known and for which accurate experimental data are available, but they try in many cases to extend beyond the normal Booster parameter range so as to understand the limitations on operations and suggest areas for improvement.
Nuclear implementation of the density functional theory (DFT) is at present the only microscopic framework applicable to the whole nuclear landscape. The extension of DFT to superfluid systems in the spirit of the Kohn-Sham approach, the superfluid local density approximation (SLDA) and its extension to time-dependent situations, time-dependent superfluid local density ap- proximation (TDSLDA), have been extensively used to describe various static and dynamical problems in nuclear physics, neutron star crust, and cold atom systems. In this paper, we present the codes that solve the static and time-dependent SLDA equations in three-dimensional coordinate space without any symmetry restriction. These codes are fully parallelized with the message passing interface (MPI) library and take advantage of graphic processing units (GPU) for accelerating execution. The dynamic codes have checkpoint/restart capabilities and for initial conditions one can use any generalized Slater determinant type of wave function. The code can describe a large number of physical problems: nuclear fission, collisions of heavy ions, the interaction of quantized vor- tices with nuclei in the nuclear star crust, excitation of superfluid fermion systems by time dependent external fields, quantum shock waves, domain wall generation and propagation, the dynamics of the Anderson-Bogoliubov-Higgs mode, dynamics of fragmented condensates, vortex rings dynamics, generation and dynamics of quantized vortices, their crossing and recombinations and the incipient phases of quantum turbulence.
Prompt loss of beam injected fast ions approaches 100% in lithium tokamak experiment-beta (LTX-β) discharges, though significantly improved confinement is expected for the higher current plasmas made available by a recent upgrade to the Ohmic heating power supply. Modeling of fast ions using TRANSP/NUBEAM finds a maximum coupled beam fraction of 76% at the near-term limits of the LTX-β operating space. The full ion orbit code POET is employed to validate NUBEAM results against possible non-adiabatic effects on fast ion orbits, but corrections to the prompt loss fraction due to collisionless transport are found to be small. The graphical method code CONBEAM is used to investigate the topology of fast ion phase space as it relates to neutral beam deposition, and counter-injected NBI is considered as a way to access a region of high field side beam deposition. A metric is developed within the CONBEAM using a beam filament model to estimate the prompt loss fraction and shown to agree well with both POET and NUBEAM, enabling near real-time analysis and potential feedback to operators between plasma discharges.
The Proton Improvement Plan phase II (PIP-II) project currently under construction at FNAL will replace the existing 400 MeV normal conducting linac with a new 800 MeV superconducting linac. The beam power in the downstream rapid-cycling Booster synchrotron will be doubled by raising the machine cycle frequency from 15 to 20 Hz and by increasing the injected beam intensity by a factor 1.5. This has to be accomplished without raising uncontrolled losses beyond the administrative limit of 500 W. In addition, slip-stacking efficiency in the Recycler, the next machine in the accelerator chain, sets an upper limit on the longitudinal emittance of the beam delivered by the Booster. As part of an effort to better understand potential losses and emittance blow-up in the Booster, we have been conducting full cycle 6D simulations using the code PyORBIT. The simulations include space charge, wall impedance effects and transition crossing. In this paper, we discuss our experience with the code and present representative results for possible operational scenarios.
The Proton Improvement Plan phase II (PIP-II) project currently under construction at FNAL will replace the existing 400 MeV normal conducting linac with a new 800 MeV superconducting linac. The beam power in the downstream rapid-cycling Booster synchrotron will be doubled by raising the machine cycle frequency from 15 to 20 Hz and by increasing the injected beam intensity by a factor 1.5. This has to be accomplished without raising uncontrolled losses beyond the administrative limit of 500 W. In addition, slip-stacking efficiency in the Recycler, the next machine in the accelerator chain, sets an upper limit on the longitudinal emittance of the beam delivered by the Booster. As part of an effort to better understand potential losses and emittance blow-up in the Booster, we have been conducting full cycle 6D simulations using the code PyORBIT. The simulations include space charge, wall impedance effects and transition crossing. In this paper, we discuss our experience with the code and present representative results for possible operational scenarios.
Optical constants characterize the interaction of materials with light and are important properties in material design. Here we present a Python-based Corvus workflow for simulations of full spectrum optical constants from the visible and ultraviolet to hard x-ray wavelengths based on the real-space Green’s function code FEFF10 and structural data from the Materials Project (MP). The Corvus workflow manager and its associated tools provide an interface to FEFF10 and the MP database. The workflow parallelizes the FEFF computations of optical constants over all absorption edges for each material in the MP database specified by a unique MP-ID. The workflow tools determine the distribution of computational resources needed for that case. Similarly, the optical constants for selected sets of materials can be computed in a single-shot. Additionally, to illustrate the approach, we present results for several elemental solids in the periodic table, as well as a sample compound, and compare our predictions with experimental results. In addition, we provide a database of calculated results for all elements for which there is a stable elemental solid at standard conditions available in the Materials Project database. As in x-ray absorption spectra, these results are interpreted in terms of an atomic-like background and fine-structure contributions.
Here, a series of Yb-substituted Zintl phases in the Ca 3–x Yb x AlSb 3 (0 ≤ x ≤ 0.81(1)) system has been synthesized by initial arc melting and post-heat treatment, and their isotypic crystal structures were characterized by both powder and single crystal X-ray diffraction analysis. All four title compounds adopted the Ca 3 AlAs 3 -type structure (space group Pnma, Pearson code oP28, Z = 4). The overall structure can be described as a combination of the 1-dimensional (1D) infinite chain of ∞ 1 [Al(Sb 2 Sb 2/2 )] formed by two vertices sharing [AlSb 4 ] tetrahedral moieties and three Ca 2+ /Yb 2+ mixed sites located in between these 1D chains. The charge balance and the resultant independency of the 1D chains in the title system were explained by the Zintl-Klemm formalism [Ca 2+ /Yb 2+ ] 3 [(4b-Al 1– )(1b-Sb 2– ) 2 (2b-Sb 1– ) 2/2 ]. A series of DFT calculations proved that (1) the band overlap between the d-orbital states from two types of cations and the p-orbital states from Sb at the high symmetry Γ point implied a heavily doped degenerate semiconducting behavior of the quaternary Ca 2 YbAlSb 3 model and (2) the site preference of Yb for the M1 site was due to the electronic-factor criterion based on the Q values of each atomic site. The electron localization function calculations also proved that the two different shapes of lone pairs of the Sb atoms—the “umbrella-shape” and the “C-shape”—are determined by local geometry and the coordination environment on the anionic frameworks. Thermoelectric measurements of the quaternary title compound Ca 2.19(1) Yb 0.81 AlSb 3 showed an approximately two times larger ZT value than that of ternary Ca 3 AlSb 3 at 623 K due to increased electrical conductivity and ultralow thermal conductivity originated from Yb substitution for Ca.
The aim of the study is to model and characterize the soft X-ray emissivity on the Earth magnetosphere for different space weather conditions (SWC), providing information to interpret the soft X-ray measurements of the Solar wind Magnetosphere Ionosphere Link Explorer space mission. The MHD code pluto in spherical coordinates is used to perform parametric studies with respect to the solar wind (SW) dynamic pressure (considering density and velocity effects independently) as well as the IMF intensity and orientation, predicting the soft X-ray emissivity for different SWC. The integrated soft X-ray emissivity inside the magnetosheath is calculated as a proxy of the soft X-ray emission dependencies with the SWC independently of the satellite orbit and camera line of sight. The analysis indicates fluctuations of the interplanetary magnetic field (IMF) orientation and magnitude may significantly affect the measured soft X-ray emission although changes in the SW dynamic pressure should be the main source of variability. The southward IMF orientation leads to the configuration with the largest soft X-ray emissivity and northward to the lowest. Strongly distorted magnetospheres explored in configurations showing SW and IMF parameters comparable to the impact of interplanetary coronal mass ejections may show a decrease of the soft X-ray emissivity as the IMF magnitude increases, explained by the strong magnetosphere compression and constriction of the magnetosheath region where the soft X-ray emissivity maximum is located. The simulations also indicate large excursions of the soft X-ray emissivity maximum inside the magnetosheath as the IMF magnitude and SW dynamic pressure fluctuate particularly for radial and ecliptic IMF orientations.
The real-space Green's function code FEFF has been extensively developed and used for calculations of X-ray and related spectra, including X-ray absorption (XAS), X-ray emission (XES), inelastic X-ray scattering, and electron energy-loss spectra. The code is particularly useful for the analysis and interpretation of the XAS fine-structure (EXAFS) and the near-edge structure (XANES) in materials throughout the periodic table. Nevertheless, many applications, such as non-equilibrium systems, and the analysis of ultra-fast pump–probe experiments, require extensions of the code including finite-temperature and auxiliary calculations of structure and vibrational properties. To enable these extensions, we have developed in tandem a new version FEFF10 and new FEFF -based workflows for the Corvus workflow manager, which allow users to easily augment the capabilities of FEFF10 via auxiliary codes. This coupling facilitates simplified input and automated calculations of spectra based on advanced theoretical techniques. The approach is illustrated with examples of high-temperature behavior, vibrational properties, many-body excitations in XAS, super-heavy materials, and fits of calculated spectra to experiment.
Remote operation of small modular reactor (SMR) facilities is an appealing prospect for streamlined operation costs, safety concerns, and user convenience. The ability to guarantee security of communication channels between offsite users and nuclear reactor facilities is critical to enabling remote operations, given the sensitive nature of reactor state data. This report proposes a configuration implementing quantum key distribution (QKD) and advanced encryption standard (AES) protocols to establish secure reactor data feeds. The model outlined, which utilizes the Argonne-developed SeQUeNCe software package and custom MATLAB code, determines free-space channel losses for microwave-band communications in a variety of weather conditions and simulates QKD via ground-satellite links. This report demonstrates the feasibility of repeatable, real-time key generation and distribution for this configuration towards enabling secure wireless reactor communications.
This work details the synthesis and the crystal structures of the quaternary Zintl phases Na2CaCdSb2, Na2SrCdSb2 and Na2EuCdSb2. They are isostructural and their noncentrosymmetric structure is with the space group Pmc21 (Pearson code oP12). All structural work is carried out via single-crystal X-ray diffraction methods. The structure features [CdSb2]4– layers of corner-shared CdSb4 tetrahedra, which are stacked along the b-crystallographic axis and are separated by cations. The results from the structure refinements suggest that in addition to full cation ordering, which is typical for this structure, there also exists a possibility for an accommodation of a small degree of cation disorder.
Accurate optical surfaces are a primary driver of concentrated solar power plant performance. Errors in pointing and tracking mirrors, the canting of individual mirror facets, and the surface slope of the mirror itself can be caused by errors during assembly, transportation, wind loading, gravity, and many other sources. The tools that exist to measure these error sources today largely rely on fringe deflectometry (SOFAST, QDec, others), or photogrammetry with targets attached to the mirror surface. Since 2022, NREL has been developing a measurement method called the Reflected Target Non-intrusive Assessment (ReTNA) system. This system differs from most established methods in that we perform deflectometry with a pattern of coded targets, identified in space with photogrammetry. Reflected target systems have several advantages over traditional fringe deflectometry systems. Firstly, they can be operated in bright or ambient lighting, a challenge for fringe systems that use a projector and screen. Reflected target systems also can use a much lighter and less expensive target than projector-based systems. Lastly, 2D slope measurement can be solved from a single image, which leads to several advantages for accommodating faster measurements and smaller sized targets. These advantages make ReTNA particularly well-suited for applications where there are space or lighting constraints, like performing heliostat quality assurance on an assembly line. It's also useful when a lightweight, flexible system is needed, like for heliostat developers to quickly measure a new heliostat design at different orientations, to observe gravitational effects on the mirror surface shape. In the last year, significant improvements were made to this tool to make it more useful for these applications. These improvements were focused around validation of the ReTNA measurement system, and automation of the setup and measurement process. First, we present an improved ReTNA layout, for use on the heliostat assembly line. Next, we detail the various changes to the ReTNA software and computer vision methods to automate data collection in this new setup, and lessons learned from this process. The goal with this new setup is to perform a full heliostat surface characterization without removing the mirror from the assembly line. Lastly, we share results from several ReTNA validation studies undertaken over the last year. These include repeated ReTNA measurement on demonstration mirror facets, comparisons with other optical measurement tools, and some studies aimed at quantifying the uncertainty of ReTNA measurement under various constraints (mirror-target spacing, camera resolution, etc.). These results are compared with 2024 HelioCon performance targets, and our planned next steps for the ReTNA measurement system are presented.
A classical reduced order model (ROM) for dynamical problems typically involves only the spatial reduction of a given problem. Recently, a novel space–time ROM for linear dynamical problems has been developed [Choi et al., Space–tume reduced order model for large-scale linear dynamical systems with application to Boltzmann transport problems, Journal of Computational Physics, 2020], which further reduces the problem size by introducing a temporal reduction in addition to a spatial reduction without much loss in accuracy. The authors show an order of a thousand speed-up with a relative error of less than 10−5 for a large-scale Boltzmann transport problem. In this work, we present for the first time the derivation of the space–time least-squares Petrov–Galerkin (LSPG) projection for linear dynamical systems and its corresponding block structures. Utilizing these block structures, we demonstrate the ease of construction of the space–time ROM method with two model problems: 2D diffusion and 2D convection diffusion, with and without a linear source term. For each problem, we demonstrate the entire process of generating the full order model (FOM) data, constructing the space–time ROM, and predicting the reduced-order solutions, all in less than 120 lines of Python code. We compare our LSPG method with the traditional Galerkin method and show that the space–time ROMs can achieve O(10−3) to O(10−4) relative errors for these problems. Depending on parameter–separability, online speed-ups may or may not be achieved. For the FOMs with parameter–separability, the space–time ROMs can achieve O(10) online speed-ups. Finally, we present an error analysis for the space–time LSPG projection and derive an error bound, which shows an improvement compared to traditional spatial Galerkin ROM methods.
Generalized parton distributions (GPDs) characterize the 3-dimensional structure of hadrons, combining information about their internal quark and gluon longitudinal momentum distributions and transverse position within the hadron. The dependence of GPDs on the factorization scale Q 2 allows one to connect hard exclusive processes involving GPDs at disparate energy and momentum scales, which is needed in global analyses of experimental data. Here, in this work, we explore how finite element methods can be used to construct fast and differentiable Q 2 evolution codes for GPDs in momentum space, which can be used in a machine learning framework. We show numerical benchmarks of the methods' accuracy, including a comparison to an existing evolution code from PARTONS/APFEL++, and provide a repository where the code can be accessed.
In vitro directed evolution techniques, such as mRNA display, enable peptide ligand discovery and optimization. However, physical libraries that rely on a genetic code can only search a small fraction of sequence space due to inherent biases in the genetic code and experimental limitations. To address this challenge, denoising diffusion implicit models (DDIMs) are applied to generate novel peptide ligands against B‐cell lymphoma extra‐large (Bcl‐x L ), a key cancer target. Starting with high‐throughput sequencing data from previous selections, a DDIM is trained to produce novel sequences with high affinity binding. Experimental validation confirms that most generated sequences are functionally equivalent to the original library members for Bcl‐x L binding and demonstrated comparable binding kinetics and affinity relative to the wildtype and nearest original neighbors. Importantly, this approach generated rare sequences not easily accessible via mutation and directed evolution. These results indicate that DDIMs can complement and expand directed evolution data, efficiently exploring underrepresented regions of sequence space. This approach provides a broadly applicable framework for accelerating ligand discovery and optimizing molecular properties across diverse targets.
The proposed code, named FROG (Fusion neutron Generator) is built upon the open-source particle transport Monte Carlo toolkit Geant4. Geant4 provides C++ classes that can be leveraged to build application-specific codes dealing with the transport of particles through matter. Geant4-based codes are applied in high-energy particle physics experiments, medical applications, shielding, and space applications for example. The FROG code allows the user to define the geometry of a neutron converter device shaped as a hollow cylinder, where a neutron breeding material such as lithium deuteride (LiD) is cladded by two concentric cylinders. Such neutron converter is then placed inside a regular nuclear fission reactor, where thermal neutrons will react with the neutron breeder material (typically, Lithium 6), and through a series of reactions, will generate high-energy neutrons – neutrons whose kinetic energy are around 14 MeV. The hollowed central portion can hold a specimen that will be bombarded by high-energy neutrons created inside the neutron breeding material. Figuratively speaking, this type of device transforms neutrons from thermal (~0.625 eV) to fusion (~14 MeV) energies and is sometimes termed “fusion-to-thermal neutron converters” in the literature. The code consists of C++ source file compiled and linked to generate an executable. The user can select the dimensions of the converter (radius, length, and thickness of the breeder material), the breeder material type, the cladding material, and the specimen material that will be activated or irradiated. As input, the neutron flux for a specific location inside a reactor, for instance, positions in ATR, is required. As output, the code predicts the number of high-energy neutrons produced, the total neutron flux and fluence as well as its detailed spectrum. The physics involved in such device is very complex, as it requires modeling neutron transport, light-ion (tritons) transport, as well as fusion reactions. The Geant4 toolkit provides the required physical models.