Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Numerical simulation of Jet-A combustion approximated by improved propane chemical kinetics

Through the effort devoted to the chemical kinetics for propane air combustion, three mechanisms are developed. The full mechanism consists of 131 reactions. This mechanism is used as a guide for the evaluation of other mechanisms, but because of the long expected cpu time, it is not to be incorporated into the computer code KIVA-II for actual simulation. Through the sensitivity analysis, a reduced mechanism of 45 reactions is produced. But the calculated results from the 45 reaction mechanism are always low in temperature. Some efforts are devoted to correct this situation and details are included in this report. A simplified mechanism of reactions is successfully improved and computed results are compared with experimental data. Contour plots of physical parameters and species concentrations and results for emission indices of CO and NOx are presented.

Ying, Shuh-Jing↗
NASA NTRS1992

A variational approach to three-dimensional elasticity solutions of laminated composite plates

The displacements in a laminated composite are represented as products of two sets of unknown functions, one of which is only a function of the thickness coordinate and the other is a function of the in-plane coordinates (i.e., separation of variables approach), and the minimization of the total potential energy is reduced to a sequence of iterative linear problems. Analytical solutions are developed for cross-ply and angle-ply laminated composite rectangular plates. The solution for simply-supported cross-ply plates under sinusoidal transverse load reduces to that of Pagano. Numerical results for stresses and is placements for antisymmetric angle-ply laminates are presented. The three-dimensional elasticity solutions developed are important because they can be used to study the behavior of composite laminates, in addition to serving as reference for approximate solutions by numerical methods and two-dimensional theories.

Savoia, M.↗
NASA NTRS1966

An orthonormalization procedure for multivariable function approximation

Where a function of several variables is given numerically in tabular form, an orthonormalization technique allows an approximation of the numerical data to be determined in a convenient functional form. In this technique, the speed and accuracy of coefficient computation are much improved.

Ingram, H. L.↗
NASA NTRS1989

Recursive linearization of multibody dynamics equations of motion

The equations of motion of a multibody system are nonlinear in nature, and thus pose a difficult problem in linear control design. One approach is to have a first-order approximation through the numerical perturbations at a given configuration, and to design a control law based on the linearized model. Here, a linearized model is generated analytically by following the footsteps of the recursive derivation of the equations of motion. The equations of motion are first written in a Newton-Euler form, which is systematic and easy to construct; then, they are transformed into a relative coordinate representation, which is more efficient in computation. A new computational method for linearization is obtained by applying a series of first-order analytical approximations to the recursive kinematic relationships. The method has proved to be computationally more efficient because of its recursive nature. It has also turned out to be more accurate because of the fact that analytical perturbation circumvents numerical differentiation and other associated numerical operations that may accumulate computational error, thus requiring only analytical operations of matrices and vectors. The power of the proposed linearization algorithm is demonstrated, in comparison to a numerical perturbation method, with a two-link manipulator and a seven degrees of freedom robotic manipulator. Its application to control design is also demonstrated.

Lin, Tsung-Chieh↗
NASA NTRS1989

A model for the formation of the Local Group

Observational tests of a model for the formation of the Local Group are presented and analyzed in which the mass concentration grows by gravitational accretion of local-pressure matter onto two seed masses in an otherwise homogeneous initial mass distribution. The evolution of the mass distribution is studied in an analytic approximation and a numerical computation. The initial seed mass and separation are adjusted to produce the observed present separation and relative velocity of the Andromeda Nebula and the Galaxy. If H(0) is adjusted to about 80 km/s/Mpc with density parameter Omega = 1, then the model gives a good fit to the motions of the outer members of the Local Group. The same model gives particle orbits at radius of about 100 kpc that reasonably approximate the observed distribution of redshifts of the Galactic satellites.

Peebles, P. J. E.↗
NASA NTRS1981

Mixed spectral/finite difference approximations for slightly viscous flows

A numerical scheme is presented, which employs a standard finite difference approximation for the viscous terms in high Reynolds number flows, and resorts to pseudo-spectral methods for the advection terms to greatly reduce the major source of numerical error without dramatically increasing computational cost. The spectral/finite difference (SFD) method evaluates the advection term and second-order differences to evaluate the diffusion term. The fully finite difference (FFD) method with second-order central differences on both terms is also used. The SFD method can handle strong shocks and can outperform the FFD method at moderate viscosity.

Zang, T. A.↗
NASA NTRS1978

Theoretical simulation of solar spectra in the middle ultraviolet and visible for atmospheric trace constituent measurements

Two balloon flights reaching float altitudes of approximately 30 and 40 km respectively, were used to obtain scans of the ultraviolet and visible solar spectra. Both flights covered the UV (2800-3500A) at approximately 0.3A resolution and the visible at approximately 0.6A. Numerous scans were obtained during ascent and from float for both flights. All spectral scans obtained at float, from high sun to low sun, were calibrated in wavelength by using several standard solar spectra for line position references. Comparisons of low sun scans and high sun scans show significant atmospheric continuum extinction and have the potential of being used to identify atmospheric lines superimposed on the attenuated solar spectrum. The resolution was mathematically degraded to approximately 5A to better see the broad band atmospheric extinction. This low resolution is also appropriate for the available low resolution absorption coefficients of NO2 and O3, allowing the identification of NO2 and O3 features on the sunset spectra.

Goldman, A.↗
NASA NTRS1995

An extension of the Lighthill theory of jet noise to encompass refraction and shielding

A formalism for jet noise prediction is derived that includes the refractive 'cone of silence' and other effects; outside the cone it approximates the simple Lighthill format. A key step is deferral of the simplifying assumption of uniform density in the dominant 'source' term. The result is conversion to a convected wave equation retaining the basic Lighthill source term. The main effect is to amend the Lighthill solution to allow for refraction by mean flow gradients, achieved via a frequency-dependent directional factor. A general formula for power spectral density emitted from unit volume is developed as the Lighthill-based value multiplied by a squared 'normalized' Green's function (the directional factor), referred to a stationary point source. The convective motion of the sources, with its powerful amplifying effect, also directional, is already accounted for in the Lighthill format: wave convection and source convection are decoupled. The normalized Green's function appears to be near unity outside the refraction dominated 'cone of silence', this validates our long term practice of using Lighthill-based approaches outside the cone, with extension inside via the Green's function. The function is obtained either experimentally (injected 'point' source) or numerically (computational aeroacoustics). Approximation by unity seems adequate except near the cone and except when there are shrouding jets: in that case the difference from unity quantifies the shielding effect. Further extension yields dipole and monopole source terms (cf. Morfey, Mani, and others) when the mean flow possesses density gradients (e.g., hot jets).

Ribner, Herbert S.↗
NASA NTRS1995

Risk cross sections and their application to risk estimation in the galactic cosmic-ray environment

Radiation risk cross sections (i.e. risks per particle fluence) are discussed in the context of estimating the risk of radiation-induced cancer on long-term space flights from the galactic cosmic radiation outside the confines of the earth's magnetic field. Such quantities are useful for handling effects not seen after low-LET radiation. Since appropriate cross-section functions for cancer induction for each particle species are not yet available, the conventional quality factor is used as an approximation to obtain numerical results for risks of excess cancer mortality. Risks are obtained for seven of the most radiosensitive organs as determined by the ICRP [stomach, colon, lung, bone marrow (BFO), bladder, esophagus and breast], beneath 10 g/cm2 aluminum shielding at solar minimum. Spectra are obtained for excess relative risk for each cancer per LET interval by calculating the average fluence-LET spectrum for the organ and converting to risk by multiplying by a factor proportional to R gamma L Q(L) before integrating over L, the unrestricted LET. Here R gamma is the risk coefficient for low-LET radiation (excess relative mortality per Sv) for the particular organ in question. The total risks of excess cancer mortality obtained are 1.3 and 1.1% to female and male crew, respectively, for a 1-year exposure at solar minimum. Uncertainties in these values are estimated to range between factors of 4 and 15 and are dominated by the biological uncertainties in the risk coefficients for low-LET radiation and in the LET (or energy) dependence of the risk cross sections (as approximated by the quality factor). The direct substitution of appropriate risk cross sections will eventually circumvent entirely the need to calculate, measure or use absorbed dose, equivalent dose and quality factor for such a high-energy charged-particle environment.

NASA Discipline Radiation Health↗
NASA NTRS1989

Numerical integration of asymptotic solutions of ordinary differential equations

Classical asymptotic analysis of ordinary differential equations derives approximate solutions that are numerically stable. However, the analysis also leads to tedious expansions in powers of the relevant parameter for a particular problem. The expansions are replaced with integrals that can be evaluated by numerical integration. The resulting numerical solutions retain the linear independence that is the main advantage of asymptotic solutions. Examples, including the Falkner-Skan equation from laminar boundary layer theory, illustrate the method of asymptotic analysis with numerical integration.

Thurston, Gaylen A.↗
NASA NTRS1972

On propagation of long waves in curved ducts.

A two-dimensional detailed study of the behavior of long waves in curved ducts and in junctions between straight and curved ducts will be given. The mathematical treatment of the problem utilizes the method of separation variables. Solutions and expressions for principal mode of the wave are obtained by using the linearized equation of motion solved for its characteristic values. The unavoidable approximations in the numerical solutions of the cylindrical functions are due to use of series expansion of Bessel functions and from restrictions necessary to solve infinite matrices.

Rostafinski, W.↗
NASA NTRS1972

The transmission of sonic boom signals into rooms through open windows.

This paper is the second of two companion papers. In the first paper, expressions for the acoustic field generated inside a room with an open window due to incoming transient or periodic signals have been presented. In this paper, the technique has been applied to the specific example of an N-wave type signal. Detailed mathematical analysis is followed by its approximate version and numerical computation. Results of controlled experiments, using simulated sonic booms, have been compared with theoretical predictions.

Vaidya, P. G.↗
NASA NTRS1975

Flapping response characteristics of hingeless rotor blades by a gereralized harmonic balance method

Linearized equations of motion for the flapping response of flexible rotor blades in forward flight are derived in terms of generalized coordinates. The equations are solved using a matrix form of the method of linear harmonic balance, yielding response derivatives for each harmonic of the blade deformations and of the hub forces and moments. Numerical results and approximate closed-form expressions for rotor derivatives are used to illustrate the relationships between rotor parameters, modeling assumptions, and rotor response characteristics. Finally, basic hingeless rotor response derivatives are presented in tabular and graphical form for a wide range of configuration parameters and operating conditions.

Peters, D. A.↗
NASA NTRS1975

Detailed gravimetric geoid confirmation of short wavelength features of sea surface topography detected by the Skylab S-193 altimeter in the Atlantic Ocean

A detailed gravimetric geoid was computed for the Northwest Atlantic Ocean and Caribbean Sea area in support of the calibration and evaluation of the GEOS-C altimeter. This geoid, computed on a 15 ft. x 15 ft. grid was based upon a combination of surface gravity data with the GSFC GEM-6 satellite derived gravity data. A comparison of this gravimetric geoid with 10 passes of SKYLAB altimeter data is presented. The agreement of the two data types is quite good with the differences generally less than 2 meters. Sea surface manifestations of numerous short wavelength (approximately 100 km) oceanographic features are now indicated in the gravimetric geoid and are also confirmed by the altimetry data.

Marsh, J. G.↗
NASA NTRS1975

Magnetic attitude control system for dual-spin satellites

A closed-loop control law is developed for a dual-spin satellite control system which utilizes the interaction of the geomagnetic field with the satellite dipole parallel to the spin axis. The control law consists of the linear combination of the pitch axis component of the rate of change of the geomagnetic field and the product of the roll angle and roll axis component of the geomagnetic field. Application of the method of multiple time scales yields approximate solutions for the feedback gains in terms of the system parameters. Approximate solutions are also obtained for the response of the system to disturbance torques. A comparison of the approximate solutions and numerical solutions obtained by numerical integration of the exact equations of motion is then given.

Alfriend, K. T.↗
NASA NTRS1981

Three-dimensional turbulent boundary layer development and separation in V/STOL engine inlets at incidence with small-cross flow and curvature influences

The study concerns the influence of the three-dimensional cross flows on the compressible turbulent boundary layer development and flow separation prediction around V/STOL engine inlets at high incidence. The governing equations for the three-dimensional boundary layer flow with small-cross approximation are solved numerically on an intrinsic streamline coordinate system. Results are presented to illustrate the effects of small cross-flow, compressibility and streamline curvatures on the flow. Comparisons of the results with the wind tunnel data for scaled model and with data obtained from another existing compressible axisymmetric turbulent boundary layer scheme are included in the analysis.

Chou, D. C.↗
NASA NTRS1981

Simulation of charge exchange plasma propagation near an ion thruster propelled spacecraft

A model describing the charge exchange plasma and its propagation is discussed, along with a computer code based on the model. The geometry of an idealized spacecraft having an ion thruster is outlined, with attention given to the assumptions used in modeling the ion beam. Also presented is the distribution function describing charge exchange production. The barometric equation is used in relating the variation in plasma potential to the variation in plasma density. The numerical methods and approximations employed in the calculations are discussed, and comparisons are made between the computer simulation and experimental data. An analytical solution of a simple configuration is also used in verifying the model.

Robinson, R. S.↗
NASA NTRS1982

Numerical simulation of one-dimensional heat transfer in composite bodies with phase change

A numerical simulation was developed to investigate the one dimensional heat transfer occurring in a system composed of a layered aircraft blade having an ice deposit on its surface. The finite difference representation of the heat conduction equations was done using the Crank-Nicolson implicit finite difference formulation. The simulation considers uniform or time dependent heat sources, from heaters which can be either point sources or of finite thickness. For the ice water phase change, a numerical method which approximates the latent heat effect by a large heat capacity over a small temperature interval was applied. The simulation describes the temperature profiles within the various layers of the de-icer pad, as well as the movement of the ice water interface. The simulation could also be used to predict the one dimensional temperature profiles in any composite slab having different boundary conditions.

Dewitt, K. J.↗