Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical integration”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

A 4-node Bilinear Isoparametric Element in Rockwell NASTRAN

The development and evaluation of the Rockwell NASTRAN four node quadrilateral (QUAD4) element is presented. The element derivation utilizes bilinear isoparametric techniques both for membrane and bending characteristics. The QUAD4 element coordinate system, membrane properties, lumped mass matrix, and treatment of warping are based upon the COSMIC/NASTRAN QDMEMI element while the bending characteristics are based on a paper by T. J. R. Hughes. The effects of warping on the bending stiffness, consistent mass, and geometric stiffness are based upon a paper by R. H. MacNeal. Numerical integration is accomplished by Gaussian quadrature on a 2 x 2 grid. Practical user support features include variable element thickness, thermal analysis and layered composite material definitions.

Liao, C.↗
NASA NTRS1985

Vortex studies relating to boundary layer turbulence and noise

The present study considers the two-dimensional case of an array of N rectilinear, like-sign vortices above an infinite flat boundary. The method of images can be employed with this configuration to reduce the problem to that of 2N vortices in free space, constrained by 2N symmetry relations. This system is Hamiltonian and therefore certain invariants of the motion are known. Futher, from the Hamiltonian constant, the equations of motion are readily derived and may be integrated numerically to determine the vortex trajectories. This knowledge of the time-dependent vortex motion then allows the resulting noise radiation to be computed by standard aeroacoustic techniques.

Adelman, J. L.↗
NASA NTRS1985

Preliminary Results from Numerical Experiments on the Summer 1980 Heat Wave and Drought

During the summer of 1980, a prolonged heat wave and drought affected the United States. A preliminary set of experiments has been conducted to study the effect of varying boundary conditions on the GLA model simulation of the heat wave. Five 10-day numerical integrations with three different specifications of boundary conditions were carried out: a control experiment which utilized climatological boundary conditions, an SST experiment which utilized summer 1980 sea-surface temperatures in the North Pacific, but climatological values elsewhere, and a Soil Moisture experiment which utilized the values of Mintz-Serafini for the summer, 1980. The starting dates for the five forecasts were 11 June, 7 July, 21 July, 22 August, and 6 September of 1980. These dates were specifically chosen as days when a heat wave was already established in order to investigate the effect of soil moistures or North Pacific sea-surface temperatures on the model's ability to maintain the heat wave pattern. The experiments were evaluated in terms of the heat wave index for the South Plains, North Plains, Great Plains and the entire U.S. In addition a subjective comparison of map patterns has been performed.

Wolfson, N.↗
NASA NTRS1985

On Fofonoff's Mode

The classical barotropic ocean model of Veronis (1966) is reexamined to verify the possible existence of Fofonoff's (1954) steady free inertial mode by considering the balance where J is the Jacobian, eta is the streamfunction of the horizontal velocity, and R is a parameter of nonlinearity. A linear functional relationship between the potential vorticity and the streamfunction and inertial boundary layers along the ocean basin was suggested by Fofonoff. However, the numerical integration showed that the system is barotropically unstable. The typical instantaneous field is given. The statistical equilibrium state eta is also given. Although it resembles Fofonoff's solution, it is maintained by the Reynolds stress field induced by the eddies in addition to the wind stress and dissipation. The plot of eta versus R delta sq +y is given. It is shown that the Fofonoff's linear relation between the potential vorticity and the streamfunction does not hold uniformly over the entire basin. The above results show that the Fofonoff's solution can not be realized physically.

Merkine, L. O.↗
NASA NTRS1985

SPAR improved structural-fluid dynamic analysis capability

The results of a study whose objective was to improve the operation of the SPAR computer code by improving efficiency, user features, and documentation is presented. Additional capability was added to the SPAR arithmetic utility system, including trigonometric functions, numerical integration, interpolation, and matrix combinations. Improvements were made in the EIG processor. A processor was created to compute and store principal stresses in table-format data sets. An additional capability was developed and incorporated into the plot processor which permits plotting directly from table-format data sets. Documentation of all these features is provided in the form of updates to the SPAR users manual.

Pearson, M. L.↗
NASA NTRS1985

Application of closed-form solutions to a mesh point field in silicon solar cells

A computer simulation method is discussed that provides for equivalent simulation accuracy, but that exhibits significantly lower CPU running time per bias point compared to other techniques. This new method is applied to a mesh point field as is customary in numerical integration (NI) techniques. The assumption of a linear approximation for the dependent variable, which is typically used in the finite difference and finite element NI methods, is not required. Instead, the set of device transport equations is applied to, and the closed-form solutions obtained for, each mesh point. The mesh point field is generated so that the coefficients in the set of transport equations exhibit small changes between adjacent mesh points. Application of this method to high-efficiency silicon solar cells is described; and the method by which Auger recombination, ambipolar considerations, built-in and induced electric fields, bandgap narrowing, carrier confinement, and carrier diffusivities are treated. Bandgap narrowing has been investigated using Fermi-Dirac statistics, and these results show that bandgap narrowing is more pronounced and that it is temperature-dependent in contrast to the results based on Boltzmann statistics.

Lamorte, M. F.↗
NASA NTRS1985

Electrophoretic mobility of cells in a vertical Ficoll gradiant

The upward migration of living cells and test particles under the influence of a constant electric field in a low conductivity Ficoll gradient occurs at nearly constant velocity. Viscosity and neutral polymer concentration affect migration rate. Decreasing viscosity speeds up the particle migration, decreasing neutral polymer (Ficoll) concentration, slows particle migration, since electrophoretic mobility increases approximately linearly with neutral polymer concentration. Neutral polymers interact with the cell surface to effectively raise its zeta potential. An analytic function was developed from the known dependence of these physical variables on migration distance; the analysis expresses migration velocity as an explicit function of position in the density gradient. It predicts an almost linear increase in velocity of about 12 to 16% over the working region of the gradient. It was numerically integrated and correctly predicts cell migration distance vs time curves without the use of any fitted parameters. The resulting migration curves follow the expected slowly varying exponential form that closely resembles a straight line. The ability to determine standard electrophoretic mobilities from such curves depends on knowledge of the effect of Ficoll on the zeta potential of the cell type that is separated.

Plank, L. D.↗
NASA NTRS1985

Pulsars as cosmic ray particle accelerators: Proton orbits

Proton orbits are calculated in the electromagnetic vacuum field of a magnetic point dipole rotating with its angular velocity omega perpendicular to its dipole moment mu by numerical integration of the Lorentz-Dirac equation. Trajectories are shown and discussed for various initial conditions. A critical surface is shown separating initial positions of protons which finally hit the pulsar in the polar region from those which finally recede to infinity.

Thielheim, K. O.↗
NASA NTRS1985

Pulsars as cosmic ray accelerators: Energy development of protons

Results of numerical integrations of the Lorentz-Dirac-equation with Landau approximation for protons in the electromagnetic vacuum field of a magnetic dipole rotating with its angular velocity omega perpendicular to its magnetic moment mu will be discussed with special attention to the energy development of protons.

Thielheim, K. O.↗
NASA NTRS1985

A new world survey expression for cosmic ray vertical intensity vs. depth in standard rock

The cosmic ray data on vertical intensity versus depth below 10 to the 5th power g sq cm is fitted to a 5 parameter empirical formula to give an analytical expression for interpretation of muon fluxes in underground measurements. This expression updates earlier published results and complements the more precise curves obtained by numerical integration or Monte Carlo techniques in which the fit is made to an energy spectrum at the top of the atmosphere. The expression is valid in the transitional region where neutrino induced muons begin to be important, as well as at great depths where this component becomes dominant.

Crouch, M.↗
NASA NTRS1984

Final stages of transition to turbulence in plane channel flow

This paper involves a numerical simulation of the final stages of transition to turbulence in plane channel flow at a Reynolds number of 1500. Three-dimensional incompressible Navier-Stokes equations are numerically integrated to obtain the time evolution of two- and three-dimensional finite-amplitude disturbances. Computations are performed on the CYBER-203 vector processor for a 32 x 51 x 32 grid. Solutions indicate the existence of structures similar to those observed in the laboratory and characteristics of the various stages of transition that lead to final breakdown. In particular, evidence points to the formation of a upside-down-V-shaped vortex and the subsequent system of horseshoe vortices inclined to the main flow direction as the primary elements of transition. Details of the resulting flow field after breakdown indicate the evolution of streaklike formations found in turbulent flows. Although the flow field does approach a steady state (turbulent channel flow), the introduction of subgrid-scale terms seems necessary to obtain fully developed turbulence statistics.

Biringen, S.↗
NASA NTRS1985

Bounds and approximations for rapid evaluation of coherent MPSK error probabilities

The error probability for MPSK detection is expressed here as a sum of a familiar upper-bound approximation and a correction term which is well suited for numerical integration. An asymptotic expansion of the correction term is derived, from which simple and extremely tight upper and lower bounds are obtained. Relative errors in using these bounds and approximations of them are graphed for cases of practical interest.

Chie, C. M.↗
NASA NTRS1985

A comparison of observed and forecast energetics over North America

The observed kinetic energy balance is calculated over North America and compared with that computed from forecast fields for the 13-15 January 1979 cyclone. The FGGE upper-air rawinsonde network serves as the observational database while the forecast energetics are derived from a numerical integration with the GLAS fourth-order general circulation model initialized at 00 GMT 13 January. Maps of the observed and predicted kinetic energy and eddy conversion are in good qualitative agreement, although the model eddy conversion tends to be 2 to 3 times stronger than the observed values. Both the forecast and observations exhibit the lower and upper tropospheric maxima in vertical profiles of kinetic energy generation and dissipation typically found in cyclonic disturbances. An interesting time lag is noted in the observational analysis with the maximum observed kinetic energy occurring 12 h later than the maximum eddy conversion over the same region.

Baker, W. E.↗
NASA NTRS1985

Transition probabilities of the B-prime 1Sigma(u)(+) to X 1Sigma(g)(+) system of molecular hydrogen

From published potential energy curves and transition dipole moments, there are obtained by numerical integration of the equations of nuclear motion the vibrational eigenfunctions of the X 1Sigma(g)(+) and B-prime 1Sigma(u)(+) states of H2. The probabilities of radiative transitions from the discrete vibrational levels of the excited B-prime 1Sigma(u)(+) electronic state of H2 to the discrete and continuum vibrational levels of the ground X 1Sigma(g)(+) electronic state are calculated. The Franck-Condon factors are also presented.

Kwok, T. L.↗
NASA NTRS1985

Meteorites may follow a chaotic route to earth

A numerical integration was employed to show that the 3/1 Kirkwood gap region between Jupiter and the sun can serve as the source for asteroidal debris that strikes the earth as chondrites. A previous study by Wetherill (1968) considered planetary close encounters, secular resonances, collisional destruction in space and subsequent atmospheric abrasion of the chondrites. The factors would have resulted in chaotic trajectories for asteroids in the asteroid belt, i.e., eccentricities large enough to make the objects earth-crossing. A sampling of trajectories which could reach large eccentricities is integrated in terms of deterministic solutions of Hamiltonian systems of equations which behave erratically and are dependent on the initial conditions. The examinations successfully predicted a trajectory with a peak eccentricity of 0.652, a time interval of 340,000 yr, and a semimajor axis between 2.4-2.53 AU that was earth-crossing.

Wisdom, J.↗
NASA NTRS1985

Rotordynamic coefficients for compressible flow in tapered annular seals

Derivation of the governing equations for compressible flow in a tapered annular seal is based on Hirs' turbulent bulk-flow model. Zeroth and first-order perturbation equations are developed by an expansion in the eccentricity ratio. These equations are numerically integrated to obtain the leakage, and the direct and cross-coupled stiffness and damping coefficients. Seal parameters similar to the Space Shuttle Main Engine High Pressure Oxidizer Turbopump are used to demonstrate output from the analysis procedure. The effects of preswirl and seal taper are shown for three different length-to-diameter ratios. Generally the results indicate that prerotating the fluid significantly increases the cross-coupled stiffness but has little effect on the other coefficients, and increasing the convergent taper increases the direct stiffness while decreasing the direct damping and cross-coupled stiffness.

Nelson, C. C.↗
NASA NTRS1985

Tolerance bounds for log gamma regression models

The present procedure for finding lower confidence bounds for the quantiles of Weibull populations, on the basis of the solution of a quadratic equation, is more accurate than current Monte Carlo tables and extends to any location-scale family. It is shown that this method is accurate for all members of the log gamma(K) family, where K = 1/2 to infinity, and works well for censored data, while also extending to regression data. An even more accurate procedure involving an approximation to the Lawless (1982) conditional procedure, with numerical integrations whose tables are independent of the data, is also presented. These methods are applied to the case of failure strengths of ceramic specimens from each of three billets of Si3N4, which have undergone flexural strength testing.

Jones, R. A.↗
NASA NTRS1985

Ground based concept for time control of aircraft entering the terminal area

This paper summarizes current research for increasing the efficiency of traffic flow at hub airports by application of four dimensional (4D) guidance techniques. A method for generating 4D guidance commands to control the landing time of aircraft that are not equipped with on-board 4D guidance systems is described. In one possible implementation the commands are generated in an auxiliary processor linked to the ATC host computer and displayed on a controller's monitor in the form of profile descent advisories. Accurate time control is achieved by generating air traffic controller advisories in a ground-based algorithm that combines aerodynamic, thrust and atmospheric models with an efficient numerical integration method. The time accuracy and fuel efficiency achieved when the pilot responds to these advisories were evaluated in a piloted simulation of a transport aircraft.

Erzberger, H.↗