Search NASA⌕ Search

SEARCH · Search NASA

Results for “iterative thresholding”

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 163 records · Page 9

Calibration of coronary calcium scores determined using iterative image reconstruction (AIDR 3D) at 120, 100, and 80 kVp

Purpose: Computed tomography (CT) radiation dose reduction is frequently achieved by applying lower tube voltages and using iterative reconstruction (IR). For calcium scoring, the reference protocol at 120 kVp with filtered back projection (FBP) is still used, because kVp and IR may influence the Agatston score (AS) and volume score (VS). The authors present a two-step method to optimize dose: first, to determine the lowest feasible exposure and highest noise thresholds; second, to define a calibration method that ensures that the AS and VS are similar to the reference protocol. Methods: AS and VS were measured for an anthropomorphic thoracic phantom that includes a calcium calibration module. The phantom was scanned on a 320-row CT scanner, at tube voltages of 120 kVp using FBP, and 120, 100, and 80 kVp using adaptive iterative dose reduction (AIDR 3D) reconstruction. The minimum CTDIs were determined based on three objective quality criteria. Calibration was performed to estimate adjusted CT number thresholds for the lower kVp acquisitions. Finally, the accuracies of the total and individual insert scores at dose level close to the minimum CTDI level were investigated and compared to low (FBP{sub LD} − 120) and high (FBP{sub HD} − 120) dose reference protocols (based on ten repeated acquisitions for each group). Results: IR allows the exposure to be reduced by 69% at 120 kVp, with no significant effect on the total scores when averaged over all included dose steps and compared to FBP-120 (AS: 693 vs 699, p = 0.182; VS: 588 vs 587 mm{sup 3}, p = 0.569). Also when averaged over ten repeated scans and compared to FBP{sub HD} − 120 (AS: 709 vs 704, p = 0.435; VS: 604 vs 601 mm{sup 3}, p = 0.479), there is no statistical significant effect. Reducing the peak tube voltage allows even greater dose reductions: 73% at 100 kVp and 76% at 80 kVp. The calibrated CT number thresholds for analysis at 120, 100, and 80 kVp were, respectively, 130, 133, and 160 HU for the Agatston score, and 130, 132, and 140 HU for the volume score. Following the calibration, the mean scores of the four groups with dose variation were not significantly different from the reference protocol, at 100 kVp (AS: 698 vs 699, p = 0.818; VS: 584 vs 587 mm{sup 3}, p = 0.365) or at 80 kVp (AS: 698 vs 699, p = 0.996; VS: 586 vs 587 mm{sup 3}, p = 0.827). Similarly, there was no significant score difference with FBP{sub LD} − 120 during repeated scanning: 100 kVp (AS: 690 vs 694, p = 0.394; VS: 579 vs 585 mm{sup 3}, p = 0.168) and 80 kVp (AS: 703 vs 694, p = 0.115; VS: 588 vs 585 mm{sup 3}, p = 0.613). Compared to FBP{sub HD} − 120 group, the relative score deviation for the accuracy of the 400 and 800 mg/cm{sup 3} HA inserts with 3 and 5 mm diameter is less than 7%. However, the relative deviation of the smaller 1 mm inserts is poorer (up to 41% deviations for scores <3). Conclusions: With iterative reconstruction using AIDR 3D, deviations of the total Agatston and volume scores remain within 4% of the reference protocol. The 1 mm inserts were detected as calcification, but scores less than ten tend to be underestimated. Following the calibration process, the application of IR in combination with reduced tube voltages allows up to 76% lower radiation dose.

60 APPLIED LIFE SCIENCES↗

Highly Scalable Matching Pursuit Signal Decomposition Algorithm

Matching Pursuit Decomposition (MPD) is a powerful iterative algorithm for signal decomposition and feature extraction. MPD decomposes any signal into linear combinations of its dictionary elements or atoms . A best fit atom from an arbitrarily defined dictionary is determined through cross-correlation. The selected atom is subtracted from the signal and this procedure is repeated on the residual in the subsequent iterations until a stopping criterion is met. The reconstructed signal reveals the waveform structure of the original signal. However, a sufficiently large dictionary is required for an accurate reconstruction; this in return increases the computational burden of the algorithm, thus limiting its applicability and level of adoption. The purpose of this research is to improve the scalability and performance of the classical MPD algorithm. Correlation thresholds were defined to prune insignificant atoms from the dictionary. The Coarse-Fine Grids and Multiple Atom Extraction techniques were proposed to decrease the computational burden of the algorithm. The Coarse-Fine Grids method enabled the approximation and refinement of the parameters for the best fit atom. The ability to extract multiple atoms within a single iteration enhanced the effectiveness and efficiency of each iteration. These improvements were implemented to produce an improved Matching Pursuit Decomposition algorithm entitled MPD++. Disparate signal decomposition applications may require a particular emphasis of accuracy or computational efficiency. The prominence of the key signal features required for the proper signal classification dictates the level of accuracy necessary in the decomposition. The MPD++ algorithm may be easily adapted to accommodate the imposed requirements. Certain feature extraction applications may require rapid signal decomposition. The full potential of MPD++ may be utilized to produce incredible performance gains while extracting only slightly less energy than the standard algorithm. When the utmost accuracy must be achieved, the modified algorithm extracts atoms more conservatively but still exhibits computational gains over classical MPD. The MPD++ algorithm was demonstrated using an over-complete dictionary on real life data. Computational times were reduced by factors of 1.9 and 44 for the emphases of accuracy and performance, respectively. The modified algorithm extracted similar amounts of energy compared to classical MPD. The degree of the improvement in computational time depends on the complexity of the data, the initialization parameters, and the breadth of the dictionary. The results of the research confirm that the three modifications successfully improved the scalability and computational efficiency of the MPD algorithm. Correlation Thresholding decreased the time complexity by reducing the dictionary size. Multiple Atom Extraction also reduced the time complexity by decreasing the number of iterations required for a stopping criterion to be reached. The Course-Fine Grids technique enabled complicated atoms with numerous variable parameters to be effectively represented in the dictionary. Due to the nature of the three proposed modifications, they are capable of being stacked and have cumulative effects on the reduction of the time complexity.

Christensen, Daniel↗

Power and isotope effects in the ITER baseline scenario with tungsten and tungsten-equivalent radiators in DIII-D

Abstract Experiments in DIII-D document the ITER Baseline Scenario (IBS) at q 95 ∼ 3 and P IN / P LH ∼ 1–2, in both deuterium and hydrogen utilizing Kr and Xe as Tungsten-equivalent radiators. The power threshold for H-mode operation ( P LH ) was determined experimentally without added impurities and found to be about a factor of two higher than the scaling law. In recent IBS experiments in deuterium, intrinsic levels of metals such as Tungsten (W) or molybdenum and inconel are present that reduce the pedestal pressure by 20%–25%. A complete radiative collapse of deuterium IBS plasmas occurs at W core concentrations C W = 10 −5 . Simulations show that for core temperatures expected for ITER, the plasmas would not have a radiative collapse at C W = 1 × 10 −5 , moreover Q = 8–10 would still be achieved for C W up to 3 × 10 −5 . In contrast to deuterium, the IBS in hydrogen is not affected by intrinsic high-Z impurities, indicating that hydrogen H-modes in ITER may not inform the D-T phase with respect to W accumulation and discharge survival. Compared to deuterium, the pedestal pressure in hydrogen is ∼25% lower, with much higher ELM frequency of 150 Hz, decreasing with input power. Krypton was injected in a matrix scan of input power and impurity flow in IBS hydrogen discharges. Krypton impurity density profiles in hydrogen are similar to deuterium plasmas, but at Kr flows that are 2–3 times higher for the same input power. Krypton is transported into the core and affects the whole radius; at the highest injection rates a radiative collapse occurs at core radiation fractions of 0.3–0.35, consistent with the expected maximum W radiation fraction for ITER core plasmas. Comparing the results with previous International Tokamak Physics Activity database studies of the IBS confirms that at higher radiation fraction due to high-Z impurities, a drop in H 98 of >10% is observed. On the other hand, the results using Kr as a W-equivalent radiator indicate that metal (W) devices at lower core temperatures than ITER may provide overly pessimistic performance extrapolations to ITER for deuterium-tritium operation. The new DIII-D results support a more attractive option for the ITER Research Plan with a short hydrogen phase for system commissioning, transitioning to deuterium operations as soon as possible to provide relevant conditions for deuterium-tritium operations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluating iterative reconstruction performance in computed tomography

Purpose: Iterative reconstruction (IR) offers notable advantages in computed tomography (CT). However, its performance characterization is complicated by its potentially nonlinear behavior, impacting performance in terms of specific tasks. This study aimed to evaluate the performance of IR with both task-specific and task-generic strategies. Methods: The performance of IR in CT was mathematically assessed with an observer model that predicted the detection accuracy in terms of the detectability index (d′). d′ was calculated based on the properties of the image noise and resolution, the observer, and the detection task. The characterizations of image noise and resolution were extended to accommodate the nonlinearity of IR. A library of tasks was mathematically modeled at a range of sizes (radius 1–4 mm), contrast levels (10–100 HU), and edge profiles (sharp and soft). Unique d′ values were calculated for each task with respect to five radiation exposure levels (volume CT dose index, CTDI{sub vol}: 3.4–64.8 mGy) and four reconstruction algorithms (filtered backprojection reconstruction, FBP; iterative reconstruction in imaging space, IRIS; and sinogram affirmed iterative reconstruction with strengths of 3 and 5, SAFIRE3 and SAFIRE5; all provided by Siemens Healthcare, Forchheim, Germany). The d′ values were translated into the areas under the receiver operating characteristic curve (AUC) to represent human observer performance. For each task and reconstruction algorithm, a threshold dose was derived as the minimum dose required to achieve a threshold AUC of 0.9. A task-specific dose reduction potential of IR was calculated as the difference between the threshold doses for IR and FBP. A task-generic comparison was further made between IR and FBP in terms of the percent of all tasks yielding an AUC higher than the threshold. Results: IR required less dose than FBP to achieve the threshold AUC. In general, SAFIRE5 showed the most significant dose reduction potentials (11–54 mGy, 77%–84%), followed by SAFIRE3 (7–36 mGy, 50%–61%) and IRIS (6–26 mGy, 37%–50%). The dose reduction potentials highly depended on task size and task contrast, with tasks of lower contrasts and smaller sizes, i.e., more challenging tasks, indicating higher dose reductions. Softer edge profile showed higher dose reduction potentials with SAFIRE3 and SAFIRE5, but not with IRIS. The task-generic comparison between IR and FBP demonstrated the overall superiority of IR performance, as IR allowed a larger percent of tasks to exceed the threshold AUC: IRIS, 8%–12%; SAFIRE3, 10%–16%; and SAFIRE5, 20%–33%. The improvement with IR was generally more pronounced at lower dose levels. Conclusions: Expanding beyond traditional contrast and noise based assessments of IR, we performed both task-specific and task-generic evaluations of IR performance. The task-specific evaluation demonstrated the dependency of IR’s dose reduction potential on task attributes, which can be employed to optimize IR for clinical indications with specific range of size and contrast. The task-generic evaluation demonstrated IR’s overall superiority over FBP in terms of the range of tasks exceeding a threshold performance level, which can be employed for general comparisons between algorithms.

60 APPLIED LIFE SCIENCES↗

Higher Order Bases in a 2D Hybrid BEM/FEM Formulation

The advantages of using higher order, interpolatory basis functions are examined in the analysis of transverse electric (TE) plane wave scattering by homogeneous, dielectric cylinders. A boundary-element/finite-element (BEM/FEM) hybrid formulation is employed in which the interior dielectric region is modeled with the vector Helmholtz equation, and a radiation boundary condition is supplied by an Electric Field Integral Equation (EFIE). An efficient method of handling the singular self-term arising in the EFIE is presented. The iterative solution of the partially dense system of equations is obtained using the Quasi-Minimal Residual (QMR) algorithm with an Incomplete LU Threshold (ILUT) preconditioner. Numerical results are shown for the case of an incident wave impinging upon a square dielectric cylinder. The convergence of the solution is shown versus the number of unknowns as a function of the completeness order of the basis functions.

Fink, Patrick W.↗

Chevron beam dump for ITER edge Thomson scattering system

This paper contains the design of the beam dump for the ITER edge Thomson scattering system and mainly concerns its lifetime under the harsh thermal and electromagnetic loads as well as tight space allocation. The lifetime was estimated from the multi-pulse laser-induced damage threshold. In order to extend its lifetime, the structure of the beam dump was optimized. A number of bent sheets aligned parallel in the beam dump form a shape called a chevron which enables it to avoid the concentration of the incident laser pulse energy. The chevron beam dump is expected to withstand thermal loads due to nuclear heating, radiation from the plasma, and numerous incident laser pulses throughout the entire ITER project with a reasonable margin for the peak factor of the beam profile. Structural analysis was also carried out in case of electromagnetic loads during a disruption. Moreover, detailed issues for more accurate assessments of the beam dump's lifetime are clarified. Variation of the bi-directional reflection distribution function (BRDF) due to erosion by or contamination of neutral particles derived from the plasma is one of the most critical issues that needs to be resolved. In this paper, the BRDF was assumed, and the total amount of stray light and the absorbed laser energy profile on the beam dump were evaluated.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-E-T-111: An Iterative Approach to Modelling Dosimetric Properties of a Varian MLC

Purpose: In effort to efficiently and accurately model MLC parameters for use in Varian Eclipse treatment planning system (TPS), an iterative technique for characterizing the dosimetric properties of MLCs has been developed. Methods: The initial characterization of the dosimetric leaf gap (DLG) and MLC transmission was performed using varied sizes of moving gaps and transmission measurements through central axis of the radiation beam. These values were determined during the initial data collection. After the initial values were determined and entered into the TPS, two plans were calculated to test the veracity of DLG and transmission independently. The plan used to determine the DLG has openings of varied sizes which are affected differently by the DLG. A large field uniform fluence that required a MLC carriage shift was used to determine the transmission value independent of DLG. After characterizing the dosimetric properties of the MLC independently, clinical IMRT plans were tested and used to tweak the final values. Results: This approach to modelling MLC dosimetric properties was employed for: standard 120 leaf MLC (SDMLC) and 120 leaf high definition MLC (HDMLC). Both types were successfully commissioned, tested, and approved for clinical use on 4 SDMLC machines and 2 HDMLC machines. The pass rate of clinical IMRT plans was 96.9% for 6MV and 98.3% for 10MV for SDMLC and 96.3% for 6MV, 96.5% for 6FFF, and 97.5% for 10FFF for HDMLCs using MapCheck2 with gamma analysis of 2%/2mm and 10% threshold. The final DLG and transmission parameters for 6X were 0.5mm and 0.01, and 1.35mm and 0.02 for HD and SD MLCs respectively. Conclusion: This iterative approach of determining dosimetric properties of MLCs is viable for both SDMLC and HDMLC and multiple energies. This approach accurately and efficiently determines the values needed for commissioning MLCs when they are employed for modulated therapies.

60 APPLIED LIFE SCIENCES↗

Vapor shielding models and the energy absorbed by divertor targets during transient events

The erosion of divertor targets caused by high heat fluxes during transients is a serious threat to ITER operation, as it is going to be the main factor determining the divertor lifetime. Under the influence of extreme heat fluxes, the surface temperature of plasma facing components can reach some certain threshold, leading to an onset of intense material evaporation. The latter results in formation of cold dense vapor and secondary plasma cloud. This layer effectively absorbs the energy of the incident plasma flow, turning it into its own kinetic and internal energy and radiating it. This so called vapor shielding is a phenomenon that may help mitigating the erosion during transient events. In particular, the vapor shielding results in saturation of energy (per unit surface area) accumulated by the target during single pulse of heat load at some level E{sub max}. Matching this value is one of the possible tests to verify complicated numerical codes, developed to calculate the erosion rate during abnormal events in tokamaks. The paper presents three very different models of vapor shielding, demonstrating that E{sub max} depends strongly on the heat pulse duration, thermodynamic properties, and evaporation energy of the irradiated target material. While its dependence on the other shielding details such as radiation capabilities of material and dynamics of the vapor cloud is logarithmically weak. The reason for this is a strong (exponential) dependence of the target material evaporation rate, and therefore the “strength” of vapor shield on the target surface temperature. As a result, the influence of the vapor shielding phenomena details, such as radiation transport in the vapor cloud and evaporated material dynamics, on the E{sub max} is virtually completely masked by the strong dependence of the evaporation rate on the target surface temperature. However, the very same details define the amount of evaporated particles, needed to provide an effective shielding to the target, and, therefore, strongly influence resulting erosion rate. Thus, E{sub max} cannot be used for validation of shielding models and codes, aimed at the target material erosion calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transport at high ${\beta_p}$ and development of candidate steady state scenarios for ITER

On DIII-D, the high β p scenario has an internal transport barrier (ITB), β N ~β p ~3,q 95 ~10, and very high normalized confinement H 98,y2 ~1.6. Recently, plasmas starting with these conditions have been dynamically driven to q 95 ~6 and β p ~2, where we find the ITB and high performance persist for five energy confinement times. These conditions are projected to meet the ITER steady-state goal of Q = 5. The ITB is maintained at lower β p with a strong reverse shear, consistent with predictions that negative central shear can lower the β p threshold for the ITB. There are two observed confinement states in the high β p scenario: H-mode confinement state with a high edge pedestal, and an enhanced confinement state with a low pedestal and an ITB. It has been observed in a scan of external resonant magnetic perturbation amplitude that when there are no large type-I ELMs, there is no transition to enhanced confinement. This is consistent with the proposed mechanism for ITB formation being a type-I ELM. Quasilinear gyro-Landau fluid predictive modeling of ITER suggests that only a modest reverse shear is required to achieve the ITB formation necessary for Q=5 when electromagnetic physics including the kinetic ballooning mode (KBM) is incorporated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Method of up-front load balancing for local memory parallel processors

In a parallel processing computer system with multiple processing units and shared memory, a method is disclosed for uniformly balancing the aggregate computational load in, and utilizing minimal memory by, a network having identical computations to be executed at each connection therein. Read-only and read-write memory are subdivided into a plurality of process sets, which function like artificial processing units. Said plurality of process sets is iteratively merged and reduced to the number of processing units without exceeding the balance load. Said merger is based upon the value of a partition threshold, which is a measure of the memory utilization. The turnaround time and memory savings of the instant method are functions of the number of processing units available and the number of partitions into which the memory is subdivided. Typical results of the preferred embodiment yielded memory savings of from sixty to seventy five percent.

Baffes, Paul Thomas↗

Shatter Thresholds and Fragment Size Distributions of Deuterium–Neon Mixture Cryogenic Pellets for Tokamak Thermal Mitigation

Reliable mitigation is necessary to eliminate the detrimental effects of a disruption event in large high-current tokamaks such as ITER. To avoid serious damage to plasma-facing components during the thermal quench phase of a disruption, material is injected to radiate the plasma energy over the inner surface of the machine. The most promising method of material injection is a process known as shattered pellet injection (SPI). SPI utilizes cryogenic cooling to desublimate gas into the barrel of a pipe gun to form a solid pellet. High-pressure gas or a mechanical punch is used to dislodge the pellet and accelerate it into a bent tube to intentionally fracture it. Pellets made of a mixture of deuterium and neon are likely candidates for thermal mitigation. The survivability of these pellets throughout their flight path, before striking the shatter tube, is essential for reliable SPI operation. Experiments were conducted to determine intact speed limits for various mixtures. In this paper, we outline the details of brittle fracture theory and compare a theory-based model to experimental results from various mixtures of deuterium and neon pellets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Diagnostic of fusion neutrons on JET tokamak using diamond detector

In 2011-2012, an experimental campaign with a significant yield of fusion neutrons was carried out on the JET tokamak. During this campaign the facility was equipped with two diamond detectors based on natural and artificial CVD diamond. These detectors were designed and manufactured in State Research Center of Russian Federation TRINITI. The detectors measure the flux of fast neutrons with energies above 0.2 MeV. They have been installed in the torus hall and the distance from the center of plasma was about 3 m. For some of the JET pulses in this experiment, the neutron flux density corresponded to the operational conditions in collimator channels of ITER Vertical Neutron Camera. The main objective of diamond monitors was the measurement of total fast neutron flux at the detector location and the estimation of the JET total neutron yield. The detectors operate as threshold counters. Additionally a spectrometric measurement channel has been configured that allowed us to distinguish various energy components of the neutron spectrum. In this paper we describe the neutron signal measuring and calibration procedure of the diamond detector. Fluxes of DD and DT neutrons at the detector location were measured. It is shown that the signals of total neutron yield measured by the diamond detector correlate with signals measured by the main JET neutron diagnostic based on fission chambers with high accuracy. This experiment can be considered as a successful test of diamond detectors in ITER-like conditions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Modelling and experiment to stabilize disruptive tearing modes in the ITER baseline scenario in DIII-D

Abstract The achievement of high gain, stationary conditions in a tokamak scenario aimed at producing fusion energy in the ITER Project is crucial to the demonstration that this form of energy can be used in future reactors to provide cheap and clean energy globally. Disruptions are a challenge for the fusion energy field, in particular for the ‘ITER Baseline Scenario’ (IBS), as reproduced in the DIII-D tokamak. This work shows that a solution has been found for the m = 2/ n = 1 tearing modes that have consistently caused disruptions in the IBS: stable operation down to zero input torque was achieved by modifying the current density profile at the beginning of the pressure flattop and the ELM character later in the discharges, guided by previous results showing that the most likely cause of these instabilities is the current density profile. The coupling between sawteeth, n>2 modes and the 2/1 TMs is shown to not be statistically significant, nor the leading origin for the evolution towards instability. Ideal and resistive MHD modeling provide positive verification that a steeper ‘well’ in the region of the q = 2 rational surface leads to worse ideal stability, higher tearing index Δ’ and lower threshold Δ’ c for resistive instabilities, consistent with the experimental results. This provides confidence that the methods used in this work can be extrapolated to other devices and applied to avoid disruptions in ITER and pulsed fusion devices worldwide.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Disruption prediction with artificial intelligence techniques in tokamak plasmas

In nuclear fusion reactors, plasmas are heated to very high temperatures of more than 100 million kelvin and, in so-called tokamaks, they are confined by magnetic fields in the shape of a torus. Light nuclei, such as deuterium and tritium, undergo a fusion reaction that releases energy, making fusion a promising option for a sustainable and clean energy source. Tokamak plasmas, however, are prone to disruptions as a result of a sudden collapse of the system terminating the fusion reactions. As disruptions lead to an abrupt loss of confinement, they can cause irreversible damage to present-day fusion devices and are expected to have a more devastating effect in future devices. Disruptions expected in the next-generation tokamak, ITER, for example, could cause electromagnetic forces larger than the weight of an Airbus A380. Furthermore, the thermal loads in such an event could exceed the melting threshold of the most resistant state-of-the-art materials by more than an order of magnitude. In this work, to prevent disruptions or at least mitigate their detrimental effects, empirical models obtained with artificial intelligence methods, of which an overview is given here, are commonly employed to predict their occurrence—and ideally give enough time to introduce counteracting measures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

60 Co γ Irradiation of a CMS CROC v2 to 1.6 Grad

This paper reports the irradiation of the production version of the readout chip for the CMS HL-LHC pixel detector, CROC~v2, to 1.6~Grad using γ s from a 60 Co source at Sandia National Laboratory. The CROC~v2 was kept powered, configured, and maintained as close as possible to its nominal operating temperature of -20~ o C during irradiation. Internal chip voltages, pixel response, and diagnostic ring oscillators were continuously monitored, and thresholds were optimized frequently. The CROC~v2 performed according to specifications at 1.6~Grad. Where available, these results will be compared with those for this and earlier iterations of the RD53 readout chip irradiated to a few Grad with lower-energy X-ray sources.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Climate-Water Adaptation for Future US Electricity Infrastructure

Future climate-water conditions are anticipated to increase electricity demand, reduce transmission capacity, and limit power production. Yet, typical electricity capacity expansion planning does not consider climate-water constraints. We project four alternative U.S. power system configurations using an iterative modeling and data exchange platform that integrates climate-driven hydrological, thermal power plant, and capacity expansion models. Through a comparison with traditional modeling approaches, we show that this novel approach provides greater confidence in electricity capacity projections by incorporating feasibility checks that adjust infrastructure development to reach grid reliability thresholds under climate-water constraints. Initial projections without climate-water impacts on electricity generation show future power systems become less vulnerable, independent of climate-water adaptation, as economic drivers increase renewable and natural gas-based capacity, while water-intensive coal and nuclear plants retire. However, power systems may face reliability challenges without climate-water adaptation, revealing the significance of incorporating climate-water impacts into power system planning. Climate-adjusted (Iterative approach) projections require a 5.3-12.0% increase in national-level capacity, relative to Initial projections, leading to an additional $125-143 billion (5.0-7.0%) in infrastructure costs. Variable renewable and natural gas technologies account for nearly all the additional capacity and, together with regional trade-offs in electricity generation, enhance grid performance to reach reliability thresholds. These adaptation transitions also lower water use and emissions, contributing to climate change mitigation, and highlight the trade-offs and impacts of both near and long-term electricity generation planning decisions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Fast Active-Set Thresholding Method for Nonnegative Least Squares

Nonnegative Least Squares (NNLS) is a fundamental constrained optimization problem encountered in many applications such as image deblurring, signal processing, nonnegative matrix factorization, magnetic microscopy, and hyperspectral imaging. Active-set based methods are a common class of algorithms for solving NNLS which identify the optimal variable set of the NNLS solution. They do so by iteratively solving a series of unconstrained least squares problems, identifying which variables violate the nonnegativity constraints, and then swapping variables in/out of consideration until the optimal set of variables is found. Several variations improving upon this method exist in the literature. In this work, we propose an active-set swap heuristic which further improves upon existing active-set based methods for NNLS. Our optimizations are based upon adding multiple variables to the passive set within a threshold of the smallest gradient value and removing variables within a similar threshold of the closest boundary constraint. We leverage these optimizations to yield a Fast Active-Set Thresholding NNLS (FAST-NNLS) algorithm which significantly outperforms the existing state-of-the-art NNLS algorithms for a wide range of problems. Rigorous convergence guarantees are proven for the proposed method. We demonstrate the effectiveness of our proposed method on multiple synthetic datasets and two realworld text analysis applications. In doing so, we present the most comprehensive NNLS solver comparison in the literature to date.

Cobb, Benjamin [Georgia Institute of Technology]↗

Extending Helium Partial Pressure Measurement Technology to JET DTE2 and ITER

The detection limit for helium (He) partial pressure monitoring via the Penning discharge optical emission diagnostic, mainly used for tokamak divertor effluent gas analysis, is shown here to be possible for He concentrations down to 0.1% in predominantly deuterium effluents. This result from a dedicated laboratory study means that the technique can now be extended to intrinsically (non-injected) He produced as fusion reaction ash in deuterium-tritium experiments. The paper also examines threshold ionization mass spectroscopy as a potential backup to the optical technique, but finds that further development is needed to attain with plasma pulse-relevant response times. Both these studies are presented in the context of continuing development of plasma pulse-resolving, residual gas analysis for the upcoming JET deuterium-tritium campaign (DTE-2) and for ITER.

ITER↗