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 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1985

A fast, time-accurate unsteady full potential scheme

The unsteady form of the full potential equation is solved in conservation form by an implicit method based on approximate factorization. At each time level, internal Newton iterations are performed to achieve time accuracy and computational efficiency. A local time linearization procedure is introduced to provide a good initial guess for the Newton iteration. A novel flux-biasing technique is applied to generate proper forms of the artificial viscosity to treat hyperbolic regions with shocks and sonic lines present. The wake is properly modeled by accounting not only for jumps in phi, but also for jumps in higher derivatives of phi, obtained by imposing the density to be continuous across the wake. The far field is modeled using the Riemann invariants to simulate nonreflecting boundary conditions. The resulting unsteady method performs well which, even at low reduced frequency levels of 0.1 or less, requires fewer than 100 time steps per cycle at transonic Mach numbers. The code is fully vectorized for the CRAY-XMP and the VPS-32 computers.

Shankar, V.↗

NASA NTRS1974

The radiative effect of aerosols in the earth's atmosphere

A modified two-flux approximation is employed to compute the transfer of radiation in a finite, inhomogeneous, turbid atmosphere. A perturbation technique is developed to allow the treatment of nongray gaseous absorption with multiple scattering. The perturbation method, which employs a backscatter factor as a parameter, can be used with anisotropic particle scattering as well as Rayleigh scattering. This method is used to study the effect of aerosols on radiative solar heating and infrared cooling as well as the radiative-convective temperature distribution in the earth's atmosphere. It is found that the effect of aerosols in the infrared cannot be neglected; while in the visible, the effect can be of the same order as that due to absorption by water vapor. For a high surface albedo (greater than 0.30) heating of the earth-atmosphere system results due to the presence of aerosols. The aerosols also reduce the amount of convection needed to maintain a stable atmosphere. For the case of a dense haze a temperature inversion is found to exist close to the ground.

Wang, W.-C.↗

NASA NTRS1952

Radiant-interchange Configuration Factors

A study is presented of the geometric configuration factors required for computing radiant heat transfer between opaque surfaces separated by a nonabsorbing medium and various methods of determining the configuration factors are discussed. Configuration-factor solutions available in the literature have been checked and the more complicated equations are presented as families of curves. Cases for point, line, and finite-area sources are worked out over a wide range of geometric proportions. These cases include several new configurations involving rectangles, triangles, and cylinders of finite length which are integrated and tabulated. An analysis is presented, in which configuration factors are employed of the radiant heat transfer to the rotor blades of a typical gas turbine under different conditions of temperature and pressure. (author)

Hamilton, D C :↗

NASA NTRS1992

Organic crystal growth experiment facility (13-IML-1)

The interesting nature of metal-like organic compounds composed of charge transfer complexes has been recently realized. Crystals of these complexes can usually be grown by the solution crystallization method. It is difficult to grow such organic crystals on Earth, especially from the chemical reactions through diffusion controlled process in the solutions, because of gravitational disturbances, or sedimentation. The International Microgravity Lab. (IML-1) Organic Crystal Growth with G-Gitter Preventive Measure (OCGP) experiment is expected to grow a single crystal large enough to allow its intrinsic physical properties to be measured and its detailed crystal structure to be determined. This experiment also attempts to assess the experimental conditions including the microgravity environment for further study of the fundamental process of solution crystallization, nucleation, and growth from supersaturated phases including chemical reactions. Microgravity disturbances, G-jitter, may be an important environmental factor in the experimental method to assess. The vibration damping effects on organic crystal growth can be carefully studied.

Kanbayashi, Akio↗

NASA NTRS1975

Satellite animal tracking feasibility studies

A study was initiated in Tsavo National Park to determine movements and home ranges of individual elephants and their relations to overall distribution patterns and environmental factors such as rainfall. Methods used were radio tracking and observations of visually identifiable individuals. Aerial counts provided data on overall distribution. Two bulls and two cows were radio-tagged in Tsavo West and two bulls and four cows in Tsavo East, providing home range and movement data. The movements of individuals were useful in interpreting relatively major shifts in elephant distribution. Results point to the following preliminary conclusions: (1) elephants in the Tsavo area undertook long distance movements in fairly direct response to localized rainfall; (2) a subdivision of the overall population into locally distinct units may exist during the dry season but did not occur after significant rainfall; and (3) food appears to be the primary factor governing movements and distribution of elephants in the area.

Buechner, H. K.↗

NASA NTRS1987

Numerical study of laminar and turbulent flows inside a turnaround duct with and without guide vanes

The purpose of this study is to examine in detail incompressible laminar and turbulent flows inside a turnaround duct with and without guide vanes and to investigate the effects of vanes on the flow characteristics. To perform this study, an implicit finite difference code cast in general curvilinear coordinates is further developed. The code is based on the method of pseudo-compressibility and utilize ADI or implicit approximate factorization algorithm to achieve computational efficiency. Method of segmental sweeping is developed to overcome the multiple-zone problem due to imposition of guide vanes. In the present study, several test cases have been computed. These test cases include laminar and turbulent flows inside a turnaround duct without and with two or three guide vanes. The study reveals that: (1) there exists large recirculation zones inside the duct if no vanes are present; (2) properly shaped and positional guide vanes are effective in eliminating flow separation; and (3) laminar and turbulent flows have similar flow features. But turbulent flow has less total pressure drop.

Lin, S.-J.↗

NASA NTRS1987

A Navier-Stokes flow simulation of the Space Shuttle Main Engine Hot Gas Manifold

Incompressible viscous flow inside the turnaround duct, the fuel bowl, the transfer duct and the racetrack of the Space Shuttle Main Engine (SSME) Hot Gas Manifold (HGM) has been computed using the method of pseudo-compressibility together with an implicit, approximate-factorization algorithm. A multiple-zone method is used to make solution of flows in complex geometries easy. A model which predicts the pressure loading for the shield and the injector post arrangement without solving the complex flow field in the main injector region is proposed. The computed results show good qualitative agreement with experimental data.

Yang, Ruey-Jen↗

NASA NTRS1995

Design and Testing of an Automated System using Thermochromatic Liquid Crystals to Determine Local Heat Transfer Coefficients for an Impinging Jet

Using thermochromatic liquid crystal to measure surface temperature, an automated transient method with time-varying free-stream temperature is developed to determine local heat transfer coefficients. By allowing the free-stream temperature to vary with time, the need for complicated mechanical components to achieve a step temperature change is eliminated, and by using the thermochromatic liquid crystals as temperature indicators, the labor intensive task of installing many thermocouples is omitted. Bias associated with human perception of the transition of the thermochromatic liquid crystal is eliminated by using a high speed digital camera and a computer. The method is validated by comparisons with results obtained by the steady-state method for a circular Jet impinging on a flat plate. Several factors affecting the accuracy of the method are evaluated.

Tan, Benjamin↗

NASA NTRS1989

Efficient computer algebra algorithms for polynomial matrices in control design

The theory of polynomial matrices plays a key role in the design and analysis of multi-input multi-output control and communications systems using frequency domain methods. Examples include coprime factorizations of transfer functions, cannonical realizations from matrix fraction descriptions, and the transfer function design of feedback compensators. Typically, such problems abstract in a natural way to the need to solve systems of Diophantine equations or systems of linear equations over polynomials. These and other problems involving polynomial matrices can in turn be reduced to polynomial matrix triangularization procedures, a result which is not surprising given the importance of matrix triangularization techniques in numerical linear algebra. Matrices with entries from a field and Gaussian elimination play a fundamental role in understanding the triangularization process. In the case of polynomial matrices, matrices with entries from a ring for which Gaussian elimination is not defined and triangularization is accomplished by what is quite properly called Euclidean elimination. Unfortunately, the numerical stability and sensitivity issues which accompany floating point approaches to Euclidean elimination are not very well understood. New algorithms are presented which circumvent entirely such numerical issues through the use of exact, symbolic methods in computer algebra. The use of such error-free algorithms guarantees that the results are accurate to within the precision of the model data--the best that can be hoped for. Care must be taken in the design of such algorithms due to the phenomenon of intermediate expressions swell.

Baras, J. S.↗

NASA NTRS1988

Storage and computationally efficient permutations of factorized covariance and square-root information matrices

A unified method to permute vector-stored upper-triangular diagonal factorized covariance (UD) and vector stored upper-triangular square-root information filter (SRIF) arrays is presented. The method involves cyclical permutation of the rows and columns of the arrays and retriangularization with appropriate square-root-free fast Givens rotations or elementary slow Givens reflections. A minimal amount of computation is performed and only one scratch vector of size N is required, where N is the column dimension of the arrays. To make the method efficient for large SRIF arrays on a virtual memory machine, three additional scratch vectors each of size N are used to avoid expensive paging faults. The method discussed is compared with the methods and routines of Bierman's Estimation Subroutine Library (ESL).

Muellerschoen, R. J.↗

NASA NTRS1946

Stress Analysis of Columns and Beam Columns by the Photoelastic Method

Principles of similarity and other factors in the design of models for photoelastic testing are discussed. Some approximate theoretical equations, useful in the analysis of results obtained from photoelastic tests are derived. Examples of the use of photoelastic techniques and the analysis of results as applied to uniform and tapered beam columns, circular rings, and statically indeterminate frames, are given. It is concluded that this method is an effective tool for the analysis of structures in which column action is present, particularly in tapered beam columns, and in statically indeterminate structures in which the distribution of loads in the structures is influenced by bending moments due to axial loads in one or more members.

BEAM-COLUMNS - STRESSES↗

NASA NTRS1991

An iterative method for obtaining the optimum lightning location on a spherical surface

A brief introduction to the basic principles of an eigen method used to obtain the optimum source location of lightning is presented. The location of the optimum source is obtained by using multiple direction finders (DF's) on a spherical surface. An improvement of this method, which takes the distance of source-DF's as a constant, is presented. It is pointed out that using a weight factor of signal strength is not the most ideal method because of the inexact inverse signal strength-distance relation and the inaccurate signal amplitude. An iterative calculation method is presented using the distance from the source to the DF as a weight factor. This improved method has higher accuracy and needs only a little more calculation time. Some computer simulations for a 4DF system are presented to show the improvement of location through use of the iterative method.

Chao, Gao↗

NASA NTRS1984

Methods for analysis of cracks in three-dimensional solids

Analytical and numerical methods evaluating the stress-intensity factors for three-dimensional cracks in solids are presented, with reference to fatigue failure in aerospace structures. The exact solutions for embedded elliptical and circular cracks in infinite solids, and the approximate methods, including the finite-element, the boundary-integral equation, the line-spring models, and the mixed methods are discussed. Among the mixed methods, the superposition of analytical and finite element methods, the stress-difference, the discretization-error, the alternating, and the finite element-alternating methods are reviewed. Comparison of the stress-intensity factor solutions for some three-dimensional crack configurations showed good agreement. Thus, the choice of a particular method in evaluating the stress-intensity factor is limited only to the availability of resources and computer programs.

Raju, I. S.↗

NASA NTRS1983

Vibration Analysis Reduces Computer Time

Improved calculation method promises to reduce computer time by factor of 10 for vibration analysis of complex structures. Method originally developed to predict mechanical response of spacecraft to maneuvers in various cases in which same booster is used with different payloads. Method is more exact than conventional normal-mode methods; it involves no initial approximations of assumptions. Calculation cost is reduced without loss of accuracy.

Engels, R. C.↗

NASA NTRS1984

Three-dimensional finite-element analysis of chevron-notched fracture specimens

Stress-intensity factors and load-line displacements were calculated for chevron-notched bar and rod fracture specimens using a three-dimensional finite-element analysis. Both specimens were subjected to simulated wedge loading (either uniform applied displacement or uniform applied load). The chevron-notch sides and crack front were assumed to be straight. Crack-length-to-specimen width ratios (a/w) ranged from 0.4 to 0.7. The width-to-thickness ratio (w/B) was 1.45 or 2. The bar specimens had a height-to-width ratio of 0.435 or 0.5. Finite-element models were composed of singularity elements around the crack front and 8-noded isoparametric elements elsewhere. The models had about 11,000 degrees of freedom. Stress-intensity factors were calculated by using a nodal-force method for distribution along the crack front and by using a compliance method for average values. The stress intensity factors and load-line displacements are presented and compared with experimental solutions from the literature. The stress intensity factors and load-line displacements were about 2.5 and 5 percent lower than the reported experimental values, respectively.

Raju, I. S.↗

NASA NTRS1995

Sensitivity analysis for aeroacoustic and aeroelastic design of turbomachinery blades

A new method for computing the effect that small changes in the airfoil shape and cascade geometry have on the aeroacoustic and aeroelastic behavior of turbomachinery cascades is presented. The nonlinear unsteady flow is assumed to be composed of a nonlinear steady flow plus a small perturbation unsteady flow that is harmonic in time. First, the full potential equation is used to describe the behavior of the nonlinear mean (steady) flow through a two-dimensional cascade. The small disturbance unsteady flow through the cascade is described by the linearized Euler equations. Using rapid distortion theory, the unsteady velocity is split into a rotational part that contains the vorticity and an irrotational part described by a scalar potential. The unsteady vorticity transport is described analytically in terms of the drift and stream functions computed from the steady flow. Hence, the solution of the linearized Euler equations may be reduced to a single inhomogeneous equation for the unsteady potential. The steady flow and small disturbance unsteady flow equations are discretized using bilinear quadrilateral isoparametric finite elements. The nonlinear mean flow solution and streamline computational grid are computed simultaneously using Newton iteration. At each step of the Newton iteration, LU decomposition is used to solve the resulting set of linear equations. The unsteady flow problem is linear, and is also solved using LU decomposition. Next, a sensitivity analysis is performed to determine the effect small changes in cascade and airfoil geometry have on the mean and unsteady flow fields. The sensitivity analysis makes use of the nominal steady and unsteady flow LU decompositions so that no additional matrices need to be factored. Hence, the present method is computationally very efficient. To demonstrate how the sensitivity analysis may be used to redesign cascades, a compressor is redesigned for improved aeroelastic stability and two different fan exit guide vanes are redesigned for reduced downstream radiated noise. In addition, a framework detailing how the two-dimensional version of the method may be used to redesign three-dimensional geometries is presented.

Lorence, Christopher B.↗

NASA NTRS1984

The crater density of a linear feature

The standard method of crater counting presently employed permits determination of the relative ages of areal geologic features on cratered planets and satellites. A new method is presented which calculates the areal density of the craters superposed upon linear features. The method produces an effective synthetic area around the linear feature for each crater bin size, permitting the line counts to be compared directly with standard area counts. The shape of the synthetic area produced by the method is an oval, with the calculation of this area for each in being dependent on the median crater size for that bin, as well as the length and width of the linear feature. Incorporated into the method's equation is a factor to normalize the count to one million square kilometers for each bin. Because of the dependence of the method on the crater sizes, the largest possible number of bins should be counted for the method to be the most effective.

Spong, L.↗

NASA NTRS2023

Usability Evaluation Technical Report: A Usability Evaluation of the Communication Tool Prototype Prepared in Adobe XD for High Density Vertiplex (HDV) Team

Purpose: As cities get more crowded, the roadway infrastructure cannot keep up with the travel demands. Aviation can be a solution. Organizations supporting NASA’s Urban Air Mobility (UAM) concept are conducting studies on feasible concepts of operations for the new air traffic management system required to implement UAM. NASA’s High Density Vertiplex (HDV) team conducts studies using both live flights and simulations. Significance: The goal of this study is to evaluate the effectiveness of a communication tool developed for this team to enhance team situation awareness. Specifically, the study aims to determine if the communication tool is effective in allowing NASA-Langley and NASA-Ames team members playing different roles to understand what others are doing so that completed steps can be communicated to the whole team. Methods: This study asks participants to conduct four distinct scenarios using an Adobe XD prototype. Participants included five members of the target audience, who each received up to 15 minutes of training on how to use the Adobe XD prototype. The participants were NASA-Ames or NASA-Langley researchers, developers or tech leads who had 0 to 15 years of experience. Their ages ranged from 20 to 40 years old, and they were of both genders. After training, each participant was asked to do 4 tasks of varying difficulty, using the prototype testbed. After each task they reported on task difficulty and after all tasks were completed they were asked to fill out the System Usability Survey (SUS). Performance will be assessed through task completion time and accuracy rating. Results: User perception was evaluated by the SUS which showed a mean score over the six participants of 84. All the participants were above the acceptable range which was 70 and nobody was in the marginal usability zone between 50 and 70. In addition, most of the participants did well in terms of time on task, success rate and task difficulty on all the tasks except for Task 3. Participants had two minutes to do the task correctly to succeed at the task, otherwise it was marked as a failure. The geomean for Task 1 for completion time was 22.20; Task 2 was 58.28; Task 3 was 84.68 and Task 4 was 19.61. In Task 3, three out of six people failed this task but all five scored this with an average score 4.2 out of 5 for self-rating on task difficulty.

human factors↗