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 325 records · Page 18

Chemistry-split techniques for viscous reactive blunt body flow computations

The weak-coupling structure between the fluid and species equations has been exploited and resulted in three, closely related, time-iterative implicit techniques. While the primitive variables are solved in two separated groups and each by an Alternating Direction Implicit (ADI) factorization scheme, the rate-species Jacobian can be treated in either full or diagonal matrix form, or simply ignored. The latter two versions render the split technique to solving for species as scalar rather than vector variables. The solution is completed at the end of each iteration after determining temperature and pressure from the flow density, energy and species concentrations. Numerical experimentation has shown that the split scalar technique, using partial rate Jacobian, yields the best overall stability and consistency. Satisfactory viscous solutions were obtained for an ellipsoidal body of axis ratio 3:1 at Mach 35 and an angle of attack of 20 degrees.

Li, C. P.↗

Convenient stability criteria for difference approximations of hyperbolic initial-boundary value problems. II

The results of Goldberg and Tadmor (1985) are extended to achieve improved stability criteria for a large class of approximations to the initial boundary value problem associated with a particular hyperbolic system in a quarter plane. In a stability analysis, it is shown that the entire approximation is stable if and only if the scalar outflow components of its principal part are stable. Thus, the global stability question is reduced to that of a scalar, homogeneous outflow problem. The stability criteria for the reduced problem, which depend both on the basic scheme and the boundary conditions, but very little on the interaction between the two, are stated and used to establish previous examples and new ones, including a host of dissipative and nondissipative examples. There is no difficulty in extending the stability criteria to two-boundary problems and initial-boundary value problems with variable coefficients.

Goldberg, Moshe↗

Analysis and performance prediction of scramjet inlets utilizing a three-dimensional Navier-Stokes code

A series of inlet analysis codes (2-D, axisymmetric, 3-D) were developed which can analyze complicated flow through complex inlet geometries in a reasonably efficient manner. The codes were verified and are being used extensively to analyze practical inlet geometries both at Langley as well as industries. Newly installed VPS 32 computer will allow more complex configurations to be analyzed. Scalar FORTRAN versions are available to increase transportability of the codes for use on other Scalar computers and on the Cray vector processing computer.

Kumar, Ajay↗

An internal variable constitutive model for the large deformation of metals at high temperatures

The advent of large deformation finite element methodologies is beginning to permit the numerical simulation of hot working processes whose design until recently has been based on prior industrial experience. Proper application of such finite element techniques requires realistic constitutive equations which more accurately model material behavior during hot working. A simple constitutive model for hot working is the single scalar internal variable model for isotropic thermal elastoplasticity proposed by Anand. The model is recalled and the specific scalar functions, for the equivalent plastic strain rate and the evolution equation for the internal variable, presented are slight modifications of those proposed by Anand. The modified functions are better able to represent high temperature material behavior. The monotonic constant true strain rate and strain rate jump compression experiments on a 2 percent silicon iron is briefly described. The model is implemented in the general purpose finite element program ABAQUS.

Brown, Stuart↗

Cyclic uniaxial and biaxial hardening of type 304 stainless steel modeled by the viscoplasticity theory based on overstress

The isotropic theory of viscoplasticity based on overstress does not use a yield surface or a loading and unloading criterion. The inelastic strain rate depends on overstress, the difference between the stress and the equilibrium stress, and is assumed to be rate dependent. Special attention is paid to the modeling of elastic regions. For the modeling of cyclic hardening, such as observed in annealed Type 304 stainless steel, and additional growth law for a scalar quantity which represents the rate independent asymptotic value of the equilibrium stress is added. It is made to increase with inelastic deformation using a new scalar measure which differentiates between nonproportional and proportional loading. The theory is applied to correlate uniaxial data under two step amplitude loading including the effect of further hardening at the high amplitude and proportional and nonproportional cyclic loadings. Results are compared with corresponding experiments.

Yao, David↗

Analysis of non-premixed turbulent reacting flows

Studies of chemical reactions occurring in turbulent flows are important in the understanding of combustion and other applications. Current numerical methods are limited in their applications due to the numerical resolution required to completely capture all length scales, but, despite the fact that realistic combustion cannot be solved completely, numerical simulations can be used to give insight into the interaction between the processes of turbulence and chemical reaction. The objective was to investigate the effects of turbulent motion on the effects of chemical reaction to gain some insight on the interaction of turbulence, molecular diffusion, and chemical reaction to support modeling efforts. A direct turbulence simulation spectral code was modified to include the effects of chemical reaction and applied to an initial value problem of chemical reaction between non-premixed species. The influence of hydrodynamics on the instantaneous structure of the reaction was investigated. The local scalar dissipation rates and the local reaction rates were examined to determine the influence of vorticity or rate of strain on the reaction and the structure of the scalar field.

Leonard, Andy D.↗

Hypercluster - Parallel processing for computational mechanics

An account is given of the development status, performance capabilities and implications for further development of NASA-Lewis' testbed 'hypercluster' parallel computer network, in which multiple processors communicate through a shared memory. Processors have local as well as shared memory; the hypercluster is expanded in the same manner as the hypercube, with processor clusters replacing the normal single processor node. The NASA-Lewis machine has three nodes with a vector personality and one node with a scalar personality. Each of the vector nodes uses four board-level vector processors, while the scalar node uses four general-purpose microcomputer boards.

Blech, Richard A.↗

Aerothermal modeling program, phase 2. Element B: Flow interaction experiment

The design process was improved and the efficiency, life, and maintenance costs of the turbine engine hot section was enhanced. Recently, there has been much emphasis on the need for improved numerical codes for the design of efficient combustors. For the development of improved computational codes, there is a need for an experimentally obtained data base to be used at test cases for the accuracy of the computations. The purpose of Element-B is to establish a benchmark quality velocity and scalar measurements of the flow interaction of circular jets with swirling flow typical of that in the dome region of annular combustor. In addition to the detailed experimental effort, extensive computations of the swirling flows are to be compared with the measurements for the purpose of assessing the accuracy of current and advanced turbulence and scalar transport models.

Nikjooy, M.↗

A class of high resolution explicit and implicit shock-capturing methods

An attempt is made to give a unified and generalized formulation of a class of high resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock wave computations. Included is a systematic review of the basic design principle of the various related numerical methods. Special emphasis is on the construction of the basis nonlinear, spatially second and third order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and the flux vector splitting approaches. Generalization of these methods to efficiently include equilibrium real gases and large systems of nonequilibrium flows are discussed. Some issues concerning the applicability of these methods that were designed for homogeneous hyperbolic conservation laws to problems containing stiff source terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for 1-, 2- and 3-dimensional gas dynamics problems.

Yee, H. C.↗

A dual potential formulation of the Navier-Stokes equations

A dual potential formulation for numerically solving the Navier-Stokes equations is developed and presented. The velocity field is decomposed using a scalar and vector potential. Vorticity and dilatation are used as the dependent variables in the momentum equations. Test cases in two dimensions verify the capability to solve flows using approximations from potential flow to full Navier-Stokes simulations. A three-dimensional incompressible flow formulation is also described. An interesting feature of this approach to solving the Navier-Stokes equations is the decomposition of the velocity field into a rotational part (vector potential) and an irrotational part (scalar potential). The Helmholtz decomposition theorem allows this splitting of the velocity field. This approach has had only limited use since it increases the number of dependent variables in the solution. However, it has often been used for incompressible flows where the solution scheme is known to be fast and accurate. This research extends the usage of this method to fully compressible Navier-Stokes simulations by using the dilatation variable along with vorticity. A time-accurate, iterative algorithm is used for the uncoupled solution of the governing equations. Several levels of flow approximation are available within the framework of this method. Potential flow, Euler and full Navier-Stokes solutions are possible using the dual potential formulation. Solution efficiency can be enhanced in a straightforward way. For some flows, the vorticity and/or dilatation may be negligible in certain regions (e.g., far from a viscous boundary in an external flow). It is possible to drop the calculation of these variables then and optimize the solution speed. Also, efficient Poisson solvers are available for the potentials. The relative merits of non-primitive variables versus primitive variables for solution of the Navier-Stokes equations are also discussed.

Gegg, S. G.↗

Climatological aspects of Nimbus-7 SMMR data

An algorithm was developed for calculating simultaneously SSTs, SWs, and TAUs on a global basis using only the 6.6 and 18 GHz channels of the SMMR. Samples of the retrievals were calculated in each of eight of the SMMR years and found to produce independent results, consistent with weather charts and climatic records, even in the presence of high winds. Another new algorithm for calculating high-latitude scalar winds from Nimbus-7 Scanning Multichannel Microwave Radiometer (SMMR) data was devised and tuned with surface observations from a number of documented Arctic Polar Low events. The algorithm utilizes the horizontally and vertically polarized radiances from the 0.8 and 1.7 cm wavelength channels of the ten-channel SMMR to retrieve near surface oceanic scalar winds and cloud water in the column, and takes advantage of the relatively small fluctuations in atmospheric water vapor at high latitudes. An advantage of this algorithm for high latitude winds from SMMR over the global algorithm is an inherently better spatial resolution as a result of the shorter wavelengths used.

Gloersen, P.↗

Advances in turbulent transport modeling based on direct simulations of turbulence

Recent work using numerical simulations to study the modeling of scalar transport is outlined. A model developed by reference to simulations for homogeneous shear flow is presented. This model predicts cross-gradient transports that can give rise to the illusion of counter-gradient transport. The model is found to work remarkably well when compared to simulations for channel flow, but fails to predict important features of the scalar fluxes in mixing layers. Reasons for this failure are discussed and suggestions are made for future work.

Reynolds, W. C.↗

Investigation of heat and momentum transport in turbulent flows via numerical simulations

Turbulent transport of heat is studied by examining the flow fields obtained from a direct simulation of a turbulent channel flow. The turbulence structures associated with the velocity and scalar fields are presented using air (Pr = 0.71) as the medium. A comparison is made between the wall-layer structures identified by the temperature field and the structures found in the velocity field. Consideration is also given to the role of the organized turbulence structures in scalar transport.

Kim, John↗

On thick domain walls in general relativity

Planar scalar field configurations in general relativity differ considerably from those in flat space. It is shown that static domain walls of finite thickness in curved space-time do not possess a reflection symmetry. At infinity, the space-time tends to the Taub vacuum on one side of the wall and to the Minkowski vacuum (Rindler space-time) on the other. Massive test particles are always accelerated towards the Minkowski side, i.e., domain walls are attractive on the Taub side, but repulsive on the Minkowski side (Taub-vacuum cleaner). It is also proved that the pressure in all directions is always negative. Finally, a brief comment is made concerning the possibility of infinite, i.e., bigger than horizon size, domain walls in our universe. All of the results are independent of the form of the potential V(phi) greater than or equal to 0 of the scalar field phi.

Goetz, Guenter↗

Potential theory of radiation

A theoretical method is being developed by which the structure of a radiation field can be predicted by a radiation potential theory, similar to a classical potential theory. The introduction of a scalar potential is justified on the grounds that the spectral intensity vector is irrotational. The vector is also solenoidal in the limits of a radiation field in complete radiative equilibrium or in a vacuum. This method provides an exact, elliptic type equation that will upgrade the accuracy and the efficiency of the current CFD programs required for the prediction of radiation and flow fields. A number of interesting results emerge from the present study. First, a steady state radiation field exhibits an optically modulated inverse square law distribution character. Secondly, the unsteady radiation field is structured with two conjugate scalar potentials. Each is governed by a Klein-Gordon equation with a frictional force and a restoring force. This steady potential field structure and the propagation of radiation potentials are consistent with the well known results of classical electromagnetic theory. The extension of the radiation potential theory for spray combustion and hypersonic flow is also recommended.

Chiu, Huei-Huang↗

Thermal relics: Do we know their abundances

The relic abundance of a particle species that was once in thermal equilibrium in the expanding Universe depends upon a competition between the annihilation rate of the species and the expansion rate of the Universe. Assuming that the Universe is radiation dominated at early times the relic abundance is easy to compute and well known. At times earlier than about 1 sec after the bang there is little or no evidence that the Universe had to be radiation dominated, although that is the simplest and standard assumption. Because early-Universe relics are of such importance both to particle physics and to cosmology, three nonstandard possibilities are considered in detail for the Universe at the time a species' abundance froze in: energy density dominated by shear (i.e., anisotropic expansion), energy density dominated by some other nonrelativistic species, and energy density dominated by the kinetic energy of the scalar field that sets the gravitational constant in a Brans-Dicke-Jordan cosmological mode. In the second case the relic abundance is less than the standard value, while in the other two cases it can be enhanced by a significant factor. Two other more exotic possibilities for enhancing the relic abundance of a species are also mentioned--a larger value of Newton's constant at early times (e.g., as might occur in superstring or Kaluza-Klein theories) or a component of the energy density at early times with a very stiff equation of state (p greater than rho/3), e.g., a scalar field phi with potential V(phi) = Beta /phi/ (exp n) with n greater than 4. Results have implications for dark matter searches and searches for particle relics in general.

Kamionkowski, Marc↗

Microwave background distortions from domain walls

Domain walls arising in a cosmic phase transition after decoupling were recently proposed as seeds for the formation of large scale structure. The distortion induced in the microwave background radiation is calculated in dependence of the wall thickness, surface density, scalar field potential, cosmic redshift and the velocity of the wall. It was found that the maximal redshift distortion for both spherical and planar walls is of the order pi G sigma H(sup -1)(sub 0), where sigma is the surface energy density and H(sup -1)(sub 0) the Hubble parameter. It was also found that, for a wall thickness smaller than the horizon, walls can be treated as infinitely thin, i.e., the redshift distortion is independent of the wall thickness and the specific form of the scalar potential. For planar walls moving with a Lorentz-factor gamma the redshift distortion is enhanced by gamma cubed.

Goetz, Guenter↗

Turbulence spectrum of strongly conductive temperature field in a rapidly stirred fluid

A simulation is performed of a passive scalar field convected by a rapidly fluctuating velocity field whose correlation time approaches zero. By using a code proposed in a previous study (Chasnov et al., 1988), the turbulence spectrum of the passive temperature field in the conductive subrange is determined. A theoretical model is proposed which explains the result obtained by representing the transfer of scalar variance by an eddy conductivity, whose correlation time is limited by the correlation time of the velocity field.

Chasnov, J.↗