SEARCH · Engineering Papers

Results for “SCALARS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Krylov subspace methods on supercomputers

A short survey of recent research on Krylov subspace methods with emphasis on implementation on vector and parallel computers is presented. Conjugate gradient methods have proven very useful on traditional scalar computers, and their popularity is likely to increase as three-dimensional models gain importance. A conservative approach to derive effective iterative techniques for supercomputers has been to find efficient parallel/vector implementations of the standard algorithms. The main source of difficulty in the incomplete factorization preconditionings is in the solution of the triangular systems at each step. A few approaches consisting of implementing efficient forward and backward triangular solutions are described in detail. Polynomial preconditioning as an alternative to standard incomplete factorization techniques is also discussed. Another efficient approach is to reorder the equations so as to improve the structure of the matrix to achieve better parallelism or vectorization. An overview of these and other ideas and their effectiveness or potential for different types of architectures is given.

Saad, Youcef↗

NASA NTRS1989

Higher order dispersion in the propagation of a gravity wave packet

To the first order of approximation, the complex amplitude of a wave packet in an anisotropic and dispersive medium is convected with the group of velocity. However, a gravity wave is a vector wave. Its wave packet must be formed by superposition of various wave numbers with corresponding frequencies, as is the case for scalar waves, and additionally by superposing many eigenmodes which also depend on the wave number. To represent the vector wave packet self-consistently, it is found that a gradient term must be included in the expansion. For a Guassian wave packet, this gradient term is shown to have important implications on the velocity vector as represented by its hodograph. Numerical results show that the hodograph is influenced by the location of the relative position of interest from the center of a Gaussian pulse. Higher order expansion shows that an initial Gaussian wave packet will retain its Gaussian shape as it propagates, but the pulse will spread in all directions with its major axis undergoing a rotation. Numerical results indicate that these higher order dispersive effects may be marginally observable in the atmosphere.

Yeh, K. C.↗

NASA NTRS1988

Effects of bearing offset and flexibility on the mesh force distribution of spiral bevel gears

For straight or spiral bevel gears of nominal spherical involute design, the resultant total force vector transmitted by the gear mesh, in the absence of friction, lies in the plane of tooth contact. This force vector can be characterized by three scalar components, two orthogonal force components lying in the plane of contact and the resultant moment taken about the nominal center of the zone of contact. Equations for these three generalized force components are derived. The equations are expressed in terms of tooth pair/gear body stiffnesses, bearing/bearing support flexibility influence coefficients, the shaft input torque, deviations of the tooth running surfaces from perfect spherical involute surfaces, and bearing centerline offsets from the positions occupied by the base cone axes of the perfect involute bevel gear counterparts to the actual gears under consideration. Inertial forces arising from transverse and axial vibrations of the gear bodies are assumed to be negligible in comparison with the bearing support reaction forces.

Mark, W. D.↗

NASA NTRS1988

A new method of imposing boundary conditions in pseudospectral approximations of hyperbolic equations

A new method to impose boundary conditions for pseudospectral approximations to hyperbolic equations is suggested. This method involves the collocation of the equation at the boundary nodes as well as satisfying boundary conditions. Stability and convergence results are proven for the Chebyshev approximation of linear scalar hyperbolic equations. The eigenvalues of this method applied to parabolic equations are shown to be real and negative.

Funaro, D.↗

NASA NTRS1989

A multigrid and upwind viscous flow solver on 3-D embedded and overlapped grids

A numerically efficient method is presented for solving the three-dimensional governing equations of the viscous compressible flow about complex configurations with topologically different components. The physical domain is decomposed into regions for which the grid generation is relatively simple and virtually with no significant restrictions. The Navier-Stokes equations are solved by an implicit, approximately factored, upwind, finite-volume scheme. The block inversions and the diagonalized scalar inversions of the coefficient matrices are modified to allow the holes created in the computational domain by the embedded and overlapped grids. The convergence is accelerated by a multigrid algorithm despite the existence of such holes. The solution for s supersonic flow past a blunt-nose-cylinder at high angle-of-attack is obtained using a C-O grid embedded in a global Cartesian grid.

Baysal, Oktay↗

NASA NTRS1988

Eddy viscosity measurements in a rectangular jet

The flow field of a rectangular jet with a 2:1 aspect ratio was studied at a Reynolds number of 100,000 (Mach number 0.09) using three-dimensional laser Doppler velocimetry (LDV). Velocity gradients, Reynolds stress tensor components, and scalar eddy viscosities are presented for the major and minor axis planes of the jet. The eddy viscosity model was found to be applicable only in the direction of maximum mean velocity gradient.

Swan, David H.↗

NASA NTRS1988

Second order accurate finite difference approximations for the transonic small disturbance equation and the full potential equation

New shock-capturing finite difference approximations for solving two scalar conservation law nonlinear partial differential equations describing inviscid, isentropic, compressible flows of aerodynamics at transonic speeds are presented. A global linear stability theorem is applied to these schemes in order to derive a necessary and sufficient condition for the finite element method. A technique is proposed to render the described approximations total variation-stable by applying the flux limiters to the nonlinear terms of the difference equation dimension by dimension. An entropy theorem applying to the approximations is proved, and an implicit, forward Euler-type time discretization of the approximation is presented. Results of some numerical experiments using the approximations are reported.

Mostrel, M. M.↗

NASA NTRS1989

Field structure of collapsing wave packets in 3D strong Langmuir turbulence

A simple model is constructed for the electric fields in the collapsing wave packets found in 3D simulations of driven and damped isotropic strong Langmuir turbulence. This model, based on a spherical-harmonic decomposition of the electrostatic potential, accounts for the distribution of wave-packet shapes observed in the simulations, particularly the predominance of oblate wave packets. In contrast with predictions for undamped and undriven subsonic collapse of scalar fields, oblate vector-field wave packets do not flatten during collapse but, instead, remain approximately self-similar and rigid.

Newman, D. L.↗

NASA NTRS1989

The stability analysis of magnetohydrodynamic equilibria - Comparing the thermodynamic approach with the energy principle

This paper is a contribution to the stability analysis of current-carrying plasmas, i.e., plasma systems that are forced by external mchanisms to carry a nonrelaxing electrical current. Under restriction to translationally invariant configurations, the thermodynamic stability criterion for a multicomponent plasma is rederived within the framework of nonideal MHD. The chosen dynamics neglects scalar resistivity, but allows for other types of dissipation effects both in Ohm's law and in the equation of motion. In the second section of the paper, the thermodynamic stability criterion is compared with the ideal MHD based energy principle of Bernstein et al. With the help of Schwarz's inequality, it is shown that the former criterion is always more 'pessimistic' than the latter, i.e., that thermodynamic stability implies stability according to the MHD principle, but not vice versa. This reuslt confirms the physical plausible idea that dissipational effects tend to weaken the stability properties of current-carrying plasma equilibria by breaking the constraints of ideal MHD and allowing for possibly destabilizing effects such as magnetic field line reconfiguration.

Brinkmann, R. P.↗

NASA NTRS1989

The diffusion approximation and transport theory for cosmic rays in relativistic flows

Equations describing the transport of cosmic rays in relativistic flows in the diffusion approximation are obtained. The analysis is based on the zeroth, first, and second differential moment equations of the relativistic Boltzmann equation with a BGK collision term. A perturbation solution of the moment equations in the diffusion approximation yields both the co-moving frame particle current and viscous stresses. The resultant cosmic-ray continuity equation contains three readily recognized energy change terms: the adiabatic energy change term; the viscous shear energy change term; and a term proportional to the scalar product of the acceleration vector of the scattering frame and the heat flux.

Webb, G. M.↗

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

A fully-coupled implicit method for thermo-chemical nonequilibrium air at sub-orbital flight speeds

A CFD technique is described in which the finite-rate chemistry in thermal and chemical nonequilibrium air is fully and implicitly coupled with the fluid motion. Developed for use in the suborbital hypersonic flight speed range, the method accounts for nonequilibrium vibrational and electronic excitation and dissociation, but not ionization. The steady-state solution to the resulting system of equations is obtained by using a lower-upper factorization and symmetric Gauss-Seidel sweeping technique through Newton iteration. Inversion of the left-hand-side matrices is replaced by scalar multiplications through the use of the diagonal dominance algorithm. The code, named CENS2H (Compressible-Euler-Navier-Stokes Two-Dimensional Hypersonic), is fully vectorized and requires about 8.8 x 10 to the -5th sec per node point per iteration using a Cray X-MP computer. Converged solutions are obtained after about 2400 iterations. Sample calculations are made for a circular cylinder and a 10 percent airfoil at 5 deg angle of attack. The calculated cylinder flow field agrees with that obtained experimentally. The code predicts a 10 percent change in lift, drag, and pitching moment for the airfoil due to the thermochemical phenomena.

Park, Chul↗

NASA NTRS1989

The numerical solution of the Navier-Stokes equations for 3-dimensional, unsteady, incompressible flows by compact schemes

The present numerical method for the solution of unsteady, incompressible three-dimensional flow Navier-Stokes equations using velocity-vorticity variables and irregular Cartesian grids proceeds by solving: (1) equations of Cauchy-Riemann type for the velocity; and (2) transport-diffusion equations for the vorticity, whose solenoidal vorticity components are generated by a Poisson equation for an appropriate scalar potential. Iterations are used to solve the finite difference equations, facilitating the use of vector and parallel-computing methods; numerical experiments have verified the method's second-order spatial and temporal accuracy.

Gatski, T. B.↗

NASA NTRS1988

An invariant asymptotic formula for solutions of second-order linear ODE's

An invariant-matrix technique for the approximate solution of second-order ordinary differential equations (ODEs) of form y-double-prime = phi(x)y is developed analytically and demonstrated. A set of linear transformations for the companion matrix differential system is proposed; the diagonalization procedure employed in the final stage of the asymptotic decomposition is explained; and a scalar formulation of solutions for the ODEs is obtained. Several typical ODEs are analyzed, and it is shown that the Liouville-Green or WKB approximation is a special case of the present formula, which provides an approximation which is valid for the entire interval (0, infinity).

Gingold, H.↗

NASA NTRS1989

Light scattering by randomly oriented crystals

The scattering phase function and the degree of linear polarization for small crystals oriented randomly in space have been computed using the geometric ray tracing theory and assuming that the crystals are homogeneous and isotropic. Calculations have been carried out for the main crystal geometries. Detection of halos from crystals other than hexagonal water ice is briefly discussed. The crystal size and shape parameters have also been averaged over some simple distributions in order to examine general light scattering properties of sharp-edged particles. A scalar physical optics correction has been developed for the geometric optics phase functions. Results can be applied to light scattering from regoliths and planetary rings, and possibly also to atmospheric halos. Retroreflecting crystals in the regolith would cause an opposition spike, a phenomenon observed for many bright satellites.

Muinonen, Karri↗

NASA NTRS1989

Pdf - Transport equations for chemically reacting flows

The closure problem for the transport equations for pdf and the characteristic functions of turbulent, chemically reacting flows is addressed. The properties of the linear and closed equations for the characteristic functional for Eulerian and Lagrangian variables are established, and the closure problem for the finite-dimensional case is discussed for pdf and characteristic functions. It is shown that the closure for the scalar dissipation term in the pdf equation developed by Dopazo (1979) and Kollmann et al. (1982) results in a single integral, in contrast to the pdf, where double integration is required. Some recent results using pdf methods obtained for turbulent flows with combustion, including effects of chemical nonequilibrium, are discussed.

Kollmann, W.↗

NASA NTRS1989

Adaptation of boundary elements to the dynamics of elastic solids

The method of boundary elements is adapted to the dynamics of elastic axisymmetric solids. Three scalar potentials are used, each satisfying the Helmholtz equation and each represented by its own source distribution. The sources are distributed over a surface enclosing the solid boundary. Stress or displacement conditions on the boundary are realized by third-order derivatives of the potentials. Discrete ring elements on the source enclosure combine their influences at discrete ring elements on the solid boundary.

El-Raheb, Michael↗

NASA NTRS1989

Computational analysis of a stability robustness margin for structured real-parameter perturbations

An efficient computational method is presented for stability robustness analysis with structured real-parameter perturbations. A generic model of a class of uncertain dynamical systems is used as an example. The parameter uncertainty is characterized by a real scalar, epsilon. Multilinearity of the closed-loop characteristic polynomial is exploited to permit application of the mapping theorem to calculate the stability robustness margin. It is found that sensitive geometry of the stability boundary in the epsilon, omega-plane renders problematic the calculation of the minimum epsilon as a function of omega. This difficulty is avoided by calculating the minimum distance to the image of the uncertainty domain over omega as a function of epsilon. It is also shown that a certain class of uncertain dynamical systems has the required multilinearity property and are thus amenable to the proposed technique.

Wedell, Evan↗