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 145 records · Page 8

A Practical Approach to Uncertainty Quantification Using Probability Boxes

To date, while the use of CFD for aerospace vehicle design and development is prevalent, the documentation of uncertainties associated with the simulations are rare. Instead, the current state-of-the-art relies heavily on the experience of the CFD practitioner to estimate the uncertainty associated with their simulations through simple sensitivity studies or subject matter expertise. This practice will have to be replaced with a formal uncertainty quantification (UQ) process if CFD is to play an expanded role in the research and engineering design community, test and evaluation community, and ultimately certification for flight. Accounting for uncertainties in a formal manner is a tedious process. Moreover, the typical CFD practitioner is not likely to be familiar with formal UQ methods. These factors have prevented the adoption of UQ methods in the engineering design and development cycle. This presentation will outline a credible approach to UQ using Probability Boxes that is straightforward to apply, and can readily be automated using existing UQ tool sets such as the DAKOTA packaged developed at Sandia. The added expense incurred when moving away from a deterministic CFD process to a stochastic one that captures uncertainties to enable risk-informed decision making will be discussed, as well as effective ways to reduce the computational costs.

Uncertainty Quantification↗

The factorization of large composite numbers on the MPP

The continued fraction method for factoring large integers (CFRAC) was an ideal algorithm to be implemented on a massively parallel computer such as the Massively Parallel Processor (MPP). After much effort, the first 60 digit number was factored on the MPP using about 6 1/2 hours of array time. Although this result added about 10 digits to the size number that could be factored using CFRAC on a serial machine, it was already badly beaten by the implementation of Davis and Holdridge on the CRAY-1 using the quadratic sieve, an algorithm which is clearly superior to CFRAC for large numbers. An algorithm is illustrated which is ideally suited to the single instruction multiple data (SIMD) massively parallel architecture and some of the modifications which were needed in order to make the parallel implementation effective and efficient are described.

Mckurdy, Kathy J.↗

Krylov methods preconditioned with incompletely factored matrices on the CM-2

The performance is measured of the components of the key interative kernel of a preconditioned Krylov space interative linear system solver. In some sense, these numbers can be regarded as best case timings for these kernels. Sweeps were timed over meshes, sparse triangular solves, and inner products on a large 3-D model problem over a cube shaped domain discretized with a seven point template. The performance of the CM-2 is highly dependent on the use of very specialized programs. These programs mapped a regular problem domain onto the processor topology in a careful manner and used the optimized local NEWS communications network. The rather dramatic deterioration in performance was documented when these ideal conditions no longer apply. A synthetic workload generator was developed to produce and solve a parameterized family of increasingly irregular problems.

Berryman, Harry↗

Analysis of Yttrium-Barium-Copper-Oxide by x ray diffraction and mechanical characterization

The efforts in developing high-temperature superconductor (HTSC) YBa2Cu3O7 electrical leads are to benefit future NASA missions that will carry payloads with sensitive instruments operating at cryogenic temperatures. Present-day leads made of copper or magnesium are responsible for as much as 50 percent of the parasitic heat load on cryogenic systems. A reduction of this load could be achieved by replacing the conventional materials with HTSC ceramic electrical leads. Superconductor quality has become a concern in the industry, as has the development of effective evaluation methods. The factors that need to be examined for these materials include material purity, mechanical properties, and superconducting ability below the critical temperature. We applied several methods to study these factors: thermogravimetric analysis, x-ray diffraction, tensile testing, and laser-generated ultrasound. Our objectives were to determine the average tensile strength and Young's modulus of the HTSC material and to compare them to those values for copper and manganin.

Arsenovic, Petar↗

A review of nonstatistical techniques for the estimation of vertical atmospheric structure from remote infrared measurements

Mathematical techniques used in atmospheric profile inversion are reviewed, with the discussion limited to methods that do not require extensive information on the statistical properties of the profiles. The basic problem of the calculation of vertical profiles from noise contaminated infrared spectral measurements is examined. Three specific examples of methods of solution are given: the iterative minimum estimation method, a second linear iterative method, and the algorithm of Chahine. Satellite data are used in an empirical comparison of the methods. Several factors are considered, including the stability of solutions, convergence behavior of the iterative methods, and dependence of the solutions on the first guess, a factor that is particularly important in the treatment of data from other planets.

Conrath, B. J.↗

TRASYS form factor matrix normalization

A method has been developed for adjusting a TRASYS enclosure form factor matrix to unity. This approach is not limited to closed geometries, and in fact, it is primarily intended for use with open geometries. The purpose of this approach is to prevent optimistic form factors to space. In this method, nodal form factor sums are calculated within 0.05 of unity using TRASYS, although deviations as large as 0.10 may be acceptable, and then, a process is employed to distribute the difference amongst the nodes. A specific example has been analyzed with this method, and a comparison was performed with a standard approach for calculating radiation conductors. In this comparison, hot and cold case temperatures were determined. Exterior nodes exhibited temperature differences as large as 7 C and 3 C for the hot and cold cases, respectively when compared with the standard approach, while interior nodes demonstrated temperature differences from 0 C to 5 C. These results indicate that temperature predictions can be artificially biased if the form factor computation error is lumped into the individual form factors to space.

Tsuyuki, Glenn T.↗

Optical methods of stress analysis applied to cracked components

After briefly describing the principles of frozen stress photoelastic and moire interferometric analyses, and the corresponding algorithms for converting optical data from each method into stress intensity factors (SIF), the methods are applied to the determination of crack shapes, SIF determination, crack closure displacement fields, and pre-crack damage mechanisms in typical aircraft component configurations.

Smith, C. W.↗

Evaluation of factors determining the accuracy of linearized subsonic panel methods

A systematic evaluation of the factors determining the accuracy of linearized subsonic panel methods is presented. In particular, the constant and quadratically varying doublet panel methods are compared for thin and thick surface modelings in two and three dimensions. The sensitivity of results to panel edge and control point locations is studied for both of the methods. The first order convergence of the quadratic doublet method near network edges and the subsequent effect on the Kutta condition is investigated. Results from a quadratic doublet method specifically designed for a vector processing computer are shown.

Thomas, J. L.↗

Analysis of corner cracks at hole by a 3-D weight function method with stresses from finite element method

Stress intensity factors for quarter-elliptical corner cracks emanating from a circular hole are determined using a 3-D weight function method combined with a 3-D finite element method. The 3-D finite element method is used to analyze uncracked configuration and provide stress distribution in the region where crack is to occur. Using this stress distribution as input, the 3-D weight function method is used to determine stress intensity factors. Three different loading conditions, i.e. remote tension, remote bending and wedge loading, are considered for a wide range in geometrical parameters. The significance in using 3-D uncracked stress distribution and the difference between single and double corner cracks are studied. Typical crack opening displacements are also provided. Comparisons are made with solutions available in the literature.

Zhao, W.↗

Correction factory techniques for improving aerodynamic prediction methods

A method for correcting discrete element lifting surface theory to reflect given experimental data is presented. Theoretical pressures are modified such that imposed constraints are satisfied while minimizing the changes to the pressures. Several types of correction procedures are presented and correlated; (1) scaling of pressures; (2) scaling of downwash values; and (3) addition of an increment to the downwash that is proportioned to pressure. Some special features are included in these methods and they include: (1) consideration of experimental data from multiple deflection modes, (2) limitation of the amplitudes of the correction factors, and (3) the use of correction factor mode shapes. These methods are correlated for cases involving all three Mach Number ranges using a FORTRAN IV computer program. Subsonically, a wing with an oscillating partial span control surface and a wing with a leading edge droop are presented. Transonically a two-dimensional airfoil with an oscillating flap is considered. Supersonically an arrow wing with and without camber is analyzed. In addition to correction factor methods an investigation is presented dealing with a new simplified transonic modification of the two-dimensional subsonic lifting surface theory. Correlations are presented for an airfoil with an oscillating flap.

Giesing, J. P.↗

Finite Element Analysis in Concurrent Processing: Computational Issues

The purpose of this research is to investigate the potential application of new methods for solving large-scale static structural problems on concurrent computers. It is well known that traditional single-processor computational speed will be limited by inherent physical limits. The only path to achieve higher computational speeds lies through concurrent processing. Traditional factorization solution methods for sparse matrices are ill suited for concurrent processing because the null entries get filled, leading to high communication and memory requirements. The research reported herein investigates alternatives to factorization that promise a greater potential to achieve high concurrent computing efficiency. Two methods, and their variants, based on direct energy minimization are studied: a) minimization of the strain energy using the displacement method formulation; b) constrained minimization of the complementary strain energy using the force method formulation. Initial results indicated that in the context of the direct energy minimization the displacement formulation experienced convergence and accuracy difficulties while the force formulation showed promising potential.

Sobieszczanski-Sobieski, Jaroslaw↗

Solution of elliptic partial differential equations by fast Poisson solvers using a local relaxation factor. 2: Two-step method

A two-step semidirect procedure is developed to accelerate the one-step procedure described in NASA TP-2529. For a set of constant coefficient model problems, the acceleration factor increases from 1 to 2 as the one-step procedure convergence rate decreases from + infinity to 0. It is also shown numerically that the two-step procedure can substantially accelerate the convergence of the numerical solution of many partial differential equations (PDE's) with variable coefficients.

Chang, S. C.↗

Solution of elliptic partial differential equations by fast Poisson solvers using a local relaxation factor. 1: One-step method

An algorithm for solving a large class of two- and three-dimensional nonseparable elliptic partial differential equations (PDE's) is developed and tested. It uses a modified D'Yakanov-Gunn iterative procedure in which the relaxation factor is grid-point dependent. It is easy to implement and applicable to a variety of boundary conditions. It is also computationally efficient, as indicated by the results of numerical comparisons with other established methods. Furthermore, the current algorithm has the advantage of possessing two important properties which the traditional iterative methods lack; that is: (1) the convergence rate is relatively insensitive to grid-cell size and aspect ratio, and (2) the convergence rate can be easily estimated by using the coefficient of the PDE being solved.

Chang, S. C.↗

Single Phase Passive Rectification Versus Active Rectification Applied to High Power Stirling Engines

Stirling engine converters are being considered as potential candidates for high power energy conversion systems required by future NASA explorations missions. These types of engines typically contain two major moving parts, the displacer and the piston, in which a linear alternator is attached to the piston to produce a single phase sinusoidal waveform at a specific electric frequency. Since all Stirling engines perform at low electrical frequencies (less or equal to 100 Hz), space explorations missions that will employ these engines will be required to use DC power management and distribution (PMAD) system instead of an AC PMAD system to save on space and weight. Therefore, to supply such DC power an AC to DC converter is connected to the Stirling engine. There are two types of AC to DC converters that can be employed, a passive full bridge diode rectifier and an active switching full bridge rectifier. Due to the inherent line inductance of the Stirling Engine-Linear Alternator (SE-LA), their sinusoidal voltage and current will be phase shifted producing a power factor below 1. In order to keep power the factor close to unity, both AC to DC converters topologies will implement power factor correction. This paper discusses these power factor correction methods as well as their impact on overall mass for exploration applications. Simulation results on both AC to DC converters topologies with power factor correction as a function of output power and SE-LA line inductance impedance are presented and compared.

Santiago, Walter↗

An extension of A-stability to alternating direction implicit methods

An alternating direction implicit (ADI) scheme was constructed by the method of approximate factorization. An A-stable linear multistep method (LMM) was used to integrate a model two-dimensional hyperbolic-parabolic partial differential equation. Sufficient conditions for the A-stability of the LMM were determined by applying the theory of positive real functions to reduce the stability analysis of the partial differential equations to a simple algebraic test. A linear test equation for partial differential equations is defined and then used to analyze the stability of approximate factorization schemes. An ADI method for the three-dimensional heat equation is also presented.

Warming, R. F.↗

A parallel-vector algorithm for rapid structural analysis on high-performance computers

A fast, accurate Choleski method for the solution of symmetric systems of linear equations is presented. This direct method is based on a variable-band storage scheme and takes advantage of column heights to reduce the number of operations in the Choleski factorization. The method employs parallel computation in the outermost DO-loop and vector computation via the 'loop unrolling' technique in the innermost DO-loop. The method avoids computations with zeros outside the column heights, and as an option, zeros inside the band. The close relationship between Choleski and Gauss elimination methods is examined. The minor changes required to convert the Choleski code to a Gauss code to solve non-positive-definite symmetric systems of equations are identified. The results for two large-scale structural analyses performed on supercomputers, demonstrate the accuracy and speed of the method.

Storaasli, Olaf O.↗