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.
This code provides fast Gaussian process prediction algorithms based on the MuyGPs scalable hyperparameter optimization algorithm. This code is the companion to a research paper preprint soon to be made publicly available.
An optimized plasma current ramp-down strategy is critical for safe and fast termination of plasma discharges in a tokamak demonstration fusion reactor (DEMO), both in planned and emergency scenarios, avoiding plasma disruptions and excessive heat loads to the first wall. Plasma stability limits and machine-specific technical requirements constrain the stable envelope through which the plasma must be navigated. Large amounts of auxiliary heating are required throughout the ramp-down phase, to avoid a radiative collapse in the presence of intrinsic tungsten and seeded xenon impurities, as quantitatively estimated in this work. As the plasma current is reduced, the current density becomes increasingly peaked, reflected by a growing value of the internal inductance $\ell$ $i3$ , resulting in reduced controllability of the vertical position of the plasma. The feasibility of different plasma current ramp-down rates is tested by applying an automated optimization framework embedding the RAPTOR core transport solver. Optimal time traces for plasma current $I$ $p$ $(t)$ and plasma elongation $κ(t)$ are proposed, to satisfy an $I$ $p$ -dependent upper limit on the plasma internal inductance, as obtained from vertical stability studies using the CREATE-NL code, as well as a constraint on the time evolution of $q$ 95 , to avoid an ideal MHD mode. A negative current density near the plasma edge is observed in our simulations, even for the most conservative $I$ $p$ ramp-down rate, indicating significant transient dynamics due to a large resistive time.
Extreme DC fast charging for electric vehicles (EVs) could be competitive with the internal combustion engine refueling experience and enable longer-distance travel, which could help with EV adoption and decarbonization, but these systems have high capital costs and extremely variable high-power demands. Behind-the-meter systems (BTMS) could support extreme-fast-charging (XFC) stations to increase nationwide adoption of EVs. Here, this study examines the optimal break-even levelized cost of charging (LCOC) across 96 BTMS scenarios to enable low-wait XFC stations providing 200 miles of charge in 10 min. This research simulates LCOC via synthetic XFC-capable EV loads, machine-learned battery life models from testing data, and nonlinear optimal controls, co-minimizing complex utility costs and battery replacements. An aggregate optimal BTMS design treating each EV load as equal likely gives an optimal LCOC per utility rate, the average of which is $\$$0.59/kWh. In addition, the sensitivity of optimal and off-optimal design factors, the long-life LMO/LTO chemistry, and optimized controls are analyzed. The battery control model, based on battery stressors to compare chemistries, optimizes LMO/LTO resting state of charge and cycle depth without compromising cost reduction, which enables greater flexibility in operation. The LCOC savings due to replacement reduction are small, up to $\$$0.035/kWh (6%), with an average of $\$$0.02/kWh (3.5%). Compared with gasoline stations, the aggregate XFC station design achieves comparable speed, experience of service, and cost at $\$$3.81/gal gasoline, showing that EVs can replace gasoline vehicles even for longer-distance travel.
Initial computational plasma physics scoping and a finite element method antenna modeling design workflow for a symmetric center-fed high-field side high harmonic fast wave traveling wave array (TWA) antenna are reported here. The TWA is designed to generate a test population of fast deuterium ions in an existing D–D tokamak by heating neutral beam deuterium ions, accelerating them from 80 keV to several hundred keV. The resulting fast particles are tailored to mimic key reactor energetic particle parameters with regards to exciting Alfven eigenmode instabilities, allowing for a D–D tokamak like DIII-D or ASDEX-U to replicate reactor-relevant conditions experimentally. Initial scenario scoping for high single-pass absorption as well as good preferential fast ion damping relative to electron damping was completed using the ray-tracing/Fokker–Planck codes GENRAY and CQL3D. Python RF network analysis packages were used to create a custom TWA optimization tool to inform a COMSOL flat antenna design, and Petra-M was used to study cold plasma effects. The TWA produced by this workflow has several novel features when compared to previous TWA studies, including symmetric center feeding, and passive end straps for image current cancellation for reduced impurity production. We show here that the antenna design workflow can readily produce TWA antennas optimized for reflection coefficient, image current cancellation, and launched power spectrum shape; and that a population of fast ions can be generated in the correct region of parameter space, warranting future more detailed studies.
Abstract We present the first results from the HPSoC ASIC designed for readout of Ultra-fast Silicon Detectors. The 4-channel ASIC manufactured in 65 nm CMOS by TSMC has been optimized for 50 μm thick AC-LGAD. The evaluation of the analog front end with β-particles impinging on 3 × 3 AC-LGAD arrays (500 μm pitch, 200 × 200 μm 2 metal) confirms a fast output rise time of 600 ps and good timing performance with a jitter of 45 ps. Further calibration experiments and TCT laser studies indicate some gain limitations that are being investigated and are driving the design of the second-generation pre-amplification stages to reach a jitter of 15 ps.
Electrification of heavy-duty (HD) vehicles marks an important milestone and technical challenge in the electric vehicle (EV) industry and the public grid. However, implementing EV charging at this scale will necessitate that traditional truck stops be updated with EV charging infrastructure that could represent 10's of MW in electricity consumption. Furthermore, as the transportation sector is represented as critical infrastructure, supporting resiliency considerations in EV charging infrastructure will be critical. This paper proposes an optimization-based approach for optimally designing a MW-scale microgrid charging network. This approach transforms conventional designed truck stops into a reliable HDEV charging stations capable of overnight slow charging and 30-minute to 1 hour fast charging. Using a mixed-integer linear program formulation blending capacity planning and reliability constraints, an optimal network configuration can be solved for a proposed EV charging station that includes photovoltaic and battery energy storage capabilities.
Understanding and controlling atomic-level processes at solid-liquid interfaces is key to advancing technologies in energy storage, carbon capture, critical element recovery, and materials synthesis. Many of these processes are dominated by the formation of short-lived intermediate precipitates that determine the final properties of synthesized materials. However, studying these intermediates is challenging due to their sensitivity and the reliance on trial-and-error methods. To address this, we developed an automated variable-volume mixed-flow reactor (MFR) to optimize metastable material synthesis and investigate rapid kinetic processes. This state-of-the-art MFR system, paired with an automated modeling framework, enables efficient synthesis and real-time analysis of transient phases. Benchmarking with advanced capabilities, such as wide-/small-angle X-ray scattering, allows us to resolve fast nucleation and growth dynamics that were previously inaccessible. By combining automation, ML-guided optimization, and tailored kinetic modeling, this approach provides a robust platform for improving material design and achieving precise control over solid-liquid reactions.
Here, the Comprehensive Nuclear-Test-Ban Treaty bans all nuclear tests worldwide. As part of treaty compliance, the concentration of radioactive nuclides in the atmosphere is monitored to detect nuclear weapons tests. Radioactive noble gas fission products, specifically radioxenon, can vent into the atmosphere after a nuclear weapons test, even if the test is well contained underground or underwater. Radioxenon thus serves as a signal for nuclear weapons tests. All atmospheric monitoring systems require samples of radioxenon isotopes for detector calibration, quality control, and certification. Here, we present a novel, improved method for creating samples of 133m Xe via neutron irradiation of 132 Xe in the Washington State University TRIGA reactor. 132 Xe neutron absorption results in either 133 Xe or 133m Xe—thermal neutron absorption results in 133m Xe 12% of the time, while fast neutron absorption (above ~1 MeV) results in 133m Xe ~50% of the time. To optimize the production of 133m Xe via neutron absorption in 132 Xe in the thermal TRIGA reactor, spectral tuning using an irradiation chamber is required to maximize the fraction of fast neutrons being absorbed and minimize the number of thermal neutrons interacting with the 132 Xe. We used MCNP to tally 132 Xe absorptions with the isotopic tally function, flux tallies and neutron attenuation to estimate the number of neutrons reaching the 132 Xe through the irradiation chamber, and the adjoint importance function to improve the source strength estimate. Additionally, we performed a heat transfer analysis for safety considerations. It was determined that the use of a 96% enriched 10 B boron carbide chamber, placed next to the fuel elements in reactor position D8, increases the 133m Xe/ 133 Xe activity ratio from a baseline value of 0.3 to 1.0, a 233% increase. Additionally, it was determined that the alpha heating produced in the boron does not become an unmanageable problem in the Washington State University reactor.
Renewables and especially photovoltaics (PV) have benefitted significantly from a host of incentives and policies targeted toward enhanced integration and adoption of specific energy technologies. However, with the push to move forward into a subsidy-free market framework, behind-the-meter residential PV applications have generally struggled to retain their value (unlike utility scale and commercial projects) [1]. This project focused on developing control-theoretic solutions aimed at improving the integration and interaction of behind-the-meter residential PV with other distribution system assets (controllable and non-controllable) to enhance the integrated value of residential PV. To this end, a suite of decentralized control methodologies have been developed to enable effective coordination and control of behind-the-meter residential load customers’ PV, battery storage systems (BSS), controllable loads and other similar assets within a distribution feeder. This interaction aims at procuring energy savings and, thus, energy bill savings. The main source of savings is drawn from reducing the effect of demand charge pricing and is realized at the feeder level, assuming community level interaction and management among the aforementioned assets. Optimal control of the assets is implemented with a distributed optimization methodology, leveraging consensus-based algorithms. The results gathered from the optimal control simulations demonstrates that the savings can be duly achieved and the algorithm decision times (to dynamically control asset set points, for example) are fast. As for the overall efficiency of PV+BSS systems, to procure energy savings from curtailment of the demand charge pricing effects, the optimal control is set up so as to minimize the variance of the load for all customers, throughout a feeder and throughout time in a rolling horizon scheduling with model predictive control. The control takes into account inter-temporal electrochemical storage (battery) degradation costs: specifically, we have developed a long-term lifetime model for the BSS that weighs in the effect of the degradation factor in the dispatch formulations, thus, a considerable operating cost that affects energy decision making. The levelized cost of energy (LCOE – redefined for the purpose of quantifying asset integration effectiveness through the customers’ energy cost) is shown to be below the threshold set for the combined PV+BSS topology of $ 0.14/kWh for multiple cases of PV penetration all the way up to 50%, provided that a policy of shared ownership of and savings is in place. Further, the LCOE calculated for the case before the deployment PV+BSS systems is also achievable, i.e. the deployment of PV+BSS, if planned and scheduled optimally. will have no effect on customers’ energy costs. From the control methodology viewpoint, the developed consensus-based algorithms are shown to converge for a wide range of problem cases (spanning normal operating scenarios and contingencies), guaranteeing dispatch solutions under forecasting errors, communication break-downs and cyber-security attacks. The proposed control solutions are scalable and real-time implementable, with dispatch computations and device set-point updates converging in less than 2s in most practical instances of the above events.
Recent trends in chemical engineering research point towards an increasing reliance on data-driven modeling approaches. Neural networks, for instance, have proven to be accurate when data is plentiful and high-dimensional, but in many cases, they require computationally-intensive training procedures. Here, in this work, we describe hyperplane decision trees (HT) as a highly expressive and low-compute machine learning model architecture. These models are locally linear and have linear decision boundaries, resulting in a piecewise linear model of the data. This property allows them to be converted into mixed-integer linear constraints which can be globally optimized. Our open-source PyTorch implementation of this method is a fast, flexible, and accessible way to build accurate piecewise linear models of data.
Fast secondary neutrons created at the Brookhaven Linac Isotope Producer (BLIP) facility following proton irradiation were characterized by the foil activation technique and compared with FLUKA Monte Carlo simulations. The FLUKA-simulated neutron flux was spectrally adjusted following the maximum entropy formalism using the International Reactor Dosimetry and Fusion File (IRDFF-II), with predictions agreeing with experimental measurements to within 9 % following the adjustment procedure. Here, a multitude of degrader configurations were simulated to assess the feasibility of improving the fast (E n > 20 MeV) secondary neutron yield at the proposed neutron target position (“N-slot”). A configuration where the N-slot is closest to the proton degrader produced the highest fast neutron yield, with tungsten degraders achieving the best performance. Assuming the optimized target-degrader configuration proposed in this work, we discuss potential isotope production opportunities with secondary neutrons. In most cases the yields are in the order of several mCi.
GX is a code designed to solve the nonlinear gyrokinetic system for low-frequency turbulence in magnetized plasmas, particularly tokamaks and stellarators. In GX, our primary motivation and target is a fast gyrokinetic solver that can be used for fusion reactor design and optimization along with wide-ranging physics exploration. Here, this has led to several code and algorithm design decisions, specifically chosen to prioritize time to solution. First, we have used a discretization algorithm that is pseudospectral in the entire phase space, including a Laguerre–Hermite pseudospectral formulation of velocity space, which allows for smooth interpolation between coarse gyrofluid-like resolutions and finer conventional gyrokinetic resolutions and efficient evaluation of a model collision operator. Additionally, we have built GX to natively target graphics processors (GPUs), which are among the fastest computational platforms available today. Finally, we have taken advantage of the reactor-relevant limit of small $\rho _*$ by using the radially local flux-tube approach. In this paper we present details about the gyrokinetic system and the numerical algorithms used in GX to solve the system. We then present several numerical benchmarks against established gyrokinetic codes in both tokamak and stellarator magnetic geometries to verify that GX correctly simulates gyrokinetic turbulence in the small $\rho _*$. Moreover, we show that the convergence properties of the Laguerre–Hermite spectral velocity formulation are quite favourable for nonlinear problems of interest. Coupled with GPU acceleration, which we also investigate with scaling studies, this enables GX to be able to produce useful turbulence simulations in minutes on one (or a few) GPUs and higher fidelity results in a few hours using several GPUs. GX is open-source software that is ready for fusion reactor design studies.
Advancement of solid-state electrolytes (SSEs) for all solid-state batteries typically focuses on modification of a known structural framework to improve conductivity, e.g., cation substitution for an immobile ion or varying the concentration of the mobile ions. Novel frameworks can be disruptive by enabling fast ion conduction aided by different structure and diffusion mechanisms, thereby unlocking optimal conductors with different properties. Herein, we perform a high-throughput survey of a structural framework for sodium ion conduction, Na8–x A x P2O9 (NAP), to understand the family’s thermodynamic stability, synthesizability, and ionic conduction. We show that the parent phase Na4TiP2O9 (NTP) undergoes a structural distortion (with accompanying conductivity transition) due to unstable phonons arising from pseudo-Jahn–Teller mode in the 1D titanium chains. Screening compounds in which Ti is substituted by other metals computationally reveal a number of candidates that are predicted to be low in formation energy and have high predicted ionic conductivities. High-throughput experimental and subsequent methodology optimization trials deliver one new compound, Na4SnP2O9 (NSP). X-ray diffraction (XRD), microscopy, and spectroscopy characterization indicate that the room-temperature structure of NSP is similar to the high-temperature, orthorhombic NTP phase but with some small unresolved structural differences. These uncharacterized structural details are speculated to limit the ion conductivity. Temperature-dependent XRD and electrochemical impedance spectroscopy indicate multiple coupled conductivity–structure transitions at a high temperature. We demonstrate the challenges with synthesis development and a priori identification of promising SSE phases as a major bottleneck in new (energy) materials development.
The Flexible Resource Scheduler is a hierarchical controller that manages the distributed energy resources in a distribution substation or distribution feeder to provide a firm commitment of power flow at the substation or feeder head to be scheduled in transmission-level markets as an aggregated demand resource. It is the reference controller for the FAST-DERMS Architecture, developed in tandem with the architecture under the DOE FAST-DERMS project. It is comprised of a day-ahead stochastic optimization, which schedules substation power flow and reserves, a intra-hour MPC, which generates dispatch base points for DER, and a real-time PID controller maintaining that dispatches DER to maintain the substation power around the base points. The repository also includes a representative aggregator controller, and all of the necessary components to run a simulation using PNNL's GridAPPS-D software with the controller.
The objective of the project is to design, fabricate, and optimize Reversible Metal Electrodeposition (RME) windows that simultaneously achieve color-neutral, fast, uniform switching over a wide dynamic range using low-cost, high-throughput materials and processes.
Here, a conceptual design of a scintillator-based fast-ion loss detector (FILD) has been developed for the Wendelstein 7-X stellarator (W7-X). Simulations using the Monte Carlo codes ASCOT5 and BEAMS3D have been performed to calculate the expected flux of neutral beam injection (NBI)-generated fast hydrogen ions onto the conceptual detector probe head. These fast-ion loss fluxes have been calculated for several magnetic field configurations as well as probe insertion positions. At the maximum insertion position, both co-and counter-going losses with high incident pitch angles are observed; however, at retracted positions, only co-going fast ions reach the probe head. The FILDSIM code has been used to optimize the geometry of the detector entrance and collimating elements to achieve a wide velocity space coverage as well as a high velocity–space resolution. A synthetic FILD signal is calculated for the expected loss distribution via forward modeling using the instrument response function. The synthetic signal is found to largely retain the velocity space features of the loss distribution.
High energy colliders provide a critical tool in nuclear physics study by probing the fundamental structure and dynamics of matter. To maximize the potential of scientific discovery in nuclear physics study, it is important to optimize the parameters of these colliders to attain the best performance. The performance of a collider is typically measured by its integrated luminosity of colliding beams since the probability of a new event is proportional to the integrated luminosity. However, the achievable luminosity is limited by the electromagnetic interactions (beam-beam effects) of two colliding beams at higher energy, and the interplay between the space-charge effects and the beam-beam effects at lower energy. To achieve the best performance of a collider means to attain the highest luminosity of the collider with optimized collider parameters. Optimizing the collider’s machine parameters is both computationally and experimentally expensive. A fast and robust computational framework including beam-beam and space-charge effects will be critical to attaining the best performance of the collider. In this project, we will study the beam dynamics challenges, specifically the interplay of the space-charge and the beam-beam effects, and the machine tuning models for maximizing the performance of RHIC experiments. We will develop an advanced modeling framework based on first-principles physical simulations, lattice models and the state-of-the-art machine learning methods and apply this framework to performance improvement of the RHIC in operation. We will build data manipulation packages to connect the simulation data and the experimental data with the framework, develop a self-consistent hybrid model of space-charge and beam-beam effects, study underlying physics mechanisms, build surrogate models using the labeled data, integrate the models into the advanced modeling framework, and apply the framework to RHIC luminosity (STAR and sPHENIX) optimization. The success of this project would substantially improve the performance of existing and future colliders and increase the opportunity for scientific discovery.
In the last 50 years there has been a tremendous progress in solving numerically the Navier-Stokes equations using finite differences, finite elements, spectral, and even meshless methods. Yet, in many real cases, we still cannot incorporate seamlessly (multi-fidelity) data into existing algorithms, and for industrial-complexity applications the mesh generation is time consuming and still an art. Moreover, solving ill-posed problems (e.g., lacking boundary conditions) or inverse problems is often prohibitively expensive and requires different formulations and new computer codes. Here, we employ physics-informed neural networks (PINNs), encoding the governing equations directly into the deep neural network via automatic differentiation, to overcome some of the aforementioned limitations for simulating incompressible laminar and turbulent flows. We develop the Navier-Stokes flow nets (NSFnets) by considering two different mathematical formulations of the Navier-Stokes equations: the velocity-pressure (VP) formulation and the vorticity-velocity (VV) formulation. Since this is a new approach, we first select some standard benchmark problems to assess the accuracy, convergence rate, computational cost and flexibility of NSFnets; analytical solutions and direct numerical simulation (DNS) databases provide proper initial and boundary conditions for the NSFnet simulations. The spatial and temporal coordinates are the inputs of the NSFnets, while the instantaneous velocity and pressure fields are the outputs for the VP-NSFnet, and the instantaneous velocity and vorticity fields are the outputs for the VV-NSFnet. This is unsupervised learning and, hence, no labeled data are required beyond boundary and initial conditions and the fluid properties. The residuals of the VP or VV governing equations, together with the initial and boundary conditions, are embedded into the loss function of the NSFnets. No data is provided for the pressure to the VP-NSFnet, which is a hidden state and is obtained via the incompressibility constraint without extra computational cost. Unlike the traditional numerical methods, NSFnets inherit the properties of neural networks (NNs), hence the total error is composed of the approximation, the optimization, and the generalization errors. Here, we empirically attempt to quantify these errors by varying the sampling (“residual”) points, the iterative solvers, and the size of the NN architecture. For the laminar flow solutions, we show that both the VP and the VV formulations are comparable in accuracy but their best performance corresponds to different NN architectures. The initial convergence rate is fast but the error eventually saturates to a plateau due to the dominance of the optimization error. For the turbulent channel flow, we show that NSFnets can sustain turbulence at , but due to expensive training we only consider part of the channel domain and enforce velocity boundary conditions on the subdomain boundaries provided by the DNS data base. We also perform a systematic study on the weights used in the loss function for balancing the data and physics components, and investigate a new way of computing the weights dynamically to accelerate training and enhance accuracy. In the last part, we demonstrate how NSFnets should be used in practice, namely for ill-posed problems with incomplete or noisy boundary conditions as well as for inverse problems. We obtain reasonably accurate solutions for such cases as well without the need to change the NSFnets and at the same computational cost as in the forward well-posed problems. As a result, we also present a simple example of transfer learning that will aid in accelerating the training of NSFnets for different parameter settings.