Search NASA⌕ Search

SEARCH · Search NASA

Results for “ITERATION”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Volumetric quantification of lung nodules in CT with iterative reconstruction (ASiR and MBIR)

Purpose: Volume quantifications of lung nodules with multidetector computed tomography (CT) images provide useful information for monitoring nodule developments. The accuracy and precision of the volume quantification, however, can be impacted by imaging and reconstruction parameters. This study aimed to investigate the impact of iterative reconstruction algorithms on the accuracy and precision of volume quantification with dose and slice thickness as additional variables.Methods: Repeated CT images were acquired from an anthropomorphic chest phantom with synthetic nodules (9.5 and 4.8 mm) at six dose levels, and reconstructed with three reconstruction algorithms [filtered backprojection (FBP), adaptive statistical iterative reconstruction (ASiR), and model based iterative reconstruction (MBIR)] into three slice thicknesses. The nodule volumes were measured with two clinical software (A: Lung VCAR, B: iNtuition), and analyzed for accuracy and precision.Results: Precision was found to be generally comparable between FBP and iterative reconstruction with no statistically significant difference noted for different dose levels, slice thickness, and segmentation software. Accuracy was found to be more variable. For large nodules, the accuracy was significantly different between ASiR and FBP for all slice thicknesses with both software, and significantly different between MBIR and FBP for 0.625 mm slice thickness with Software A and for all slice thicknesses with Software B. For small nodules, the accuracy was more similar between FBP and iterative reconstruction, with the exception of ASIR vs FBP at 1.25 mm with Software A and MBIR vs FBP at 0.625 mm with Software A.Conclusions: The systematic difference between the accuracy of FBP and iterative reconstructions highlights the importance of extending current segmentation software to accommodate the image characteristics of iterative reconstructions. In addition, a calibration process may help reduce the dependency of accuracy on reconstruction algorithms, such that volumes quantified from scans of different reconstruction algorithms can be compared. The little difference found between the precision of FBP and iterative reconstructions could be a result of both iterative reconstruction's diminished noise reduction at the edge of the nodules as well as the loss of resolution at high noise levels with iterative reconstruction. The findings do not rule out potential advantage of IR that might be evident in a study that uses a larger number of nodules or repeated scans.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2017

Neutronics Analysis for ITER Diagnostic Generic Upper Port Plug

This paper presents new results of neutronics analysis performed for the updated design of the ITER Diagnostic Generic Upper Port Plug (DGUPP). Such analysis is performed regularly during the development of the ITER Ports design. Neutronics analysis is an important step in support of the development of ITER components irradiated by 14 MeV neutrons generated in D-T fusion thermonuclear reactions. The neutrons are born in the ITER tokamak plasma chamber and then impinge on the surrounding blanket and port plugs, after that they diffuse through shielding material of vacuum vessel components or could stream through the diagnostic channels and assembly gaps far away from plasma up to the bio-shield and in a small amount beyond it. While ITER machine is in operation it is inaccessible, after the shut-down it is planned to provide personnel access to the Upper Port Inter-Space Structure (UP ISS) behind the Port Plug via the corridors. Neutrons activate the ITER materials and those radioactive materials emit decay gammas during the shut-down causing a Shut-Down Dose Rate (SDDR) field distributed all around the activated materials. Decay gammas could also stream through the void channels. That emphasizes the need for calculation of decay gamma transport in evaluations of SDDR, which is reflected in the use of transport - activation interface code described in Section of Computational Modelling. The assessment of radiation environment (fluxes and SDDR) at DGUPP was focused on the port back-side area of UP ISS, where personnel access is planned along the ISS maintenance corridors. Minimization of SDDR at the front part of the ISS maintenance corridors was the target value for optimization of the DGUPP shielding and analysis of related neutronics effects. Following the dose ALARA (As Low As Reasonably Achievable) principle, every possibility to reduce SDDR should be investigated, emphasizing the importance of this analysis. The challenges of such analysis are inextricably connected with the engineering complexity of the ITER machine, with integration of many diagnostics systems having channels and gaps inside the ITER ports. Those gaps result in radiation streaming and negatively affect the shielding performance of the port plugs. The computational complication is caused by the need to model precisely the ITER geometry and then to run deep penetrating radiation transport in the substantially heterogenic toroidal geometry of ITER large-size tokamak, with 30 m height and diameter of its bio-shield. The whole mass of Upper Port Plug (UPP) could be up to 25 tons and it must be very precisely positioned inside the port opening in blanket and Vacuum Vessel. Facing the engineering challenges for the assembly of DGUPP, there were proposed two Project Change Requests (PCRs) No. 439 and 709 concerning the design updates of UPP/blanket interface and its remote handling.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Methodological Aspects of Shielding and Activation Analyses for ITER Diagnostics Systems - Paper 125

It is expected that the diagnostics systems of the ITER tokamak will provide valuable measurement data about the first in its kind D-T fusion plasma with emission of 14 MeV neutrons. These systems must withstand the neutron and secondary photon radiation along the ITER operation. Under the neutron irradiation the ITER materials are getting radioactive and therefore an adequate shielding must be provided along all the time of ITER operation and shut-down for the reliable machine operation and for the safe personnel access to its components, such as the diagnostics ports. These are the prerequisites for the shielding and activation analysis of the ITER diagnostics systems discussed in this paper. This paper provides insight into the current tasks of neutronics support for the development of ITER diagnostics systems. The tasks themselves range from relatively straightforward ones such as Monte Carlo radiation transport calculations for shielding analyses of the diagnostics systems, to more advanced tasks of Shut-Down Dose Rate (SDDR) analyses for which the dedicated Rigorous 2-Step mesh-based (R2Smesh) and Direct 1-Step (D1S) methodologies were applied. The accuracy and precision of the ITER neutronics analysis is dependent on the methodological approaches and the implying modeling assumptions. In addition to the SDDR challenges, another distinguishing feature of fusion neutronics is a wide use of CAD-based neutronics modeling to design components of fusion devices such as ITER and DEMO. Basically, there are two possibilities of CAD geometry settings in Monte Carlo (MC) radiation transport models, either by the geometry format conversion of CAD models into the format of MC models by means of the interface programs such as McCad and MCAM or by the direct use of CAD geometry engine in MC calculations such as it was realized in the DAGMC code approach. Several conversion interfaces and the DAGMC software library have been tested on the 3D ITER model developed for the CAD-MCNP benchmark. It was concluded that the considered approaches have reached the maturity level to allow their application to real ITER design analyses.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

ITER Central Solenoid Module Fabrication

The fabrication of the modules for the ITER Central Solenoid (CS) has started in a dedicated production facility located in Poway, California, USA. The necessary tools have been designed, built, installed, and tested in the facility to enable the start of production. The current schedule has first module fabrication completed in 2017, followed by testing and subsequent shipment to ITER. The Central Solenoid is a key component of the ITER tokamak providing the inductive voltage to initiate and sustain the plasma current and to position and shape the plasma. The design of the CS has been a collaborative effort between the US ITER Project Office (US ITER), the international ITER Organization (IO) and General Atomics (GA). GA’s responsibility includes: completing the fabrication design, developing and qualifying the fabrication processes and tools, and then completing the fabrication of the seven 110 tonne CS modules. The modules will be shipped separately to the ITER site, and then stacked and aligned in the Assembly Hall prior to insertion in the core of the ITER tokamak. A dedicated facility in Poway, California, USA has been established by GA to complete the fabrication of the seven modules. Infrastructure improvements included thick reinforced concrete floors, a diesel generator for backup power, along with, cranes for moving the tooling within the facility. The fabrication process for a single module requires approximately 22 months followed by five months of testing, which includes preliminary electrical testing followed by high current (48.5 kA) tests at 4.7K. The production of the seven modules is completed in a parallel fashion through ten process stations. The process stations have been designed and built with most stations having completed testing and qualification for carrying out the required fabrication processes. The final qualification step for each process station is achieved by the successful production of a prototype coil. Fabrication of the first ITER module is in progress. The seven modules will be individually shipped to Cadarache, France upon their completion. This paper describes the processes and status of the fabrication of the CS Modules for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Final report on the MAGARC software consulting support contract (WFO No. 15904) for ITER magnet safety assessment

he ITER (not an acronym but Latin for "the Way”) International Organization (IO) is building a large-scale scientific experiment that is now under construction in Cadarache, France. The ITER Tokamak fusion reactor aims to demonstrate that it is possible to produce energy from fusion. The ITER’s power output goal is to deliver ten times the power it consumes. During its operational lifetime, ITER will test key technologies necessary for the next step: the demonstration fusion power plant that will prove that it is possible to capture fusion energy for commercial use. General information on the ITER project can be found at http://www.iter.org. The ITER reactor’s superconducting magnet systems create ~41GJ of stored magnetic field energy when fully energized. The consequences of such energy dissipating via electrical arcs under faulted conditions gives rise to safety concerns even though these fault events are extremely unlikely. As part of the ITER licensing procedure, and as a result of such safety concerns, this contract was established to benchmark a developing magnet safety analysis tool at Culham Centre for Fusion Energy against the MAGARC computer code. MAGARC was developed by the U.S. Idaho National Laboratory (INL) under the ITER Engineering Design Activity (EDA) Project and was applied to produce previous analyses, in order to support magnet fault and safety assessment for the ITER design description document. This contract report documents the completion of requested modifications to the MAGARC code by the ITER IO, a time-step size sensitivity analysis with MAGARC during an unmitigated quench event, a comparison of the modeling approaches taken in developing MAGARC and the CCFE analysis package, and the data required for the benchmarking activity for these codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Reducing the effects of acoustic heterogeneity with an iterative reconstruction method from experimental data in microwave induced thermoacoustic tomography

Purpose: An iterative reconstruction method has been previously reported by the authors of this paper. However, the iterative reconstruction method was demonstrated by solely using the numerical simulations. It is essential to apply the iterative reconstruction method to practice conditions. The objective of this work is to validate the capability of the iterative reconstruction method for reducing the effects of acoustic heterogeneity with the experimental data in microwave induced thermoacoustic tomography. Methods: Most existing reconstruction methods need to combine the ultrasonic measurement technology to quantitatively measure the velocity distribution of heterogeneity, which increases the system complexity. Different to existing reconstruction methods, the iterative reconstruction method combines time reversal mirror technique, fast marching method, and simultaneous algebraic reconstruction technique to iteratively estimate the velocity distribution of heterogeneous tissue by solely using the measured data. Then, the estimated velocity distribution is used subsequently to reconstruct the highly accurate image of microwave absorption distribution. Experiments that a target placed in an acoustic heterogeneous environment are performed to validate the iterative reconstruction method. Results: By using the estimated velocity distribution, the target in an acoustic heterogeneous environment can be reconstructed with better shape and higher image contrast than targets that are reconstructed with a homogeneous velocity distribution. Conclusions: The distortions caused by the acoustic heterogeneity can be efficiently corrected by utilizing the velocity distribution estimated by the iterative reconstruction method. The advantage of the iterative reconstruction method over the existing correction methods is that it is successful in improving the quality of the image of microwave absorption distribution without increasing the system complexity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Iterative metal artifact reduction for x-ray computed tomography using unmatched projector/backprojector pairs

Purpose: Metal artifact reduction (MAR) is a major problem and a challenging issue in x-ray computed tomography (CT) examinations. Iterative reconstruction from sinograms unaffected by metals shows promising potential in detail recovery. This reconstruction has been the subject of much research in recent years. However, conventional iterative reconstruction methods easily introduce new artifacts around metal implants because of incomplete data reconstruction and inconsistencies in practical data acquisition. Hence, this work aims at developing a method to suppress newly introduced artifacts and improve the image quality around metal implants for the iterative MAR scheme. Methods: The proposed method consists of two steps based on the general iterative MAR framework. An uncorrected image is initially reconstructed, and the corresponding metal trace is obtained. The iterative reconstruction method is then used to reconstruct images from the unaffected sinogram. In the reconstruction step of this work, an iterative strategy utilizing unmatched projector/backprojector pairs is used. A ramp filter is introduced into the back-projection procedure to restrain the inconsistency components in low frequencies and generate more reliable images of the regions around metals. Furthermore, a constrained total variation (TV) minimization model is also incorporated to enhance efficiency. The proposed strategy is implemented based on an iterative FBP and an alternating direction minimization (ADM) scheme, respectively. The developed algorithms are referred to as “iFBP-TV” and “TV-FADM,” respectively. Two projection-completion-based MAR methods and three iterative MAR methods are performed simultaneously for comparison. Results: The proposed method performs reasonably on both simulation and real CT-scanned datasets. This approach could reduce streak metal artifacts effectively and avoid the mentioned effects in the vicinity of the metals. The improvements are evaluated by inspecting regions of interest and by comparing the root-mean-square errors, normalized mean absolute distance, and universal quality index metrics of the images. Both iFBP-TV and TV-FADM methods outperform other counterparts in all cases. Unlike the conventional iterative methods, the proposed strategy utilizing unmatched projector/backprojector pairs shows excellent performance in detail preservation and prevention of the introduction of new artifacts. Conclusions: Qualitative and quantitative evaluations of experimental results indicate that the developed method outperforms classical MAR algorithms in suppressing streak artifacts and preserving the edge structural information of the object. In particular, structures lying close to metals can be gradually recovered because of the reduction of artifacts caused by inconsistency effects.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

ITER Neutronic Challenges for Upper Port 14

The future of magnetically confined nuclear fusion facilities is extremely promising, yet there are extensive challenges that still have to be addressed. These challenges range significantly, but one of the more difficult challenges is nuclear shielding and nuclear safety. The challenge of adequately shielding neutrons is especially prevalent at ITER. Throughout the lifetime of ITER there will be large quantities of high energy (14.1- MeV) neutrons generated, with the expected neutron lifetime fluence to be around 0.3 MW.a/m{sup 2}. For such a large neutron fluence, shielding is a very important factor for the operation of ITER. Not only the shielding of the people that will be supporting the facility, but also to the components that surround the plasma. Without the many diagnostics that are planned for ITER operation, ITER would not be able to function. But with the inclusion of the diagnostics and the many apertures required, an interesting neutronics environment is created. This environment is like nothing we have seen or are currently able to create on the earth. Due to this unknown, significant time and effort is being placed on the research, design, and analysis of the ITER neutronics environment. A neutron has an interesting life after birth within the plasma on ITER. There are many different paths and reactions that can take place; so many, that trying to simplify the problem has brought significant unknowns into the analysis. Due to unique material compositions, neutron streaming paths, and neutronics requirements, the simplest of ITER neutronics calculations have to be thought out in great detail. Placing a new hole for a diagnostic creates a new streaming path; internal components must survive a high fast neutron flux, while external components face a completely different neutron energy spectrum of a high flux; material activation within human accessible areas must be maintained to an allowable level. As the complexity of the ITER device grows, so does the complexity of the neutronics challenges.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2022

Developments and Challenges in the Design of the ITER DRGA

We report the ITER Diagnostic Residual Gas Analyzer (DRGA) will measure the distribution of gas species, i.e., deuterium (D), tritium (T), and impurities, in the divertor exhaust stream and in the plasma periphery, with time resolution relevant to fusion plasma–wall particle dynamics. The uniqueness of the DRGA, over previous implementations of plasma dynamics residual gas analysis, is an integrated approach, combining mass and low-temperature plasma-activated optical spectroscopy, in a differentially pumped analysis station. A further unique feature of the ITER divertor-specific DRGA is an ~8-m separation of the analysis station from the sampled pumping duct, while still maintaining a ~1-s response time for hydrogen isotopic concentrations. ITER DRGA final design activities are strongly benefiting from testing of prototypical DRGA components and methods on present fusion devices, most currently on JET and W7-X. DRGA systems are implemented on both these devices and include sensors (and pumping methods) that are directly relevant to the ITER DRGA design. The recent JET-DTE2 campaign has provided the first experience on operating the combined ITER DRGA sensors with D-T plasmas. While enhancing system design for ITER, this experience on operating devices has also revealed additional engineering challenges, which further guide the continuing final design project. Meanwhile, the recent determination that the ITER DRGA, with slight optimization, will resolve the helium isotopes well enough to support an ITER pre-DT, He-3-based heating scheme, has greatly increased ITER Research Program interest in the DRGA and its implementation well ahead of the DT phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

xSDK-batched Subcontract - Ginkgo Batched Iterative Solver Development (Final Report)

Iterative solvers are fundamentally different from direct solvers in terms of execution as they generally do not execute a pre-defined sequence of operations or steps, but adapt the number of iterations to the specific problem and the preset solution quality. Generally, the adaptation of the iteration count to the problem is realized by monitoring the solver convergence and stopping the iteration process once the monitored metric, e.g., the residual norm, hits a pre-defined threshold. When addressing a set of problems with different properties, it is necessary to monitor the threshold for each problem individually and break up the SIMD execution style to avoid excess iterations for “easier” problems. Ginkgo integrates a simple but customizable stopping criterion for the residual norm and generally uses a pre-defined (relative or absolute) residual norm as the stopping criterion. In order to avoid the overhead of launching a kernel at every iteration, the iteration convergence and iteration control is part of the solver kernel. Each thread maintains its own copy of the iteration count.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2012

Influence of Primary Gage Sensitivities on the Convergence of Balance Load Iterations

The connection between the convergence of wind tunnel balance load iterations and the existence of the primary gage sensitivities of a balance is discussed. First, basic elements of two load iteration equations that the iterative method uses in combination with results of a calibration data analysis for the prediction of balance loads are reviewed. Then, the connection between the primary gage sensitivities, the load format, the gage output format, and the convergence characteristics of the load iteration equation choices is investigated. A new criterion is also introduced that may be used to objectively determine if the primary gage sensitivity of a balance gage exists. Then, it is shown that both load iteration equations will converge as long as a suitable regression model is used for the analysis of the balance calibration data, the combined influence of non linear terms of the regression model is very small, and the primary gage sensitivities of all balance gages exist. The last requirement is fulfilled, e.g., if force balance calibration data is analyzed in force balance format. Finally, it is demonstrated that only one of the two load iteration equation choices, i.e., the iteration equation used by the primary load iteration method, converges if one or more primary gage sensitivities are missing. This situation may occur, e.g., if force balance calibration data is analyzed in direct read format using the original gage outputs. Data from the calibration of a six component force balance is used to illustrate the connection between the convergence of the load iteration equation choices and the existence of the primary gage sensitivities.

Ulbrich, Norbert Manfred↗
DOE OSTI2015

Novel aspects of plasma control in ITER

ITER plasma control design solutions and performance requirements are strongly driven by its nuclear mission, aggressive commissioning constraints, and limited number of operational discharges. In addition, high plasma energy content, heat fluxes, neutron fluxes, and very long pulse operation place novel demands on control performance in many areas ranging from plasma boundary and divertor regulation to plasma kinetics and stability control. Both commissioning and experimental operations schedules provide limited time for tuning of control algorithms relative to operating devices. Although many aspects of the control solutions required by ITER have been well-demonstrated in present devices and even designed satisfactorily for ITER application, many elements unique to ITER including various crucial integration issues are presently under development. We describe selected novel aspects of plasma control in ITER, identifying unique parts of the control problem and highlighting some key areas of research remaining. Novel control areas described include control physics understanding (e.g., current profile regulation, tearing mode (TM) suppression), control mathematics (e.g., algorithmic and simulation approaches to high confidence robust performance), and integration solutions (e.g., methods for management of highly subscribed control resources). We identify unique aspects of the ITER TM suppression scheme, which will pulse gyrotrons to drive current within a magnetic island, and turn the drive off following suppression in order to minimize use of auxiliary power and maximize fusion gain. The potential role of active current profile control and approaches to design in ITER are discussed. Issues and approaches to fault handling algorithms are described, along with novel aspects of actuator sharing in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Radiation dose reduction in medical x-ray CT via Fourier-based iterative reconstruction

Purpose: A Fourier-based iterative reconstruction technique, termed Equally Sloped Tomography (EST), is developed in conjunction with advanced mathematical regularization to investigate radiation dose reduction in x-ray CT. The method is experimentally implemented on fan-beam CT and evaluated as a function of imaging dose on a series of image quality phantoms and anonymous pediatric patient data sets. Numerical simulation experiments are also performed to explore the extension of EST to helical cone-beam geometry. Methods: EST is a Fourier based iterative algorithm, which iterates back and forth between real and Fourier space utilizing the algebraically exact pseudopolar fast Fourier transform (PPFFT). In each iteration, physical constraints and mathematical regularization are applied in real space, while the measured data are enforced in Fourier space. The algorithm is automatically terminated when a proposed termination criterion is met. Experimentally, fan-beam projections were acquired by the Siemens z-flying focal spot technology, and subsequently interleaved and rebinned to a pseudopolar grid. Image quality phantoms were scanned at systematically varied mAs settings, reconstructed by EST and conventional reconstruction methods such as filtered back projection (FBP), and quantified using metrics including resolution, signal-to-noise ratios (SNRs), and contrast-to-noise ratios (CNRs). Pediatric data sets were reconstructed at their original acquisition settings and additionally simulated to lower dose settings for comparison and evaluation of the potential for radiation dose reduction. Numerical experiments were conducted to quantify EST and other iterative methods in terms of image quality and computation time. The extension of EST to helical cone-beam CT was implemented by using the advanced single-slice rebinning (ASSR) method. Results: Based on the phantom and pediatric patient fan-beam CT data, it is demonstrated that EST reconstructions with the lowest scanner flux setting of 39 mAs produce comparable image quality, resolution, and contrast relative to FBP with the 140 mAs flux setting. Compared to the algebraic reconstruction technique and the expectation maximization statistical reconstruction algorithm, a significant reduction in computation time is achieved with EST. Finally, numerical experiments on helical cone-beam CT data suggest that the combination of EST and ASSR produces reconstructions with higher image quality and lower noise than the Feldkamp Davis and Kress (FDK) method and the conventional ASSR approach. Conclusions: A Fourier-based iterative method has been applied to the reconstruction of fan-bean CT data with reduced x-ray fluence. This method incorporates advantageous features in both real and Fourier space iterative schemes: using a fast and algebraically exact method to calculate forward projection, enforcing the measured data in Fourier space, and applying physical constraints and flexible regularization in real space. Our results suggest that EST can be utilized for radiation dose reduction in x-ray CT via the readily implementable technique of lowering mAs settings. Numerical experiments further indicate that EST requires less computation time than several other iterative algorithms and can, in principle, be extended to helical cone-beam geometry in combination with the ASSR method.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Accelerated iterative beam angle selection in IMRT

Purpose: Iterative methods for beam angle selection (BAS) for intensity-modulated radiation therapy (IMRT) planning sequentially construct a beneficial ensemble of beam directions. In a naïve implementation, the nth beam is selected by adding beam orientations one-by-one from a discrete set of candidates to an existing ensemble of (n − 1) beams. The best beam orientation is identified in a time consuming process by solving the fluence map optimization (FMO) problem for every candidate beam and selecting the beam that yields the largest improvement to the objective function value. This paper evaluates two alternative methods to accelerate iterative BAS based on surrogates for the FMO objective function value. Methods: We suggest to select candidate beams not based on the FMO objective function value after convergence but (1) based on the objective function value after five FMO iterations of a gradient based algorithm and (2) based on a projected gradient of the FMO problem in the first iteration. The performance of the objective function surrogates is evaluated based on the resulting objective function values and dose statistics in a treatment planning study comprising three intracranial, three pancreas, and three prostate cases. Furthermore, iterative BAS is evaluated for an application in which a small number of noncoplanar beams complement a set of coplanar beam orientations. This scenario is of practical interest as noncoplanar setups may require additional attention of the treatment personnel for every couch rotation. Results: Iterative BAS relying on objective function surrogates yields similar results compared to naïve BAS with regard to the objective function values and dose statistics. At the same time, early stopping of the FMO and using the projected gradient during the first iteration enable reductions in computation time by approximately one to two orders of magnitude. With regard to the clinical delivery of noncoplanar IMRT treatments, we could show that optimized beam ensembles using only a few noncoplanar beam orientations often approach the plan quality of fully noncoplanar ensembles. Conclusions: We conclude that iterative BAS in combination with objective function surrogates can be a viable option to implement automated BAS at clinically acceptable computation times.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

TGLF Recalibration for ITER Standard Case Parameters FY2015: Theory and Simulation Performance Target Final Report

This work was motivated by the observation, as early as 2008, that GYRO simulations of some ITER operating scenarios exhibited nonlinear zonal-flow generation large enough to effectively quench turbulence inside r /a ~ 0.5. This observation of flow-dominated, low-transport states persisted even as more accurate and comprehensive predictions of ITER profiles were made using the state-of-the-art TGLF transport model. This core stabilization is in stark contrast to GYRO-TGLF comparisons for modern-day tokamaks, for which GYRO and TGLF are typically in very close agreement. So, we began to suspect that TGLF needed to be generalized to include the effect of zonal-flow stabilization in order to be more accurate for the conditions of reactor simulations. While the precise cause of the GYRO-TGLF discrepancy for ITER parameters was not known, it was speculated that closeness to threshold in the absence of driven rotation, as well as electromagnetic stabilization, created conditions more sensitive the self-generated zonal-flow stabilization than in modern tokamaks. Need for nonlinear zonal-flow stabilization: To explore the inclusion of a zonal-flow stabilization mechanism in TGLF, we started with a nominal ITER profile predicted by TGLF, and then performed linear and nonlinear GYRO simulations to characterize the behavior at and slightly above the nominal temperature gradients for finite levels of energy transport. Then, we ran TGLF on these cases to see where the discrepancies were largest. The predicted ITER profiles were indeed near to the TGLF threshold over most of the plasma core in the hybrid discharge studied (weak magnetic shear, q > 1). Scanning temperature gradients above the TGLF power balance values also showed that TGLF overpredicted the electron energy transport in the low-collisionality ITER plasma. At first (in Q3), a model of only the zonal-flow stabilization (Dimits shift) was attempted. Although we were able to construct an ad hoc model of the zonal flows that fit the GYRO simulations, the parameters of the model had to be tuned to each case. A physics basis for the zonal flow model was lacking. Electron energy transport at short wavelength: A secondary issue – the high-k electron energy flux – was initially assumed to be independent of the zonal flow effect. However, detailed studies of the fluctuation spectra from recent multiscale (electron and ion scale) GYRO simulations provided a critical new insight into the role of zonal flows. The multiscale simulations suggested that advection by the zonal flows strongly suppressed electron-scale turbulence. Radial shear of the zonal E×B fluctuation could not compete with the large electron-scale linear growth rate, but the k x -mixing rate of the E×B advection could. This insight led to a preliminary new model for the way zonal flows saturate both electron- and ion-scale turbulence. It was also discovered that the strength of the zonal E×B velocity could be computed from the linear growth rate spectrum. The new saturation model (SAT1), which replaces the original model (SAT0), was fit to the multiscale GYRO simulations as well as the ion-scale GYRO simulations used to calibrate the original SAT0 model. Thus, SAT1 captures the physics of both multiscale electron transport and zonal-flow stabilization. In future work, the SAT1 model will require significant further testing and (expensive) calibration with nonlinear multiscale gyrokinetic simulations over a wider variety of plasma conditions – certainly more than the small set of scans about a single C-Mod L-mode discharge. We believe the SAT1 model holds great promise as a physics-based model of the multiscale turbulent transport in fusion devices. Correction to ITER performance predictions: Finally, the impact of the SAT1model on the ITER hybrid case is mixed. Without the electron-scale contribution to the fluxes, the Dimits shift makes a significant improvement in the predicted fusion power as originally posited. Alas, including the high-k electron transport reduces the improvement, yielding a modest net increase in predicted fusion power compared to the TGLF prediction with the original SAT0 model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

First ITER CS module test results

The ITER Central Solenoid (CS) is under fabrication by the US ITER organization and its subcontractors. US ITER will supply seven modules to ITER IO, six of which will be assembled in a stack that forms the ITER Central Solenoid, with one as a spare. The first module fabrication has been completed by General Atomics (GA) at their facility and has begun testing including high voltage testing, Paschen testing in the vacuum and then testing at 4.5 K and up to 40 kA in order to demonstrate compliance of the coil to ITER requirements. In the paper we present the Test Plan and results of the CS Module performance tests, especially at 40 kA current. Herein, AC losses, joint resistances and hydraulic characteristics of the coil are all measured. Displacements of the coil height and hoop strain of the CS Module are also measured to verify structural and mechanical characteristics of the coil along with cooldown shrinkage of the coil. This information is used for verification of the stack behavior of CS in ITER operation. The test results and preliminary analyses results are presented, compared to expectations, and discussed. All measured parameters suggest that the CS module will perform well in ITER machine.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Test of the ITER-like resonant magnetic perturbation configurations for edge-localized mode crash suppression on KSTAR

KSTAR has demonstrated divertor heat flux broadening during edge-localized-mode (ELM)-crash-suppression using ITER-like 3-row resonant magnetic perturbation (RMP) for the first time. To address a couple of critical issues in ITER RMP, robust ELM-crash-suppression methodology has been explored at low q95 and established in KSTAR using low-nRMPs. Taking full advantage of the ITER-like 3-row in-vessel control coils (IVCC) in KSTAR, a set of intentionally misaligned RMP configurations (IMC) was tested to investigate whether or not IMC could be compatible with ELM-crash-suppression, while minimizing electromagnetic loads on RMP coils. As a result, the ITER-like 3-row IMCs were found not only to have been compatible with the ELM-crash-suppression, but also to have broadened the divertor heat fluxesin the vicinity of the outer strike point. In comparison, the 2-row RMPs have rarely affected the near scrape-off-layer (SOL) heat flux despite slightly broadened profile change in the far-SOL. Since the divertor heat flux broadening reflects the dispersal of the peaked near-SOL heat flux, the experimental outcome is quite favorable to the ITER choice of 3-rows, instead of 2-rows. Nonetheless, since the IMC-driven broadening observed in the attached plasmas in KSTAR might appear substantially different in the partially detached plasmas in ITER, additional investigation has been conducted to see if RMP-driven, ELM-crash-suppression could be compatible with detached plasmas. Although no detached plasmas have been identified with ELM-crash-suppression yet, significantly reduced divertor heat flux was confirmed in high density, ELM-crash-suppressed plasmas at q 95 =3.4using n=2 RMPs. These new findings elevate the confidence level about the ITER RMP system, while the remaining uncertainties need to be further clarified using the 3-row IVCCs in KSTAR. As long as mode-locking percussion is minimized along with a quick recovery of wall conditions, the access to the targeted q95(~ 3) for ITER is foreseen to be feasible in KSTAR

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Development of iterative techniques for the solution of unsteady compressible viscous flows

Efficient iterative solution methods are being developed for the numerical solution of two- and three-dimensional compressible Navier-Stokes equations. Iterative time marching methods have several advantages over classical multi-step explicit time marching schemes, and non-iterative implicit time marching schemes. Iterative schemes have better stability characteristics than non-iterative explicit and implicit schemes. Thus, the extra work required by iterative schemes can also be designed to perform efficiently on current and future generation scalable, missively parallel machines. An obvious candidate for iteratively solving the system of coupled nonlinear algebraic equations arising in CFD applications is the Newton method. Newton's method was implemented in existing finite difference and finite volume methods. Depending on the complexity of the problem, the number of Newton iterations needed per step to solve the discretized system of equations can, however, vary dramatically from a few to several hundred. Another popular approach based on the classical conjugate gradient method, known as the GMRES (Generalized Minimum Residual) algorithm is investigated. The GMRES algorithm was used in the past by a number of researchers for solving steady viscous and inviscid flow problems with considerable success. Here, the suitability of this algorithm is investigated for solving the system of nonlinear equations that arise in unsteady Navier-Stokes solvers at each time step. Unlike the Newton method which attempts to drive the error in the solution at each and every node down to zero, the GMRES algorithm only seeks to minimize the L2 norm of the error. In the GMRES algorithm the changes in the flow properties from one time step to the next are assumed to be the sum of a set of orthogonal vectors. By choosing the number of vectors to a reasonably small value N (between 5 and 20) the work required for advancing the solution from one time step to the next may be kept to (N+1) times that of a noniterative scheme. Many of the operations required by the GMRES algorithm such as matrix-vector multiplies, matrix additions and subtractions can all be vectorized and parallelized efficiently.

Sankar, Lakshmi N.↗