Search NASA⌕ Search

SEARCH · Search NASA

Results for “iterative analysis”

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

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

At least 1,081 records · Page 60

Analysis of a multi-machine database on divertor heat fluxes

A coordinated effort to measure divertor heat flux characteristics in fully attached, similarly shaped H-mode plasmas on C-Mod, DIII-D, and NSTX was carried out in 2010 in order to construct a predictive scaling relation applicable to next step devices including ITER, FNSF, and DEMO. Few published scaling laws are available and those that have been published were obtained under widely varying conditions and divertor geometries, leading to conflicting predictions for this critically important quantity. This study was designed to overcome these deficiencies. Analysis of the combined data set reveals that the primary dependence of the parallel heat flux width is robustly inverse with I{sub p}, which all three tokamaks independently demonstrate. An improved Thomson scattering system on DIII-D has yielded very accurate scrape off layer (SOL) profile measurements from which tests of parallel transport models have been made. It is found that a flux-limited model agrees best with the data at all collisionalities, while a Spitzer resistivity model agrees at higher collisionality where it is more valid. The SOL profile measurements and divertor heat flux scaling are consistent with a heuristic drift based model as well as a critical gradient model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A protection system for the JET ITER-like wall based on imaging diagnostics

The new JET ITER-like wall (made of beryllium and tungsten) is more fragile than the former carbon fiber composite wall and requires active protection to prevent excessive heat loads on the plasma facing components (PFC). Analog CCD cameras operating in the near infrared wavelength are used to measure surface temperature of the PFCs. Region of interest (ROI) analysis is performed in real time and the maximum temperature measured in each ROI is sent to the vessel thermal map. The protection of the ITER-like wall system started in October 2011 and has already successfully led to a safe landing of the plasma when hot spots were observed on the Be main chamber PFCs. Divertor protection is more of a challenge due to dust deposits that often generate false hot spots. In this contribution we describe the camera, data capture and real time processing systems. We discuss the calibration strategy for the temperature measurements with cross validation with thermal IR cameras and bi-color pyrometers. Most importantly, we demonstrate that a protection system based on CCD cameras can work and show examples of hot spot detections that stop the plasma pulse. The limits of such a design and the associated constraints on the operations are also presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

First Simultaneous Extraction of Spin-Dependent Parton Distributions and Fragmentation Functions from a Global QCD Analysis

In this paper, we perform the first global QCD analysis of polarized inclusive and semi-inclusive deep-inelastic scattering and single-inclusive $e^+e^-$ annihilation data, simultaneously fitting the parton distribution and fragmentation functions using the iterative Monte Carlo method. Without imposing SU(3) symmetry relations, we find the strange polarization to be very small, consistent with zero for both inclusive and semi-inclusive data, which provides a resolution to the strange quark polarization puzzle. Finally, the combined analysis also allows the direct extraction from data of the isovector and octet axial charges, and is consistent with a small SU(2) flavor asymmetry in the polarized sea.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Sensitivity Analysis of X-ray Spectra from Scanning Electron Microscopes

The primary goal of this project is to evaluate x-ray spectra generated within a scanning electron microscope (SEM) to determine elemental composition of small samples. This will be accomplished by performing Monte Carlo simulations of the electron and photon interactions in the sample and in the x-ray detector. The elemental inventories will be determined by an inverse process that progressively reduces the difference between the measured and simulated x-ray spectra by iteratively adjusting composition and geometric variables in the computational model. The intended benefit of this work will be to develop a method to perform quantitative analysis on substandard samples (heterogeneous phases, rough surfaces, small sizes, etc.) without involving standard elemental samples or empirical matrix corrections (i.e., true standardless quantitative analysis).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scalable, In-situ Data Clustering Data Analysis for Extreme Scale Scientific Computing (Final Report)

The objective of this project is to address challenges in the design and development of scalable in-situ data clustering and analytics algorithms and software. Our goal is to develop parallel software consisting of a set of spatio-temporal data clustering and anomaly detection functions, both of which are very important for large-scale analysis and have wide applicability for in-situ runs as well as post-processing analysis. Our design principles for in-situ analysis consider the following: (1) identify parts of the computation can be done close to the data within the nodes, while it is still in memory; (2) extract analysis components can (and should) be performed in remote staging and analysis nodes; (3) develop error-bound approximation methods for applications tolerable for small errors; (4) identify the type of derived distributions and statistics, for spatio-temporal data, that can be kept locally in order to both accelerate computations and meet energy constraints in subsequent iterations and phases; (5) use a self-describing data format so that data can be consistent and understood among local storage (memory and SSDs) and at staging and analysis nodes, thereby providing portability and flexibility; (6) develop service-oriented functions that can schedule in-situ and post-hoc analysis tasks based on the dynamic requirements of applications. Our development focus is to produce the parallel data analysis software/library that will be scalable, reusable, extensible, and generic for applications in different disciplines. The software will be able to run in-situ with the simulations as well as post-hoc analysis. This approach will satisfy many synergistic requirements for data intensive applications executed on data coming from instruments and experiments. In particular, the proposed multilevel approach is directly applicable to perform design tradeoffs for running part of the algorithms near the instruments and the rest on remote (analysis) systems.

97 MATHEMATICS AND COMPUTING↗

SU-F-207-02: Use of Postmortem Subjects for Subjective Image Quality Assessment in Abdominal CT Protocols with Iterative Reconstruction

Purpose: New radiation dose reduction technologies are emerging constantly in the medical imaging field. The latest of these technologies, iterative reconstruction (IR) in CT, presents the ability to reduce dose significantly and hence provides great opportunity for CT protocol optimization. However, without effective analysis of image quality, the reduction in radiation exposure becomes irrelevant. This work explores the use of postmortem subjects as an image quality assessment medium for protocol optimizations in abdominal CT. Methods: Three female postmortem subjects were scanned using the Abdomen-Pelvis (AP) protocol at reduced minimum tube current and target noise index (SD) settings of 12.5, 17.5, 20.0, and 25.0. Images were reconstructed using two strengths of iterative reconstruction. Radiologists and radiology residents from several subspecialties were asked to evaluate 8 AP image sets including the current facility default scan protocol and 7 scans with the parameters varied as listed above. Images were viewed in the soft tissue window and scored on a 3-point scale as acceptable, borderline acceptable, and unacceptable for diagnosis. The facility default AP scan was identified to the reviewer while the 7 remaining AP scans were randomized and de-identified of acquisition and reconstruction details. The observers were also asked to comment on the subjective image quality criteria they used for scoring images. This included visibility of specific anatomical structures and tissue textures. Results: Radiologists scored images as acceptable or borderline acceptable for target noise index settings of up to 20. Due to the postmortem subjects’ close representation of living human anatomy, readers were able to evaluate images as they would those of actual patients. Conclusion: Postmortem subjects have already been proven useful for direct CT organ dose measurements. This work illustrates the validity of their use for the crucial evaluation of image quality during CT protocol optimization, especially when investigating the effects of new technologies.

60 APPLIED LIFE SCIENCES↗

Iterative sampling of expensive simulations for faster deep surrogate training

Deep neural network (DNN) surrogates of expensive physics simulations are enabling a rapid change in the way that common experimental design and analysis tasks are approached. Surrogate models allow simulations to be performed in parallel and separately from downstream tasks, thereby enabling analyses that would be impossible with the simulation in-the-loop; surrogates based on DNNs can effectively emulate diverse non-scalar data of the types collected in fusion and laboratory-astrophysics experiments. The challenge is in training the surrogate model, for which large ensembles of physics simulations must be run, preferably without wasting computational effort on uninteresting simulations. Here, in this paper, we present an iterative sampling scheme that can preferentially propose simulations in interesting regions of parameter space without neglecting unexplored regions, allowing high-quality and wide-ranging surrogate models to be trained using 2–3 times fewer simulations compare to space-filling designs. Our approach uses an explicit importance function defined on the simulation output space, balanced against a measure of simulation density which serves as a proxy for surrogate accuracy. It is easy to implement and can be tuned to find interesting simulations early in the study, allowing surrogates to be trained quickly and refined as new simulations become available; this represents an important step towards the routine generation of deep surrogate models quickly enough to be truly relevant to experimental work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Provable Convergence of Plug-and-Play Priors With MMSE Denoisers

Plug-and-play priors (PnP) is a methodology for regularized image reconstruction that specifies the prior through an image denoiser. While PnP algorithms are well understood for denoisers performing maximum a posteriori probability (MAP) estimation, they have not been analyzed for the minimum mean squared error (MMSE) denoisers. Here we address this gap by establishing the first theoretical convergence result for the iterative shrinkage/thresholding algorithm (ISTA) variant of PnP for MMSE denoisers. We show that the iterates produced by PnP-ISTA with an MMSE denoiser converge to a stationary point of some global cost function. We validate our analysis on sparse signal recovery in compressive sensing by comparing two types of denoisers, namely the exact MMSE denoiser and the approximate MMSE denoiser obtained by training a deep neural net.

97 MATHEMATICS AND COMPUTING↗

Machine learning inversion of interatomic force constants from single-crystal inelastic neutron scattering

Atomic vibrations govern many macroscopic properties of materials, but experiments to comprehensively probe them remain challenging. Inelastic neutron scattering (INS) is a powerful technique to map phonon dispersions in crystals, especially when leveraging modern time-of-flight (ToF) spectrometers with large detectors. However, efficiently and robustly extracting interatomic force constants (FCs) parameterizing phonon dynamics from experimental spectra remains a bottleneck due to the complexity and high dimensionality of ToF INS datasets. Here, we present a machine learning approach for the direct inversion of FCs from single-crystal INS measurements. The framework leverages synthetic training data generated using universal machine-learned force fields and an efficient physics-based forward model. We benchmark two neural architectures–one emphasizing structured latent representation learning and the other direct, supervised spectral regression–across simulated datasets for two materials under idealized and noisy conditions. The latent-representation model is subsequently applied to experimental single-crystal INS data on germanium. The model is shown to reproduce FCs derived from both first-principles simulations and from iterative optimization, and furthermore achieves reliable inference even from sparse, single-orientation measurements representing short data acquisitions. Analysis of the learned latent space reveals semantically continuous and physically interpretable encodings that support strong cross-domain generalization. By bridging theoretical and experimental domains, we establish a path toward rapid inversion of experimental spectra and data-driven interpretation of temperature-dependent lattice dynamics.

42 ENGINEERING↗

Quantification Of Cesium In Negative Hydrogen Ion Sources By Laser Absorption Spectroscopy

The use of cesium in negative hydrogen ion sources and the resulting cesium dynamics caused by the evaporation and redistribution in the vacuum and plasma phase makes a reliable and on-line monitoring of the cesium amount in the source highly desirable. For that purpose, a robust and compact laser absorption setup suitable for the ion source environment has been developed utilizing the Cs D{sub 2} resonance line at 852.1 nm. First measurements are taken in a small laboratory plasma chamber with cesium evaporation. A detection limit of {approx_equal}5x10{sup 13} m{sup -3} at a typical path length of 15 cm has been obtained with a dynamic range of more than three orders of magnitude, limited by line saturation at high densities. For on-line monitoring an automatic data analysis is established achieving a temporal resolution of 100 ms. The setup has then been applied to the ITER prototype ion sources developed at IPP. It is been shown that the method is well suited for routine measurements revealing a new insight into the cesium dynamics during source operation and cesium conditioning.

43 PARTICLE ACCELERATORS↗

Extracted current saturation in negative ion sources

The extraction of negatively charged particles from a negative ion source is one of the crucial issues in the development of the neutral beam injector system for future experimental reactor ITER. Full 3D electrostatic particle-in-cell Monte Carlo collision code - ONIX [S. Mochalskyy et al., Nucl. Fusion 50, 105011 (2010)] - is used to simulate the hydrogen plasma behaviour and the extracted particle features in the vicinity of the plasma grid, both sides of the aperture. It is found that the contribution to the extracted negative ion current of ions born in the volume is small compared with that of ions created at the plasma grid walls. The parametric study with respect to the rate of negative ions released from the walls shows an optimum rate. Beyond this optimum, a double layer builds-up by the negative ion charge density close to the grid aperture surface reducing thus extraction probability, and therefore the extracted current. The effect of the extraction potential and magnetic field magnitudes on the extraction is also discussed. Results are in good agreement with available experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ion source development for a photoneutralization based NBI system for fusion reactors

The next step after ITER is to demonstrate the viability and generation of electricity by a future fusion reactor (DEMO). The specifications required to operate an NBI system on DEMO are very demanding. The system has to provide a very high level of power and energy, ~100MW of D° beam at 1MeV, including high wall-plug efficiency (η > 60%). For this purpose, a new injector concept, called Siphore, is under investigation between CEA and French universities. Siphore is based on the stripping of the accelerated negative ions by photo-detachment provided by several Fabry-Perot cavities (3.5MW of light power per cavity) implemented along the D{sup −} beam. The beamline is designed to be tall and narrow in order that the photon flux overlaps the entire negative ion beam. The paper will describe the present R and D at CEA which addresses the development of an ion source and pre-accelerator prototypes for Siphore, the main goal being to produce an intense negative ion beam sheet. The negative ion source Cybele is based on a magnetized plasma column where hot electrons are emitted from the source center. Parametric studies of the source are performed using Langmuir probes in order to characterize the plasma and to compare with numerical models being developed in French universities.

43 PARTICLE ACCELERATORS↗

Balancing Power Absorption Against Structural Loads With Viscous Drag and Power-Takeoff Efficiency Considerations

The focus of this paper is to balance power absorption against structural loading for a novel fixed-bottom oscillating surge wave energy converter in both regular and irregular wave environments. The power-to-load ratio will be evaluated using pseudospectral control (PSC) to determine the optimum power-takeoff (PTO) torque based on a multiterm objective function. This paper extends the pseudospectral optimal control problem to not just maximize the time-averaged absorbed power but also include measures for the surge-foundation force and PTO torque in the optimization. The objective function may now potentially include three competing terms that the optimizer must balance. Separate weighting factors are attached to the surge-foundation force and PTO control torque that can be used to tune the optimizer performance to emphasize either power absorption or load shedding. To correct the pitch equation of motion, derived from linear hydrodynamic theory, a quadratic-viscous-drag torque has been included in the system dynamics; however, to continue the use of quadratic programming solvers, an iteratively obtained linearized drag coefficient was utilized that provided good accuracy in the predicted pitch motion. Furthermore, the analysis considers the use of a nonideal PTO unit to more accurately evaluate controller performance. The PTO efficiency is not directly included in the objective function but rather the weighting factors are utilized to limit the PTO torque amplitudes, thereby reducing the losses resulting from the bidirectional energy flow through a nonideal PTO. Furthermore, results from PSC show that shedding a portion of the available wave energy can lead to greater reductions in structural loads, peak-to-average power ratio, and reactive power requirement.

16 TIDAL AND WAVE POWER↗

Dynamics of vortex dipoles in anisotropic Bose-Einstein condensates

We study the motion of a vortex dipole in a Bose-Einstein condensate confined to an anisotropic trap. We focus on a system of ODEs describing the vortices' motion, which is in turn a reduced model of the Gross-Pitaevskii equation describing the condensate's motion. Using a sequence of canonical changes of variables, we reduce the dimension and simplify the equations of motion. In this study, we uncover two interesting regimes. Near a family of periodic orbits known as guiding centers, we find that the dynamics is essentially that of a pendulum coupled to a linear oscillator, leading to stochastic reversals in the overall direction of rotation of the dipole. Near the separatrix orbit in the isotropic system, we find other families of periodic, quasi-periodic, and chaotic trajectories. In a neighborhood of the guiding center orbits, we derive an explicit iterated map that simplifies the problem further. Numerical calculations are used to illustrate the phenomena discovered through the analysis. Using the results from the reduced system, we are able to construct complex periodic orbits in the original, PDE, mean-field model for Bose-Einstein condensates, which corroborates the phenomenology observed in the reduced dynamical equations.

74 ATOMIC AND MOLECULAR PHYSICS↗

System design considerations for free-fall materials processing

The design constraints for orbiting materials processing systems are dominated by the limitations of the flight vehicle/crew and not by the processes themselves. Although weight, size and power consumption are all factors in the design of normal laboratory equipment, their importance is increased orders of magnitude when the equipment must be used in an orbital facility. As a result, equipment intended for space flight may have little resemblance to normal laboratory apparatus although the function to be performed may be identical. The same considerations influence the design of the experiment itself. The processing requirements must be carefully understood in terms of basic physical parameters rather than defined in terms of equipment operation. Preliminary experiments and analysis are much more vital to the design of a space experiment than they are on earth where iterative development is relatively easy. Examples of these various considerations are illustrated with examples from the M518 and MA-010 systems. While these are specific systems, the conclusions apply to the design of flight materials processing systems both present and future.

Seidensticker, R. G.↗

Inference of variations in the gravity field from satellite-to-satellite range rate

An analytic scheme for inferring variations of the gravity field from satellite-to-satellite range rate (low-low) is developed. As a test, it is applied to a pair of satellites in polar orbit, at altitude 160 km and spacing 100 km, with 72 data points per revolution. An assumed gravity field of tesseral spherical harmonics up to the eighth degree is completely recovered in three iterations over 64 revolutions. It is apparent that data points at regular intervals enable the use of data analysis techniques that avoid massive matrix inversions.

Kaula, W. M.↗

Infrared band extinctions and complex refractive indices of crystalline C2H2 and C4H2

Thermal IR absorption intensities are obtained for thin films of crystalline C2H2 and C4H2 at 70 K, and their n and k complex refractive indices are ascertained by separating true film absorption from interface reflection on the basis of an analysis of the transmission spectrum ratio for two sample thicknesses. This method significantly simplifies the n and k iteration process. The n and k values determined in the laboratory will in most cases reproduce a given sample thickness' observed transmission to within + or - 5 percent.

Khanna, R. K.↗