Engineering Papers⌕ Search

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.

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

At least 739 records · Page 41

Euler solutions using an implicit multigrid technique

A coarse-grid correction algorithm has been implemented into an implicit upwind Euler solver and tested for transonic airfoil problems. The Euler solver uses split-flux formulation and penta-diagonal scalar equations, respectively, for the explicit and implicit operators. The multigrid sequence starts at the fine grid level, then steps down to each coarse grid level to smooth error components using implicit operators. Estimate of residuals can be obtained by two approaches, which differ in the level where the residuals are collected. Both approaches will lead to a work reduction factor of 12 for a Mach 0.75 flow at 2 degrees incidence on a 65x26 grid. The work reduction factor is found to increase proportional to the number of grid levels.

Li, Chien-Peng↗

A high resolution shock capturing scheme for high Mach number internal flow

An accurate shock capturing scheme is developed for high speed flow. The flux limited dissipation scheme can be converted to a symmetric TVD scheme in the case of a scalar conservation law. Application to an oblique shock wave problem shows the improved resolution and accuracy when compared to the results obtained by the adaptive dissipation method, while exhibiting comparable computational efficiency.

Yoon, Seokkwan↗

An LU-SSOR scheme for the Euler and Navier-Stokes equations

A new multigrid relaxation scheme, lower-upper symmetric successive overrelaxation (LU-SSOR) is developed for the steady-state solution of the Euler and Navier-Stokes equations. The scheme, which is based on central differences, does not require flux splitting for approximate Newton iteration. Application to transonic flow shows that the new method is efficient and robust. The vectorizable LU-SSOR scheme needs only scalar diagonal inversions.

Yoon, Seokkwan↗

Effective medium approximation for effective propagation constant calculation in a dense random medium

The effective medium approximation (EMA) formalism developed for scalar wave calculations in solid state physics is generalized to electromagnetic wave scattering in a dense random medium. Results are applied to compute the effective propagation constant in a dense medium involving discrete spherical scatterers. When compared with a common quasicrystalline approximation (QCA), it is found that EMA accounts for backward scattering and the effect of correlation among three scatterers which are not available in QCA. It is also found that there is not much difference in the calculated normalized phase velocity between the use of these two approximations. However, there is a significant difference in the computed effective loss tangent in a nonabsorptive random medium. The computed effective loss tangent using EMA and measurements from a snow medium are compared, showing good agreement.

Zhu, P. Y.↗

A Constitutive Model for the Inelastic Multiaxial Cyclic Response of a Nickel Base Superalloy Rene 80

The objective was to develop unified constitutive equations which can model a variety of nonlinear material phenomena observed in Rene 80 at elevated temperatures. A constitutive model was developed based on back stress and drag stress. The tensorial back stress was used to model directional effects; whereas, the scalar drag stress was used to model isotropic effects and cyclic hardening or softening. A flow equation and evolution equations for the state variables were developed in multiaxial form. Procedures were developed to generate the material parameters. The model predicted very well the monotonic tensile, cyclic, creep, and stress relaxation behavior of Rene 80 at 982 C. The model was then extended to 871, 760, and 538 C. It was shown that strain rate dependent behavior at high temperatures and strain rate independent behavior at the lower temperatures could be predicted very well. A large number of monotonic tensile, creep, stress relation, and cyclic experiments were predicted. The multiaxial capabilities of the model were verified extensively for combined tension/torsion experiments. The prediction of the model agreed very well for proportional, nonproportional, and pure shear cyclic loading conditions at 982 and 871 C.

Ramaswamy, V. G.↗

Automating the parallel processing of fluid and structural dynamics calculations

The NASA Lewis Research Center is actively involved in the development of expert system technology to assist users in applying parallel processing to computational fluid and structural dynamic analysis. The goal of this effort is to eliminate the necessity for the physical scientist to become a computer scientist in order to effectively use the computer as a research tool. Programming and operating software utilities have previously been developed to solve systems of ordinary nonlinear differential equations on parallel scalar processors. Current efforts are aimed at extending these capabilities to systems of partial differential equations, that describe the complex behavior of fluids and structures within aerospace propulsion systems. This paper presents some important considerations in the redesign, in particular, the need for algorithms and software utilities that can automatically identify data flow patterns in the application program and partition and allocate calculations to the parallel processors. A library-oriented multiprocessing concept for integrating the hardware and software functions is described.

Arpasi, Dale J.↗

An optimal output feedback gain variation scheme for the control of plants exhibiting gross parameter changes

A concept for optimally designing output feedback controllers for plants whose dynamics exhibit gross changes over their operating regimes was developed. This was to formulate the design problem in such a way that the implemented feedback gains vary as the output of a dynamical system whose independent variable is a scalar parameterization of the plant operating point. The results of this effort include derivation of necessary conditions for optimality for the general problem formulation, and for several simplified cases. The question of existence of a solution to the design problem was also examined, and it was shown that the class of gain variation schemes developed are capable of achieving gain variation histories which are arbitrarily close to the unconstrained gain solution for each point in the plant operating range. The theory was implemented in a feedback design algorithm, which was exercised in a numerical example. The results are applicable to the design of practical high-performance feedback controllers for plants whose dynamics vary significanly during operation. Many aerospace systems fall into this category.

Moerder, Daniel D.↗

A practical assessment of spectral accuracy for hyperbolic problems with discontinuities

Physical and transform space filtering has been applied to the Fourier spectral collocation solution of the constant coefficient scalar wave equation with a discontinuous initial condition. High order accuracy can be extracted from the unfiltered solution. Smooth, high order Fourier space filtering gives expected polynomial order solutions away from the discontinuity. Spectral accuracy is observed with the physical space filter of Gottlieb and Tadmor.

Kopriva, David A.↗

Finite element models of wire rope for vibration analysis

The usefulness of wire rope in shock and vibration isolation is briefly reviewed and its modeling for the purpose of vibration analysis is addressed. A model of a nominally straight segment of wire rope is described in which the rope structure is represented by a maiden, or central, strand of wire with one (or more) strand(s) wrapped around it in a helix (helices). The individual strands are modeled using finite elements and MSC NASTRAN. Small linear segments of each wire are modeled mathematically by dividing them lengthwise into triangular prisms representing each prism by a solid NASTRAN element. To model pretensioning and allow for extraction of internal force information from the NASTRAN model, the wound strands are connected to the maiden strand and each other using spring (scalar elastic) elements. Mode shapes for a length of wire rope with one and fixed to a moving base and the other attached to a point mass, are presented. The use of the NASTRAN derived mode shapes to approximate internal normal forces in equations of motion for vibration analyses is considered.

Cochran, J. E., Jr.↗

A three-dimensional dual potential procedure with applications to wind tunnel inlets and interacting boundary layers

A dual potential decomposition of the velocity field into a scalar and a vector potential function is extended to three dimensions and used in the finite-difference simulation of steady three-dimensional inviscid rotational flows and viscous flow. The finite-difference procedure was used to simulate the flow through the 80 by 120 ft wind tunnel at NASA Ames Research Center. Rotational flow produced by the stagnation pressure drop across vanes and screens which are located at the entrance of the inlet is modeled using actuator disk theory. Results are presented for two different inlet vane and screen configurations. The numerical predictions are in good agreement with experimental data. The dual potential procedure was also applied to calculate the viscous flow along two and three dimensional troughs. Viscous effects are simulated by injecting vorticity which is computed from a boundary layer algorithm. For attached flow over a three dimensional trough, the present calculations are in good agreement with other numerical predictions. For separated flow, it is shown from a two dimensional analysis that the boundary layer approximation provides an accurate measure of the vorticity in regions close to the wall; whereas further away from the wall, caution has to be exercised in using the boundary-layer equations to supply vorticity to the dual potential formulation.

Rao, K. V.↗

A generalized procedure for constructing an upwind based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some linear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to scalar nonlinear equations and a constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance in numerical experiments. The formulation is extended to the nonlinear system. It is demonstrated that the present procedure can easily convert existing central or upwind, and second- or higher-order differencing schemes to preserve monotonicity and yield physically admissible solutions. The formulation is simple mathematically as well as numerically; both matrix-vector multiplication and Riemann solver are avoided. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made.

Liou, Meng-Sing↗

Smoothing and the second law

The technique of obtaining second order, oscillation free, total variation diminishing (TVD), scalar difference schemes by adding a limited diffusion flux (smoothing) to a second order centered scheme is explored. It is shown that such schemes do not always converge to the correct physical answer. The approach presented here is to construct schemes that numerically satisfy the second law of thermodynamics on a cell by cell basis. Such schemes can only converge to the correct physical solution and in some cases can be shown to be TVD. An explicit scheme with this property and second order spatial accuracy was found to have an extremely restrictive time step limitation (Delta t less than Delta x squared). Switching to an implicit scheme removed the time step limitation.

Merriam, Marshal L.↗

On the coalescence-dispersion modeling of turbulent molecular mixing

The general coalescence-dispersion (C/D) closure provides phenomenological modeling of turbulent molecular mixing. The models of Curl and Dopazo and O'Brien appear as two limiting C/D models that bracket the range of results one can obtain by various models. This finding is used to investigate the sensitivtiy of the results to the choice of the model. Inert scalar mixing is found to be less model-sensitive than mixing accompanied by chemical reaction. Infinitely fast chemistry approximation is used to relate the C/D approach to Toor's earlier results. Pure mixing and infinite rate chemistry calculations are compared to study further a recent result of Hsieh and O'Brien who found that higher concentration moments are not sensitive to chemistry.

Givi, Peyman↗

A new method of imposing boundary conditions for 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.↗

Addition of the AUTOSPC feature to NASTRAN

A new capability called the AUTOSPC feature has been incorporated into the April 1987 release of NASTRAN. It gives the user the option of automatically applying single-point constraints for the purpose of removing potential grid and scalar point singularities that have not been otherwise already constrained out. Details of this implementation are given. Its usage is described and it is illustrated with an example problem.

Pamidi, P. R.↗

ICL: The Image Composition Language

The Image Composition Language (ICL) provides a convenient way for programmers of interactive graphics application programs to define how the video look-up table of a raster display system is to be loaded. The ICL allows one or several images stored in the frame buffer to be combined in a variety of ways. The ICL treats these images as variables, and provides arithematic, relational, and conditional operators to combine the images, scalar variables, and constants in image composition expressions. The objective of ICL is to provide programmers with a simple way to compose images, to relieve the tedium usually associated with loading the video look-up table to obtain desired results.

Foley, James D.↗

Grid generation on surfaces in 3 dimensions

The development of a surface grid generation algorithm was initiated. The basic adaptive movement technique of mean-value-relaxation was extended from the viewpoint of a single coordinate grid over a surface described by a single scalar function to that of a surface more generally defined by vector functions and covered by a collection of smoothly connected grids. Within the multiconnected assemblage, the application of control was examined in several instances.

Eiseman, Peter R.↗

The convergence of spectral methods for nonlinear conservation laws

The convergence of the Fourier method for scalar nonlinear conservation laws which exhibit spontaneous shock discontinuities is discussed. Numerical tests indicate that the convergence may (and in fact in some cases must) fail, with or without post-processing of the numerical solution. Instead, a new kind of spectrally accurate vanishing viscosity is introduced to augment the Fourier approximation of such nonlinear conservation laws. Using compensated compactness arguments, it is shown that this spectral viscosity prevents oscillations, and convergence to the unique entropy solution follows.

Tadmor, Eitan↗