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 883 records · Page 49

Iterative solution of large, sparse linear systems on a static data flow architecture - Performance studies

The applicability of static data flow architectures to the iterative solution of sparse linear systems of equations is investigated. An analytic performance model of a static data flow computation is developed. This model includes both spatial parallelism, concurrent execution in multiple PE's, and pipelining, the streaming of data from array memories through the PE's. The performance model is used to analyze a row partitioned iterative algorithm for solving sparse linear systems of algebraic equations. Based on this analysis, design parameters for the static data flow architecture as a function of matrix sparsity and dimension are proposed.

Reed, D. A.↗

Domain decomposition algorithms and computation fluid dynamics

In the past several years, domain decomposition was a very popular topic, partly motivated by the potential of parallelization. While a large body of theory and algorithms were developed for model elliptic problems, they are only recently starting to be tested on realistic applications. The application of some of these methods to two model problems in computational fluid dynamics are investigated. Some examples are two dimensional convection-diffusion problems and the incompressible driven cavity flow problem. The construction and analysis of efficient preconditioners for the interface operator to be used in the iterative solution of the interface solution is described. For the convection-diffusion problems, the effect of the convection term and its discretization on the performance of some of the preconditioners is discussed. For the driven cavity problem, the effectiveness of a class of boundary probe preconditioners is discussed.

Chan, Tony F.↗

Orbiter Trajectory Analysis for a Two-Stage Reusable Launch Vehicle

Trajectory analysis performed on NASA's reference two-stage-to-orbit launch vehicle upper stage will be presented. The work was completed in support of the Hypersonics Multidisciplinary Analysis and Optimization effort for the NASA-Air Force Joint System Study. Three degree-of-freedom (3-DOF) untrimmed trajectory analysis was performed for the orbiter ascent, closure and re-entry. An iterative closure process resulted in a 333,000 lb initial mass for the orbiter. The re-entry trajectory satisfied heating constraints for all payload out cases and met the constraints with reduced margins for payload in cases. Abort trajectories for engine out at staging, engine out during ascent, and failure to circularize in orbit, gave insight to the robustness of the orbiter. A trimmed ascent trajectory defined an engine gimbal location and the body flap angle best suited for maximizing injected mass. A trimmed re-entry trajectory revealed a need to update the trim routine to accommodate full flap aerodynamic data.

Cowling, Adam L.↗

TPSAS-NF1676L-33433-DND

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction, and are characterized by a complex valued frequency-dependent impedance. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic scattering of geometric bodies treated with acoustic liners. This work considers using either an impedance boundary condition or admittance boundary condition to allow acoustic scattering problems to be modeled with geometries consisting of both un-lined and lined surfaces; admittance is defined to be the inverse of impedance. Three different acoustic liner models will be discussed: the Extended Helmholtz Resonator Model, the Three-Parameter Impedance Model, and a Broadband Model. In the Extended Helmholtz Resonator Model and Three-Parameter Impedance Model, the liner impedance is specified at a single frequency; the Broadband Model allows for the investigation of multiple frequencies simultaneously. The impedance and admittance boundary conditions for each model will be derived and coupled with a time-domain boundary integral equation. The scattering solution will be obtained iteratively using spatial and temporal basis functions. Stability will be demonstrated through eigenvalue analysis. Moreover, fast algorithms and high performance computing are considered to both assess and reduce the computational cost of the simulation.

Michelle E Rodio↗

Feasibility Analysis of Unmanned Aerial Vehicle Based on Tensegrity Structure

For future unmanned aerial vehicles to be ubiquitous, they need to be safe to other aircraft and people on the ground even in the event of unintended collisions. Unfortunately, this is difficult for current rigid designs as a crash of all but the smallest aircraft will have a significant impact on its target. As an alternative, this pa-per explores using a soft tensegrity structure (interconnected cables and rods) as a lightweight, robust chassis for a UAS and analyze its ability to achieve stable, controllable flight and mitigate damage during a collision. Tensegrity structures are three dimensional geometric constructions which are able to undergo severe deformations without permanent damage, and can immediately return to their original shape. Experiments and analysis are performed using a simple physics simulation environment, allowing for rapid design iteration and thorough exploration of system behavior. We test multiple configurations of the vehicle to understand impacts of changing various physical characteristics, and take limitations such as modern propeller thrust coefficients and motor speeds into account to ensure our results are inline with those of a real world vehicle. We ultimately find that such a vehicle has considerable real world potential, and that a physical prototype could theoretically be constructed and flown with current technology.

Bjorn Johnson↗

Front Range Wildland Fires: Evaluating the Efficacy of Remote Sensing Imagery in Monitoring Forest Fuels Treatment Methods

Over the last several decades, wildfire frequency and severity in forested areas along Colorado’s Front Range have increased due to a buildup of fuels. This has led to an increase in forest treatments, as well as an increased need to evaluate the success of these treatments. Remote sensing products offer an efficient and cost-effective way to monitor forest treatments; however, not all remote sensing products and analysis techniques have been explored by Coloradan land managers. Specifically, project partners at the Colorado State Forest Service (CSFS) and the Colorado Forest Restoration Institute (CFRI) were interested in using an effective and streamlined method of mapping canopy cover to better monitor forest treatment success. To support their needs, the NASA DEVELOP Front Range Wildland Fires team explored National Agricultural Imagery Program (NAIP) imagery at different spatial resolutions and numbers of training points with NASA’s Shuttle Radar Topography Mission (SRTM) Data Elevation Model (DEM) as a predictor in addition to NAIP imagery spectral predictors. From this analysis, we created classified canopy cover rasters, and compared accuracy metrics across model iterations. We also determined that the best performing model, with an overall accuracy of 0.900 uses 2021 NAIP imagery at 2-meter resolution, 800 training points, 200 testing points, does not use topographic predictors, and reclassifies shadow pixels via a pre-selected NDVI threshold.

Remote Sensing↗

Adaptive time scaling for accelerating explicit finite element analysis

A method for accelerating an explicit finite element analysis (FEA) simulation of a modeled system or process includes performing an initial iteration of the FEA simulation according to a baseline time interval via an FEA computing network, and calculating a criteria ratio of a predetermined set of scaling criteria for the modeled system or process. The method includes determining a time-scaling factor using the criteria ratio via the FEA computing network as a function of the criteria ratio, and then applying the time-scaling factor to the baseline time interval to generate a scaled time interval. The scaled time interval accelerates simulation time of the FEA simulation. The method includes performing a subsequent iteration of the explicit FEA simulation at the scaled time interval using the FEA computing network. The process continues for subsequent iterations, with the time-scaling factor adapting with each iteration.

Chen, Jian↗

Efficient thermo-mechanical simulation approaches for arc-and-wire based additive manufacturing

This study aims at the thermo-mechanical features inherent with the arc and wire metal deposition process to advance the structural simulation of additive manufacturing. Simulation approaches including iterative substructure method (ISM), dynamic mesh refining method (DMRM), and graphics processing unit (GPU) based explicit FEM code (GPU code) were developed for accelerating stress analysis and compared in terms of efficiency and accuracy. ISM and DMRM both take advantage of the strong nonlinearity phenomenon near the moving heat source, by reducing the global iterations and refining the local mesh respectively. The explicit FEM code was accelerated by a mass scaling technique based on the inherent strain concept. Two large parts by wire and arc additive manufacturing were analyzed, and the numerical results agreed very well with each other. The predicted distortion and residual stresses in the deposited layers were validated by experimental measurement. Compared with commercial software ABAQUS, the novel approaches reduced the computation cost substantially without compromising on accuracy. GPU code shows the highest computation efficiency (x30~70), while DMRM and ISM has an acceleration factor of x9 and x3 respectively. Owing to the minimum model simplification, such advanced simulation approaches will be very useful for structural analysis of additive manufacturing with complex deposition paths and product geometry.

Huang, Hui↗

TRANSP Tests Of TGLF and Predictions For ITER

Gyro kinetic simulations of turbulence capture some of the features observed in transport, fluctuations, and correlations measured in tokamak plasmas. These codes calculations are CPU intensive, and are not practical for incorporation in present time-dependant transport codes, so reduced models based on these gyro kinetic codes are being used. An example is the TGLF model [1] which is a quasilinear gyrofluid model calibrated to nonlinear results from the GYRO code [2]. Recently TGLF has been incorporated into TRANSP. Analysis of experimental data using TRANSP with such models provides fundamental understanding of turbulent transport. Predictions of ITER performance with various plasma scenarios using such models are useful for optimizing design and for exposing issues that can be addressed in present experiments and theory. For instance, which combinations of heating, torquing, and current drive are optimal. Another application is for nuclear licensing (e.g. system integrity, neutron rates). Others are generating inputs for design of diagnostic systems and for theoretical studies. An example of the later is Alfv´en Eigenmode and AE-induced loss of fast ions. The beam ion distribution can either enhance or reduce the alpha pressure drive of the AE instability. The AE instability can cause dangerous amounts of fast ion losses, as was seen in TFTR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Advanced applications of BEM to inelastic analysis of solids

Newly evolved boundary element solution algorithms are utilized for solving a range of elastoplastic and thermoplastic problems. The first algorithm is the typical iterative procedure in which the unknown boundary solution and the initial stress (or strain) rates are found together in an incremental iterative fashion. The second one is a variable stiffness type approach in which the incremental boundary solution is obtained in a direct (noniterative) manner. The third procedure differs from the two prior algorithms in that volume integration is not needed to incorporate the nonlinear effects in the analysis. Instead, initial stress rates are introduced in the boundary element system by way of particular integrals. The methodology is used in a general purpose, multiregion system that employs quadratic isoparametric shape functions to model the geometry and field variables of the body and can admit up to 15 substructured areas of different material properties.

Banerjee, P. K.↗

CT head-scan dosimetry in an anthropomorphic phantom and associated measurement of ACR accreditation-phantom imaging metrics under clinically representative scan conditions

Purpose: To measure radiation absorbed dose and its distribution in an anthropomorphic head phantom under clinically representative scan conditions in three widely used computed tomography (CT) scanners, and to relate those dose values to metrics such as high-contrast resolution, noise, and contrast-to-noise ratio (CNR) in the American College of Radiology CT accreditation phantom.Methods: By inserting optically stimulated luminescence dosimeters (OSLDs) in the head of an anthropomorphic phantom specially developed for CT dosimetry (University of Florida, Gainesville), we measured dose with three commonly used scanners (GE Discovery CT750 HD, Siemens Definition, Philips Brilliance 64) at two different clinical sites (Walter Reed National Military Medical Center, National Institutes of Health). The scanners were set to operate with the same data-acquisition and image-reconstruction protocols as used clinically for typical head scans, respective of the practices of each facility for each scanner. We also analyzed images of the ACR CT accreditation phantom with the corresponding protocols. While the Siemens Definition and the Philips Brilliance protocols utilized only conventional, filtered back-projection (FBP) image-reconstruction methods, the GE Discovery also employed its particular version of an adaptive statistical iterative reconstruction (ASIR) algorithm that can be blended in desired proportions with the FBP algorithm. We did an objective image-metrics analysis evaluating the modulation transfer function (MTF), noise power spectrum (NPS), and CNR for images reconstructed with FBP. For images reconstructed with ASIR, we only analyzed the CNR, since MTF and NPS results are expected to depend on the object for iterative reconstruction algorithms.Results: The OSLD measurements showed that the Siemens Definition and the Philips Brilliance scanners (located at two different clinical facilities) yield average absorbed doses in tissue of 42.6 and 43.1 mGy, respectively. The GE Discovery delivers about the same amount of dose (43.7 mGy) when run under similar operating and image-reconstruction conditions, i.e., without tube current modulation and ASIR. The image-metrics analysis likewise showed that the MTF, NPS, and CNR associated with the reconstructed images are mutually comparable when the three scanners are run with similar settings, and differences can be attributed to different edge-enhancement properties of the applied reconstruction filters. Moreover, when the GE scanner was operated with the facility's scanner settings for routine head exams, which apply 50% ASIR and use only approximately half of the 100%-FBP dose, the CNR of the images showed no significant change. Even though the CNR alone is not sufficient to characterize the image quality and justify any dose reduction claims, it can be useful as a constancy test metric.Conclusions: This work presents a straightforward method to connect direct measurements of CT dose with objective image metrics such as high-contrast resolution, noise, and CNR. It demonstrates that OSLD measurements in an anthropomorphic head phantom allow a realistic and locally precise estimation of magnitude and spatial distribution of dose in tissue delivered during a typical CT head scan. Additional objective analysis of the images of the ACR accreditation phantom can be used to relate the measured doses to high contrast resolution, noise, and CNR.

61 RADIATION PROTECTION AND DOSIMETRY↗

Coupled model for liquid lithium plasma facing components

Numerical analysis provides the design choice and operating window of liquid metal Plasma Facing Components (PFC) concepts. Coupled analysis of boundary plasma together with the surrounding boundary structures is required. Here, to achieve this goal, PPPL is developing a comprehensive multi-physics model for modeling of PFCs in fusion devices. The model includes the fluid-kinetic code SOLPS-ITER and the flow and heat transfer code CFX from ANSYS. SOLPS-ITER was augmented with a liquid metal boundary condition algorithm, allowing direct two-way coupling of the plasma analysis with the two-dimensional analytical slab flow model which includes heat convection in the liquid metal PFC. The target heat flux resulting from this coupled analysis is used as a boundary condition for detailed 3D Computational Fluid Dynamics (CFD) Magneto Hydro Dynamics (MHD) and heat transfer analysis. A new formulation of MHD equations is introduced in the numerical procedure ensuring current conservation of the discretized equations. Results of the 3D analysis are used for final validation of the coupled model. A PFC design where a porous wall is used to stabilize the liquid metal surface, while MHD drive is used to push the liquid metal flow inside the PFC, will be investigated in the regimes where vapor shielding is created for enhanced volumetric plasma heat dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Diverse Power Iteration Embeddings and Its Applications

Abstract—Spectral Embedding is one of the most effective dimension reduction algorithms in data mining. However, its computation complexity has to be mitigated in order to apply it for real-world large scale data analysis. Many researches have been focusing on developing approximate spectral embeddings which are more efficient, but meanwhile far less effective. This paper proposes Diverse Power Iteration Embeddings (DPIE), which not only retains the similar efficiency of power iteration methods but also produces a series of diverse and more effective embedding vectors. We test this novel method by applying it to various data mining applications (e.g. clustering, anomaly detection and feature selection) and evaluating their performance improvements. The experimental results show our proposed DPIE is more effective than popular spectral approximation methods, and obtains the similar quality of classic spectral embedding derived from eigen-decompositions. Moreover it is extremely fast on big data applications. For example in terms of clustering result, DPIE achieves as good as 95% of classic spectral clustering on the complex datasets but 4000+ times faster in limited memory environment.

97 MATHEMATICS AND COMPUTING↗

An Online Plug-and-Play Algorithm for Regularized Image Reconstruction

Plug-and-play priors (PnP) is a powerful framework for regularizing imaging inverse problems by using advanced denoisers within an iterative algorithm. Recent experimental evidence suggests that PnP algorithms achieve state-of-the-art performance in a range of imaging applications. Here, we introduce a new online PnP algorithm based on the proximal gradient method (PGM). The proposed algorithm uses only a subset of measurements at every iteration, which makes it scalable to very large datasets. We present a new theoretical convergence analysis, for both batch and online variants of PnP-PGM, for denoisers that do not necessarily correspond to proximal operators. We also present simulations illustrating the applicability of the algorithm to image reconstruction in diffraction tomography. Finally, the results in this paper have the potential to expand the applicability of the PnP framework to very large datasets.

97 MATHEMATICS AND COMPUTING↗

Ion beam analysis of fusion plasma-facing materials and components: facilities and research challenges

Following the IAEA Technical Meeting on 'Advanced Methodologies for the Analysis of Materials in Energy Applications Using Ion Beam Accelerators', our work reviews the current status of ion beam analysis (IBA) techniques and some aspects of ion-induced radiation damage in materials for the field of materials relevant to fusion. Available facilities, apparatus development, future research options and challenges are presented and discussed. The analysis of beryllium and radioactivity-containing samples from future experiments in JET or ITER represents not only an analytical but also a technical challenge. A comprehensive list of the facilities, their current status, and analytical capabilities comes alongside detailed descriptions of the labs. A discussion of future issues of sample handling and the current status of facilities at JET complete the technical section. To prepare the international IBA community for these challenges, the IAEA technical meeting finalizes the necessity for determining new nuclear reaction cross-sections and improving the inter-laboratory comparability by defining international standards and testing these via a round-robin test.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Self-driving thin film laboratory: autonomous epitaxial atomic-layer synthesis via real-time computer vision analysis of electron diffraction

Emerging materials science platforms with the ability to make autonomous decisions on the fly are fundamentally changing the outlook and protocols for materials optimization and discovery. Because AI-driven self-navigating schemes can effectively reduce the total number of iterations needed to arrive at the "answer" (i.e. the best stochiometric composition for a desired physical property, optimum materials processing parameters, etc.) by significant margins, they have the potential to revolutionize materials and chemical manufacturing processes at large in research laboratory settings as well as in industrial plants. Here, we demonstrate a successful implementation of real-time closed-loop autonomous navigation of a multi-dimensional materials synthesis parameter space for fabricating phase-pure epitaxial films of a metastable phase of a functional oxide in a combinatorial pulsed laser deposition chamber. Sequential epitaxial growth iterations in search of the optimized recipe to stabilize the desired crystal phase were performed using frame-by-frame quantitative computer vision analysis of reflection high-energy electron diffraction (RHEED) images of the unit-cell level film being deposited. The autonomous scheme regularly resulted in > 30-fold reduction in the number of required experiments compared to a comprehensive mapping of the parameter space. The real-time workflow developed here can be readily extended to a variety of thin film synthesis platforms opening the door for self-driving atomic-level materials design as well as autonomous optimization of semiconductor manufacturing.

36 MATERIALS SCIENCE↗

SU-D-12A-06: A Comprehensive Parameter Analysis for Low Dose Cone-Beam CT Reconstruction

Purpose: There is always a parameter in compressive sensing based iterative reconstruction (IR) methods low dose cone-beam CT (CBCT), which controls the weight of regularization relative to data fidelity. A clear understanding of the relationship between image quality and parameter values is important. The purpose of this study is to investigate this subject based on experimental data and a representative advanced IR algorithm using Tight-frame (TF) regularization. Methods: Three data sets of a Catphan phantom acquired at low, regular and high dose levels are used. For each tests, 90 projections covering a 200-degree scan range are used for reconstruction. Three different regions-of-interest (ROIs) of different contrasts are used to calculate contrast-to-noise ratios (CNR) for contrast evaluation. A single point structure is used to measure modulation transfer function (MTF) for spatial-resolution evaluation. Finally, we analyze CNRs and MTFs to study the relationship between image quality and parameter selections. Results: It was found that: 1) there is no universal optimal parameter. The optimal parameter value depends on specific task and dose level. 2) There is a clear trade-off between CNR and resolution. The parameter for the best CNR is always smaller than that for the best resolution. 3) Optimal parameters are also dose-specific. Data acquired under a high dose protocol require less regularization, yielding smaller optimal parameter values. 4) Comparing with conventional FDK images, TF-based CBCT images are better under a certain optimally selected parameters. The advantages are more obvious for low dose data. Conclusion: We have investigated the relationship between image quality and parameter values in the TF-based IR algorithm. Preliminary results indicate optimal parameters are specific to both the task types and dose levels, providing guidance for selecting parameters in advanced IR algorithms. This work is supported in part by NIH (1R01CA154747-01)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Uncertainty propagation in a stratospheric model. I - Development of a concise stratospheric model. II - Monte Carlo analysis of imprecisions due to reaction rates

The paper discusses the development of a concise stratospheric model which uses iteration to obtain coupling between interacting species. The one-dimensional, steady-state, diurnally-averaged model generates diffusion equations with appropriate sources and sinks for species odd oxygen, H2O, H2, CO, N2O, odd nitrogen, CH4, CH3Cl, CCl4, CF2Cl2, CFCl3, and odd chlorine. The model evaluates steady-state perturbations caused by injections of chlorine and NO(x) and may be used to predict ozone depletion. The model is used in a Monte Carlo study of the propagation of reaction-rate imprecisions by calculating an ozone perturbation caused by the addition of chlorine. Since the model is sensitive to only 10 of the more than 50 reaction rates considered, only about 1000 Monte Carlo cases are required to span the space of possible results.

Rundel, R. D.↗