Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “and Operations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,657 records · Page 92

On the effectiveness of neural operators at zero-shot weather downscaling

Machine-learning (ML) methods have shown great potential for weather downscaling. These data-driven approaches provide a more efficient alternative for producing high-resolution weather datasets and forecasts compared to physics-based numerical simulations. Neural operators, which learn solution operators for a family of partial differential equations, have shown great success in scientific ML applications involving physics-driven datasets. Neural operators are grid-resolution-invariant and are often evaluated on higher grid resolutions than they are trained on, i.e., zero-shot super-resolution. Given their promising zero-shot super-resolution performance on dynamical systems emulation, we present a critical investigation of their zero-shot weather downscaling capabilities, which is when models are tasked with producing high-resolution outputs using higher upsampling factors than are seen during training. To this end, we create two realistic downscaling experiments with challenging upsampling factors (e.g., 8x and 15x) across data from different simulations: the European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5) and the Wind Integration National Dataset Toolkit. While neural operator-based downscaling models perform better than interpolation and a simple convolutional baseline, we show the surprising performance of an approach that combines a powerful transformer-based model with parameter-free interpolation at zero-shot weather downscaling. We find that this Swin-Transformer-based approach mostly outperforms models with neural operator layers in terms of average error metrics, whereas an Enhanced Super-Resolution Generative Adversarial Network-based approach is better than most models in terms of capturing the physics of the ground truth data. We suggest their use in future work as strong baselines.

17 WIND ENERGY↗

Conservative discontinuous Galerkin schemes for nonlinear Dougherty–Fokker–Planck collision operators

Here, we present a novel discontinuous Galerkin algorithm for the solution of a class of Fokker–Planck collision operators. These operators arise in many fields of physics, and our particular application is for kinetic plasma simulations. In particular, we focus on an operator often known as the ‘Lenard–Bernstein’ or ‘Dougherty’ operator. Several novel algorithmic innovations, based on the concept of weak equality, are reported. These weak equalities are used to define weak operators that compute primitive moments, and are also used to determine a reconstruction procedure that allows an efficient and accurate discretization of the diffusion term. We show that when two integrations by parts are used to construct the discrete weak form, and finite velocity-space extents are accounted for, a scheme that conserves density, momentum and energy exactly is obtained. One novel feature is that the requirements of momentum and energy conservation lead to unique formulas to compute primitive moments. Careful definition of discretized moments also ensure that energy is conserved in the piecewise linear case, even though the kinetic-energy term, $v^2$ is not included in the basis set used in the discretization. A series of benchmark problems is presented and shows that the scheme conserves momentum and energy to machine precision. Empirical evidence also indicates that entropy is a non-decreasing function. The collision terms are combined with the Vlasov equation to study collisional Landau damping and plasma heating via magnetic pumping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Collision operator for electron runaway in cold weakly ionised plasmas

In cold weakly ionised plasmas, the Dreicer generation mechanism can be non-diffusive as demonstrated in a recent study. By expanding the previous letter, we present the detailed description of a proper collision operator to precisely account for the non-diffusive electron kinetics. The operator appropriately combines the Fokker–Planck operator and Boltzmann operator where free-bound collision cross-sections are valid in the low-energy region. The proposed operator is envisaged to predict runaway electron generation in cold weakly ionised plasmas, particularly to design a runaway-free reactor tokamak start-up.

fusion plasma↗

Internal Photoemission Spectroscopy Measurements of Interfacial Energy Barriers in Operating TaN/Hf 0.5 Zr 0.5 O 2 /TaN Metal/Ferroelectric/Metal (MFM) Devices

The effect of the "waking" and subsequent "poling" operations on the electron barriers at both top and bottom electrode interfaces in operating ferroelectric hafnium zirconium oxide (Hf 0.5 Zr 0.5 O 2 , HZO) metal/ferroelectric/metal (MFM) devices are measured for the first time via internal photoemission (IPE) spectroscopy. Top and bottom (TaN/HZO and HZO/TaN) barriers for pristine devices were measured at 2.6 eV and 2.9 eV, respectively. The waking operation (10 kHz bipolar voltage cycling above the coercive field) increased the top barrier to 2.8 eV while leaving the bottom barrier essentially unchanged. Poling operations (application of a longer 10 msec unipolar pulse) were found to significantly decrease both top and bottom barriers. Here, the poling direction (polarity) had relatively little impact. The barrier for the top interface under positive (P↓) and negative poling (P↑) was 2.1 and 2.2 eV, respectively, while the bottom barrier was 2.3 eV for P↓ and 2.4 eV for P↑. All barrier heights remained unchanged after several months at room temperature. Potential physical mechanisms responsible for the changes observed under electrical operation are consistent with movement and/or creation of charged oxygen vacancy defects.

36 MATERIALS SCIENCE↗

The Molybdenum Oxide Interface Limits the High-Temperature Operational Stability of Unencapsulated Perovskite Solar Cells

Here, we report on the improved operational stability of unencapsulated perovskite solar cells (PSCs) aged in an ambient atmosphere at elevated temperatures (70 °C) for >1000 h under constant illumination and bias at 30–50% relative humidity. We identify a previously unseen interfacial degradation mechanism concerning the use of a MoO x interlayer, which was originally added to increase operational stability. Specifically, the hole-transport layer/MoO x interface buckles under illumination at 70 °C, which leads to delamination and rapid losses of short-circuit current density corresponding to an average $t_{80}$ of ~55 h. By judiciously evaluating various hole-transport layers, interlayers, and contacts, we find that replacing the MoO x with a VO x interlayer, regardless of the other components in the solar cell, alleviates this buckling issue due to its higher activation barrier toward crystallization, leading to significant gains in PSC operational stability. Unencapsulated devices aged in an ambient atmosphere with a VO x interlayer retain 71% of their initial PCE on average after constant illumination and bias at 70 °C for 1100 h ($t_{80}$ ~ 645 h). Currently, this is the highest temperature reported for the operational stability of unencapsulated n-i-p PSCs aged in air. Identification of a new facet of the complex degradation mechanisms in PSCs will allow for targeted acceleration testing to speed the deployment of low-cost, long-lasting electricity generation under realistic operating temperatures.

14 SOLAR ENERGY↗

Optimizing Desalination Operations for Energy Flexibility

Despite the value of energy optimization in desalination processes, modeling dynamic operations for monthly billing periods has remained a computational challenge. This work proposes a framework for energy flexibility optimization, which includes new modeling features for independent operation of parallel skids, start-up delays associated with chemical stabilization, the consideration of industrial energy tariff structures, and inclusion of hourly electrical carbon intensities. This is done using a modular and computationally efficient formulation that guarantees a globally optimal solution with standard optimization solvers. In this study, the approach is demonstrated in two distinct case studies: a seawater desalination plant in Santa Barbara, CA, and an indirect potable reuse facility in San Jose, CA. Trends predicted from the model are validated against operational facility measurements from a demand response shutdown event. Preliminary results show that optimizing energy flexibility can result in 18.51% monthly cost savings over energy efficiency-optimized operation. The value extracted from a facility-wide shutdown during peak electricity price hours is hampered by start-up delays in post-treatment chemical stabilization. In cases in which a facility does not have much excess capacity, using a flow equalization tank or operating over a wide recovery range may be cost-effective.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Universal programmable quantum circuit schemes to emulate an operator

Unlike fixed designs, programmable circuit designs support an infinite number of operators. The functionality of a programmable circuit can be altered by simply changing the angle values of the rotation gates in the circuit. Here, we present a new quantum circuit design technique resulting in two general programmable circuit schemes. The circuit schemes can be used to simulate any given operator by setting the angle values in the circuit. This provides a fixed circuit design whose angles are determined from the elements of the given matrix-which can be non-unitary-in an efficient way. We also give both the classical and quantum complexity analysis for these circuits and show that the circuits require a few classical computations. For the electronic structure simulation on a quantum computer, one has to perform the following steps: prepare the initial wave function of the system; present the evolution operator U=e{sup -iHt} for a given atomic and molecular Hamiltonian H in terms of quantum gates array and apply the phase estimation algorithm to find the energy eigenvalues. Thus, in the circuit model of quantum computing for quantum chemistry, a crucial step is presenting the evolution operator for the atomic and molecular Hamiltonians in terms of quantum gate arrays. Since the presented circuit designs are independent from the matrix decomposition techniques and the global optimization processes used to find quantum circuits for a given operator, high accuracy simulations can be done for the unitary propagators of molecular Hamiltonians on quantum computers. As an example, we show how to build the circuit design for the hydrogen molecule.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characteristic operator functions for quantum input-plant-output models and coherent control

We introduce the characteristic operator as the generalization of the usual concept of a transfer function of linear input-plant-output systems to arbitrary quantum nonlinear Markovian input-output models. This is intended as a tool in the characterization of quantum feedback control systems that fits in with the general theory of networks. The definition exploits the linearity of noise differentials in both the plant Heisenberg equations of motion and the differential form of the input-output relations. Mathematically, the characteristic operator is a matrix of dimension equal to the number of outputs times the number of inputs (which must coincide), but with entries that are operators of the plant system. In this sense, the characteristic operator retains details of the effective plant dynamical structure and is an essentially quantum object. We illustrate the relevance to model reduction and simplification definition by showing that the convergence of the characteristic operator in adiabatic elimination limit models requires the same conditions and assumptions appearing in the work on limit quantum stochastic differential theorems of Bouten and Silberfarb [Commun. Math. Phys. 283, 491-505 (2008)]. This approach also shows in a natural way that the limit coefficients of the quantum stochastic differential equations in adiabatic elimination problems arise algebraically as Schur complements and amounts to a model reduction where the fast degrees of freedom are decoupled from the slow ones and eliminated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The quiescent H-mode regime for high performance edge localized mode-stable operation in future burning plasmas

For the first time, DIII-D experiments have achieved stationary quiescent H-mode (QH-mode) operation for many energy confinement times at simultaneous ITER-relevant values of beta, confinement, and safety factor, in an ITER-like shape. QH-mode provides excellent energy confinement, even at very low plasma rotation, while operating without edge localized modes (ELMs) and with strong impurity transport via the benign edge harmonic oscillation (EHO). By tailoring the plasma shape to improve the edge stability, the QH-mode operating space has also been extended to densities exceeding 80% of the Greenwald limit, overcoming the long-standing low-density limit of QH-mode operation. In the theory, the density range over which the plasma encounters the kink-peeling boundary widens as the plasma cross-section shaping is increased, thus increasing the QH-mode density threshold. The DIII-D results are in excellent agreement with these predictions, and nonlinear magnetohydrodynamic analysis of reconstructed QH-mode equilibria shows unstable low n kink-peeling modes growing to a saturated level, consistent with the theoretical picture of the EHO. Furthermore, high density operation in the QH-mode regime has opened a path to a new, previously predicted region of parameter space, named “Super H-mode” because it is characterized by very high pedestals that can be more than a factor of two above the peeling-ballooning stability limit for similar ELMing H-mode discharges at the same density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Operational Modes of a 2.0 MWth Chloride Molten-Salt Pilot-Scale System

The limit of traditional solar-salt thermal stability is around 600 °C with ambient air as the cover gas. Nitrate molten salt concentrating solar power (CSP) systems are currently deployed globally and are considered to be state-of the art heat transfer fluids (HTFs) for present day high-temperature operation. However, decomposition challenges occur with these salts for operation beyond 600 C. Although slightly higher limits may be possible with solar salt, to fully realize SunShot efficiency goals of $15/kWhth HTFs and an LCOE of 6¢/kWh, molten-salt technologies working at higher temperatures (e.g., 650 °C to 750 °C) will require an alternative salt chemistry composition, such as chlorides. In this investigation a 2.0MWth Pilot-scale CSP plant design is developed to assess thermodynamic performance potential for operation up to 720 C. Here, an Engineering Equation Solver (EES) model is developed with respect to 14 state-points from the base of a solar tower at the Sandia National Laboratories, National Solar Thermal Test Facility (NSTTF), to solar receiver mounted 120 ft. above the ground. The system design considers a ternary chloride ternary chloride (20%NaCl/40%MgCl/40%KCl by mol%) salt as the HTF, with 6 hrs. of storage and a 1 MWth primary salt to sCO2 heat exchanger. Preliminary system modelling results indicate a minimum non-dimensional Cv of 60 required for both cold and hot-side throttle recirculation valves for the operational pump operating between speeds of 1800 and 2400 RPM. Further receiver comparison study results suggest that the ternary salt requires an average 15.2% higher receiver flux with a slightly lower calculated receiver efficiency when compared to a binary carnelite salt to achieve a 2.0 MWth desired input power design.

associated liquids↗

Thin THz QCL active regions for improved continuous-wave operating temperature

We compare the performance of 10 and 5 μm thick metal–metal waveguide terahertz quantum-cascade laser ridges operating around 2.7 THz and based on a 4-well phonon depopulation active region design. Thanks to reduced heat dissipation and lower thermal resistance, the 5 μm thick material shows an 18 K increase in continuous wave operating temperature compared to the 10 μm material, despite a lower maximum pulsed-mode operating temperature and a larger input power density. A maximum continuous wave operating temperature of 129 K is achieved using the 5 μm thick material and a 15 μm wide ridge waveguide, which lased up to 155 K in the pulsed mode. The use of thin active regions is likely to become increasingly important to address the increasing input power density of emerging 2- and 3-well active region designs that show the highest pulsed operating temperatures.

47 OTHER INSTRUMENTATION↗

3D vacuum magnetic field modelling of the ITER ELM control coil during standard operating scenarios

In-vessel, non-axisymmetric, control coils have proven to be an important option for mitigating and suppressing edge-localized modes (ELMs) in high performance operating regimes on a growing number of tokamaks. Additionally, an in-vessel non-axisymmetric ELM control coil is being considered in the ITER baseline design. For this paper, in preparing for the initial operation of this coil set, a comprehensive study was carried out to characterize the linear superposition of the 3D vacuum magnetic field, produced by the ELM coil, on a series of equilibria representing nine standard ITER operating scenarios. Here, the spatial phase angle of toroidally distributed currents, specified with a cosine waveform, in the upper and lower rows of the ITER ELM coil (IEC) set is varied in 2° steps while holding the current in the equatorial row of coils constant. The peak current in each of the three toroidal rows of window-frame coils making up the IEC is scanned between 5 kAt and 90 kAt in 5 kAt steps and the width of the edge region covered by overlapping vacuum field magnetic islands is calculated. This width is compared to a vacuum field ELM suppression correlation criterion found in DIII-D. A minimum coil current satisfying the DIII-D criterion, along with an associated set of phase angles, is identified for each ITER operating scenario. These currents range from 20 kAt to 75 kAt depending on the operating scenario being used and the toroidal mode number (n) of the cosine waveform. Comparisons between the scaling of the divertor footprint area in cases with n = 3 perturbation fields versus those with n = 4 show significant advantages when using n = 3. In addition, it is found that the DIII-D correlation criterion can be satisfied in the event that various combinations of individual IEC window-frame coils need to be turned off due to malfunctioning components located inside the vacuum vessel. Finally, details of these results for both the full set of 27 window-frame coils and various reduced sets, using either n = 3 and n = 4 perturbation fields, are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A dimensionality reduction algorithm for mapping tokamak operational regimes using a variational autoencoder (VAE) neural network

Abstract A variational autoencoder (VAE) is a type of unsupervised neural network which is able to learn meaningful data representations in a reduced dimensional space. We present an application of VAE in identifying the operational stability boundary of tokamak plasma discharges. This model was implemented using a dataset of over 3000 discharges from the high beta tokamak-extended pulse (HBT-EP) device. We found the VAE model to be capable of forming a continuous low-dimensional operational space map and identifying the operational boundaries using a specified warning time window. By projecting the operational parameters onto the same reduced space, this provides an intuitive way for the machine operator or an automated control system to perform disruption avoidance using a relevant control actuator as a discharge approaches a boundary. Pre-programmed GPU control experiments were conducted to demonstrate this control technique using HBT-EP’s saddle control coils as a horizontal position actuator, showing the ability to avoid the oncoming disruptive event and extend the duration of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma–surface interaction in the stellarator W7-X: conclusions drawn from operation with graphite plasma-facing components

W7-X completed its plasma operation in hydrogen with island divertor and inertially cooled test divertor unit (TDU) made of graphite. A substantial set of plasma-facing components (PFCs), including in particular marker target elements, were extracted from the W7-X vessel and analysed post-mortem. The analysis provided key information about underlying plasma–surface interactions (PSI) processes, namely erosion, transport, and deposition as well as fuel retention in the graphite components. The net carbon (C) erosion and deposition distribution on the horizontal target (HT) and vertical target (VT) plates were quantified and related to the plasma time in standard divertor configuration with edge transform ι = 5/5, the dominant magnetic configuration of the two operational phases (OP) with TDU. The operation resulted in integrated high net C erosion rate of 2.8 mg s -1 in OP1.2B over 4809 plasma seconds. Boronisations reduced the net erosion on the HT by about a factor 5.4 with respect to OP1.2A owing to the suppression of oxygen (O). In the case of the VT, high peak net C erosion of 11μm at the strike line was measured during OP1.2B which converts to 2.5 nm s -1 or 1.4 mg s -1 when related to the exposed area of the target plate and the operational time in standard divertor configuration. PSI modelling with ERO2.0 and WallDYN-3D is applied in an interpretative manner and reproduces the net C erosion and deposition pattern at the target plates determined by different post-mortem analysis techniques. This includes also the 13 C tracer deposition from the last experiment of OP1.2B with local 13 CH 4 injection through a magnetic island in one half module. The experimental findings are used to predict the C erosion, transport, and deposition in the next campaigns aiming in long-pulse operation up to 1800 s and utilising the actively cooled carbon-fibre composite (CFC) divertor currently being installed. The CFC divertor has the same geometrical design as the TDU and extrapolation depends mainly on the applied plasma boundary. Extrapolation from campaign averaged information obtained in OP1.2B reveals a net erosion of 7.6 g per 1800 s for a typical W7-X attached divertor plasma in hydrogen.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of Wendelstein 7-X high-performance operation

The Wendelstein 7-X (W7-X) stellarator has completed two consecutive experimental campaigns OP 2.2 (Sep.-Dec. 2024) and OP 2.3 (Feb.-May 2025) under a new operational strategy enabling more than one year of uninterrupted device availability. This approach, supported by exceptionally high subsystem reliability, allowed sustained high-efficiency plasma operations with up to 80–100 discharges per day across a broad range of magnetic configurations. Several key technical upgrades-most notably the first operation of a 1.5 MW class steady-state gyrotron, a new steady-state pellet injector, and advanced real-time feedback control systems significantly enhanced heating, fueling, and plasma control capabilities. Together, these improvements enabled major advances in long-pulse performance, high-β operation, and confinement optimization. Long-pulse discharges achieved 1.8 GJ of injected energy under fully detached divertor conditions, while reduced-field scenarios facilitated record volume-averaged β values approaching 3%. High-performance plasmas with centrally peaked density profiles, created via neutral beam injection (NBI) or sustained pellet fueling, demonstrated strongly reduced turbulent transport and stellarator-record fusion triple products. Complementary studies of power exhaust and divertor heat loads revealed the role of scrape-off-layer drift physics in shaping strike-line patterns under attached conditions. Together, the results from OP 2.2 and OP 2.3 significantly expand the operational space of W7-X and strengthen its role as a leading platform for steady-state stellarator research and reactor-relevant plasma scenarios.

Wendelstein 7-X↗

Uncertainty quantification in the analyses of operational wind power plant performance

In the present work, we examine the variation introduced in the evaluation of an operating plant's wind power production as a result of the choices analysts make in the processing of the operational data. For this study, an idealized power production for individual turbines over an operational period was predicted by fitting power curves to the turbine production data collected during expected operation (that is, without curtailment or availability losses). A set of 240 possible methods were developed for (a) defining what data represented expected operation and (b) modeling the power curve. The spread in the idealized power production as predicted by the different methods was on average almost 3% for the 100 turbines considered. Such significant variation places a lower bound on the precision with which analysts may employ such data as benchmarks for calibration of their energy estimation processes and limits the potential for identification of refinements to the energy estimation models for improved accuracy.

17 WIND ENERGY↗

Progress towards 28 GHz operations at the Facility for Rare Isotope Beams (FRIB)

This paper provides an overview on the design, operational parameters as well as recent development and recent performance, of the superconducting electron cyclotron resonance (ECR) ion source at the Facility for Rare Isotope Beams (FRIB). The ion source has been used in support of beam operations since the end of 2022 and was designed to operate at 28 GHz and be capable of delivering high intensity uranium beams to help support beam delivery of up to a final beam power of 400 kW. Currently operating at up to 20 kW many primary beams have been developed and routinely used with the superconducting ECR including uranium. Initially operated at 18 GHz, the ion source has now been coupled to a 28 GHz gyrotron and high intensity beams have been demonstrated ahead of the facility planned power ramp up. text with the article abstract. Replace this text with the article abstract.

Machicoane, G. A. [Michigan State University, East↗

Neutral gas analysis for JET DT operation

Neutral gas analysis, the measurement and evaluation of total and partial pressures, is a key technique to study the impact of neutral gas dynamics on retention, recycling and release processes of fuel or impurity species in fusion devices. At JET, the experiment closest to ITER in terms of operating parameters and size, various detectors and techniques for partial pressure and total pressure measurements are deployed together to characterise neutral gas dynamics during and after plasma operation on various toroidal and poloidal locations. An extensive modification of JET's sub-divertor neutral gas diagnostic system aims at retaining and extending established measurement capabilities in the forthcoming Deuterium-Tritium (DT) experiments (DTE2). To achieve DT compatibility, a separation of radiation-sensitive electronics from the sensor and adequate radiation shielding is required, as well as utilisation of a DT compatible differential pumping system with adjustable throughput to account for the strong pressure variation in the sub-divertor region. Finally, the sub-divertor neutral gas diagnostic will be equipped with multiple Residual Gas Analysers (RGAs), utilising quadrupole mass spectrometry and electrostatic ion-trap-principles, all operating with remote electronics located behind the biological radiation shield. These RGAs will record data in a fast selected discrete mass mode during plasma pulses (cycle ~2 s) and will automatically switch back to continuous data recording (cycle ~100 s) afterwards. They will be complemented by a newly improved Penning gauge spectroscopy configuration in particular supporting the He and D 2 separation relevant for DT operation. The distance between these devices and their associated control unit is typically 15m . A newly developed RGA with a cable length of 80 m, compatible with the ITER environment, will also be employed for the first time. Finally, this set-up and its operation in DTE2 will provide vital input to the development of the ITER divertor RGA in the most relevant environment currently available.

47 OTHER INSTRUMENTATION↗