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 883 records · Page 49

Method and models for R-curve instability calculations

This paper presents a simple method for performing elastic R-curve instability calculations. For a single material-structure combination, the calculations can be done on some pocket calculators. On microcomputers and larger, it permits the development of a comprehensive program having libraries of driving force equations for different configurations and R-curve model equations for different materials. The paper also presents several model equations for fitting to experimental R-curve data, both linear elastic and elastoplastic. The models are fit to data from the literature to demonstrate their viability.

Orange, Thomas W.↗

Calculating Flows With Interfering Shock Waves

Equilibrium Air Shock Interference, EASI, program takes account of dissociation of air molecules. Revives and updates older computational methods for calculating inviscid flow field and maximum heating from interference of shock waves. Expands methods to solve problems involving six shock-wave interference patterns on two-dimensional cylindrical leading edge with equilibrium, chemically-reacting gas mixture. Written in FORTRAN 77.

Glass, Christopher E.↗

Implementation of equivalent domain integral method in the two-dimensional analysis of mixed mode problems

An equivalent domain integral (EDI) method for calculating J-intergrals for two-dimensional cracked elastic bodies is presented. The details of the method and its implementation are presented for isoparametric elements. The total and product integrals consist of the sum of an area of domain integral and line integrals on the crack faces. The line integrals vanish only when the crack faces are traction free and the loading is either pure mode 1 or pure mode 2 or a combination of both with only the square-root singular term in the stress field. The EDI method gave accurate values of the J-integrals for two mode I and two mixed mode problems. Numerical studies showed that domains consisting of one layer of elements are sufficient to obtain accurate J-integral values. Two procedures for separating the individual modes from the domain integrals are presented. The procedure that uses the symmetric and antisymmetric components of the stress and displacement fields to calculate the individual modes gave accurate values of the integrals for all problems analyzed. The EDI method when applied to a problem of an interface crack in two different materials showed that the mode 1 and mode 2 components are domain dependent while the total integral is not. This behavior is caused by the presence of the oscillatory part of the singularity in bimaterial crack problems. The EDI method, thus, shows behavior similar to the virtual crack closure method for bimaterial problems.

Raju, I. S.↗

Charts for determining the pitching moment of tapered wings with sweepback and twist

This report presents a convenient method for calculating the pitching-moment characteristics of tapered wings with sweepback and twist. The method is based on the fact that the pitching-moment characteristics of a wing may be specified by giving the value of the pitching moment at zero lift and the location of the axis about which the axis is constant. Data for calculating these characteristics are presented by curves which apply to wings having a linear distribution of twist along the span and which cover a large range of aspect ratios. The curves are given for wings having straight taper and distorted elliptical plan forms. The characteristics of wings of other shapes may be determined by interpolation.

Anderson, Raymond F↗

Signal Quality Influence on Averaging Procedure for DIAL Pollution Monitoring

Presented here is an analysis of the relative effect of each stochastic phenomenon on the accuracy of LIDAR measurements. The preferred averaging procedure to be applied is noted. The analysis takes into account the temporal correlation of the transmitted power variations which are associated with laser modal noise, short term temporal correlation of atmospheric turbulence, and 'white' additive noise. A distinction is made between two types of systems, one that transmits a sequence of pulses tuned into two wavelengths, but with a relatively large time separation, and another that transmits pulses in two wavelengths almost coincidentally. This distinction is important because the effect of speckle noise can be minimized when the pulses at two wavelengths are sent within the correlation time of the atmospheric turbulence. The presented analysis compares the signal to noise ratio in the calculation of the ratio beta (sub 1)/beta (sub 2), as achieved by using two methods of calculations: first averaging the data related to each wavelength and then taking the ratio, and first taking the ratio and then averaging. For both methods, normalization with respect to the transmitted power is considered, either prior to or after the averaging of the signals in each wavelength.

Egert, S.↗

Aerodynamic Loads on Tails at High Angles of Attack and Sideslip

Results are presented for the loads and moments acting on the individual tail surfaces of a body-tail combination over a wide range of angles of attack and sideslip. The effects of forebody length and panel-panel interference on the characteristics are included. It is shown that large nonlinear variations in these loads and moments, which occur at some combinations of angle of attack and sideslip, cannot be predicted by low-angle theory. A relatively simple, but general, theoretical method for calculating these load and moment characteristics is described, and the results from this method are found to be in good agreement with experiment provided the initial positions of the forebody vortices are known. It is shown that a simple application of slender-body theory can be used to predict the side loads due to sideslip that are contributed by a vertical tail on a wide variety of wing-body-tail combinations at low angles of attack. For several configurations, changes are indicated which reduced the vertical-tail loads per unit yawing moment of each complete configuration at large angles of attack. Some results are presented on the effect of high angle of attack on the induced-flow field and tail loads due to a wing at supersonic speed.

SLENDER BODY↗

Collisional redistribution of radiation. I - The density matrix

As a first step in deriving accurate, calculable expressions for collisional redistribution of radiation, valid outside the impact regime a method for calculating an accurate binary-collision density operator for an atom in the presence of a driving field and perturbers is presented. Projection-operator techniques to establish first the precise validity of the Markoff approximation for radiative relaxation are used. Using the same techniques, it is shown how the collisional-relaxation problem may be analyzed in a practical manner outside the Markoff approximation. The effect of correlations between radiative and collisional events is included in a consistent way and the physical implications of these correlation effects are demonstrated for a simple example, the two-level atom with nondegenerate levels.

Burnett, K.↗

Loran-C plotting program for plotting lines of position on standard charts

A set of programs designed to be run on the IBM System/370 computer is described. These programs are used to plot Loran-C lines of position (LOP) on any common map or standard aviation sectional chart. The Loran-C plotting program JRPLOT FORTRAN uses a standard Calcomp-compatible plotting subroutine package for the Hewlett-Packard 7203A graphic plotter. A general description of the features of the Loran-C plotting program. This program involves a simple add/substrate method to calculate the LOP. A description on how to use the program and some methods of operation are included.

Roman, J. P.↗

Pole of rotating analysis of present-day Juan de Fuca plate motion

Convergence rates between the Juan de Fuca and North American plates are calculated by means of their relative, present-day pole of rotation. A method of calculating the propagation of errors in addition to the instantaneous poles of rotation is also formulated and applied to determine the Euler pole for Pacific-Juan de Fuca. This pole is vectorially added to previously published poles for North America-Pacific and 'hot spot'-Pacific to obtain North America-Juan de Fuca and 'hot spot'-Juan de Fuca, respectively. The errors associated with these resultant poles are determined by propagating the errors of the two summed angular velocity vectors. Under the assumption that hot spots are fixed with respect to a mantle reference frame, the average absolute velocity of the Juan de Puca plate is computed at approximately 15 mm/yr, thereby making it the slowest-moving of the oceanic plates.

Nishimura, C.↗

Fatigue life estimates for helicopter loading spectra

Helicopter loading histories applied to notch metal samples are used as examples, and their fatigue lives are calculated by using a simplified version of the local strain approach. This simplified method has the advantage that it requires knowing the loading history in only the reduced form of ranges and means and number of cycles from the rain-flow cycle counting method. The calculated lives compare favorably with test data.

Khosrovaneh, A. K.↗

Fatigue life estimates for helicopter loading spectra

Helicopter loading histories applied to notch metal samples are used as examples, and their fatigue lives are calculated by using a simplified version of the local strain approach. This simplified method has the advantage that it requires knowing the loading history in only the reduced form of ranges and means and number of cycles from the rain-flow cycle counting method. The calculated lives compare favorably with test data.

Khosrovaneh, A. K.↗

Wall catalytic recombination and boundary conditions in nonequilibrium hypersonic flows - With applications

The meaning of catalysis and its relation to aerodynamic heating in nonequilibrium hypersonic flows are discussed. The species equations are described and boundary conditions for them are derived for a multicomponent gas and for a binary gas. Slip effects are included for application of continuum methods to low-density flows. Measurement techniques for determining catalytic wall recombination rates are discussed. Among them are experiments carried out in arc jets as well as flow reactors. Diagnostic methods for determining the atom or molecule concentrations in the flow are included. Results are given for a number of materials of interest to the aerospace community, including glassy coatings such as the RCG coating of the Space Shuttle and for high temperature refractory metals such as coated niobium. Methods of calculating the heat flux to space vehicles in nonequilibrium flows are described. These methods are applied to the Space Shuttle, the planned Aeroassist Flight Experiment, and a hypersonic slender vehicle such as a transatmospheric vehicle.

Scott, Carl D.↗

Characterization on defect solids

The purpose of this research has been to develop new semi-empirical techniques to describe ceramics, alloys, and metal/ceramic interfaces with applications in mind that support the materials aspects of the high speed civil transport program (HSCT). HSCT requires methods that aid in the design of alloys and ceramics for new high strength, high temperature applications. Current theoretical methods are not capable of carrying out this mission. Hence, new accurate and more efficient theoretical techniques are needed to facilitate the design of new materials and the examination of their properties. This program concentrated on modeling ceramics, but also dealt with alloy properties to a lesser degree. The primary accomplishment of this research was the development of a new equation of state (EOS) that models the energy/bond length relationship and includes terms that represent charge transfer between cations and anions. The new EOS has been used by researchers at the Naval Research Laboratory to develop a new semi-empirical method for calculating properties of ceramics and metal ceramic interfaces based on the embedded atom method. The results of the director's discretionary fund (DDF) research describing the binding energy relation was the enabling theoretical basis. Also the workers at Cleveland State and NASA Lewis have detailed an approach based on equivalent crystal theory as part of the DDF which is currently under development. In addition, initial results for using the Harris functional method have been developed. Finally, the techniques developed by us have proved useful for the study of high pressure properties of solids and have been used extensively by researchers in this field.

Schlosser, Herbert↗

Metafitting: Weight optimization for least-squares fitting of PTTI data

For precise time intercomparisons between a master frequency standard and a slave time scale, we have found it useful to quantitatively compare different fitting strategies by examining the standard uncertainty in time or average frequency. It is particularly useful when designing procedures which use intermittent intercomparisons, with some parameterized fit used to interpolate or extrapolate from the calibrating intercomparisons. We use the term 'metafitting' for the choices that are made before a fitting procedure is operationally adopted. We present methods for calculating the standard uncertainty for general, weighted least-squares fits and a method for optimizing these weights for a general noise model suitable for many PTTI applications. We present the results of the metafitting of procedures for the use of a regular schedule of (hypothetical) high-accuracy frequency calibration of a maser time scale. We have identified a cumulative series of improvements that give a significant reduction of the expected standard uncertainty, compared to the simplest procedure of resetting the maser synthesizer after each calibration. The metafitting improvements presented include the optimum choice of weights for the calibration runs, optimized over a period of a week or 10 days.

Douglas, Rob J.↗

Finite-field Calculations of Molecular Polarizabilities Using Field-Induced Polarization Functions: Second- and Third-order Perturbation Correlation Corrections to the Coupled Hartree-Fock Polarizability of H2O

Ordinary Rayleigh-Schroudinger perturbation theory with Moller-Plesset (RSMP) partitioning is used to calculate second- and third-order correlation corrections to the CHF polarizability and dipole moment of the water molecule by a finite-field procedure. Pade approximants are found to be useful in accelerating the convergence of the property perturbation expansions. Field-induced polarization functions suitable for polarizability calculations are determined. The average polarizability calculated, neglecting vibrational averaging, with Dunning's (9s5p/4s-4s2p/2s) contracted GTO basis set augmented by field-induced lslp2d/lp polarization functions is within 3 per cent of the experimental result. Correlation corrections to the dipole moment and polarizability of the water molecule calculated by the finite-field RSMP and single + double excitation CI(SDCI) methods for the same basis set are found to be in close agreement. The RSMP approach has the advantages of being size-consistent and of being capable of greater efficiency than the SCDI method. Comparative calculations carried out using Epstein-Nesbet partitioning show that through third order RSEN correlation perturbation expansions for the dipole moment and polarizability are less rapidly convergent than RSMP expansions. However, reasonable accord with RSMP results can be achieved by using Pade approximants to accelerate the convergence of RSEN energy perturbation expansions. The convergence of RSMP property correlation expansions based on the zeroth-order uncoupled-Hartree-Fock (UCHF) and coupled-Hartree-Fock (CHF) approximations are compared through third order. Whereas the CHF + RSMP expansions are for practical purposes fully converged, the UCHF + RSMP expansions are not adequately converged.

Langhoff, S. R.↗

Comparison of Artificial Compressibility Methods

Various artificial compressibility methods for calculating the three-dimensional incompressible Navier-Stokes equations are compared. Each method is described and numerical solutions to test problems are conducted. A comparison based on convergence behavior, accuracy, and robustness is given.

Kiris, Cetin↗

Prediction Interval Development for Wind-Tunnel Balance Check-Loading

Results from the Facility Analysis Verification and Operational Reliability project revealed a critical gap in capability in ground-based aeronautics research applications. Without a standardized process for check-loading the wind-tunnel balance or the model system, the quality of the aerodynamic force data collected varied significantly between facilities. A prediction interval is required in order to confirm a check-loading. The prediction interval provides an expected upper and lower bound on balance load prediction at a given confidence level. A method has been developed which accounts for sources of variability due to calibration and check-load application. The prediction interval method of calculation and a case study demonstrating its use is provided. Validation of the methods is demonstrated for the case study based on the probability of capture of confirmation points.

Landman, Drew↗

Calculating Free Energies Using Average Force

A new, general formula that connects the derivatives of the free energy along the selected, generalized coordinates of the system with the instantaneous force acting on these coordinates is derived. The instantaneous force is defined as the force acting on the coordinate of interest so that when it is subtracted from the equations of motion the acceleration along this coordinate is zero. The formula applies to simulations in which the selected coordinates are either unconstrained or constrained to fixed values. It is shown that in the latter case the formula reduces to the expression previously derived by den Otter and Briels. If simulations are carried out without constraining the coordinates of interest, the formula leads to a new method for calculating the free energy changes along these coordinates. This method is tested in two examples - rotation around the C-C bond of 1,2-dichloroethane immersed in water and transfer of fluoromethane across the water-hexane interface. The calculated free energies are compared with those obtained by two commonly used methods. One of them relies on determining the probability density function of finding the system at different values of the selected coordinate and the other requires calculating the average force at discrete locations along this coordinate in a series of constrained simulations. The free energies calculated by these three methods are in excellent agreement. The relative advantages of each method are discussed.

Darve, Eric↗