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 667 records · Page 37

Shape reanalysis and sensitivities utilizing preconditioned iterative boundary solvers

The computational advantages associated with the utilization of preconditined iterative equation solvers are quantified for the reanalysis of perturbed shapes using continuum structural boundary element analysis (BEA). Both single- and multi-zone three-dimensional problems are examined. Significant reductions in computer time are obtained by making use of previously computed solution vectors and preconditioners in subsequent analyses. The effectiveness of this technique is demonstrated for the computation of shape response sensitivities required in shape optimization. Computer times and accuracies achieved using the preconditioned iterative solvers are compared with those obtained via direct solvers and implicit differentiation of the boundary integral equations. It is concluded that this approach employing preconditioned iterative equation solvers in reanalysis and sensitivity analysis can be competitive with if not superior to those involving direct solvers.

Guru Prasad, K.↗

3D neutronic/thermal-hydraulic coupled analysis of MYRRHA

The current tendency in multiphysics calculations applied to reactor physics is the use of already validated computer codes, coupled by means of an iterative approach. In this paper such an approach is explained concerning neutronics and thermal-hydraulics coupled analysis with MCNPX and COBRA-IV codes using a driver program and file exchange between codes. MCNPX provides the neutronic analysis of heterogeneous nuclear systems, both in critical and subcritical states, while COBRA-IV is a subchannel code that can be used for rod bundles or core thermal-hydraulics analysis. In our model, the MCNP temperature dependence of nuclear data is handled via pseudo-material approach, mixing pre-generated cross section data set to obtain the material with the desired cross section temperature. On the other hand, COBRA-IV has been updated to allow for the simulation of liquid metal cooled reactors. The coupled computational tool can be applied to any geometry and coolant, as it is the case of single fuel assembly, at pin-by-pin level, or full core simulation with the average pin of each fuel-assembly. The coupling tool has been applied to the critical core layout of the SCK-CEN MYRRHA concept, an experimental LBE cooled fast reactor presently in engineering design stage. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Three-dimensional Finite Element Formulation and Scalable Domain Decomposition for High Fidelity Rotor Dynamic Analysis

This paper has two objectives. The first objective is to formulate a 3-dimensional Finite Element Model for the dynamic analysis of helicopter rotor blades. The second objective is to implement and analyze a dual-primal iterative substructuring based Krylov solver, that is parallel and scalable, for the solution of the 3-D FEM analysis. The numerical and parallel scalability of the solver is studied using two prototype problems - one for ideal hover (symmetric) and one for a transient forward flight (non-symmetric) - both carried out on up to 48 processors. In both hover and forward flight conditions, a perfect linear speed-up is observed, for a given problem size, up to the point of substructure optimality. Substructure optimality and the linear parallel speed-up range are both shown to depend on the problem size as well as on the selection of the coarse problem. With a larger problem size, linear speed-up is restored up to the new substructure optimality. The solver also scales with problem size - even though this conclusion is premature given the small prototype grids considered in this study.

Datta, Anubhav↗

Low heat transfer oxidizer heat exchanger design and analysis

The RL10-IIB engine, a derivative of the RLIO, is capable of multi-mode thrust operation. This engine operates at two low thrust levels: tank head idle (THI), which is approximately 1 to 2 percent of full thrust, and pumped idle (PI), which is 10 percent of full thrust. Operation at THI provides vehicle propellant settling thrust and efficient engine thermal conditioning; PI operation provides vehicle tank pre-pressurization and maneuver thrust for log-g deployment. Stable combustion of the RL10-IIB engine at THI and PI thrust levels can be accomplished by providing gaseous oxygen at the propellant injector. Using gaseous hydrogen from the thrust chamber jacket as an energy source, a heat exchanger can be used to vaporize liquid oxygen without creating flow instability. This report summarizes the design and analysis of a United Aircraft Products (UAP) low-rate heat transfer heat exchanger concept for the RL10-IIB rocket engine. The design represents a second iteration of the RL10-IIB heat exchanger investigation program. The design and analysis of the first heat exchanger effort is presented in more detail in NASA CR-174857. Testing of the previous design is detailed in NASA CR-179487.

Kanic, P. G.↗

Computational Investigation of a Boundary-Layer Ingestion Propulsion System for the Common Research Model

This thesis will examine potential propulsive and aerodynamic benefits of integrating a boundary-layer ingestion (BLI) propulsion system with a typical commercial aircraft using the Common Research Model geometry and the NASA Tetrahedral Unstructured Software System (TetrUSS). The Numerical Propulsion System Simulation (NPSS) environment will be used to generate engine conditions for CFD analysis. Improvements to the BLI geometry will be made using the Constrained Direct Iterative Surface Curvature (CDISC) design method. Previous studies have shown reductions of up to 25% in terms of propulsive power required for cruise for other axisymmetric geometries using the BLI concept. An analysis of engine power requirements, drag, and lift coefficients using the baseline and BLI geometries coupled with the NPSS model are shown. Potential benefits of the BLI system relating to cruise propulsive power are quantified using a power balance method and a comparison to the baseline case is made. Iterations of the BLI geometric design are shown and any improvements between subsequent BLI designs presented. Simulations are conducted for a cruise flight condition of Mach 0.85 at an altitude of 38,500 feet and an angle of attack of 2deg for all geometries. A comparison between available wind tunnel data, previous computational results, and the original CRM model is presented for model verification purposes along with full results for BLI power savings. Results indicate a 14.3% reduction in engine power requirements at cruise for the BLI configuration over the baseline geometry. Minor shaping of the aft portion of the fuselage using CDISC has been shown to increase the benefit from boundary-layer ingestion further, resulting in a 15.6% reduction in power requirements for cruise as well as a drag reduction of eighteen counts over the baseline geometry.

Blumenthal, Brennan↗

Update to the Phase II Corrective Action Investigation Plan for Corrective Action Units 101 and 102, Central and Western Pahute Mesa, Nevada National Security Site, Nye County, Nevada

This report is an update to the Phase II Corrective Action Investigation Plan for Corrective Action Units 101 and 102: Central and Western Pahute Mesa, Nevada Test Site, Nye County, Nevada (NNSA/NSO, 2009) activities resulting from data collection and analysis conducted between 2009 and 2019. This report contains two sections: (1) a summary of how the new data have informed the updated approach to complete the characterization activities under the corrective action investigation (CAI), and (2) a table documenting the specific activities identified in the CAIP with a status and expected outcome of each activity (see Table A-1 in Appendix A). As stated in Sections 1.5 and 1.5.1 of the CAIP, the CAIP was conceptualized to be iterative, with continuous review and assessment of the data, data analysis, and modeling to provide a guide for decisions regarding additional work. This report provides the revised plan, based on an assessment of the data, to complete the characterization activities identified in the CAIP.

54 ENVIRONMENTAL SCIENCES↗

The solution of linear systems of equations with a structural analysis code on the NAS CRAY-2

Two methods for solving linear systems of equations on the NAS Cray-2 are described. One is a direct method; the other is an iterative method. Both methods exploit the architecture of the Cray-2, particularly the vectorization, and are aimed at structural analysis applications. To demonstrate and evaluate the methods, they were installed in a finite element structural analysis code denoted the Computational Structural Mechanics (CSM) Testbed. A description of the techniques used to integrate the two solvers into the Testbed is given. Storage schemes, memory requirements, operation counts, and reformatting procedures are discussed. Finally, results from the new methods are compared with results from the initial Testbed sparse Choleski equation solver for three structural analysis problems. The new direct solvers described achieve the highest computational rates of the methods compared. The new iterative methods are not able to achieve as high computation rates as the vectorized direct solvers but are best for well conditioned problems which require fewer iterations to converge to the solution.

Poole, Eugene L.↗

Nonlinear Structural Analysis Methodology and Dynamics Scaling of Inflatable Parabolic Reflector Antenna Concepts

Ultra-light weight and ultra-thin membrane inflatable antenna concepts are fast evolving to become the state-of-the-art antenna concepts for deep-space applications. NASA Langley Research Center has been involved in the structural dynamics research on antenna structures. One of the goals of the research is to develop structural analysis methodology for prediction of the static and dynamic response characteristics of the inflatable antenna concepts. This research is focused on the computational studies to use nonlinear large deformation finite element analysis to characterize the ultra-thin membrane responses of the antennas. Recently, structural analyses have been performed on a few parabolic reflector antennas of varying size and shape, which are referred in the paper as 0.3 meters subscale, 2 meters half-scale, and 4 meters full-scale antenna. The various aspects studied included nonlinear analysis methodology and solution techniques, ways to speed convergence in iterative methods, the sensitivities of responses with respect to structural loads, such as inflation pressure, gravity, and pretension loads in the ground and in-space conditions, and the ultra-thin membrane wrinkling characteristics. Several such intrinsic aspects studied have provided valuable insight into evaluation of structural characteristics of such antennas. While analyzing these structural characteristics, a quick study was also made to assess the applicability of dynamics scaling of the half-scale antenna. This paper presents the details of the nonlinear structural analysis results, and discusses the insight gained from the studies on the various intrinsic aspects of the analysis methodology. The predicted reflector surface characteristics of the three inflatable ultra-thin membrane parabolic reflector antenna concepts are presented as easily observable displacement fringe patterns with associated maximum values, and normal mode shapes and associated frequencies. Wrinkling patterns are presented to show how surface wrinkle progress with increasing tension loads. Antenna reflector surface accuracies were found to be very much dependent on the type and size of the antenna, the reflector surface curvature, reflector membrane supports in terms of spacing of catenaries, as well as the amount of applied load.

Sreekantamurthy, Tham↗

Software issues in three dimensional continuum shape optimization employing boundary formulations

The paper addresses the issue of how individual computational techniques for the accurate and economical calculation of information required for large-scale 3D continuum structural shape optimization via boundary element analysis (BEA) formulations can be prudently selected and incorporated in large-scale BEA programs, employing either direct or iterative equation solvers. A series of techniques that allows for the incorporation of economical shape design sensitivity analysis capability with a minimal amount of software development is described. The use of reanalysis is shown to be the key ingredient associated with each of these techniques. It is concluded that, from a software engineering perspective, capabilities that facilitate 3D shape optimization can be implemented in BEA programs with only modest investments in program development.

Kane, James H.↗

Local convergence analysis of an inexact trust-region method for nonsmooth optimization

In Baraldi, we introduced an inexact trust-region algorithm for minimizing the sum of a smooth nonconvex function and a nonsmooth convex function in Hilbert space—a class of problems that is ubiquitous in data science, learning, optimal control, and inverse problems. Furthermore, this algorithm has demonstrated excellent performance and scalability with problem size. In this paper, we enrich the convergence analysis for this algorithm, proving strong convergence of the iterates with guaranteed rates. In particular, we demonstrate that the trust-region algorithm recovers superlinear, even quadratic, convergence rates when using a second-order Taylor approximation of the smooth objective function term.

97 MATHEMATICS AND COMPUTING↗

A STATISTICAL RECONSTRUCTION OF THE PLANET POPULATION AROUND KEPLER SOLAR-TYPE STARS

Using the cumulative catalog of planets detected by the NASA Kepler mission, we reconstruct the intrinsic occurrence of Earth- to Neptune-size (1-4 R {sub ⊕}) planets and their distributions with radius and orbital period. We analyze 76,711 solar-type (0.8 < R {sub *}/R {sub ☉} < 1.2) stars with 430 planets on 20-200 day orbits, excluding close-in planets that may have been affected by the proximity to the host star. Our analysis considers errors in planet radii and includes an ''iterative simulation'' technique that does not bin the data. We find a radius distribution that peaks at 2-2.8 Earth radii, with lower numbers of smaller and larger planets. These planets are uniformly distributed with logarithmic period, and the mean number of such planets per star is 0.46 ± 0.03. The occurrence is ∼0.66 if planets interior to 20 days are included. We estimate the occurrence of Earth-size planets in the ''habitable zone'' (defined as 1-2 R {sub ⊕}, 0.99-1.7 AU for solar-twin stars) as 6.4{sub −1.1}{sup +3.4}%. Our results largely agree with those of Petigura et al., although we find a higher occurrence of 2.8-4 Earth-radii planets. The reasons for this excess are the inclusion of errors in planet radius, updated Huber et al. stellar parameters, and also the exclusion of planets that may have been affected by proximity to the host star.

79 ASTRONOMY AND ASTROPHYSICS↗

Real-time confinement regime detection in fusion plasmas with convolutional neural networks and high-bandwidth edge fluctuation measurements

Abstract A real-time detection of the plasma confinement regime can enable new advanced plasma control capabilities for both the access to and sustainment of enhanced confinement regimes in fusion devices. For example, a real-time indication of the confinement regime can facilitate transition to the high-performing wide-pedestal (WP) quiescent H-mode, or avoid unwanted transitions to lower confinement regimes that may induce plasma termination. To demonstrate real-time confinement regime detection, we use the 2D beam emission spectroscopy (BES) diagnostic system to capture localized density fluctuations of long wavelength turbulent modes in the edge region at a 1 MHz sampling rate. BES data from 330 discharges in either L-mode, H-mode, quiescent H (QH)-mode, or WP QH-mode were collected from the DIII-D tokamak and curated to develop a high-quality database to train a deep-learning classification model for real-time confinement detection. We utilize the 6×8 spatial configuration with a time window of 1024 µ s and recast the input to obtain spectral-like features via fast Fourier transform preprocessing. We employ a shallow 3D convolutional neural network for the multivariate time-series classification task and utilize a softmax in the final dense layer to retrieve a probability distribution over the different confinement regimes. Our model classifies the global confinement state on 44 unseen test discharges with an average F 1 score of 0.94, using only ∼1 ms snippets of BES data at a time. This activity demonstrates the feasibility for real-time data analysis of fluctuation diagnostics in future devices such as ITER, where the need for reliable and advanced plasma control is urgent.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NF-κB perturbation reveals unique immunomodulatory functions in Prx1 + fibroblasts that promote development of atopic dermatitis

Skin is composed of diverse cell populations that cooperatively maintain homeostasis. Up-regulation of the nuclear factor κB (NF-κB) pathway may lead to the development of chronic inflammatory disorders of the skin, but its role during the early events remains unclear. Here, through analysis of single-cell RNA sequencing data via iterative random forest leave one out prediction, an explainable artificial intelligence method, we identified an immunoregulatory role for a unique paired related homeobox-1 (Prx1) + fibroblast subpopulation. Disruption of Ikkb–NF-κB under homeostatic conditions in these fibroblasts paradoxically induced skin inflammation due to the overexpression of C-C motif chemokine ligand 11 (CCL11; or eotaxin-1) characterized by eosinophil infiltration and a subsequent T H 2 immune response. Because the inflammatory phenotype resembled that seen in human atopic dermatitis (AD), we examined human AD skin samples and found that human AD fibroblasts also overexpressed CCL11 and that perturbation of Ikkb–NF-κB in primary human dermal fibroblasts up-regulated CCL11. Monoclonal antibody treatment against CCL11 was effective in reducing the eosinophilia and T H 2 inflammation in a mouse model. Together, the murine model and human AD specimens point to dysregulated Prx1 + fibroblasts as a previously unrecognized etiologic factor that may contribute to the pathogenesis of AD and suggest that targeting CCL11 may be a way to treat AD-like skin lesions.

60 APPLIED LIFE SCIENCES↗

Predictive Maturity of Multi-Scale Simulation Models for Fuel Performance

The project proposed to provide a Predictive Maturity Framework with its companion metrics that (1) introduce a formalized, quantitative means to communicate information between interested parties, (2) provide scientifically dependable means to claim completion of Validation and Uncertainty Quantification (VU) activities, and (3) guide the decision makers in the allocation of Nuclear Energy’s resources for code development and physical experiments. The project team proposed to develop this framework based on two complimentary criteria: (1) the extent of experimental evidence available for the calibration of simulation models and (2) the sophistication of the physics incorporated in simulation models. The proposed framework is capable of quantifying the interaction between the required number of physical experiments and degree of physics sophistication. The project team has developed this framework and implemented it with a multi-scale model for simulating creep of a core reactor cladding. The multi-scale model is composed of the viscoplastic self-consistent (VPSC) code at the meso-scale, which represents the visco-plastic behavior and changing properties of a highly anisotropic material and a Finite Element (FE) code at the macro-scale to represent the elastic behavior and apply the loading. The framework developed takes advantage of the transparency provided by partitioned analysis, where independent constituent codes are coupled in an iterative manner. This transparency allows model developers to better understand and remedy the source of biases and uncertainties, whether they stem from the constituents or the coupling interface by exploiting separate-effect experiments conducted within the constituent domain and integral-effect experiments conducted within the full-system domain. The project team has implemented this procedure with the multi- scale VPSC-FE model and demonstrated its ability to improve the predictive capability of the model. Within this framework, the project team has focused on optimizing resource allocation for improving numerical models through further code development and experimentation. Related to further code development, we have developed a code prioritization index (CPI) for coupled numerical models. CPI is implemented to effectively improve the predictive capability of the coupled model by increasing the sophistication of constituent codes. In relation to designing new experiments, we investigated the information gained by the addition of each new experiment used for calibration and bias correction of a simulation model. Additionally, the variability of ‘information gain’ through the design domain has been investigated in order to identify the experiment settings where maximum information gain occurs and thus guide the experimenters in the selection of the experiment settings. This idea was extended to evaluate the information gain from each experiment can be improved by intelligently selecting the experiments, leading to the development of the Batch Sequential Design (BSD) technique. Additionally, we evaluated the importance of sufficiently exploring the domain of applicability in experiment-based validation of high-consequence modeling and simulation by developing a new metric to quantify coverage. This metric has also been incorporated into the design of new experiments. Finally, we have proposed a data-aware calibration approach for the calibration of numerical models. This new method considers the complexity of a numerical model (the number of parameters to be calibrated, parameter uncertainty, and form of the model) and seeks to identify the number of experiments necessary to calibrate the model based on the level of sophistication of the physics. The final component in the project team’s work to improve model calibration and validation methods is the incorporation of robustness to non-probabilistic uncertainty in the input parameters. This is an improvement to model validation and uncertainty quantification stemming beyond the originally proposed scope of the project. We have introduced a new metric for incorporating the concept of robustness into experiment-based validation of numerical models. This project has accounted for the graduation of two Ph.D. students (Kendra Van Buren and Josh Hegenderfer) and two M.S. students (Matthew Egeberg and Parker Shields). One of the doctoral students is now working in the nuclear engineering field and the other one is a post-doctoral fellow at the Los Alamos National Laboratory. Additionally, two more Ph.D. students (Garrison Stevens and Tunc Kulaksiz) who are working towards graduation have been supported by this project.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Calculations of molecular ionization energies using a self-consistent-charge Hartree-Fock-Slater method

A numerical-variational method for performing self-consistent molecular calculations in the Hartree-Fock-Slater (HFS) model is presented. Molecular wavefunctions are expanded in terms of basis sets constructed from numerical HFS solutions of selected one-center atomlike problems. Binding energies and wavefunctions for the molecules are generated using a discrete variational method for a given molecular potential. In the self-consistent-charge (SCC) approximation to the complete self-consistent-field (SCF) method, results of a Mulliken population analysis of the molecular eigenfunctions are used in each iteration to produce 'atomic' occupation numbers. The simplest SCC potential is then obtained from overlapping spherical atomlike charge distributions. Molecular ionization energies are calculated using the transition-state procedure; results are given for CO, H2O, H2S, AlCl, InCl, and the Ni5O surface complex. Agreement between experimental and theoretical ionization energies for the free-molecule valence levels is generally within 1 eV. The simple SCC procedure gives a reasonably good approximation to the molecular potential, as shown by comparison with experiment, and with complete SCF calculations for CO, H2O, and H2S.

Rosen, A.↗

Low cost solar array project: Experimental process system development unit for producing semiconductor-grade silicon using silane-to-silicon process

The design, fabrication, and installation of an experimental process system development unit (EPSDU) were analyzed. Supporting research and development were performed to provide an information data base usable for the EPSDU and for technological design and economical analysis for potential scale-up of the process. Iterative economic analyses were conducted for the estimated product cost for the production of semiconductor grade silicon in a facility capable of producing 1000-MT/Yr.

Source record↗

Simplified solution for stresses and deformation

Conventional contact deformation analysis for ball bearings, gears, and cams involves tedious iterative procedures or the use of design charts. A simplified approach that makes the elastic deformation at the center of contact easy to calculate was previously reported. The range of validity in which these equations can be used is extended. A simplified approach to the calculation of the location and magnitude of subsurface stresses developed in machine element applications is included.

Hamrock, B. J.↗

Temperature retrievals from TIROS-N

A nonstatistical, nonlinear iterative scheme, has been developed, based on the methods of Chahine (1970), to retrieve surface temperature, vertical temperature profile and cloud distribution from TIROS-N HIRS 2/MSU observations. Fundamental to the technique is the ability to accurately calculate clear column radiances given atmospheric and surface characteristics. The cloud filtering algorithm, making simultaneous use of IR 15 micron channels and microwave channels, allows accurate infrared retrievals in the presence of up to 80 percent cloud cover. Details for the iterative scheme are presented together with preliminary results from analysis of TIROS-N data.

Susskind, J.↗