Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Methods”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

A comparison of the efficiency of numerical methods for integrating chemical kinetic rate equations

The efficiency of several algorithms used for numerical integration of stiff ordinary differential equations was compared. The methods examined included two general purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes were applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code available for the integration of combustion kinetic rate equations. It is shown that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient then evaluating the temperature by integrating its time-derivative.

Radhakrishnan, K.↗
NASA NTRS1975

Numerical methods for the calculation of three-dimensional nozzle exhaust flow fields

Numerical codes developed for the calculation of three-dimensional nozzle exhaust flow fields associated with hypersonic airbreathing aircraft are described. Both codes employ reference plane grid networks with respect to three coordinate systems. Program CHAR3D is a characteristic code utilizing a new wave preserving network within the reference planes, while program BIGMAC is a finite difference code utilizing conservation variables and a one-sided difference algorithm. Secondary waves are numerically captured by both codes, while the underexpansion shock and plume boundary are treated discretely. The exhaust gas properties consist of hydrogen-air combustion product mixtures in local chemical equilibrium. Nozzle contours are treated by a newly developed geometry package based on dual cubic splines. Results are presented for simple configurations demonstrating two- and three-dimensional multiple wave interactions.

Dash, S. M.↗
NASA NTRS1984

Finite-analytic numerical method for unsteady two-dimensional Navier-Stokes equations

A finite analytic (FA) numerical solution is developed for unsteady two-dimensional Navier-Stokes equations. The FA method utilizes the analytic solution in a small local element to formulate the algebraic representation of partial differential equations. The combination of linear and exponential functions that satisfy the governing equation is adopted as the boundary function, thereby improving the accuracy of the finite analytic solution. Two flows, one a starting cavity flow and the other a vortex shedding flow behind a rectangular block, are solved by the FA method. The starting square cavity flow is solved for Reynolds number of 400, 1000, and 2000 to show the accuracy and stability of the FA solution. The FA solution for flow over a rectangular block (H x H/4) predicts the Strouhal number for Reynolds numbers of 100 and 500 to be 0.156 and 0.125. Details of the flow patterns are given. In addition to streamlines and vorticity distribution, rest-streamlines are given to illustrate the vortex motion downstream of the block.

Chen, C.-J.↗
NASA NTRS1997

Comparison of Several Numerical Methods for Simulation of Compressible Shear Layers

An investigation is conducted on several numerical schemes for use in the computation of two-dimensional, spatially evolving, laminar variable-density compressible shear layers. Schemes with various temporal accuracies and arbitrary spatial accuracy for both inviscid and viscous terms are presented and analyzed. All integration schemes use explicit or compact finite-difference derivative operators. Three classes of schemes are considered: an extension of MacCormack's original second-order temporally accurate method, a new third-order variant of the schemes proposed by Rusanov and by Kutier, Lomax, and Warming (RKLW), and third- and fourth-order Runge-Kutta schemes. In each scheme, stability and formal accuracy are considered for the interior operators on the convection-diffusion equation U(sub t) + aU(sub x) = alpha U(sub xx). Accuracy is also verified on the nonlinear problem, U(sub t) + F(sub x) = 0. Numerical treatments of various orders of accuracy are chosen and evaluated for asymptotic stability. Formally accurate boundary conditions are derived for several sixth- and eighth-order central-difference schemes. Damping of high wave-number data is accomplished with explicit filters of arbitrary order. Several schemes are used to compute variable-density compressible shear layers, where regions of large gradients exist.

Kennedy, Christopher A.↗
NASA NTRS1972

A direct numerical method for predicting concentration profiles in a turbulent boundary layer over a flat plate

A numerical solution of the turbulent mass transport equation utilizing the concept of eddy diffusivity is presented as an efficient method of investigating turbulent mass transport in boundary layer type flows. A FORTRAN computer program is used to study the two-dimensional diffusion of ammonia, from a line source on the surface, into a turbulent boundary layer over a flat plate. The results of the numerical solution are compared with experimental data to verify the results of the solution. Several other solutions to diffusion problems are presented to illustrate the versatility of the computer program and to provide some insight into the problem of mass diffusion as a whole.

Dow, J. W.↗
NASA NTRS1982

A new numerical method for the simulation of three dimensional flow in a pipe

A new numerical technique for simulating three dimensional, unsteady, incompressible pipe flows is presented and its utility and accuracy is shown. Each vector function in the expansion of the velocity field is divergence free and satisfies the boundary conditions for viscous flow. Some of the benefits of the expansion technique are that pressure is eliminated from the dynamics, only two unknowns per mesh point are required, implicit treatment of the viscous terms is provided at no extra computational cost, and no fractional time steps are required. The method uses spectral expansions: Fourier series in the azimuthal and streamwise directions, and Jacobi polynominals in the radial direction.

Leonard, A.↗
NASA NTRS1982

A new numerical method for the simulation of three-dimensional flow in a pipe

A new numerical technique for simulating three dimensional, unsteady, incompressible pipe flows is presented and its utility and accuracy is shown. Each vector function in the expansion of the velocity field is divergence free and satisfies the boundary conditions for viscous flow. Some of the benefits of the expansion technique are that pressure is eliminated from the dynamics, only two unknowns per mesh point are required, implicit treatment of the viscous terms is provided at no extra computational cost, and no fractional time steps are required. The method uses spectral expansions: Fourier series in the azimuthal and streamwise directions, and Jacobi polynominals in the radial direction. Previously announced in STAR as N82-31644

Leonard, A.↗
NASA NTRS1999

Acoustic Treatment Design Scaling Methods: Numerical Simulation of the Nonlinear Acoustic Impedance of a Perforated Plate Single-Degree-of-Freedom Resonator Using a Time-Domain Finite Difference Method - Volume 4

Single-degree-of-freedom resonators consisting of honeycomb cells covered by perforated facesheets are widely used as acoustic noise suppression liners in aircraft engine ducts. The acoustic resistance and mass reactance of such liners are known to vary with the intensity of the sound incident upon the panel. Since the pressure drop across a perforated liner facesheet increases quadratically with the flow velocity through the facesheet, this is known as the nonlinear resistance effect. In the past, two different empirical frequency domain models have been used to predict the Sound Pressure Level effect of the incident wave on the perforated liner impedance, one that uses the incident particle velocity in isolated narrowbands, and one that models the particle velocity as the overall velocity. In the absence of grazing flow, neither frequency domain model is entirely accurate in predicting the nonlinear effect that is measured for typical perforated sheets. The time domain model is developed in an attempt to understand and improve the model for the effect of spectral shape and amplitude of multi-frequency incident sound pressure on the liner impedance. A computer code for the time-domain finite difference model is developed and predictions using the models are compared to current frequency-domain models.

Kraft, R. E.↗
NASA NTRS1972

Numerical methods for solving the Boltzmann equation (a review)

The methods are reviewed which are utilized in principal attempts to obtain the numerical solution or modeling of the Boltzmann equation over a broad range of Knudsen numbers. The primary methods considered are the Monte Carlo and the discrete velocities methods. The conculsions drawn from the analysis include the following: (1) The Monte Carlo methods are not well suited in the area of small Knudsen numbers. (2) Among the Monte Carlo methods, the Bird method appears to be the most attractive, since it is more directly related to the Boltzmann equation. (3) The deterministic methods, which include the discrete ordinate technique, offer great possibilities but require exceedingly large computer times. (4) The use of approximating equations in combination with the discrete velocities method will possibly improve computation time and reduce the required memory volume.

Limar, Y. F.↗
NASA NTRS1984

A comparison of the efficiency of numerical methods for integrating chemical kinetic rate equations

A comparison of the efficiency of several algorithms recently developed for the efficient numerical integration of stiff ordinary differential equations is presented. The methods examined include two general-purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D, and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes are applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code currently available for the integration of combustion kinetic rate equations. An important finding is that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient than evaluating the temperature by integrating its time-derivative.

Radhakrishnan, K.↗
NASA NTRS1972

Numerical method and FORTRAN program for the solution of an axisymmetric electrostatic collector design problem

The numerical calculation is described of the steady-state flow of electrons in an axisymmetric, spherical, electrostatic collector for a range of boundary conditions. The trajectory equations of motion are solved alternately with Poisson's equation for the potential field until convergence is achieved. A direct (noniterative) numerical technique is used to obtain the solution to Poisson's equation. Space charge effects are included for initial current densities as large as 100 A/sq cm. Ways of dealing successfully with the difficulties associated with these high densities are discussed. A description of the mathematical model, a discussion of numerical techniques, results from two typical runs, and the FORTRAN computer program are included.

Reese, O. W.↗
NASA NTRS1978

Numerical method for boundary layers with blowing - The exponential box scheme

The paper describes a new numerical scheme based on exponential difference operator concepts combined with Keller's (1968) box scheme approach to produce a stable second-order accurate finite-difference scheme for convection-diffusion problems arising in boundary layer flows in the presence of massive injection through a porous surface. The technique is demonstrated by application to the self-similar boundary layer equations with massive blowing at the surface.

El-Mistikawy, T. M.↗
NASA NTRS1987

Aerothermal modeling program. Phase 2, element A: Improved numerical methods for turbulent viscous recirculating flows

The objective of this effort is to develop improved numerical schemes for predicting combustor flow fields. Various candidate numerical schemes were evaluated, and promising schemes were selected for detailed assessment. The criteria for evaluation included accuracy, computational efficiency, stability, and ease of extension to multidimensions. The candidate schemes were assessed against a variety of simple one- and two-dimensional problems. These results led to the selection of the following schemes for further evaluation: flux spline schemes (linear and cubic) and controlled numerical diffusion with internal feedback (CONDIF). The incorporation of the flux spline scheme and direct solution strategy in a computer program for three-dimensional flows is in progress.

Karki, K. C.↗
NASA NTRS1981

Forebody and base region real gas flow in severe planetary entry by a factored implicit numerical method. II - Equilibrium reactive gas

The factored-implicit, gas-dynamic algorithm has been adapted to the numerical simulation of equilibrium reactive flows. Changes required in the perfect gas version of the algorithm are developed, and the method of coupling gas-dynamic and chemistry variables is discussed. A flow-field solution that approximates a Jovian entry case was obtained by this method and compared with the same solution obtained by HYVIS, a computer program much used for the study of planetary entry. Comparison of surface pressure distribution and stagnation line shock-layer profiles indicates that the two solutions agree well.

Davy, W. C.↗
NASA NTRS1979

Numerical methods for problems involving the Drazin inverse

The objective was to try to develop a useful numerical algorithm for the Drazin inverse and to analyze the numerical aspects of the applications of the Drazin inverse relating to the study of homogeneous Markov chains and systems of linear differential equations with singular coefficient matrices. It is felt that all objectives were accomplished with a measurable degree of success.

Meyer, C. D., Jr.↗
NASA NTRS1978

Some recent advances in numerical methods in structural dynamics

The purpose of this paper is to present some recent developments in the area of numerical dynamic analysis of large structural systems. The first part of the paper is devoted to the finite dynamic element method, a new concept in the discretization of a continuum, the adoption of which effects considerable reduction in the number of degrees of freedom for a required solution accuracy, when compared with the usual finite element method. This is followed by a brief description of eigenproblem solution techniques, for the free vibration analysis of both rotating and nonrotating structures, that fully exploit the associated matrix sparsity. Together, these techniques constitute a powerful tool for the free vibration and subsequent dynamic response analysis of complex practical structures.

Gupta, K. K.↗
NASA NTRS1986

A numerical method for vortex sheet roll-up

The problem of computing vortex sheet roll-up from periodic analytic initial data is studied. Previous theoretical and numerical work is reviewed. Computational difficulties arising from ill posedness and singularity formation are discussed. A desingularization method is proposed to diminish these difficulties. Computations indicate that this approach converges past the time at which previous numerical investigations have failed to converge.

Krasny, R.↗