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 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS1979

Stress-intensity factors and crack-opening displacements for round compact specimens

A two dimensional, boundary collocation stress analysis was used to analyze various round compact specimens. The influence of the round external boundary and of pin-loaded holes on stress intensity factors and crack opening displacements was determined as a function of crack-length-to-specimen-width ratios. A wide-range equation for the stress intensity factors was developed. Equations for crack-surface displacements and load-point displacements were also developed. In addition, stress intensity factors were calculated from compliance methods to demonstrate that load-displacement records must be made at the loading points and not along the crack line for crack-length-to-specimen-width ratios less than about 0.4.

Newman, J. C., Jr.↗

NASA NTRS1982

Climatic influence of background and volcanic stratosphere aerosol models

A simple modelization of the earth atmosphere system including tropospheric and stratospheric aerosols has been derived and tested. Analytical expressions are obtained for the albedo variation due to a thin stratospheric aerosol layer. Also outlined are the physical procedures and the respective influence of the main parameters: aerosol optical thickness, single scattering albedo and asymmetry factor, and sublayer albedo. The method is applied to compute the variation of the zonal and planetary albedos due to a stratospheric layer of background H2SO4 particles and of volcanic ash.

Deschamps, P. Y.↗

NASA NTRS1986

Broadcasting satellite service synthesis using gradient and cyclic coordinate search procedures

Two search techniques are considered for solving satellite synthesis problems. Neither is likely to find a globally optimal solution. In order to determine which method performs better and what factors affect their performance, we design an experiment and solve the same problem under a variety of starting solution configuration-algorithm combinations. Since there is no randomization in the experiment, we present results of practical, rather than statistical, significance. Our implementation of a cyclic coordinate search procedure clearly finds better synthesis solutions than our implementation of a gradient search procedure does with our objective of maximizing the minimum C/I ratio computed at test points on the perimeters of the intended service areas. The length of the available orbital arc and the configuration of the starting solution are shown to affect the quality of the solutions found.

Reilly, C. H.↗

NASA NTRS1986

Inclusion of elastically connected members in exact buckling and frequency calculations

A standard stiffness matrix procedure which permits any combination of rigid, elastic, pinned or sliding connections of the degrees of freedom at the ends of a member to the nodes of its parent structure is described, in order to show how easily it can be extended to allow an existing algorithm to be used to ensure that no eigenvalues of the parent structure can be missed even when 'exact' member theory is used. The eigenvalues are the natural frequencies of undamped free vibration analyses or the critical load factors of buckling problems. The method preserves the exactness of the member theory and an efficient method for computer application is indicated. The theory also permits any combination of rigid, elastic, pinned or sliding connections between the freedoms of a substructure and those of its parent structure.

Williams, F. W.↗

NASA NTRS1986

Broadcasting satellite service synthesis using gradient and cyclic coordinate search procedures

Two search techniques are considered for solving satellite synthesis problems. Neither is likely to find a globally optimal solution. In order to determine which method performs better and what factors affect their performance, an experiment is designed and the same problem is solved under a variety of starting solution configuration-algorithm combinations. Since there is no randomization in the experiment, results of practical, rather than statistical, significance are presented. Implementation of a cyclic coordinate search procedure clearly finds better synthesis solutions than implementation of a gradient search procedure does with the objective of maximizing the minimum C/I ratio computed at test points on the perimeters of the intended service areas. The length of the available orbital arc and the configuration of the starting solution are shown to affect the quality of the solutions found.

Reilly, C. H.↗

NASA NTRS1986

Further development of the dynamic gas temperature measurement system

A compensated thermocouple measurement method was experimentally verified. Dynamic signal content from an atmospheric pressure laboratory burner was measured by the dynamic temperature sensor and a relatively delicate fine-wire resistance thermometer. Compensated data from the two dynamic temperature sensor thermoelements were compared with the compensated fine-wire data in the frequency domain. Absolute differences between spectral line amplitudes measured with different sensors are small relative to the mean temperature and verify the compenation method. Increases in precision of the measurement method require optimization of several factors, and directions for further work are identified.

Elmore, D. L.↗

NASA NTRS1986

Further development of the dynamic gas temperature measurement system. Volume 1: Technical efforts

A compensated dynamic gas temperature thermocouple measurement method was experimentally verified. Dynamic gas temperature signals from a flow passing through a chopped-wheel signal generator and an atmospheric pressure laboratory burner were measured by the dynamic temperature sensor and other fast-response sensors. Compensated data from dynamic temperature sensor thermoelements were compared with fast-response sensors. Results from the two experiments are presented as time-dependent waveforms and spectral plots. Comparisons between compensated dynamic temperature sensor spectra and a commercially available optical fiber thermometer compensated spectra were made for the atmospheric burner experiment. Increases in precision of the measurement method require optimization of several factors, and directions for further work are identified.

Elmore, D. L.↗

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.↗