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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

A domain method for shape design sensitivity analysis of built-up structures

A numerical method of obtaining accurate shape design sensitivity information for built-up structures is developed and demonstrated through examples. The method utilizes the basic character of the finite element method that gives accurate information not on the boundary, but in the domain. The method is shown to drastically simplify the derivation of shape design sensitivity formulas for complex built-up structures.

Choi, Kyung K.↗
NASA NTRS2015

Higher-Order Methods for Compressible Turbulent Flows Using Entropy Variables

Turbulent flows have a large range of spatial and temporal scales which need to be resolved in order to obtain accurate predictions. Higher-order methods can provide greater efficiency for simulations requiring high spatial and temporal resolution, allowing for solutions with fewer degrees of freedom and lower computational cost than traditional second-order computational fluid dynamics (CFD) methods.1 Higher-order methods have been widely used for turbulent flows. However, the reduced numerical stabilization present in higher-order schemes implies that special care needs to be taken in the development of numerical methods to suppress nonlinear instabilities.2–6 In this work we present the development of a higher-order space-time discontinuous Galerkin method with a focus on the aspects of our numerical scheme required for ensuring nonlinear stability for turbulent simulations at high Reynolds numbers.

Diosady, Laslo T.↗
NASA NTRS1987

On direct numerical simulations of turbulent reacting flows

A description of the emerging field of direct numerical simulations of turbulent, chemically reacting flows is presented. The types of direct numerical simulations, physical issues related to implementing the simulations, as well as the various numerical methods used are described. Examples are presented of recent applications of direct numerical simulations to a variety of problems, displaying both the potential of the method and also some of its limitations. Finally, our view of the potential role of direct numerical simulations in future research on turbulent, chemically reacting flows is presented.

Jou, W.-H.↗
NASA NTRS2003

Numerical Simulation of One- and Two-Phase Flows in Propulsion Systems

Four subprojects were conducted using analytical methods, numerical simulation and experimental tests: (A) Shock wave mitigation by spike-shaped blunt bodies with application for the purpose of drag, lift and longitudinal momentum optimization. The main result in this subproject is: application of a single needle against a supersonic flow provides higher benefits for blunt body drag reduction and heat transfer to the body than the application of multiple needles. (B) Solid particles, liquid and air jet injection through the front of a blunt body against a supersonic flow. In this case, the research conducted and analysis of multiple previous investigations in this area have shown essential benefits and preferable application of solid particle injection. (C) Comparison of different methods of fuel injection into supersonic duct flows. Preliminary numerical simulations and theoretical analysis show promising results for Telescope-shaped inlet applications in SCRAMJET; and (D) Development of an acoustic source location method for different applications including propulsion systems.

Gilinsky, Mikhail↗
NASA NTRS1975

Recent progress in numerical techniques for flow simulation

Recent developments in the use of numerical methods for fluid flow simulation show an increasing tendency to use numerical operators that can, in various ways, be factored. Use of methods having this property often increases the accuracy and efficiency of computer codes. Tracing the factorization property provides a unity to the basic concepts involved in the development of cyclic reduction, predictor-corrector, splitting, fast Fourier transform, and pseudospectral methods.-

Lomax, H.↗
NASA NTRS1997

An Investigation of High-Order Shock-Capturing Methods for Computational Aeroacoustics

This project is motivated by the desire to develop numerical methods that will be useful in the study of compressible flows that exhibit aeroacoustic phenomena. Solutions to linear problems have been investigated through the development of a computer code based on the recent dispersion-relation-preserving (DRP) methodology. In regard to nonlinear problems, the class of essentially nonoscillatory (ENO) schemes have been considered as the primary candidates for solving aeroacoustic problems in which discontinuities are involved. Discontinuities in the solution itself (e.g. shocks) as well as in the geometry on which the problem is defined have been studied. Two-dimensional nonlinear problems were considered in order to determine if the one-dimensional results obtained in the first phase of this project were extendable to more realistic problems. Conclusions have been drawn in regard to the ability to numerically predict solutions of nonlinear problems with shocks to high-order accuracy.

Casper, Jay↗
NASA NTRS1992

Numerical study on mixing of sprayed liquid in an LNG storage tank

This paper presents a numerical method to simulate the mixing of heavier LNG sprayed on lighter layer. Numerical results for evolutions of flow field and density field are obtained in a rectangular computational domain which includes the vicinity of the liquid surface. At the surface boundary, uniform distributions of the fluid velocity and the density are assumed. Detail structure of flow caused by impingements of liquid drops are neglected. But, to trigger a realistic motion, a series of random numbers is employed. It is used as an initial distribution of the density near the surface. This method successfully gives a realistic simulation of the mixing process. Numerical results for mixing velocity shows good agreement with experimental data.

Uchida, Hiroyuki↗
NASA NTRS1972

Nuclear radiation analysis

A study program of radiation shielding against the deleterious effects of nuclear radiation on man and equipment is reported. The methods used to analyze the radiation environment from bremsstrahlung photons are discussed along with the methods employed by transport code users. The theory and numerical methods used to solve transport of neutrons and gammas are described, and the neutron and cosmic fluxes that would be present on the gamma-ray telescope were analyzed.

Knies, R. J.↗
NASA NTRS1972

Monte Carlo investigation of transient acoustic fields in partially or completely bounded medium

A simple repetitive calculation was used to investigate what happens to the field in terms of the signal paths of disturbances originating from the energy source. The computation allowed the field to be reconstructed as a function of space and time on a statistical basis. The suggested Monte Carlo method is in response to the need for a numerical method to supplement analytical methods of solution which are only valid when the boundaries have simple shapes, rather than for a medium that is bounded. For the analysis, a suitable model was created from which was developed an algorithm for the estimation of acoustic pressure variations in the region under investigation. The validity of the technique was demonstrated by analysis of simple physical models with the aid of a digital computer. The Monte Carlo method is applicable to a medium which is homogeneous and is enclosed by either rectangular or curved boundaries.

Thanedar, B. D.↗
NASA NTRS1984

A method for computing the leading-edge suction in a higher-order panel method

Experimental data show that the phenomenon of a separation induced leading edge vortex is influenced by the wing thickness and the shape of the leading edge. Both thickness and leading edge shape (rounded rather than point) delay the formation of a vortex. Existing computer programs used to predict the effect of a leading edge vortex do not include a procedure for determining whether or not a vortex actually exists. Studies under NASA Contract NAS1-15678 have shown that the vortex development can be predicted by using the relationship between the leading edge suction coefficient and the parabolic nose drag. The linear theory FLEXSTAB was used to calculate the leading edge suction coefficient. This report describes the development of a method for calculating leading edge suction using the capabilities of the higher order panel methods (exact boundary conditions). For a two dimensional case, numerical methods were developed using the double strength and downwash distribution along the chord. A Gaussian quadrature formula that directly incorporates the logarithmic singularity in the downwash distribution, at all panel edges, was found to be the best method.

Ehlers, F. E.↗
NASA NTRS1987

Analysis of non-premixed turbulent reacting flows

Studies of chemical reactions occurring in turbulent flows are important in the understanding of combustion and other applications. Current numerical methods are limited in their applications due to the numerical resolution required to completely capture all length scales, but, despite the fact that realistic combustion cannot be solved completely, numerical simulations can be used to give insight into the interaction between the processes of turbulence and chemical reaction. The objective was to investigate the effects of turbulent motion on the effects of chemical reaction to gain some insight on the interaction of turbulence, molecular diffusion, and chemical reaction to support modeling efforts. A direct turbulence simulation spectral code was modified to include the effects of chemical reaction and applied to an initial value problem of chemical reaction between non-premixed species. The influence of hydrodynamics on the instantaneous structure of the reaction was investigated. The local scalar dissipation rates and the local reaction rates were examined to determine the influence of vorticity or rate of strain on the reaction and the structure of the scalar field.

Leonard, Andy D.↗
NASA NTRS1989

A hybrid perturbation-Galerkin method for differential equations containing a parameter

A two-step hybrid perturbation-Galerkin method to solve a variety of differential equations which involve a parameter is presented and discussed. The method consists of: (1) the use of a perturbation method to determine the asymptotic expansion of the solution about one or more values of the parameter; and (2) the use of some of the perturbation coefficient functions as trial functions in the classical Bubnov-Galerkin method. This hybrid method has the potential of overcoming some of the drawbacks of the perturbation method and the Bubnov-Galerkin method when they are applied by themselves, while combining some of the good features of both. The proposed method is illustrated first with a simple linear two-point boundary value problem and is then applied to a nonlinear two-point boundary value problem in lubrication theory. The results obtained from the hybrid method are compared with approximate solutions obtained by purely numerical methods. Some general features of the method, as well as some special tips for its implementation, are discussed. A survey of some current research application areas is presented and its degree of applicability to broader problem areas is discussed.

Geer, James F.↗
NASA NTRS1984

Methods for analysis of cracks in three-dimensional solids

Various analytical and numerical methods used to evaluate the stress intensity factors for cracks in three-dimensional (3-D) solids are reviewed. Classical exact solutions and many of the approximate methods used in 3-D analyses of cracks are reviewed. The exact solutions for embedded elliptic cracks in infinite solids are discussed. The approximate methods reviewed are the finite element methods, the boundary integral equation (BIE) method, the mixed methods (superposition of analytical and finite element method, stress difference method, discretization-error method, alternating method, finite element-alternating method), and the line-spring model. The finite element method with singularity elements is the most widely used method. The BIE method only needs modeling of the surfaces of the solid and so is gaining popularity. The line-spring model appears to be the quickest way to obtain good estimates of the stress intensity factors. The finite element-alternating method appears to yield the most accurate solution at the minimum cost.

Raju, I. S.↗
NASA NTRS1984

A method for modeling finite-core vortices in wake-flow calculations

A numerical method for computing nonplanar vortex wakes represented by finite-core vortices is presented. The approach solves for the velocity on an Eulerian grid, using standard finite-difference techniques; the vortex wake is tracked by Lagrangian methods. In this method, the distribution of continuous vorticity in the wake is replaced by a group of discrete vortices. An axially symmetric distribution of vorticity about the center of each discrete vortex is used to represent the finite-core model. Two distributions of vorticity, or core models, are investigated: a finite distribution of vorticity represented by a third-order polynomial, and a continuous distribution of vorticity throughout the wake. The method provides for a vortex-core model that is insensitive to the mesh spacing. Results for a simplified case are presented. Computed results for the roll-up of a vortex wake generated by wings with different spanwise load distributions are presented; contour plots of the flow-field velocities are included; and comparisons are made of the computed flow-field velocities with experimentally measured velocities.

Stremel, P. M.↗
NASA NTRS1985

Computational methods for hypersonic viscous flow over finite ellipsoid-cones at incidence

A numerical method, which is simpler than others currently in use, is proposed for determining the full viscous flow over a finite body in hypersonic stream at high altitude. It treats the shock layer surrounding the blunt foebody and the near wake behind the base simultaneously by formulating the Navier-Stokes equations in conformal and azimuthal-angle coordinates. The computational domain is confined to the body wall, outflow surface and the bow shock, which is adjusted along the coordinate normal to the wall in the course of iterations. Because of the optimal grid and a well developed alternating direction implicit factorization technique for the governing equations, reasonably accurate results can be obtained on a 30 by 36 by 6 grid with 400 time-marching iterations. Results for body shapes belonging to the ellipsoid-cone family are compared with the experimental data for the Apollo command module and the Viking aeroshell. Validation of the method based on self-consistency is also discussed.

Li, C. P.↗
NASA NTRS1994

Theoretical Investigation of Wave-Vector-Dependent Analytical and Numerical Formulations of the Interband Impact-Ionization Transition Rate for Electron in Bulk Silicon and GaAs

The electron interband impact-ionization rate for both silicon and gallium arsenide is calculated using an ensemble Monte Carlo simulation with the expressed purpose of comparing different formulations of the interband ionization transition rate. Specifically, three different treatments of the transition rate are examined: the traditional Keldysh formula, a new k-dependent analytical formulation first derived by W. Quade, E Scholl, and M. Rudan, and a more exact, numerical method of Y. Wang and K. F. Brennan. Although the completely numerical formulation contains no adjustable parameters and as such provides a very reliable result, it is highly computationally intensive. Alternatively, the Keldysh formula, although inherently simple and computationally efficient, fails to include the k dependence as well as the details of the energy band structure. The k-dependent analytical formulation of Quade and co-workers overcomes the limitations of both of these models but at the expense of some new parameterization. It is found that the k-dependent analytical method of Quade and co-workers produces very similar results to those obtained with (he completely numerical model for some quantities. Specifically, both models predict that the effective threshold for impact ionization in GaAs and silicon is quite soft, that the majority of ionization events originate from the second conduction band in both materials, and that the transition rate is k dependent. Therefore, it is concluded that the k-dependent analytical model can qualitatively reproduce results similar to those obtained with the numerical model yet with far greater computational efficiency. Nevertheless, there exist some important drawbacks to the k-dependent analytical model of Quade and co-workers: These are that it does not accurately reproduce the quantum yield data for bulk silicon, it requires determination of a new parameter, related physically to (he overlap integrals of the Bloch state which can only be adjusted by comparison to experiment, and fails to account for any wave-vector dependence of the overlap integrals. As such [he transition rate may be overestimated at those points for which "near vertical," small change in k, transitions occur.

Kolnik, Jan↗
NASA NTRS1991

A Runge-Kutta discontinuous finite element method for high speed flows

A Runge-Kutta discontinuous finite element method is developed for hyperbolic systems of conservation laws in two space variables. The discontinuous Galerkin spatial approximation to the conservation laws results in a system of ordinary differential equations which are marched in time using Runge-Kutta methods. Numerical results for the two-dimensional Burger's equation show that the method is (p+1)-order accurate in time and space, where p is the degree of the polynomial approximation of the solution within an element and is capable of capturing shocks over a single element without oscillations. Results for this problem also show that the accuracy of the solution in smooth regions is unaffected by the local projection and that the accuracy in smooth regions increases as p increases. Numerical results for the Euler equations show that the method captures shocks without oscillations and with higher resolution than a first-order scheme.

Bey, Kim S.↗