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,063 records · Page 59

Phase Diagram of the Two-Chain Hubbard Model

We have calculated the charge gap and spin gap for the two-chain Hubbard model as a function of the on-site Coulomb interaction and the interchain hopping amplitude. We used the density matrix renormalization group method and developed a method to calculate separately the gaps numerically for the symmetric and antisymmetric modes with respect to the exchange of the chain indices. We have found very different behaviors for the weak and strong interaction cases. Our calculated phase diagram is compared to the one obtained by Balents and Fisher using the weak coupling renormalization group technique.

Park, Youngho↗

An equivalent domain integral for analysis of two-dimensional mixed mode problems

An equivalent domain integral (EDI) method for calculating J-integrals for two-dimensional cracked elastic bodies subjected to mixed mode loading is presented. The total and product integrals consist of the sum of an area or domain integral and line integrals on the crack faces. 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 the problems analyzed.

Raju, I. S.↗

Lift distribution and longitudinal stability of an airplane

The preliminary calculation of the airplane polar and hence of the flight performances and characteristics rests on the assumption of an elliptical lift distribution at all altitudes. For large angles of attack below C(sub a (sub max)), this method of calculation yields no satisfactory agreement with measurements made in flight. An attempt is made to eliminate the errors in the preliminary calculation by the assumption of a disturbance of the lift distribution in this angle-of-attack range, which is so important for the constructor. An explanation is also given of the great differences found in flight with and without propeller slipstream.

Topfer, Carl↗

Numerical computation of transonic flows with shock waves

Some recent developments in numerical methods for calculating solutions to the transonic potential flow equation are reviewed, including (1) the construction of the stable coordinate independent difference schemes; (2) the use of conservation form to insure proper shock jump conditions; (3) analysis of the relaxation method by the time dependent analogy; (4) accelerated iterative schemes.

Jameson, A.↗

Static and dynamic structural-sensitivity derivative calculations in the finite-element-based Engineering Analysis Language (EAL) system

The implementation of static and dynamic structural-sensitivity derivative calculations in a general purpose, finite-element computer program denoted the Engineering Analysis Language (EAL) System is described. Derivatives are calculated with respect to structural parameters, specifically, member sectional properties including thicknesses, cross-sectional areas, and moments of inertia. Derivatives are obtained for displacements, stresses, vibration frequencies and mode shapes, and buckling loads and mode shapes. Three methods for calculating derivatives are implemented (analytical, semianalytical, and finite differences), and comparisons of computer time and accuracy are made. Results are presented for four examples: a swept wing, a box beam, a stiffened cylinder with a cutout, and a space radiometer-antenna truss.

Camarda, C. J.↗

A turbulence model for nonequilibrium adverse pressure gradient flows

This paper presents a method for calculating compressible adverse pressure gradient boundary layers by using an algebraic eddy viscosity turbulence model that has been modified for variable pressure gradient and turbulence memory effects. The pressure gradient corrections are based on previous incompressible data correlations. Several methods for including the effects of turbulence memory are evaluated. A new lag model, which gives good agreement with available experimental data, is developed. Finally, a correlation is developed for the lag length parameter employed in the model as a function of the known experimental flow variables.

Horstman, C. C.↗

Estimation of additive forces and moments for supersonic inlets

A technique for estimating the additive forces and moments associated with supersonic, external compression inlets as a function of mass flow ratio has been developed. The technique makes use of a low order supersonic paneling method for calculating minimum additive forces at maximum mass flow conditions. A linear relationship between the minimum additive forces and the maximum values for fully blocked flow is employed to obtain the additive forces at a specified mass flow ratio. The method is applicable to two-dimensional inlets at zero or nonzero angle of attack, and to axisymmetric inlets at zero angle of attack. Comparisons with limited available additive drag data indicate fair to good agreement.

Perkins, Stanley C., Jr.↗

High frequency scattering from trihedral corner reflectors and other benchmark targets - SBR versus experiment

A general method for calculating the radar cross section (RCS) from a three-dimensional target is described. The target is first constructed by using a solid-geometry-modeling computer-aided design (CAD) package. Following the shooting and bouncing ray (SBR) method, a very dense grid of rays is launched from the incident direction toward the target. Each ray is traced according to the geometrical optics theory including the effect of ray tube divergence, polarization, and material reflection coefficient. At the point where the ray exits the target, a physical optics-type integration is performed to obtain the scattered far fields. This method is tested using several simple examples involving interaction among plates, cylinders, and spheres. The theoretical results are generally in good agreement with measured data.

Baldauf, John↗

Theoretical evaluation of high-speed aerodynamics for arrow-wing configurations

The use of the theoretical methods to calculate the high-speed aerodynamic characteristics of arrow-wing supersonic cruise configurations was studied. Included are correlations of theoretical predictions with wind-tunnel data at Mach numbers from 0.8 to 2.7, examples of the use of theoretical methods to extrapolate the wind-tunnel data to full-scale flight condition, and presentation of a typical supersonic data package for an advanced supersonic transport application. A brief description of the methods and their application is given.

Dollyhigh, S. M.↗

Velocities of guided ultrasonic waves in heterogeneous medium

Experimental and theoretical examinations were performed of the longitudinal velocity characteristics of waves in trilaminar and encapsulated waveguides. The study was confined to waveguides with core material that featured transverse wave velocities much worse than the longitudinal wave velocities. The velocities were obtained using a dispersion equation, with consideration given to both the core and encapsulant. Asymptotic velocities were also calculated for bending and twisting in trilaminar waveguides. Trials were run with bimetallic waveguides for comparison with the theoretical predictions. Good agreement was found between the predicted velocity of the propagation of the fundamental mode and the measured velocities. The method was calculated valid for modes above four, confirming that the data were contained in either the core or outer layer, and were unsensitive to the encapsulant.

Touratier, M.↗

Crippling Strength of Axially Loaded Rods

A new empirical formula was developed that holds good for any length and any material of a rod, and agrees well with the results of extensive strength tests. To facilitate calculations, three tables are included, giving the crippling load for solid and hollow sectioned wooden rods of different thickness and length, as well as for steel tubes manufactured according to the standards of Army Air Services Inspection. Further, a graphical method of calculation of the breaking load is derived in which a single curve is employed for determination of the allowable fiber stress. Finally, the theory is discussed of the elastic curve for a rod subject to compression, according to which no deflection occurs, and the apparent contradiction of this conclusion by test results is attributed to the fact that the rods under test are not perfectly straight, or that the wall thickness and the material are not uniform. Under the assumption of an eccentric rod having a slight initial bend according to a sine curve, a simple formula for the deflection is derived, which shows a surprising agreement with test results. From this a further formula is derived for the determination of the allowable load on an eccentric rod. The resulting relations are made clearer by means of a graphical representation of the relation of the moments of the outer and inner forces to the deflection.

Natalis, FR↗

Ice Accretion Calculations for a Commercial Transport Using the LEWICE3D, ICEGRID3D and CMARC Programs

Collection efficiency and ice accretion calculations were made for a commercial transport using the NASA Lewis LEWICE3D ice accretion code, the ICEGRID3D grid code and the CMARC panel code. All of the calculations were made on a Windows 95 based personal computer. The ice accretion calculations were made for the nose, wing, horizontal tail and vertical tail surfaces. Ice shapes typifying those of a 30 minute hold were generated. Collection efficiencies were also generated for the entire aircraft using the newly developed unstructured collection efficiency method. The calculations highlight the flexibility and cost effectiveness of the LEWICE3D, ICEGRID3D, CMARC combination.

Bidwell, Colin S.↗

Research in Computational Astrobiology

We report on several projects in the field of computational astrobiology, which is devoted to advancing our understanding of the origin, evolution and distribution of life in the Universe using theoretical and computational tools. Research projects included modifying existing computer simulation codes to use efficient, multiple time step algorithms, statistical methods for analysis of astrophysical data via optimal partitioning methods, electronic structure calculations on water-nuclei acid complexes, incorporation of structural information into genomic sequence analysis methods and calculations of shock-induced formation of polycylic aromatic hydrocarbon compounds.

Chaban, Galina↗

A numerical method for determining the temperature structure of planetary atmospheres.

A numerical method for calculating the time-average, vertical temperature structure of planetary atmospheres is presented. It is assumed that the atmospheres are in radiative-convective equilibrium, which is a good first approximation to many situations. Numerical tests of the rate of convergence and accuracy of the answer are presented. The method can readily handle molecular sources of opacity. Accurate results can be obtained with a minimum of computer time, because the number of iterations needed (about 4) is small and the number of pressure levels at which the net flux needs to be evaluated (about 10) is small. As an application of this procedure, some model atmospheres of Jupiter are calculated.

Pollack, J. B.↗

Vector sweep

Program is numerical procedure for calculating radiation/geometrical configuration factors and/or space-environmental heat-flux histories. System can accomodate complex systems of diverse geometric shapes residing in unlimited variety of mutual shadow conditions. Basic analytical device is finite-difference employment of unit-sphere method of calculating configuration factors.

Moutrie, C. L.↗

The desorptivity model of bulk soil-water evaporation

Available models of bulk evaporation from a bare-surfaced soil are difficult to apply to field conditions where evaporation is complicated by two main factors: rate-limiting climatic conditions and redistribution of soil moisture following infiltration. Both factors are included in the "desorptivity model', wherein the evaporation rate during the second stage (the soil-limiting stage) of evaporation is related to the desorptivity parameter, A. Analytical approximations for A are presented. The approximations are independent of the surface soil moisture. However, calculations using the approximations indicate that both soil texture and soil moisture content at depth significantly affect A. Because the moisture content at depth decreases in time during redistribution, it follows that the A parameter also changes with time. Consequently, a method to calculate a representative value of A was developed. When applied to field data, the desorptivity model estimated cumulative evaporation well. The model is easy to calculate, but its usefulness is limited because it requires an independent estimate of the time of transition between the first and second stages of evaporation. The model shows that bulk evaporation after the transition to the second stage is largely independent of climatic conditions.

Clapp, R. B.↗