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 775 records · Page 43

A Performant, Scalable Processing Pipeline for High‐Quality and FAIR Environmental Sensor Data

High-resolution environmental monitoring is necessary to record, understand, and predict biogeochemical and ecological changes particularly in coastal systems but brings significant challenges in processing and making rapidly available the resulting data. The COMPASS-FME project established a network of coastal observational sites across the Chesapeake Bay and western Lake Erie regions extensively instrumented with soil, vegetation, and weather sensors logging data every 15 min. Our data processing framework, written in R and completely open source, prioritizes rapid model-experiment iteration and makes biogeochemical data rapidly available for quality assurance/quality control, analysis, and model ingestion. This pipeline is distinguished by a standardized and modular approach to data curation, extensive metadata and documentation, and its high performance. These attributes combine to make biogeochemical data rapidly accessible across COMPASS-FME and the broader community. Flexible, powerful, and reproducible approaches to handling high-volume environmental data are crucial for accelerating biogeosciences research.

Pennington, Stephanie C. [Pacific Northwest Nation↗

Exact and inexact subsampled Newton methods for optimization

Abstract The paper studies the solution of stochastic optimization problems in which approximations to the gradient and Hessian are obtained through subsampling. We first consider Newton-like methods that employ these approximations and discuss how to coordinate the accuracy in the gradient and Hessian to yield a superlinear rate of convergence in expectation. The second part of the paper analyzes an inexact Newton method that solves linear systems approximately using the conjugate gradient (CG) method, and that samples the Hessian and not the gradient (the gradient is assumed to be exact). We provide a complexity analysis for this method based on the properties of the CG iteration and the quality of the Hessian approximation, and compare it with a method that employs a stochastic gradient iteration instead of the CG method. We report preliminary numerical results that illustrate the performance of inexact subsampled Newton methods on machine learning applications based on logistic regression.

Mathematics↗

SU-F-T-575: Verification of a Monte-Carlo Small Field SRS/SBRT Dose Calculation System

Purpose: To develop a methodology for validation of a Monte-Carlo dose calculation model for robotic small field SRS/SBRT deliveries. Methods: In a robotic treatment planning system, a Monte-Carlo model was iteratively optimized to match with beam data. A two-part analysis was developed to verify this model. 1) The Monte-Carlo model was validated in a simulated water phantom versus a Ray-Tracing calculation on a single beam collimator-by-collimator calculation. 2) The Monte-Carlo model was validated to be accurate in the most challenging situation, lung, by acquiring in-phantom measurements. A plan was created and delivered in a CIRS lung phantom with film insert. Separately, plans were delivered in an in-house created lung phantom with a PinPoint chamber insert within a lung simulating material. For medium to large collimator sizes, a single beam was delivered to the phantom. For small size collimators (10, 12.5, and 15mm), a robotically delivered plan was created to generate a uniform dose field of irradiation over a 2×2cm{sup 2} area. Results: Dose differences in simulated water between Ray-Tracing and Monte-Carlo were all within 1% at dmax and deeper. Maximum dose differences occurred prior to dmax but were all within 3%. Film measurements in a lung phantom show high correspondence of over 95% gamma at the 2%/2mm level for Monte-Carlo. Ion chamber measurements for collimator sizes of 12.5mm and above were within 3% of Monte-Carlo calculated values. Uniform irradiation involving the 10mm collimator resulted in a dose difference of ∼8% for both Monte-Carlo and Ray-Tracing indicating that there may be limitations with the dose calculation. Conclusion: We have developed a methodology to validate a Monte-Carlo model by verifying that it matches in water and, separately, that it corresponds well in lung simulating materials. The Monte-Carlo model and algorithm tested may have more limited accuracy for 10mm fields and smaller.

60 APPLIED LIFE SCIENCES↗

Perspectives on AI Architectures and Codesign for Earth System Predictability

Abstract Recently, the U.S. Department of Energy (DOE), Office of Science, Biological and Environmental Research (BER), and Advanced Scientific Computing Research (ASCR) programs organized and held the Artificial Intelligence for Earth System Predictability (AI4ESP) workshop series. From this workshop, a critical conclusion that the DOE BER and ASCR community came to is the requirement to develop a new paradigm for Earth system predictability focused on enabling artificial intelligence (AI) across the field, laboratory, modeling, and analysis activities, called model experimentation (ModEx). BER’s ModEx is an iterative approach that enables process models to generate hypotheses. The developed hypotheses inform field and laboratory efforts to collect measurement and observation data, which are subsequently used to parameterize, drive, and test model (e.g., process based) predictions. A total of 17 technical sessions were held in this AI4ESP workshop series. This paper discusses the topic of the AI Architectures and Codesign session and associated outcomes. The AI Architectures and Codesign session included two invited talks, two plenary discussion panels, and three breakout rooms that covered specific topics, including 1) DOE high-performance computing (HPC) systems, 2) cloud HPC systems, and 3) edge computing and Internet of Things (IoT). We also provide forward-looking ideas and perspectives on potential research in this codesign area that can be achieved by synergies with the other 16 session topics. These ideas include topics such as 1) reimagining codesign, 2) data acquisition to distribution, 3) heterogeneous HPC solutions for integration of AI/ML and other data analytics like uncertainty quantification with Earth system modeling and simulation, and 4) AI-enabled sensor integration into Earth system measurements and observations. Such perspectives are a distinguishing aspect of this paper. Significance Statement This study aims to provide perspectives on AI architectures and codesign approaches for Earth system predictability. Such visionary perspectives are essential because AI-enabled model-data integration has shown promise in improving predictions associated with climate change, perturbations, and extreme events. Our forward-looking ideas guide what is next in codesign to enhance Earth system models, observations, and theory using state-of-the-art and futuristic computational infrastructure.

54 ENVIRONMENTAL SCIENCES↗

Perspectives on AI Architectures and Co-design for Earth System Predictability

Recently, the U.S. Department of Energy (DOE), Office of Science, Biological and Environmental Research (BER), and Advanced Scientific Computing Research (ASCR) programs organized and held the Artificial Intelligence for Earth System Predictability (AI4ESP) workshop series. From this workshop, a critical conclusion that the DOE BER and ASCR community came to is the requirement to develop a new paradigm for Earth system predictability focused on enabling artificial intelligence (AI) across the field, laboratory, modeling, and analysis activities, called model experimentation (ModEx). BER’s ModEx is an iterative approach that enables process models to generate hypotheses. The developed hypotheses inform field and laboratory efforts to collect measurement and observation data, which are subsequently used to parameterize, drive, and test model (e.g., process based) predictions. A total of 17 technical sessions were held in this AI4ESP workshop series. This paper discusses the topic of the AI Architectures and Codesign session and associated outcomes. The AI Architectures and Codesign session included two invited talks, two plenary discussion panels, and three breakout rooms that covered specific topics, including 1) DOE high-performance computing (HPC) systems, 2) cloud HPC systems, and 3) edge computing and Internet of Things (IoT). We also provide forward-looking ideas and perspectives on potential research in this codesign area that can be achieved by synergies with the other 16 session topics. These ideas include topics such as 1) reimagining codesign, 2) data acquisition to distribution, 3) heterogeneous HPC solutions for integration of AI/ML and other data analytics like uncertainty quantification with Earth system modeling and simulation, and 4) AI-enabled sensor integration into Earth system measurements and observations. Such perspectives are a distinguishing aspect of this paper.

58 GEOSCIENCES↗

Reassessing energy deposition for the ITER 5 MA vertical displacement event with an improved DINA model

The beryllium (Be) main chamber wall interaction during a 5 MA/1.8 T upward, unmitigated VDE scenario, first analysed in [J. Coburn et al., Phys. Scr. T171 (2020) 014076] for ITER, has been re-evaluated using the latest energy deposition analysis software. Updates to the DINA disruption model are summarized, including an improved numerical convergence for the 0D power balance, limitations on the safety factor within the plasma core, and the choice to maintain a constant plasma +halo poloidal cross-section. Such updates result in a broad halo region and higher radiated power fractions compared to previous models. The new scenario lasts for ~75 ms and deposits ~29 MJ of energy, with the radial distribution of parallel heat flux q‖(r)resembling an exponential falloff with an effective λ_q=75-198 mm. A maximum halo width w_h of 0.52 m at the outboard midplane is observed. SMITER field line tracing and energy deposition simulations calculate a q_{⊥,max} of ~83 MW/m^2 on the upper first wall panels (FWP). Heat transfer calculations with the MEMOS-U code show that the FWP surface temperature reaches ~1000 K, well below the Be melt threshold. Variations of this 5 MA scenario with Be im-purity densities from 0 to 3∙10^19 m^-3 also remain below the melt threshold despite differences in energy deposition and duration. These results are in contrast to the early study which predicted melt damage to the first wall [J. Coburn et al., Phys. Scr. T171 (2020) 014076], and emphasize the importance of accurate models for the halo width w_h and the heat flux distribution q‖(r)within that halo width. The 2020 halo model in DINA has been compared with halo current experiments on COMPASS, JET, and Alcator C-Mod, and the preliminary results build confidence in the broad halo width predictions. Results for the 5 MA VDE are compared with those for a 15 MA equivalent, generated using the new DINA model. At the higher current, significant melting of the upper FWP is to be expected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Forward marching procedure for separated boundary-layer flows

A forward-marching procedure for separated boundary-layer flows which permits the rapid and accurate solution of flows of limited extent is presented. The streamwise convection of vorticity in the reversed flow region is neglected, and this approximation is incorporated into a previously developed (Carter, 1974) inverse boundary-layer procedure. The equations are solved by the Crank-Nicolson finite-difference scheme in which column iteration is carried out at each streamwise station. Instabilities encountered in the column iterations are removed by introducing timelike terms in the finite-difference equations. This provides both unconditional diagonal dominance and a column iterative scheme, found to be stable using the von Neumann stability analysis.

Carter, J. E.↗

Numerical solution of nonlinear partial differential equations of mixed type

A review is presented of some recently developed numerical methods for the solution of nonlinear equations of mixed type. The methods considered use finite difference approximations to the differential equation. Central difference formulas are employed in the subsonic zone and upwind difference formulas are used in the supersonic zone. The relaxation method for the small disturbance equation is discussed and a description is given of difference schemes for the potential flow equation in quasi-linear form. Attention is also given to difference schemes for the potential flow equation in conservation form, the analysis of relaxation schemes by the time dependent analogy, the accelerated iterative method, and three-dimensional calculations.

Jameson, A.↗

Simplified contact analysis

A simple and straightforward method for calculating the elastic deformation at the center of contact between two solids with different radii of curvature is proposed instead of the conventional contact deformation analysis for ball bearings, gears, and cams, which usually involves tedious iterative procedures or the use of design charts. Ellipticity of the contact is approximated from a least-squares power fit. Simplified expressions are derived which allow rapid calculation of deformation from the material properties and geometry of the contacting elements.

Brewe, D. E.↗

Direct solution for thermal stresses in a nose cap under an arbitrary axisymmetric temperature distribution

The design of a nose cap for a hypersonic vehicle is an iterative process requiring a rapid, easy to use and accurate stress analysis. The objective of this paper is to develop such a stress analysis technique from a direct solution of the thermal stress equations for a spherical shell. The nose cap structure is treated as a thin spherical shell with an axisymmetric temperature distribution. The governing differential equations are solved by expressing the stress solution to the thermoelastic equations in terms of a series of derivatives of the Legendre polynomials. The process of finding the coefficients for the series solution in terms of the temperature distribution is generalized by expressing the temperature along the shell and through the thickness as a polynomial in the spherical angle coordinate. Under this generalization the orthogonality property of the Legendre polynomials leads to a sequence of integrals involving powers of the spherical shell coordinate times the derivative of the Legendre polynomials. The coefficients of the temperature polynomial appear outside of these integrals. Thus, the integrals are evaluated only once and their values tabulated for use with any arbitrary polynomial temperature distribution.

Davis, Randall C.↗

Deducing Electronic Unit Internal Response During a Vibration Test Using a Lumped Parameter Modeling Approach

During random vibration testing of electronic boxes there is often a desire to know the dynamic response of certain internal printed wiring boards (PWBs) for the purpose of monitoring the response of sensitive hardware or for post-test forensic analysis in support of anomaly investigation. Due to restrictions on internally mounted accelerometers for most flight hardware there is usually no means to empirically observe the internal dynamics of the unit, so one must resort to crude and highly uncertain approximations. One common practice is to apply Miles Equation, which does not account for the coupled response of the board in the chassis, resulting in significant over- or under-prediction. This paper explores the application of simple multiple-degree-of-freedom lumped parameter modeling to predict the coupled random vibration response of the PWBs in their fundamental modes of vibration. A simple tool using this approach could be used during or following a random vibration test to interpret vibration test data from a single external chassis measurement to deduce internal board dynamics by means of a rapid correlation analysis. Such a tool might also be useful in early design stages as a supplemental analysis to a more detailed finite element analysis to quickly prototype and analyze the dynamics of various design iterations. After developing the theoretical basis, a lumped parameter modeling approach is applied to an electronic unit for which both external and internal test vibration response measurements are available for direct comparison. Reasonable correlation of the results demonstrates the potential viability of such an approach. Further development of the preliminary approach presented in this paper will involve correlation with detailed finite element models and additional relevant test data.

random vibration response↗

Digital Human Modeling

The development of models to represent human characteristics and behaviors in human factors is broad and general. The term "model" can refer to any metaphor to represent any aspect of the human; it is generally used in research to mean a mathematical tool for the simulation (often in software, which makes the simulation digital) of some aspect of human performance and for the prediction of future outcomes. This section is restricted to the application of human models in physical design, e.g., in human factors engineering. This design effort is typically human interface design, and the digital models used are anthropometric. That is, they are visual models that are the physical shape of humans and that have the capabilities and constraints of humans of a selected population. They are distinct from the avatars used in the entertainment industry (movies, video games, and the like) in precisely that regard: as models, they are created through the application of data on humans, and they are used to predict human response; body stresses workspaces. DHM enable iterative evaluation of a large number of concepts and support rapid analysis, as compared with use of physical mockups. They can be used to evaluate feasibility of escape of a suited astronaut from a damaged vehicle, before launch or after an abort (England, et al., 2012). Throughout most of human spaceflight, little attention has been paid to worksite design for ground workers. As a result of repeated damage to the Space Shuttle which adversely affected flight safety, DHM analyses of ground assembly and maintenance have been developed over the last five years for the design of new flight systems (Stambolian, 2012, Dischinger and Dunn Jackson, 2014). The intent of these analyses is to assure the design supports the work of the ground crew personnel and thereby protect the launch vehicle. They help the analyst address basic human factors engineering questions: can a worker reach the task site from the work platform provided; can she or he see the task site; can she or he control tools, which, if dropped, might damage the system? Figure 7.3.1 provides an example of such analysis for a future NASA launch vehicle. [figure 7.3.1 here] In-space systems for operation by astronauts have long been targets for DHM analysis, given the focus on mission success and concerns for astronaut safety. Figure 7.3.2 illustrates the analysis of the design to support astronaut tasks for an International Space Station glovebox. [Figure 7.3.2 here] Use by

Dischinger, H. Charles, Jr.↗

The Effects of Flux Spectrum Perturbation on Transmutation of Actinides: Optimizing the Production of Transcurium Isotopes

This research presented herein involves the optimization of transcurium production in the High Flux Isotope Reactor at Oak Ridge National Laboratory. Due to the dependence of isotope cross sections on incoming neutron energy, the efficiency with which an isotope is transmuted is highly dependent upon the flux spectrum. There are certain energy bands in which the rate of fission of transcurium production feedstock materials is minimized, relative to the rate of non-fission absorptions. It is proposed that by perturbing the flux spectrum, it is possible to increase the amount of key isotopes, such as 249Bk and 252Cf, that are produced during a transmutation cycle, relative to the consumption of feedstock material. This optimization process is carried out by developing an iterative objective framework involving problem definition, flux spectrum and cross section analysis, simulated transmutation, and analysis of final yields and transmutation parameters. It is shown that it is possible to perturb the local flux spectrum in the transcurium target by perturbing the composition of the target. It is further shown that these perturbations are able to alter the target yields in a non-negligible way. Future work is necessary to develop the optimization framework, and identify the necessary algorithms to update the problem definition based upon progress towards the optimization goals.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

URANIUM ACCUMULATIONS IN CASTING OPERATIONS

In December 2017, unexpected accumulations of uranium were discovered in several locations throughout the casting process area at the Y-12 National Security Complex that subsequently placed casting operations on hold. The Y-12 uranium casting process takes place in enclosures and hoods, many of which are open to a large floor surface area that has limited visibility and access. At the time of discovery, there were no Nuclear Criticality Safety (NCS) analyses or limits for material accumulations in the casting processing areas. Due to changes in the operation of the casting process and invalid assumptions in the criticality safety evaluation, significant material accumulations were occurring inside the enclosures. The material accumulations were identified as a result of an extent of condition review from another process, and an occurrence report was filed. The NCS organization worked with Production and Operations personnel to identify additional areas of uranium accumulation in the casting area, which resulted in several iterations of cleaning and material collection. Following material recovery, a causal analysis was performed that identified several contributing factors that allowed the accumulations to go unnoticed. In order to resume operations, the criticality safety evaluation was revised to analyze and control uranium accumulations in the casting area. Additionally, improvements in nuclear material control and accountability (NMC&A) controls were required. Casting operations have since resumed with frequently required material collections, which have limited the amount of uranium material accumulation.

42 ENGINEERING↗

URANIUM ACCUMULATIONS IN CASTING OPERATIONS

In December 2017, unexpected accumulations of uranium were discovered in several locations throughout the casting process area at the Y-12 National Security Complex that subsequently placed casting operations on hold. The Y-12 uranium casting process takes place in enclosures and hoods, many of which are open to a large floor surface area that has limited visibility and access. At the time of discovery, there were no Nuclear Criticality Safety (NCS) analyses or limits for material accumulations in the casting processing areas. Due to changes in the operation of the casting process and invalid assumptions in the criticality safety evaluation, significant material accumulations were occurring inside the enclosures. The material accumulations were identified as a result of an extent of condition review from another process, and an occurrence report was filed. The NCS organization worked with Production and Operations personnel to identify additional areas of uranium accumulation in the casting area, which resulted in several iterations of cleaning and material collection. Following material recovery, a causal analysis was performed that identified several contributing factors that allowed the accumulations to go unnoticed. In order to resume operations, the criticality safety evaluation was revised to analyze and control uranium accumulations in the casting area. Additionally, improvements in nuclear material control and accountability (NMC&A) controls were required. Casting operations have since resumed with frequently required material collections, which have limited the amount of uranium material accumulation.

61 RADIATION PROTECTION AND DOSIMETRY↗

Error field predictability and consequences for ITER

Abstract ITER coil tolerances are re-evaluated using the modern understanding of coupling to least-stable plasma modes and an updated center-line-traced model of ITER’s coil windings. This reassessment finds the tolerances to be conservative through a statistical, linear study of n = 1 error fields (EFs) due to tilted, shifted misplacements and nominal windings of central solenoid and poloidal field coils within tolerance. We also show that a model-based correction scheme remains effective even when metrology quality is sub-optimal, and compare this to projected empirical correction schemes. We begin with an analysis of the necessity of error field correction (EFC) for daily operation in ITER using scalign laws for the EF penetration threshold. We then consider the predictability of EF dominant mode overlap across early planned ITER scenarios and, as measuring EFs in high power scenarios can pose risks to the device, the potential for extrapolation to the ITER Baseline Scenario (IBS). We find that carefully designing a scenario matching currents proportionally to those of the IBS is far more important than plasma shape or profiles in accurately measuring an optimal correction current set.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mathematical Model for the Emission Infrared Tomography of the Temperature Field in an Isotropic Layer

We consider a mathematical model for the determination of the temperature field in a layer that emits, absorbs, and scatters infrared radiation both in the volume and on the surface. Within the framework of the model, we formulate nonlinear direct and inverse problems of emission infrared tomography of the temperature field in the layer according to known space-and-angular distributions of the flows of infrared radiation emitted into the ambient medium. The iterative-variational methods are proposed for the solution of the indicated direct and inverse problems. As a specific example, we perform the numerical analysis of the developed algorithms.

97 MATHEMATICS AND COMPUTING↗