SEARCH · Search NASA

Results for “Methods: numerical”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Hypersonic shock structure with Burnett terms in the viscous stress and heat flux

The continuum Navier-Stokes and Burnett equations are solved for one-dimensional shock structure in various monatomic gases. A new numerical method is employed which utilizes the complete time-dependent continuum equations and obtains the steady-state shock structure by allowing the system to relax from arbitrary initial conditions. Included is discussion of numerical difficulties encountered when solving the Burnett equations. Continuum solutions are compared to those obtained utilizing the Direct Simulation Monte Carlo method. Shock solutions are obtained for a hard sphere gas and for argon from Mach 1.3 to Mach 50. Solutions for a Maxwellian gas are obtained from Mach 1.3 to Mach 3.8. It is shown that the Burnett equations yield shock structure solutions in much closer agreement to both Monte Carlo and experimental results than do the Navier-Stokes equations. Shock density thickness, density asymmetry, and density-temperature separation are all more accurately predicted by the Burnett equations than by the Navier-Stokes equations.

Chapman, Dean R.↗

NASA NTRS1986

Methods for the solution of radiative transfer equation

To obtain an exact solution of the radiative-transfer equation in media where both absorption and scattering are significant, the usual approach is to use a numerical method. Three methods are known in the literature: invariant imbedding, eigenvalue-eigenfunction, and matrix doubling. This paper examines the practical application of these methods to the problem of emission from an inhomogeneous (Rayleigh) layer, the effects of layer parameters on the stability. It is found that invariant imbedding is most suitable for computing emission from an inhomogeneous layer with a temperature profile but tends to be unstable as the optical thickness of the layer increases beyond 0.5. On the other hand, the matrix-doubling method is stable for arbitrary optical thickness but is not suitable for handling multilayers. The eigenvalue-eigenfunction method is more stable than the invariant imbedding as optical thickness increases up to 2.0. It also permits temperature profile in the layer, but the computation is much more complicated. It is less stable than the matrix-doubling method when optical thickness is larger than 2.0. In general, the choice of a method is dependent on the nature of the problem.

Chen, M. F.↗

NASA NTRS1978

Transonic wing analysis using advanced computational methods

This paper discusses the application of three-dimensional computational transonic flow methods to several different types of transport wing designs. The purpose of these applications is to evaluate the basic accuracy and limitations associated with such numerical methods. The use of such computational methods for practical engineering problems can only be justified after favorable evaluations are completed. The paper summarizes a study of both the small-disturbance and the full potential technique for computing three-dimensional transonic flows. Computed three-dimensional results are compared to both experimental measurements and theoretical results. Comparisons are made not only of pressure distributions but also of lift and drag forces. Transonic drag rise characteristics are compared. Three-dimensional pressure distributions and aerodynamic forces, computed from the full potential solution, compare reasonably well with experimental results for a wide range of configurations and flow conditions.

Henne, P. A.↗

NASA NTRS1983

Unsteady laminar boundary-layer calculations on oscillating configurations including backflow. Part 1: Flat plate, oscillating in its own plane

A finite difference method has been developed to calculate the unsteady boundary layer over an oscillating flat plate. Low- and high frequency approximations were used for comparison with numerical results. Special emphasis was placed on the behavior of the flow and on the numerical calculation procedure as soon as reversed flow has occurred over part of the oscillation cycle. The numerical method displayed neither problems nor singular behavior at the beginning of or within the reversed flow region. Calculations, however, came to a limit where the back-flow region reached the plate's leading edge in the case of high oscillation amplitudes. It is assumed that this limit is caused by the special behavior of the flow at the plate's leading edge where the boundary layer equations are not valid.

Geissler, W.↗

NASA NTRS1986

Partial spline models for the inclusion of tropopause and frontal boundary information in otherwise smooth two- and three-dimensional objective analysis

A new method, based on partial spline models, is developed for including specified discontinuities in otherwise smooth two- and three-dimensional objective analyses. The method is appropriate for including tropopause height information in two- and three-dimensinal temperature analyses, using the O'Sullivan-Wahba physical variational method for analysis of satellite radiance data, and may in principle be used in a combined variational analysis of observed, forecast, and climate information. A numerical method for its implementation is described and a prototype two-dimensional analysis based on simulated radiosonde and tropopause height data is shown. The method may also be appropriate for other geophysical problems, such as modeling the ocean thermocline, fronts, discontinuities, etc.

Shiau, Jyh-Jen↗

NASA NTRS1992

Schroedinger's radial equation - Solution by extrapolation

A high-accuracy numerical method for the solution of a 1D Schroedinger equation that is suitable for a diatomic molecule, obtained by combining a finite-difference method with iterative extrapolation to the limit, is presently shown to have several advantages over more conventional methods. Initial guesses for the term values are obviated, and implementation of the algorithm is straightforward. The method is both less sensitive to round-off error, and faster than conventional methods for equivalent accuracy. These advantages are illustrated through the solution of Schroedinger's equation for a Morse potential function suited for HCl and a numerically derived Rydberg-Klein-Rees potential function for the X 1Sigma(+) state of CO.

Goorvitch, D.↗

NASA NTRS2002

Discontinuous Galerkin Methods for Turbulence Simulation

A discontinuous Galerkin (DG) method is formulated, implemented, and tested for simulation of compressible turbulent flows. The method is applied to turbulent channel flow at low Reynolds number, where it is found to successfully predict low-order statistics with fewer degrees of freedom than traditional numerical methods. This reduction is achieved by utilizing local hp-refinement such that the computational grid is refined simultaneously in all three spatial coordinates with decreasing distance from the wall. Another advantage of DG is that Dirichlet boundary conditions can be enforced weakly through integrals of the numerical fluxes. Both for a model advection-diffusion problem and for turbulent channel flow, weak enforcement of wall boundaries is found to improve results at low resolution. Such weak boundary conditions may play a pivotal role in wall modeling for large-eddy simulation.

Collis, S. Scott↗

NASA NTRS1977

Error behavior of multistep methods applied to unstable differential systems

The problem of modeling a dynamic system described by a system of ordinary differential equations which has unstable components for limited periods of time is discussed. It is shown that the global error in a multistep numerical method is the solution to a difference equation initial value problem, and the approximate solution is given for several popular multistep integration formulas. Inspection of the solution leads to the formulation of four criteria for integrators appropriate to unstable problems. A sample problem is solved numerically using three popular formulas and two different stepsizes to illustrate the appropriateness of the criteria.

Brown, R. L.↗

NASA NTRS1978

Error behaviour of multistep methods applied to unstable differential systems

The problem of modelling a dynamic system described by a system of ordinary differential equations which has unstable components for limited periods of time is discussed. It is shown that the global error in a multistep numerical method is the solution to a difference equation initial value problem, and the approximate solution is given for several popular multistep integration formulae. Inspection of the solution leads to the formulation of four criteria for integrators appropriate to unstable problems. A sample problem is solved numerically using three popular formulae and two different stepsizes to illustrate the appropriateness of the criteria.

Brown, R. L.↗

DOE OSTI2023

Variational Monte Carlo calculations of n + H 3 scattering

A paramount goal in nuclear physics is to unify ab initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound-state calculations in light systems but have seen minimal implementation in unbound systems. Here, we introduce a numerical method to improve the efficiency and accuracy of unbound-state calculations in QMC. As an initial application, we compute scattering observables for the smallest system available to probe three-body forces, the neutron-triton system, using variational Monte Carlo (VMC) wave functions. The method involves inferring long-range amplitudes in the wave function from integrals over the short-range region where all the particles interact. This approach using integral relations is well established in the literature; here, we develop it for the QMC framework. We validate our approach with a consistency check between short-range spectroscopic overlap functions computed from direct evaluation and from the integral relations; scattering amplitudes are long-range asymptotics of those overlaps. Comparison against published benchmark calculations using the same potential demonstrates that when applied to the current VMC wave functions, the integral method produces more accurate scattering observables than direct evaluation from the same variational wave function. However, it still differs noticeably from the exact results. Using additional interactions, we then present phase shifts and mixing parameters for the n + 3 H system. In particular, we present one of the first applications of the Norfolk family of local coordinate-space chiral potentials in unbound systems of A > 2. The Norfolk results accurately describe s-wave scattering but predict p-wave cross sections too large. Compared with previous QMC scattering calculations, the integral method avoids difficulties associated with the precise computation of energy differences and with convergence outside the interaction region, which is particularly severe in the variational calculation. Application of the integral method here paves the way for its use in Green's function Monte Carlo (GFMC) calculations. In GFMC, the wave functions are more accurate, but the high-precision convergence of their tails is slow, and there are additional difficulties in reading out amplitudes. Here, the integral methods will address both of those remaining problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1987

Response of joint dominated space structures

An approximate method is developed for estimating the transient response of nonlinear systems in terms of linearized modes of response. Its advantages are that it is computationally more efficient than the time integration method and that it is possible to view the design problem in the more traditional physical terms of modal response. The major drawback of the approximate method is loss of accuracy. It seems that both approximate methods and time integration have their roles in design. Approximate methods provide efficient tools for performing parametric studies and they supply physical insights into how to optimize system performance that are not easily inferred from strictly numerical methods. Time integration provides a method for assessing the accuracy of the approximate solution for key simulations and for fine tuning the final design. In the procedure presented the nonlinear system is approximated by an equivalent linear system in which the system parameters are constant over the range of transient response.

Source record↗

NASA NTRS1987

Numerical simulations of unsteady airfoil-vortex interactions

Numerical methods based on the thin layer Navier-Stokes, Euler and transonic small disturbance equations, all three of which use the perturbation form of an implicit numerical algorithm, are used to formulate the unsteady interactions of a moving vortex with a stationary airfoil. Results are presented for both transonic and subsonic flows, as well as for weak and strong interactions. In general, the results exhibit an overwhelming influence of the vortex on the flow field around the airfoil, especially when the convecting vortex is within one chord of the airfoil. It is noted that even the ostensibly weak interactions produced vortex/shock-induced boundary layer separation.

Srinivasan, G. R.↗

NASA NTRS1990

Shape optimization for maximum stability and dynamic stiffness

Any optimization of structures for maximum stability or for maximum dynamic stiffness deals with an eigenvalue problem. The goal of this optimization is to raise the lowest eigenvalue (or eigenvalues) of the problem to its highest (optimal) level at a constant volume of the structure. Likely the lowest eigenvalue may be either inherently multi-modal or it can become multi-modal as a result of the optimization process. The multimodeness introduces some ambiguity to the eigenvalue problem and make the optimization difficult to handle. Thus far, only the simplest cases of multi-modal structures have been effectively optimized using rather elaborate analytical methods. Numerous publications report design of a minimum volume structure with different eigenvalues constraints, in which, however, the modality of the problem is assumed a priori. The method presented here utilizes a multi-modal optimality criteria and allows for inclusion of an arbitrary number of buckling or vibrations modes which might influence the optimization process. The real multi-modality of the problem, that is the number of modes participating in the final optimal design is determined iteratively. Because of a natural use of the FEM technique the method is easy to program and might be helpful in design of large flexible space structures.

Szyszkowski, W.↗

NASA NTRS1997

Multiresolution and Explicit Methods for Vector Field Analysis and Visualization

This is a request for a second renewal (3d year of funding) of a research project on the topic of multiresolution and explicit methods for vector field analysis and visualization. In this report, we describe the progress made on this research project during the second year and give a statement of the planned research for the third year. There are two aspects to this research project. The first is concerned with the development of techniques for computing tangent curves for use in visualizing flow fields. The second aspect of the research project is concerned with the development of multiresolution methods for curvilinear grids and their use as tools for visualization, analysis and archiving of flow data. We report on our work on the development of numerical methods for tangent curve computation first.

Nielson, Gregory M.↗

NASA NTRS2013

Why Is the Overheating Problem Difficult: the Role of Entropy

The development of computational fluid dynamics over the last few decades has yielded enormous successes and capabilities being routinely employed today; however there remain some open problems to be properly resolved-some are fundamental in nature and some resolvable by operational changes. These two categories are distinguished and broadly explored previously. One, that belongs to the former, is the so-called overheating problem, especially in rarefying flow. This problem up to date still dogs every method known to the author; a solution to it remains elusive. The study in this paper concludes that: (1) the entropy increase is quantitatively linked to the increase in the temperature increase, (2) it is argued that the overheating is inevitable in the current shock capturing or traditional finite difference framework, and (3) a simple hybrid method is proposed that removes the overheating problem in the rarefying problems, but also retains the property of accurate shock capturing. This remedy (enhancement of current numerical methods) can be included easily in the present Eulerian codes.

Liou, Meng-Sing↗

NASA NTRS2017

Improving Estimation of Ground Casualty Risk From Reentering Space Objects

A recent improvement to the long-term estimation of ground casualties from reentering space debris is the further refinement and update to the human population distribution. Previous human population distributions were based on global totals with simple scaling factors for future years, or a coarse grid of population counts in a subset of the world's countries, each cell having its own projected growth rate. The newest population model includes a 5-fold refinement in both latitude and longitude resolution. All areas along a single latitude are combined to form a global population distribution as a function of latitude, creating a more accurate population estimation based on non-uniform growth at the country and area levels. Previous risk probability calculations used simplifying assumptions that did not account for the ellipsoidal nature of the Earth. The new method uses first, a simple analytical method to estimate the amount of time spent above each latitude band for a debris object with a given orbit inclination and second, a more complex numerical method that incorporates the effects of a non-spherical Earth. These new results are compared with the prior models to assess the magnitude of the effects on reentry casualty risk.

Ostrom, Chris L.↗

NASA NTRS2017

Improving Estimation of Ground Casualty Risk from Reentering Space Objects

A recent improvement to the long-term estimation of ground casualties from reentering space debris is the further refinement and update to the human population distribution. Previous human population distributions were based on global totals with simple scaling factors for future years, or a coarse grid of population counts in a subset of the world's countries, each cell having its own projected growth rate. The newest population model includes a 5-fold refinement in both latitude and longitude resolution. All areas along a single latitude are combined to form a global population distribution as a function of latitude, creating a more accurate population estimation based on non-uniform growth at the country and area levels. Previous risk probability calculations used simplifying assumptions that did not account for the ellipsoidal nature of the earth. The new method uses first, a simple analytical method to estimate the amount of time spent above each latitude band for a debris object with a given orbit inclination, and second, a more complex numerical method that incorporates the effects of a non-spherical Earth. These new results are compared with the prior models to assess the magnitude of the effects on reentry casualty risk.

Ostrom, C.↗

NASA NTRS1989

Projection methods for the numerical solution of Markov chain models

Projection methods for computing stationary probability distributions for Markov chain models are presented. A general projection method is a method which seeks an approximation from a subspace of small dimension to the original problem. Thus, the original matrix problem of size N is approximated by one of dimension m, typically much smaller than N. A particularly successful class of methods based on this principle is that of Krylov subspace methods which utilize subspaces of the form span(v,av,...,A(exp m-1)v). These methods are effective in solving linear systems and eigenvalue problems (Lanczos, Arnoldi,...) as well as nonlinear equations. They can be combined with more traditional iterative methods such as successive overrelaxation, symmetric successive overrelaxation, or with incomplete factorization methods to enhance convergence.

Saad, Youcef↗