Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CALCULATION METHODS”

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 1,009 records · Page 56

LSENS, The NASA Lewis Kinetics and Sensitivity Analysis Code

A general chemical kinetics and sensitivity analysis code for complex, homogeneous, gas-phase reactions is described. The main features of the code, LSENS (the NASA Lewis kinetics and sensitivity analysis code), are its flexibility, efficiency and convenience in treating many different chemical reaction models. The models include: static system; steady, one-dimensional, inviscid flow; incident-shock initiated reaction in a shock tube; and a perfectly stirred reactor. In addition, equilibrium computations can be performed for several assigned states. An implicit numerical integration method (LSODE, the Livermore Solver for Ordinary Differential Equations), which works efficiently for the extremes of very fast and very slow reactions, is used to solve the "stiff" ordinary differential equation systems that arise in chemical kinetics. For static reactions, the code uses the decoupled direct method to calculate sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of dependent variables and/or the rate coefficient parameters. Solution methods for the equilibrium and post-shock conditions and for perfectly stirred reactor problems are either adapted from or based on the procedures built into the NASA code CEA (Chemical Equilibrium and Applications).

Radhakrishnan, K.↗

A Reduced-frequency Approach for Calculating Dynamic Derivatives

Computational Fluid Dynamics (CFD) is increasingly being used to both augment and create an aerodynamic performance database for aircraft configurations. This aerodynamic database contains the response of the aircraft to varying flight conditions and control surface deflections. The current work presents a novel method for calculating dynamic stability derivatives which reduces the computational cost over traditional unsteady CFD approaches by an order of magnitude, while still being applicable to arbitrarily complex geometries over a wide range of flow regimes. The primary thesis of this work is that the response to a forced motion can often be represented with a small, predictable number of frequency components without loss of accuracy. By resolving only those frequencies of interest, the computational effort is significantly reduced so that the routine calculation of dynamic derivatives becomes practical. The current implementation uses this same non-linear, frequency-domain approach and extends the application to the 3-D Euler equations. The current work uses a Cartesian, embedded-boundary method to automate the generation of dynamic stability derivatives.

Murman, Scott M.↗

Stability of matrix equations arising in electromagnetics

Many problems in electromagnetics can be formulated as a Fredholm equation of the first kind. In these problems it is desired to determine some function f(x) from its integral transform g(y). Except for special cases for which an explicit inverse operator is known, the equation is usually discretized by the method of moments into a matrix equation Af = g, which is then solved on a computer. A method of calculating the sensitivity of the solution of Af = g due to errors in either the vector g or the matrix A is developed. The usefulness of this method is demonstrated with problems from antenna theory, scattering theory, and remote sensing.

Klein, C.↗

A discussion of dynamic stability measurement techniques

Techniques for the measurement of the dynamic stability of linear systems are discussed. Particular attention is given to an analysis of the errors in the procedures, and to methods for calculating the system damping from the data. The techniques discussed include: transient decay, moving block analysis, spectral analysis, random decrement signatures, transfer function analysis, and parameter identification methods. The special problems of rotorcraft dynamic stability testing are discussed.

Johnson, W.↗

Computation of wave drag for transonic flow

The paper develops a method for calculating wave drag for two-dimensional transonic flow, with particular application to the prediction of the drag rise Mach number of a supercritical wing section. The method is based on a transonic similarity model which is defined by a normalized small perturbation equation and represents shock waves by the addition of an artificial viscosity term in the region of supersonic flow to the partial differential equation. The drag formula obtained allows the computer simulation of transonic wind tunnel data taking account of boundary layer and wall effects.

Garabedian, P. R.↗

Calculation of the planar supercritical flow over a NASA supercritical profile

An unsteady finite difference procedure was used to calculate the steady inviscid flow over an 11% thickness ratio NASA supercritical profile of LWP 505 at M sub infinity = 0.80 and alpha = 0. An attempt is made to include the viscous effects using a modified form of Head's entrainment method to calculate the turbulent boundary layer. The attempt to predict the viscous effects using the compressible form of Head's integral method with a modified auxiliary equation for the form factors was unsatisfactory. Though a reasonably separation bubble was established on the lower surface, a grossly exaggerated displacement effect resulted downstream of the shock on the upper surface. There clearly is substantial further effort required to evolve a satisfactory boundary layer procedure, which must then be coupled in a still unproven manner with the inviscid procedure.

Yoshihara, H.↗

Raman determination of layer stresses and strains for heterostructures and its application to the cubic SiC/Si system

A set of formulas for a generalized axial stress in diamond and zinc-blende semiconductors under axial stress is derived to calculate stress-related Raman shifts. By analyzing known Raman data on cubic SiC under hydrostatic pressures, one of the Raman-stress coefficients was obtained, and the existing elastic stiffness constants of cubic SiC were optimized. A method for calculating the stress and strain in SiC films on (100) Si is proposed. It is suggested that the stress and strain expressions and the method of the stress and strain determinations in heterostructures are quite general and may be used for other systems.

Feng, Z. C.↗

Regional surface fluxes from satellite-derived surface temperatures (AVHRR) and radiosonde profiles

Radiometric surface temperatures, derived from measurements by the AVHRR instrument aboard the NOAA-9 and the NOAA-11 polar orbiting satellites, were used in combination with wind velocity and temperature profiles measured by radiosondes, to calculate surface fluxes of sensible heat. The measurements were made during FIFE, the First ISLSCP (International Satellite Land Surface Climatology Project) Field Experiment, in a hilly tall grass prairie area of northeastern Kansas. The method of calculation was based on turbulent similarity formulations for the atmospheric boundary layer. Good agreement (r = 0.7) was obtained with reference values of sensible heat flux, taken as arithmetric means of measurements with the Bowen ratio method at six ground stations. The values of evaporation (latent heat fluxes), derived from these sensible heat fluxes by means of the energy budget, were also in good agreement (r = 0.94) with the corresponding reference values from the ground stations.

Brutsaert, Wilfried↗

Mn-53 and Al-26 evidence for solar cosmic ray constancy - An improved model for interpretation

A general, numerical method for calculating activities of solar cosmic ray produced radionuclides at any point within an irregularly shaped lunar rock of known surface contour and lunar surface orientation is described. This method is then used to predict the activities of Mn-53 andAl-26 as a function of postion within lunar rock 68815 for various assumed values of solar cosmic ray flux (J), rigidity (R sub 0), and rock erosion rate (ER). The predicted activities agree with the measured activities of Kohl et al. (1978) when values of R sub 0 = 100 MV, J = 70 p/sq cm-sec (4 pi, E greater than 10 MeV), ER not greater than 1 mm/m.y. and a total exposure time of 2 m.y. are assumed. These values are in agreement with those found for rocks exposed for not less than 10 m.y. and provide no evidence for variation of the average solar cosmic ray parameters between the last 2 and 10 m.y. intervals. When interpreted with the improved model the activity vs. depth profiles for three faces of 68815 show no evidence of SCR anisotropy or differential erosion.

Russ, G. P., III↗

Advances in computational structural and solid mechanics; Proceedings of the First World Congress on Computational Mechanics, University of Texas, Austin, Sept. 22-26, 1986

Papers are presented on the nonlinear dynamic analysis of quasi-symmetric anisotropic structures, the analysis of large deformations of membrane shells by the generalized finite difference method, the vibration of shear-deformable laminated plate structures by the finite strip method, and Coon's surface method for the formulation of finite elements for plates and shells. Other topics include an evaluation of higher-order modal methods for calculating the transient structural response, a nonlinear dynamic analysis of frame structures, and a nonlinear analysis of a truss by energy minimization. Also considered are an incremental Galerkin method for plates and stiffened plates, a comparison of the postbuckling behavior of plates and shells, and a general solution of bending in a cylindrical shell.

Noor, Ahmed K.↗

On the integration of image sources in exact image method of field analysis

Convergence conditions of image integration in the exact image method of field calculation have been investigated, and the method is extended to include more general media than previously considered. It is demonstrated that the integral is well behaved and the method works best when the field is calculated in a medium with less loss. If the medium has more loss, it is shown that the image line may enter the improper half-space and still produce valid results. Means of correctly selecting the integration path branch for the case of crossing branch cuts of the Green function in complex integration planes are proposed.

Lindell, I. V.↗

Discrete-vortex model for the symmetric-vortex flow on cones

A relatively simple but accurate potential flow model was developed for studying the symmetric vortex flow on cones. The model is a modified version of the model first developed by Bryson, in which discrete vortices and straight-line feeding sheets were used to represent the flow field. It differs, however, in the zero-force condition used to position the vortices and determine their circulation strengths. The Bryson model imposed the condition that the net force on the feeding sheets and discrete vortices must be zero. The proposed model satisfies this zero-force condition by having the vortices move as free vortices, at a velocity equal to at the local crossflow velocity at their centers. When the free-vortex assumption is made, a solution is obtained in the form of two nonlinear algebraic equations that relate the vortex center coordinates and vortex strengths to the cone angle and angle of attack. The vortex center locations calculated using the model are in good agreement with experimental values. The cone normal forces as well as center locations are in good agreement with the vortex cloud method of calculating symmetric flow fields.

Gainer, Thomas G.↗

An accelerated subspace iteration for eigenvector derivatives

An accelerated subspace iteration method for calculating eigenvector derivatives has been developed. Factors affecting the effectiveness and the reliability of the subspace iteration are identified, and effective strategies concerning these factors are presented. The method has been implemented, and the results of a demonstration problem are presented.

Ting, Tienko↗

Radiation transfer in gamma-ray bursts

We have calculated gamma-ray radiative transport in regions of high-energy density, such as gamma-ray burst source regions, using a discrete ordinate, discrete energy group method. The calculations include two-photon pair production and annihilation, as well as three-photon annihilation. The radiation fields itself acts as an absorbing medium, and the optical depth depends on its intensity, so the problem is intrinsically nonlinear. Spherical divergence produces effective collimation of the flux. At high optical depth the high energy (E is greater than 1 MeV) portion of the emergent spectrum assumes a nearly universal form. An approximate limit is derived for the high-energy flux from a gamma-ray burst source region of given size, and the implications of this limit for the distance to the 1979 March 5 event are briefly discussed. We discuss more generally the problem of very luminous bursts, and implications of Galactic halo distances for flare models.

Carrigan, B. J.↗

Predicting the Rotor-Stator Interaction Acoustics of a Ducted Fan Engine

A Navier-Stokes computation is performed for a ducted-fan configuration with the goal of predicting rotor-stator noise generation without having to resort to heuristic modeling. The calculated pressure field in the inlet region is decomposed into classical infinite-duct modes, which are then used in either a hybrid finite-element/Kirchhoff surface method or boundary integral equation method to calculate the far field noise. Comparisons with experimental data are presented, including rotor wake surveys and far field sound pressure levels for two blade passage frequency (BPF) tones.

Biedron, Robert T.↗

Predicting the Rotor-Stator Interaction Acoustics of a Ducted Fan Engine

A Navier-Stokes computation is performed for a ducted-fan configuration with the goal of predicting rotor-stator noise generation without having to re- to heuristic modeling. The calculated pressure field in the inlet region is decomposed into classical infinite-duct, modes, which are then used in either a hybrid finite-element /Kirchhoff surface method or boundary integral equation method to calculate the far field noise. Comparisons with experimental data are presented, including rotor wake surveys and far field sound pressure levels for 2 blade passage frequency (BPF) tones.

Biedron, Robert T.↗

Compatibility Condition in Theory of Solid Mechanics (Elasticity, Structures, and Design Optimization)

The strain formulation in elasticity and the compatibility condition in structural mechanics have neither been understood nor have they been utilized. This shortcoming prevented the formulation of a direct method to calculate stress. We have researched and understood the compatibility condition for linear problems in elasticity and in finite element analysis. This has lead to the completion of the method of force with stress (or stress resultant) as the primary unknown. The method in elasticity is referred to as the completed Beltrami-Michell formulation (CBMF), and it is the integrated force method (IFM) in structures. The dual integrated force method (IFMD) with displacement as the primary unknown has been formulated. IFM and IFMD produce identical responses. The variational derivation of the CBMF yielded the new boundary compatibility conditions. The CBMF can be used to solve stress, displacement, and mixed boundary value problems. The IFM in structures produced high-fidelity response even with a modest finite element model. The IFM has influenced structural design considerably. A fully utilized design method for strength and stiffness limitation has been developed. The singularity condition in optimization has been identified. The CBMF and IFM tensorial approaches are robust formulations because of simultaneous emphasis on the equilibrium equation and the compatibility condition.

Patnaik, Surya N.↗