Search NASA⌕ Search

SEARCH · Search NASA

Results for “higher-order integration”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Higher-order numerical methods derived from three-point polynomial interpolation

Higher-order collocation procedures resulting in tridiagonal matrix systems are derived from polynomial spline interpolation and Hermitian finite-difference discretization. The equations generally apply for both uniform and variable meshes. Hybrid schemes resulting from different polynomial approximations for first and second derivatives lead to the nonuniform mesh extension of the so-called compact or Pade difference techniques. A variety of fourth-order methods are described and this concept is extended to sixth-order. Solutions with these procedures are presented for the similar and non-similar boundary layer equations with and without mass transfer, the Burgers equation, and the incompressible viscous flow in a driven cavity. Finally, the interpolation procedure is used to derive higher-order temporal integration schemes and results are shown for the diffusion equation.

Rubin, S. G.↗
NASA NTRS1979

A progress report on a NASA research program for embedded computer systems software

The paper presents the results of the second stage of the Multipurpose User-oriented Software Technology (MUST) program. Four primary areas of activities are discussed: programming environment, HAL/S higher-order programming language support, the Integrated Verification and Testing System (IVTS), and distributed system language research. The software development environment is provided by the interactive software invocation system. The higher-order programming language (HOL) support chosen for consideration is HAL/S mainly because at the time it was one of the few HOLs with flight computer experience and it is the language used on the Shuttle program. The overall purpose of IVTS is to provide a 'user-friendly' software testing system which is highly modular, user controlled, and cooperative in nature.

Foudriat, E. C.↗
NASA NTRS1987

A computational study of thrust augmenting ejectors based on a viscous-inviscid approach

A viscous-inviscid interaction technique is advocated as both an efficient and accurate means of predicting the performance of two-dimensional thrust augmenting ejectors. The flow field is subdivided into a viscous region that contains the turbulent jet and an inviscid region that contains the ambient fluid drawn into the device. The inviscid region is computed with a higher-order panel method, while an integral method is used for the description of the viscous part. The strong viscous-inviscid interaction present within the ejector is simulated in an iterative process where the two regions influence each other en route to a converged solution. The model is applied to a variety of parametric and optimization studies involving ejectors having either one or two primary jets. The effects of nozzle placement, inlet and diffuser shape, free stream speed, and ejector length are investigated. The inlet shape for single jet ejectors is optimized for various free stream speeds and Reynolds numbers. Optimal nozzle tilt and location are identified for various dual-ejector configurations.

Lund, Thomas S.↗
NASA NTRS1983

Calculations of atmospheric refraction for spacecraft remote-sensing applications

Analytical solutions to the refraction integrals appropriate for ray trajectories along slant paths through the atmosphere are derived in this paper. This type of geometry is commonly encountered in remote-sensing applications utilizing an occultation technique. The solutions are obtained by evaluating higher-order terms from expansion of the refraction integral and are dependent on the vertical temperature distributions. Refraction parameters such as total refraction angles, air masses, and path lengths can be accurately computed. It is also shown that the method can be used for computing refraction parameters in astronomical refraction geometry for large zenith angles.

Chu, W. P.↗
NASA NTRS1989

Advanced development of the boundary element method for steady-state heat conduction

Considerable progress has been made in recent years toward advancing the state-of-the-art in solid mechanics boundary element technology. In the present work, much of this new technology is applied in the development of a general-purpose boundary element method (BEM) for steady-state heat conduction. In particular, the BEM implementation involves the use of higher-order conforming elements, self-adaptive integration and multi-region capability. Two- and three-dimensional, as well as axisymmetric analysis, are incorporated within a unified framework. In addition, techniques are introduced for the calculation of boundary flux, and for the inclusion of thermal resistance across interfaces. As a final extension, an efficient formulation is developed for the analysis of solid three-dimensional bodies with embedded holes. For this last class of problems, the new BEM formulation is particularly attractive, since use of the alternatives (i.e. finite element or finite difference methods) is not practical. A number of detailed examples illustrate the suitability and robustness of the present approach for steady-state heat conduction.

Dargush, G. F.↗
NASA NTRS1989

Derivation of a Tappered p-Version Beam Finite Element

A tapered p-version beam finite element suitable for dynamic applications is derived. The taper in the element is represented by allowing the area moments of inertia to vary as quartic polynomials along the length of the beam, and the cross-sectional area to vary as a quadratic polynomial. The p-version finite-element characteristics are implemented through a set of polynomial shape functions. The lower-order shape functions are identical to the classical cubic and linear shape functions normally associated with a beam element. The higher-order shape functions are a hierarchical set of polynomials that are integrals of orthogonal polynomials. Explicit expressions for the mass and stiffness matrices are presented for an arbitrary value of p. The element has been verified to be numerically stable using shape functions through 22nd order.

Hinnant, Howard E.↗
NASA NTRS1987

A viscous-inviscid interaction method for ejectors having two primary jets

A viscous-inviscid interaction methodology is developed for evaluating the performance of two-dimensional incompressible flow ejectors that utilize two primary jets. In this method, the inviscid portion of the flow field, which contains the ambient fluid drawn into the device, is modeled using a higher-order panel method. The viscous region, which contains the turbulent jets, is described by an integral method. The viscous and inviscid solutions are matched together in an iterative process that simulates the interaction between the jets and the ambient fluid. The model is applied to a variety of parametric and optimization studies that illustrate the connection between the details of the ejector geometry and the level of thrust augmentation. The advantages of a dual-jet configuration are assessed through comparison with similar calculations for an otherwise identical single-jet ejector.

Lund, T. S.↗
NASA NTRS1987

An evaluation of higher-order modal methods for calculating transient structural response

The present study evaluates a higher-order modal method proposed by Leung for transient structural analysis entitled the force-derivative method. This method repeatedly integrates by parts with respect to time the convolution-integral form of the structural response to produce successively better approximations to the contribution of the higher modes which are neglected in the modal summation. Comparisons are made of the force-derivative, the mode-displacement, and the mode-acceleration methods for several numerical example problems for various times, levels of damping, and forcing functions. The example problems include a tip-loaded cantilevered beam and a simply-supported multispan beam. The force-derivative method is shown to converge to an accurate solution in fewer modes than either the mode-displacement or the mode-acceleration methods. In addition, for problems in which there are a large number of closely-spaced frequencies whose mode shapes have a negligible contribution to the response, the force derivative method is very effective in representing the effect of the important, but otherwise neglected, higher modes.

Camarda, Charles J.↗
NASA NTRS1987

An evaluation of higher-order model methods for calculating transient structural response

A higher-order modal method proposed by Leung for transient structural analysis entitled the force-derivative method is evaluated. This method repeatedly integrates by parts with respect to time the convolution-integral form of the structural response to produce successively better approximations to the contribution of the higher modes which are neglected in the modal summation. Comparisons are made of the force-derivative, the mode-displacement, and the mode-acceleration methods for several numerical example problems for various times, levels of damping, and forcing functions. The example problems include a tip-loaded cantilevered beam and a simply-supported multispan beam. The force-derivative method is shown to converge to an accurate solution in fewer modes than either the mode-displacement or the mode-acceleration methods. In addition, for problems in which there are a large number of closely-spaced frequencies whose mode shapes have a negligible contribution to the response, the force-derivative method is very effective in representing the effect of the important, but otherwise neglected, higher modes.

Camarda, Charles J.↗
NASA NTRS1988

Tunneling near the peaks of potential barriers - Consequences of higher-order Wentzel-Kramers-Brillouin corrections

The quantum-mechanical tunneling transmission coefficient is derived for energies near the peak of a general potential barrier, to fourth order in the WKB approximation, using a method based on an eighth-order Taylor expansion of the potential near the peak. The result agrees with contour-integral formulations of WKB tunneling. For the Poeschl-Teller potential the WKB series converges rapidly to the known exact transmission coefficient. Higher-order WKB corrections may be important in attempts to determine the internuclear potential by inversion using low-energy fusion cross-section data.

Will, Clifford M.↗
NASA NTRS1987

Actuators for actively controlled space structures

A free-flying platform of about 4 x 17 m overall dimensions, carrying a variety of imaging and sounding payloads, calls for an intelligent structure with active dynamic control of structural resonances. The actuators for such a structure must be lightweight, require low power, and allow integration into the structure without degradation of its integrity; the dc-to-100 Hz dynamic range required may entail several types of actuators, as is presently emphasized. Broadband damping of higher-order modes requires modeling of the structure with a distributed array of sensors and actuators.

Studer, P.↗
NASA NTRS1973

Auroral ion velocity distributions using a relaxation model.

Calculation of ion velocity distributions for a weakly-ionized plasma subjected to crossed electric and magnetic fields for application to the auroral ionosphere. By replacing the Boltzmann collision integral with a simple relaxation model, an exact solution to Boltzmann's equation could be obtained. This solution has the advantage over a series expansion in that all the higher-order velocity moments are inherent in it. The exact solution is particularly advantageous when studying large departures of the distribution from its Maxwellian form, because these departures are caused by the higher velocity moments. In general, however, a simple relaxation model can only be used to obtain qualitative information on the distribution function. Consequently, it is possible to determine when the higher-order velocity moments affect the ion velocity distribution and the nature of their effect, but it is not possible to obtain accurate quantitative results.

St-Maurice, J.-P.↗
NASA NTRS1976

Kinetic theory of turbulent compressible flows and comparison with classical theory

A complete set of turbulent correlations is given in terms of expansion coefficients in a double series of Hermite polynomials of a two-particle correlation in phase space. Only two of these coefficients, corresponding to Reynolds stress and turbulent heat flux, are shown to appear in gasdynamic equations of turbulence. This is accomplished by incorporating the phase-space correlation in the collision integral of Boltzmann's equation and by deriving generalized Navier-Stokes and Fourier transport relations. In contrast to the classical formalism, the simple expression remains invariant in form whether the flow is compressible or not. Nevertheless, a bilinear transformation of fluctuating quantities shows that the two formalisms are identical with regard to terms of double correlations. The kinetic theory justifies nonexistence of the higher-order correlations characteristic of compressible-turbulence equations in the classical regime.

Tsuge, S.↗
NASA NTRS1989

Euler and Navier-Stokes solutions for hypersonic flows

An upwind finite-element technique that uses cell-centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured grids in two and three dimensions. A perfect gas model as well as an equilibrium air model is implemented for solving high-speed flows. A first-order basic scheme and a higher-order flux-corrected transport (FCT) scheme have been implemented. This technique has been used to predict 'Type III and IV' shock interactions on a cylinder in two dimensions and a swept cylinder in three dimensions, with a view to determine the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively-generated grids. Three-dimensional corner flow, typically encountered in engine inlets due to compression of the flow by ramps in the walls, is also modeled. This procedure is the first step in developing an integrated fluid, thermal, structural analysis capability for hypersonic flight vehicles like the National Aero-Space Plane.

Thareja, Rajiv R.↗
NASA NTRS1977

Use of finite element techniques in the determination of the acoustic properties of turbofan inlets

This paper describes the application of the finite-element method in combination with Galerkin's method in the determination of the acoustic properties of turbofan inlets containing high-subsonic-Mach-number flows. An approximate solution to the steady inviscid flow field is obtained using an integral method for calculating the potential flow field in the inlet with a correction to account for compressibility effects. The accuracy of the finite-element technique in predicting the acoustic properties of annular ducts has been checked by comparison with available analytical solutions for the problems of planeand spinning-wave propagation through a hard-walled annular duct with a constant mean flow. Results are presented comparing low-frequency plane-wave propagation through a hard-walled turbofan inlet containing a one-dimensional flow with the same inlet containing a fully two-dimensional axisymmetric mean flow. It is shown that when one-dimensional mean flow is assumed to exist in the duct, the plane wave propagates with relatively little distortion. However, propagation of a plane wave through the fully two-dimensional flow field in the inlet produces severe distortion due to the excitation of higher-order modes.

Sigman, R. K.↗
NASA NTRS1990

A real-time approximate optimal guidance law for flight in a plane

A real-time guidance scheme is presented for the problem of maximizing the payload into orbit subject to the equations of motion of a rocket over a nonrotating spherical earth. The flight is constrained to a path in the equatorial plane while reaching an orbital altitude at orbital injection speeds. The dynamics of the problem can be separated into primary and perturbation effects by a small parameter, epsilon, which is the ratio of the atmospheric scale height to the radius of the earth. The Hamilton-Jacobi-Bellman or dynamic programming equation is expanded in an asymptotic series where the zeroth-order term (epsilon = 0) can be obtained in closed form. The neglected perturbation terms are included in the higher-order terms of the expansion, which are determined from the solution of first-order linear partial differential equations requiring only integrations which are quadratures. The quadratures can be performed rapidly with emerging computer capability, so that real-time approximate optimization can be used to construct the launch guidance law. The application of this technique to flight in three-dimensions is made apparent from the solution presented.

Feeley, Timothy S.↗
NASA NTRS1942

On the use of residue theory for treating the subsonic flow of a compressible fluid

A new mathematical technique, due to Milne-Thomson, is used to obtain an improved form of the method of Poggi for calculating the effect of compressibility on the subsonic flow past an obstacle. By means of this new method, the difficult surface integrals of the original Poggi method can be replaced by line integrals. These line integrals are then solved by the use of residue theory. In this way an equation is obtained giving the second-order effect of compressibility on the velocity of the fluid. The method is practicable for obtaining the higher-order effects of compressibility on the velocity field. As an illustration of the general result, the flow past an elliptic cylinder is discussed.

Kaplan, Carl↗