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 361 records · Page 20

Low-Rank Incremental Methods for Computing Dominant Singular Subspaces

Computing the singular values and vectors of a matrix is a crucial kernel in numerous scientific and industrial applications. As such, numerous methods have been proposed to handle this problem in a computationally efficient way. This paper considers a family of methods for incrementally computing the dominant SVD of a large matrix A. Specifically, we describe a unification of a number of previously disparate methods for approximating the dominant SVD via a single pass through A. We tie the behavior of these methods to that of a class of optimization-based iterative eigensolvers on A'*A. An iterative procedure is proposed which allows the computation of an accurate dominant SVD via multiple passes through A. We present an analysis of the convergence of this iteration, and provide empirical demonstration of the proposed method on both synthetic and benchmark data.

99 GENERAL AND MISCELLANEOUS↗

Statistical analysis of the Nb 3 Sn strand production for the ITER toroidal field coils

The ITER toroidal field (TF) strand procurement initiated the largest Nb 3 Sn superconducting strand production hitherto. The industrial-scale production started in Japan in 2008 and finished in summer 2015. Six ITER partners (so-called Domestic Agencies, or DAs) are in charge of the procurement and involved eight different strand suppliers all over the world, of which four are using the bronze route (BR) process and four the internal-tin (IT) process. In total more than 500 tons have been produced including excess material covering losses during the conductor manufacturing process, in particular the cabling. The procurement is based on a functional specification where the main strand requirements like critical current, hysteresis losses, Cu ratio and residual resistance ratio are specified but not the strand production process or layout. We introduce the analysis on the data acquired during the quality control (QC) process that was carried out to ensure the same conductor performance requirements are met by the different strand suppliers regardless of strand design. The strand QC is based on 100% billet testing and on applying statistical process control (SPC) limits. Throughout the production, samples adjacent to the strand pieces tested by the suppliers are cross-checked ('verified') by their respective DAs reference labs. The level of verification was lowered from 100% at the beginning of the procurement progressively to approximately 25% during the final phase of production. Based on the complete dataset of the TF strand production, an analysis of the SPC limits of the critical strand parameters is made and the related process capability indices are calculated. In view of the large-scale production and costs, key manufacturing parameters such as billet yield, number of breakages and piece-length distribution are also discussed. Furthermore, the results are compared among all the strand suppliers, focusing on the difference between BR and IT processes. Following the completion of the largest Nb 3 Sn strand production, our experience gained from monitoring the execution of the QC activities and from auditing the results from the measurements is summarised for future superconducting strand material procurement activities.

Vostner, A.↗

ITER diagnostic systems in development in Ioffe Institute

Three diagnostic systems are being developed in Ioffe Institute for ITER. Those are Neutral Particle Analysis (NPA), Thomson Scattering in Divertor (TSD) and Gamma Spectroscopy (GS). The main objective of NPA in ITER is to measure D/T fuel ration in plasma on the basis of measurement of neutralized fluxes of D and T ions [1]. Fuel ratio is one of the key parameters needed by ITER control system to provide the optimal conditions in plasma and the most effective plasma burning. Another objective is to measure the distribution function of fast ions (including alpha particles) generated as a result of the additional heating and nuclear fusion reactions. Thomson Scattering in Divertor (TSD) [2] will be used to measure electron temperature and density in the scrape-off layer in outer leg of ITER divertor. The main task of TSD is to protect the machine from divertor overloading. Gamma Spectroscopy (GS) [3] is based on the measurement of spectral lines of MeV range gammas generated in nuclear reactions in plasma. 2-D gamma-ray emission measurements give valuable information on the confined alpha particles in DT plasma. They also provide important information on the location of MeV range runaway electron beams in ITER plasma. For all three cases the physical basis and instrumentation are presented. The simple NPA version for measurements of D/T ratio in DEMO is also briefly described.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A diffusion synthetic acceleration scheme for rectangular geometries based on bilinear discontinuous finite elements

A DSA technique to accelerate the iterative convergence of S{sub n} transport solves is derived for bilinear discontinuous (BLD) finite elements on rectangular grids. The diffusion synthetic acceleration equations are discretized using BLD elements by adapting the Modified Interior Penalty technique, introduced in [4] for triangular grids. The MIP-DSA equations are SPD and thus are solved using a preconditioned CG technique. Fourier analyses and implementation of the technique in a BLD S{sub n} transport code show that the technique is stable is effective. (authors)

97 MATHEMATICS AND COMPUTING↗

Graphical and Numerical Description of Strain-Gage Balance Interactions

A new approach for the graphical and numerical description of balance interactions is presented. The approach uses data from single-component loads as input. This choice has the two advantages. First, the number of applied loads is at the minimum needed for interactions to be observed. In addition, the loads used for the description of interactions can easily be repeated at different sites. Output differences relative to the outputs of the zero load point of a load series are used for the description of interactions. Similarly, load differences relative to the loads of the zero load point of a load series are used for the description of the loads. Interactions are plotted versus the load differences for each load component while omitting the outputs of the primary gage of the chosen load component. The resulting plots have a characteristic star pattern as interactions are zero at zero load. Numerical estimates of the slopes of the interactions can be reverse-engineered from the load prediction equations of a balance if calibration data is examined. The slopes are the off-diagonal coefficients of the inverse of the matrix that has the coefficients of the linear terms of the fitted loads if the Non-Iterative Method is used for the analysis. Similarly, the slopes are the off-diagonal coefficients of the inverse of the matrix that is the non-iterative part of the primary load iteration equation if the Iterative Method is used for the analysis. Data from the calibration of a force balance is processed to illustrate the proposed graphical and numerical description of interactions.

wind tunnel test↗

Graphical and Numerical Description of Strain-Gage Balance Interactions

A new approach for the graphical and numerical description of strain-gage balance interactions is presented. The approach uses data from single-component loads as input. This choice has the two advantages. First, the number of applied loads is at the minimum needed for interactions to be observed. In addition, the loads used for the description of interactions can easily be repeated at different sites. Output differences relative to the outputs of the zero load point of a load series are used for the description of interactions. Similarly, load differences relative to the loads of the zero load point of a load series are used for the description of the loads. Interactions are plotted versus the load differences for each load component while omitting outputs of the primary gage of the chosen load component. The resulting plots have a star pattern as all interactions are zero at zero load. Estimates of the slopes of the interactions can be obtained from the load prediction equations of a balance if calibration data is examined. The slopes are the off-diagonal coefficients of the inverse of the matrix that has the coefficients of the linear terms of the fitted loads if the Non-Iterative Method is used for the analysis. Similarly, the slopes are the off-diagonal coefficients of the inverse of the matrix that is the non-iterative part of the primary load iteration equation if the Iterative Method is used for the analysis. Data sets from a manual and a machine calibration of a force balance are processed to illustrate the proposed description of interactions.

strain-gage balance↗

Elevated temperature mechanical properties of Inconel 617 surface oxide using nanoindentation

Inconel 617 is a principal candidate material for helium gas cooled very-high-temperature reactors with outlet temperatures of 700–950 °C. Recent findings showed that this alloy develops unique surface oxide layers especially at high temperature (HT) helium environment with distinctive wear, friction and contact properties. This study investigates the elevated temperature mechanical properties of Inconel 617 top surface layers aged in HT helium environment. Nanoindentation technique is used to obtain load-displacement graphs of the alloy top surface oxide in temperatures ranging from 25 °C up to 600 °C. In addition, using finite element analysis along with an iterative regression technique, a semi-numerical method is developed to further measure and quantify the material parameters and, in particular, time-independent and creep characteristics of the oxide. Although Young's modulus of the oxide is found to be relatively close to the bulk for the tested temperatures, the yield strength and hardness, in comparison to the bulk material, increase significantly as the material is oxidized after aging. The oxide exhibits significant softening as the temperature increases to 600 °C. Unlike the bulk material, diffusion through the grains is found to be the dominant creep mechanism for the oxide. Considerable difference between the mechanical properties of the oxide and the bulk material shows the need for accurate measurements of near surface mechanical properties, if reliable predictive contact and tribological models are sought at elevated temperatures.

36 MATERIALS SCIENCE↗

Experimental study of ECH pre-ionization on J-TEXT

An experimental study on electron cyclotron heating (ECH) pre-ionization has been conducted on J-TEXT in support of the joint experiment research for ITER plasma initiation. In this experiment, ECH power was injected to the vessel before the application of loop voltage ionize the neutral gas and form the initial plasma or so-called pre-plasma. The impact of several significant factors, such as magnetic field configuration, pre-fill gas pressure, ECH toroidal injection angle and ECH power on the evolution of pre-plasma are systematically studied, aiming to identify shared features, clarify their potential relationship and optimize the discharge parameters to generate a rather high pre-plasma density. By separating ECH power from the inductive start-up, the effect of pre-plasma on tokamak start-up can be observed. A dynamic magnetic configuration facilitates the transition of pre-plasma to tokamak plasma. To assess the influence of pre-plasma density on tokamak start-up, two kinds of magnetic field configurations are examined. While the effect of different pre-plasma densities on tokamak start-up is negligible, a significant difference is observed between pure ohmic start-up and start-up with pre-ionization. The studies presented here show evolutionary trends and threshold values needed to optimize ECH pre-ionization and also a feasible way to improve pre-plasma density and a configuration to stabilize pre-plasma for transition. Eventually, these results may contribute to multi-machine research and physics analysis that can assist ITER with its optimal preparation for first plasma operation.

ECH↗

Development of a Fast Spectrum Test Loop at the Advanced Test Reactor

The United States currently has no infrastructure to provide a fast neutron flux for research purposes. With the increased interest in Sodium Fast Reactor (SFR) technology, this type of infrastructure is necessary for further development of key components such as Versatile Test Reactor. The purpose of this research is to characterize the thermal hydraulics of a test loop capable of providing SFR like operating conditions. This includes showcasing feasibility in mechanical design, reaching a high fast to thermal neutron flux ratio ideally near 15, an inlet temperature of 600 K, a temperature increase of 50 K from inlet to outlet, and test chamber radial temperature profiles similar to SFR fuel rods. Existing infrastructure must be used to provide the neutron flux for the test loop, and it was decided that the Advanced Test Reactor (ATR) was the optimal choice due to its availability and previous work conducted similar to this project. To achieve the goals of this project, several mechanical designs were put forward with the most feasible design being chosen. The neutronics and thermal hydraulics analysis were conducted in tandem, with an iterative process used to maximize efficiency of both systems. Analysis of the neutronics system was done using Attila and MCNP while the thermal hydraulic analysis was done using RELAP5. This document outlines the thermal hydraulic analysis, though it includes the results from the mechanical design and neutronics analysis. The loop design includes a primary system with molten sodium, a set of “booster fuel plates” to increase the fast flux within the test section, a water system to cool the booster fuel plates, a helium system to separate the sodium and water systems, and a thermal neutron shield to decrease the thermal flux. Each of these systems is modeled in the neutronics and thermal hydraulic analyses. The RELAP5 model showcased the feasibility of the design proposed by this report. For all designs tested, the desired 50 K temperature increase from inlet to outlet was reached. The fuel rod radial temperature profile was lower than was expected, likely due to a lower rod thermal power than expected. A parametric study was also conducted to characterize the system with three fuel pins. This study showed the effect of thermal power, fluid pressure, fluid velocity, and inlet temperature on the system. ©Copyright by Anton Higgins December 11, 2019 All Rights Reserved

99 GENERAL AND MISCELLANEOUS↗

Computational methods of robust controller design for aerodynamic flutter suppression

The development of Riccati iteration, a tool for the design and analysis of linear control systems is examined. First, Riccati iteration is applied to the problem of pole placement and order reduction in two-time scale control systems. Order reduction, yielding a good approximation to the original system, is demonstrated using a 16th order linear model of a turbofan engine. Next, a numerical method for solving the Riccati equation is presented and demonstrated for a set of eighth order random examples. A literature review of robust controller design methods follows which includes a number of methods for reducing the trajectory and performance index sensitivity in linear regulators. Lastly, robust controller design for large parameter variations is discussed.

Anderson, L. R.↗

ESTIMATING THE DEEP SOLAR MERIDIONAL CIRCULATION USING MAGNETIC OBSERVATIONS AND A DYNAMO MODEL: A VARIATIONAL APPROACH

We show how magnetic observations of the Sun can be used in conjunction with an axisymmetric flux-transport solar dynamo model in order to estimate the large-scale meridional circulation throughout the convection zone. Our innovative approach rests on variational data assimilation, whereby the distance between predictions and observations (measured by an objective function) is iteratively minimized by means of an optimization algorithm seeking the meridional flow that best accounts for the data. The minimization is performed using a quasi-Newton technique, which requires knowledge of the sensitivity of the objective function to the meridional flow. That sensitivity is efficiently computed via the integration of the adjoint flux-transport dynamo model. Closed-loop (also known as twin) experiments using synthetic data demonstrate the validity and accuracy of this technique for a variety of meridional flow configurations, ranging from unicellular and equatorially symmetric to multicellular and equatorially asymmetric. In this well-controlled synthetic context, we perform a systematic study of the behavior of our variational approach under different observational configurations by varying their spatial density, temporal density, and noise level, as well as the width of the assimilation window. We find that the method is remarkably robust, leading in most cases to a recovery of the true meridional flow to within better than 1%. These encouraging results are a first step toward using this technique to (i) better constrain the physical processes occurring inside the Sun and (ii) better predict solar activity on decadal timescales.

79 ASTRONOMY AND ASTROPHYSICS↗

VERA and VERA-EDU 3.5 Release Notes

The Virtual Environment for Reactor Applications components included in this distribution include selected computational tools and supporting infrastructure that solve neutronics, thermal-hydraulics, fuel performance, and coupled neutronics-thermal hydraulics problems. The infrastructure components provide a simplified common user input capability and provide for the physics integration with data transfer and coupled-physics iterative solution algorithms. Neutronics analysis can be performed for 2D lattices, 2D core and 3D core problems for pressurized water reactor geometries that can be used to calculate criticality and fission rate distributions by pin for input fuel compositions. MPACT uses the Method of Characteristics transport approach for 2D problems. For 3D problems, MPACT uses the 2D/1D method which uses 2D MOC in a radial plane and diffusion or SPn in the axial direction. MPACT includes integrated cross section capabilities that provide problem-specific cross sections generated using the subgroup methodology. The code can be executed both 2D and 3D problems in parallel to reduce overall run time. A thermal-hydraulics capability is provided with CTF (an updated version of COBRA-TF) that allows thermal-hydraulics analyses for single and multiple assemblies using the simplified VERA common input. This distribution also includes coupled neutronics/thermal-hydraulics capabilities to allow calculations using MPACT coupled with CTF. The VERA fuel rod performance component BISON calculates, on a 2D or 3D basis, fuel rod temperature, fuel rod internal pressure, free gas volume, clad integrity and fuel rod waterside diameter. These capabilities allow simulation of power cycling, fuel conditioning and deconditioning, high burnup performance, power uprate scoping studies, and accident performance. Input/Output capabilities include the VERA Common Input (VERAIn) script which converts the ASCII common input file to the intermediate XML used to drive all of the physics codes in the VERA Core Simulator (VERA-CS). VERA component codes either input the VERA XML format directly, or provide a preprocessor which can convert the XML into native input. VERAView is an interactive graphical interface for the visualization and engineering analyses of output data from VERA. The python-based software is easy to install and intuitive to use, and provides instantaneous 2D and 3D images, 1D plots, and alpha-numeric data from VERA multi-physics simulations. Testing within CASL has focused primarily on Westinghouse four-loop reactor geometries and conditions with example problems included in the distribution.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

RAMI Analyses of Heating Neutral Beam and Diagnostic Neutral Beam Systems for ITER

A RAMI (Reliability, Availability, Maintainability, Inspectability) analysis has been performed for the heating (and current drive) neutral beam (HNB) and diagnostic neutral beam (DNB) systems of the ITER device. The objective of these analyses is to implement RAMI engineering requirements for design and testing to prepare a reliability-centred plan for commissioning, operation, and maintenance of the system in the framework of technical risk control to support the overall ITER Project. These RAMI requirements will correspond to the RAMI targets for the ITER project and the compensating provisions to reach them as deduced from the necessary actions to decrease the risk level of the function failure modes. The RAMI analyses results have to match with the procurement plan of the systems.

43 PARTICLE ACCELERATORS↗

Steady state thermo-mechanics and material property definition framework for analyzing DCLL blanket in the fusion nuclear science facility

In this work, a thermo-mechanics model that relies on creating the material property definition framework (MPDF) and multiphysics coupling of the heat transfer and the solid mechanics modules is developed to determine the structural integrity of the recently designed dual cooled lead lithium (DCLL) inboard blanket (IB) for the Fusion Nuclear Science Facility under steady state loads. The MPDF is called to supply fusion relevant neutron irradiation and temperature induced changes in material properties during multiphysics finite element runs, and PbLi temperature profiles are used to approximate Magnetohydrodynamics effect and the nuclear volumetric heating on the PbLi. Neutron irradiation and temperature induced reduction of the yield and ultimate strengths of F82H steel at the first wall (FW) are quantified for one year. A blanket in an assembly with gaps between blanket sectors and another blanket in an assembly with no gaps between blanket sectors, both exposed to radiation damage that lasted for one year are analyzed. Analysis using the elastic ITER structural design criteria for in-vessel components (ITER SDC-IC) design rules and a linear isotropic-hardening-type elastoplastic material model are used where most appropriate. The IB blanket with gaps between blanket sectors will withstand the steady state combined thermal and coolant loads for one year operational period but will fail if no gaps are allowed between blanket sectors. It is recommended that a gap of about 7.62 mm should be provided between IB blanket sectors during assembly which would close up during service, stop neutron streaming, reduce stresses and reduce bending of the FW into the scrape-off layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Determining Partial Atomic Charges for Liquid Water: Assessing Electronic Structure and Charge Models

Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, natural population analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree–Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 Å of the central water) in the quantum chemical calculation converges the charges of the central water molecule. Furthermore, our final pragmatic quantum chemical charge-assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 Å.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluating Water Assets Using Water Efficiency Framework for Infrastructure Intensive Agency

Water infrastructure-intensive organizations should incorporate sustainability principles when making decisions on infrastructure investments to ensure the long-term resilience of mission-critical equipment and components. A water efficiency evaluation framework to inform infrastructure investments is discussed and applied to three federal campuses. The methodology considered emphasizes the iterative process between data analysis and decision-making; evaluates cost-sharing solution alternatives when making final recommendations; considers both internal and external stakeholder perceptions of project outcomes and future planning; and makes explicit government and community preferences when evaluating projects. Case study findings underscore condition evaluations for infrastructure that are critical to setting a baseline for system performance and determining desired potential upgrades and sustainment strategies. Such a holistic approach has demonstrated effectiveness for understanding, classifying, and recommending innovative water asset solutions that are necessary for the public sector to consider in today’s complex, modern planning environment.

Asset Management, Alternative Water, evaluation, P↗

A high-speed, high-efficiency phase controller for coherent beam combining based on SPGD algorithm

A phase controller for coherent beam combining (CBC) of fibre lasers has been designed and manufactured based on a stochastic parallel gradient descent (SPGD) algorithm and a field programmable gate array (FPGA). The theoretical analysis shows that the iteration rate is higher than 1.9 MHz, and the average compensation bandwidth of CBC for 5 or 20 channels is 50 kHz or 12.5 kHz, respectively. The tests show that the phase controller ensures reliable phase locking of lasers: When the phases of five lasers are locked by the improved control strategy with a variable gain, the energy encircled in the target is increased by 23 times than that in the single output, the phase control accuracy is better than λ/20, and the combining efficiency is 92%. (control of laser radiation parameters)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Elucidating the G″ overshoot in soft materials with a yield transition via a time-resolved experimental strain decomposition

Significance The yielding transition in soft materials manifests macroscopically as a reversible transition from solid-like deformation to liquid-like flow under sufficiently large loading or deformation. Despite the ubiquity of materials that exhibit yielding, an understanding of the transition has remained elusive due to experimental limitations of the simple rheological tests typically used to study it. We present here a more thorough understanding of yielding, obtained via an iterative test protocol and analysis that allow for the separation of solid-like and liquid-like behaviors in the response of these materials. We show that yielding is a gradual transition, occurring to an increased extent as the deformation increases.

Donley, Gavin J.↗