Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical algorithm computing”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

A three-dimensional Navier-Stokes/Euler code for blunt-body flow computations

The formulation computation method of an improved version of the three-dimensional Navier-Stokes/Euler computation algorithm of Li (1981) for the numerical simulation of blunt-body reentry flows are discussed, and results for five sample problems are presented graphically. The vector notation of the coordinate systems is defined; the governing equations are presented in full; the Jacobian matrices and damping terms of the factorization technique (based on the alternating-direction implicit procedure of Beam and Warming, 1978) are explained; and the capabilities, limitations and proper use of the code are summarized. Examples presented include the equilibrium-flow problem for Shuttle-orbiter reentry at Mach 22 and angle of attack 40.8 deg and the perfect-gas problem of an aerobraking orbital-transfer vehicle with an ellipsoidal/60-deg cone and a toroidal sonic shoulder. The advantages of the improved code in terms of accuracy and computation time are indicated.

Li, C. P.↗
NASA NTRS1984

Research in applied mathematics, numerical analysis, and computer science

Research conducted at the Institute for Computer Applications in Science and Engineering (ICASE) in applied mathematics, numerical analysis, and computer science is summarized and abstracts of published reports are presented. The major categories of the ICASE research program are: (1) numerical methods, with particular emphasis on the development and analysis of basic numerical algorithms; (2) control and parameter identification; (3) computational problems in engineering and the physical sciences, particularly fluid dynamics, acoustics, and structural analysis; and (4) computer systems and software, especially vector and parallel computers.

Source record↗
DOE OSTI2019

On Xenon Fission Product Poisoning

We all appreciate the special role Xenon-135 plays in reactor transients. Unless there is sufficient reactivity override, on shutdown a conventional reactor cannot immediately restart because of the buildup of strongly absorbing Xenon from Iodine decay and the absence of neutron flux to continue its elimination. We are also aware of how the flux imbalance in a large reactor core gives rise to Xenon induced power precession around the core in unpredictable ways. In addition, with regard to computational methods, we lack readily accessible numerical algorithms enabling reliable characterization of Xenon oscillations with simple models to explain this fundamentally important phenomenon. Nowhere is this more apparent than in the classroom. In this presentation, we first consider a new 12-element comprehensive decay to model Xenon and its corresponding numerical implementation and verification to establish a predictive model of Xenon buildup in a constant flux. We then follow with a reduced self-consistent flux model in an infinite reactor- a model not found in textbooks. This presentation should appeal to those interested in developing high order numerical methods for reactor transients.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2011

High-Performance Data Analysis Tools for Sun-Earth Connection Missions

The data analysis tool of choice for many Sun-Earth Connection missions is the Interactive Data Language (IDL) by ITT VIS. The increasing amount of data produced by these missions and the increasing complexity of image processing algorithms requires access to higher computing power. Parallel computing is a cost-effective way to increase the speed of computation, but algorithms oftentimes have to be modified to take advantage of parallel systems. Enhancing IDL to work on clusters gives scientists access to increased performance in a familiar programming environment. The goal of this project was to enable IDL applications to benefit from both computing clusters as well as graphics processing units (GPUs) for accelerating data analysis tasks. The tool suite developed in this project enables scientists now to solve demanding data analysis problems in IDL that previously required specialized software, and it allows them to be solved orders of magnitude faster than on conventional PCs. The tool suite consists of three components: (1) TaskDL, a software tool that simplifies the creation and management of task farms, collections of tasks that can be processed independently and require only small amounts of data communication; (2) mpiDL, a tool that allows IDL developers to use the Message Passing Interface (MPI) inside IDL for problems that require large amounts of data to be exchanged among multiple processors; and (3) GPULib, a tool that simplifies the use of GPUs as mathematical coprocessors from within IDL. mpiDL is unique in its support for the full MPI standard and its support of a broad range of MPI implementations. GPULib is unique in enabling users to take advantage of an inexpensive piece of hardware, possibly already installed in their computer, and achieve orders of magnitude faster execution time for numerically complex algorithms. TaskDL enables the simple setup and management of task farms on compute clusters. The products developed in this project have the potential to interact, so one can build a cluster of PCs, each equipped with a GPU, and use mpiDL to communicate between the nodes and GPULib to accelerate the computations on each node.

Messmer, Peter↗
NASA NTRS1991

An Explicit Upwind Algorithm for Solving the Parabolized Navier-Stokes Equations

An explicit, upwind algorithm was developed for the direct (noniterative) integration of the 3-D Parabolized Navier-Stokes (PNS) equations in a generalized coordinate system. The new algorithm uses upwind approximations of the numerical fluxes for the pressure and convection terms obtained by combining flux difference splittings (FDS) formed from the solution of an approximate Riemann (RP). The approximate RP is solved using an extension of the method developed by Roe for steady supersonic flow of an ideal gas. Roe's method is extended for use with the 3-D PNS equations expressed in generalized coordinates and to include Vigneron's technique of splitting the streamwise pressure gradient. The difficulty associated with applying Roe's scheme in the subsonic region is overcome. The second-order upwind differencing of the flux derivatives are obtained by adding FDS to either an original forward or backward differencing of the flux derivative. This approach is used to modify an explicit MacCormack differencing scheme into an upwind differencing scheme. The second order upwind flux approximations, applied with flux limiters, provide a method for numerically capturing shocks without the need for additional artificial damping terms which require adjustment by the user. In addition, a cubic equation is derived for determining Vegneron's pressure splitting coefficient using the updated streamwise flux vector. Decoding the streamwise flux vector with the updated value of Vigneron's pressure splitting improves the stability of the scheme. The new algorithm is applied to 2-D and 3-D supersonic and hypersonic laminar flow test cases. Results are presented for the experimental studies of Holden and of Tracy. In addition, a flow field solution is presented for a generic hypersonic aircraft at a Mach number of 24.5 and angle of attack of 1 degree. The computed results compare well to both experimental data and numerical results from other algorithms. Computational times required for the upwind PNS code are approximately equal to an explicit PNS MacCormack's code and existing implicit PNS solvers.

Korte, John J.↗
NASA NTRS1974

Eigenproblem solution of damped structural systems

This article presents an efficient numerical algorithm and a complete listing of the associated computer program, developed for the eigenproblem solution of discrete damped structures, including spinning ones. The numerically stable procedure is based on a combined Sturm sequence and inverse iteration technique, which fully exploits the banded form of the relevant matrices and proves to be most substantially economical when compared to similar existing softwares. Numerical results are presented for representative structures, solved by the present computer program written in FORTRAN V for the JPL UNIVAC 1108 computer.

Gupta, K. K.↗
DOE OSTI2014

Robust parallel iterative solvers for linear and least-squares problems (Final Technical Report)

The primary goal of this project is to study and develop robust iterative methods for solving linear systems of equations and least squares systems. The focus of the Minnesota team is on algorithms development, robustness issues, and on tests and validation of the methods on realistic problems. 1. The project begun with an investigation on how to practically update a preconditioner obtained from an ILU-type factorization, when the coefficient matrix changes. 2. We investigated strategies to improve robustness in parallel preconditioners in a specific case of a PDE with discontinuous coefficients. 3. We explored ways to adapt standard preconditioners for solving linear systems arising from the Helmholtz equation. These are often difficult linear systems to solve by iterative methods. 4. We have also worked on purely theoretical issues related to the analysis of Krylov subspace methods for linear systems. 5. We developed an effective strategy for performing ILU factorizations for the case when the matrix is highly indefinite. The strategy uses shifting in some optimal way. The method was extended to the solution of Helmholtz equations by using complex shifts, yielding very good results in many cases. 6. We addressed the difficult problem of preconditioning sparse systems of equations on GPUs. 7. A by-product of the above work is a software package consisting of an iterative solver library for GPUs based on CUDA. This was made publicly available. It was the first such library that offers complete iterative solvers for GPUs. 8. We considered another form of ILU which blends coarsening techniques from Multigrid with algebraic multilevel methods. 9. We have released a new version on our parallel solver - called pARMS [new version is version 3]. As part of this we have tested the code in complex settings - including the solution of Maxwell and Helmholtz equations and for a problem of crystal growth.10. As an application of polynomial preconditioning we considered the problem of evaluating f(A)v which arises in statistical sampling. 11. As an application to the methods we developed, we tackled the problem of computing the diagonal of the inverse of a matrix. This arises in statistical applications as well as in many applications in physics. We explored probing methods as well as domain-decomposition type methods. 12. A collaboration with researchers from Toulouse, France, considered the important problem of computing the Schur complement in a domain-decomposition approach. 13. We explored new ways of preconditioning linear systems, based on low-rank approximations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Effective construction of eigenvectors for a class of singular sparse matrices

Fundamental results and an efficient algorithm for constructing eigenvectors corresponding to non-zero eigenvalues of matrices with zero rows and/or columns are developed. The formulation is built on the relation between eigenvectors of such matrices and the eigenvectors of their submatrices after removing all zero rows and columns. While being easily utilized, the algorithm decreases the computation time needed for numerical eigenanalysis, and resolves potential numerical eigensolver instabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A conservative Galerkin solver for the quasilinear diffusion model in magnetized plasmas

We propose a conservative Galerkin scheme for the quasilinear model in three-dimensional momentum space and three-dimensional spectral space, with cylindrical symmetry. We construct an unconditionally conservative weak form and use a discretization that preserves conservation properties independent of the wave emission probability. The discrete operators, combined with a consistent quadrature rule, preserve all the conservation laws rigorously. The proposed scheme is quite general: it works for both relativistic and non-relativistic systems, for both magnetized and unmagnetized plasmas, and even for problems with time-dependent dispersion relations. We represent the particle distribution by continuous basis functions and use discontinuous basis functions for the wave spectral energy density, which enables the application of a positivity-preserving technique. We adopt the marching simplex algorithm, designed initially for computer graphics, for numerical integration on the resonance manifold. Furthermore, the numerical examples with a bump-on-tail initial configuration show how the unstable waves produce strong momentum space diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Highly stable evolution of Earth's future orbit despite chaotic behavior of the solar system

Due to the chaotic nature of the solar system, the question of its dynamic long-term stability can only be answered in a statistical sense, for instance, based on numerical ensemble integrations of nearby orbits. Destabilization of the inner planets, including catastrophic encounters and/or collisions involving the Earth, has been suggested to be initiated through a large increase in Mercury’s eccentricity (e{sub M}), with an estimated probability of ∼1%. However, it has recently been shown that the statistics of numerical solar system integrations are sensitive to the accuracy and type of numerical algorithm. Here, I report results from computationally demanding ensemble integrations (N = 1600 with slightly different initial conditions) at unprecedented accuracy based on the full equations of motion of the eight planets and Pluto over 5 Gyr, including contributions from general relativity. The standard symplectic algorithm used for long-term integrations produced spurious results for highly eccentric orbits and during close encounters, which were hence integrated with a suitable Bulirsch–Stoer algorithm, specifically designed for these situations. The present study yields odds for a large increase in Mercury’s eccentricity that are less than previous estimates. Strikingly, in two solutions, Mercury continued on highly eccentric orbits (after reaching e{sub M} values >0.93) for 80–100 Myr before colliding with Venus or the Sun. Most importantly, none of the 1600 solutions led to a close encounter involving the Earth or a destabilization of Earth’s orbit in the future. I conclude that Earth’s orbit will be dynamically highly stable for billions of years in the future, despite the chaotic behavior of the solar system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

BoBa

BoBa is a C++ software library for working with large matrices, tensors, and tensor decompositions. The library provides tools for dense matrix and tensor operations, tensor decompositions, and tensor decomposition methods that support modern CPU and GPU architectures. It includes portable abstractions for linear algebra, tensor algebra, and multidimensional computation. BoBa is intended for scientific computing applications that involve large multidimensional data sets or high dimensional mathematical models. Its capabilities support tasks such as data compression, linear algebra, efficient numerical computation, and the development of scalable algorithms for heterogeneous hardware. Tutorials, tests, and example applications are included to help users learn and apply the library.

Yao, Jin [Lawrence Livermore National Laboratory (↗
NASA NTRS1971

Theoretical buckling loads of boron/aluminum and graphite/resin fiber composite anisotropic plates

Theoretical results are presented for the buckling of anisotropic plates. The plates are subjected to simple and combined in-plane loading. The plates are made from fiber composite material of boron/aluminum or high-modulus graphite/resin. The results are presented in nondimensional form as buckling load against fiber orientation angle for various plate aspect ratios. The results indicate that buckling loads of boron/aluminum plates are independent of fiber direction if the plate aspect ratios are greater than about 1, and moderately dependent when this ratio is less than about 1. In addition, the results indicate that the buckling loads are independent of aspect ratio for plates with aspect ratios greater than about 2. Boron/ aluminum composite plates can resist buckling loads more efficiently than graphite/resin composites on a specific buckling stress basis. The numerical algorithm and a listing of the computer code used to obtain the results are included.

Chamis, C. C.↗
NASA NTRS1981

Studies in astronomical time series analysis. I - Modeling random processes in the time domain

Several random process models in the time domain are defined and discussed. Attention is given to the moving average model, the autoregressive model, and relationships between and combinations of these models. Consideration is then given to methods for investigating pulse structure, procedures of model construction, computational methods, and numerical experiments. A FORTRAN algorithm of time series analysis has been developed which is relatively stable numerically. Results of test cases are given to study the effect of adding noise and of different distributions for the pulse amplitudes. A preliminary analysis of the light curve of the quasar 3C 272 is considered as an example.

Scargle, J. D.↗
NASA NTRS1995

Spurious Numerical Solutions Of Differential Equations

Paper presents detailed study of spurious steady-state numerical solutions of differential equations that contain nonlinear source terms. Main objectives of this study are (1) to investigate how well numerical steady-state solutions of model nonlinear reaction/convection boundary-value problem mimic true steady-state solutions and (2) to relate findings of this investigation to implications for interpretation of numerical results from computational-fluid-dynamics algorithms and computer codes used to simulate reacting flows.

Lafon, A.↗
NASA NTRS2009

Analysis of Mass Averaged Tissue Doses in CAM, CAF, MAX, and FAX

To estimate astronaut health risk due to space radiation, one must have the ability to calculate exposure-related quantities averaged over specific organs and tissue types. In this study, we first examine the anatomical properties of the Computerized Anatomical Man (CAM), Computerized Anatomical Female (CAF), Male Adult voXel (MAX), and Female Adult voXel (FAX) models by comparing the masses of various tissues to the reference values specified by the International Commission on Radiological Protection (ICRP). Major discrepancies are found between the CAM and CAF tissue masses and the ICRP reference data for almost all of the tissues. We next examine the distribution of target points used with the deterministic transport code HZETRN to compute mass averaged exposure quantities. A numerical algorithm is used to generate multiple point distributions for many of the effective dose tissues identified in CAM, CAF, MAX, and FAX. It is concluded that the previously published CAM and CAF point distributions were under-sampled and that the set of point distributions presented here should be adequate for future studies involving CAM, CAF, MAX, or FAX. It is concluded that MAX and FAX are more accurate than CAM and CAF for space radiation analyses.

Slaba, Tony C.↗
NASA NTRS2007

NASA Tech Briefs, November 2007

Topics include: Wireless Measurement of Contact and Motion Between Contact Surfaces; Wireless Measurement of Rotation and Displacement Rate; Portable Microleak-Detection System; Free-to-Roll Testing of Airplane Models in Wind Tunnels; Cryogenic Shrouds for Testing Thermal-Insulation Panels; Optoelectronic System Measures Distances to Multiple Targets; Tachometers Derived From a Brushless DC Motor; Algorithm-Based Fault Tolerance for Numerical Subroutines; Computational Support for Technology- Investment Decisions; DSN Resource Scheduling; Distributed Operations Planning; Phase-Oriented Gear Systems; Freeze Tape Casting of Functionally Graded Porous Ceramics; Electrophoretic Deposition on Porous Non- Conductors; Two Devices for Removing Sludge From Bioreactor Wastewater; Portable Unit for Metabolic Analysis; Flash Diffusivity Technique Applied to Individual Fibers; System for Thermal Imaging of Hot Moving Objects; Large Solar-Rejection Filter; Improved Readout Scheme for SQUID-Based Thermometry; Error Rates and Channel Capacities in Multipulse PPM; Two Mathematical Models of Nonlinear Vibrations; Simpler Adaptive Selection of Golomb Power-of- Two Codes; VCO PLL Frequency Synthesizers for Spacecraft Transponders; Wide Tuning Capability for Spacecraft Transponders; Adaptive Deadband Synchronization for a Spacecraft Formation; Analysis of Performance of Stereoscopic-Vision Software; Estimating the Inertia Matrix of a Spacecraft; Spatial Coverage Planning for Exploration Robots; and Increasing the Life of a Xenon-Ion Spacecraft Thruster.

Source record↗
DOE OSTI2017

Investigation of Xenon Oscillations via Point Kinetics

We all appreciate the special role Xenon-135 plays in reactor transients. Unless there is sufficient reactivity override, on shutdown a conventional reactor cannot immediately restart because of the buildup of strongly absorbing Xenon from Iodine decay and the absence of neutron flux to continue its elimination. We are also aware of how the flux imbalance in a large reactor core gives rise to Xenon induced power precession around the core in unpredictable ways. In addition, with regard to computational methods, we lack readily accessible numerical algorithms enabling reliable characterization of Xenon oscillations with simple models to explain this fundamentally important phenomenon. Nowhere is this more apparent than in the classroom. In this presentation, we first consider a new 12-element comprehensive decay model to Xenon and its corresponding numerical implementation and verification to establish a predictive model of Xenon buildup in a constant flux. We then follow with a reduced self-consistent flux model in an infinite reactor- a model not found in textbooks. This presentation should appeal to those interested in developing high order numerical methods for reactor transients.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1982

Mixed time integration methods for transient thermal analysis of structures

The computational methods used to predict and optimize the thermal structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a different yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. The computer implementation aspects and numerical evaluation of these mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally useful method of estimating the critical time step for linear quadrilateral element is also given. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.

Liu, W. K.↗