SEARCH · Engineering Papers

Results for “factorization 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Two-color method for optical astrometry - Theory and preliminary measurements with the Mark III stellar interferometer

The two-color method for interferometric astrometry provides a means of reducing the error in a stellar position measurement attributable to atmospheric turbulence. The primary limitation of the method is shown to be turbulent water vapor fluctuations. Secondary atmospheric effects caused by diffraction from small refractive-index inhomogeneities and differential refraction for the observation of stars away from zenith are shown to introduce errors that behave as white noise and which should not be significant. Other potential error sources due to photon noise, systematic instrumental effects, and imperfect data reduction are also considered. The improvement in accuracy possible with the two-color method is estimated as a factor of 5-10 over the corresponding one-color measurement. Some preliminary two-color measurements with the Mark III stellar interferometer at Mt. Wilson are presented, which demonstrate a factor of about 5 reduction in the amplitude of the atmospheric fluctuations in a stellar position measurement.

Colavita, M. Mark↗

NASA NTRS1987

Fracture of composites using singular finite elements

Six noded quarter point natural isoparametric triangular elements are employed to obtain displacement and stress distributions around the crack tip in unidirectional graphite epoxy. The performance of the quarter point element for predicting crack growth direction in unidirectional composite materials is studied for different loading conditions and fiber orientations. The finite element results are compared to analytical solutions of anisotropic elasticity where available. Notched Iosipescu shear specimens are analyzed to study fracture behavior under shear loading. Two methods for determining stress intensity factors are presented.

Choksi, G. N.↗

NASA NTRS1988

Developing fluid flow in a curved duct of square cross-section and its fully developed dual solutions

Developing fluid flow in a curved duct of square cross-section is studied numerically by a factored ADI finite-difference method on a staggered grid. A central-difference scheme with primitive variables is used inside the computational domain to reduce numerical diffusion. Two Reynolds numbers, 574 and 790, based upon a bulk velocity and hydraulic diameter are chosen for curvature ratios of 1/6.45 and 1/2.3, respectively. It is found that the secondary flow is far more complicated than expected, with the appearance of at least two pairs of vortices. Main-flow separation is also observed for the higher curvature ratio. Furthermore, it is observed that the flow develops into two quite different states downstream, depending upon the inlet conditions. Solutions of the fully developed Navier-Stokes equations is shown to be not unique beyond a certain critical Reynolds number. Developing flow seems to evolve into the fully developed state along a particular branch into which the fully developed solution bifurcates.

Soh, W. Y.↗

NASA NTRS1989

Algorithm For Solution Of Navier-Stokes Equations

Advantages of two previous algorithms combined. Steady-state Navier-Stokes equations put in implicit finite-difference form solved by approximate Newton iteration. LU-SSOR scheme is new relaxation method combining advantages of LU factorization with Gauss-Seidel relaxation. Vectorizable LU-SSOR scheme, based on central differences, requires scalar diagonal inversions. Application of scheme to approximate Newton iteration of Navier-Stokes equations yields set of equations requiring no implicit smoothing on left side. Only adaptive, total-variation-diminishing, flux-limited dissipation terms added to right side. Algorithm used to predict laminar, turbulent, and hypersonic flows.

Yoon, Seokkwan↗

NASA NTRS1990

Extension of a streamwise upwind algorithm to a moving grid system

A new streamwise upwind algorithm was derived to compute unsteady flow fields with the use of a moving-grid system. The temporally nonconservative LU-ADI (lower-upper-factored, alternating-direction-implicit) method was applied for time marching computations. A comparison of the temporally nonconservative method with a time-conservative implicit upwind method indicates that the solutions are insensitive to the conservative properties of the implicit solvers when practical time steps are used. Using this new method, computations were made for an oscillating wing at a transonic Mach number. The computed results confirm that the present upwind scheme captures the shock motion better than the central-difference scheme based on the beam-warming algorithm. The new upwind option of the code allows larger time-steps and thus is more efficient, even though it requires slightly more computational time per time step than the central-difference option.

Obayashi, Shigeru↗

NASA NTRS1989

Advances in upwind relaxation methods

Numerical techniques for solving the compressible Euler and Navier-Stokes equations are discussed with an emphasis on characteristic-based schemes. Two popular approaches, flux difference splitting and flux vector splitting, are described in one-dimensional Cartesian coordinates and then extended to three-dimensional generalized coordinates. A technique for increasing the spatial accuracy is presented, followed by a discussion of numerical dissipation mechanisms. An introduction to the use of implicit time integration schemes for accelerating the convergence rate to steady-state solutions including Newton's method, relaxation strategies, and approximate factorization techniques and their implementation on a vector processor concludes the chapter.

Walters, R. W.↗

NASA NTRS1990

Extension of transonic flow computational concepts in the analysis of cavitated bearings

An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings.

Vijayaraghavan, D.↗

NASA NTRS1990

Extension of transonic flow computational concepts in the analysis of cavitated bearings

Vijayaraghavan, D.↗

NASA NTRS1991

Laboratory simulation of cratering on small bodies

A new technique using external pressure was developed to simulate the lithostatic pressure due to self-gravity of small bodies. A 13-in. diameter cylindrical test chamber with L/D of 1 was fabricated to accommodate firing explosive charges with gas overpressures of up to 6000 psi. The chamber was hydrotested to 9000 psi. The method allows much larger scale factors that can be obtained with existing centrifuges and has the correct spherical geometry of self gravity. A simulant for jointed rock to be used in this fixture was developed using weakly cemented basalt. Various strength/pressure scaling theories can now be examined and tested.

Schmidt, Robert M.↗

NASA NTRS1992

A Comparison Between Raman Lidar and Conventional Contact Measurements of Atmospheric Temperature

Described here are the results of comparison between lidar and conventional contact measurements of the vertical temperature profile of the atmosphere. The lidar measurements are based on the method of temperature dependence of pure rotational Raman scattering of nitrogen and oxygen molecules. The presented results show that, as a whole, the motion of the lidar and conventional profiles coincide in the confidence intervals. Still there are districts in which the differences lay out of the confidence intervals. In analysis and estimation of the coincidence and the distinction in moving the profiles, we must take into account the differences between the methods for measuring the temperature profiles by lidar and contact methods, i.e. the next additional factors. The tied balloon meter gives the momentary values of the temperature in particular points of the profile, as these momentary values are the results of consequent, but not simultaneous measurements. In the case of free flying balloon, the time of measuring the lower several hundred meters is little because of high vertical velocity (about 300 m/min). This leads to an indefinite increase of the measuring error. The complete space coincidence between the lidar's and conventional profiles isn't possible. The contribution of the factors mentioned above about the error of comparison could increase because of nonstable layers in the planetary boundary layer of the atmosphere as well as by the influence of the mountains and the city situated nearby.

Mitev, V. M.↗

NASA NTRS1992

A boundary element alternating method for two-dimensional mixed-mode fracture problems

A boundary element alternating method, denoted herein as BEAM, is presented for two dimensional fracture problems. This is an iterative method which alternates between two solutions. An analytical solution for arbitrary polynomial normal and tangential pressure distributions applied to the crack faces of an embedded crack in an infinite plate is used as the fundamental solution in the alternating method. A boundary element method for an uncracked finite plate is the second solution. For problems of edge cracks a technique of utilizing finite elements with BEAM is presented to overcome the inherent singularity in boundary element stress calculation near the boundaries. Several computational aspects that make the algorithm efficient are presented. Finally, the BEAM is applied to a variety of two dimensional crack problems with different configurations and loadings to assess the validity of the method. The method gives accurate stress intensity factors with minimal computing effort.

Raju, I. S.↗

NASA NTRS1992

Recent Advances in Statistical Energy Analysis

Statistical Energy Analysis (SEA) has traditionally been developed using modal summation and averaging approach, and has led to the need for many restrictive SEA assumptions. The assumption of 'weak coupling' is particularly unacceptable when attempts are made to apply SEA to structural coupling. It is now believed that this assumption is more a function of the modal formulation rather than a necessary formulation of SEA. The present analysis ignores this restriction and describes a wave approach to the calculation of plate-plate coupling loss factors. Predictions based on this method are compared with results obtained from experiments using point excitation on one side of an irregular six-sided box structure. Conclusions show that the use and calculation of infinite transmission coefficients is the way forward for the development of a purely predictive SEA code.

K H Heron↗

NASA NTRS1992

Standards for space automation and robotics

The AIAA's Committee on Standards for Space Automation and Robotics (COS/SAR) is charged with the identification of key functions and critical technologies applicable to multiple missions that reflect fundamental consideration of environmental factors. COS/SAR's standards/practices/guidelines implementation methods will be based on reliability, performance, and operations, as well as economic viability and life-cycle costs, simplicity, and modularity.

Kader, Jac B.↗

NASA NTRS1992

A study of rain estimation methods from space using dual-wavelength radar measurements at near-nadir incidence over ocean

A question arising from the recent interest in spaceborne weather radar is what methods can be used to estimate precipitation parameters from space. In this paper, dual-wavelength airborne radar data obtained from flights conducted during 1988 and 1989 are used to compare rain rates derived from backscattering and attenuation methods. To help interpret the results the surface reference methods are studied by means of scatterplots of the surface cross sections at the two frequencies under rain and no-rain conditions. Approximate criteria are given on combining attenuation and backscattering methods to increase the effective dynamic range of the radar. The dual-wavelength capability of the radar is also used to examine the vertical structure of the precipitation. Another factor affecting the accuracy of the methods is the drop-size distribution. In the final section of the paper a procedure to estimate the profiled drop-size distribution is applied to the measured radar data.

Meneghini, R.↗

NASA NTRS1993

Aerospace devices for magnetic replicas

Retained persistent magnetic field has been studied and improved in the superconductor YBa2Cu3O7 (Y123). During the study, trapped magnetic field, B(t), has been increased by over a factor of 10(exp 5). Methods used to improve magnetic field trapping were principally: (1) the adoption of the Melt Texturing process to increase grain size; (2) the addition of excess Y to disperse deposits of Y2BaCuO5 (Y211) and again increase grain size; (3) irradiation with high energy particles including 1H+, 3He++, 4He++, and fission fragments; and (4) utilizing temperatures below 77 K has also been quantified as a way to increase trapped field. In addition, in our study of B(t), we have found laws governing creep, activation, temperature dependence, creep vs. current flow, etc. In the range 20 K less than or equal to T less than or equal to 65 K, and for B less than 10 Tesla, a simple empirical relationship was found: B(trap) (T2) = B(trap) (T1) ((Tc - T2)/(Tc - T1))squared where Tc is the critical temperature. The highest experimental trapped field was B(trap) = 3.96 Tesla, at 65 K. We believe this to be the highest persistent field ever produced, by any method. A two component model of the persistent currents has been developed. This accurately reproduces the data, using as parameters only the magnitude of a constant surface current, J(s), and a constant volume current J(v). The model successfully predicts B(t) (xyz) for the case of maximum trapped field, for all samples observed. It has also been extended to describe the unsaturated case either zero field cooled, or field cooled. Loss of strap with time has been studied for the critical state (Bt,max), and non critical state (Bt less than Bt,max), for times from a few minutes to a few months, for unirradiated material, for irradiation by 1H+, 3He++, 4He++, high z projectiles, and neutrons, and for all materials used in the overall study. We conclude that: (1) multi Tesla trapped fields are attained; (2) fields over 10 T are achievable; (3) creep is not a large problem; (4) application is feasible to motors, generators, magnets for particle beam optics, separators, levitating bearings, energy storage, shielding, and transportation.

Weinstein, Roy↗

NASA NTRS1952

An engineering method for estimating notch-size effect in fatigue tests on steel

Neuber's proposed method of calculating a practical factor of stress concentration for parts containing notches of arbitrary size depends on the knowledge of a "new material constant" which can be established only indirectly. In this paper, the new constant has been evaluated for a large variety of steels from fatigue tests reported in the literature, attention being confined to stresses near the endurance limit. Reasonably satisfactory results were obtained with the assumption that the constant depends only on the tensile strength of the steel. Even in cases where the notches were cracks of which only the depth was known, reasonably satisfactory agreement was found between calculated and experimental factors. It is also shown that the material constant can be used in an empirical formula to estimate the size effect on unnotched specimens tested in bending fatigue.

Kuhn, Paul↗

NASA NTRS1994

Fluorescence imaging to quantify crop residue cover

Crop residues, the portion of the crop left in the field after harvest, can be an important management factor in controlling soil erosion. Methods to quantify residue cover are needed that are rapid, accurate, and objective. Scenes with known amounts of crop residue were illuminated with long wave ultraviolet (UV) radiation and fluorescence images were recorded with an intensified video camera fitted with a 453 to 488 nm band pass filter. A light colored soil and a dark colored soil were used as background for the weathered soybean stems. Residue cover was determined by counting the proportion of the pixels in the image with fluorescence values greater than a threshold. Soil pixels had the lowest gray levels in the images. The values of the soybean residue pixels spanned nearly the full range of the 8-bit video data. Classification accuracies typically were within 3(absolute units) of measured cover values. Video imaging can provide an intuitive understanding of the fraction of the soil covered by residue.

Daughtry, C. S. T.↗

NASA NTRS1995

Nozzle flow with vibrational nonequilibrium

This research concerns the modeling and numerical solutions of the coupled system of compressible Navier-Stokes equations in cylindrical coordinates under conditions of equilibrium and nonequilibrium thermodynamics. The problem considered was the modeling of a high temperature diatomic gas N2 flowing through a converging-diverging high expansion nozzle. The problem was modeled in two ways. The first model uses a single temperature with variable specific heats as functions of this temperature. For the second model we assume that the various degrees of freedom all have a Boltzmann distribution and that there is a continuous redistribution of energy among the various degrees of freedom as the gas passes through the nozzle. Each degree of freedom is assumed to have its own temperature and, consequently, each system state can be characterized by these temperatures. This suggests that formulation of a second model with a vibrational degree of freedom along with a rotational-translation degree of freedom, each degree of freedom having its own temperature. Initially the vibrational degree of freedom is excited by heating the gas to a high temperature. As the high temperature gas passes through the nozzle throat there is a sudden drop in temperature along with a relaxation time for the vibrational degree of freedom to achieve equilibrium with the rotational-translation degree of freedom. That is, we assume that the temperature change upon passing through the throat is so great that the changes in the vibrational degree of freedom occur at a much slower pace and consequently lags behind the rotational-translational energy changes. This lag results in a finite relaxation time. In this context the term nonequilibrium is used to denote the fact that the energy content of the various degrees of freedom are characterized by two temperatures. We neglect any chemical reactions which could also add nonequilibrium effects. We develop the energy equations for the nonequilibrium model from first principles. The resulting equations, which model the nozzle flow, can be expressed in various forms. In most forms the resulting equations are coupled systems of nonlinear partial differential equations subject to certain boundary conditions. To solve the resulting coupled system of nonlinear partial differential equations, several numerical techniques were investigated: (1) the explicit MacCormack method, (2) the explicit-implicit MacCormack method, (3) the method of operator splitting, (4) factorization schemes, and (5) the Steger-Warming scheme.

Heinbockel, J. H.↗