SEARCH · Engineering Papers
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.
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.
The spatial structure of turbulent motion in incompressible isotropic turbulence is investigated using a conditional average in which the conditional event specifies the local deformation tensor in addition to the local velocity vector. This average gives the best estimate of the flow field around a fixed point x given the kinematic state at x. Estimates are calculated for various kinematic states in isotropic turbulence, including pure translation, pure shear, plane strain, axisymmetric strain, and pure rotation. It is demonstrated that the large-scale motion is dominated by a vortex ring structure associated with the translational component, except at critical points of the velocity vector field.
The use of the symbolic manipulation program MACSYMA for the automatic generation of constitutive equations describing the thermomechanical behavior of complex materials such as metal-matrix composites is demonstrated. Problem-specific algorithms are developed; e.g., for partial differentiation by the chain rule, automatic generation of tree structures, special tensor analysis, and the simplification of expressions. The resulting subprogram package, SDICE, is shown to be fully applicable to both isotropic and anisotropic materials; its potential usefulness as a basis for FEM analyses is indicated.
The mechanisms involved in damage and failure of laminated orthotropic composites are investigated theoretically. The continuum model developed accounts for both matrix cracks and interply delamination using second-order tensor-valued internal-state variables based on the locally averaged microcrack dynamics. The derivation of the model is given in detail, and numerical results for sample problems are presented in extensive graphs and tables. The model is shown to be effective in predicting stresses at the ply level, and significant damage-induced decreases in laminate stress states are found.
The effect of a standing wave field in a rectangular chamber on the shape and location of an acoustically positioned drop or bubble is calculated. The sample deformation and equilibrium position are obtained from an analysis of the spherical harmonic projections of the total surface stress tensor. The method of calculation relies on the assumed condition that the sample is only slightly distorted from a spherical form. The equilibrium location of a levitated drop is combined with a formula introduced by Hasegawa (1979) to calcualte the ka dependence of the radiation force function. The present theory is valid for large as well as small ka values. Calculations in the small ka limit agree with previous theories and experimental results. Examples are presented for nonplane-wave modes as well as plane-wave rectangular modes.
The stochastic estimation provides a means of estimating conditional averages from unconditional statistics, namely, the two-point correlation tensor. This technique is applied to three-dimensional data obtained in a fully turbulent boundary layer. The results of linear and nonlinear estimates are compared to the corresponding conditional averages and the agreement is found to be excellent, with little improvement for nonlinear estimates. Extensions of the technique to space-time estimates of one- and two-point conditional averages are also presented, the latter providing a means to examine structures associated with a given length scale.
A damage-dependent constitutive model for laminated composites has been developed for the combined damage modes of matrix cracks and delaminations. The model is based on the concept of continuum damage mechanics and uses second-order tensor valued internal state variables to represent each mode of damage. The internal state variables are defined as the local volume average of the relative crack face displacements. Since the local volume for delaminations is specified at the laminate level, the constitutive model takes the form of laminate analysis equations modified by the internal state variables. Model implementation is demonstrated for the laminate engineering modulus E(x) and Poisson's ratio nu(xy) of quasi-isotropic and cross-ply laminates. The model predictions are in close agreement to experimental results obtained for graphite/epoxy laminates.
The development of the Reynolds stress field was studied for flows in which an initially two-dimensional boundary layer was skewed sideways by a spanwise pressure gradient ahead of an upstream-facing wedge. Two different wedges were used, providing a variation in the boundary-layer skewing. Measurements of all components of the Reynolds stress tensor and all ten triple products were measured using a rotatable cross-wire anemometer. The results show the expected lag of the shear stress vector behind the strain rate. Comparison of the two present experiments with previous data suggests that the lag can be estimated if the radius of curvature of the free-stream streamline is known. The magnitude of the shear stress vector in the plane of the wall is seen to decrease rapidly as the boundary-layer skewing increases. The amount of decrease is apparently related to the skewing angle between the wall and the free stream. The triple products evolve rapidly and profiles in the three-dimensional boundary layer are considerably different than two-dimensional profiles, leaving little hope for gradient transport models for the Reynolds stresses. The simplified model presented by Rotta (1979) performs reasonably well providing that an appropriate value of the T-parameter is chosen.
Flow field measurements are presented of 3 subsonic rectangular cold air jets. The 3 cases presented had aspect ratios of 1 x 2, 1 x 4 at a Mach number of 0.09 and an aspect ratio of 1 x 2 at a Mach number of 0.9. All measurements were made using a 3-D laser Doppler anemoneter system. The presented data includes the mean velocity vector, all Reynolds stress tensor components, turbulent kinetic energy and velocity correlation coefficients. The data is presented in tabular and graphical form. No analysis of the measured data or comparison to other published data is made. All tabular data are available in ASCII format on MS-DOS compatible disks.
A thermodynamic foundation using the concept of internal state variables is given for a general theory of viscoplasticity, as it applies to initially isotropic materials. Three fundamental internal state variables are admitted. They are: a tensor valued back stress for kinematic effects, and the scalar valued drag and yield strengths for isotropic effects. All three are considered to phenomenologically evolve according to competitive processes between strain hardening, strain induced dynamic recovery, and time induced static recovery. Within this phenomenological framework, a thermodynamically admissible set of evolution equations is put forth. This theory allows each of the three fundamental internal variables to be composed as a sum of independently evolving constituents.
The concept of the temperature of a surface is introduced from the viewpoint of the physical chemistry of surfaces. The surface, near surface and microlayer regions of the interface are defined. Most methods measure the temperature of the microlayer or at best the near surface region and may err in representing the surface temperature. Methods based on capillary ripples actually measure the surface temperature since surface tension (or surface tension tensor when a monolayer has been spread or absorbed at the interface) is the main restoring force that controls their propagation. Light scattering methods are described for determining the elevation of very small amplitude capillary waves through the computation of various correlation functions from which the surface tension can be estimated. Procedures for estimating the surface temperature are described.
The equivalence principle prohibits the distinction of gravity from acceleration by a local measurement. However, by making a differential measurement of acceleration over a baseline, platform accelerations can be cancelled and gravity gradients detected. In an in-line superconducting gravity gradiometer, this differencing is accomplished with two spring-mass accelerometers in which the proof masses are confined to motion in a single degree of freedom and are coupled together by superconducting circuits. Platform motions appear as common mode accelerations and are cancelled by adjusting the ratio of two persistent currents in the sensing circuit. The sensing circuit is connected to a commercial SQUID amplifier to sense changes in the persistent currents generated by differential accelerations, i.e., gravity gradients. A three-axis gravity gradiometer is formed by mounting six accelerometers on the faces of a precision cube, with the accelerometers on opposite faces of the cube forming one of three in-line gradiometers. A dedicated satellite mission for mapping the earth's gravity field is an important one. Additional scientific goals are a test of the inverse square law to a part in 10(exp 10) at 100 km, and a test of the Lense-Thirring effect by detecting the relativistic gravity magnetic terms in the gravity gradient tensor for the earth.
The highlights are presented of an Engineering Memorandum, Dynamic Estimation for Floated Gradiometers. The original impetus for the work was that gradiometers, in principle, measure components of the gravity gradient tensor, plus rotation effects, similar to centrifugal and Coriolis effects in accelerometers. The problem is that the rotation effects are often quite large, compared to the gradient, and that available inertial instruments can't measure them to adequate accuracy. The paper advances the idea that, if the instruments can be floated in a package subject to very low disturbances, a dynamic estimation, based on the Euler and translational equations of motion, plus models of all the instruments, can be used to greatly strengthen the estimates of the gradient and the rotation parameters. Moreover, symmetry constraints can be imposed directly in the filter, further strengthening the solution.
Development of new material models for describing the high temperature constitutive behavior of real materials represents an important area of research in engineering disciplines. Derivation of mathematical expressions (constitutive equations) which describe this high temperature material behavior can be quite time consuming, involved and error prone; thus intelligent application of symbolic systems to facilitate this tedious process can be of significant benefit. A computerized procedure (SDICE) capable of efficiently deriving potential based constitutive models, in analytical form is presented. This package, running under MACSYMA, has the following features: partial differentiation, tensor computations, automatic grouping and labeling of common factors, expression substitution and simplification, back substitution of invariant and tensorial relations and a relational data base. Also limited aspects of invariant theory were incorporated into SDICE due to the utilization of potentials as a starting point and the desire for these potentials to be frame invariant (objective). Finally not only calculation of flow and/or evolutionary laws were accomplished but also the determination of history independent nonphysical coefficients in terms of physically measurable parameters, e.g., Young's modulus, was achieved. The uniqueness of SDICE resides in its ability to manipulate expressions in a general yet predefined order and simplify expressions so as to limit expression growth. Results are displayed when applicable utilizing index notation.
Although the amount of data yielded by the initial airborne and surface tests was modest, it was sufficient to demonstrate that the full gravity gradient tensor was successfully measured from moving platforms both in the air and on the surface. The measurements were effectively continuous with spatial along-track resolution limited only by choice of integration lengths taken to reduce noise. The airborne data were less noisy (800 E squared/Hz typical) than were the Gravity Gradiometer Survey System (GGSS) measurements taken at the surface (5000 E squared/Hz typical). Single tracks of surface gravity disturbances recovered from airborne data were accurate to 3 to 4 mgal in each component of gravity when compared to 5 x 5 mean gravity anomalies over a 90 km track. Multitrack processing yielded 2 to 3 mgal when compared to 5 x 5 mean anomalies. Deflection of the vertical recovery over a distance of 150 km was about one arcsecond.
The European ARISTOTELES mission aims at the determination of the Earth's gravity field at short wavelength with a global coverage. Gravity gradient measurements will be achieved during six months by the GRADIO instrument onboard a dedicated satellite in a near dawn-dusk sun-synchronous orbit at an altitude of 200 km. The objective is an accuracy of better than 5 mgals for gravity anomalies, at ground level for blocks of 1 x 1 deg. According to present knowledge of the potential, the recovery of higher spherical harmonics (degree and order greater than 30) is of main importance. This leads to focus on the variations of the measured components T(sub ij) of the gravity gradient tensor, at frequencies greater than 5 x 10(exp -3) Hz. The resolution, required for the gradiometer is 10(exp -2) Eotvos (i.e., 10(exp -11)/s squared) with an averaging time of 4 s.
An essentially nonoscillatory (ENO) formulation is described for hyperbolic systems of conservation laws. ENO approaches are based on smart interpolation to avoid spurious numerical oscillations. ENO schemes are a superset of Total Variation Diminishing (TVD) schemes. In the recent past, TVD formulations were used to construct shock capturing finite difference methods. At extremum points of the solution, TVD schemes automatically reduce to being first-order accurate discretizations locally, while away from extrema they can be constructed to be of higher order accuracy. The new framework helps construct essentially non-oscillatory finite difference methods without recourse to local reductions of accuracy to first order. Thus arbitrarily high orders of accuracy can be obtained. The basic general ideas of the new approach can be specialized in several ways and one specific implementation is described based on: (1) the integral form of the conservation laws; (2) reconstruction based on the primitive functions; (3) extension to multiple dimensions in a tensor product fashion; and (4) Runge-Kutta time integration. The resulting method is fourth-order accurate in time and space and is applicable to uniform Cartesian grids. The construction of such schemes for scalar equations and systems in one and two space dimensions is described along with several examples which illustrate interesting aspects of the new approach.
Flow field measurements of three subsonic rectangular cold air jets are presented. The three cases had aspect ratios of 1x2, 1x4 at a Mach number of 0.09 and an aspect ratio of 1x2 at a Mach number of 0.9. All measurements were made using a 3-D laser Doppler anemometer system. The data includes the mean velocity vector, all Reynolds stress tensor components, turbulent kinetic energy and velocity correlation coefficients. The data are presented in tabular and graphical form. No analysis of the measured data or comparison to other published data is made.