SEARCH · Engineering Papers
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.
We report the results of simulations and experiments on the formation of light distributions by holograms and diffraction elements displayed by digital micromirror devices. Examples of successful experiments on the implementation of optical-digital systems based on them are demonstrated. (paper)
In May 2016, the U.S. Nuclear Regulatory Commission (NRC) staff provided a digital instrumentation and control (I&C) regulatory infrastructure integrated action plan to the NRC for approval. One of the objectives of that plan was to establish a clear regulatory structure with reduced regulatory uncertainty to enable the expanded safe use of digital I&C in commercial nuclear reactors while continuing to ensure safety and security. To achieve this end, the NRC, with collaboration from industry, developed a streamlined License Amendment Request Alternate Review (AR) process for safety-related (SR) digital I&C upgrades. In spite of this effort, the industry has remained reluctant to perform such I&C upgrades because of perceived regulatory and financial risks associated with being the first or an early adopter of the AR process for SR I&C upgrades. The U.S. Department of Energy Light Water Reactor Sustainability Program at the Idaho National Laboratory performed Initial Scoping Phase research to help break this impasse by supporting a SR I&C Pilot Upgrade, working with MPR Associates, ScottMadden Inc., and Exelon Generation. Exelon’s Limerick Generating Station (LGS) was selected as the target for this research. This paper summarizes the Initial Scoping Phase engineering and operations, licensing, and project management activities necessary to bound the scope, schedule, and estimated cost of the project sufficiently to enable utility management authorization of Conceptual Design Phase activities. These efforts and associated products are intended to provide a template to support larger industry efforts to perform similar upgrades as a foundation stone for a digital transformation that will improve plant safety, reliability, and operational performance while lowering plant total cost of ownership. As a result of the combined effort of Exelon Generation and research participants, Conceptual Design Phase activities for the subject upgrade at LGS were approved by Exelon. Further, the U.S. Department of Energy also awarded a $50 million cost share award to Exelon in order to pave the way for SR I&C modernization and associated control room upgrades across the U.S. nuclear fleet. Additional research reports are planned for the Conceptual Design Phase, Detailed Design Phase, and the Implementation Phase of the LGS project to document the process followed and promulgate lessons learned to industry.
The DIII-D has upgraded its upper divertor to a modular system using copper alloy pedestals to alter the divertor geometry without changing the vessel structure. Six new graphite tiles were designed for the shape and volume rise (SVR) divertor. Here, these new SVR tiles facilitated the formation of a poloidal array of 27 surface eroding thermocouples (SETCs) in the upper divertor region. At one location, a specialized recessed SETC, paired with a standard flush SETC, was installed in one of the ceiling tiles to provide comprehensive heat flux measurements, distinguishing between charged and noncharged particle contributions. Upgrades were made to the SETC system in the small-angle slot (SAS) divertor to improve overall performance. These upgrades included optimizing the feedthrough system to double the thermocouple cable capacity and reallocating cables from the SAS area to the SVR divertor. A compact isolation amplifier system with a fixed gain of 41 was employed to improve the signal level and minimize interference. Additionally, two analog-to-analog fiber systems were implemented for transmitting thermocouple signals over a single fiber, significantly reducing both noise levels and costs. The newly installed SETCs in the SVR divertor successfully completed initial commissioning testing. The SETCs captured the in-out asymmetry in the power distribution between the inner and outer strike points and demonstrated the dependency of the heat flux profile on the outer strike point location. During divertor detachment, heat flux mitigation was noted at the outer strike point location, while significant heat flux contributions from neutral particles were measured in the SVR divertor.
Abstract A computational approach by an implementation of the principle component analysis (PCA) with K -means and Gaussian mixture (GM) clustering methods from machine learning algorithms to identify structural and dynamical heterogeneities of supercooled liquids is developed. In this method, a collection of the average weighted coordination numbers ( W C N s ‾ ) of particles calculated from particles’ positions are used as an order parameter to build a low-dimensional representation of feature (structural) space for K -means clustering to sort the particles in the system into few meso-states using PCA. Nano-domains or aggregated clusters are also formed in configurational (real) space from a direct mapping using associated meso-states’ particle identities with some misclassified interfacial particles. These classification uncertainties can be improved by a co-learning strategy which utilizes the probabilistic GM clustering and the information transfer between the structural space and configurational space iteratively until convergence. A final classification of meso-states in structural space and domains in configurational space are stable over long times and measured to have dynamical heterogeneities. Armed with such a classification protocol, various studies over the thermodynamic and dynamical properties of these domains indicate that the observed heterogeneity is the result of liquid–liquid phase separation after quenching to a supercooled state.
Plasma accelerators can generate high-energy, high-brightness electron beams over centimeter-scale distances, offering novel pathways to compact x-ray free-electron lasers. Generating beams pre-bunched at the desired radiation wavelength would significantly enhance longitudinal coherence and reduce saturation length. Plasma density-modulated downramp injection offers an in-situ way to generate such beams with nanometer-scale bunching. Here we report the first experimental evidence of this mechanism in a laser-driven wakefield accelerator, showing that modulated density downramps generate modulated electron energy spectra absent in unmodulated cases. Particle-in-cell simulations reproduce these observations and reveal bunching factors of 0.05 at 0.4 μm, with a compression factor of approximately 7. Building on this demonstration, we propose a beam-driven implementation for FACET-II to generate multi-GeV beams pre-bunched at hundreds of nanometers wavelength, with sub-micrometer emittance, kiloampere peak current, and sub-percent slice energy spread. Two-stage magnetic compression enables tunable bunching from optical to extreme ultraviolet wavelengths while achieving peak currents exceeding 100 kA. Coherent transition radiation calculations confirm diagnostic feasibility. This approach offers a promising path towards compact, high-energy pre-bunched electron sources for advanced photon science applications.
Abstract Using the quantum information model of dynamical diffraction we consider a neutron cavity composed of two perfect crystal silicon blades capable of containing the neutron wavefunction. We show that the internal confinement of the neutrons through Bragg diffraction can be modelled by a quantum random walk. Furthermore, we introduce a toolbox for modelling crystal imperfections such as surface roughness and defects. Good agreement is found between the simulation and the experimental implementation, where leakage beams are present, modelling of which is impractical with the conventional theory of dynamical diffraction. Analysis of the standing neutron waves is presented in regards to the crystal geometry and parameters; and the conditions required for well-defined bounces are derived. The presented results enable new approaches to studying the setups utilizing neutron confinement, such as the experiments to measure neutron magnetic and electric dipole moments.
We develop an analytic method for implementing the IR-resummation of [1], which allows one to correctly and consistently describe the imprint of baryon acoustic oscillations (BAO) on statistical observables in large-scale structure. We show that the final IR-resummed correlation function can be computed analytically without relying on numerical integration, thus allowing for an efficient and accurate use of these predictions on real data in cosmological parameter fitting. In this work we focus on the one-loop correlation function and the BAO peak. We show that, compared with the standard numerical integration method of IR-resummation, the new method is accurate to better than 0.2%, and is quite easily improvable. We also give an approximate resummation scheme which is based on using the linear displacements of a fixed fiducial cosmology, which when combined with the method described above, is about six times faster than the standard numerical integration. Finally, we show that this analytic method is generalizable to higher loop computations.
ICARUS T600 is the far detector of the Short Baseline Neutrino program at Fermilab (U.S.A.), which foresees three Liquid Argon Time Projection Chambers along the Booster Neutrino Beam line to search for LSND-like sterile neutrino signal. The T600 detector underwent a significant overhauling process at CERN, introducing new technological developments while maintaining the already achieved performances. The realization of a new liquid argon scintillation light detection system is a primary task of the detector overhaul. As the detector will be subject to a huge flux of cosmic rays, the light detection system should allow the 3D reconstruction of events contributing to the identification of neutrino interactions in the beam spill gate. The design and implementation of the new scintillation light detection system of ICARUS T600 is described.
The Auger Muon Infill Ground Array (AMIGA) is part of the AugerPrime upgrade of the Pierre Auger Observatory. It consists of particle counters buried 2.3 m underground next to the water-Cherenkov stations that form the 23.5 km 2 large infilled array. The reduced distance between detectors in this denser area allows the lowering of the energy threshold for primary cosmic ray reconstruction down to about 1017 eV. At the depth of 2.3 m the electromagnetic component of cosmic ray showers is almost entirely absorbed so that the buried scintillators provide an independent and direct measurement of the air showers muon content. This work describes the design and implementation of the AMIGA embedded system, which provides centralized control, data acquisition and environment monitoring to its detectors. The presented system was firstly tested in the engineering array phase ended in 2017, and lately selected as the final design to be installed in all new detectors of the production phase. The system was proven to be robust and reliable and has worked in a stable manner since its first deployment.
The LINAC at the Los Alamos Neutron Science Center (LANSCE) has been utilizing the Delta-t method to match the RF cavities to the design acceleration parameters since its commissioning in 1972. The differences in time-of-flight between two subsequent Beam Position and Phase Monitors (BPPMs) are measured with both accelerated and drifting beams, depending on the whether the module is set to on or off. The algorithm optimizes the module amplitude and phase via iterative measurements if the initial phase is in the vicinity of the design value. With an upgrade to a faster readout system, a scan over the whole RF cavity phase range requires relatively less time than the classical optimization procedure. The Phase Scan Signature Matching (PSSM) method provides a time-efficient method that ensures the phase selection lands on the bunching side and empowers future analyses to build module-specific models. The PSSM also utilizes a direct model to determine the correct amplitude to sub-percent level instead of using linearized matrices. Furthermore, lacking a reliable energy measurement method in the LINAC, we measure the beam phases at two downstream locations to increase the precision of energy measurements. In this letter, we also discuss the sensitivities of PSSM, error propagation, and the implementation results for the 2019 and 2020 beam cycles.
The PTOLEMY transverse drift filter is a new concept to enable precision analysis of the energy spectrum of electrons near the tritium β-decay endpoint. Here, we detail the implementation and optimization methods for successful operation of the filter for electrons with a known pitch angle. We present the first demonstrator that produces the required magnetic field properties with an iron return-flux magnet. Two methods for the setting of filter electrode voltages are detailed. The challenges of low-energy electron transport in cases of low field are discussed, such as the growth of the cyclotron radius with decreasing magnetic field, which puts a ceiling on filter performance relative to fixed filter dimensions. Additionally, low pitch angle trajectories are dominated by motion parallel to the magnetic field lines and introduce non-adiabatic conditions and curvature drift. To minimize these effects and maximize electron acceptance into the filter, we present a three-potential-well design to simultaneously drain the parallel and transverse kinetic energies throughout the length of the filter. These optimizations are shown, in simulation, to achieve low-energy electron transport from a 1 T iron core (or 3 T superconducting) starting field with initial kinetic energy of 18.6 keV drained to < 10 eV (< 1 eV) in about 80 cm. This result for low field operation paves the way for the first demonstrator of the PTOLEMY spectrometer for measurement of electrons near the tritium endpoint to be constructed at the Gran Sasso National Laboratory (LNGS) in Italy.
SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.
A glow-curve analysis code was previously developed in C++ to analyze thermoluminescent dosimeter glow curves using automated peak detection while applying a first-order kinetics model. A newer version of this code was implemented to improve the automated peak detection and curve fitting models. The Stochastic Gradient Descent Algorithm was introduced to replace the prior approach of taking first and second-order derivatives for peak detection. Additionally, early stopping mechanisms were invoked to improve the previously used Levenberg-Marquardt Algorithm employed for curve fitting. The two software versions were compared through glow curve analysis of different thermoluminescent dosimeter materials and calculation of the corresponding figures of merit. Altogether improvements were shown, namely an increase in the number of peaks detected and a reduction of the mean figure of merit by approximately 46%.
Frenkel excitons are of interest for ultrafast switching applications due to the femtosecond timescale at which coherent exciton transfer between chromophores takes place. The Frenkel Hamiltonian, which governs the dynamics of Frenkel excitons, belongs to a class of Hamiltonians that enable universal quantum computation. It is thus of interest to determine how a complete set of exciton-based gates for quantum computing could be constructed as aggregates of chromophores. Here, we demonstrate that a controlled phase shifter can be constructed as a pair of exciton transmission lines that employ two types of chromophores and that are cross coupled by a two-exciton interaction. This element facilitates the construction of controlled basis-change gates, thereby enabling the implementation of a complete set of exciton-based gates for universal quantum computation.
Here, we propose an implementation of bivariate bicycle codes [S. Bravyi et al., Nature (London) 627, 778 (2024)] based on long-range Rydberg gates between stationary neutral atom qubits. An optimized layout of data and ancilla qubits reduces the maximum Euclidean communication distance needed for nonlocal parity-check operators. An optimized Rydberg gate pulse design enables 𝖢𝖹 entangling operations with fidelity $\mathscr{F}$ >0.999 at a distance greater than 12 µm. The combination of optimized layout and gate design leads to a quantum error correction cycle time of ∼1.28 ms for a [[144,12,12]] code, which is nearly a factor-of-two improvement over previous designs.
Computing the GW quasiparticle band structure and Bethe-Salpeter equation (BSE) absorption spectra for materials with spin-orbit coupling have commonly been done by treating GW corrections and spin-orbit coupling (SOC) as separate perturbations to density-functional theory. However, accurate treatment of materials with strong spin-orbit coupling (such as many topological materials of recent interest, and thermoelectrics) often requires a nonperturbative approach using spinor wave functions in the Kohn-Sham equation and GW/BSE. Such calculations have only recently become available, in particular for the BSE. Here, we have implemented this approach in the plane-wave pseudopotential GW/BSE code BerkeleyGW, which is highly parallelized and widely used in the electronic-structure community. We present reference results for quasiparticle band structures and optical absorption spectra of solids with different strengths of spin-orbit coupling, including Si, Ge, GaAs, GaSb, CdSe, Au, and Bi 2 Se 3 . The calculated quasiparticle band gaps of these systems are found to agree with experiment to within a few tens of meV. SOC splittings are found to be generally in better agreement with experiment, including quasiparticle corrections to band energies. The absorption spectrum of GaAs is not significantly impacted by the inclusion of spin-orbit coupling due to its relatively small value (0.2 eV) in the Λ direction, while the absorption spectrum of GaSb calculated with the spinor GW/BSE captures the large spin-orbit splitting of peaks in the spectrum. For the prototypical topological insulator Bi 2 Se 3 , we find a drastic change in the low-energy band structure compared to that of DFT, with the spinorial treatment of the GW approximation correctly capturing the parabolic nature of the valence and conduction bands after including off-diagonal self-energy matrix elements. We present the detailed methodology, approach to spatial symmetries for spinors, comparison against other codes, and performance compared to spinless GW/BSE calculations and perturbative approaches to SOC. This work aims to spur further development of spinor GW/BSE methodology in excited-state research software and enables a more accurate and detailed exploration of electronic and optical properties of materials containing elements with large atomic numbers.
Predicting the properties of the matter ejected during and after a neutron star merger is crucial to our ability to use electromagnetic observations of these mergers to constrain the masses of the neutron stars, the equation of state of dense matter, and the role of neutron star mergers in the enrichment of the Universe in heavy elements. Furthermore, our ability to reliably provide such predictions is however limited by a broad range of factors, including the finite resolution of numerical simulations, their treatment of magnetic fields, neutrinos, and neutrino-matter interactions, and the approximate modeling of the equation of state of dense matter. In this manuscript, we study specifically the role that a small residual eccentricity and different implementations of the same equation of state have on the matter ejected during the merger of a 1.3M ⊙ –1.4M ⊙ binary neutron star system. We find that a residual eccentricity e ~ 0.01, as measured ~ 4–6 orbits before merger, causes O(25%–30%) changes in the amount of ejected mass, mainly due to changes in the amount of matter ejected as a result of core bounces during merger. We note that O(1%) residual eccentricities have regularly been used in binary neutron star merger simulations as proxy for circular binaries, potentially creating an additional source of error in predictions for the mass of the dynamical ejecta.
The speed of elementary quantum gates, particularly two-qubit gates, ultimately sets the limit on the speed at which quantum circuits can operate. In this work, we experimentally demonstrate commonly used two-qubit gates at nearly the fastest possible speed allowed by the physical interaction strength between two superconducting transmon qubits. We achieve this quantum speed limit by implementing experimental gates designed using a machine-learning-inspired optimal control method. Importantly, our method only requires the single-qubit drive strength to be moderately larger than the interaction strength to achieve an arbitrary two-qubit gate close to its analytical speed limit with high fidelity. Thus the method is applicable to a variety of platforms, including those with comparable single-qubit and two-qubit gate speeds, or those with always-on interactions. We expect our method to offer significant speedups for non-native two-qubit gates that are typically achieved with a long sequence of single-qubit and native two-qubit gates.