Engineering Papers⌕ Search

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.

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

At least 433 records · Page 24

A Constant Pressure Bomb

This report describes a new optical method of unusual simplicity and of good accuracy suitable to study the kinetics of gaseous reactions. The device is the complement of the spherical bomb of constant volume, and extends the applicability of the relationship, pv=rt for gaseous equilibrium conditions, to the use of both factors p and v. The method substitutes for the mechanical complications of a manometer placed at some distance from the seat of reaction the possibility of allowing the radiant effects of reaction to record themselves directly upon a sensitive film. It is possible the device may be of use in the study of the photoelectric effects of radiation. The method makes possible a greater precision in the measurement of normal flame velocities than was previously possible. An approximate analysis shows that the increase of pressure and density ahead of the flame is negligible until the velocity of the flame approaches that of sound.

Stevens, F W↗

Human factors process failure modes and effects analysis (HF PFMEA) software tool

Methods, computer-readable media, and systems for automatically performing Human Factors Process Failure Modes and Effects Analysis for a process are provided. At least one task involved in a process is identified, where the task includes at least one human activity. The human activity is described using at least one verb. A human error potentially resulting from the human activity is automatically identified, the human error is related to the verb used in describing the task. A likelihood of occurrence, detection, and correction of the human error is identified. The severity of the effect of the human error is identified. The likelihood of occurrence, and the severity of the risk of potential harm is identified. The risk of potential harm is compared with a risk threshold to identify the appropriateness of corrective measures.

Chandler, Faith T.↗

Stress-intensity factors for internal surface cracks in cylindrical pressure vessels

The stress intensity factors were calculated by a three dimensional finite element method. The finite element models employed singularity elements along the crack front and linear strain elements elsewhere. The models had about 6500 degrees of freedom. The stress intensity factors were evaluated from a nodal force method. An equation for the stress intensity factors was obtained form the results of the present analysis. The equation applies over a wide range of configuration parameters and was within about 5 percent of the present results. A comparison was made between the present results and other analyses of internal surface cracks in cylinders. The results from a boundary integral equation method were in agreement (+ or - 2 percent) and those from another finite element were in fair agreement (+ or - 8 percent) with the present results.

Newman, J. C., Jr.↗

Three-dimensional stress intensity factor analysis of a surface crack in a high-speed bearing

The boundary element method is applied to calculate the stress intensity factors of a surface crack in the rotating inner raceway of a high-speed roller bearing. The three-dimensional model consists of an axially stressed surface cracked plate subjected to a moving Hertzian contact loading. A multidomain formulation and singular crack-tip elements were employed to calculate the stress intensity factors accurately and efficiently for a wide range of configuration parameters. The results can provide the basis for crack growth calculations and fatigue life predictions of high-performance rolling element bearings that are used in aircraft engines.

Ballarini, Roberto↗

Determination of Parachute Joint Factors using Seam and Joint Testing

This paper details the methodology for determining the joint factor for all parachute components. This method has been successfully implemented on the Capsule Parachute Assembly System (CPAS) for the NASA Orion crew module for use in determining the margin of safety for each component under peak loads. Also discussed are concepts behind the joint factor and what drives the loss of material strength at joints. The joint factor is defined as a "loss in joint strength...relative to the basic material strength" that occurs when "textiles are connected to each other or to metals." During the CPAS engineering development phase, a conservative joint factor of 0.80 was assumed for each parachute component. In order to refine this factor and eliminate excess conservatism, a seam and joint testing program was implemented as part of the structural validation. This method split each of the parachute structural joints into discrete tensile tests designed to duplicate the loading of each joint. Breaking strength data collected from destructive pull testing was then used to calculate the joint factor in the form of an efficiency. Joint efficiency is the percentage of the base material strength that remains after degradation due to sewing or interaction with other components; it is used interchangeably with joint factor in this paper. Parachute materials vary in type-mainly cord, tape, webbing, and cloth -which require different test fixtures and joint sample construction methods. This paper defines guidelines for designing and testing samples based on materials and test goals. Using the test methodology and analysis approach detailed in this paper, the minimum joint factor for each parachute component can be formulated. The joint factors can then be used to calculate the design factor and margin of safety for that component, a critical part of the design verification process.

Mollmann, Catherine↗

Evolving Multi Rover Systems in Dynamic and Noisy Environments

In this chapter, we address how to evolve control strategies for a collective: a set of entities that collectively strives to maximize a global evaluation function that rates the performance of the full system. Addressing such problems by directly applying a global evolutionary algorithm to a population of collectives is unworkable because the search space is prohibitively large. Instead, we focus on evolving control policies for each member of the collective, where each member is trying to maximize the fitness of its own population. The main difficulty with this approach is creating fitness evaluation functions for the members of the collective that induce the collective to achieve high performance with respect to the system level goal. To overcome this difficulty, we derive member evaluation functions that are both aligned with the global evaluation function (ensuring that members trying to achieve high fitness results in a collective with high fitness) and sensitive to the fitness of each member (a member's fitness depends more on its own actions than on actions of other members). In a difficult rover coordination problem in dynamic and noisy environments, we show how to construct evaluation functions that lead to good collective behavior. The control policy evolved using aligned and member-sensitive evaluations outperforms global evaluation methods by up to a factor of four. in addition we show that the collective continues to perform well in the presence of high noise levels and when the environment is highly dynamic. More notably, in the presence of a larger number of rovers or rovers with noisy sensors, the improvements due to the proposed method become significantly more pronounced.

Tumer, Kagan↗

Comparative Mirror Cleaning Study: 'A Study on Removing Particulate Contamination'

The cleanliness of optical surfaces is recognized as an industry-wide concern for the performance of optical devices such as mirrors and telescopes, microscopes and lenses, lasers and interferometers, and prisms and optical filters. However, no standard has been established for optical cleaning and there is no standard definition of a 'clean' optical element. This study evaluates the effectiveness of commonly used optical cleaning techniques based on wafer configuration, contamination levels, and the number and size of removed particles. It is concluded that cleaning method and exposure time play a significant factor in obtaining a high removal percentage. The detergent bath and solvent rinse method displayed an increase in effective removal percentage as the contamination exposure increased. Likewise, CO2 snow cleaning showed a relatively consistent cleaning effectiveness. The results can help ensure mission success to flight projects developed for the NASA Origins Program. Advantages and disadvantages of each of the optical cleaning methods are described.

Houston, Karrie↗

Investigation of stress concentration factor - Stress intensity factor interaction for flaws in filleted rods

To predict the fatigue life of machine parts which contain flaws in regions of stress concentration, the engineer needs to treat a series of crack sizes and/or shapes quickly and economically. In this paper the stress intensity factor distributions for circular surface flaws in the roots of fillets in rods loaded by remote tension are presented and examined. The examination reveals insight into the interaction of cracks with other stress concentrating discontinuities, and based on this it is suggested that weight function methods used in conjunction with known crack solutions may provide needed, efficient stress intensity factor prediction methods for such three-dimensional geometries.

Springfield, C. W., Jr.↗

Structural Deterministic Safety Factors Selection Criteria and Verification

Though current deterministic safety factors are arbitrarily and unaccountably specified, its ratio is rooted in resistive and applied stress probability distributions. This study approached the deterministic method from a probabilistic concept leading to a more systematic and coherent philosophy and criterion for designing more uniform and reliable high-performance structures. The deterministic method was noted to consist of three safety factors: a standard deviation multiplier of the applied stress distribution; a K-factor for the A- or B-basis material ultimate stress; and the conventional safety factor to ensure that the applied stress does not operate in the inelastic zone of metallic materials. The conventional safety factor is specifically defined as the ratio of ultimate-to-yield stresses. A deterministic safety index of the combined safety factors was derived from which the corresponding reliability proved the deterministic method is not reliability sensitive. The bases for selecting safety factors are presented and verification requirements are discussed. The suggested deterministic approach is applicable to all NASA, DOD, and commercial high-performance structures under static stresses.

Verderaime, V.↗

Stress analysis and stress intensity factors for finite geometry solids containing rectangular surface cracks

The line method of analysis is applied to the Navier-Cauchy equations of elastic equilibrium to calculate the displacement field in a finite geometry bar containing a variable depth rectangular surface crack under extensionally applied uniform loading. The application of this method to these equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. Using the obtained displacement field, normal stresses and the stress intensity factor variation along the crack periphery are calculated for different crack depth to bar thickness ratios. Crack opening displacements and stress intensity factors are also obtained for a through-thickness, center cracked bar with variable thickness. The reported results show a considerable potential for using this method in calculating stress intensity factors for commonly encountered surface crack geometries in finite solids.

Gyekenyesi, J. P.↗

A finite element tensor approach to plate buckling and postbuckling.

A practical finite element method is presented for geometrically nonlinear (von Karman) plate problems. Symmetric strain energy tensors provide an efficient formulation and solution, and allow the effects of initial imperfections to be considered. An unbalanced force iteration is employed, with the advantage that the time consuming Gauss decomposition needs to be performed only once. The convergence is controlled by an automatically computed relaxation factor, and this makes the method applicable in advanced stages of non-linearity. Several different triangular elements are used in explicitly deriving the required tensors.

Vos, R. G.↗

Extending substructure based iterative solvers to multiple load and repeated analyses

Direct solvers currently dominate commercial finite element structural software, but do not scale well in the fine granularity regime targeted by emerging parallel processors. Substructure based iterative solvers--often called also domain decomposition algorithms--lend themselves better to parallel processing, but must overcome several obstacles before earning their place in general purpose structural analysis programs. One such obstacle is the solution of systems with many or repeated right hand sides. Such systems arise, for example, in multiple load static analyses and in implicit linear dynamics computations. Direct solvers are well-suited for these problems because after the system matrix has been factored, the multiple or repeated solutions can be obtained through relatively inexpensive forward and backward substitutions. On the other hand, iterative solvers in general are ill-suited for these problems because they often must restart from scratch for every different right hand side. In this paper, we present a methodology for extending the range of applications of domain decomposition methods to problems with multiple or repeated right hand sides. Basically, we formulate the overall problem as a series of minimization problems over K-orthogonal and supplementary subspaces, and tailor the preconditioned conjugate gradient algorithm to solve them efficiently. The resulting solution method is scalable, whereas direct factorization schemes and forward and backward substitution algorithms are not. We illustrate the proposed methodology with the solution of static and dynamic structural problems, and highlight its potential to outperform forward and backward substitutions on parallel computers. As an example, we show that for a linear structural dynamics problem with 11640 degrees of freedom, every time-step beyond time-step 15 is solved in a single iteration and consumes 1.0 second on a 32 processor iPSC-860 system; for the same problem and the same parallel processor, a pair of forward/backward substitutions at each step consumes 15.0 seconds.

Farhat, Charbel↗

Simplified methods for calculating photodissociation rates

Simplified methods for calculating the transmission of solar UV radiation and the dissociation coefficients of various molecules are compared. A significant difference sometimes appears in calculations of the individual band, but the total transmission and the total dissociation coefficients integrated over the entire SR (solar radiation) band region agree well between the methods. The ambiguities in the solar flux data affect the calculated dissociation coefficients more strongly than does the method. A simpler method is developed for the purpose of reducing the computation time and computer memory size necessary for storing coefficients of the equations. The new method can reduce the computation time by a factor of more than 3 and the memory size by a factor of more than 50 compared with the Hudson-Mahle method, and yet the result agrees within 10 percent (in most cases much less) with the original Hudson-Mahle results, except for H2O and CO2. A revised method is necessary for these two molecules, whose absorption cross sections change very rapidly over the SR band spectral range.

Shimazaki, T.↗

Assessment of Computational Fluid Dynamics (CFD) Models for Shock Boundary-Layer Interaction

A workshop on the computational fluid dynamics (CFD) prediction of shock boundary-layer interactions (SBLIs) was held at the 48th AIAA Aerospace Sciences Meeting. As part of the workshop numerous CFD analysts submitted solutions to four experimentally measured SBLIs. This paper describes the assessment of the CFD predictions. The assessment includes an uncertainty analysis of the experimental data, the definition of an error metric and the application of that metric to the CFD solutions. The CFD solutions provided very similar levels of error and in general it was difficult to discern clear trends in the data. For the Reynolds Averaged Navier-Stokes methods the choice of turbulence model appeared to be the largest factor in solution accuracy. Large-eddy simulation methods produced error levels similar to RANS methods but provided superior predictions of normal stresses.

DeBonis, James R.↗

Compressible turbulent flow simulation with a multigrid multiblock method

We describe a multigrid multiblock method for compressible turbulent flow simulations and present results obtained from calculations on a two-element airfoil. A vertex-based spatial discretization method and explicit multistage Runge-Kutta time-stepping are used. The slow convergence of a single grid method makes the multigrid method, which yields a speed factor of about 20, indispensable. The numerical predictions are in good agreement with experimental results. It is shown that the convergence of the multigrid process depends considerably on the ordering of the various loops. If the block loop is put inside the stage loop, the process converges more rapidly than if the block loop is situated outside the stage loop when a three-stage Runge-Kutta method is used. If a five-stage scheme is used, the process does not converge in the latter block ordering. Finally, the process based on the five-stage method is about 60 percent more efficient than with the three-stage method, if the block loop is inside the stage loop.

Kuerten, Hans↗

Effect of impregnation method on cycle life of the nickel electrode

The nickel electrode has been identified as the life limiting component for individual pressure vessel (IPV) nickel-hydrogen cells when cycled under a low earth orbit (LEO) cycle regime at deep depths of discharge. As a part of an overall program to develop a long life nickel electrode for nickel-hydrogen cells, the effect of two different methods of electrochemical impregnation on the cycle life of the nickel electrode was investigated. One method was the Pickett (aqueous/ethanolic) process. The other was the modified Bell (aqueous) process. The plaques for both impregnation methods were made by sintering dry carbonyl nickel powder in a reducing atmosphere. The plaques contain a nickel screen substrate. Electrodes made from both processes were cycle tested in Air Force design IPV nickel-hydrogen cells. The only factor different for this test was the method of plaque impregnation; all other factors were the same. The cells were cycled to failure under a 90 min LEO cycle regime at a deep depth of discharge (80 percent DOD). Failure for this test was defined to occur when the cell voltage degraded to 1.0 V prior to the completion of the 35 min discharge.

Smithrick, J. J.↗

Effect of impregnation method on cycle life of the nickel electrode

The nickel electrode has been identified as the life limiting component for individual pressure vessel (IPV) nickel-hydrogen cells when cycled under a low earth orbit (LEO) cycle regime at deep depths of discharge. As a part of an overall program to develop a long life nickel electrode for nickel-hydrogen cells, the effect of two different methods of electrochemical impregnation on the cycle life of the nickel electrode was investigated. One method was the Pickett (aqueous/ethanolic) process. The other was the modified Bell (aqueous) process. The plaques for both impregnation methods were made by sintering dry carbonyl nickel powder in a reducing atmosphere. The plaques contain a nickel screen substrate. Electrodes made from both processes were cycle tested in Air Force design IPV nickel-hydrogen cells. The only factor different for this test was the method of plaque impregnation; all other factors were the same. The cells were cycled to failure under a 90 min LEO cycle regime at a deep depth of discharge (80 percent DOD). Failure for this test was defined to occur when the cell voltage degraded to 1.0 V prior to the completion of the 35 min discharge.

Smithrick, J. J.↗

On the Dome Effect of Flux Radiometers to Radiative Forcing

Since the introduction of thermopile, pyranometers (solar, e.g., 0.3-3.0 microns) and pyrgeometers (terrestrial, e.g., 4-50 microns) have become instruments commonly used for measuring the broadband hemispherical irradiances at the surface in a long-term, monitoring mode for decades. These commercially available radiometers have been manufactured in several countries such as from the United States, Asia, and Europe, and are generally reliable and economical. These worldwide distributions of surface measurements become even more important in the era of Earth remote sensing in studying climate forcing. However, recent studies from field campaigns have pointed out that erroneous factors (e.g., temperature gradients between the filter dome and detector, emissivity of the thermopile) are responsible for the unacceptable level of uncertainty (e.g., 10-20 W/square Meter). Using a newly developed instrument of Quantum Well Infrared Photodetector (QWIP), we have characterized the brightness temperature fields of pyranometers and pyrgeometers under various sky conditions. The QWIP is based on the superlattice (GaAs/AlGaAs) technology and has a noise equivalent temperature (NE delta T) less than 0.1 K. The quality of pyranometer and pyrgeometer measurements can be improved largely by applying proper knowledge of the thermal parameters affecting the operation of the thermopile systems. For example, we show a method to determine the "dome factor" (the longwave emission divided by the longwave transmission of a pyrgeometer dome) from field measurements. The results show, and are verified independently by the QWIP, that our dome factors of 0.59 and 0.90 are much smaller than the value of 4.0 assumed by the WMO. Data correction procedure and algorithm will be presented and discussed.

Tsay, S.-C.↗