SEARCH · Search NASA

Results for “numerics”

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 505 records · Page 28

All sources NASA NTRS DOE OSTI

NASA NTRS1977

A numerical investigation of the finite element method in compressible primitive variable Navier-Stokes flow

The results of a comprehensive numerical investigation of the basic capabilities of the finite element method (FEM) for numerical solution of compressible flow problems governed by the two-dimensional and axis-symmetric Navier-Stokes equations in primitive variables are presented. The strong and weak points of the method as a tool for computational fluid dynamics are considered. The relation of the linear element finite element method to finite difference methods (FDM) is explored. The calculation of free shear layer and separated flows over aircraft boattail afterbodies with plume simulators indicate the strongest assets of the method are its capabilities for reliable and accurate calculation employing variable grids which readily approximate complex geometry and capably adapt to the presence of diverse regions of large solution gradients without the necessity of domain transformation.

Cook, C. H.↗

NASA NTRS1977

Ray trajectories in a torus: An application of MACSYMA to a complex numerical computation

The study of ray trajectories of plasma waves in a torodial geometry using MACSYMA is an example of how symbolic, numerical, and graphical facilities can be used in concert to accomplish a complex computational goal. Computational features of this study which are of particular significance include: the derivation of code (i.e. writing functions to generate program fragments), the use of array functions to simplify the specification of a numerical iteration scheme, and the graphical presentation of the results. Mathematically, this study originates in the solution of a linear inhomogeneous partial differential equation in 3 dimensions by the method of characteristics. While it is possible to describe this equation compactly by using vector notation, and by specifying the spatial variation of the coefficients in terms of intermediate parameters, the transformation of the equation into a form amenable to solution is very tedious. A MACSYMA program is presented for obtaining description of the rf field structure excited by a waveguide located at the edge of a toroidal plasma confinement device.

Kulp, J. L.↗

NASA NTRS1976

Sensitivity of inelastic response to numerical integration of strain energy

The exact solution to the quasi-static, inelastic response of a cantilever beam of rectangular cross section subjected to a bending moment at the tip is obtained. The material of the beam is assumed to be linearly elastic-linearly strain-hardening. This solution is then compared with three different numerical solutions of the same problem obtained by minimizing the total potential energy using Gaussian quadratures of two different orders and a Newton-Cotes scheme for integrating the strain energy of deformation. Significant differences between the exact dissipative strain energy and its numerical counterpart are emphasized. The consequence of this on the nonlinear transient responses of a beam with solid cross section and that of a thin-walled beam on elastic supports under impulsive loads are examined.

Kamat, M. P.↗

NASA NTRS1976

A numerical comparison of discrete Kalman filtering algorithms - An orbit determination case study

An improved Kalman filter algorithm based on a modified Givens matrix triangularization technique is proposed for solving a nonstationary discrete-time linear filtering problem. The proposed U-D covariance factorization filter uses orthogonal transformation technique; measurement and time updating of the U-D factors involve separate application of Gentleman's fast square-root-free Givens rotations. Numerical stability and accuracy of the algorithm are compared with those of the conventional and stabilized Kalman filters and the Potter-Schmidt square-root filter, by applying these techniques to a realistic planetary navigation problem (orbit determination for the Saturn approach phase of the Mariner Jupiter-Saturn Mission, 1977). The new algorithm is shown to combine the numerical precision of square root filtering with the efficiency of the original Kalman algorithm.

Thornton, C. L.↗

NASA NTRS1977

Reynolds number effects on shock-wave turbulent boundary-layer interactions - A comparison of numerical and experimental results

An experiment is described that tests and guides computations of a shock-wave turbulent boundary-layer interaction flow over a 20-deg compression corner at Mach 2.85. Numerical solutions of the time-averaged Navier-Stokes equations for the entire flow field, employing various turbulence models, are compared with the data. Each model is critically evaluated by comparisons with the details of the experimental data. Experimental results for the extent of upstream pressure influence and separation location are compared with numerical predictions for a wide range of Reynolds numbers and shock-wave strengths.

Horstman, C. C.↗

NASA NTRS1977

Numerical solution of axisymmetric boattail flow fields with plume simulators

Turbulent separating flows over axisymmetric afterbody-boattail configurations with solid sting plume simulators are computed with a time-dependent finite-difference method to solve the compressible Navier-Stokes equations. The Reynolds stress terms are replaced with a two-layer eddy viscosity model including a relaxation formula to model the nonequilibrium effects of the separated flow. The mesh is alined with the boattail body through an analytic transformation which accommodates a wide variety of boattail geometries. Numerical results for a series of boattail geometries over a wide range of Reynolds number (140,000 to 140 million) are presented and, when possible, compared with experimental data or independent numerical results.

Holst, T. L.↗

NASA NTRS1976

A numerical study of flow phenomena within a vortex sink rate sensor

A numerical investigation of a viscous incompressible fluid with constant properties within a vortex sink rate sensor is undertaken. A computer program was developed, for obtaining the numerical solution of the Navier-Stokes equations, assuming laminar flow, having general prescribed inlet conditions and axisymmetric boundary conditions. Computational results for various viscous flows and assorted boundary conditions have been obtained. In the vortex chamber, the radial momentum is conserved, whereas the axial momentum is negligible. In the sink region, the radial and tangential velocities, as well as the tangential vorticities, have significantly high values. In the sink tube, the radial momentum is negligible. The radial momentum within the vortex chamber is converted to the axial momentum in the sink tube. A change in the flow rate had an appreciable effect within the sensor and particularly near the sink tube entrance. As the flow rate was increased, both the tangential velocity and tangential vorticity increased rapidly. At the higher flow rates, vortices were produced at the corner of the entrance section of the sink tube.

Patel, D. K.↗

NASA NTRS1976

Dynamic analysis of a class of mixed lumped-distributed parameter systems via numerical techniques

The modeling via transfer matrices of mixed lumped-distributed parameter systems with feedback control is discussed and novel methods of obtaining information on the controlled system's dynamic response via numerical methods is presented. These methods utilize numerical searches in the complex plane for the roots of the transcendental characteristic equation to determine the closed loop system eigenvalues. The fast Fourier transform (FFT) algorithm is utilized but in the inverse fashion (from frequency domain to time domain). The information obtained is identical to that contained in a model in modal coordinates, and can be used to construct such a model. An interactive computer program to do the required calculation is described, as well as its application to a manipulator arm design problem.

Book, W. J.↗

NASA NTRS1976

Numerical optimization of gasdynamic lasers

A sophisticated numerical optimization procedure has been coupled to a detailed theoretical model of the carbon dioxide gasdynamic laser (GDL) for the purpose of systematically evaluating the performance of these systems. The integrated model was used to predict the performance of a wide variety of combustion-driven gasdynamic lasers and was then used to maximize the theoretical available specific energy of these lasers. Up to six design variables and three nonlinear constraints were used in these optimizations, although the method is not limited to this number. Results are presented to show the performance of typical lasers before and after optimization. Improvements approaching a factor of two in the available specific energy are found, suggesting that numerical optimization techniques are valuable in improving the design of gasdynamic lasers.

Saunders, R. C., III↗

NASA NTRS1976

Numerical problems in the advection of pollutants

Prediction of pollutant transport involves solution of the advection equation numerically. Truncation errors arising from integrating the advection equation by using finite difference schemes are discussed. Numerical errors manifest themselves in artificial dissipation and dispersion of pollutants. Criteria for measuring the dissipation and dispersion are presented, and minimization of those errors is discussed. An alternative procedure by integration in the wavenumber space is examined. Examples illustrating quantitative discrepancy in actual forecasting are also shown.

Liu, C. Y.↗

NASA NTRS1977

A critical study of higher-order numerical methods for solving the boundary-layer equations

A fourth-order box method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method is the natural extension of the second-order Keller Box Scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary-layer equations. The efficiency of the present method is compared with other two-point and three-point higher-order methods; namely, the Keller Box Scheme with Richardson extrapolation, the method of deferred corrections, and the three-point spline methods. For equivalent accuracy, numerical results show the present method to be more efficient than the other higher-order methods for both laminar and turbulent flows.

Wornom, S. F.↗

NASA NTRS1977

Application of numerical optimization to the design of supercritical airfoils without drag-creep

Recent applications of numerical optimization to the design of advanced airfoils for transonic aircraft have shown that low-drag sections can be developed for a given design Mach number without an accompanying drag increase at lower Mach numbers. This is achieved by imposing a constraint on the drag coefficient at an off-design Mach number while the drag at the design Mach number is the objective function. Such a procedure doubles the computation time over that for single design-point problems, but the final result is worth the increased cost of computation. The ability to treat such multiple design-point problems by numerical optimization has been enhanced by the development of improved airfoil shape functions. Such functions permit a considerable increase in the range of profiles attainable during the optimization process.

Hicks, R. M.↗

NASA NTRS1976

Comparison of numeric keyboard and CRT line-labeled buttons for information access

Test were conducted to determine whether differences in speed and accuracy are experienced when using either line-labeled index buttons or a numeric keyboard for page selection in airborne CRT-display area navigation systems. The experiment was conducted with six airline pilots, each flying the same two simulated RNAV routes. Three pilot subjects used line-labeled buttons adjacent to the CRT screen, while three used a numeric keyboard for page access. The hypothesis of no differences in response times between the two modes of access could not be rejected.

Williams, D.↗

NASA NTRS1977

Numerical solutions for viscous and potential flow about arbitrary two-dimensional bodies using body-fitted coordinate systems

A procedure for numerical solution of the time-dependent, incompressible Navier-Stokes equations for the flow about arbitrarily shaped two-dimensional bodies is given. This solution is based on a technique of automatic numerical generation of a curvilinear coordinate system having a coordinate line coincident with the body contour regardless of its shape. The implicit solution utilizes the vorticity-stream function formulation with a false-position iterative adjustment of the surface vorticity in satisfaction of the no-slip boundary condition. Excellent agreement with the Blasius boundary layer solution is obtained for a semi-infinite flat plate. Results are presented for Reynolds numbers up to 2000 for several airfoils and a cambered rock.

Thames, F. C.↗

NASA NTRS1977

Numerical simulation of spacecraft charging phenomena

A numerical simulation program is being constructed having the following features: (1) infinite circular cylindrical geometry with angle-dependence, (2) inclusion of incident particles, photoelectrons, secondary electrons, backscattered electrons, any gun emissions, and any internal current pathways including surface conductive layers, (3) quasistatic time-dependent iteration, in which sheath potential changes during particle transit times are ignored, (4) use of approximate, locally-dependent space charge density expressions in solving Poisson's equation for sheath potentials, with use of numerical orbit-following to determine surface currents, and (5) incident particle velocity distributions isotropic or beam-like, or some superposition of these. Rationales for each of these features are discussed.

Laframboise, J. G.↗

NASA NTRS1977

Stability of numerical integration techniques for transient rotor dynamics

A finite element model of a rotor bearing system was analyzed to determine the stability limits of the forward, backward, and centered Euler; Runge-Kutta; Milne; and Adams numerical integration techniques. The analysis concludes that the highest frequency mode determines the maximum time step for a stable solution. Thus, the number of mass elements should be minimized. Increasing the damping can sometimes cause numerical instability. For a uniform shaft, with 10 mass elements, operating at approximately the first critical speed, the maximum time step for the Runge-Kutta, Milne, and Adams methods is that which corresponds to approximately 1 degree of shaft movement. This is independent of rotor dimensions.

Kascak, A. F.↗

NASA NTRS1977

Numerical methods for inviscid/viscous flows

A study of numerical schemes for solving viscous fluid flow problems with sizable regions of predominantly inviscid flow is presented. Difficulties associated with the familiar central difference approach for such problems were analyzed and alternative finite difference approaches employing windward concepts are presented. In addition, difference relations based on exponential operators were developed. All such schemes were demonstrated and evaluated through application to the case of Falkner Skan flow with blowing - a problem in which a sizable region of predominantly inviscid flow develops near the injection surface that traditionally causes numerical difficulty.

El-Mistikawy, T. M.↗

NASA NTRS1977

Numerical determination of the fundamental eigenvalue for the Laplace operator on a spherical domain

Methods for obtaining approximate solutions for the fundamental eigenvalue of the Laplace-Beltrami operator (i.e., the membrane eignevalue problem for the vibration equation) on the unit spherical surface are developed. Two types of spherical surface domains are considered: the interior of a spherical triangle, and the exterior of a great circle arc extending for less than pi radians (a spherical surface with a slit). In both cases, zero boundary conditions are imposed. In order to solve the resulting second-order elliptic partial differential equations in two independent variables, a finite difference approximation is employed. The fundamental eigenvalue is approximated by iteration utilizing the power method and point successive overrelaxation. Some numerical results are given and compared, in certain special cases, with analytical solutions to the eigenvalue problem. The significance of the numerical eigenvalue results is discussed in terms of the singularities in the solution of three-dimensional boundary-value problems near a polyhedral corner of the domain.

Walden, H.↗