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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Application of a monotonic upstream-biased transport scheme to three-dimensional constituent transport calculations

The application of van Leer's scheme, a monotonic, upstream-biased differencing scheme, to three-dimensional constituent transport calculations is shown. The major disadvantage of the scheme is shown to be a self-limiting diffusion. A major advantage of the scheme is shown to be its ability to maintain constituent correlations. The scheme is adapted for a spherical coordinate system with a hybrid sigma-pressure coordinate in the vertical. Special consideration is given to cross-polar flow. The vertical wind calculation is shown to be extremely sensitive to the method of calculating the divergence. This sensitivity implies that a vertical wind formulation consistent with the transport scheme is essential for accurate transport calculations. The computational savings of the time-splitting method used to solve this equation are shown. Finally, the capabilities of this scheme are illustrated by an ozone transport and chemistry model simulation.

Allen, Dale J.↗
NASA NTRS2012

Calculation of Sensitivity Derivatives in an MDAO Framework

During gradient-based optimization of a system, it is necessary to generate the derivatives of each objective and constraint with respect to each design parameter. If the system is multidisciplinary, it may consist of a set of smaller "components" with some arbitrary data interconnection and process work ow. Analytical derivatives in these components can be used to improve the speed and accuracy of the derivative calculation over a purely numerical calculation; however, a multidisciplinary system may include both components for which derivatives are available and components for which they are not. Three methods to calculate the sensitivity of a mixed multidisciplinary system are presented: the finite difference method, where the derivatives are calculated numerically; the chain rule method, where the derivatives are successively cascaded along the system's network graph; and the analytic method, where the derivatives come from the solution of a linear system of equations. Some improvements to these methods, to accommodate mixed multidisciplinary systems, are also presented; in particular, a new method is introduced to allow existing derivatives to be used inside of finite difference. All three methods are implemented and demonstrated in the open-source MDAO framework OpenMDAO. It was found that there are advantages to each of them depending on the system being solved.

Moore, Kenneth T.↗
NASA NTRS1995

Calculating 3-D Crack-Opening Behavior Of A Fatigue Crack

Method for calculating through-the-thickness, three-dimensional (3-D) crack-opening behavior in standard fatigue-crack-growth test specimens devised. Involves determination of crack-opening loads from crack-propagation data and finite-element analyses. In method, one considers variation in rate of crack growth along crack front, relationship between rate of crack growth and effective stress intensity factor, and variation in three-dimensional stress-intensity factor along crack front. Developed to provide simpler means for obtaining fatigue-crack-opening behavior loads.

Newman, J. C., Jr.↗
NASA NTRS1975

Calculation of arbitrary-order diffraction efficiencies of thick gratings with arbitrary grating shape

A method for calculating arbitrary-order diffraction efficiencies of thick, lossless transmission gratings with arbitrary periodic grating shapes has been developed. A Fourier-series representation of the grating is employed, along with a coupled-mode theory of diffraction. For illustration, numerical values of the diffraction efficiencies at the first three Bragg angles are calculated for sinusoidal, square-wave, triangular, and saw-tooth gratings. Numerical results for the same grating shapes with the same parameters are also calculated for comparison, by extending Burckhardt's numerical method for analyzing thick sinusoidal gratings. The comparison shows that the coupled-mode theory provides results with relative computational ease and results that are in agreement with calculations obtained by extending the more-rigorous Burckhardt theory to nonsinusoidal grating shapes and to higher-order Bragg angles.

Su, S. F.↗
NASA NTRS1990

The C2 Ballik-Ramsay transitions in carbon stars

The bands of C2 are an important source of opacity in the near-infrared spectral region of carbon stars. In the 2 micron spectral region, the Delta v = -1 and -2 transitions of the Ballik-Ramsay b 3Sigma(g)(-) - a 3Pi(u) system are especially important. A general method for calculating highly accurate transition parameters of the Ballik-Ramsay band is presented. This method allows the easy calculation of the satellite lines and the forbidden lines resulting from the perturbation of the upper b 3Sigma(g)(-) state by the X 1Sigma(0)(+) state. Included in the calculation is the effect of the isotopic molecules (C-12)2, (C-12)(C-13), and (C-13)2. Three consistency tests of the calculated strengths are presented. A comparison of synthetic spectra is made with an observed spectrum of WZ Cas in the 3940/cm spectral region.

Goorvitch, David↗
NASA NTRS2010

Lightning Reporting at 45th Weather Squadron: Recent Improvements

The 45th Weather Squadron (45 WS) provides daily lightning reports to space launch customers at CCAFS/KSC. These reports are provided to assess the need to inspect the electronics of satellite payloads, space launch vehicles, and ground support equipment for induced current damage from nearby lightning strokes. The 45 WS has made several improvements to the lightning reports during 2008-2009. The 4DLSS, implemented in April 2008, provides all lightning strokes as opposed to just one stroke per flash as done by the previous system. The 45 WS discovered that the peak current was being truncated to the nearest kilo amp in the database used to generate the daily lightning reports, which led to an up to 4% underestimate in the peak current for average lightning. This error was corrected and led to elimination of this underestimate. The 45 WS and their mission partners developed lightning location error ellipses for 99% and 95% location accuracies tailored to each individual stroke and began providing them in the spring of 2009. The new procedure provides the distance from the point of interest to the best location of the stroke (the center of the error ellipse) and the distance to the closest edge of the ellipse. This information is now included in the lightning reports, along with the peak current of the stroke. The initial method of calculating the error ellipses could only be used during normal duty hours, i.e. not during nights, weekends, or holidays. This method was improved later to provide lightning reports in near real-time, 24/7. The calculation of the distance to the closest point on the ellipse was also significantly improved later. Other improvements were also implemented. A new method to calculate the probability of any nearby lightning stroke. being within any radius of any point of interest was developed and is being implemented. This may supersede the use of location error ellipses. The 45 WS is pursuing adding data from nine NLDN sensors into 4DLSS in real-time. This will overcome the problem of 4DLSS missing some of the strong local strokes. This will also improve the location accuracy, reduce the size and eccentricity of the location error ellipses, and reduce the probability of nearby strokes being inside the areas of interest when few of the 4DLSS sensors are used in the stroke solution. This will not reduce 4DLSS performance when most of the 4DLSS sensors are used in the stroke solution. Finally, several possible future improvements were discussed, especially for improving the peak current estimate and the error estimate for peak current, and upgrading the 4DLSS. Some possible approaches for both of these goals were discussed.

Finn, Frank C.↗
NASA NTRS1983

Homogeneous crystal nucleation in binary metallic melts

A method for calculating the homogeneous crystal nucleation frequency in binary metallic melts is developed. The free energy of crystallization is derived from regular solution models for the liquid and solid and is used, together with model-based estimates of the interfacial tension, to calculate the nucleation frequency from the classical theory. The method can account for the composition dependence of the maximum undercooling observed in a number of experiments on small droplet dispersions. It can also be used to calculate the driving force for crystal growth and to obtain more precise estimates of the homogeneous crystal nucleation frequency in glass-forming alloys. This method, although approximate, is simple to apply, and requires only knowledge of the phase diagram and a few readily available thermodynamic quantities as input data.

Thompson, C. V.↗
NASA NTRS1983

Calculation of symmetric vortex separation affecting subsonic bodies at high incidence

A method for calculating the longitudinal aerodynamic coefficients and the pressure distributions on a body at reasonably high angles of attack is presented. The body is represented by a combination of source elements and vortex-lattice elements, including separation of the vortices at increasing angles of attack. The method is self-consistent in that the body and the separated vortex wake are treated as an integrated interacting system. The location of the separation line can be included as an arbitrary input from experimental data or can be evaluated approximately by a pressure-dependent criterion. Calculated values of the aerodynamic coefficients and pressure distributions on cone-cylinder and ogive-cylinder bodies compare well, qualitatively and quantitatively, with experimental data, including simulation of the dependence on Reynolds number.

Almosnino, D.↗
NASA NTRS2007

Asymptotic Energies and QED Shifts for Rydberg States of Helium

This paper reviews progress that has been made in obtaining essentially exact solutions to the nonrelativistic three-body problem for helium by a combination of variational and asymptotic expansion methods. The calculation of relativistic and quantum electrodynamic corrections by perturbation theory is discussed, and in particular, methods for the accurate calculation of the Bethe logarithm part of the electron self energy are presented. As an example, the results are applied to the calculation of isotope shifts for the short-lived 'halo' nucleus He-6 relative to He-4 in order to determine the nuclear charge radius of He-6 from high precision spectroscopic measurements carried out at the Argonne National Laboratory. The results demonstrate that the high precision that is now available from atomic theory is creating new opportunities to create novel measurement tools, and helium, along with hydrogen, can be regarded as a fundamental atomic system whose spectrum is well understood for all practical purposes.

Drake, G.W.F.↗
NASA NTRS1958

Comparison of Calculated and Experimental Temperatures and Coolant Pressure Losses for a Cascade of Small Air-Cooled Turbine Rotor Blades

Average spanwise blade temperatures and cooling-air pressure losses through a small (1.4-in, span, 0.7-in, chord) air-cooled turbine blade were calculated and are compared with experimental nonrotating cascade data. Two methods of calculating the blade spanwise metal temperature distributions are presented. The method which considered the effect of the length-to-diameter ratio of the coolant passage on the blade-to-coolant heat-transfer coefficient and assumed constant coolant properties based on the coolant bulk temperature gave the best agreement with experimental data. The agreement obtained was within 3 percent at the midspan and tip regions of the blade. At the root region of the blade, the agreement was within 3 percent for coolant flows within the turbulent flow regime and within 10 percent for coolant flows in the laminar regime. The calculated and measured cooling-air pressure losses through the blade agreed within 5 percent. Calculated spanwise blade temperatures for assumed turboprop engine operating conditions of 2000 F turbine-inlet gas temperature and flight conditions of 300 knots at a 30,000-foot altitude agreed well with those obtained by the extrapolation of correlated experimental data of a static cascade investigation of these blades.

Stepka, Francis S↗
NASA NTRS1993

Prediction of nonlinear optical properties of organic materials. General theoretical considerations

The prediction of nonlinear optical properties of organic materials is geared to assist materials scientists in the selection of good candidate molecules. A brief summary of the quantum mechanical methods used for estimating hyperpolarizabilities will be presented. The advantages and limitations of each technique will be discussed. Particular attention will be given to the finite-field method for calculating first and second order hyperpolarizabilities, since this method is better suited for large molecules. Corrections for dynamic fields and bulk effects will be discussed in detail, focusing on solvent effects, conformational isomerization, core effects, dispersion, and hydrogen bonding. Several results will be compared with data obtained from third-harmonic-generation (THG) and dc-induced second harmonic generation (EFISH) measurements. These comparisons will demonstrate the qualitative ability of the method to predict the relative strengths of hyperpolarizabilities of a class of compounds. The future application of molecular mechanics, as well as other techniques, in the study of bulk properties and solid state defects will be addressed. The relationship between large values for nonlinear optical properties and large conjugation lengths is well known, and is particularly important for third-order processes. For this reason, the materials with the largest observed nonresonant third-order properties are conjugated polymers. An example of this type of polymer is polydiacetylene. One of the problems in dealing with polydiacetylene is that substituents which may enhance its nonlinear properties may ultimately prevent it from polymerizing. A model which attempts to predict the likelihood of solid-state polymerization is considered, along with the implications of the assumptions that are used. Calculations of the third-order optical properties and their relationship to first-order properties and energy gaps will be discussed. The relationship between monomeric and polymeric third-order optical properties will also be considered.

Cardelino, B.↗
NASA NTRS1940

The frequencies of cantilever wings in beam and torsional vibrations

Methods are described for calculating the period and frequency of vibration of cantilever wings and similar structures in which the weight and moment of inertia vary along the span. Both the beam and torsional frequencies may be calculated by these methods. The procedure is illustrated by examples. It is shown that a surprisingly close approximation to the beam frequency may be obtained by a very brief calculation in which the curvature of the wing in vibration is assumed to be constant. A somewhat longer computation permits taking account of the true curvature of the beam by a series of successive approximations which are shown to be strongly convergent. Analogous methods are applied to calculations of the torsional frequency. For the first approximation it is assumed that the angle of twist varies linearly alone the semispan. True variation of the twist is computed by successive approximations which are strongly convergent, as in the case of beam vibrations.

Burgess, C P↗
NASA NTRS1972

A method of limit point calculation in finite element structural analysis

An approach is presented for the calculation of limit points for structures described by discrete coordinates, and whose governing equations derive from finite element concepts. The nonlinear load-displacement path of the imperfect structure is first traced by use of a direct iteration scheme and the determinant of the governing algebraic equations is calculated at each solution point. The limit point is then established by extrapolation and imposition of the condition of zero slope of the plot of load vs. determinant. Three problems are solved in illustration of the approach and in comparison with alternative procedures and test data.

Gallagher, R. H.↗
NASA NTRS1973

Iterative methods for plasma sheath calculations: Application to spherical probe

The computer cost of a Poisson-Vlasov iteration procedure for the numerical solution of a steady-state collisionless plasma-sheath problem depends on: (1) the nature of the chosen iterative algorithm, (2) the position of the outer boundary of the grid, and (3) the nature of the boundary condition applied to simulate a condition at infinity (as in three-dimensional probe or satellite-wake problems). Two iterative algorithms, in conjunction with three types of boundary conditions, are analyzed theoretically and applied to the computation of current-voltage characteristics of a spherical electrostatic probe. The first algorithm was commonly used by physicists, and its computer costs depend primarily on the boundary conditions and are only slightly affected by the mesh interval. The second algorithm is not commonly used, and its costs depend primarily on the mesh interval and slightly on the boundary conditions.

Parker, L. W.↗
NASA NTRS1975

Studies of scramjet/airframe integration techniques for hypersonic aircraft

New design and analysis techniques for engine-airframe integration were applied in a recent hypersonic vehicle design study. A new technique was developed to design the vehicle's forebody so that uniform precompressed flow was produced at the inlet entrance. Results are verified with three-dimensional characteristic calculations. Results from a new three-dimensional method for calculating nozzle flows show that the entire lower afterbody of the vehicle can be used as a scramjet exhaust nozzle to achieve efficient, controlled, and stable flight over a wide range of flight conditions.

Edwards, C. L. W.↗