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 127 records · Page 7

Stability analysis of nonlinear two-grid method for multigroup neutron diffusion problems

We present theoretical analysis of a nonlinear acceleration method for solving multigroup neutron diffusion problems. This method is formulated with two energy grids that are defined by (i) fine-energy groups structure and (ii) coarse grid with just a single energy group. The coarse-grid equations are derived by averaging of the multigroup diffusion equations over energy. The method uses a nonlinear prolongation operator. We perform stability analysis of iteration algorithms for inhomogeneous (fixed-source) and eigenvalue neutron diffusion problems. To apply Fourier analysis the equations of the method are linearized about solutions of infinite-medium problems. The developed analysis enables us to predict convergence properties of this two-grid method in different types of problems. Numerical results of problems in 2D Cartesian geometry are presented to confirm theoretical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evaluating the reforested area for the municipality of Buri by automatic analysis of LANDSAT imagery

The author has identified the following significant results. The class of reforestation (Pinus, Eucalyptus, Araucaria) was defined using iterative image analysis (1-100) and LANDSAT MSS data. Estimates of class area by 1-100 were compared with data supplied by the forestry institute in Sao Paulo. LANDSAT channels 4 and 5 served to differentiate the Pinus, Eucalyptus, and Araucaria from the other trees. Channels 6 and 7 gave best results for differentiating between the classes. A good representative spectral response was obtained for Auraucaria on these two channels. The small relative differences obtained were +4.24% for Araucaria, -7.51% for Pinus, and -32.07% for Eucalyptus.

Dejesusparada, N.↗

Combining Static Analysis and Model Checking for Software Analysis

We present an iterative technique in which model checking and static analysis are combined to verify large software systems. The role of the static analysis is to compute partial order information which the model checker uses to reduce the state space. During exploration, the model checker also computes aliasing information that it gives to the static analyzer which can then refine its analysis. The result of this refined analysis is then fed back to the model checker which updates its partial order reduction. At each step of this iterative process, the static analysis computes optimistic information which results in an unsafe reduction of the state space. However we show that the process converges to a fired point at which time the partial order information is safe and the whole state space is explored.

Brat, Guillaume↗

Comparing direct and iterative equation solvers in a large structural analysis software system

Two direct Choleski equation solvers and two iterative preconditioned conjugate gradient (PCG) equation solvers used in a large structural analysis software system are described. The two direct solvers are implementations of the Choleski method for variable-band matrix storage and sparse matrix storage. The two iterative PCG solvers include the Jacobi conjugate gradient method and an incomplete Choleski conjugate gradient method. The performance of the direct and iterative solvers is compared by solving several representative structural analysis problems. Some key factors affecting the performance of the iterative solvers relative to the direct solvers are identified.

Poole, E. L.↗

A modal aeroelastic analysis scheme for turbomachinery blading

An aeroelastic analysis is developed which has general application to all types of axial-flow turbomachinery blades. The approach is based on linear modal analysis, where the blade's dynamic response is represented as a linear combination of contributions from each of its in-vacuum free vibrational modes. A compressible linearized unsteady potential theory is used to model the flow over the oscillating blades. The two-dimensional unsteady flow is evaluated along several stacked axisymmetric strips along the span of the airfoil. The unsteady pressures at the blade surface are integrated to result in the generalized force acting on the blade due to simple harmonic motions. The unsteady aerodynamic forces are coupled to the blade normal modes in the frequency domain using modal analysis. An iterative eigenvalue problem is solved to determine the stability of the blade when the unsteady aerodynamic forces are included in the analysis. The approach is demonstrated by applying it to a high-energy subsonic turbine blade from a rocket engine turbopump power turbine. The results indicate that this turbine could undergo flutter in an edgewise mode of vibration.

Smith, Todd E.↗

Digital Safety Analysis for Small Modular Nuclear Reactors (SMRs)

A Documented Safety Analysis (DSA) is a Department of Energy (DOE) construct that defines the extent to which a nuclear facility can be operated safely. It includes a description of hazards, safe boundaries, and hazard controls. The authors assert that a Digital Safety Analysis (DgSA) is far superior to a legacy DSA for several reasons: • The underling database is structured such that it is possible to perform a comprehensive design review and safety analysis by iterating systematically across a hierarchy of linked objects versus a redundant and spotty review by entities of various abilities under unknown resource and schedule constraints. • The analysis of a new design can discover elements that are similar to elements in previous designs. The discovery of similarities is made possible by using the same structure for the underlying database for each new DgSA. The “prior learning” from previous designs is then applied automatically to new designs. • Outputs from the DgSA are from a single source to ensure consistency among various views of the same information. After the DgSA is released, the continued use of a single source implements a configuration management program to ensure consistency between the design basis, the design, the built system, and system procedures. • The development of the DgSA is agile in that any change in a linked object triggers an analysis of impacts on other linked objects and updates of linked objects are made accordingly. After the DgSA is released, the continued maintenance of these links and objects automates the “unreviewed safety question” process.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Blade design and analysis using a modified Euler solver

An iterative method for blade design based on Euler solver and described in an earlier paper is used to design compressor and turbine blades providing shock free transonic flows. The method shows a rapid convergence, and indicates how much the flow is sensitive to small modifications of the blade geometry, that the classical iterative use of analysis methods might not be able to define. The relationship between the required Mach number distribution and the resulting geometry is discussed. Examples show how geometrical constraints imposed upon the blade shape can be respected by using free geometrical parameters or by relaxing the required Mach number distribution. The same code is used both for the design of the required geometry and for the off-design calculations. Examples illustrate the difficulty of designing blade shapes with optimal performance also outside of the design point.

Leonard, O.↗

Fourier Analysis of the Multilevel QD Method with {alpha}-Approximation For Time-Dependent Radiative Transfer Problems

Thermal radiation transfer (TRT) problems describe interaction of photon radiation with matter. They are defined by the time-dependent radiative transfer (RT) equation for the specific intensity I coupled with the energy balance (EB) equation. This class of problems is characterised by high dimensionality, multiple scales and strong nonlinearity. In this paper we analyze the multilevel quasi-diffusion (QD) method for solving TRT problems. This is a nonlinear method of moments with exact closures. The QD method enables one to reduce dimensionality of the original radiative transfer problem. It is formulated by a nonlinear system of equations consisting of (i) the high-order RT (HORT) equation for the intensity of photon radiation and (ii) low-order equations for the angular and energy moments of the intensity. The nonlinear closures are formulated by means of stable factors. The low-order problem brings the radiative transfer model to the level of multi-physics equations. The multilevel QD (MLQD) method possesses attractive features for efficient multi-physics coupling. To calculate the high-order RT solution we apply the {alpha}-approximation which assumes that the intensity changes exponentially in time over a time step. The {alpha}-approximation leads to the steady-state RT equation with a modified opacity and enables one to avoid storing the high-dimensional solution from the previous time step. The rate of change is evaluated by the solution of the time-dependent low-order QD (LOQD) equations for moments. We note that the HORT equation in the {alpha}-approximation contains an extra nonlinear term. It affects the iterative behaviour of the MLQD method. The obtained computational results indicated that the rate of convergence is slower compared to the rate when the time-dependent HORT is used to calculate the intensity. In this paper we perform an iterative stability analysis of the MLQD method in a discretized form and analyze convergence of outer (transport) iterations on a time step. We consider multigroup TRT problems in 1D slab geometry. The fully implicit scheme is used to discretize equations in time. The HORT equation is approximated with the method of step characteristics. A finite volume method is applied for spatial discretization of multigroup LOQD equations. The grey (one-group) LOQD equations are algebraically consistent with the multigroup LOQD equations. To perform Fourier analysis a model TRT problem in an infinite homogeneous medium is considered. The system of nonlinear equations of the MLQD method is linearized near an equilibrium solution at some given temperature. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An approach to probabilistic finite element analysis using a mixed-iterative formulation

An efficient algorithm for computing the response sensitivity of finite element problems based on a mixed-iterative formulation is proposed. This method does not involve explicit differentiation of the tangent stiffness array and can be used with formulations for which a consistent tangent stiffness is not readily available. The method has been successfully applied to probabilistic finite element analysis of problems using the proposed mixed formulation, and this exercise has provided valuable insights regarding the extension of the method to a more general class of problems to include material and geometric nonlinearities.

Dias, J. B.↗

Neural network acceleration of large-scale structure theory calculations

Here, we make use of neural networks to accelerate the calculation of power spectra required for the analysis of galaxy clustering and weak gravitational lensing data. For modern perturbation theory codes, evaluation time for a single cosmology and redshift can take on the order of two seconds. In combination with the comparable time required to compute linear predictions using a Boltzmann solver, these calculations are the bottleneck for many contemporary large-scale structure analyses. Here, we construct neural network-based surrogate models for Lagrangian perturbation theory (LPT) predictions of matter power spectra, real and redshift space galaxy power spectra, and galaxy-matter cross power spectra that attain ~ 0.1% (at one sigma) accuracy over a broad range of scales in a ωCDM parameter space. The neural network surrogates can be evaluated in approximately one millisecond, a factor of 1000 times faster than the full Boltzmann code and LPT computations. In a simulated full-shape redshift space galaxy power spectrum analysis, we demonstrate that the posteriors obtained using our surrogates are accurate compared to those obtained using the full LPT model. We make our surrogate models public at https://github.com/sfschen/EmulateLSS, so that others may take advantage of the speed gains they provide to enable rapid iteration on analysis settings, something that is essential in complex contemporary large-scale structure analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

From Informal Safety-Critical Requirements to Property-Driven Formal Validation

Most of the efforts in formal methods have historically been devoted to comparing a design against a set of requirements. The validation of the requirements themselves, however, has often been disregarded, and it can be considered a largely open problem, which poses several challenges. The first challenge is given by the fact that requirements are often written in natural language, and may thus contain a high degree of ambiguity. Despite the progresses in Natural Language Processing techniques, the task of understanding a set of requirements cannot be automatized, and must be carried out by domain experts, who are typically not familiar with formal languages. Furthermore, in order to retain a direct connection with the informal requirements, the formalization cannot follow standard model-based approaches. The second challenge lies in the formal validation of requirements. On one hand, it is not even clear which are the correctness criteria or the high-level properties that the requirements must fulfill. On the other hand, the expressivity of the language used in the formalization may go beyond the theoretical and/or practical capacity of state-of-the-art formal verification. In order to solve these issues, we propose a new methodology that comprises of a chain of steps, each supported by a specific tool. The main steps are the following. First, the informal requirements are split into basic fragments, which are classified into categories, and dependency and generalization relationships among them are identified. Second, the fragments are modeled using a visual language such as UML. The UML diagrams are both syntactically restricted (in order to guarantee a formal semantics), and enriched with a highly controlled natural language (to allow for modeling static and temporal constraints). Third, an automatic formal analysis phase iterates over the modeled requirements, by combining several, complementary techniques: checking consistency; verifying whether the requirements entail some desirable properties; verify whether the requirements are consistent with selected scenarios; diagnosing inconsistencies by identifying inconsistent cores; identifying vacuous requirements; constructing multiple explanations by enabling the fault-tree analysis related to particular fault models; verifying whether the specification is realizable.

Cimatti, Alessandro↗

Iterative methods for large scale static analysis of structures on a scalable multiprocessor supercomputer

A parallel Preconditioned Conjugate Gradient (PCG) iterative solver has been developed and implemented on the iPSC-860 scalable hypercube. This new implementation makes use of the Parallel Automated Runtime Toolkit at ICASE (PARTI) primitives to efficiently program irregular communications patterns that exist in general sparse matrices and in particular in the finite element sparse stiffness matrices. The iterative PCG has been used to solve the finite element equations that result from discretizing large scale aerospace structures. In particular, the static response of the High Speed Civil Transport (HSCT) finite element model is solved on the iPSC-860.

Sobh, Nahil Atef↗

Numerical Analysis of Fatigue Behavior of ITER-Like Monoblock Divertor Interlayer Under Coupled Heat Loads

High-confinement mode is a very prominent operation style for future fusion device due to its unique advantages. However, the conjuncted edge localized modes (ELMs) are very difficult to control so that divertor plates are very prone to suffer both stationary high heat flux (HHF) loads of long-pulse operating mode and transient shock loads of ELMs. Most previous researches focus on degradation of plasma facing material (PFM), however, as a layer joining PFM and cooling tube, the soft copper interlayer suffers concentrated thermal stress loads due to mismatched thermal expansion of PFM and cooling tube. Its thermal fatigue behavior under such coupled loads is also of great significance to structural safety of divertor component. With such a motivation, the reduction effects on fatigue life time of a typical interlayer of monoblock divertor under series of coupled HHF and ELMs shock loading conditions are investigated. It is found that: (1) The transient shock feature of ELMs loading is propagated into interlayer with less sharp pattern. The increase of damage induced by coupled ELMs loading is limited in single cycle, while the accumulated damage of multiple consecutive coupled loading cycles is increased nonlinearly. (2) Under the coupled HHF and ELMs loading, the fatigue life time of interlayer is generally decreasing. The magnitude of decrease is increasing nonlinearly with the magnitude of ELMs peak and averaged heat flux. (3) For three characteristic parameters of ELMs shock loading such as frequency, duration and peak heat flux, the peak heat flux and frequency are two parameters more sensitive to determine coupled reduction effects on fatigue lifetime of the interlayer, while for high frequency case, time averaged heat flux takes the lead.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

GPU code optimization using abstract kernel emulation and sensitivity analysis

In this paper, we develop a novel approach to GPU kernel optimization by focusing on identification of bottleneck resources and determining optimization parameters that can alleviate the bottleneck. Performance modeling for GPUs is done by abstract kernel emulation along with latency/gap modeling of resources. Sensitivity analysis with respect to resource latency/gap parameters is used to predict the bottleneck resource for a given kernel's execution. Here, the utility of the bottleneck analysis is demonstrated in two contexts: 1) Enhancing the OpenTuner auto-tuner with the new bottleneck-driven optimization strategy. Experimental results on all kernels from the Rodinia suite and GPU tensor contraction kernels from the NWChem computational chemistry suite demonstrate effectiveness. 2) Manual code optimization: two case studies illustrate the use of the bottleneck analysis to iteratively improve performance of code from state-of-the-art DSL code generators.

97 MATHEMATICS AND COMPUTING↗

Analysis and Evaluation of NASA's MK40B Six-Component Task Balance

The MK40B balance is NASA's newest balance of Task design. Its load range was optimized for use in the NASA Ames 11-ft Transonic Wind Tunnel. The balance was calibrated in 2019 in Calspan's Automatic Balance Calibration System. The data analysis was performed at the Ames Balance Calibration Laboratory. First, bi-directional characteristics of the six gage outputs of the balance were assessed by applying a semi-empirical test to the calibration data. These tests indicated that four of the six gages have bi-directional characteristics. Therefore, the use of absolute value terms in the regression models of the calibration data is justified. Then, the machine calibration data was analyzed using both the Non-Iterative and the Iterative Method. Analysis results for the two methods were compared using the percent contributions, the calibration load residuals, the primary sensitivities, and the principle linear terms. These results confirmed that the two methods describe bi-directional characteristics of the balance gages correctly. In addition, it was concluded that the load prediction accuracies of the two methods are the same for all practical purposes.

wind tunnel balance↗

Analysis of Landsat-4 Thematic Mapper data for classification of forest stands in Baldwin County, Alabama

A computer-implemented classification has been derived from Landsat-4 Thematic Mapper data acquired over Baldwin County, Alabama on January 15, 1983. One set of spectral signatures was developed from the data by utilizing a 3x3 pixel sliding window approach. An analysis of the classification produced from this technique identified forested areas. Additional information regarding only the forested areas. Additional information regarding only the forested areas was extracted by employing a pixel-by-pixel signature development program which derived spectral statistics only for pixels within the forested land covers. The spectral statistics from both approaches were integrated and the data classified. This classification was evaluated by comparing the spectral classes produced from the data against corresponding ground verification polygons. This iterative data analysis technique resulted in an overall classification accuracy of 88.4 percent correct for slash pine, young pine, loblolly pine, natural pine, and mixed hardwood-pine. An accuracy assessment matrix has been produced for the classification.

Hill, C. L.↗

A Modified Wall Matching Treatment to Account for Local Solid to Fluid Thermal Coupling

The wall-matching methodology of Wilcox is modified to include a solid-wall, thermal-conduction model. This coupled fluid-thermal-structure model is derived assuming that the wall thermal-structure behavior is locally one-dimensional and that structural deformations, due to thermally induced stresses, are not significant. The one-dimensional coupled fluid-thermal-structure model is derived such that the wall temperature is removed as an independent boundary condition variable. The one-dimensional coupled fluid-thermal-structure model is also derived for the general case of an arbitrary mixture of thermally prefect gases and a wall of arbitrary thickness and conductivity by using a compressible, streamwise-pressure-gradient-corrected, wall-matching function and Fourier's law of heat conduction. The resulting model was implemented in the VULCAN CFD code as a new boundary condition type. VULCAN was then used to simulate a two-dimensional Mach 6 wind tunnel facility nozzle flow to demonstrate/validate the one-dimensional coupled fluid-thermal-structure model. The nozzle internal-wall surface temperature and heat transfer distributions computed using the one-dimensional coupled fluid-thermal-structure model are compared to wall temperature and heat transfer distributions from an iterative multi-dimensional analysis obtained by coupling the VULCAN CFD code and the MSC/NASTRAN-thermal code. The one-dimensional coupled fluid-thermal-structure model analysis is shown to be very robust and in excellent agreement with the multi-dimensional iteratively coupled analysis. It is also shown that the one-dimensional analysis can be used as an initial guess for the multi-dimensional iteratively coupled analysis.

White, Jeffery A.↗

A Modified Wall Matching Treatment to Account for Local Solid to Fluid Thermal Coupling

The wall-matching methodology of Wilcox is modified to include a solid-wall, thermal-conduction model. This coupled fluid-thermal-structure model is derived assuming that the wall thermal-structure behavior is locally one-dimensional and that structural deformations, due to thermally induced stresses, are not significant. The one-dimensional coupled fluid-thermal-structure model is derived such that the wall temperature is removed as an independent boundary condition variable. The one-dimensional coupled fluid-thermal-structure model is also derived for the general case of an arbitrary mixture of thermally prefect gases and a wall of arbitrary thickness and conductivity by using a compressible, streamwise-pressure-gradient-corrected, wall-matching function and Fourier's law of heat conduction. The resulting model was implemented in the VULCAN CFD code as a new boundary condition type. VULCAN was then used to simulate a two-dimensional Mach 6 wind tunnel facility nozzle flow to demonstrate/validate the one-dimensional coupled fluid-thermal-structure model. The nozzle internal-wall surface temperature and heat transfer distributions computed using the one-dimensional coupled fluid-thermal-structure model are compared to wall temperature and heat transfer distributions from an iterative multi-dimensional analysis obtained by coupling the VULCAN CFD code and the MSC/NASTRAN-thermal code. The one-dimensional coupled fluid-thermal-structure model analysis is shown to be very robust and in excellent agreement with the multi-dimensional iteratively coupled analysis. It is also shown that the one-dimensional analysis can be used as an initial guess for the multi-dimensional iteratively coupled analysis.

White, Jeffery A.↗