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 937 records · Page 52

All sources NASA NTRS DOE OSTI

NASA NTRS1990

A fractal transition in the two dimensional shear layer

The dependence of product generation with the Peclet and Reynolds number in a numerically simulated, reacting, two dimensional, temporally growing mixing layer is used to compute the fractal dimension of passive scalar interfaces. A transition from a low dimension of 4/3 to a higher one of 5/3 is identified and shown to be associated to the kinematic distortion on the flow field during the first pairing interaction. It is suggested that the structures responsible for this transition are non-deterministic, non-random, inhomogeneous fractals. Only the large scales are involved. No further transition is found for Reynolds numbers up to 20,000.

Jimenez, Javier↗

NASA NTRS1992

A differential CDM model for fatigue of unidirectional metal matrix composites

A multiaxial, isothermal, continuum damage mechanics (CDM) model for fatigue of a unidirectional metal matrix composite volume element is presented. The model is phenomenological, stress based, and assumes a single scalar internal damage variable, the evolution of which is anisotropic. The development of the fatigue damage model, (i.e., evolutionary law) is based on the definition of an initially transversely isotropic fatigue limit surface, a static fracture surface, and a normalized stress amplitude function. The anisotropy of these surfaces and function, and therefore the model, is defined through physically meaningful invariants reflecting the local stress and material orientation. This transversely isotropic model is shown, when taken to it's isotropic limit, to directly simplify to a previously developed and validated isotropic fatigue continuum damage model. Results of a nondimensional parametric study illustrate (1) the flexibility of the present formulation in attempting to characterize a class of composite materials, and (2) the capability of the formulation in predicting anticipated qualitative trends in the fatigue behavior of unidirectional metal matrix composites. Also, specific material parameters representing an initial characterization of the composite system SiC/Ti 15-3 and the matrix material (Ti 15-3) are reported.

Arnold, S. M.↗

NASA NTRS1992

Nonlinearly stable compact schemes for shock calculations

The applications of high-order, compact finite difference methods in shock calculations are discussed. The main concern is to define a local mean which will serve as a reference for introducing a local nonlinear limiting to control spurious numerical oscillations while maintaining the formal accuracy of the scheme. For scalar conservation laws, the resulting schemes can be proven total-variation stable in one space dimension and maximum-norm stable in multiple space dimensions. Numerical examples are shown to verify accuracy and stability of such schemes for problems containing shocks. These ideas can also be applied to other implicit schemes such as the continuous Galerkin finite element methods.

Cockburn, Bernardo↗

NASA NTRS1992

Treating convection in sequential solvers

The treatment of the convection terms in the sequential solver, a standard procedure found in virtually all pressure based algorithms, to compute the flow problems with sharp gradients and source terms is investigated. Both scalar model problems and one-dimensional gas dynamics equations have been used to study the various issues involved. Different approaches including the use of nonlinear filtering techniques and adoption of TVD type schemes have been investigated. Special treatments of the source terms such as pressure gradients and heat release have also been devised, yielding insight and improved accuracy of the numerical procedure adopted.

Shyy, Wei↗

NASA NTRS1992

Combustion instability analysis for liquid propellant rocket engines

The multi-dimensional numerical model has been developed to analyze the nonlinear combustion instabilities in liquid-fueled engines. The present pressure-based approach can handle the implicit pressure-velocity coupling in a non-iterative way. The additional scalar conservation equations for the chemical species, the energy, and the turbulent transport quantities can be handled by the same predictor-corrector sequences. This method is time-accurate and it can be applicable to the all-speed, transient, multi-phase, and reacting flows. Special emphasis is given to the acoustic/vaporization interaction which may act as the crucial rate-controlling mechanism in the liquid-fueled rocket engines. The subcritical vaporization is modeled to account for the effects of variable thermophysical properties, non-unitary Lewis number in the gas-film, the Stefan flow effect, and the effect of transient liquid heating. The test cases include the one-dimenisonal fast transient non-reacting and reacting flows, and the multi-dimensional combustion instabilities encountered in the liquid-fueled rocket thrust chamber. The present numerical model successfully demonstrated the capability to simulate the fast transient spray-combusting flows in terms of the limiting-cycle amplitude phenomena, correspondence between combustion and acoustics, and the steep-fronted wave and flame propagation. The investigated parameters include the spray initial conditions, air-fuel mixture ratios, and the engine geometry. Stable and unstable operating conditions are found for the liquid-fueled combustors. Under certain conditions, the limiting cycle behavior of the combusting flowfields is obtained. The numerical results indicate that the spray vaporization processes play an important role in releasing thermal energy and driving the combustion instability.

Kim, Y. M.↗

NASA NTRS1991

Characterization of frequency standard instability by estimation of their covariance matrix

The popular 3-cornered hat method used for evaluating the noise contributions of individual frequency standards is revisited. This method is used in several cases, but sometimes the results are not consistent because one or more estimated clock variances turn out to be negative. Different causes of this unacceptable result have been conjectured: among them one regards the hypothesis of uncorrelated clocks, essential in this method. Since recently realistic cases of correlation between clocks, mainly due to the environmental conditions, have been observed, this paper proposes an entirely revisited version of the 3-cornered hat method which permits to evaluate the individual variances and also the possible covariances between clocks, by relaxing the hypothesis of uncorrelation. The uncertainty and the lack of contemporaneity of the measurement series are assumed to be negligible. The lack of the uncorrelation hypothesis calls for a more general mathematical model leading to an underdetermined linear system. The estimates of the (co)variances of the measurement series us well us those of the individual clocks are introduced by means of the scalar product of the related time series and arranged in the respective covariance matrices S and R. Since covariance matrix is positive definite by definition, the problem consists in estimating the unknown R, subject to the constraint of positive definiteness, from the known S. Unfortunately, this constraint is not sufficient to estimate R. Therefore a suitable optimization criterion is proposed, which assures the positive definiteness of R and, at the same time, minimizes the global correlation among clocks. Examples of frequency instability measurements processed by the "classical" 3-cornered hat method and the here-revisited method are presented showing that the solutions are identical only when the uncorrelation hypothesis doesn't violate the positive definiteness of R.

PatriziaTravella↗

NASA NTRS1990

Land surface parameterization for regional surface heat transfer in hilly grassland

Radiometrically obtained surface temperature data and wind and temperature data from the atmospheric boundary layer are analyzed to determine the regional surface-heat flux. Monin-Obukhov similarity theory is employed to analyze the wind speed and the relationship between the surface potential temperature and the potential temperature. The radiometric scalar roughness is considered and found to be related to solar elevation and canopy height, although the regional sensible heat flux derived by the similarity theory analysis agrees with measured flux values.

Sugita, Michiaki↗

NASA NTRS1991

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. I - The dynamics of time discretization and its implications for algorithm development in computational fluid dynamics

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.

Yee, H. C.↗

NASA NTRS1991

Probing the Higgs vacuum with general relativity

It is shown that the structure of the Higgs vacuum can be revealed in gravitational experiments which probe the Schwarzschild geometry to only one order in MG/r beyond that needed for the classical tests of general relativity. The possibility that deviations from the conventional geometry are at least theoretically conceivable is explored. The deviations obtained provide a diagnostic test for searching for the existence of macroscopic scalar fields and open up the possiblity for further exploring the Higgs mechanism.

Mannheim, Philip D.↗

NASA NTRS1992

Central difference TVD and TVB schemes for time dependent and steady state problems

We use central differences to solve the time dependent Euler equations. The schemes are all advanced using a Runge-Kutta formula in time. Near shocks, a second difference is added as an artificial viscosity. This reduces the scheme to a first order upwind scheme at shocks. The switch that is used guarantees that the scheme is locally total variation diminishing (TVD). For steady state problems it is usually advantageous to relax this condition. Then small oscillations do not activate the switches and the convergence to a steady state is improved. To sharpen the shocks, different coefficients are needed for different equations and so a matrix valued dissipation is introduced and compared with the scalar viscosity. The connection between this artificial viscosity and flux limiters is shown. Any flux limiter can be used as the basis of a shock detector for an artificial viscosity. We compare the use of the van Leer, van Albada, mimmod, superbee, and the 'average' flux limiters for this central difference scheme. For time dependent problems, we need to use a small enough time step so that the CFL was less than one even though the scheme was linearly stable for larger time steps. Using a total variation bounded (TVB) Runge-Kutta scheme yields minor improvements in the accuracy.

Jorgenson, P.↗

NASA NTRS1992

Mixing in straight shear layers

Planar laser-induced fluorescence measurements were performed in a liquid plane mixing layer to extract the probability density function (pdf) of the mixture fraction of a passive scalar across the layer. Three Reynolds number (Re) cases were studied, 10,000, 33,000 and 90,000, with Re based on velocity difference and visual thickness. The results show that a non-marching pdf (central hump invariant from edge to edge of the layer) exists for Re = 10,000 but that a marching type pdf characterizes the Re = 33,000 and Re = 90,000 cases. For all cases, a broad range of mixture fraction values is found at each location across the layer. Streamwise and spanwise ramps across the layer, and structure-to-structure variation were observed and are believed to be responsible for the above behavior of the composition field. Tripping the boundary layer on the high-speed side of the splitter plate for each of the above three cases resulted in increased three-dimensionality and a change in the composition field. Average and average mixed fluid compositions are reported for all cases.

Karasso, P. S.↗

NASA NTRS1992

'Reduced' magnetohydrodynamics and minimum dissipation rates

It is demonstrated that all solutions of the equations of 'reduced' magnetohydrodynamics approach a uniform-current, zero-flow state for long times, given a constant wall electric field, uniform scalar viscosity and resistivity, and uniform mass density. This state is the state of minimum energy dissipation rate for these boundary conditions. No steady-state turbulence is possible. The result contrasts sharply with results for full three-dimensional magnetohydrodynamics before the reduction occurs.

Montgomery, David↗

NASA NTRS1991

A viscoplastic theory with thermodynamic considerations

A thermodynamic foundation using the concept of internal state variables is given for a general theory of viscoplasticity for initially isotropic materials. Three fundamental, internal, state variables are admitted; they are: a tensorial back stress for kinematic effects, and scalar drag and yield strengths for isotropic effects. All three are considered to evolve phenomenologically according to competitive processes between strain hardening, deformation-induced dynamic recovery, and thermally induced static recovery. Within this phenomenological framework, a thermodynamically admissible set of evolution equations is proposed. The theory allows each of the three internal variables to be composed as a sum of independently evolving constituents. The evolution of internal state can also include terms that vary linearly with the external variable rates, whose presence affects the energy dissipation properties of a material.

Freed, A. D.↗

NASA NTRS1992

A staggered mesh finite difference scheme for the computation of compressible flows

A simple high resolution finite difference technique is presented to approximate weak solutions to hyperbolic systems of conservation laws. The method does not rely on Riemann problem solvers and is therefore easy to extend to a wide variety of problems. The overall performance (resolution and CPU requirements) is competitive, with other state-of-the-art techniques offering sharp nonoscillatory shocks and contacts. Theoretical results confirm the reliability of the approach for linear systems and nonlinear scalar equations.

Sanders, Richard↗

NASA NTRS1991

A viscoplastic theory with thermodynamic considerations

A thermodynamic foundation using the concept of internal state variables is given for a general theory of viscoplasticity for initially isotropic materials. Three, fundamental, internal, state variables are admitted; they are: a tensorial back stress for kinematic effects, and scalar drag and yield strengths for isotropic effects. All three are considered to evolve phenomenologically according to competitive processes between strain hardening, deformation induced dynamic recovery, and thermally induced static recovery. Within this phenomenological framework, a thermodynamically admissible set of evolution equations is proposed. The theory allows each of the three internal variables to be composed as a sum of independently evolving constituents. The evolution of internal state can also include terms that vary linearly with the external variable rates, whose presence affects the energy dissipation properties of a material.

Freed, A. D.↗

NASA NTRS1992

Diskoseismology: Probing accretion disks. II - G-modes, gravitational radiation reaction, and viscosity

A scalar potential is used to derive a single partial differential equation governing the oscillation of a disk. The eigenfunctions and eigenfrequencies of a variety of disk models are found to fall into two main classes which are analogous to the p-modes and g-modes in the sun. Specifically, the eigenfunctions and eigenfrequencies of isothermal disks are computed, and the way in which these results can be generalized to other disk models is indicated. The (assumed) relatively small rates of growth or damping of the modes due to various mechanisms, in particular gravitational radiation reaction and parameterized models of viscosity are also computed. It is found that for certain parameters the p-modes are unstable to gravitational radiation reaction (CFS instability), while both the p-modes and g-modes are unstable to viscosity unless highly anisotropic viscosity models are considered.

Nowak, Michael A.↗

NASA NTRS1992

A numerical study of Li-SF6 wick combustion - Forced and mixed convective burning

A numerical study is conducted to study Li-SF6 wick diffusion flames under mixed convective burning conditions at a pressure of 0.01 MPa. Both planar and cylindrical wicks are considered. The model is based on a conserved scalar approach. The objective of this study is to assess the effects of particular parameters on the burning rate and heat transfer. The flat-plate solution yields a fuel mass burning rate per unit surface area following the x exp -1/2 dependence of the classical similarity solution, where x is the streamwise distance. Cylindrical wick geometries yield enhanced burning rates over planar wicks. For the case of mixed convective burning, the burning rate results approach either the forced or natural convective burning limits as ambient streamwise velocity is changed. Critical Richardson numbers specifying these burning limits are determined for a given condition. Reducing gravity results in a lower burning rate because the influence of natural convection is diminished. Under reduced gravity of 1/1000 of the sea-level value, mixed convective burning nearly resembles forced convection.

Damaso, R. C.↗

NASA NTRS1992

A comparative study of advanced shock-capturing schemes applied to Burgers' equation

A systematic evaluation is conducted of all extant numerical schemes for nonlinear scalar transport problems, and several advanced shock-capturing schemes are used to solve the nonlinear Burgers' equation in order to characterize their ability to resolve the sharp discontinuity, expansion zone, and propagation and collision features of shocks. For discontinuous functions, the Warming-Beam scheme generates preshock wiggles, while the Lax-Wendroff scheme generates postshock ones. Such limiters as the MUSCL or the superbee are more compressive than minimod or monotonic limiters. The performance of such TVD schemes as the upwind, the symmetric, and the Roe-Sweby, resemble each other.

Yang, H. Q.↗