Search NASA⌕ Search

SEARCH · Search NASA

Results for “Iterative”

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 OSTI2023

Feasibility of raised inner strike point equilibria scenario in ITER for detritiation from beryllium co-deposits

Abstract In ITER, tritium retention primarily occurs through co-deposition with beryllium. To avoid exceeding the strict tritium inventory limit, efficient tritium recovery techniques are essential. Baking is the ITER baseline for tritium recovery, but its effectiveness in removing tritium from thick beryllium layers is limited. A raised strike point scenario is considered an alternative method for removing tritium from the ITER inner vertical divertor target by heating components via plasma flux. This paper presents SOLPS-ITER code simulations conducted under various conditions, assessing the divertor performance and tritium outgassing of the raised strike point scenario. As the strike point is raised, recycled neutrals are not efficiently baffled by the dome and scrape-off layer, significantly changing the neutral trajectory and ionization source distribution. This improves detachment accessibility but worsens core-edge compatibility compared to the baseline scenario. However, in the partially detached condition, the impact of raising the strike point, perpendicular transport, and q 95 on target heat flux is not significant, as it primarily scales with the input power. Target heat flux is translated to target surface temperature using a simplified heat transfer model that considers the 3D target monoblock geometry and active cooling condition, excluding Be layer thermal properties. For partially detached divertor conditions, the bulk tungsten monoblock surface temperature remains below the baking temperature, which is insufficient for efficient tritium outgassing under the actively cooled ITER divertor condition. However, considering the potential thermal contact resistance between the beryllium and tungsten layers, which may significantly impact temperature distribution, the temperature of the beryllium layer can be raised to a level sufficient for efficient tritium outgassing. Therefore, the raised strike point scenario can be considered as an alternative in-vessel tritium removal technique.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Iterative reconstruction excursions for Baryon Acoustic Oscillations and beyond

ABSTRACT The density field reconstruction technique has been widely used for recovering the baryon acoustic oscillation (BAO) feature in galaxy surveys that has been degraded due to non-linearities. Recent studies advocated adopting iterative steps to improve the recovery much beyond that of the standard technique. In this paper, we investigate the performance of a few selected iterative reconstruction techniques focusing on the BAO and the broad-band shape of the two-point clustering. We include redshift-space distortions, halo bias, and shot noise and inspect the components of the reconstructed field in Fourier space and in configuration space using both density field-based reconstruction and displacement field-based reconstruction. We find that the displacement field reconstruction becomes quickly challenging in the presence of non-negligible shot noise and therefore present surrogate methods that can be practically applied to a much more sparse field such as galaxies. For a galaxy field, implementing a debiasing step to remove the Lagrangian bias appears crucial for the displacement field reconstruction. We show that the iterative reconstruction does not substantially improve the BAO feature beyond an aggressively optimized standard reconstruction with a small smoothing kernel. However, we find taking iterative steps allows us to use a small smoothing kernel more ‘stably’, i.e. without causing a substantial deviation from the linear power spectrum on large scales. In one specific example we studied, we find that a deviation of 13 per cent in $P(k\sim 0.1\, h{\rm \,\,Mpc^{-1}})$ with an aggressive standard reconstruction can reduce to 3–4 per cent with iterative steps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Modeling iterative reconstruction and displacement field in the large scale structure

The next generation of galaxy surveys like the Dark Energy Spectroscopic Instrument and Euclid will provide datasets orders of magnitude larger than anything available to date. Our ability to model nonlinear effects in late time matter perturbations will be a key to unlock the full potential of these datasets, and the area of initial condition reconstruction is attracting growing attention. Iterative reconstruction developed by Marcel Schmittfull et al. is a technique designed to reconstruct the displacement field from the observed galaxy distribution. The nonlinear displacement field and initial linear density field are highly correlated. Therefore, reconstructing the nonlinear displacement field enables us to extract the primordial cosmological information better than from the late time density field at the level of the two-point statistics. This paper will test to what extent the iterative reconstruction can recover the true displacement field and construct a perturbation theory model for the postreconstructed field. We model the iterative reconstruction process with Lagrangian perturbation theory (LPT) up to third order for dark matter in real space and compare it with N-body simulations. Here we find that the simulated iterative reconstruction does not converge to the nonlinear displacement field, and the discrepancy mainly appears in the shift term, i.e., the term correlated directly with the linear density field. On the contrary, our 3LPT model predicts that the iterative reconstruction should converge to the nonlinear displacement field. We discuss the sources of discrepancy, including numerical noise/artifacts on small scales, and present an ad hoc phenomenological model that improves the agreement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Fast Iterative Multi-site Hosting Capacity Analysis for Distribution Systems With Search Space Pruning

Interconnection studies for distributed energy resources (DERs) is a time-intensive process, primarily due to the necessity of solving large number of power flow scenarios. Hosting capacity analysis (HCA) is a time-consuming aspect of interconnection studies that is divided into single-site HCA (SHCA) and multi-site HCA (MHCA). From a computational and understandable standpoint, the industry seeks iteration-based solutions for SHCA, although it doesn't maximize the total DER hosting capacity (DERHC) of the grid, as MHCA does. While non-iterative solutions are available for MHCA, they involve a trade-off between the modeling accuracy of the distribution system, solution quality, and ease of understanding. In this work, we present a fast iterative solution for MHCA, reducing computational complexity by eliminating the need to solve power flows for a large amount of search space, thus making iterative solutions feasible. This iterative approach guarantees both a global optimal solution with sufficient time and a fast, close-to-optimal solution through efficient search space pruning. It also easily integrates with existing utility HCA tools. The results are demonstrated on select locations in the IEEE-123 bus system for community-scale interconnection studies. We highlight the benefits of skipping the need to solve millions of power flows, all while maximizing the grid's total DERHC.

Guddanti, Kishan Prudhvi↗
DOE OSTI2023

An Iterative Bidirectional Gradient Boosting Algorithm for CVR Baseline Estimation

This paper presents a novel iterative, bidirectional, gradient boosting (bidirectional-GB) algorithm for estimating the baseline of the Conservation Voltage Reduction (CVR) program. We define the CVR baseline as the load profile during the CVR period if the substation voltage is not lowered. The proposed algorithm consists of two key steps: selection of similar days and iterative bidirectional-GB training. In the first step, preand post-event temperature profiles of the targeted CVR day are used to select similar days from historical non-CVR days. In the second step, the pre-event and post-event similar days are used to train two GBMs iteratively: a forward-GBM and a backwardGBM. After each iteration, the two generated CVR baselines are reconciled and only the first and the last points on the reconciled baseline are kept. The iteration repeats until all CVR baseline points are generated. We tested two gradient boosting methods (i.e., GBM and LighGBM) with two data resolutions (i.e., 15and 30-minute). The results demonstrate that both the accuracy and performance of the algorithm are satisfactory.

Lee, Han Pyo↗
DOE OSTI2025

Accelerating iterative ptychography with an integrated neural network

Electron ptychography is a powerful and versatile tool for high-resolution and dose-efficient imaging. Iterative reconstruction algorithms are powerful but also computationally expensive due to their relative complexity and the many hyperparameters that must be optimised. Gradient descent-based iterative ptychography is a popular method, but it may converge slowly when reconstructing low spatial frequencies. Here, in this work, we present a method for accelerating a gradient descent-based iterative reconstruction algorithm by training a neural network (NN) that is applied in the reconstruction loop. The NN works in Fourier space and selectively boosts low spatial frequencies, thus enabling faster convergence in a manner similar to accelerated gradient descent algorithms. We discuss the difficulties that arise when incorporating a NN into an iterative reconstruction algorithm and show how they can be overcome with iterative training. We apply our method to simulated and experimental data of gold nanoparticles on amorphous carbon and show that we can significantly speed up ptychographic reconstruction of the nanoparticles.

4DSTEM↗
NASA NTRS1977

An implicit-iterative solution of the heat conduction equation with a radiation boundary condition

For the problem of predicting one-dimensional heat transfer between conducting and radiating mediums by an implicit finite difference method, four different formulations were used to approximate the surface radiation boundary condition while retaining an implicit formulation for the interior temperature nodes. These formulations are an explicit boundary condition, a linearized boundary condition, an iterative boundary condition, and a semi-iterative boundary method. The results of these methods in predicting surface temperature on the space shuttle orbiter thermal protection system model under a variety of heating rates were compared. The iterative technique caused the surface temperature to be bounded at each step. While the linearized and explicit methods were generally more efficient, the iterative and semi-iterative techniques provided a realistic surface temperature response without requiring step size control techniques.

Williams, S. D.↗
NASA NTRS1984

Second-order p-iterative solution of the Lambert/Gauss problem

An algorithm is presented for efficient p-iterative solution of the Lambert/Gauss orbit-determination problem using second-order Newton iteration. The algorithm is based on a universal transformation of Kepler's time-of-flight equation and approximate inverse solutions of this equation for short-way and long-way flight paths. The approximate solutions provide both good starting values for iteration and simplified computation of the second-order term in the iteration formula. Numerical results are presented which indicate that in many cases of practical significance (except those having collinear position vectors) the algorithm produces at least eight significant digits of accuracy with just two or three steps of iteration.

Boltz, F. W.↗
NASA NTRS1985

Iterative methods for mixed finite element equations

Iterative strategies for the solution of indefinite system of equations arising from the mixed finite element method are investigated in this paper with application to linear and nonlinear problems in solid and structural mechanics. The augmented Hu-Washizu form is derived, which is then utilized to construct a family of iterative algorithms using the displacement method as the preconditioner. Two types of iterative algorithms are implemented. Those are: constant metric iterations which does not involve the update of preconditioner; variable metric iterations, in which the inverse of the preconditioning matrix is updated. A series of numerical experiments is conducted to evaluate the numerical performance with application to linear and nonlinear model problems.

Nakazawa, S.↗
NASA NTRS1992

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

The development of efficient iterative solution methods for the numerical solution of two- and three-dimensional compressible Navier-Stokes equations is discussed. 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. In this work, another approach based on the classical conjugate gradient method, known as the Generalized Minimum Residual (GMRES) algorithm is investigated. The GMRES algorithm has been used in the past by a number of researchers for solving steady viscous and inviscid flow problems. Here, we investigate the suitability of this algorithm for solving the system of non-linear equations that arise in unsteady Navier-Stokes solvers at each time step.

Sankar, Lakshmi N.↗
NASA NTRS1992

A subspace iteration for eigenvector derivatives

An iterative procedure for calculating eigenvector derivatives has been developed based on the subspace iteration concept. The basic formulation is derived directly from the first variation of the resulting equations of subspace iteration for solving eigenvalue problems. Since the basic formulation conserves the convergence properties of the original subspace iteration for eigen-problems, an overrelaxation scheme similar to Bathe's approach is employed to accelerate the subspace iteration process for calculating eigenvector derivatives.

Ting, T.↗
DOE OSTI2024

THERMAL STRUCTURAL AND SHOCK EVENT EVALUATIONS OF THE FUELING PELLET INJECTION SYSTEM FOR ITER

The US is among seven partner nations in a collaborative effort to design, build, and demonstrate fusion’s ability to be a large-scale carbon free energy source. Each country has its own Domestic Agencies (DA) that contribute directly to the ITER project. US ITER, which is a DOE Office of Science project managed by Oak Ridge National Laboratory, is developing world class engineering solutions to the design, construction, and assembly of the burning plasma experiment that can demonstrate the scientific and technological feasibility of fusion. US ITER’s scope includes completing the preliminary and final design, qualifying materials and processes for manufacture and testing, executing manufacturing, and delivering the Fueling Pellet Injection System (FPIS) to the ITER site for assembly. The US-ITER FPIS system is designed to inject cryogenically frozen pellets of deuterium-tritium (D-T) into the plasma. The FPIS has two main functions: 1. Provide a steady supply of deuterium and tritium fuel, 2. Mitigating the impact of edge localized modes on the plasma facing components. The FPIS will have the capability to reside in three port cells within a tritium second barrier containment cask. Each pellet cask contains three flight tubes linking with penetrations on the torus cryopump housing and vacuum vessel (VV): two for the magnetic high field side (HFS) pellet injection and one for the magnetic low field side (LFS) pellet injection. This paper presents results of the thermal-structural analyses of the FPIS flight tube structural components when subjected to various loads such as electromagnetic (EM) loads, nuclear heating, seismic, and operational, and dynamic shock events. The resulting temperatures and stresses under combined conditions have been found to satisfy the design criteria to ensure safe and reliable operation of the FPIS flight tubes within the vacuum vessel.

Martinez, Oscar↗
DOE OSTI2023

Convergence analysis of single rate and multirate fixed stress split iterative coupling schemes in heterogeneous poroelastic media

Recently, the accurate modeling of flow–structure interactions has gained more attention and importance for both petroleum and environmental engineering applications. Of particular interest is the coupling between subsurface flow and reservoir geomechanics. Different single rate and multirate iterative and explicit coupling schemes have been proposed and analyzed in the past. In addition, Banach fixed point contraction results were obtained for iterative coupling schemes, and conditionally stable results were obtained for explicit coupling schemes. In this work, we will consider the mathematical analysis of the single rate and multirate fixed stress split iterative coupling schemes for spatially heterogeneous poroelastic media. We will re–establish the contractivity for both schemes in the localized case, and we will show that heterogeneities come at the expense of imposing more restricted conditions on the number of fine flow time steps that can be taken within one coarse mechanics time step in the multirate case. Our mathematical analysis is supplemented by numerical simulations validating our derived upper bounds. Finally, to the best of our knowledge, this is the first rigorous mathematical analysis of the multirate fixed–stress split iterative coupling scheme in heterogeneous poroelastic media.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Status of the development and testing of in-vessel and ECH-protection components for the ITER low-field side reflectometer

The ITER Low-Field Side Reflectometer (LFSR) is a critical diagnostic system designed to measure edge electron density profiles, fluctuations, and plasma rotation in ITER. This paper presents the latest developments in the design, testing, and validation of key in-vessel and Electron Cyclotron Heating (ECH) protection components. The LFSR antenna array has been optimized to provide robust coverage over expected plasma vertical displacements and a Doppler measurement for plasma rotation. Further, a comprehensive assessment of a 170-GHz diffraction grating and a novel stray-ECH power monitor demonstrates effectiveness in mitigating the impact of stray ECH power on sensitive microwave electronics. Additionally, an advancement in anti-reflective technology for millimeter waves significantly improves transmission of LFSR's vacuum window while meeting ITER's stringent safety and operational requirements. These results support the continued integration of LFSR into ITER, ensuring its diagnostic capabilities remain resilient under reactor-relevant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Overview of JET results for optimising ITER operation

The JET 2019–2020 scientific and technological programme exploited the results of years of concerted scientific and engineering work, including the ITER-like wall (ILW: Be wall and W divertor) installed in 2010, improved diagnostic capabilities now fully available, a major neutral beam injection upgrade providing record power in 2019–2020, and tested the technical and procedural preparation for safe operation with tritium. Research along three complementary axes yielded a wealth of new results. Firstly, the JET plasma programme delivered scenarios suitable for high fusion power and alpha particle (α) physics in the coming D–T campaign (DTE2), with record sustained neutron rates, as well as plasmas for clarifying the impact of isotope mass on plasma core, edge and plasma-wall interactions, and for ITER pre-fusion power operation. The efficacy of the newly installed shattered pellet injector for mitigating disruption forces and runaway electrons was demonstrated. Secondly, research on the consequences of long-term exposure to JET-ILW plasma was completed, with emphasis on wall damage and fuel retention, and with analyses of wall materials and dust particles that will help validate assumptions and codes for design and operation of ITER and DEMO. Thirdly, the nuclear technology programme aiming to deliver maximum technological return from operations in D, T and D–T benefited from the highest D–D neutron yield in years, securing results for validating radiation transport and activation codes, and nuclear data for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Loss of energetic particles due to resistive wall mode instability in ITER

Effects of an unstable n = 1 (n is the toroidal mode number) resistive wall mode (RWM) on the energetic particle (EP) confinement and loss are numerically investigated, for an ITER steady state scenario with 10 MA plasma current and 5.3 T toroidal field. The eigenfunction of the RWM is computed, with the associated three-dimensional magnetic field perturbation superposed with the 2D equilibrium field for tracing the EP drift orbits. Considered are mono-energetic EPs at 0.5 MeV and 1 MeV for deuterium ions, and 3.5 MeV for fusion-born alphas, with a range of distribution in the particle pitch angle. Modeling finds that less than 20% of EPs can be lost to the limiting surface in ITER assuming a source distribution uniform in minor radius, due to an unstable RWM that produces 100 Gauss poloidal field perturbation at the outboard mid-plane just inside the (effective) resistive wall surface. On top of the initial prompt drift orbit loss for counter-current EPs, the RWM induced particle loss occurs on a one second time scale, which is comparable to the RWM growth time in ITER. The 'wetted' area, due to the lost EPs striking the limiting surface, is generally found to be large due to the RWM. This is a favorable prediction for ITER. Here, the loss distribution in the poloidal angle is more uniform for co-current EPs. Counter-current EPs experience outward orbit drift when launched from the low-field side, and tend to more often hit the bottom region of the limiting surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Toroidal modeling of runaway electron loss due to 3D fields in ITER

Mitigation of runaway electrons (REs) by three-dimensional (3D) magnetic field perturbations is numerically investigated for the ITER 15 MA baseline D–T scenario, utilizing the MARS-F code (Liu et al Phys. Plasmas 7 3681) with a drift orbit test particle tracing module. Considered are two types of 3D fields: the n = 3 (n is the toroidal mode number) resonant magnetic perturbation (RMP) utilized for the purpose of controlling the edge localized modes in ITER, and perturbations generated by the n = 1 magneto-hydrodynamic (MHD) instabilities in a post-disruption plasma. The RMP field, applied to a pre-disruption plasma, is found to be moderately effective in mitigating the RE seeds in ITER when vacuum field model is assumed. Up to ~40% loss fraction is possible at 90 kA-turn coil current. The mitigation efficiency is however substantially reduced, down to less than 5%, when the plasma response is taken into account. This is due to strong screening of the resonant magnetic field components by the plasma response resulting in much less field line stochasticity. On the other hand, the MARS-F modeling, based on the DINA-simulated post-disruption equilibria, shows that the n = 1 resistive kink instabilities develop in these plasmas, as the edge safety factor qa evolves and drops below integer numbers. RE mitigation by these MHD instabilities is sensitive to the eigenmode structure. The best mitigation is achieved as q a drops below 3, when a global kink instability occurs that encompasses both internal and external components. This global instability is found to be capable of mitigating over 80% MeV-level passing RE orbits at a field perturbation |δB|/B 0 that is comparable to that observed in DIII-D experiments, and full mitigation if the perturbation amplitude is doubled. The 'wetted' area on the ITER limiting surface, due to MHD instability induced RE loss, generally increases with the perturbation amplitude (together with increasing loss fraction). At the highest perturbation level assumed in this study, the wetted area reaches ~60% of the total limiting surface area. Lastly, the lost RE orbits mainly strike the outer divertor region of the limiting surface, with some fraction also hitting a wide area along the inboard side of the surface.

3D fields↗
DOE OSTI2023

Evaluation of ITER divertor shunts as a synthetic diagnostic for detachment control

Reliable diagnostics that measure the detached state of the ITER divertor plasma will be necessary to control heat flux to the divertor targets during steady state, burning plasma operation. This paper conducts an initial exploration into the feasibility of the divertor shunt diagnostic as a lightweight, robust, and real-time detachment sensor. This diagnostic is a set of shunt lead pairs that measure the voltage drop along the divertor cassette body, from which the plasma scrape-off layer (SOL) current is calculated. Using SOLPS-ITER simulations for control-relevant ITER plasma scenarios, the thermoelectric current magnitude along the SOL is shown to decrease significantly with the onset of partial detachment at the outer divertor target. Electromagnetic modelling of a simplified divertor cassette is used to develop a control-oriented inductance-resistance circuit model, from which SOL currents can be calculated from shunt pair voltage measurements. The sensitivity and frequency-response of the resulting system indicates that the diagnostic will accurately measure SOL thermoelectric currents during ITER operation. These currents will be a good measure of the detached state of the divertor plasma, making the divertor shunt diagnostic a potentially extremely valuable and physically robust sensor for real-time detachment control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗