Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gauss-Seidel”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Atmospheric transfer of radiation above an inhomogeneous non-Lambertian reflective ground. I - Theory

A method for solving the three-dimensional equation of transfer for a vertically inhomogeneous atmosphere bounded by a reflecting surface of non-uniform bidirectional reflectance is presented. The technique incorporates a two-dimensional spatial Fourier transform of the transfer equation and solution of the resulting expressions for each Fourier component of the radiation field using the method of Gauss-Seidel iteration. The intensity field in the spatial domain, which is calculated for a variety of altitudes, zenith angles, and azimuths, is reconstructed using the inverse Fourier transform. An empirical surface bidirectional reflectance function is employed, permitting the consideration of non-Lambertian reflective properties.

Diner, D. J.↗
NASA NTRS1985

Upwind relaxation algorithms for the Navier Stokes equations

The development of upwind relaxation algorithms for obtaining efficient steady-state solutions to the compressible Navier-Stokes equations is described. The method is second-order accurate spatially and naturally disipative, using third-order flux splitting of the pressure and convective terms and second-order central differencing for shear and heat flux terms. A line Gauss-Seidel relaxation approach, shown to be unconditionally stable for model convection and diffusion equations, is used. The algorithm is demonstrated for several flows using the thin-layer form of the equations, including the problem of shock-induced separation over a flat plate.

Thomas, J. L.↗
NASA NTRS1985

An efficient iteration strategy for the solution of the Euler equations

A line Gauss-Seidel (LGS) relaxation algorithm in conjunction with a one-parameter family of upwind discretizations of the Euler equations in two-dimensions is described. The basic algorithm has the property that convergence to the steady-state is quadratic for fully supersonic flows and linear otherwise. This is in contrast to the block ADI methods (either central or upwind differenced) and the upwind biased relaxation schemes, all of which converge linearly, independent of the flow regime. Moreover, the algorithm presented here is easily enhanced to detect regions of subsonic flow embedded in supersonic flow. This allows marching by lines in the supersonic regions, converging each line quadratically, and iterating in the subsonic regions, thus yielding a very efficient iteration strategy. Numerical results are presented for two-dimensional supersonic and transonic flows containing both oblique and normal shock waves which confirm the efficiency of the iteration strategy.

Walters, R. W.↗
NASA NTRS1985

Three dimensional CSCM method for the compressible Navier-Stokes equations with application to a multi-nozzle exhaust flowfield

The effective design of three-dimensional flow systems such as a clustered multinozzle shrouded exhaust system requires detailed analysis of the internal flowfield. It is pointed out that new numerical techniques on the most modern vector computers are making these analyses practicable. In this paper, the upwind differenced implicit conservative supracharacteristics method (CSCM) for three-dimensional viscous flow is presented. The new method is an implicit 'method of planes' symmetric Gauss-Seidel relaxation scheme. The data is conveniently stored on successive planes along the streamwise coordinate, and the system of equations is solved twice in each successive plane of the streamwise coordinate. The developed new three-dimensional compressible Navier-Stokes algorithm combines the best features of the storage and computationally efficient space marching procedures with the generality of time dependent techniques to solve flows with elliptic and streamwise separated flows.

Bardina, J.↗
NASA NTRS1985

Implicit flux-split schemes for the Euler equations

Recent progress in the development of implicit algorithms for the Euler equations using the flux-vector splitting method is described. Comparisons of the relative efficiency of relaxation and spatially-split approximately factored methods on a vector processor for two-dimensional flows are made. For transonic flows, the higher convergence rate per iteration of the Gauss-Seidel relaxation algorithms, which are only partially vectorizable, is amply compensated for by the faster computational rate per iteration of the approximately factored algorithm. For supersonic flows, the fully-upwind line-relaxation method is more efficient since the numerical domain of dependence is more closely matched to the physical domain of dependence. A hybrid three-dimensional algorithm using relaxation in one coordinate direction and approximate factorization in the cross-flow plane is developed and applied to a forebody shape at supersonic speeds and a swept, tapered wing at transonic speeds.

Thomas, J. L.↗
NASA NTRS1985

Numerical solution of Euler's equation by perturbed functionals

A perturbed functional iteration has been developed to solve nonlinear systems. It adds at each iteration level, unique perturbation parameters to nonlinear Gauss-Seidel iterates which enhances its convergence properties. As convergence is approached these parameters are damped out. Local linearization along the diagonal has been used to compute these parameters. The method requires no computation of Jacobian or factorization of matrices. Analysis of convergence depends on properties of certain contraction-type mappings, known as D-mappings. In this article, application of this method to solve an implicit finite difference approximation of Euler's equation is studied. Some representative results for the well known shock tube problem and compressible flows in a nozzle are given.

Dey, S. K.↗
NASA NTRS1986

A numerical and experimental investigation of electrochemical aircraft deicing

This study was composed of three parts. The first part involved the extension of an existing transient two dimensional numerical code for an electrothermal deicer so that it would simulate the situation where a variable thickness ice layer existed at the outer surface. The Enthalpy Method was used to simulate the phase change, and Gauss-Seidel iteration was used to solve the resulting system of finite difference equations. A set of criteria were developed for determining when a variable thickness ice layer had an effect on deicer performance. The second part was the acquisition and analysis of experimental data. The test model was a section of a Bell UH-1H helicopter blade equipped with an electrothermal deicer. A total of fifty-two thermocouples were utilized to document the thermal response of the blade and decier assembly. In the deicing runs, the experimental temperature response data clearly showed when melting, shedding or refreezing occurred. The tests illustrated that the criterion for shedding in the three cases where it did occur was that the abrasion shield interface temperature was 32 to 34 F. The third part concerned the validation of a one dimensional transient thermal model of an electrothermal deicer by comparison of the predictions with the experimental data. The Enthalpy Method was found to effectively model the phase change which occurred, and the ice shedding algorithm employed in the simulation was also evaluated.

Leffel, K. L.↗
NASA NTRS1986

Application of program LAURA to three-dimensional AOTV flowfields

Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is a robust, finite volume, single-level storage, implicit upwind differencing algorithm which has been documented and tested on several three-dimensional blunt-body flows. The algorithm can run at unlimited Courant numbers (relaxing the steady-state equations) but requires the inversion of only a 5 x 5 matrix per computational cell. An alternating directional sweep Gauss-Seidel substitution strategy is used to relax the governing equations. At present, the Euler and thin-layer Navier-Stokes equations using Sutherland's law for viscosity have been modeled for a perfect gas, equilibrium air, and nonequilibrium air chemistry neglecting diffusion. The equilibrium and nonequilibrium air chemistry options have been described in a companion paper. Good comparisons with experimental data and another calculation method for pressure distributions, aerodynamic coefficients, and heat-transfer distributions have been demonstrated for three-dimensional blunt-body flows.

Gnoffo, P. A.↗
NASA NTRS1986

Comparison of Euler and full potential marching techniques for flows over complex configurations

Two recently developed aerodynamic prediction techniques based on the steady full potential equation and the unsteady Euler equations have been applied to a variety of three-dimensional supersonic flow problems exhibiting embedded subsonic regions. Both techniques utilize planar Gauss-Seidel relaxation in the marching direction and approximate factorization in the cross-flow plane. A conservative switching scheme and flux bias technique are employed in the full potential method to transition from the supersonic marching procedure to a subsonic relaxation algorithm and vice versa. A new unified approach with finite volume, high accuracy (up to third order) Total Variation Diminishing formulation (based on Roe's scheme) is used in the Euler solver. In the supersonic regions of the flow an 'infinitely large' time step is employed, and a finite time step is applied in the subsonic regions of the flow to reach the steady-state as a time-asymptote. Numerical solutions are obtained for a number of complex configurations, including: (1) an elliptic waverider, (2) a realistic fighter configuration, (3) the Space Shuttle, and (4) a Shuttle-like configuration. Both the Full Potential and Euler numerical results are in good agreement with available experimental data.

Szema, K. Y.↗
NASA NTRS1986

Implicit, vectorizable schemes for the flux-difference split, three-dimensional Navier-Stokes equations

Two hybrid upwind models are defined for solving the Euler equations. The algorithms both employ approximate factorization (AF) in crossplane and symmetric block Gauss-Seidel relaxation in the third direction. One approach adds an additional factorization step to lower the number of required grid point operations for inversion of the block tridiagonal matrices; however, the move permits only one third of the operations to be vectorized. Finite difference solutions are calculated on a C-H-type grid, in this case enveloping a slender, sharp-edged delta wing. Sample data are provided for the calculated vortex flow for Re of 10,000, at a 20.5 deg angle of attack, represented in a crossflow velocity vector plot and in a spanwise pressure coefficient distribution. The AF scheme, without additional factorization, when used with a grid covering 51 x 51 x 72 points provides a convergent solution with no time step lasting longer than 0.00001 sec.

Liu, C. H.↗
NASA NTRS1986

An implicit flux-difference splitting scheme for three-dimensional, incompressible Navier-Stokes solutions to leading edge vortex flows

A new, implicit finite-difference scheme designed to solve the conservative, flux-difference split Navier-Stokes equations is used to compute incompressible vortex flows around delta wings. The completely vectorizable hybrid algorithm is constructed in delta form for steady state solutions independent of the time-step sizes. The scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. The governing equations are solved in curvilinear, body-fitted coordinates for treating complex geometries. The computed flow field results are compared with other theoretical and experimental data.

Hartwich, P.-M.↗
NASA NTRS1986

Implicit hybrid schemes for the flux-difference split, three-dimensional Navier-Stokes equations

Implicit hybrid algorithms employing symmetric planar Gauss-Seidel (SPGS) relaxation and either block-tridiagonally structured coefficient matrices (AF-SPGS) or block-triangular coefficient matrices (LU-SPGS) are derived to solve the flux-difference-split Navier-Stokes equations for three-dimensional incompressible flow in an upwind scheme. The physical basis of the approach is discussed, and results for problems involving vortex flow around a thin delta wing at Reynolds numbers 900,000 and 10,000 are presented graphically. It is found that AF-SPGS converges faster on vector computers which depend on long vector lengths to achieve optimum performance, whereas LU-SPGS is preferable on sequentially operating machines and vector computers using shorter vector lengths.

Hartwich, P.-M.↗
NASA NTRS1987

Efficient solutions to the Euler equations for supersonic flow with embedded subsonic regions

A line Gauss-Seidel (LGS) relaxation algorithm in conjunction with a one-parameter family of upwind discretizations of the Euler equations in two dimensions is described. Convergence of the basic algorithm to the steady state is quadratic for fully supersonic flows and is linear for other flows. This is in contrast to the block alternating direction implicit methods (either central or upwind differenced) and the upwind biased relaxation schemes, all of which converge linearly, independent of the flow regime. Moreover, the algorithm presented herein is easily coupled with methods to detect regions of subsonic flow embedded in supersonic flow. This allows marching by lines in the supersonic regions, converging each line quadratically, and iterating in the subsonic regions, and yields a very efficient iteration strategy. Numerical results are presented for two-dimensional supersonic and transonic flows containing oblique and normal shock waves which confirm the efficiency of the iteration strategy.

Walters, Robert W.↗
NASA NTRS1987

Incompressible Navier-Stokes solutions for a sharp-edged double-delta wing

An implicit finite-difference scheme is used to compute the incompressible laminar vortical flow around a double-delta wing with an aspect ratio of 2.06. By adding a time derivative of the pressure to the continuity equation, the unsteady incompressible Navier-Stokes equations can be integrated like a conventional parabolic time-dependent system of equations. The flux-difference split scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. The trajectory of the vortical core is well predicted in comparison with experimental data.

Hsu, Chung-Hao↗
NASA NTRS1987

Three dimensional hypersonic flow simulations with the CSCM implicit upwind Navier-Stokes method

The Bardina and Lombard (1985) three-dimensional CSCM Navier-Stokes method is presently extended to the simulation of complex hypersonic reentry vehicle external flows at angle of attack. The robust stability of the method derives from the combination of conservative implicit upwind flux difference splitting with a three-dimensional diagonally-dominant approximate factorization and relaxation scheme and characteristic-based implicit boundary approximations. The method's efficiency derives from an implicit symmetric Gauss-Seidel 'method of planes' relaxation scheme with alternating directional space marching sweeps along the flow coordinate direction.

Bardina, Jorge↗
NASA NTRS1987

Incompressible Navier-Stokes computations of vortical flows over double-delta wings

An implicit flux-difference splitting scheme is used to compute incompressible laminar vortical flows about a thin round-edged double-delta wing with 80 deg and 60 deg leading-edged sweep for the strake and the main wing respectively. The numerical scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. It is second-order accurate spatially by applying TVD-like upwind discretization to the inviscid fluxes and central differencing to the viscous shear fluxes. The interaction between the two primary vortices emanating from the apex and kink leading-edges is successfully simulated. The computed trajectories of vortical cores compare well with experimental data.

Hsu, Chung-Hao↗
NASA NTRS1987

Three-dimensional incompressible Navier-Stokes simulations of slender-wing vortices

The application of an implicit flux-difference splitting scheme to leading-edge vortex flows is described. The numerical scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. Second-order spatial accuracy is achieved by applying TVD-like upwind discretization to the inviscid fluxes and central differencing to the viscous shear fluxes. Good agreement between computed results and experimental data has been obtained for three different low-aspect-ratio wings.

Hsu, Chung-Hao↗