SEARCH · Engineering Papers

Results for “numerical approximation & analysis”

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 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS2008

Unsteady Analysis of Blade and Tip Heat Transfer as Influenced by the Upstream Momentum and Thermal Wakes

The effect of the upstream wake on the blade heat transfer has been numerically examined. The geometry and the flow conditions of the first stage turbine blade of GE s E3 engine with a tip clearance equal to 2 percent of the span was utilized. Based on numerical calculations of the vane, a set of wake boundary conditions were approximated, which were subsequently imposed upon the downstream blade. This set consisted of the momentum and thermal wakes as well as the variation in modeled turbulence quantities of turbulence intensity and the length scale. Using a one-blade periodic domain, the distributions of unsteady heat transfer rate on the turbine blade and its tip, as affected by the wake, were determined. Such heat transfer coefficient distribution was computed using the wall heat flux and the adiabatic wall temperature to desensitize the heat transfer coefficient to the wall temperature. For the determination of the wall heat flux and the adiabatic wall temperatures, two sets of computations were required. The results were used in a phase-locked manner to compute the unsteady or steady heat transfer coefficients. It has been found that the unsteady wake has some effect on the distribution of the time averaged heat transfer coefficient on the blade and that this distribution is different from the distribution that is obtainable from a steady computation. This difference was found to be as large as 20 percent of the average heat transfer on the blade surface. On the tip surface, this difference is comparatively smaller and can be as large as four percent of the average.

Ameri, Ali A.↗

NASA NTRS1982

An efficient full potential implicit method based on characteristics for analysis of supersonic flows

A nonlinear aerodynamic prediction technique based on the full potential equation in conservation form has been developed for the treatment of supersonic flows. The method uses the theory of characteristic signal propagation to accurately simulate the flow structure, which includes shock waves and mixed elliptic-hyperbolic crossflow. An implicit approximate factorization scheme is employed to solve the finite-differenced equation. The necessary body-fitted grid system in every marching plane is generated numerically, using an elliptic grid solver. Results are shown for conical and nonconical wing-body combinations and compared with experimental data and Euler calculations. The method demonstrates an enormous savings in execution time and memory requirements over Euler methods.

Shankar, V.↗

NASA NTRS1992

Two dimensional analysis of a high temperature gaseous radiation receiver

The characteristics of the Flowing Gas Radiation Receiver (FGRR), a device that absorbs solar radiation volumetrically in a gas to produce high temperatures for space propulsion and power applications, are analyzed using a two-dimensional axisymmetric numerical model of the flow and radiation fields within a diffusely reflecting channel. The results show that an FGRR system is capable of generating temperatures in excess of 3000 K with collection efficiencies of approximately 75 percent for a channel with a reflectivity of 0.9. For a collinear radiation source, outflow temperatures of 3193 and 3092 K were achieved for axial and radial flow inputs, respectively, with receiver efficiencies of 0.82 and 0.76.

Mcfall, K. A.↗

NASA NTRS1992

Navier-Stokes analysis of three-dimensional unsteady flows inside turbine stages

This study presents a numerical method for solving the 3D Navier-Stokes equations for unsteady, viscous flow through multiple turbomachinery blade rows. The method solves the fully 3D Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed. A third-order accurate upwinding scheme is used to approximate convection terms, while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flowfield of the High Pressure Fuel side Turbo-Pump (HPFTP) of the Space Shuttle Main Engine (SSME). The stage calculation is performed by coupling the stator and the rotor flowfields at each time step through an over-laid grid. Numerical results for the complete geometry with the vane trailing edge cutback are presented and compared with the available experimental data.

Hah, C.↗

NASA NTRS1975

Non-numeric computation for high eccentricity orbits

Geocentric orbits of large eccentricity (e = 0.9 to 0.95) are significantly perturbed in cislunar space by the sun and moon. The time-history of the height of perigee, subsequent to launch, is particularly critical. The determination of 'launch windows' is mostly concerned with preventing the height of perigee from falling below its low initial value before the mission lifetime has elapsed. Between the extremes of high accuracy digital integration of the equations of motion and of using an approximate, but very fast, stability criteria method, this paper is concerned with the developement of a method of intermediate complexity using non-numeric computation. The computer is used as the theory generator to generalize Lidov's theory using six osculating elements. Symbolic integration is completely automatized and the output is a set of condensed formulae well suited for repeated applications in launch window analysis. Examples of applications are given.

Sridharan, R.↗

NASA NTRS1980

An analysis of turbulent diffusion flame in axisymmetric jet

The kinetic theory of turbulent flow was employed to study the mixing limited combustion of hydrogen in axisymmetric jets. The integro-differential equations in two spatial and three velocity coordinates describing the combustion were reduced to a set of hyperbolic partial differential equations in the two spatial coordinates by a binodal approximation. The MacCormick's finite difference method was then employed for solution. The flame length was longer than that predicted by the flame-sheet analysis, and was found to be in general agreement with a recent experimental result. Increase of the turbulence energy and scale resulted in an enhancement of the combustion rate and, hence, in a shorter flame length. Details of the numerical method as well as of the physical findings are discussed.

Chung, P. M.↗

NASA NTRS1981

Instability and transition in rotating disk flow

The stability of three-dimensional rotating disk flow is investigated, including the effects of Coriolis forces and streamline curvature. The numerical results show that the critical Reynolds number for establishment of stationary vortex flow is 287. These vortices spiral outward at an angle of about 11.2 deg, and transition to turbulence occurs when their total amplification is about e to the 11th. New experimental results are also reported on the spatial growth rates of the stationary vortices. It is shown that the analysis gives growth rates that compare much better with the experimental results than do results obtained using the Orr-Sommerfeld equation. The experimental results tend to support the numerical prediction that the number of stationary vortices varies directly with the Reynolds number. The calculations also indicate the existence of weakly unstable propagating (type II) modes at low Reynolds numbers (Critical Reynolds Number being approximately equal to 49).

Malik, M. R.↗

NASA NTRS1993

Multigrid properties of upwind-biased data reconstructions

The multigrid properties of two data reconstruction methods used for achieving second-order spatial accuracy when solving the two-dimensional Euler equations are examined. The data reconstruction methods are used with an implicit upwind algorithm which uses linearized backward-Euler time-differencing. The solution of the resulting linear system is performed by an iterative procedure. In the present study only regular quadrilateral grids are considered, so a red-black Gauss-Seidel iteration is used. Although the Jacobian is approximated by first-order upwind extrapolation, two alternative data reconstruction techniques for the flux integral that yield higher-order spatial accuracy at steady state are examined. The first method, probably most popular for structured quadrilateral grids, is based on estimating the cell gradients using one-dimensional reconstruction along curvilinear coordinates. The second method is based on Green's theorem. Analysis and numerical results for the two dimensional Euler equations show that data reconstruction based on Green's theorem has superior multigrid properties as compared to the one-dimensional data reconstruction method.

Warren, Gary P.↗

NASA NTRS1974

Post-flight differential correction analysis using Vinti's spheroidal method for the small astronomy satellite orbit

The results of an intensive analysis of a differential orbit improvement method utilizing observational data for a 550-kilometer altitude, near-circular, near-equatorial satellite orbit are presented. Observations of the Small Astronomy Satellite (SAS-1) are in the form of direction cosines as measured at two ground interferometer tracking stations near the Equator during the first 22 orbital revolutions (approximately 37 hours) after launch of the spacecraft. Numerical results, in both tabular and graphical form, are displayed for numerous iterated fittings of various observational arcs by differential correction of the orbital elements. Parameters varied in these comparative cases include the time duration of the observational data block, the number of pairs of direction cosine data and the number of tracking station passes included in the solution, the distribution of such passes between the two available tracking stations, and the acceptance criterion for the observational residuals in the least squares fitting procedure. It was found that three observational pairs of direction cosine data, the minimum number possible for a uniquely determined solution in theory, are sufficient to promote covergence to an accurate solution, if properly selected.

Walden, H.↗

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 NTRS1999

Evidence for Methane Segregation at the Surface of Pluto

In May 1995, a set of spectrophotometric curves of the system Pluto-Charon was recorded with the UKIRT telescope equipped with the spectrometer CGS4. As for the previous observations, the spectra cover a part of the near infrared range, between 1.4 and 2.55 micrometers, but with a higher resolution of approximately 700. In both the 1992 and 1995 data, the existence of solid methane is confirmed by numerous absorption bands, and the carbon monoxide and the nitrogen ices are identified by their respective signatures at 2.35 and 2.15 um. The solid nitrogen seems to be the principal icy component and forms a matrix in which the CH4 and CO molecules are diluted. However a spectroscopic analysis of the 1995 observations indicates that pure methane may coexist with its diluted phase in N2. In order to derive the horizontal and vertical distribution of these different species and to obtain some quantitative information about their characteristics, we have modeled the spectrum of May 15 that corresponds to the maximum of Pluto's visible light curve. This was achieved by means of a radiative transfer algorithm dealing with compact and stratified media. Among the various representations we have tested to describe the surface of Pluto, only a geographical mixture of three distinct units explains all the significant structures of the analyzed spectrum. The first unit is a thin granular layer of pure CH4 covering a compact polycrystalline substratum of N2-CH4-CO, which are in a molecular mixture (concentrations of and CO of the order of 0.45%, 0.1-0.2% respectively). It covers about 70% of the observed area and corresponds to volatile deposits that are sublimating under solar illumination. The second unit is either (a) a single thick layer of pure granular methane or (b) a unit similar to the first unit but with the two components inverted (i.e. with CH4 forming a substratum and the N2-CH4-CO mixture a superficial layer of fine grains). Covering 20% of the surface, it represents some old surfaces that have been sublimated for a long time, and eventually recovered later by very small amounts of fresh deposits of the molecular mixture N2-CH4-CO. Finally, the third unit may result from the condensation of very fine grains of nearly pure N2. It covers the remainder of the surface (about 10%). All these results allow a better understanding of the processes of deposition, metamorphism, sublimation and transport affecting the different ices detected on Pluto during its climatic cycles.

Doute, S.↗

NASA NTRS1996

Theoretical and Numerical Investigation of Radiative Extinction of Diffusion Flames

The influence of soot radiation on diffusion flames was investigated using both analytical and numerical techniques. Soot generated in diffusion flames dominate the flame radiation over gaseous combustion products and can significantly lower the temperature of the flame. In low gravity situations there can be significant accumulation of soot and combustion products in the vicinity of the primary reaction zone owing to the absence of any convective buoyant flow. Such situations may result in substantial suppression of chemical activities in a flame, and the possibility of a radiative extinction may also be anticipated. The purpose of this work was to not only investigate the possibility of radiative extinction of a diffusion flame but also to qualitatively and quantitatively analyze the influence of soot radiation on a diffusion flame. In this study, first a hypothetical radiative loss profile of the form of a sech(sup 2) was assumed to influence a pure diffusion flame. It was observed that the reaction zone can, under certain circumstances, move through the radiative loss zone and locate itself on the fuel side of the loss zone contrary to our initial postulate. On increasing the intensity and/or width of the loss zone it was possible to extinguish the flame, and extinction plots were generated. In the presence of a convective flow, however, the movement of the temperature and reaction rate peaks indicated that the flame behavior is more complicated compared to a pure diffusional flame. A comprehensive model of soot formation, oxidation and radiation was used in a more involved analysis. The soot model of Syed, Stewart and Moss was used for soot nucleation and growth and the model of Nagle and Strickland-Constable was used for soot oxidation. The soot radiation was considered in the optically thin limit. An analysis of the flame structure revealed that the radiative loss term is countered both by the reaction term and the diffusion term. The essential balance for the soot volume fraction was found to be between the processes of soot convection and soot growth. Such a balance yielded to analytical treatment and the soot volume fraction could be expressed in the form of an integral. The integral was evaluated using two approximate methods and the results agreed very well with the numerical solutions for all cases examined.

Ray, Anjan↗

NASA NTRS1995

The analysis of control trajectories using symbolic and database computing

This final report comprises the formal semi-annual status reports for this grant for the periods June 30-December 31, 1993, January 1-June 30, 1994, and June 1-December 31, 1994. The research supported by this grant is broadly concerned with the symbolic computation, mixed numeric-symbolic computation, and database computation of trajectories of dynamical systems, especially control systems. A review of work during the report period covers: trajectories and approximating series, the Cayley algebra of trees, actions of differential operators, geometrically stable integration algorithms, hybrid systems, trajectory stores, PTool, and other activities. A list of publications written during the report period is attached.

Grossman, Robert↗

NASA NTRS1988

Analysis of long term heart rate variability: methods, 1/f scaling and implications

The use of spectral techniques to quantify short term heart rate fluctuations on the order of seconds to minutes has helped define the autonomic contributions to beat-to-beat control of heart rate. We used similar techniques to quantify the entire spectrum (0.00003-1.0 Hz) of heart rate variability during 24 hour ambulatory ECG monitoring. The ECG from standard Holter monitor recordings from normal subjects was sampled with the use of a phase locked loop, and a heart rate time series was constructed at 3 Hz. Frequency analysis of the heart rate signal was performed after a nonlinear filtering algorithm was used to eliminate artifacts. A power spectrum of the entire 24 hour record revealed power that was inversely proportional to frequency, 1/f, over 4 decades from 0.00003 to 0.1 Hz (period approximately 10 hours to 10 seconds). Displaying consecutive spectra calculated at 5 minute intervals revealed marked variability in the peaks at all frequencies throughout the 24 hours, probably accounting for the lack of distinct peaks in the spectra of the entire records.

Non-NASA Center↗

NASA NTRS1980

Note on the eigensolution of a homogeneous equation with semi-infinite domain

The 'variation-iteration' method using Green's functions to find the eigenvalues and the corresponding eigenfunctions of a homogeneous Fredholm integral equation is employed for the stability analysis of fluid hydromechanics problems with a semiinfinite (infinite) domain of application. The objective of the study is to develop a suitable numerical approach to the solution of such equations in order to better understand the full set of equations for 'real-world' flow models. The study involves a search for a suitable value of the length of the domain which is a fair finite approximation to infinity, which makes the eigensolution an approximation dependent on the length of the interval chosen. In the examples investigated y = 1 = a seems to be the best approximation of infinity; for y greater than unity this method fails due to the polynomial nature of Green's functions.

Wadia, A. R.↗

NASA NTRS1982

Heat generation in aircraft tires under free rolling conditions

A method was developed for calculating the internal temperature distribution in an aircraft tire while free rolling under load. The method uses an approximate stress analysis of each point in the tire as it rolls through the contact patch, and from this stress change the mechanical work done on each volume element may be obtained and converted into a heat release rate through a knowledge of material characteristics. The tire cross-section is then considered as a body with internal heat generation, and the diffusion equation is solved numerically with appropriate boundary conditions of the wheel and runway surface. Comparison with data obtained with buried thermocouples in tires shows good agreement.

Clark, S. K.↗

NASA NTRS1983

Heat generation in aircraft tires

Clark, S. K.↗

NASA NTRS1984

Optimum design sensitivity based on approximation concepts and dual methods

An optimum design sensitivity analysis capability is reported which exploits the approximation concepts-dual method formulation of the minimum weight structural sizing problem. An efficient iterative solution technique is used to facilitate determination of sensitivity derivatives for both primal and dual variables. Estimates on the useful range of parameter perturbations, over which the optimum design sensitivity projections can be expected to yield satisfactory revised optimum designs, are also obtained. Numerical results for several example problems will be presented to illustrate the effectiveness of the capability reported.

Schmit, L. A.↗