Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical 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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

A multistage time-stepping scheme for the thin-layer Navier-Stokes equations

A finite-volume scheme for numerical integration of the Euler equations was extended to allow solution of the thin-layer Navier-Stokes equations in two and three dimensions. The extended algorithm, which is based on a class of four-stage Runge-Kutta time-stepping schemes, was made numerically efficient through the following convergence acceleration technique: (1) local time stepping, (2) enthalpy damping, and (3) residual smoothing. Also, the high degree of vectorization possible with the algorithm has yielded an efficient program for vector processors. The scheme was evaluated by solving laminar and turbulent flows. Numerical results have compared well with either theoretical or other numerical solutions and/or experimental data.

Swanson, R. C., Jr.↗
NASA NTRS1991

Ignition and structure of a laminar diffusion flame in a compressible mixing layer with finite rate chemistry

The ignition and structure of a reacting compressible mixing layer is considered using finite rate chemistry lying between two streams of reactants with different freestream speeds and temperatures. Numerical integration of the governing equations show that the structure of the reacting flow can be quite complicated depending on the magnitude of the Zeldovich number. An analysis of both the ignition a diffusion flame regimes is presented using a combination of large Zeldovich number asymptotics and numerics. This allows to analyze the behavior of these regimes as a function of the parameters of the problem.

Grosch, C. E.↗
NASA NTRS1991

Ignition and structure of a laminar diffusion flame in a compressible mixing layer with finite rate chemistry

The ignition and structure of a reacting compressible mixing layer is considered using finite rate chemistry lying between two streams of reactants with different freestream speeds and temperatures. Numerical integration of the governing equations show that the structure of the reacting flow can be quite complicated depending on the magnitude of the Zeldovich number. An analysis of both the ignition and diffusion flame regimes is presented using a combination of large Zeldovich number asymptotics and numerics. This allows to analyze the behavior of these regimes as a function of the parameters of the problem.

Grosch, C. E.↗
NASA NTRS1992

Numerical analysis of the hot-gas-side and coolant-side heat transfer in liquid rocket engine combustors

The objectives of this paper are to develop computational methods to predict the hot-gas-side and coolant-side heat transfer, and to use these methods in parametric studies to recommend optimized design of the coolant channels for regeneratively cooled liquid rocket engine combustors. An integrated numerical model which incorporates computational fluid dynamics (CFD) for the hot-gas thermal environment, and thermal analysis for the coolant channels, was developed. The mode was validated by comparing predicted heat fluxes with those of hot-firing test and industrial design methods. Parametric studies were performed to find a strategy for optimized combustion chamber coolant channel design.

Wang, T. S.↗
NASA NTRS1975

Accurate integration of geostationary orbits with Burdet's focal elements

The theory of Burdet's focal elements (1969) is outlined. The differential equations are presented, and the initial value problem is described together with the transformation to rectangular coordinates and classical elements. The focal elements are well defined for zero eccentricity and inclination and can be adopted for the computation of elliptic, parabolic, and hyperbolic motion. For the numerical integration of near-geostationary orbits, a comparison of the efficiency is made between focal elements, Kustaanheimo-Stiefel (1971) theory, and rectangular coordinates. For this class of orbits, a higher accuracy has been obtained by integrating elements than integrating rectangular coordinates.

Flury, W.↗
NASA NTRS1948

Determination of stresses in gas-turbine disks subjected to plastic flow and creep

A finite-difference method previously presented for computing elastic stresses in rotating disks is extended to include the computation of the disk stresses when plastic flow and creep are considered. A finite-difference method is employed to eliminate numerical integration and to permit nontechnical personnel to make the calculations with a minimum of engineering supervision. Illustrative examples are included to facilitate explanation of the procedure by carrying out the computations on a typical gas-turbine disk through a complete running cycle. The results of the numerical examples presented indicate that plastic flow markedly alters the elastic-stress distribution.

Millenson, M B↗
NASA NTRS1971

Comparison of a discrete steepest ascent method with the continuous steepest ascent method for optimal programing

A discrete steepest ascent method which allows controls which are not piecewise constant (for example, it allows all continuous piecewise linear controls) was derived for the solution of optimal programming problems. This method is based on the continuous steepest ascent method of Bryson and Denham and new concepts introduced by Kelley and Denham in their development of compatible adjoints for taking into account the effects of numerical integration. The method is a generalization of the algorithm suggested by Canon, Cullum, and Polak with the details of the gradient computation given. The discrete method was compared with the continuous method for an aerodynamics problem for which an analytic solution is given by Pontryagin's maximum principle, and numerical results are presented. The discrete method converges more rapidly than the continuous method at first, but then for some undetermined reason, loses its exponential convergence rate. A comparsion was also made for the algorithm of Canon, Cullum, and Polak using piecewise constant controls. This algorithm is very competitive with the continuous algorithm.

Childs, A. G.↗
NASA NTRS1999

The Determination and Long Term Integration of the Orbits of Caliban and Sycorax

The first 2 irregular satellites of Uranus, Caliban and Sycorax, were discovered in late 1997. Subsequently, pre-discovery observations of both satellites were found on plates taken by D. Cruikshank in June of 1984. Recently, P. Nicholson, D. Tholen, and W. Offutt provided observations which they made in late 1998 at Palomar Mountain, Mauna Kea, and Cloudcroft, respectively. I fit a numerical integration perturbed by the Sun, Jupiter, Saturn, and Neptune to the set of available observations. For the 47 observations of Caliban the respective rms values of the Delta alpha (cos delta) and Delta delta residuals axe 0".60 and 0".32, and for the 103 observations of Sycorax the analogous values are 0".57 and 0".59. I extended the integration to span a 6000 year period and computed osculating orbital elements at yearly intervals. An included table contains the mean values of the elements over the 6000 years, the sidereal period, and the precession periods of the argument of periapsis and longitude of the ascending node. The osculating elements (except for a) exhibit a significant long period oscillation with a period roughly half that of the argument of periapsis.

Jacobson, R. A.↗
NASA NTRS1978

A note on obtaining an approximate value for the period of a periodic solution of x-double prime = -V prime /x/

The first integral of x-double prime = - V prime (x) yields an integral for the period of a periodic solution, if such exists. In general, this integral cannot be evaluated. By means of an approximate solution along with the minimization of a mean-square error, one can obtain an approximate value for the period in terms of the amplitude of the motion. The calculated period agrees very well with the period obtained by means of numerical integration for the case of orbit-orbit resonance involving the motion of two satellites of a planet. The same method is applied to obtain an alternative derivation of the first Krylov-Bogoliuboff averaging method in nonlinear mechanics.

Lass, H.↗
NASA NTRS2011

Shuttle Propulsion Overview - The Design Challenges

The major elements of the Space Shuttle Main Propulsion System include two reusable solid rocket motors integrated into recoverable solid rocket boosters, an expendable external fuel and oxidizer tank, and three reusable Space Shuttle Main Engines. Both the solid rocket motors and space shuttle main engines ignite prior to liftoff, with the solid rocket boosters separating about two minutes into flight. The external tank separates, about eight and a half minutes into the flight, after main engine shutdown and is safely expended in the ocean. The SSME's, integrated into the Space Shuttle Orbiter aft structure, are reused after post landing inspections. The configuration is called a stage and a half as all the propulsion elements are active during the boost phase, with only the SSME s continuing operation to achieve orbital velocity. Design and performance challenges were numerous, beginning with development work in the 1970's. The solid rocket motors were large, and this technology had never been used for human space flight. The SSME s were both reusable and very high performance staged combustion cycle engines, also unique to the Space Shuttle. The multi body side mount configuration was unique and posed numerous integration and interface challenges across the elements. Operation of the system was complex and time consuming. This paper describes the design challenges and key areas where the design evolved during the program.

Owen, James W.↗
NASA NTRS2005

Spinning Spacecraft Attitude Estimation Using Markley Variables: Filter Implementation And Results

Attitude estimation is often more difficult for spinning spacecraft than for three-axis stabilized platforms due to the need to follow rapidly-varying state vector elements and the lack of three-axis rate measurements from gyros. The estimation problem simplifies when torques are negligible and nutation has damped out, but the general case requires a sequential filter with dynamics propagation. This paper describes the implementation and test results for an extended Kalman filter for spinning spacecraft attitude and rate estimation based on a novel set of variables suggested in a paper by Markley [AAS93-3301 (referred to hereafter as Markley variables). Markley has demonstrated that the new set of variables provides a superior parameterization for numerical integration of the attitude dynamics for spinning or momentum-biased spacecraft. The advantage is that the Markley variables have fewer rapidly-varying elements than other representations such as the attitude quaternion and rate vector. A filter based on these variables was expected to show improved performance due to the more accurate numerical state propagation. However, for a variety of test cases, it has been found that the new filter, as currently implemented, does not perform significantly better than a quaternion-based filter that was developed and tested in parallel. This paper reviews the mathematical background for a filter based on Markley variables. It also describes some features of the implementation and presents test results. The test cases are based on a mission using magnetometer and Sun sensor data and gyro measurements on two axes normal to the spin axis. The orbit and attitude scenarios and spacecraft parameters are modeled after one of the THEMIS (Time History of Events and Macroscale Interactions during Substorms) probes. Several tests are presented that demonstrate the filter accuracy and convergence properties. The tests include torque-free motion with various nutation angles, large constant-torque attitude slews, sensor misalignments, large initial attitude and rate errors, and cases with low data frequency. It is found that the convergence is rapid, the radius of convergence is large, and the results are reasonably accurate even in the presence of unmodeled perturbations.

Sedlak, Joseph E.↗
NASA NTRS1972

The P prime(f) to N(h) inversion problem in ionospheric soundings

A general review is given of the inversion techniques used to derive the ionospheric electron density N as a function of altitude h from group path P versus frequency f measurements obtained by vertical incidence ionospheric sounders. The paper discusses the medium under investigation, the experimental techniques used to obtain the P(f) data, the theoretical considerations leading to the integral equation relating P(f) to N(h), and the assumptions made in the inversion process. The lamination inversion technique is then presented, with special attention given to mathematical difficulties arising from discontinuities in the P(f) function, infinities in the integrand, and in some cases unknown integration limits. Methods outlined for minimizing the uncertainties due to discontinuities include the use of redundant information - that is the two distinct P(f) functions available for a given N(h) profile - and the use of models based on statistical data. Mathematical procedures are discussed that increase significantly the efficiency and accuracy of the required numerical integrations. The accuracy of the inversion technique is deduced by comparing the resulting N(h) profile with N(h) data obtained by simultaneous but independent observations.

Jackson, J. E.↗
NASA NTRS2006

Dynamic Inversion based Control of a Docking Mechanism

The problem of position and attitude control of the Stewart platform based docking mechanism is considered motivated by its future application in space missions requiring the autonomous docking capability. The control design is initiated based on the framework of the intelligent flight control architecture being developed at NASA Ames Research Center. In this paper, the baseline position and attitude control system is designed using dynamic inversion with proportional-integral augmentation. The inverse dynamics uses a Newton-Euler formulation that includes the platform dynamics, the dynamics of the individual legs along with viscous friction in the joints. Simulation results are presented using forward dynamics simulated by a commercial physics engine that builds the system as individual elements with appropriate joints and uses constrained numerical integration,

Kulkarni, Nilesh V.↗
NASA NTRS1974

A procedure for verifying the structural integrity of an existing pressurized wind tunnel

The present work describes the approach employed in verifying the structural integrity at maximum operating pressure of a low turbulence pressurized wind tunnel built in 1939-1940. The approach utilized existing ASME Code guidelines and fracture mechanics techniques in combination with the latest nondestructive examination methods to validate the desired life of the facility while minimizing reconstruction or major modifications. Detailed stress analyses using a finite element computer program and a numerical integration program were performed on critical areas of the tunnel. Critical vessel areas were evaluated for remaining fatigue life. Fracture mechanics analyses of the shell were performed at various regions, and critical flaw sizes at each location were determined. Analysis of the test data permitted the remaining useful life of the tunnel to be estimated at 8 to 10 years.

Taylor, J. T.↗
NASA NTRS1976

Algorithms for computing the geopotential using a simple density layer

Several algorithms have been developed for computing the potential and attraction of a simple density layer. These are numerical cubature, Taylor series, and a mixed analytic and numerical integration using a singularity-matching technique. A computer program has been written to combine these techniques for computing the disturbing acceleration on an artificial earth satellite. A total of 1640 equal-area, constant surface density blocks on an oblate spheroid are used. The singularity-matching algorithm is used in the subsatellite region, Taylor series in the surrounding zone, and numerical cubature on the rest of the earth.

Morrison, F.↗
NASA NTRS1982

A study of viscous flow in stator and rotor passages

A numerical technique to predict the viscous behavior in the three-dimensional domain of rotating and stationary blade passages is presented. The analysis is based upon the numerical integration of the incompressible Navier-Stokes equations. An approach is taken to combine the three momentum equations and the continuity equation into two sets of vorticity transport-stream function equations. These equations are expressed along families of arbitrarily defined orthogonally intersecting surfaces placed within the blade passage. Unlike the stream surfaces, the arbitrary surfaces allow fluid to flow across. The numerical code employs a nonorthogonal body fitted coordinate system on each of the intersecting surfaces. The solution process involves the iterative combination of the solutions on the two families of intersecting surfaces for a complete solution. The results of flow in thin blade stationary and rotating polar ducts is presented.

Sheoran, Y.↗
NASA NTRS1988

Dynamic analysis of nonlinear rotor-housing systems

Nonlinear analysis methods are developed which will enable the reliable prediction of the dynamic behavior of the space shuttle main engine (SSME) turbopumps in the presence of bearing clearances and other local nonlinearities. A computationally efficient convolution method, based on discretized Duhamel and transition matrix integral formulations, is developed for the transient analysis. In the formulation, the coupling forces due to the nonlinearities are treated as external forces acting on the coupled subsystems. Iteration is utilized to determine their magnitudes at each time increment. The method is applied to a nonlinear generic model of the high pressure oxygen turbopump (HPOTP). As compared to the fourth order Runge-Kutta numerical integration methods, the convolution approach proved to be more accurate and more highly efficient. For determining the nonlinear, steady-state periodic responses, an incremental harmonic balance method was also developed. The method was successfully used to determine dominantly harmonic and subharmonic responses fo the HPOTP generic model with bearing clearances. A reduction method similar to the impedance formulation utilized with linear systems is used to reduce the housing-rotor models to their coordinates at the bearing clearances. Recommendations are included for further development of the method, for extending the analysis to aperiodic and chaotic regimes and for conducting critical parameteric studies of the nonlinear response of the current SSME turbopumps.

Noah, Sherif T.↗
NASA NTRS2010

An Improved Neutron Transport Algorithm for HZETRN

Long term human presence in space requires the inclusion of radiation constraints in mission planning and the design of shielding materials, structures, and vehicles. In this paper, the numerical error associated with energy discretization in HZETRN is addressed. An inadequate numerical integration scheme in the transport algorithm is shown to produce large errors in the low energy portion of the neutron and light ion fluence spectra. It is further shown that the errors result from the narrow energy domain of the neutron elastic cross section spectral distributions, and that an extremely fine energy grid is required to resolve the problem under the current formulation. Two numerical methods are developed to provide adequate resolution in the energy domain and more accurately resolve the neutron elastic interactions. Convergence testing is completed by running the code for various environments and shielding materials with various energy grids to ensure stability of the newly implemented method.

Slaba, Tony C.↗