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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Spectral ordering techniques for incomplete LU preconditoners for CG methods

The effectiveness of an incomplete LU (ILU) factorization as a preconditioner for the conjugate gradient method can be highly dependent on the ordering of the matrix rows during its creation. Detailed justification for two heuristics commonly used in matrix ordering for anisotropic problems is given. The bandwidth reduction and weak connection following heuristics are implemented through an ordering method based on eigenvector computations. This spectral ordering is shown to be a good representation of the heuristics. Analysis and test cases in two and three dimensional diffusion problems demonstrate when ordering is important, and when an ILU decomposition will be ordering insensitive. The applicability of the heuristics is thus evaluated and placed on a more rigorous footing.

Clift, Simon S.↗
NASA NTRS1974

Investigation of the ionospheric Faraday rotation for use in orbit corrections

The possibility of mapping the Faraday factors on a worldwide basis was examined as a simple method of representing the conversion factors for any possible user. However, this does not seem feasible. The complex relationship between the true magnetic coordinates and the geographic latitude, longitude, and azimuth angles eliminates the possibility of setting up some simple tables that would yield worldwide results of sufficient accuracy. Tabular results for specific stations can easily be produced or could be represented in graphic form.

Llewellyn, S. K.↗
NASA NTRS1989

An indirect method for numerical optimization using the Kreisselmeir-Steinhauser function

A technique is described for converting a constrained optimization problem into an unconstrained problem. The technique transforms one of more objective functions into reduced objective functions, which are analogous to goal constraints used in the goal programming method. These reduced objective functions are appended to the set of constraints and an envelope of the entire function set is computed using the Kreisselmeir-Steinhauser function. This envelope function is then searched for an unconstrained minimum. The technique may be categorized as a SUMT algorithm. Advantages of this approach are the use of unconstrained optimization methods to find a constrained minimum without the draw down factor typical of penalty function methods, and that the technique may be started from the feasible or infeasible design space. In multiobjective applications, the approach has the advantage of locating a compromise minimum design without the need to optimize for each individual objective function separately.

Wrenn, Gregory A.↗
NASA NTRS2017

Tortuosity Computations of Porous Materials using the Direct Simulation Monte Carlo

Low-density carbon fiber preforms, used as thermal protection systems (TPS) materials for planetary entry systems, have permeable, highly porous microstructures consisting of interlaced fibers. Internal gas transport in TPS is important in modeling the penetration of hot boundary-layer gases and the in-depth transport of pyrolysis and ablation products. The gas effective diffusion coefficient of a porous material must be known before the gas transport can be modeled in material response solvers; however, there are very little available data for rigid fibrous insulators used in heritage TPS.The tortuosity factor, which reflects the efficiency of the percolation paths, can be computed from the effective diffusion coefficient of a gas inside a porous material and is based on the micro-structure of the material. It is well known, that the tortuosity factor is a strong function of the Knudsen number. Due to the small characteristic scales of porous media used in TPS applications (typical pore size of the order of 50 micron), the transport of gases can occur in the rarefied and transitional regimes, at Knudsen numbers above 1. A proper way to model the gas dynamics at these conditions consists in solving the Boltzmann equation using particle-based methods that account for movement and collisions of atoms and molecules.In this work we adopt, for the first time, the Direct Simulation Monte Carlo (DSMC) method to compute the tortuosity factor of fibrous media in the rarefied regime. To enable realistic simulations of the actual transport of gases in the porous medium, digitized computational grids are obtained from X-ray micro-tomography imaging of real TPS materials. The SPARTA DSMC solver is used for simulations. Effective diffusion coefficients and tortuosity factors are obtained by computing the mean-square displacement of diffusing particles.We first apply the method to compute the tortuosity factors as a function of the Knudsen number for computationally designed materials such as random cylindrical fibers and packed bed of spheres with prescribed porosity. Results are compared to literature values obtained using random walk methods in the rarefied and transitional regime and a finite-volume method for the continuum regime. We then compute tortuosity factors for a real carbon fiber material with a transverse isotropic structure (FiberForm), quantifying differences between through-thickness and in-plain tortuosities at various Knudsen regimes.

Tortuosity↗
NASA NTRS1994

Stress intensity factors for surface and corner cracks emanating from a wedge-loaded hole

To assist analysis of riveted lap joints, stress intensity factors are determined for surface and corner cracks emanating from a wedge-loaded hole by using a 3-D weight function method in conjunction with a 3-D finite element method. A stress intensity factor equation for surface cracks is also developed to provide a closed-form solution. The equation covers commonly-encountered geometrical ranges and retains high accuracy over the entire range.

Zhao, W.↗
NASA NTRS1991

Combining global and local approximations

A method based on a linear approximation to a scaling factor, designated the 'global-local approximation' (GLA) method, is presented and shown capable of extending the range of usefulness of derivative-based approximations to a more refined model. The GLA approach refines the conventional scaling factor by means of a linearly varying, rather than constant, scaling factor. The capabilities of the method are demonstrated for a simple beam example with a crude and more refined FEM model.

Haftka, Raphael T.↗
NASA NTRS1972

The elastic analysis of the part-circular surface flaw problem by the alternating method.

This paper summarizes and evaluates the work done on the elastic analysis of the surface flaw problem by the application of the alternating method. An attempt is made to describe the alternating method and to present the history of its application to the surface flaw problem and to related problems in fracture mechanics. Stress intensity factors obtained by this method are summarized and compared. Results are also compared to those obtained by investigators using other methods of analysis. An evaluation of the use of the alternating method is presented with the purpose of pointing out the advantages and disadvantages in the application of this technique to the surface flaw problem.

Smith, F. W.↗
NASA NTRS1988

Methodology for matching experimental and computational aerodynamic data

Correction factor methodologies have been developed which use steady experimental or analytical pressure or force data to correct steady and unsteady aerodynamic calculations. Three methods of calculating correction factors have been developed to match steady surface pressure distributions, to match airfoil section forces and moments. Data for a rectangular supercritical wing that was previously tested in the NASA Langley Research Center Transonic Dynamics Tunnel have been used to determine correction factors to match surface pressure distributions for a range of Mach numbers.

Wieseman, Carol D.↗
NASA NTRS1988

Methodology for matching experimental and analytical aerodynamic data

Correction factor methodologies have been developed which use steady experimental or analytical pressure or force data to correct steady and unsteady aerodynamic calculations. Three methods of calculating correction factors have been developed to match steady surface pressure distributions, to match airfoil section forces and moments, and to match total forces and moments. Data for a rectangular supercritical wing that was previously tested in the NASA Langley Research Center Transonic Dynamics Tunnel have been used to determine correction factors to match surface pressure distributions for a range of Mach numbers.

Wieseman, Carol D.↗
NASA NTRS1996

Multigrid Acceleration of Time-Accurate DNS of Compressible Turbulent Flow

An efficient scheme for the direct numerical simulation of 3D transitional and developed turbulent flow is presented. Explicit and implicit time integration schemes for the compressible Navier-Stokes equations are compared. The nonlinear system resulting from the implicit time discretization is solved with an iterative method and accelerated by the application of a multigrid technique. Since we use central spatial discretizations and no artificial dissipation is added to the equations, the smoothing method is less effective than in the more traditional use of multigrid in steady-state calculations. Therefore, a special prolongation method is needed in order to obtain an effective multigrid method. This simulation scheme was studied in detail for compressible flow over a flat plate. In the laminar regime and in the first stages of turbulent flow the implicit method provides a speed-up of a factor 2 relative to the explicit method on a relatively coarse grid. At increased resolution this speed-up is enhanced correspondingly.

Broeze, Jan↗
NASA NTRS1995

Formal methods and their role in digital systems validation for airborne systems

This report is based on one prepared as a chapter for the FAA Digital Systems Validation Handbook (a guide to assist FAA certification specialists with advanced technology issues). Its purpose is to explain the use of formal methods in the specification and verification of software and hardware requirements, designs, and implementations; to identify the benefits, weaknesses, and difficulties in applying these methods to digital systems used in critical applications; and to suggest factors for consideration when formal methods are offered in support of certification. The presentation concentrates on the rationale for formal methods and on their contribution to assurance for critical applications within a context such as that provided by DO-178B (the guidelines for software used on board civil aircraft); it is intended as an introduction for those to whom these topics are new.

Rushby, John↗
NASA NTRS1974

Development of modified vibration test criteria for qualifying space vehicle components

Simplified methods are described to estimate the test criteria of primary structures at component attachment points subjected to broadband random acoustic excitations. The current method utilizes a constant smeared component mass attenuation factor across the frequency range of interest. The developed method indicates that the attenuation factor is based on a frequency dependent ratio of the mechanical impedances of both the component and primary structures. The procedures used to predict the structural responses are considered as the present state-of-the-art and provide satisfactory prediction results. Example problems are used to illustrate the application procedures of the two methods and to compare the significant difference. It was found that the lower test criteria obtained by the impedance ratio method is due to the results of considering the effects of component/primary structure interaction.

Chang, K. Y.↗
NASA NTRS1993

Matching Experimental And Analytical Aerodynamic Data

Three methods of calculating correction factors developed to match steady lifting distributions of pressure, forces and moments on airfoil sections, and total forces and moments. Data for rectangular supercritical wing tested in NASA Langley Research Center Transonic Dynamics Tunnel used to determine correction factors to match distributions of surface pressure for mach numbers ranging from 0.266 to 0.8. These correction factors also applied to calculations of unsteady aerodynamic flow, and comparisons made with oscillatory experimental data for range of reduced frequencies at several mach numbers. Methods of combining analytical and experimental data should lead to more-accurate predictions of performance of airfoils. Also applied to such airfoils as those in gas-turbine engines, fans, and wind turbines.

Wieseman, Carol D.↗
NASA NTRS2023

Single-Scattering Properties of Ellipsoidal Dust Aerosols Constrained By Measured Dust Shape Distributions

Most global aerosol models approximate dust as spherical particles, whereas most remote sensing retrieval algorithms approximate dust as spheroidal particles with a shape distribution that conflicts with measurements. These inconsistent and inaccurate shape assumptions generate biases in dust single-scattering properties. Here, we obtain dust single-scattering properties by approximating dust as triaxial ellipsoidal particles with observationally constrained shape distributions. We find that, relative to the ellipsoidal dust optics obtained here, the spherical dust optics used in most aerosol models underestimate dust single-scattering albedo, mass extinction efficiency, and asymmetry parameter for almost all dust sizes in both the shortwave and longwave spectra. We further find that the ellipsoidal dust optics are in substantially better agreement with observations of the scattering matrix and linear depolarization ratio than the spheroidal dust optics used in most retrieval algorithms. However, relative to observations, the ellipsoidal dust optics overestimate the lidar ratio by underestimating the backscattering intensity by a factor of ∼2. This occurs largely because the computational method used to simulate ellipsoidal dust optics (i.e., the improved geometric optics method) underestimates the backscattering intensity by a factor of ∼2 relative to other computational methods (e.g., the physical geometric optics method). We conclude that the ellipsoidal dust optics with observationally constrained shape distributions can help improve global aerosol models and possibly remote sensing retrieval algorithms that do not use the backscattering signal.

Dust↗
NASA NTRS2020

Farm-Level Risk Factors of Increased Abortion and Mortality in Domestic Ruminants during the 2010 Rift Valley Fever Outbreak in Central South Africa

Background:Rift Valley fever (RVF) outbreaks in domestic ruminants have severe socio-economic impacts. Climate-based continental predictions providing early warnings to regions at risk for RVF outbreaks are not of a high enough resolution for ruminant owners to assess their individual risk. (2) Methods: We analyzed risk factors for RVF occurrence and severity at the farm level using the number of domestic ruminant deaths and abortions reported by farmers in central South Africa during the 2010 RVF outbreaks using a Bayesian multinomial hurdle framework. (3) Results: We found strong support that the proportion of days with precipitation, the number of water sources, and the proportion of goats in the herd were positively associated with increased severity of RVF (the numbers of deaths and abortions). We did not find an association between any risk factors and whether RVF was reported on farms. (4) Conclusions: At the farm level we identified risk factors of RVF severity; however, there was little support for risk factors of RVF occurrence. The identification of farm-level risk factors for Rift Valley fever virus (RVFV) occurrence would support and potentially improve current prediction methods and would provide animal owners with critical information needed in order to assess their herd’s risk of RVFV infection.

Melinda K Rostel↗
NASA NTRS1988

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

A unified method to permute vector stored Upper triangular Diagonal factorized covariance and vector stored upper triangular Square Root Information arrays is presented. The method involves cyclic permutation of the rows and columns of the arrays and retriangularization with fast (slow) Givens rotations (reflections). Minimal computation is performed, and a one dimensional scratch array is required. To make the method efficient for large arrays on a virtual memory machine, computations are arranged so as to avoid expensive paging faults. This method is potentially important for processing large volumes of radio metric data in the Deep Space Network.

Muellerschoen, R. J.↗
NASA NTRS1994

Beyond maximum entropy: Fractal Pixon-based image reconstruction

We have developed a new Bayesian image reconstruction method that has been shown to be superior to the best implementations of other competing methods, including Goodness-of-Fit methods such as Least-Squares fitting and Lucy-Richardson reconstruction, as well as Maximum Entropy (ME) methods such as those embodied in the MEMSYS algorithms. Our new method is based on the concept of the pixon, the fundamental, indivisible unit of picture information. Use of the pixon concept provides an improved image model, resulting in an image prior which is superior to that of standard ME. Our past work has shown how uniform information content pixons can be used to develop a 'Super-ME' method in which entropy is maximized exactly. Recently, however, we have developed a superior pixon basis for the image, the Fractal Pixon Basis (FPB). Unlike the Uniform Pixon Basis (UPB) of our 'Super-ME' method, the FPB basis is selected by employing fractal dimensional concepts to assess the inherent structure in the image. The Fractal Pixon Basis results in the best image reconstructions to date, superior to both UPB and the best ME reconstructions. In this paper, we review the theory of the UPB and FPB pixon and apply our methodology to the reconstruction of far-infrared imaging of the galaxy M51. The results of our reconstruction are compared to published reconstructions of the same data using the Lucy-Richardson algorithm, the Maximum Correlation Method developed at IPAC, and the MEMSYS ME algorithms. The results show that our reconstructed image has a spatial resolution a factor of two better than best previous methods (and a factor of 20 finer than the width of the point response function), and detects sources two orders of magnitude fainter than other methods.

Puetter, Richard C.↗