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 361 records · Page 20

Computation of optimal low- and medium-thrust orbit transfers

This paper presents the formulation of the optimal low- and medium-thrust orbit transfer control problem, numerical methods for solution, and numerical solutions of the problem. The problem formulation is for final mass maximization and allows for second-harmonic oblateness, atmospheric drag, and 3D noncoplanar nonaligned elliptic terminal orbits. We set up examples to demonstrate the ability of two indirect methods to solve the resulting two point boundary value problems (TPBVP). The methods demonstrated are the multiple point shooting method as formulated in Oberle's (1987) subroutine BOUNDSCO, and the minimizing boundary-condition method (MBCM). We find that although both methods can converge solutions, there are tradeoffs to using either method. We present numerical solutions of planar transfers in which both the initial orbit exit and final orbit entry points have been optimized. These solutions include two- and three-burn transfers. The methods used show an ability to handle thrust down to at least T/W(o) = O(10 exp -3). They also show similar convergence abilities with or without the oblateness and drag terms. We discuss the issue of maximizing with respect to the final time and provide evidence that implies a local optimum at a maximum final time for a given number of burns.

Chuang, C.-H.↗

Numerical analysis of some supersonic viscous flows related to inlet and nozzle systems

A numerical method originally developed for three-dimensional supersonic inlet flow calculations is extended and applied to the study of two-dimensional and three-dimensional flows associated with arbitrary propulsion systems. The method is based on the forward spatial marching solution of a reduced form of the three-dimensional steady Navier-Stokes equations in which streamwise pressure gradients are retained in both the subsonic and supersonic regions. The present paper briefly describes the analysis and then shows three applications. In the first application, a wall transpiration study has been performed for the two-dimensional shock wave/turbulent boundary layer interaction flow field with application to an inlet configuration. The second application treats a rectangular high speed inlet with a swept sideplate including the effects of sideplate spillage. Finally, the method is utilized to analyze the interaction of an under-expanded supersonic jet with an ambient flow. Computed results are examined and compared with available experimental measurements. It is demonstrated that the present numerical method is capable of numerically simulating complex two- and three-dimensional flows relevant to hypersonic propulsion systems in a manner which both shows good agreement with data when such data is available, and which shows the complex flow features in the absence of data.

Kim, Y.-N.↗

Summary of research in applied mathematics, numerical analysis, and computer sciences

The major categories of current ICASE research programs addressed include: numerical methods, with particular emphasis on the development and analysis of basic numerical algorithms; control and parameter identification problems, with emphasis on effective numerical methods; computational problems in engineering and physical sciences, particularly fluid dynamics, acoustics, and structural analysis; and computer systems and software, especially vector and parallel computers.

Source record↗

Viscous Transonic Airfoil Workshop compendium of results

Results from the Viscous Transonic Airfoil Workshop held at the AIAA 25th Aerospace Sciences Meeting at Reno, NV in January 1987, are compared with each other and with experimental data. Test cases used in this workshop include attached and separated transonic flows for three different airfoils: the NACA 0012 airfoil, the RAE 2822 airfoil, and the Jones airfoil. A total of 23 sets of numerical results from 15 different author groups are included. The numerical methods used vary widely and include: 16 Navier-Stokes methods, 2 Euler/boundary-layer methods, and 5 full-potential/boundary-layer methods. The results indicate a high degree of sophistication among the numerical methods with generally good agreement between the various computed and experimental results for attached or moderately-separated cases. The agreement for cases with larger separation is only fair and suggests additional work is required in this area.

Holst, Terry L.↗

Computational fluid dynamics drag prediction: Results from the Viscous Transonic Airfoil Workshop

Results from the Viscous Transonic Airfoil Workshop are compared with each other and with experimental data. Test cases used include attached and separated transonic flows for the NACA 0012 airfoil. A total of 23 sets of numerical results from 15 different author groups are included. The numerical method used vary widely and include: 16 Navier-Stokes methods, 2 Euler boundary layer methods, and 5 potential boundary layer methods. The results indicate a high degree of sophistication among the numerical methods with generally good agreement between the various computed and experimental results for attached or moderately separated cases. The agreement for cases with larger separation is only fair and suggests additional work is required in this area.

Holst, Terry L.↗

Computational fluid dynamics drag prediction: Results from the Viscous Transonic Airfoil Workshop

Results from the Viscous Transonic Airfoil Workshop held in January 1987, are compared with each other and with experimental data. Test cases used include attached and separated transonic flows for the NACA 0012 airfoil. A total of 23 sets of numerical results from 15 different author groups are included. The numerical methods used vary widely and include: 16 Navier-Stokes methods, 2 Euler/boundary-layer methods, and 5 potential/boundary-layer methods. The results indicate a high degree of sophistication among the numerical methods with generally good agreement between the various computed and experimental results for attached or moderately separated cases. The agreement for cases with larger separation is only fair and suggests additional work is required in this area.

Holst, Terry L.↗

Cumulative reports and publications through December 31, 1989

A complete list of reports from the Institute for Computer Applications in Science and Engineering (ICASE) is presented. The major categories of the current ICASE research program are: numerical methods, with particular emphasis on the development and analysis of basic numerical algorithms; control and parameter identification problems, with emphasis on effectual numerical methods; computational problems in engineering and the physical sciences, particularly fluid dynamics, acoustics, structural analysis, and chemistry; computer systems and software, especially vector and parallel computers, microcomputers, and data management. Since ICASE reports are intended to be preprints of articles that will appear in journals or conference proceedings, the published reference is included when it is available.

Source record↗

ICASE semiannual report, April 1 - September 30, 1989

The Institute conducts unclassified basic research in applied mathematics, numerical analysis, and computer science in order to extend and improve problem-solving capabilities in science and engineering, particularly in aeronautics and space. The major categories of the current Institute for Computer Applications in Science and Engineering (ICASE) research program are: (1) numerical methods, with particular emphasis on the development and analysis of basic numerical algorithms; (2) control and parameter identification problems, with emphasis on effective numerical methods; (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. ICASE reports are considered to be primarily preprints of manuscripts that have been submitted to appropriate research journals or that are to appear in conference proceedings.

Source record↗

Semiannual report, 1 April - 30 September 1991

The major categories of the current Institute for Computer Applications in Science and Engineering (ICASE) research program are: (1) numerical methods, with particular emphasis on the development and analysis of basic numerical algorithms; (2) control and parameter identification problems, with emphasis on effective numerical methods; (3) computational problems in engineering and the physical sciences, particularly fluid dynamics, acoustics, and structural analysis; and (4) computer systems and software for parallel computers. Research in these areas is discussed.

Source record↗

Reynolds stress closure in jet flows using wave models

Numerical methods were developed that will form the computational part of the turbulence closure scheme. A wave model was developed for the two-dimensional shear layer. This configuration is being used as a test case for the closure schemes. Various numerical schemes were examined to give efficient solutions of the Rayleigh equation for this geometry. These include both spectral and finite difference methods. Secondly, numerical methods are under development to solve the non-separable Rayleigh equation. This solution is required for the closure scheme in more complex geometries. A model problem was used to assist in the algorithm development. Two-dimensional spectral methods and a hybrid spectral/finite difference technique were developed. An analytic solution of the Rayleigh equation for a basic elliptic flow was obtained. This will be used to verify the stability codes developed for arbitrary geometries. Other numerical methods for solving the Rayleigh equation based on the boundary element technique were also examined. These solutions are forming the basis of a model for the shock structure in jets of arbitrary geometry.

Morris, P. J.↗

Numerical Prediction Methods (Reynolds-Averaged Navier-Stokes Simulations of Transonic Separated Flows)

During the past five years, numerous pioneering archival publications have appeared that have presented computer solutions of the mass-weighted, time-averaged Navier-Stokes equations for transonic problems pertinent to the aircraft industry. These solutions have been pathfinders of developments that could evolve into a major new technological capability, namely the computational Navier-Stokes technology, for the aircraft industry. So far these simulations have demonstrated that computational techniques, and computer capabilities have advanced to the point where it is possible to solve forms of the Navier-Stokes equations for transonic research problems. At present there are two major shortcomings of the technology: limited computer speed and memory, and difficulties in turbulence modelling and in computation of complex three-dimensional geometries. These limitations and difficulties are the pacing items of the continuing developments, although the one item that will most likely turn out to be the most crucial to the progress of this technology is turbulence modelling. The objective of this presentation is to discuss the state of the art of this technology and suggest possible future areas of research. We now discuss some of the flow conditions for which the Navier-Stokes equations appear to be required. On an airfoil there are four different types of interaction of a shock wave with a boundary layer: (1) shock-boundary-layer interaction with no separation, (2) shock-induced turbulent separation with immediate reattachment (we refer to this as a shock-induced separation bubble), (3) shock-induced turbulent separation without reattachment, and (4) shock-induced separation bubble with trailing edge separation.

Mehta, Unmeel↗

Numerical solutions of differential equations

Various numerical methods for solving differential equations were analyzed and refined in an effort to develop a method which was adaptable to a large class of problems. The prime capabilities of the method included accuracy, numerical stability, and economic use of computer time. In multistep processes the corrector was changed at each step.

Wesson, J. R.↗

Some methods for determining the profile functions of inhomogeneous media

The problem of determining the electric permittivity of an inhomogeneous dielectric medium is considered. Two analytical methods and one numerical method for solving this problem are presented. All these methods utilize the reflection coefficient measured as a function of the angular spectrum variable. This allows measurements at a fixed frequency, and these methods, therefore, can be used when the medium being studied is dispensive. A few examples are provided and some comments are made regarding the stability of solutions obtained using the various methods.

Mittra, R.↗

Incorporation of a progressive failure analysis method in the CSM testbed software system

Analysis of the postbuckling behavior of composite shell structures pose many difficult and challenging problems in the field of structural mechanics. Current analysis methods perform well for most cases in predicting the postbuckling response of undamaged components. To predict component behavior accurately at higher load levels, the analysis must include the effects of local material failures. The CSM testbed software system is a highly modular structural analysis system currently under development at Langley Research Center. One of the primary goals of the CSM testbed is to provide a software environment for the development of advanced structural analysis methods and modern numerical methods which will exploit advanced computer architecture such as parallel-vector processors. Development of a progressive failure analysis method consists of the design and implementation of a processor which will perform the ply-level progressive failure analysis and the development of a geometrically nonlinear analysis procedure which incorporates the progressive failure processor. Regarding the development of the progressive failure processor, two components are required: failure criteria and a degradation model. For the initial implementation, the failure criteria of Hashin will be used. For a matrix failure which typically indicates the development of transverse matrix cracks, the ply properties will be degraded. Work to date includes the design of the progressive failure analysis processor and initial plans for the controlling geometrically nonlinear analysis procedure. The implementation of the progressive failure analysis has begun. Access to the model database and the Hashin failure criteria are completed. Work is in progress on the input/output operations for the processor related data and the finite element model updating procedures. In total the progressive failure processor is approximately one-third complete.

Arenburg, Robert T.↗

Magnetic Field Suppression of Flow in Semiconductor Melt

One of the most promising approaches for the reduction of convection during the crystal growth of conductive melts (semiconductor crystals) is the application of magnetic fields. Current technology allows the experimentation with very intense static fields (up to 80 KGauss) for which nearly convection free results are expected from simple scaling analysis in stabilized systems (vertical Bridgman method with axial magnetic field). However, controversial experimental results were obtained. The computational methods are, therefore, a fundamental tool in the understanding of the phenomena accounting during the solidification of semiconductor materials. Moreover, effects like the bending of the isomagnetic lines, different aspect ratios and misalignments between the direction of the gravity and magnetic field vectors can not be analyzed with analytical methods. The earliest numerical results showed controversial conclusions and are not able to explain the experimental results. Although the generated flows are extremely low, the computational task is a complicated because of the thin boundary layers. That is one of the reasons for the discrepancy in the results that numerical studies reported. Modeling of these magnetically damped crystal growth experiments requires advanced numerical methods. We used, for comparison, three different approaches to obtain the solution of the problem of thermal convection flows: (1) Spectral method in spectral superelement implementation, (2) Finite element method with regularization for boundary layers, (3) Multiquadric method, a novel method with global radial basis functions, that is proven to have exponential convergence. The results obtained by these three methods are presented for a wide region of Rayleigh and Hartman numbers. Comparison and discussion of accuracy, efficiency, reliability and agreement with experimental results will be presented as well.

Fedoseyev, A. I.↗

Electromagnetic Scattering by Discrete Random Media Illuminated by a Gaussian Beam II: Solution of the Radiative Transfer Equation

In this paper, we present numerical methods for solving the phenomenological scalar radiative transfer equation for a discrete random medium illuminated by a Gaussian beam. These rely on the Fourier transform method for the horizontal variables and the discrete ordinate method with matrix exponential for solving the underlying one-dimensional radiative transfer equation in the wavenumber domain. The problem of a Gaussian beam at oblique and normal incidence, as well as, the searchlight problem are treated. A complete description of the methods and the numerical algorithms is provided.

Gaussian beam↗