Search NASA⌕ Search

SEARCH · Search NASA

Results for “SOR”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Analysis of the SOR iteration for the 9-point Laplacian

The SOR iteration for solving linear systems of equations depends upon an overrelaxation factor omega. A theory for determining omega was given by Young (1950) for consistently ordered matrices. Here the optimal omega is determined for the 9-point stencil for the model problem of Laplace's equation on a square. Several orderings of the equations are considered, including the natural rowwise and multicolor orderings, all of which lead to non-consistently ordered matrices, and two equivalence classes of orderings are found with different convergence behavior and optimal omega's. The results for the natural rowwise ordering are compared to those of Garabedian (1956) and it is explained why both results are, in a sense, correct, even though they differ. Also analyzed is a pseudo SOR method for the model problem and it is shown that it is not as effective as the SOR methods. Finally, the point SOR methods are compared to known results for line SOR methods for this problem.

Loyce M. Adams↗
NASA NTRS1993

Optimal convolution SOR acceleration of waveform relaxation with application to semiconductor device simulation

In this paper we describe a novel generalized SOR (successive overrelaxation) algorithm for accelerating the convergence of the dynamic iteration method known as waveform relaxation. A new convolution SOR algorithm is presented, along with a theorem for determining the optimal convolution SOR parameter. Both analytic and experimental results are given to demonstrate that the convergence of the convolution SOR algorithm is substantially faster than that of the more obvious frequency-independent waveform SOR algorithm. Finally, to demonstrate the general applicability of this new method, it is used to solve the differential-algebraic system generated by spatial discretization of the time-dependent semiconductor device equations.

Reichelt, Mark↗
NASA NTRS1986

Fourier analysis of the SOR iteration

The SOR iteration for solving linear systems of equations depends upon an overrelaxation factor omega. It is shown that for the standard model problem of Poisson's equation on a rectangle, the optimal omega and corresponding convergence rate can be rigorously obtained by Fourier analysis. The trick is to tilt the space-time grid so that the SOR stencil becomes symmetrical. The tilted grid also gives insight into the relation between convergence rates of several variants.

Leveque, R. J.↗
NASA NTRS1987

SOR and MGR(nu) experiments on the Crystal multicomputer

This report describes distributed implementations of the red/black SOR algorithm and of the MGR(nu) multigrid algorithm on the Crystal multicomputer. Rates of convergence and observed efficiencies for both algorithms are compared.

Kamowitz, David↗
NASA NTRS1997

Comprehensive Analyses of Data Collected from TEREK (Solar EUV Telescope) RES-C (Solar X-Ray Spectrometer) and SORS (Solar Radio Spectrometer) on board CORONAS-1 Using Magnetohydrodynamic Models

By using the observed magnetic field data obtained from the Wilcox Solar Observatory at Stanford University as the inputs to a two-dimensional plane-of-sky magnetohydrodynamic model, via numerical relaxation method, we have deduced the plasma and magnetic field parameters for the observed coronal hole by CORONAS-1. The method for this self-consistent MHD analysis will be discussed in detail. Numerical results for the magnetic field configuration, velocity distribution, density and temperature distributions will be presented. We have converted the computed density to polarization brightness in order to directly compare the MHD outputs with observations. Also included is a summary of achievements made during the grant period. This section is summarized into three categories: 1) Visit of Co-Investigators; 2) Presentations; and 3) Papers published, accepted and submitted for publication in journals.

Wu, S. T.↗
NASA NTRS1981

Numerical algorithms for finite element computations on arrays of microprocessors

The development of a multicolored successive over relaxation (SOR) program for the finite element machine is discussed. The multicolored SOR method uses a generalization of the classical Red/Black grid point ordering for the SOR method. These multicolored orderings have the advantage of allowing the SOR method to be implemented as a Jacobi method, which is ideal for arrays of processors, but still enjoy the greater rate of convergence of the SOR method. The program solves a general second order self adjoint elliptic problem on a square region with Dirichlet boundary conditions, discretized by quadratic elements on triangular regions. For this general problem and discretization, six colors are necessary for the multicolored method to operate efficiently. The specific problem that was solved using the six color program was Poisson's equation; for Poisson's equation, three colors are necessary but six may be used. In general, the number of colors needed is a function of the differential equation, the region and boundary conditions, and the particular finite element used for the discretization.

Ortega, J. M.↗
NASA NTRS1980

Application of the multigrid method to grid generation

The multigrid method (MGM), used to numerically solve the pair of nonlinear elliptic equations commonly used to generate two dimensional boundary-fitted coordinate systems is discussed. Two different geometries are considered: one involving a coordinate system fitted about a circle and the other selected for an impinging jet flow problem. Two different relaxation schemes are tried: one is successive point overrelaxation and the other is a four-color scheme vectorizeable to take advantage of a parallel processor computer for greater computational speed. Results using MGM are compared with those using SOR (doing successive overrelaxations with the corresponding relaxation scheme on the fine grid only). It is found that MGM becomes significantly more effective than SOR as more accuracy is demanded and as more corrective grids, or more grid points, are used. For the accuracy required, it is found that MGM is two to three times faster than SOR in computing time. With the four-color relaxation scheme as applied to the impinging jet problem, the advantage of MGM over SOR is not as great. This may be due to the effect of a poor initial guess on MGM for this problem.

Ohring, S.↗
NASA NTRS2024

Short Exploration Extravehicular Mobility Unit Testing Setup: Evaluation Under Realistic Pressure and Thermal Conditions

The purpose of Short Exploration Extravehicular Mobility Unit (SxEMU) thermal vacuum testing was to verify the functionality of the Design Verification Testing (DVT) prototype xEMU (SxEMU for this test) at vacuum pressures and extreme space and lunar surface temperature conditions. The SxEMU Thermal Vacuum Test was the culmination of the DVT xEMU project. This paper’s focus is on the pre-Extravehicular Activity (EVA) test setup, and general performance of the SxEMU Portable Life Support Subsystem (xPLSS), with focus on the performance of the Primary Oxygen Assembly (POA) and Secondary Oxygen Assembly (SOA), including Secondary Oxygen Regulator (SOR) takeover and the POA and SOA low-setpoint change inhibit. The initial pre-EVA test preparation included recharging the batteries and replenishing consumables, including test-system water, oxygen assemblies (with gaseous nitrogen), and the integrated thermal loops, including the Feedwater Supply Assemblies. xPLSS functionality testing included carbon dioxide (CO2) removal via the Rapid Cycle Amine swingbed system, thermal loop temperature control, and monitoring of suit ventilation loop pressure, temperature, and CO2 percentages. Testing evaluated automatic takeover of suit pressure control by the SOR after the primary oxygen supply is depleted. The Primary Oxygen Regulator and SOR low-setpoint change inhibit function prevents the crewmember from inadvertently setting the primary regulator to a low pressure setpoint during an EVA.

xPLSS↗
NASA NTRS1979

Vector processor algorithms for transonic flow calculations

This paper discusses a number of algorithms for solving the transonic full-potential equation in conservative form on a vector computer, such as the CDC STAR-100 or the CRAY-1. Recent research with the 'artificial density' method for transonics has led to development of some new iteration schemes which take advantage of vector-computer architecture without suffering significant loss of convergence rate. Several of these more promising schemes are described and 2-D and 3-D results are shown comparing the computational rates on the STAR and CRAY vector computers, and the CYBER-175 serial computer. Schemes included are: (1) Checkerboard SOR, (2) Checkerboard Leapfrog, (3) odd-even vertical line SOR, and (4) odd-even horizontal line SOR.

South, J. C., Jr.↗
NASA NTRS1981

Approximate factorization for incompressible flow

For computational solution of the incompressible Navier-Stokes equations, the approximate factorization (AF) algorithm is used to solve the vectorized momentum equation in delta form based on the pressure calculated in the previous time step. The newly calculated velocities are substituted into the pressure equation (obtained from a linear combination of the continuity and momentum equation), which is then solved by means of line SOR. Computational results are presented for the NACA 66 sub 3 018 airfoil at Reynolds numbers of 1000 and 40,000 and attack angles of 0 and 6 degrees. Comparison with wind tunnel data for Re = 40,000 indicates good qualitative agreement between measured and calculated pressure distributions. Quantitative agreement is only fair, however, with the calculations somewhat displaced from the measurements. Furthermore, the computed velocity profiles are unrealistically thick around the airfoil, due to the excessive amount of artificial viscosity needed for stability. Based on the performance of the algorithm with regard to stability, it is concluded that AF/SOR is suitable for calculations at Reynolds numbers less than 10,000. Speedwise, the method is faster than point SOR by at least a factor of two.

Bernard, R. S.↗
NASA NTRS1983

Multi-grid solution of Neumann pressure problem for viscous flows using primitive variables

The multi-grid (MG) technique has been advanced for use with the Neumann boundary-value problem in clustered curvilinear orthogonal coordinates. This comprises an important step in the analysis of viscous flows using the velocity-pressure formulation of the Navier-Stokes equations. With successive over-relaxation (SOR) as the smoothing operator and with suitably formulated restriction and coarse-grid-correction operators, a 4-grid procedure enhances the efficiency of fine-grid solutions of the Neumann problem by a factor of 3 to 5, depending on the problem parameters. Thy influence of the smoothing operator is also examined by employing the alternating-direction implicit and the strongly implicit techniques instead of SOR.

Ghia, U.↗
NASA NTRS1984

A vectorized solution for incompressible flow

An algorithm is developed to obtain solutions to the unsteady Reynolds-averaged incompressible Navier-Stokes equations in general curvilinear coordinates on a vector processor. The governing equations are in nonconservative form with the velocity and pressure as dependent variables. Two momentum equations and the Poisson equation for pressure form a set of three governing equations for three flow field unknowns: u, v, and p. The governing equations and boundary conditions are expressed in terms of boundary-conforming curvilinear coordinates, and a checkerboard SOR iteration is used to solve the governing equations. Several possible sequences for a checkerboard SOR iteration are investigated for finding the best overall convergence rate. The efficiency and capability of the present algorithm was assessed using the example of an 18 percent thick NACA 66(3)018 airfoil at zero degree angle of attack for chord Reynolds number range 1000-40,000.

Patel, N. R.↗
NASA NTRS1986

Transformation of two and three-dimensional regions by elliptic systems

Grid smoothing and orthogonalization procedures were developed and implemented in the construction of two and three dimensional grids. The procedures are based on the variational methods of grid generation. The two-dimensional examples were computed using the MSU IRIS Graphics Workstation. It was demonstrated that the elliptic grid generation equations, with arbitrary forcing functions, can be solved, in their variational formulation, using a gradient method. Since gradient methods have a global convergence property, the divergence problems often encountered when using SOR iterative methods can be avoided. It is not to be concluded, however, that SOR methods should be abandoned, since gradient methods tend to converge very slowly. In fact, slow convergence was the major problem encountered in the three-dimensional grids. Further progress was made on the continuing effort to develop conservative interpolation formulas for overlapping grids.

Mastin, C. Wayne↗
NASA NTRS1988

On the convergence of difference approximations to scalar conservation laws

A unified treatment is given for time-explicit, two-level, second-order-resolution (SOR), total-variation-diminishing (TVD) approximations to scalar conservation laws. The schemes are assumed only to have conservation form and incremental form. A modified flux and a viscosity coefficient are introduced to obtain results in terms of the latter. The existence of a cell entropy inequality is discussed, and such an equality for all entropies is shown to imply that the scheme is an E scheme on monotone (actually more general) data, hence at most only first-order accurate in general. Convergence for TVD-SOR schemes approximating convex or concave conservation laws is shown by enforcing a single discrete entropy inequality.

Osher, Stanley↗
NASA NTRS2003

A Sounding Rocket Investigation of the Fine Structure in the Mesopause Region in Conjunction with High Spatial Resolution Lidar Measurements

The grant funds were provided to carry out chemical tracer wind and turbulence measurements in the mesosphere and lower thermosphere as part of the Turbulent Oxygen Mixing Experiment (TOMEX) for which Dr. James Hecht of the Aerospace Corp. was the Principal Investigator. Clemson University designed, built, and tested two chemical tracer release payload sections for the 21.126 and 21.127 payloads which, in addition, had photometer and ionization gauge instrumentation. The tracer chemical was trimethyl aluminum (TMA). The experiment was carried out at White Sands Missile Range in New Mexico on October 26,2000. The location was chosen because of the proximity to the Starfire Optical Range (SOR) which hosted the powerful University of Illinois sodium lidar for an extended period prior to and also during the launch window. Since the SOR telescope is fully steerable, lidar measurements in the same volume sampled by the rocket were possible. The primary objective of the experiment was to measure the turbulent diffusion and mixing in the mesosphere and lower thermosphere, especially in layers characterized by convective and/or dynamical instabilities. The lidar, which is capable of measuring the sodium density, temperatures, and winds with good range and time resolution, provided the launch criteria, as well as context measurements for the in situ rocket observations.

Larsen, M. F.↗
NASA NTRS1997

Results of the Compensated Earth-Moon-Earth Retroreflector Laser Link (CEMERLL) Experiment

Adaptive optics techniques can be used to realize a robust low bit-error-rate link by mitigating the atmosphere-induced signal fades in optical communications links between ground-based transmitters and deep-space probes. Phase I of the Compensated Earth-Moon-Earth Retroreflector Laser Link (CEMERLL) experiment demonstrated the first propagation of an atmosphere-compensated laser beam to the lunar retroreflectors. A 1.06-micron Nd:YAG laser beam was propagated through the full aperture of the 1.5-m telescope at the Starfire Optical Range (SOR), Kirtland Air Force Base, New Mexico, to the Apollo 15 retroreflector array at Hadley Rille. Laser guide-star adaptive optics were used to compensate turbulence-induced aberrations across the transmitter's 1.5-m aperture. A 3.5-m telescope, also located at the SOR, was used as a receiver for detecting the return signals. JPL-supplied Chebyshev polynomials of the retroreflector locations were used to develop tracking algorithms for the telescopes. At times we observed in excess of 100 photons returned from a single pulse when the outgoing beam from the 1.5-m telescope was corrected by the adaptive optics system. No returns were detected when the outgoing beam was uncompensated. The experiment was conducted from March through September 1994, during the first or last quarter of the Moon.

Wilson, K. E.↗
NASA NTRS1976

LSI arrays for space stations

Two approaches have been taken to study CCD's and some of their fundamental limitations. First a numerical analysis approach has been developed to solve the coupled transport and Poisson's equation for a thorough analysis of charge transfer in a CCD structure. The approach is formulated by treating the minority carriers as a surface distribution at the Si-SiO2 interface and setting up coupled difference equations for the charge and the potential. The SOR method is proposed for solving the two dimensional Poisson's equation for the potential. Methods are suggested for handling the discontinuities to improve convergence. Second, CCD shift registers were fabricated with parameters which should allow complete charge transfer independent of the transfer electrode gap width. A test instrument was designed and constructed which can be used to test this, or any similar, three phase CCD shift register.

Gassaway, J. D.↗
NASA NTRS1977

Numerical solution of the compressible Navier-Stokes equations using density gradients as additional dependent variables

Numerical solution of two dimensional, time dependent, compressible viscous Navier-Stokes equations about arbitrary bodies was treated using density gradients as additional dependent variables. Thus, six dependent variables were computed with the SOR iteration method. Besides formulation for pressure gradient terms, a formulation for computing the body density was presented. To approximate the governing equations, an implicit finite difference method was employed. In computing the solution for the flow about a circular cylinder, a problem arose near the wall at both stagnation points. Thus, computations with various conditions were tried to examine the problem. Also, computations with and without formulations are compared. The flow variables were computed on 37 by 40 field first, then on an 81 by 40 field.

Kwon, J. H.↗