Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “precondition”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Modelling the Antarctic lower stratosphere

Results form modeling studies of the Antarctic lower stratosphere which have attempted to simulate the large springtime ozone losses and corresponding changes in other trace constituents are given. These studies were carried out in a photochemical box model, a one-dimensional model without transport and in a two-dimensional photochemical-dynamical-radiation model. The photochemical studies have investigated inter alia the sensitivity of ozone to inclusion in the model of heterogeneous chemistry, and to the inclusion of the ClO dimer. When both of these are incorporated in the model, ozone depletions resembling whose found in Halley Bay in 1987 (J.C. Farman, Nature, 329, 1987) can be reproduced. The temporal variations (both diurnal and during the August to October period) of a number of important tracers including HCl, ClONO2, OClO and BrO are discussed. The two-dimensional study concentrated on the difficulty of establishing in the model the dynamical preconditioning of the lower polar stratosphere - low temperatures, low N2O, etc., high ClOx. Calculations are presented to show: (1) the depletion of ozone during the springtime season, (2) the effect of large ozone losses on lower latitudes, and (3) the longer term (multi-year) variations of ozone in Antarctica, assuming realistic increases in the atmospheric halogen burden.

Chipperfield, M. P.↗
NASA NTRS1987

An investigation of the causes of failure of flexible thermal protection materials in an aerodynamic environment

Tests of small panels of advanced flexible reusable surface insulation (AFRSI) were conducted using a small wind tunnel that was designed to simulate Space Shuttle Orbiter entry mean-flow and pulsating aerodynamic loads. The wind tunnel, with a 3 inch wide by 1.75 inch high by 7.5 inch long test section, proved to be capable of continuous flow at dynamic pressures q near 580 psf with fluctuating pressures over 2 psi RMS at an excitation frequency f sub E of 200 Hz. For this investigation, however, the wind tunnel was used to test entry-temperature preconditioned and heat-cleaned AFRSI at q = 280 psf, Prms was nearly equal to 1.2 psi and f sub E = 200 Hz. The objective of these tests was to determine the mechanism of failure of AFRSI at Orbiter entry conditions. Details of the test apparatus and test results are presented.

Coe, Charles F.↗
NASA NTRS1989

Strain energy release rate as a function of temperature and preloading history utilizing the edge delamination fatique test method

Static laminate and tension-tension fatigue tests of IM7/8551-7 composite materials was performed. The Edge Delamination Test (EDT) was utilized to evaluate the temperature and preloading history effect on the critical strain energy release rate. Static and fatigue testing was performed at room temperature and 180 F (82 C). Three preloading schemes were used to precondition fatigue test specimens prior to performing the normal tension-tension fatigue EDT testing. Computer software was written to perform all fatigue testing while monitoring the dynamic modulus to detect the onset of delamination and record the test information for later retrieval and reduction.

Zimmerman, Richard S.↗
NASA NTRS1989

A genuinely multi-dimensional upwind cell-vertex scheme for the Euler equations

The solution of the two-dimensional Euler equations is based on the two-dimensional linear convection equation and the Euler-equation decomposition developed by Hirsch et al. The scheme is genuinely two-dimensional. At each iteration, the data are locally decomposed into four variables, allowing convection in appropriate directions. This is done via a cell-vertex scheme with a downwind-weighted distribution step. The scheme is conservative, and third-order accurate in space. The derivation and stability analysis of the scheme for the convection equation, and the derivation of the extension to the Euler equations are given. Preconditioning techniques based on local values of the convection speeds are discussed. The scheme for the Euler equations is applied to two channel-flow problems. It is shown to converge rapidly to a solution that agrees well with that of a third-order upwind solver.

Powell, Kenneth G.↗
NASA NTRS1988

On conjugate gradient type methods and polynomial preconditioners for a class of complex non-Hermitian matrices

Conjugate gradient type methods are considered for the solution of large linear systems Ax = b with complex coefficient matrices of the type A = T + i(sigma)I where T is Hermitian and sigma, a real scalar. Three different conjugate gradient type approaches with iterates defined by a minimal residual property, a Galerkin type condition, and an Euclidian error minimization, respectively, are investigated. In particular, numerically stable implementations based on the ideas behind Paige and Saunder's SYMMLQ and MINRES for real symmetric matrices are proposed. Error bounds for all three methods are derived. It is shown how the special shift structure of A can be preserved by using polynomial preconditioning. Results on the optimal choice of the polynomial preconditioner are given. Also, some numerical experiments for matrices arising from finite difference approximations to the complex Helmholtz equation are reported.

Freund, Roland↗
NASA NTRS1989

Parallel linear equation solvers for finite element computations

The overall objective of this research is to develop efficient methods for the solution of linear and nonlinear systems of equations on parallel and supercomputers, and to apply these methods to the solution of problems in structural analysis. Attention has been given so far only to linear equations. The methods considered for the solution of the stiffness equation Kx=f have been Choleski factorization and the conjugate gradient iteration with SSOR and Incomplete Choleski preconditioning. More detail on these methods will be given on subsequent slides. These methods have been used to solve for the static displacements for the mast and panel focus problems in conjunction with the CSM testbed system based on NICE/SPAR.

Ortega, James M.↗
NASA NTRS1988

An explicit Runge-Kutta method for 3D turbulent incompressible flows

A computer code has been developed to solve for the steady-state solution of the 3D incompressible Reynolds-averaged Navier-Stokes equations. The approach is based on the cell-center, central-difference, finite-volume formulation and an explicit one-step, multistage Runge-Kutta time-stepping scheme. The Baldwin-Lomax turbulence model is used. Techniques to accelerate the rate of convergence to a steady-state solution include the preconditioned method, the local time stepping, and the implicit residual smoothing. Improvements in computational efficiency have been demonstrated in several areas. This numerical procedure has been used to simulate the turbulent horseshoe vortex flow around an airfoil/flat-plate juncture.

Sung, Chao-Ho↗
NASA NTRS1989

A genuinely multi-dimensional upwind cell-vertex scheme for the Euler equations

A scheme of solving the two-dimensional Euler equations is developed. The scheme is genuinely two-dimensional. At each iteration, the data are locally decomposed into four variables, allowing convection in appropriate directions. This is done via a cell-vertex scheme with a downwind-weighted distribution step. The scheme is conservative and third-order accurate in space. The derivation and stability analysis of the scheme for the convection equation, and the derivation of the extension to the Euler equations are given. Preconditioning techniques based on local values of the convection speeds are discussed. The scheme for the Euler equations is applied to two channel-flow problems. It is shown to converge rapidly to a solution that agrees well with that of a third-order upwind solver.

Powell, Kenneth G.↗
NASA NTRS1989

Galerkin/Runge-Kutta discretizations for parabolic equations with time-dependent coefficients

A new class of fully discrete Galerkin/Runge-Kutta methods is constructed and analyzed for linear parabolic initial boundary value problems with time dependent coefficients. Unlike any classical counterpart, this class offers arbitrarily high order convergence while significantly avoiding what has been called order reduction. In support of this claim, error estimates are proved, and computational results are presented. Additionally, since the time stepping equations involve coefficient matrices changing at each time step, a preconditioned iterative technique is used to solve the linear systems only approximately. Nevertheless, the resulting algorithm is shown to preserve the original convergence rate while using only the order of work required by the base scheme applied to a linear parabolic problem with time independent coefficients. Furthermore, it is noted that special Runge-Kutta methods allow computations to be performed in parallel so that the final execution time can be reduced to that of a low order method.

Keeling, Stephen L.↗
NASA NTRS1989

Design of optimally smoothing multi-stage schemes for the Euler equations

In this paper, a method is developed for designing multi-stage schemes that give optimal damping of high-frequencies for a given spatial-differencing operator. The objective of the method is to design schemes that combine well with multi-grid acceleration. The schemes are tested on a nonlinear scalar equation, and compared to Runge-Kutta schemes with the maximum stable time-step. The optimally smoothing schemes perform better than the Runge-Kutta schemes, even on a single grid. The analysis is extended to the Euler equations in one space-dimension by use of 'characteristic time-stepping', which preconditions the equations, removing stiffness due to variations among characteristic speeds. Convergence rates independent of the number of cells in the finest grid are achieved for transonic flow with and without a shock. Characteristic time-stepping is shown to be preferable to local time-stepping, although use of the optimally damping schemes appears to enhance the performance of local time-stepping. The extension of the analysis to the two-dimensional Euler equations is hampered by the lack of a model for characteristic time-stepping in two dimensions. Some results for local time-stepping are presented.

Van Leer, Bram↗
NASA NTRS1989

Convergence acceleration of viscous and inviscid hypersonic flow calculations

The convergence of inviscid and viscous hypersonic flow calculations using a two-dimensional flux-splitting code is accelerated by applying a Richardson-type overrelaxation method. Successful results are presented for various cases; and a 50 percent savings in computer time is usually achieved. An analytical formula for the overrelaxation factor is derived, and the performance of this scheme is confirmed numerically. Moreover, application of this overrelaxation scheme produces a favorable preconditioning for Wynn's epsilon-algorithm. Both techniques have been extended to viscous three-dimensional flows and applied to accelerate the convergence of the compressible Navier-Stokes code. A savings of 40 percent in computer time is achieved in this case.

Cheer, A.↗
NASA NTRS1989

The stratospheric major warming of early December 1987

The NMC observations of temperature and geopotential carried out during the stratospheric major warming event of early December 1987 are analyzed. The December 1987 warming event displayed features typical of other major warmings, including vortex displacement, erosion, and splitting and, as such, was the earliest major warming yet observed in the Northern Hemisphere winter. The December 1987 major warming was unusual in that it took place during the (deep) westerly phase of the equatorial quasi-biennial oscillation and was not preceded by any significant preconditioning of the extratropical vortex. It is suggested that this unusually early major warming was mainly attributable to an anomalously large tropospheric forcing.

Baldwin, Mark P.↗
NASA NTRS1990

Concurrent Finite-Element Analysis On Hypercube Computers

Improved approach to execution of finite-element codes on hypercube and similar concurrent data processors increases efficiency of computation for many different types of problems. Based on flexible general model of computation on, and communication among, large-node parallel processors. Hybrid combining direct methods within subdomains and preconditioned-conjugate-gradient (PCG) iteration on remaining boundary system to obtain method of solution both robust and efficient. Preserves general structure and function of conventional sequential finite-element-method software. Generalized to both distributed- and shared-memory multicomputers, eliminating degree of machine specificity restricting general usefulness.

Lyzenga, G. A.↗
NASA NTRS1989

Gravity and animal embryos

Out of more than 4,500 rat hours in space there was only one experimental attempt (Cosmos 1129) at mating with an apparent absence of fertilization, implantation and subsequent development to term and partuition. Portions of this process were successfully flown, however, including the major portion of organogenesis in the rat (Cosmos 1524). These observations show that the cellular and molecular events underlying morphogenesis and differentiation in a small mammal can proceed normally in-utero under microgravity and other conditions encountered during short-duration flight. However, it is not known whether this situation will hold for larger mammals over several generations during extended missions that venture outside of near Earth. Furthermore, it is not understood why the previous attempt at obtaining copulation, fertilization and implantation in orbit failed but may be related to limitations of the rat habitat for meeting the preconditions for reproductive behavior. With respect to mammalian development it is important to appreciate that fertilization and development occur internally within the female and take a long time to complete and their success will, therefore, be contingent upon the maternal response to the space environment. One process central to development (the establishment of cell lines) is initiated prior to implantation by environmental asymmetries preceived by progenitor cells. These asymmetries appear to result from the formation of asymmetric cell-cell contacts and the concommitant development of an electrical axis across the progenitor cells. Other asymmetries were also documented. It is not known whether any of the known asymmetries perceived by progenitor cells are influenced by gravity vectors and/or by the maternal response to microgravity and other conditions encountered in space.

Wiley, Lynn M.↗
NASA NTRS1990

Application of the TRANAIR rectangular grid approach to the aerodynamic analysis of complex configurations

A numerical method is described which uses a rectangular grid to solve the nonlinear full potential equation about complex configurations. The grid is locally refined to resolve high velocity gradients arising from leading edge expansions or shock waves. The grid penetrates the boundary (described by networks of quadrilateral panels) and is generated automatically. Discrete operators are constructed using the finite element method. The system of nonlinear discrete equations is solved iteratively using a Krylov subspace method preconditioned by an exterior Poisson solver and a direct sparse solver. The primary emphasis is to provide design engineers with an aerodynamic analysis tool (the TRANAIR code) which is accurate, reliable, economical, and flexible to use. Computational results for many interesting configurations are presented.

Johnson, Forrester T.↗
NASA NTRS1989

Parallel conjugate gradient algorithms for manipulator dynamic simulation

Parallel conjugate gradient algorithms for the computation of multibody dynamics are developed for the specialized case of a robot manipulator. For an n-dimensional positive-definite linear system, the Classical Conjugate Gradient (CCG) algorithms are guaranteed to converge in n iterations, each with a computation cost of O(n); this leads to a total computational cost of O(n sq) on a serial processor. A conjugate gradient algorithms is presented that provide greater efficiency using a preconditioner, which reduces the number of iterations required, and by exploiting parallelism, which reduces the cost of each iteration. Two Preconditioned Conjugate Gradient (PCG) algorithms are proposed which respectively use a diagonal and a tridiagonal matrix, composed of the diagonal and tridiagonal elements of the mass matrix, as preconditioners. Parallel algorithms are developed to compute the preconditioners and their inversions in O(log sub 2 n) steps using n processors. A parallel algorithm is also presented which, on the same architecture, achieves the computational time of O(log sub 2 n) for each iteration. Simulation results for a seven degree-of-freedom manipulator are presented. Variants of the proposed algorithms are also developed which can be efficiently implemented on the Robot Mathematics Processor (RMP).

Fijany, Amir↗
NASA NTRS1990

Phase compatibilities of YBa2Cu3O(9-delta) type structure in quintenary systems Y-Ba-Cu-O-X (impurity)

Electrical transport properties of the oxidic high T(sub c) superconductors are significantly affected by the presence of minor amounts of various elements adventing as impurities, e.g., from the chemical environment during manufacturing. YBa2Cu3O(9-delta) is prone to an extinction of the superconductivity on (partial) substitution of all four elemental components. E.g., Pr (for Y), La (for Ba), Zn (for Cu) or peroxygroup (for O) substituents will alter some of the superconductivity preconditions, like mixed valence state in Cu3O7/O(9-delta) network or structural distortion of the network. Although various pseudoternary chemical equilibrium phase diagrams of the Y(O)-Ba(O)-Cu(O) system now are available, no consensus is generally shown, however, this is partly due to lack of compatible definitions of the equilibrium conditions. Less information is available about the phase compatibilities in the appropriate quaternary phase diagram (including oxygen) and virtually no information exists about any pentenary phase diagrams (including one impurity). Unfortunately, complexity of such systems, stemming both from number of quaternary or pentenary compounds and from visualizing the five-component phase system, limits this presentation to more or less close surroundings of the YBa2Cu3O(9-delta) type phase in appropriate pseudoquaternary or pseudopseudoternary diagrams, involving Y-Ba-Cu and O, O-CO2, alkaline metals, Mg and alkaline earths, and Sc and most of the 3-d and 4-f elements. The systems were investigated by means of x ray diffraction, neutron diffraction and chemical analytical methods on samples prepared by sol-gel technique from citrates. The superconductivity was characterized by measuring the diamagnetic susceptibility by SQUID.

Karen, P.↗
NASA NTRS1989

Efficient conjugate gradient algorithms for computation of the manipulator forward dynamics

The applicability of conjugate gradient algorithms for computation of the manipulator forward dynamics is investigated. The redundancies in the previously proposed conjugate gradient algorithm are analyzed. A new version is developed which, by avoiding these redundancies, achieves a significantly greater efficiency. A preconditioned conjugate gradient algorithm is also presented. A diagonal matrix whose elements are the diagonal elements of the inertia matrix is proposed as the preconditioner. In order to increase the computational efficiency, an algorithm is developed which exploits the synergism between the computation of the diagonal elements of the inertia matrix and that required by the conjugate gradient algorithm.

Fijany, Amir↗