SEARCH · Engineering Papers

Results for “Numerical approximations”

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 1,153 records · Page 64

All sources NASA NTRS DOE OSTI

NASA NTRS1980

Solution of Euler's Equations of Motion and Eulerian Angles for near symmetric rigid bodies subject to constant moments

Analytic expressions are found for Euler's Equations of Motion and for the Eulerian Angles for both symmetric and near symmetric rigid bodies under the influence of arbitrary constant body-fixed torques. These solutions provide the body-fixed angular velocities and the attitude of the body, respectively, as functions of time. They are of special interest in applications to spinning spacecraft (such as the Galileo Spacecraft to be launched in 1984) because they include the effect of time-varying spin rate. Thus they can be applied to spin-up and spin-down maneuvers as well as to error analysis for thruster misalignments. The solutions are given for arbitrary initial conditions in terms of Fresnel, Sine and Cosine Integrals. Numerical integration of the governing differential equations has verified that the approximate analytic solutions are very accurate in many physical situations of interest.

Longuski, J. M.↗

NASA NTRS1982

Comparison of starting values for iterative solutions to a universal Kepler's equation

General starting values for the iterative numerical solution of a universal Kepler's equation for position in a conic orbit at a specified time are investigated. Three starting values based on recent refinements of previously obtained bounds on the solution are derived and tested numerically. Of these, a simple starting value based on a cubic approximation to Kepler's equation provides the most rapid convergence using both first and second order Newton algorithms. The performance of the starting values are compared with similar studies which used the restricted case of elliptical orbits with the initial epoch at periapse.

Bergam, M. J.↗

NASA NTRS1982

Solar particle propagation and modulation theory

Developments related to solar flare particle propagation theory are considered, taking into account the transport equation currently employed to describe the propagation of solar flare particles. This equation does not take into consideration convection with the solar wind, drifts, and changes in particle energy. The equation was solved numerically to generate typical solar flare particle time profiles at 1 AU. The results, however, were not satisfactory. It was found that the numerical method used was apparently in error. New advances which provide approximate analytical solutions to the equation are reported. Questions of galactic cosmic ray propagation in interplanetary space are also discussed.

Lee, M. A.↗

NASA NTRS1983

Review of pseudoinverse control for use with kinematically redundant manipulators

Kinematically redundant manipulators have a number of potential advantages over current manipulator designs. For this type of arm, velocity control through pseudoinverses has been suggested. Questions associated with pseudoinverse control are examined in detail and show that in some cases this control leads to undesired arm configurations. A method for distributing joint angles of a redundant arm in a good approximation to a true minimax criterion is described. In addition several numerical considerations are discussed.

Klein, C. A.↗

NASA NTRS1983

Steady flow approximations to the helium r-process

A steady flow approximation to the r-process is presented and used for numerical experiments with physical quantities to determine the sensitivity of the process to variations in those quantities. The effect of neutron capture cross sections along the capture path and of recently available improved beta decay rates on the r-process are discussed. The peaks in the observed r-process yield curve near mass numbers 80 and 130 are roughly characterized by a neutron number density of 10 to the 20th per cu/cm; the mean beta decay rates are about 10/s, and the freezing time is comparable to or less than 0.1 s. The peak near mass number 195 is roughly characterized by a neutron number density of 10 to the 21st/cm, the mean beta decay rates are about 100/s, and the freezing time is comparable to or less than 0.01 s. The flow path of the steady state r-process is sensitively dependent upon the neutron capture cross sections in the flow network and on the values of the beta decay rates.

Cameron, A. G. W.↗

NASA NTRS1985

A three-dimensional viscous flow analysis for the helicopter tip vortex generation problem

The tip vortex flow field occurring in the vicinity of the tip region of a a helicopter rotor blade is a very complicated three-dimensional, viscous flow phenomenon. The details of the flow in the tip region can have a major effect in determining the generated rotor noise and can significantly affect the performance and dynamic loading of the rotor blade. The three-dimensional viscous subsonic tip vortex generation processes is investigated by a numerical procedure which allows spatial forward-marching integration, utilizing flow approximations from the velocity-decomposition approach of Briley and McDonald. The approach has been applied to compute the laminar and turbulent tip vortex flows for a constant thickness slab airfoil with a square tip, a constant thickness slab airfoil with a half round tip and a NACA 0012 airfoil with a half round tip. The basic mechanism of the tip vortex generation process as well as the prediction of vortex appearance, strength and secondary flow shown by the calculations are in qualitative agreement with experimental results.

Lin, S. J.↗

NASA NTRS1985

Recent improvements in efficiency, accuracy, and convergence for implicit approximate factorization algorithms

In 1977 and 1978, general purpose centrally space differenced implicit finite difference codes in two and three dimensions have been introduced. These codes, now called ARC2D and ARC3D, can run either in inviscid or viscous mode for steady or unsteady flow. Since the introduction of the ARC2D and ARC3D codes, overall computational efficiency could be improved by making use of a number of algorithmic changes. These changes are related to the use of a spatially varying time step, the use of a sequence of mesh refinements to establish approximate solutions, implementation of various ways to reduce inversion work, improved numerical dissipation terms, and more implicit treatment of terms. The present investigation has the objective to describe the considered improvements and to quantify advantages and disadvantages. It is found that using established and simple procedures, a computer code can be maintained which is competitive with specialized codes.

Pulliam, T. H.↗

NASA NTRS1985

Refined numerical solution of the transonic flow past a wedge

A numerical procedure combining the ideas of solving a modified difference equation and of adaptive mesh refinement is introduced. The numerical solution on a fixed grid is improved by using better approximations of the truncation error computed from local subdomain grid refinements. This technique is used to obtain refined solutions of steady, inviscid, transonic flow past a wedge. The effects of truncation error on the pressure distribution, wave drag, sonic line, and shock position are investigated. By comparing the pressure drag on the wedge and wave drag due to the shocks, a supersonic-to-supersonic shock originating from the wedge shoulder is confirmed.

Liang, S.-M.↗

NASA NTRS1985

Nonlinear acoustic wave propagation in atmosphere

A model problem that simulates an atmospheric acoustic wave propagation situation that is nonlinear is considered. The model is derived from the basic Euler equations for the atmospheric flow and from the regular perturbations for the acoustic part. The nonlinear effects are studied by obtaining two successive linear problems in which the second one involves the solution of the first problem. Well posedness of these problems is discussed and approximations of the radiation boundary conditions that can be used in numerical simulations are presented.

Hariharan, S. I.↗

NASA NTRS1986

Nonlinear acoustic wave propagation in atmosphere

In this paper a model problem is considered that simulates an atmospheric acoustic wave propagation situation that is nonlinear. The model is derived from the basic Euler equations for the atmospheric flow and from the regular perturbations for the acoustic part. The nonlinear effects are studied by obtaining two successive linear problems in which the second one involves the solution of the first problem. Well-posedness of these problems is discussed and approximations of the radiation boundary conditions that can be used in numerical simulations are presented.

Hariharan, S. I.↗

NASA NTRS1986

Estimation of discontinuous coefficients and boundary parameters for hyperbolic systems

The problem of estimating discontinuous coefficients, including locations of discontinuities, that occur in second order hyperbolic systems typical of those arising in I-D surface seismic problems is discussed. In addition, the problem of identifying unknown parameters that appear in boundary conditions for the system is treated. A spline-based approximation theory is presented, together with related convergence findings and representative numerical examples.

Lamm, P. K.↗

NASA NTRS1986

A three-dimensional viscous flow analysis for the helicopter tip vortex generation process about square and round tipped blades

The tip vortex flow field occurring in the vicinity of the tip region of a helicopter rotor blade is a very complicated three-dimensional, viscous flow phenomenon. The details of the flow in the tip region can have a major effect in determining the generated rotor noise and can significantly effect the performance and dynamic loading of the rotor blade. The three-dimensional viscous subsonic tip vortex generation processes is investigated by a numerical procedure which allows spatial forward-marching integration, utilizing flow approximations from the velocity-decomposition approach of Briley and McDonald. The approach has been applied to compute the laminar and turbulent tip vortex flows for a constant thickness slab airfoil with a square tip, a constant thickness slab airfoil with a half round tip and a NACA 0012 airfoil with a half round tip. The basic mechanism of the tip vortex generation process as well as the prediction of vortex appearance, strength and secondary flow shown by the calculations are in qualitative agreement with experimental results.

Lin, S.-J.↗

NASA NTRS1987

Calculating Sonic-Boom Propagation

Nonlinear effects included, enabling more-realistic modeling. Modified Method of Characteristics Sonic Boom Extrapolation Program (MMOC) is computer program for sonic-boom propagation that includes shock coalescence and incorporates effects of asymmetries due to volume and lift. Numerically integrates nonlinear governing equations using data on initial data line approximately one body length from aircraft and yields sonic-boom pressure at ground as function of time or of position at given time. MMOC written in FORTRAN IV.

Darden, C. M.↗

NASA NTRS1987

Optical communication with semiconductor laser diodes

A 25 megabit/sec direct detection optical communication system that used Q=4 PPM signalling was constructed and its performance measured under laboratory conditions. The system used a single-mode AlGaAs laser diode transmitter and low noise silicon avalanche photodiode (APD) photodetector. Comparison of measured performance with the theoretical revealed that modeling the APD output as a Gaussian process under conditions of negligible background radiation and low (less than 10 to the -12 power A) APD bulk leakage currents leads to substantial underestimates of optimal APD gain to use and overestimates of system bit error probability. A procedure is given to numerically compute system performance which uses the more accurate Webb's Approximation of the exact Conradi distribution for the APD ouput signal that does not require excessive amounts of computer time (a few minutes of VAX 8600 CPU time per system operating point). Examples are given which illustrate the breakdown of the Gaussian approximation in assessing system performance. This system achieved a bit error probability of 10 to the -6 power at a received signal energy corresponding to an average of 60 absorbed photons/bit and optimal APD gain of 700.

Davidson, F.↗

NASA NTRS1986

Nonlinear flexure and torsion of rotating beams, with application to helicopter rotor blades. I - Formulation. II - Response and stability results

The dynamic response and aeroelastic stability of rotating beams such as helicopter blades is investigated analytically. The Hamilton principle is used to formulate the equations of motion for extensional and inextensional beams with precone angles and variable pitch angles, taking higher-order nonlinearities into account. The derivation of the equations and their approximate solution by a Galerkin procedure are explained in detail, and numerical results of equilibrium solutions and stability analyses are presented graphically.

Crespo Da Silva, M. R. M.↗

NASA NTRS1988

A diagonally inverted LU implicit multigrid scheme

A new Diagonally Inverted LU Implicit scheme is developed within the framework of the multigrid method for the 3-D unsteady Euler equations. The matrix systems that are to be inverted in the LU scheme are treated by local diagonalizing transformations that decouple them into systems of scalar equations. Unlike the Diagonalized ADI method, the time accuracy of the LU scheme is not reduced since the diagonalization procedure does not destroy time conservation. Even more importantly, this diagonalization significantly reduces the computational effort required to solve the LU approximation and therefore transforms it into a more efficient method of numerically solving the 3-D Euler equations.

Yokota, Jeffrey W.↗

NASA NTRS1988

Improved 3-D turbomachinery CFD algorithm

The building blocks of a computer algorithm developed for the time-accurate flow analysis of rotating machines are described. The flow model is a finite volume method utilizing a high resolution approximate Riemann solver for interface flux definitions. This block LU implicit numerical scheme possesses apparent unconditional stability. Multi-block composite gridding is used to orderly partition the field into a specified arrangement. Block interfaces, including dynamic interfaces, are treated such as to mimic interior block communication. Special attention is given to the reduction of in-core memory requirements by placing the burden on secondary storage media. Broad applicability is implied, although the results presented are restricted to that of an even blade count configuration. Several other configurations are presently under investigation, the results of which will appear in subsequent publications.

Janus, J. Mark↗

NASA NTRS1987

An analytical solution to obtain the optimum source location using multiple direction finders on a spherical surface

An analytical solution is presented for determining the optimum location of a radiating source on the surface of a sphere, given multiple bearings. The bearings are assumed to have small errors of the order of 0-10 deg. The optimum location is found by minimizing the sum of the squares of the perpendicular great-circle distances from the source to the bearing lines. This is achieved analytically through an eigenvalue approach, rather than the usual iterative, numerical approach. Bearings of different weight are taken into account by approximating the distance from each direction finder to the source. The result is general and may have wide application. Since it is simple and nearly as fast as the triangulation technique for source location, it is now used in the SUNY-Albany East Coast Lightning Detection Network to compute the optimum location for lightning in real time.

Orville, Richard E., Jr.↗