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.
Final Technical Report for DOE SC0019287 award titled Radio Frequency Algorithms. This project supported the design, simulation, analysis, and optimization of robust algorithms to configure, optimize, and control Radio Frequency and Low-Level RF systems for several modern accelerator systems.
Block Encoding (BE) is a crucial subroutine in many modern quantum algorithms, including those with near-optimal scaling for simulating quantum many-body systems, which often rely on Quantum Signal Processing (QSP). Currently, the primary methods for constructing BEs are the Linear Combination of Unitaries (LCU) and the sparse oracle approach. In this work, we demonstrate that QSP-based techniques, such as Quantum Singular Value Transformation (QSVT) and Quantum Eigenvalue Transformation for Unitary Matrices (QETU), can themselves be efficiently utilized for BE implementation. Specifically, we present several examples of using QSVT and QETU algorithms, along with their combinations, to block encode Hamiltonians for lattice bosons, an essential ingredient in simulations of high-energy physics. We also introduce a straightforward approach to BE based on the exact implementation of Linear Operators Via Exponentiation and LCU (LOVE-LCU). We find that, while using QSVT for BE results in the best asymptotic gate count scaling with the number of qubits per site, LOVE-LCU outperforms all other methods for operators acting on up to qubits, highlighting the importance of concrete circuit constructions over mere comparisons of asymptotic scalings. Using LOVE-LCU to implement the BE, we simulate the time evolution of single-site and two-site systems in the lattice theory using the Generalized QSP algorithm and compare the gate counts to those required for Trotter simulation.
The growth of the aviation industry coupled with its dependence on energy dense, liquid fuels has brought sustainable aviation fuel (SAF) research to the forefront of the biofuels community. Petroleum refineries will need to decide how to satisfy the projected increase in jet fuel demand with either capital investments to debottleneck current operations or by integrating bio-blendstocks. This work seeks to compare jet production strategies on a risk-adjusted, economic performance basis using Monte-Carlo simulation and refinery optimization models. Additionally, incentive structures aiming to de-risk initial SAF production from the refiner’s perspective are explored. Results show that market sensitive incentives can reduce the financial risks associated with producing SAFs and deliver marginal abatement costs ranging between 136-182 $/Ton-CO2e.
Modern utility-scale wind turbines are evolving toward larger, lighter, and more flexible designs to meet the growing demand for renewable energy while minimizing logistical costs. However, these advancements in lightweight design result in heightened aeroelastic sensitivity, leading to complex interactions which affect the rotor’s capacity to withstand aerodynamic loading and the cascading effects that manifest in the wake’s vortex-structure evolution under variable atmospheric conditions. In this paper, we analyze the influence of stream-wise fluctuating atmospheric flow conditions on wind turbines with large, flexible rotors through simulations of the National Rotor Testbed (NRT) turbine, located at Sandia National Labs’ Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas. The Common Ordinary Differential Equation Framework (CODEF) modeling suite is used to simulate wind turbine aeroelastic oscillatory behavior and wind farm vortex–wake interactions for a range of conditions with spatially variant atmospheric flow. CODEF solutions for turbine operation in wind conditions featuring only one parameter fluctuation are compared to wind conditions with several wind parameter variations in combination. By isolating individual inflow variations and comparing them to multi-parameter scenarios, we determine the contributions of each atmospheric factor to rotor dynamics, wake evolution, and downstream wind farm interactions. The purpose of this paper is to analyze the effects of spatial variations in atmospheric flow on the topological evolution of wind turbine vortex wakes, which constitutes a gap in the current understanding of wind turbine wake dynamics. The insights gained from this study are particularly valuable for the development of wind farm control strategies aimed at mitigating the adverse effects of wake interactions, enhancing energy capture, and improving the overall stability of wind farm operations. With these insights, we aim to contribute to the development of modeling and simulation tools to optimize utility-scale wind power plants operating in diverse atmospheric environments.
Passive surface acoustic wave (SAW) devices are attractive candidates for continuous wireless monitoring of corrosion in large infrastructures. However, acoustic loss in the aqueous medium and limited read range usually create challenges in their widespread use for monitoring large systems such as oil and gas (O&G) pipelines, aircraft, and processing plants. This paper presents the investigation of impedance-loaded reflective delay line (IL-RDL) SAW devices for monitoring metal corrosion under O&G pipeline-relevant conditions. Specifically, we studied the effect of change in resistivity of a reflector on the backscattered signal of an RDL and investigated an optimal range through simulation. This was followed by the experimental demonstrations of real-time monitoring of Fe film corrosion in pressurized (550 psi) humid CO 2 conditions. Additionally, remote monitoring of Fe film corrosion in an acidic solution inside a 70 m carbon steel pipe was demonstrated using guided waves. This paper also suggests potential ways to improve the sensing response of IL-RDLs.
In this paper we apply ideas from information theory to create a method for the design of optimal filters for photometric redshift estimation. We show the method applied to a series of simple example filters in order to motivate an intuition for how photometric redshift estimators respond to the properties of photometric passbands. We then design a realistic set of six filters covering optical wavelengths that optimize photometric redshifts for z <= 2.3 and i < 25.3. We create a simulated catalog for these optimal filters and use our filters with a photometric redshift estimation code to show that we can improve the standard deviation of the photometric redshift error by 7.1% overall and improve outliers 9.9% over the standard filters proposed for the Large Synoptic Survey Telescope (LSST). We compare features of our optimal filters to the LSST and find that the LSST filters incorporate key features for optimal photometric redshift estimation. Finally, we describe how information theory can be applied to a range of optimization problems in astronomy.
The emergence of connected and automated vehicles and smart cities technologies create the opportunity for new mobility modes and routing decision tools, among many others. To achieve maximum mobility and minimum energy consumption, it is critical to understand the energy cost of decisions and optimize accordingly. The Route Energy prediction model (RouteE) enables accurate estimation of energy consumption for a variety of vehicle types over trips or sub-trips where detailed drive cycle data are unavailable. Applications include vehicle route selection, energy accounting and optimization in transportation simulation, and corridor energy analyses, among others. The software is a Python package that includes a variety of pre-trained models from the National Renewable Energy Laboratory (NREL). However, RouteE also enables users to train custom models using their own data sets, making it a robust and valuable tool for both fast calculations and rigorous, data-rich research efforts. The pre-trained RouteE models are established using NREL's Future Automotive Systems Technology Simulator paired with approximately 1 million miles of drive cycle data from the Transportation Secure Data Center, resulting in energy consumption behavior estimates over a representative sample of driving conditions for the United States. Validations have been performed using on-road fuel consumption data for conventional and electrified vehicle powertrains. Transferring the results of the on-road validation to a larger set of real-world origin-destination pairs, it is estimated that implementing the present methodology in a green-routing application would accurately select the route that consumes the least fuel 90% of the time. The novel machine learning techniques used in RouteE make it a flexible and robust tool for a variety of transportation applications.
Distributed scientific workflows are becoming more important with the interest in incorporating AI into their loops. A critical programming and performance question is how to compose workflow tasks when data is produced on one system but must be consumed on another. Since the dominant technique is composition with remote I/O, this paper explores its performance expectations. We describe BigFlowSim, a workflow I/O simulator that captures key implementation choices for remote I/O, including intensity, reuse, locality, access pattern, and data movement.With BigFlowSim, we generate a synthetic benchmark. We quantify the effects of each parameter with a performance sensitivity study. We explain trends in terms of data movement reduction and show that, under certain conditions, it is possible to establish a total order among most parameters. We apply these insights to a high energy physics workflow, Belle II Monte Carlo and simulate several I/O optimizations. Speedups range from 5% to 2×, without changing compute time.
Particle physics is an ever-changing field seeking to explain phenomena that lie at the limits of our current knowledge. Among discovered and hypothesized interactions, Coherent Elastic Neutrino-Nucleus Scattering (CE?NS) stands out as a fascinating and elusive phenomenon, offering a unique window into the nature of neutrinos and their interactions with matter. CE?NS, a process predicted by the Standard Model, represents one of the subtle ways in which neutrinos can interact with atomic nuclei. Unlike other neutrinonucleus interactions, CE?NS occurs coherently, meaning that the entire nucleus collectively responds to the neutrino-mediated weak force, resulting in small energy transfers to the nucleus as a whole. The extremely rare nature of CE?NS events presents a formidable challenge to experimenters seeking to detect and study these interactions. Recent advances in detector technologies and the implementation of large-scale experiments, such as coherent neutrino scattering experiments at research facilities around the world, have opened up new avenues for exploring this elusive phenomenon. These experiments aim not only to observe CE?NS directly but also to extract valuable information on neutrino properties, nuclear structure, and potential deviations from the Standard Model. This thesis aims to study the feasibility of an experiment to study CE?NS at The Thomas Jefferson National Accelerator Facility, Newport News, VA, USA. This would be an experiment that exploits the production of neutrinos by an electron beam dumped on a thick target (the beam dump). The large quantity of neutrinos produced and their energy profile are well suited to the purpose of studying the CE?NS. In my thesis, I will start from previous studies on CE?NS and from theoretical considerations regarding the process, then quantify the expected events for a hypothetical experiment and possible backgrounds through simulations aimed at optimizing the geometry of the future detector. I will also show some tests on selected detector components that will be part of a first prototype.
In this research three artificial intelligent (AI)-based techniques are proposed to regulate the voltage and frequency of a grid-connected inverter. The increase in the penetration of renewable energy sources (RESs) into the power grid has led to the increase in the penetration of fast-responding inertia-less power converters. The increase in the penetration of these power electronics converters changes the nature of the conventional grid, in which the existing kinetic inertia in the rotating parts of the enormous generators plays a vital role. The concept of virtual inertia control scheme is proposed to make the behavior of grid connected inverters more similar to the synchronous generators, by mimicking the mechanical behavior of a synchronous generator. Conventional control techniques lack to perform optimally in nonlinear, uncertain, inaccurate power grids. Besides, the decoupled control assumption in conventional VSGs makes them nonoptimal in resistive grids. The neural network predictive controller, the heuristic dynamic programming, and the dual heuristic dynamic programming techniques are presented in this research to overcome the draw backs of conventional VSGs. The nonlinear characteristics of neural networks, and the online training enable the proposed methods to perform as robust and optimal controllers. The simulation and the experimental laboratory prototype results are provided to demonstrate the effectiveness of the proposed techniques.
Shock wave electrodialysis (SWED) is a highly promising technique for energy-efficient ion separation in the context of a circular economy. This paper presents a approach way of modeling and improving SWED using a two-dimensional multicell model combined with the COMSOL program and response surface methodology. The model integrates the Nernst-Planck equation, Darcy's law, and first-order electroosmosis to examine the local concentration, flux of ionic species, distribution of current, and velocity of flow in SWED cells under various operating conditions. We first illustrate the clear depiction of concentration, velocity, and electric potential distribution through contours which aids in identifying optimal operating conditions and designing scalable SWED systems. The results emphasize the significance of surface charge density and voltage in influencing the features of shock waves for obtaining effective ion separation while optimizing energy consumption and improving current efficiency by controlling the retention time of feed flow. Here, this study defines two crucial characteristics of shock waves, namely the length of the flat depletion zone of a fully developed shock wave (shock wave height) and the distance of shock wave propagation (shock wave length). These properties significantly impact separation performance, as determined by the simulation results. Additionally, the response surface methodology is incorporated with the COMSOL models to develop predictive models and graph responses, enabling a more comprehensive understanding of the interactions between parameters and performance indicators, such as removal ratio, energy consumption, and water recovery. Finally, this work suggests design tactics for expanding SWED processes and outlines potential areas for further research. This research provides valuable insights into the prospective applications, design optimization, and scalability of SWED in the field of electrokinetic separation technologies for green chemistry and a circular economy.
Uranium silicide (U 3 Si 2 ) is a candidate material for the high-density nuclear fuel in commercial light water reactors [1], [2]. Its higher uranium density, 11.3 g-U/cm3, compared to that of uranium dioxide (UO 2 ), 9.7 g-U/cm3, can improve the performance of a nuclear reactor while using low enriched uranium (LEU) and diversify the choice of cladding materials [1]–[3]. It also has a higher thermal conductivity than UO 2 , which can reduce the thermal stress on the material caused by a temperature gradient across the fuel pellet and provide a larger margin for some postulated accidents [1], [2], [4]. Furthermore, compared to U3Si, another high-density fuel candidate, it has better resistance to in-pile swelling due to less irradiation-induced rapid amorphization [1], [3]. Corresponding to its importance in nuclear engineering, many previous studies have reported the properties of U3Si2. Experiments showed that U 3 Si 2 is a paramagnetic (PM) metal, where a slight linear increase in magnetic susceptibility was measured with increasing temperature [5], [6]. In addition, thermodynamic quantities such as thermal expansion coefficient, heat capacity, and thermal conductivity were experimentally determined over a wide temperature range [1], [7], [8]. In several computational studies, ab initio atomistic simulations based on density functional theory (DFT) were performed to calculate various properties including elastic constants, electronic density of states (DOS), and phonon dispersion curves [9]–[12]. Nevertheless, thermal neutron scattering cross sections, which are critical to the prediction of the parameters in reactor physics that are ultimately related to reactor criticality, have not yet been evaluated for U3Si2. The scattering cross section can be calculated from the phonon DOS, or the energy spectrum of lattice vibrations, of the crystalline system [13], [14]. However, there is also no experimental data available for the phonon DOS of U 3 Si 2 . While some computational studies reported the phonon DOS and/or dispersion curves from ab initio simulations [9]–[12], the accuracy cannot be guaranteed because it is unclear whether the spin-polarization behavior of PM U 3 Si 2 was properly described. In the present study, the thermal neutron scattering cross section for U 3 Si 2 is evaluated for the first time by calculating the phonon DOS for U3Si2 from ab initio lattice dynamics (AILD) simulations based on DFT. First, U 3 Si 2 is modeled based on the experimental structure, and AILD simulations are performed on the modeled U3Si2 to optimize the structure. Next, AILD simulations are performed for supercells with atomic displacement to calculate Hellmann-Feynman forces. Based on the calculated forces, partial phonon DOSs for U and Si are obtained, and the thermal neutron scattering law (TSL) for U 3 Si 2 is finally evaluated. To verify the accuracy of the calculations in the present study, the calculation results are compared with experimental data on the structure and heat capacity of U3Si2 [1], [7], [8], [15].
New powertrain solutions are needed to address societal challenges stemming in part from a transportation sector which relies heavily on combustion of conventional hydrocarbon fuels like gasoline. One way to advance new solutions is through predictive simulations. Predictive simulations can potentially limit the extent of experimental validations, reduce time spent in engineering design cycles, unlock new strategies for high efficiency combustion with power density, and minimize tailpipe emissions. However, the simulation tools currently available are either too computationally expensive, inadequate in their accuracy, or a combination of the two. A DOE-funded consortium of national laboratories called the Partnership to Advance Combustion Engines (PACE) seeks to address this gap by rapidly delivering new knowledge and tools which enable market-competitive powertrain solutions for light-duty vehicles. Approaches to modeling combustion in powertrain systems typically incorporate computational fluid dynamics (CFD) simulations with chemical kinetic models. These CFD simulations generally use surrogate fuels featuring a limited number of existing components in a reduced kinetic model to limit computational costs. Such surrogate fuels often sacrifice matching several or more combustion and physical properties of the target gasoline fuel. Validation of surrogate fuels is also often only pursued for a small subset of standard metrics such as research or motor octane number (RON, MON). The purpose of this project is to develop an optimal E10 gasoline surrogate fuel suitable for tasks across PACE to facilitate rapid common analysis and progression toward consortium goals. This project received contributions from tasks led by eight PACE principal investigators. The project funding here reflects all tasks under Pitz, whose tasks also contributed to additional PACE projects. Based on the optimal gasoline surrogate fuel composition, a reduced chemical kinetic model will be shared with PACE researchers, industry, and the broader combustion community.
Abstract Scenario development in tokamaks is an open area of investigation that can be approached in a variety of different ways. Experimental trial and error has been the traditional method, but this required a massive amount of experimental time and resources. As high fidelity predictive models have become available, offline development and testing of proposed scenarios has become an option to reduce the required experimental resources. The use of predictive models also offers the possibility of using a numerical optimization process to find the controllable inputs that most closely achieve the desired plasma state. However, this type of optimization can require as many as hundreds or thousands of predictive simulation cases to converge to a solution; many of the commonly used high fidelity models have high computational burdens, so it is only reasonable to run a handful of predictive simulations. In order to make use of numerical optimization approaches, a compromise needs to be found between model fidelity and computational burden. This compromise can be achieved using neural networks surrogates of high fidelity models that retain nearly the same level of accuracy as the models they are trained to replicate while reducing the computation time by orders of magnitude. In this work, a model-based numerical optimization tool for scenario development is described. The predictive model used by the optimizer includes neural network surrogate models integrated into the fast Control-Oriented Transport simulation framework. This optimization scheme is able to converge to the optimal values of the controllable inputs that produce the target plasma scenario by running thousands of predictive simulations in under an hour without sacrificing too much prediction accuracy.
Concentrated solar power (CSP) with thermal energy storage (TES) has the potential to achieve grid parity. This can be realized by operating CSP systems at temperatures above 700 °C with high-efficiency sCO 2 power cycles. However, operating CSP systems at such temperatures poses several challenges, among which the design of solar receivers to accommodate increased thermal loads is critical. To this end, this work explores and optimizes various swirl-inducing internal fin designs for solar receiver tubes. These fin designs not only improve the thermal performance of receiver tubes but also levelize temperature unevenness caused by non-uniform thermal loading. In this work, the geometric parameters of the fin designs are optimized to maximize the Nusselt number with a constraint on the friction factor. This optimization, however, is computationally intensive, requiring hundreds of simulation calls to computational fluid dynamics (CFD) models. Here, to circumvent this problem, surrogate models are used to approximate the simulation outputs needed during the optimization. In addition, this study also examines the fin designs from an entropy generation perspective. To this end, the entropy contributions from thermal and viscous effects are quantitatively compared while varying the operational Reynolds number.
The bulk metallic electromagnetic levitation technique has a promising application prospect in condensed physics and materials science. In order to resolve the problem of weak levitation force and the small size of samples, an optimized solution based on finite element analysis (FEA) is developed to seek the best electromagnetic field distribution in the levitation zone. The influence of material properties on the distribution of magnetic flux density, surface pressure, and the Lorentz force were studied by the FEA method with various sample sizes of titanium and copper balls. The levitation capability of cylindrical coil system which consists of a single-layer coil with an inner diameter of 30 mm and the levitated object with a diameter of 20 mm was analyzed. To keep a balance between the turn number and increment efficiency of force, the best number of 7 turns is confirmed. Then, two types of optimized coils were predicted by numerical simulations, and the levitation force is improved drastically. The double-layer coil displays an outstanding levitation capability owing to the strong electromagnetic field generated by multilayer turns. In addition, 5 to 50 g aluminum, titanium, copper, and nickel are levitated by the above three kinds of coils and the levitation capacities are validated from the recorded data of the minimum operating current that can levitate the samples. The experimental results show good agreements with the simulation results, and the minimum current for optimized coils is reduced efficiently.
In this paper, a three dimensional optimization strategy is presented to improve the performance of a subsonic axial turbine exposed to high Mach number unsteady inflows encountered in pressure gain combustors. The optimization aims at altering the stator hub and shroud end wells. A multi-objective differential evolution optimizer was used to enhance specific work and efficiency. The optimization routine consisted of a multi-step approach where in a first step an optimization was carried out where the turbine was subjected to a steady high Mach inflow. The turbine stage was parametrized via an in-house turbine design tool based on Bezier curves. Seven design parameters were selected to contour the stator end walls and rotor hub end wall. The parametrized turbine stage was meshed and the steady Reynolds-Averaged Navier-Stokes equations were solved via the commercial solver CFD++ and calculation time of each stage was around 10 hours. Results of the steady optimization identified three major design parameters of stator end walls, which were taken as inputs for the next phase unsteady optimization. For the unsteady simulations, a large inlet fluctuation was imposed and the turbine was modeled with the mixing plane approach. Total calculation time of each individual was 48 hours. An efficiency improvement of 10% from the steady optimization and 13% from the unsteady optimization was found compared to the baseline case. Finally, full unsteady simulations with two stator and three rotor passages were performed on selected steady and unsteady optimized configurations.
Turbulent transport is regarded as one of the key issues in magnetic confinement nuclear fusion, both for tokamaks and stellarators. Here, in this work, we show that a significant decrease in a microstability-based proxy, as opposed to a geometric one, for the turbulent heat flux, namely the quasilinear heat flux, can be obtained in an efficient manner by coupling stellarator optimization with linear gyrokinetic simulations. This is accomplished by computing the quasilinear heat flux at each step of the optimization process, as well as the deviation from quasisymmetry, and minimizing their sum, leading to a balance between neoclassical and the turbulent transport proxy.