Search NASA⌕ Search

SEARCH · Search NASA

Results for “split accelerating”

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

Multigrid acceleration of the flux split Euler equations

Multigrid acceleration is applied to a flux-split algorithm for solving the Euler equations in two and three dimensions. The basic algorithm is an implicit spatially-split approximate factorization method. The stability of the scheme in comparison to other factorization is examined. Results are presented for two-dimensional airfoil flows and three-dimensional wing flows which demonstrate substantially improved convergence with the multigrid algorithm. An asymptotic spectral radius of 0.89 and 0.93 is attained for a 97 x 17 x 17 wing solution at subcritical and supercritical conditions, respectively.

Anderson, W. K.↗

Flux-vector splitting and Runge-Kutta methods for the Euler equations

Runge-Kutta schemes have been used as a method of solving the Euler equations exterior to an airfoil. In the past this has been coupled with central differences and an artificial vesocity in space. In this study the Runge-Kutta time-stepping scheme is coupled with an upwinded space approximation based on flux-vector splitting. Several acceleration techniques are also considered including a local time step, residual smoothing and multigrid.

Turkel, E.↗

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.↗

Time-split finite-volume method for three-dimensional blunt-body flow

An efficient numerical method for calculating plane, axisymmetric, and fully three-dimensional blunt-body flow is presented. It is a second-order-accurate, time-dependent finite-volume procedure that solves the Euler equations in integral conservation-law form. These equations are written with respect to a Cartesian coordinate system in which an embedded mesh adjusts in time to the motion of the bow shock that is automatically captured as part of the weak solution. With such an adjusting mesh, oscillations in flow properties near the shock are shown to be virtually eliminated. The scheme uses a time-splitting concept that accelerates the convergence appreciably. Comparisons are made between computed and experimental results.

Rizzi, A. W.↗

The problem of split comets in review

Cometary splitting is investigated from the dynamical and physical standpoints. A simple two-parameter model is proposed in which the rate of recession of the fragments is determined by the momentum from outgassing, so that the net differential force is of the same nature as the nongravitational perturbations detected in the motions of unsplit comets; the two parameters of the model are the differential radial acceleration and the time of splitting. It is shown that the model successfully represents the positional observations of nearly all the 21 known split comets. The following candidate triggering mechanisms for cometary splitting are considered: tidal forces, rotation, dust-mantle dumping, and radioactive heating. A model that fits the dynamical data and physical characteristics is presented which suggests that most fragments must be appreciably nonspherical and rapidly precessing.

Sekanina, Z.↗

A highly parallel multigrid-like method for the solution of the Euler equations

We consider a highly parallel multigrid-like method for the solution of the two dimensional steady Euler equations. The new method, introduced as filtering multigrid, is similar to a standard multigrid scheme in that convergence on the finest grid is accelerated by iterations on coarser grids. In the filtering method, however, additional fine grid subproblems are processed concurrently with coarse grid computations to further accelerate convergence. These additional problems are obtained by splitting the residual into a smooth and an oscillatory component. The smooth component is then used to form a coarse grid problem (similar to standard multigrid) while the oscillatory component is used for a fine grid subproblem. The primary advantage in the filtering approach is that fewer iterations are required and that most of the additional work per iteration can be performed in parallel with the standard coarse grid computations. We generalize the filtering algorithm to a version suitable for nonlinear problems. We emphasize that this generalization is conceptually straight-forward and relatively easy to implement. In particular, no explicit linearization (e.g., formation of Jacobians) needs to be performed (similar to the FAS multigrid approach). We illustrate the nonlinear version by applying it to the Euler equations, and presenting numerical results. Finally, a performance evaluation is made based on execution time models and convergence information obtained from numerical experiments.

Tuminaro, Ray S.↗

Hubble Space Telescope Servicing Mission 3A Rendezvous Operations

The Hubble Space Telescope (HST) hardware complement includes six gas bearing, pulse rebalanced rate integrating gyros, any three of which are sufficient to conduct the science mission. After the loss of three gyros between April 1997 and April 1999 due to a known corrosion mechanism, NASA decided to split the third HST servicing mission into SM3A, accelerated to October 1999, and SM3B, scheduled for November 2001. SM3A was developed as a quick turnaround 'Launch on Need' mission to replace all six gyros. Loss of a fourth gyro in November 1999 caused HST to enter Zero Gyro Sunpoint (ZGSP) safemode, which uses sun sensors and magnetometers for attitude determination and momentum bias to maintain attitude stability during orbit night. Several instances of large attitude excursions during orbit night were observed, but ZGSP performance was adequate to provide power-positive sun pointing and to support low gain antenna communications. Body rates in ZGSP were estimated to exceed the nominal 0.1 deg/sec rendezvous limit, so rendezvous operations were restructured to utilize coarse, limited life, Retrieval Mode Gyros (RMGs) under Hardware Sunpoint (HWSP) safemode. Contingency procedures were developed to conduct the rendezvous in ZGSP in the event of RMGA or HWSP computer failure. Space Shuttle Mission STS-103 launched on December 19, 1999 after a series of weather and Shuttle-related delays. After successful rendezvous and grapple under HWSP/RMGA, the crew changed out all six gyros. Following deploy and systems checkout, HST returned to full science operations.

Lee, S.↗

Sources of the ultraheavy cosmic rays

The suggestions that the source abundances of cosmic ray nuclei heavier then Fe differ significantly from Solar System abundances are not well supported by the data without assuming preferential acceleration. The Solar System abundances of Pb and Bi are split into r-, standard s-, and cyclic 8-process components; the apprarent deficiency of Pb seen in the HEAO-3 Heavy Nuclei Experiment data might indicate an absence of Pb from the recycling 8-process.

Margolis, S. H.↗

The heaviest cosmic-ray nuclei

The implications of abundance measurements of cosmic-ray nuclei heavier than Fe are examined in light of standard propagation models and nucleosynthetic theory. Suggestions that the cosmic-ray source abundances differ significantly from solar system abundances are not well supported by the data without assuming ad hoc models of preferential acceleration. The solar system abundances of Pb and Bi are split into r-, standard s-, and cyclic s-process components to assess the contributions of these nucleosynthetic processes to the cosmic-ray source. The apparent deficiency of Pb seen in the HEAO 3 HNE data might indicate an absence of Pb from the recycling s-process.

Margolis, S. H.↗

Numerical Modeling of a Vortex Stabilized Arcjet

Arcjet thrusters are being actively considered for use in Earth orbit maneuvering applications. Experimental studies are currently the chief means of determining an optimal thruster configuration. Earlier numerical studies have failed to include all of the effects found in typical arcjets including complex geometries, viscosity, and swirling flow. Arcjet geometries are large area ratio converging nozzles with centerbodies in the subsonic portion of the nozzle. The nozzle walls serve as the anode while the centerbody functions as the cathode. Viscous effects are important because the Reynolds number, based on the throat radius, is typically less than 1,000. Experimental studies have shown that a swirl or circumferential velocity component stabilizes a constricted arc. This dissertation describes the equations governing flow through a constricted arcjet thruster. An assumption that the flowfield is in local thermodynamic equilibrium leads to a single fluid plasma temperature model. An order of magnitude analysis reveals the governing fluid mechanics equations are uncoupled from the electromagnetic field equations. A numerical method is developed to solve the governing fluid mechanics equations, the Thin Layer Navier-Stokes equations. A coordinate transformation is employed in deriving the governing equations to simplify the application of boundary conditions in complex geometries. An axisymmetric formulation is employed to include the swirl velocity component as well as the axial and radial velocity components. The numerical method is an implicit finite-volume technique and allows for large time steps to reach a converged steady-state solution. The inviscid fluxes are flux-split, and Gauss-Seidel line relaxation is used to accelerate convergence. Converging-diverging nozzles with exit-to-throat area ratios up to 100:1 and annular nozzles were examined. Quantities examined included Mach number and static wall pressure distributions, and oblique shock structures. As the level of swirl and viscosity in the flowfield increased the mass flow rate and thrust decreased. The technique was used to predict the flow through a typical arcjet thruster geometry. Results indicate swirl and viscosity play an important role in the complex geometry of an arcjet.

Pawlas, Gary E.↗

Issues concerning the Orion gamma ray line observations: Line splitting and LiBeB origin

The phenomenon of broad gamma ray line splitting is discussed, together with a scenario for the suppression of the blue wings of these broad lines due to the geometry of the accelerated particle interaction region. The broad and narrow gamma ray line emissions are compared taking into account the line splitting effect. It is concluded that the observed gamma ray lines from Orion are most likely broad, implying that the low energy cosmic rays which produce this line emission consist mostly of C and heavier ions. The suppression of the proton and alpha particle abundances requires acceleration conditions such as the acceleration of the ejecta of the supernovae before mixing with the interstellar medium. Similar conditions are implied by observations of the B and Be in low metallicity stars formed during the first Gyr of galactic evolution.

Ramaty, R.↗

Geomagnetic field distortions and their effects on radiation belt particles.

One of the quadrupole terms in the main geomagnetic field is found to contribute to a north-south 'shear distortion' of the particle-drift shells, whereas one of the octupole components causes a longitude-dependent radial deformation and associated 'drift-shell splitting.' The collective action of all higher multipoles on trapped-particle motion is then used to analyze the 'true' anomalies or distortions of the internal geomagnetic field that are independent of the quadrupole-related eccentricity of the main dipole. These 'true' anomalies must originate in upper-mantle or crustal perturbations that lie relatively near the earth's surface on both sides of the mid-Atlantic ridge; they influence trapped-particle drift shells only where the latter have their closest approach to the earth (South Atlantic and South African areas). The quadrupole and octupole perturbations, on the other hand, obviously originate deep in the earth's core. In the final part of this review, we discuss the effects of external magnetospheric currents. A time-dependent symmetric ring current causes drift shells to be displaced radially, with associated particle acceleration; magnetopause currents introduce a day-night asymmetry, causing shell splitting.

Roederer, J. G.↗

An efficient solution technique for shockwave-boundary layer interactions with flow separation and slot suction effects

An efficient method for computing two-dimensional compressible Navier-Stokes flow fields is presented. The solution algorithm is a fully-implicit approximate factorization technique based on an unsymmetric line Gauss-Seidel splitting of the equation system Jacobian matrix. Convergence characteristics are improved by the addition of acceleration techniques based on Shamanskii's method for nonlinear equations and Broyden's quasi-Newton update. Characteristic-based differencing of the equations is provided by means of Van Leer's flux vector splitting. In this investigation, emphasis is placed on the fast and accurate computation of shock-wave-boundary layer interactions with and without slot suction effects. In the latter context, a set of numerical boundary conditions for simulating the transpiration flow in an open slot is devised. Both laminar and turbulent cases are considered, with turbulent closure provided by a modified Cebeci-Smith algebraic model. Comparisons with computational and experimental data sets are presented for a variety of interactions, and a fully-coupled simulation of a plenum chamber/inlet flowfield with shock interaction and suction is also shown and discussed.

Edwards, Jack R.↗

Multigrid acceleration of the isenthalpic form of the compressible flow equations

A numerical method for solving the isenthalpic form of the governing equations for compressible inviscid flows was developed. The method is based on the concept of flux vector splitting in its implicit form and was tested on several demanding configurations. Time marching to steady state was accelerated by the implementation of the multigrid procedure which very effectively increased the rate of convergence. High quality steady-state results were obtained for various test cases and required only short computational times due to the relative efficiency of the basic method.

Melson, N. Duane↗

Multigrid acceleration of the isenthalpic form of the compressible flow equations

A numerical method for solving the isenthalpic form of the governing equations for compressible inviscid flows was developed. The method is based on the concept of flux vector splitting in its implicit form and was tested on several demanding configurations. Time marching to steady state was accelerated by the implementation of the multigrid procedure which very effectively increased the rate of convergence. High quality steady-state results were obtained for various test cases and required only short computational times due to the relative efficiency of the basic method.

Melson, N. Duane↗

Efficient solution methods for the Navier-Stokes equations

Implicit finite difference schemes for solving two-dimensional and three-dimensional Euler and thin layer Navier-Stokes equations are addressed. The methods are demonstrated in fully vectorized codes for a Cray type architecture. The Beam and Warming implicit approximate factorization algorithm in generalized coordinates is used. The methods are either time accurate or accelerated non-time accurate steady state schemes. Acceleration and efficiency modifications such as matrix reduction, diagonalization, and flux split schemes are presented. Two dimensional inviscid and viscous calculations (e.g., airfoils with a deflected spoiler, circulation control airfoils, and unsteady buffeting) and of three dimensional viscous elliptical bodies, exhausting boattails, and generic oblique wing computations are discussed.

Pulliam, T. H.↗

Implicit solution methods in computational fluid dynamics

A class of implicit finite difference schemes of the Beam and Warming approximate factorization type will be addressed. The development and analysis of various aspects of this class of schemes will be given along with the motivations behind many of the choices. Various acceleration and efficiency modifications such as matrix reduction, diagonalization and flux split schemes will be presented. The methods are demonstrated in fully vectorized codes for a CRAY type architecture. The emphasis will be on the Euler equations in generalized coordinates.

Pulliam, Thomas H.↗