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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

SU-D-9A-02: Relative Effects of Threshold Choice and Spatial Resolution Modeling On SUV and Volume Quantification in F18-FDG PET Imaging of Anal Cancer Patients

Purpose: PET imaging with F18-FDG is utilized for treatment planning, treatment assessment, and prognosis. A region of interest (ROI) encompassing the tumor may be determined on the PET image, often by a threshold T on the PET standard uptake values (SUVs). Several studies have shown prognostic value for relevant ROI properties including maximum SUV value (SUVmax), metabolic tumor volume (MTV), and total glycolytic activity (TGA). The choice of threshold T may affect mean SUV value (SUVmean), MTV, and TGA. Recently spatial resolution modeling (SRM) has been introduced on many PET systems. SRM may also affect these ROI properties. The purpose of this work is to investigate the relative influence of SRM and threshold choice T on SUVmean, MTV, TGA, and SUVmax. Methods: For 9 anal cancer patients, 18F-FDG PET scans were performed prior to treatment. PET images were reconstructed by 2 iterations of Ordered Subsets Expectation Maximization (OSEM), with and without SRM. ROI contours were generated by 5 different SUV threshold values T: 2.5, 3.0, 30%, 40%, and 50% of SUVmax. Paired-samples t tests were used to compare SUVmean, MTV, and TGA (a) for SRM on versus off and (b) between each pair of threshold values T. SUVmax was also compared for SRM on versus off. Results: For almost all (57/60) comparisons of 2 different threshold values, SUVmean, MTV, and TGA showed statistically significant variation. For comparison of SRM on versus off, there were no statistically significant changes in SUVmax and TGA, but there were statistically significant changes in MTV for T=2.5 and T=3.0 and in SUVmean for all T. Conclusion: The near-universal statistical significance of threshold choice T suggests that, regarding harmonization across sites, threshold choice may be a greater concern than choice of SRM. However, broader study is warranted, e.g. other iterations of OSEM should be considered.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2018

Implementation and Use of a Convergence Test for Iterative Wind Tunnel Balance Load Predictions

A modified and extended version of a convergence test for wind tunnel strain-gage balance load iterations was implemented. The test uses an upper bound of the Lipschitz constant to assess convergence characteristics of balance load predictions if the Iterative Method is applied. Convergence is expected within the use envelope of the balance whenever this upper bound is less than the threshold of one. It is explained in great detail how the convergence test can be applied to the two load iteration equation types that are currently being used in the aerospace testing community. In addition, the application of the test to balances with bi-directional output characteristics is discussed. It is also shown how numerical differentiation can be used in order to obtain partial derivatives that are needed for the calculation of the Lipschitz constant. Finally, machine calibration data of NASA's MC60E six-component force balance is selected to demonstrate both implementation and use of the convergence test.

calibration analysis↗
DOE OSTI2026

Pedestal origin and extrapolation of high-density small edge-localised-modes peak parallel energy fluence in ITER and SPARC

Experimental analysis and simulations with the BOUT++ code show that small edge-localised modes (ELMs) in reactor-relevant high-density regimes originate in a region close to the separatrix and only marginally perturb the pedestal structure. The measured divertor peak parallel energy fluence (ε ∥,peak ) for a database of small ELM scenarios in DIII-D and ASDEX Upgrade can be reproduced, within 40 % accuracy on average, if an ad hoc modification of the Eich peak parallel ELM energy fluence model is applied to account for the small ELM pedestal birth location. This allows for first-order extrapolation of small-ELM divertor ε ∥,peak to ITER and SPARC, resulting in values that satisfy the nominal melting threshold of tungsten monoblocks of 12 MJ m −2 . The findings reported in this study, both via modelling and direct measurements, constitute a step forward in assessing small ELMs in high edge-collisionality scenarios as a viable plasma regime for the operation of next-generation fusion machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Divertor-safe nonlinear burn control based on a SOLPS parameterized core-edge model for ITER

Abstract For ITER operations, the range of desirable burning-plasma regimes with high fusion power output will be restricted by various operational constraints. These constraints include the saturation of ITER’s various heating and fueling actuators such as the neutral beam injectors, the ion and electron cyclotron heating systems, the gas puffing system, and the deuterium–tritium pellet injectors. In addition to these actuator constraints, the H-mode power threshold, divertor detachment, and the heat load on the divertor targets may apply limitations to ITER’s operational space. In this work, Plasma Operation Contour (POPCON) plots that map the aforementioned constraints to the temperature-density space are used to investigate which constraints are most limiting towards accessing regimes with high fusion power output. The presented POPCON plots are based on a control-oriented core-edge model that couples the nonlinear density and energy response models for the core-plasma region with SOLPS4.3 parameterizations for conditions in the edge-plasma regions (scrape-off-layer and divertor). Using this control-oriented core-edge model, a nonlinear burn controller, which aims to regulate the plasma temperature and density in the core-plasma region, is constructed in this work. This controller is augmented with an online optimization scheme that governs the control references such that the plasma can be guided towards regimes with high fusion powers while protecting the divertor targets from dangerously high heat loads. A closed-loop simulation study illustrates the capability of this burn control scheme.

Physics↗
DOE OSTI2019

The effect of plasma shape and neutral beam mix on the rotation threshold for RMP-ELM suppression

Varying the neutral beam injection (NBI) mix reveals a clear pedestal-top rotation threshold for edge localized mode (ELM) suppression by resonant magnetic perturbations. Guided by expectations for the RMP penetration mechanism, the rotation threshold is found to correspond to a critical radius for the ExB rotation zero-crossing. No such critical radius is observed for the electron perpendicular rotation zero-crossing. Varying the amount and ratio of power in different NBI source geometries (termed the NBI mix) also reveals that the rotation threshold can be crossed at widely varying total injected NBI torques. Computing the local torque density at the edge, the rotation threshold is found to be crossed when the local edge NBI torque is negative in nearly all discharges. Reducing the upper triangularity (delta_up ) from the ITER-similar value of 0.3 to 0.1 significantly reduces the pedestal height and width. This in turn: 1) Increases the rotation threshold and yields a more outward critical ExB rotation zero-crossing location. 2) Decreases the density threshold, consistent with a comparable collisionality range at lower pedestal temperatures. 3) Increases the input torque requirement, due to observed lower confinement and smaller intrinsic torque. These findings represent an important step along the road to predicting ELM suppression access conditions in future tokamaks such as ITER, where the toroidal rotation is expected to be low and consequently the rotation zero-crossing far from the pedestal-top.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

Higher Order, Hybrid BEM/FEM Methods Applied to Antenna Modeling

In this presentation, the authors address topics relevant to higher order modeling using hybrid BEM/FEM formulations. The first of these is the limitation on convergence rates imposed by geometric modeling errors in the analysis of scattering by a dielectric sphere. The second topic is the application of an Incomplete LU Threshold (ILUT) preconditioner to solve the linear system resulting from the BEM/FEM formulation. The final tOpic is the application of the higher order BEM/FEM formulation to antenna modeling problems. The authors have previously presented work on the benefits of higher order modeling. To achieve these benefits, special attention is required in the integration of singular and near-singular terms arising in the surface integral equation. Several methods for handling these terms have been presented. It is also well known that achieving ~he high rates of convergence afforded by higher order bases may als'o require the employment of higher order geometry models. A number of publications have described the use of quadratic elements to model curved surfaces. The authors have shown in an EFIE formulation, applied to scattering by a PEC .sphere, that quadratic order elements may be insufficient to prevent the domination of modeling errors. In fact, on a PEC sphere with radius r = 0.58 Lambda(sub 0), a quartic order geometry representation was required to obtain a convergence benefi.t from quadratic bases when compared to the convergence rate achieved with linear bases. Initial trials indicate that, for a dielectric sphere of the same radius, - requirements on the geometry model are not as severe as for the PEC sphere. The authors will present convergence results for higher order bases as a function of the geometry model order in the hybrid BEM/FEM formulation applied to dielectric spheres. It is well known that the system matrix resulting from the hybrid BEM/FEM formulation is ill -conditioned. For many real applications, a good preconditioner is required to obtain usable convergence from an iterative solver. The authors have examined the use of an Incomplete LU Threshold (ILUT) preconditioner . to solver linear systems stemming from higher order BEM/FEM formulations in 2D scattering problems. Although the resulting preconditioner provided aD excellent approximation to the system inverse, its size in terms of non-zero entries represented only a modest improvement when compared with the fill-in associated with a sparse direct solver. Furthermore, the fill-in of the preconditioner could not be substantially reduced without the occurrence of instabilities. In addition to the results for these 2D problems, the authors will present iterative solution data from the application of the ILUT preconditioner to 3D problems.

Fink, P. W.↗
DOE OSTI2016

Recursive Factorization of the Inverse Overlap Matrix in Linear Scaling Quantum Molecular Dynamics Simulations

We present a reduced complexity algorithm to compute the inverse overlap factors required to solve the generalized eigenvalue problem in a quantum-based molecular dynamics (MD) simulation. Our method is based on the recursive iterative re nement of an initial guess Z of the inverse overlap matrix S. The initial guess of Z is obtained beforehand either by using an approximate divide and conquer technique or dynamically, propagated within an extended Lagrangian dynamics from previous MD time steps. With this formulation, we achieve long-term stability and energy conservation even under incomplete approximate iterative re nement of Z. Linear scaling performance is obtained using numerically thresholded sparse matrix algebra based on the ELLPACK-R sparse matrix data format, which also enables e cient shared memory parallelization. As we show in this article using selfconsistent density functional based tight-binding MD, our approach is faster than conventional methods based on the direct diagonalization of the overlap matrix S for systems as small as a few hundred atoms, substantially accelerating quantum-based simulations even for molecular structures of intermediate size. For a 4,158 atom water-solvated polyalanine system we nd an average speedup factor of 122 for the computation of Z in each MD step.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simulation of pellet ELM triggering in low-collisionality, ITER-like discharges

3D nonlinear, as well as 2D linear M3D-C1 simulations are used to model ELM triggering by small pellets in DIII-D discharges in the ITER relevant, peeling-limited pedestal stability regime. A critical pellet size threshold is found in both experiment and modeling depending on pedestal conditions, pellet velocity and injection direction. Using radial injection at the outboard midplane, the threshold is determined by M3D-C1 for multiple time slices of a DIII-D low-collisionality discharge that has pellet ELM triggering. Experimental observations show that a larger pellet size than the standard 1.3 mm diameter is necessary for ELM triggering; 1.8 mm pellets triggered several ELMs in cases where a smaller pellet failed. The M3D-C1 simulations are in good agreement with these observations. While the 2D linear simulations give insight into the change of growth rates for various toroidal modes with pellet size, the 3D nonlinear simulations apply a pellet ablation model that mimics the actual injection with good match to the experiment. The 3D nonlinear simulation confirms the pellet ELM triggering for a pellet size larger than the threshold found by the linear simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Testing predictions of electron scale turbulent pedestal transport in two DIII-D ELMy H-modes

In this study, we present linear and nonlinear gyrokinetic analyses in the pedestal region of two DIII-D ELMy H-mode discharges using the CGYRO code. The otherwise matched discharges employ different divertor configurations to investigate the impact of varying recycling and particle source on pedestal profiles. Linear gyrokinetic simulations find electrostatic ion-scale instabilities (ion temperature gradient and trapped electron modes, ITG-TEM) are present just inside the top of the pedestal with growth rates that are enhanced significantly by parallel velocity shear. In the sharp gradient region, E x B shearing rates are comparable or larger than ion scale growth rates, suggesting the suppression of ITG-TEM modes in this region. Instead, the electron temperature profiles are found to be correlated with and just above the electron temperature gradient (ETG) instability thresholds. Using gradients varied within experimental uncertainties, nonlinear electron-scale gyrokinetic simulations predict electron heat fluxes from ETG turbulence, that when added to neoclassical ion thermal transport simulated by NEO, account for 30-60% of the total experimental heat flux. In addition, the neoclassical electron particle flux is found to contribute significantly to the experimental fluxes inferred from SOLPS-ITER analysis. Additional nonlinear gyrokinetic simulations are run varying input gradients to develop a threshold-based reduced model for ETG transport, finding a relatively simple dependence on η e = L ne /L Te . Predictive transport simulations are used to validate this pedestal-specific ETG model, in conjunction with a model for neoclassical particle transport. In both discharges, the predicted electron temperatures are always overpredicted, indicative of the insufficient stiffness in the ETG pedestal model to account for all of the experimental electron thermal transport. In the case of the closed divertor discharge with lower particle source, the predicted electron density is close to the experiment, consistent with the magnitude of neoclassical particle transport in that discharge. However, the density profiles are overpredicted in the open divertor discharge (larger particle source), due to insufficient model transport. In conclusion, the implications for other mechanisms accounting for the remainder of transport in the sharp gradient region in the two discharges are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

An ITPA joint experiment to study runaway electron generation and suppression

Recent results from an ITPA joint experiment to study the onset, growth, and decay of relativistic electrons (REs) indicate that loss mechanisms other than collisional damping may play a dominant role in the dynamics of the RE population, even during the quiescent I{sub p} flattop. Understanding the physics of RE growth and mitigation is motivated by the theoretical prediction that disruptions of full-current (15 MA) ITER discharges could generate up to 10 MA of REs with 10–20 MeV energies. The ITPA MHD group is conducting a joint experiment to measure the RE detection threshold conditions on a number of tokamaks under quasi-steady-state conditions in which V{sub loop}, n{sub e}, and REs can be well-diagnosed and compared to collisional theory. Data from DIII-D, C-Mod, FTU, KSTAR, and TEXTOR have been obtained so far, and the consensus to date is that the threshold E-field is significantly higher than predicted by relativistic collisional theory, or conversely, the density required to damp REs is significantly less than predicted, which could have significant implications for RE mitigation on ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

L-H Transition Dynamics in ITER-Similar Plasmas with Applied n=3 Magnetic Perturbations*

In ITER-similar shape (ISS), low rotation plasmas [1] in DIII-D [ =1.5-5x1019m-3, Bt=1.9-2T, Ip=1.5MA, q95~3.6, balanced neutral beam injection (NBI)], the L-H threshold power PLH with applied n=3 Resonant Magnetic Perturbations (RMP) [2,3] is found to increase strongly with decreasing collisionality. This is a concern for H-mode access in primarily ECH-heated ITER plasmas since RMP may be applied before the L-H transition in ITER to safely suppress the first ELM. The L-H power threshold increase with RMP is determined by initially accessing H-mode via applying balanced NBI and increasing steps in ECH power at high RMP field (applied via the DIII-D I-coil system consisting of six internal coils above and six coils below the outboard tokamak midplane [2]). NBI and ECH power are then removed to let the plasma return to L-mode. A second L-H transition is induced later in the same shot by applying NBI and stepping up ECH power again to obtain the threshold power with no applied RMP as a reference. Using this method to bracket PLH, a scan of PLH vs. RMP perturbation strength is then obtained via stepping down the I-coil current (RMP perturbation strength) at constant beam and ECH power, as shown in Fig. 1.

Schmitz, L.↗
NASA NTRS1993

An Unipolar Terminal-Attractor Based Neural Associative Memory with Adaptive Threshold and Perfect Convergence

For the first time, a unipolar terminal-attractor based neural associative memory (TABAM)system with adaptive threshold and perfect convergence is presented. By adaptively setting thethreshold values for the dynamic iteration for the unipolar binary neuron states with terminal-attractors and inner-product approach, we demonstrate via computer simulation the achievement ofperfect convergence and correct retrieval. The simulation is completed with a small number of storedstates (M) and a small number of neurons (N) but a large M/N ratio. An experiment with exclusive-or logic operation using LCTV SLMs is used to show feasibility of the optoelectronic implementationof the models.

Associative↗
DOE OSTI2013

Reduction of edge localized mode intensity on DIII-D by on-demand triggering with high frequency pellet injection and implications for ITER

The injection of small deuterium pellets at high repetition rates up to 12× the natural edge localized mode (ELM) frequency has been used to trigger high-frequency ELMs in otherwise low natural ELM frequency H-mode deuterium discharges in the DIII-D tokamak [J. L. Luxon and L. G. Davis, Fusion Technol. 8, 441 (1985)]. The resulting pellet-triggered ELMs result in up to 12× lower energy and particle fluxes to the divertor than the natural ELMs. The plasma global energy confinement and density are not strongly affected by the pellet perturbations. The plasma core impurity density is strongly reduced with the application of the pellets. These experiments were performed with pellets injected from the low field side pellet in plasmas designed to match the ITER baseline configuration in shape and normalized β operation with input heating power just above the H-mode power threshold. Nonlinear MHD simulations of the injected pellets show that destabilization of ballooning modes by a local pressure perturbation is responsible for the pellet ELM triggering. This strongly reduced ELM intensity shows promise for exploitation in ITER to control ELM size while maintaining high plasma purity and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Soft syndrome iterative decoding of quantum LDPC codes and hardware architectures

In practical quantum error correction implementations, the measurement of syndrome information is an unreliable step—typically modeled as a binary measurement outcome flipped with some probability. However, the measured syndrome is in fact a discretized value of the continuous voltage or current values obtained in the physical implementation of the syndrome extraction. In this paper, we use this “soft” or analog information to benefit iterative decoders for decoding quantum low-density parity-check (QLDPC) codes. Syndrome-based iterative belief propagation decoders are modified to utilize the soft syndrome to correct both data and syndrome errors simultaneously. We demonstrate the advantages of the proposed scheme not only in terms of comparison of thresholds and logical error rates for quasi-cyclic lifted-product QLDPC code families but also with faster convergence of iterative decoders. Additionally, we derive hardware (FPGA) architectures of these soft syndrome decoders and obtain similar performance in terms of error correction to the ideal models even with reduced precision in the soft information. The total latency of the hardware architectures is about 600 ns (for the QLDPC codes considered) in a 20 nm CMOS process FPGA device, and the area overhead is almost constant—less than 50% compared to min-sum decoders with noisy syndromes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1993

Dynamically reconfigurable optical morphological processor and its applications

An innovative optically implemented morphological processor is introduced. With the use of a large space-bandwidth-product Dammann grating and a high-speed shutter spatial light modulator, effective structuring element with large size and arbitrary shape can be constructed with dynamic reconfigurability. This reconfigurability is a major improvement over the conventional correlator-based morphological processor in which fixed holographic filters are used as structuring elements (Casasent and Botha, 1988). A novel two-dimensional thresholding photodetector array, capable of performing parallel thresholding and feedback, is utilized in this system and makes possible the implementation of many complex morphological operations requiring iterative feedbacks and full programmability. The optical architecture and the principle of operation are presented. Experimental demonstration of binary image morphological erosion, dilation, opening, and closing are also demonstrated. A technique for extending this technique to gray-scale image using thresholding decomposition technique is also discussed.

Chao, Tien-Hsin↗
DOE OSTI2014

An algorithm to estimate the object support in truncated images

Purpose: Truncation artifacts in CT occur if the object to be imaged extends past the scanner field of view (SFOV). These artifacts impede diagnosis and could possibly introduce errors in dose plans for radiation therapy. Several approaches exist for correcting truncation artifacts, but existing correction algorithms do not accurately recover the skin line (or support) of the patient, which is important in some dose planning methods. The purpose of this paper was to develop an iterative algorithm that recovers the support of the object. Methods: The authors assume that the truncated portion of the image is made up of soft tissue of uniform CT number and attempt to find a shape consistent with the measured data. Each known measurement in the sinogram is interpreted as an estimate of missing mass along a line. An initial estimate of the object support is generated by thresholding a reconstruction made using a previous truncation artifact correction algorithm (e.g., water cylinder extrapolation). This object support is iteratively deformed to reduce the inconsistency with the measured data. The missing data are estimated using this object support to complete the dataset. The method was tested on simulated and experimentally truncated CT data. Results: The proposed algorithm produces a better defined skin line than water cylinder extrapolation. On the experimental data, the RMS error of the skin line is reduced by about 60%. For moderately truncated images, some soft tissue contrast is retained near the SFOV. As the extent of truncation increases, the soft tissue contrast outside the SFOV becomes unusable although the skin line remains clearly defined, and in reformatted images it varies smoothly from slice to slice as expected. Conclusions: The support recovery algorithm provides a more accurate estimate of the patient outline than thresholded, basic water cylinder extrapolation, and may be preferred in some radiation therapy applications.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Effect of magnetic perturbations on turbulence-flow dynamics at the L-H transition on DIII-D

Detailed 2D turbulence measurements from the DIII-D tokamak provide an explanation for how resonant magnetic perturbations (RMPs) raise the L-H power threshold P LH [P. Gohil et al. , Nucl. Fusion 51 , 103020 (2011)] in ITER-relevant, low rotation, ITER-similar-shape plasmas with favorable ion ∇ B direction. RMPs simultaneously raise the turbulence decorrelation rate Δ ω D and reduce the flow shear rate ω shear in the stationary L-mode state preceding the L-H transition, thereby disrupting the turbulence shear suppression mechanism. RMPs also reduce the Reynolds stress drive for poloidal flow, contributing to the reduction of ω shear . On the ~100 μ s timescale of the L-H transition, RMPs reduce Reynolds-stress-driven energy transfer from turbulence to flows by an order of magnitude, challenging the energy depletion theory for the L-H trigger mechanism. In contrast, non-resonant magnetic perturbations, which do not significantly affect P LH , do not affect Δ ω D and only slightly reduce ω shear and Reynolds-stress-driven energy transfer.<

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗