Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tensor”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

On the Computation of Sound by Large-Eddy Simulations

The effect of the small scales on the source term in Lighthill's acoustic analogy is investigated, with the objective of determining the accuracy of large-eddy simulations when applied to studies of flow-generated sound. The distribution of the turbulent quadrupole is predicted accurately, if models that take into account the trace of the SGS stresses are used. Its spatial distribution is also correct, indicating that the low-wave-number (or frequency) part of the sound spectrum can be predicted well by LES. Filtering, however, removes the small-scale fluctuations that contribute significantly to the higher derivatives in space and time of Lighthill's stress tensor T(sub ij). The rms fluctuations of the filtered derivatives are substantially lower than those of the unfiltered quantities. The small scales, however, are not strongly correlated, and are not expected to contribute significantly to the far-field sound; separate modeling of the subgrid-scale density fluctuations might, however, be required in some configurations.

Piomelli, Ugo↗
NASA NTRS1998

The Wavelet Element Method: Realization and Additional Features in 2D and 3D - Part 2

The Wavelet Element Method (WEM) provides a construction of multiresolution systems and biorthogonal wavelets on fairly general domains. These are split into subdomains that are mapped to a single reference hypercube. Tensor products of scaling functions and wavelets defined on the unit interval are used on the reference domain. By introducing appropriate matching conditions across the interelement boundaries, a globally continuous biorthogonal wavelet basis on the general domain is obtained. This construction does not uniquely define the basis functions but rather leaves some freedom for fulfilling additional features. In this paper we detail the general construction principle of the WEM to the 1D, 2D and 3D cases. We address additional features such as symmetry, vanishing moments and minimal support of the wavelet functions in each particular dimension. The construction is illustrated by using biorthogonal spline wavelets on the interval.

Canuto, Claudio↗
NASA NTRS1998

Heat Conduction in Ceramic Coatings: Relationship Between Microstructure and Effective Thermal Conductivity

Analysis of the effective thermal conductivity of ceramic coatings and its relation to the microstructure continued. Results (obtained in Task 1) for the three-dimensional problem of heat conduction in a solid containing an inclusion (or, in particular, cavity - thermal insulator) of the ellipsoidal shape, were further advanced in the following two directions: (1) closed form expressions of H tensor have been derived for special cases of ellipsoidal cavity geometry: spheroid, crack-like spheroidal cavity and needle shaped spheroidal cavity; (2) these results for one cavity have been incorporated to construct heat energy potential for a solid with many spheroidal cavities (in the approximation of non-interacting defects). This problem constitutes a basic building block for further analyses.

Kachanov, Mark↗
NASA NTRS1998

Determination of Stress Coefficient Terms in Cracked Solids for Monoclinic Materials with Plane Symmetry at x3 = 0

Determination of all the coefficients in the crack tip field expansion for monoclinic materials under two-dimensional deformation is presented in this report. For monoclinic materials with a plane of material symmetry at x(sub 3) = 0, the in-plane deformation is decoupled from the anti-plane deformation. In the case of in-plane deformation, utilizing conservation laws of elasticity and Betti's reciprocal theorem, together with selected auxiliary fields, T-stress and third-order stress coefficients near the crack tip are evaluated first from path-independent line integrals. To determine the T-stress terms using the J-integral and Betti's reciprocal work theorem, auxiliary fields under a concentrated force and moment acting at the crack tip are used respectively. Through the use of Stroh formalism in anisotropic elasticity, analytical expressions for all the coefficients including the stress intensity factors are derived in a compact form that has surprisingly simple structure in terms of the Barnett-Lothe tensors, L. The solution forms for degenerated materials, orthotropic, and isotropic materials are presented.

Yuan, F. G.↗
NASA NTRS1998

Numerical Study of Rotating Turbulence with External Forcing

Direct numerical simulation at 256(exp 3) resolution have been carried out to study the response of isotropic turbulence to the concurrent effects of solid-body rotation and numerical forcing at the large scales. Because energy transfer to the smaller scales is weakened by rotation, energy input from forcing gradually builds up at the large scales, causing the overall kinetic energy to increase. At intermediate wavenumbers the energy spectrum undergoes a transition from a limited k(exp -5/3) inertial range to k(exp -2) scaling recently predicted in the literature. Although the Reynolds stress tensor remains approximately isotropic and three-components, evidence for anisotropy and quasi- two-dimensionality in length scales and spectra in different velocity components and directions is strong. The small scales are found to deviate from local isotropy, primarily as a result of anisotropic transfer to the high wavenumbers. To understand the spectral dynamics of this flow we study the detailed behavior of nonlinear triadic interactions in wavenumber space. Spectral transfer in the velocity component parallel to the axis of rotation is qualitatively similar to that in non-rotating turbulence; however the perpendicular component is characterized by a greatly suppressed energy cascade at high wavenumber and a local reverse transfer at the largest scales. The broader implications of this work are briefly addressed.

Yeung, P. K.↗
NASA NTRS1996

Seismic Excitation of the Polar Motion, 1977-1993

The mass redistribution in the earth as a result of an earthquake faulting changes the earth's inertia tensor, and hence its rotation. Using the complete formulae developed by CHAO and GROSS (1987) based on the normal mode theory, we calculated the earthquake-induced polar motion excitation for the largest 11,015 earthquakes that occurred during 1977.0-1993.6. The seismic excitations in this period are found to be two orders of magnitude below the detection threshold even with today's high precision earth rotation measurements. However, it was calculated that an earthquake of only one tenth the size of the great 1960 Chile event, if happened today, could be comfortably detected in polar motion observations. Furthermore, collectively these seismic excitations have a strong statistical tendency to nudge the pole towards approximately 140deg E, away from the actual observed polar drift direction. This non-random behavior, similarly found in other earthquake-induced changes in earth rotation and low-degree gravitational field by CHAO and GROSS (1987), manifests some geodynamic behavior yet to be explored.

Chao, Benjamin Fong↗
NASA NTRS1982

High-Resolution LDA Measurements of Reynolds Stress in Boundary Layers and Wakes

The turbulent character of the boundary layer and wake associated with an airfoil has been studied at a Reynolds number of 10(exp 6) and a Mach number of 0.1. To accomplish these measurements, a unique laser Doppler anemometer (LDA) has been developed that is capable of sensing two velocity components from a remote distance of 2.13 m. Using special simultaneity logic and counter-type signal processors, the geometrical features of the LDA have been exploited to provide variable spatial resolution as low as 0.2 mm. By combining the LDA with an on-line computerized data acquisition and display system, it has been possible to measure mean velocity and Reynolds stress tensor distribution at several locations along the upper surface of a 0.9 m chord, flapped airfoil installed in the Ames 7- by 10-Foot Wind Tunnel.

Orloff, Kenneth↗
NASA NTRS1999

Automated Fluid Feature Extraction from Transient Simulations

In the past, feature extraction and identification were interesting concepts, but not required to understand the underlying physics of a steady flow field. This is because the results of the more traditional tools like iso-surfaces, cuts and streamlines were more interactive and easily abstracted so they could be represented to the investigator. These tools worked and properly conveyed the collected information at the expense of much interaction. For unsteady flow-fields, the investigator does not have the luxury of spending time scanning only one "snap-shot" of the simulation. Automated assistance is required in pointing out areas of potential interest contained within the flow. This must not require a heavy compute burden (the visualization should not significantly slow down the solution procedure for co-processing environments like pV3). And methods must be developed to abstract the feature and display it in a manner that physically makes sense. The following is a list of the important physical phenomena found in transient (and steady-state) fluid flow: (1) Shocks, (2) Vortex cores, (3) Regions of recirculation, (4) Boundary layers, (5) Wakes. Three papers and an initial specification for the (The Fluid eXtraction tool kit) FX Programmer's guide were included. The papers, submitted to the AIAA Computational Fluid Dynamics Conference, are entitled : (1) Using Residence Time for the Extraction of Recirculation Regions, (2) Shock Detection from Computational Fluid Dynamics results and (3) On the Velocity Gradient Tensor and Fluid Feature Extraction.

Haimes, Robert↗
NASA NTRS1999

Towards a Rational Model for the Triple Velocity Correlations of Turbulence

This paper presents a rational approach to modelling the triple velocity correlations that appear in the transport equations for the Reynolds stresses. All existing models of these correlations have largely been formulated on phenomenological grounds and are defective in one important aspect: they all neglect to allow for the dependence of these correlations on the local gradients of mean velocity. The mathematical necessity for this dependence will be demonstrated in the paper. The present contribution lies in the novel use of Group Representation Theory to determine the most general tensorial form of these correlations in terms of all the second- and third-order tensor quantities that appear in the exact equations that govern their evolution. The requisite representation did not exist in the literature and therefore had to be developed specifically for this purpose by Professor G. F. Smith. The outcome of this work is a mathematical framework for the construction of algebraic, explicit, and rational models for the triple velocity correlations that are theoretically consistent and include all the correct dependencies. Previous models are reviewed, and all are shown to be an incomplete subset of this new representation, even to lowest order.

Younis, B. A.↗
NASA NTRS1998

Self-Consistent Simulation of the Magnetotail

Kinetic effects may be critically important at various stages of magnetotail evolution prior to and after the break up of substorms. The quasi-adiabaticity of ion dynamics results in the appearance of pronounced off diagonal terms in the pressure tensor (an effect that cannot be reproduced using the standard MHD (magnetohydrodynamic) approach). Only the appearance of these terms could support the formation of thin current sheets that recent experimental data have repeatedly indicated. We present in this paper the results of a self-consistent large-scale kinetic simulation of the Earth's magnetotail. We launch ions from the plasma mantle and self-consistently adjust the background magnetic field by computing the perturbations magnetic field from the particle current. We find that a thin current sheet with a half-thickness of approximately 0.3 R(sub E) quickly develops in the simulation. Also, the magnetotail evolves into a periodic state where the location of the X-line in the tail oscillates between x approximately 40 R(sub E) and x approximately 70 R(sub E) on a 5-minute time scale. Evidence of this periodic motion is also seen in the auroral precipitation profile of ions.

Peroomian, Vahe↗
NASA NTRS1997

Molecular Static Third-Order Polarizabilities of Carbon-Cage Fullerenes and their Correlation with Three Geometric Parameters: Group Order, Aromaticity, and Size

Calculations were performed on the valence contribution to the static molecular third-order polarizabilities (gamma) of thirty carbon-cage fullerenes (C60, C70, five isomers of C78, and twenty-three isomers of C84). The molecular structures were obtained from B3LYP/STO-3G calculations. The values of the tensor elements and an associated numerical uncertainty were obtained using the finite-field approach and polynomial expansions of orders four to eighteen of polarization versus static electric field data. The latter information was obtained from semiempirical calculations using the AM1 hamiltonian.

Moore, Craig E.↗
NASA NTRS1989

Second-Order Modeling of Low-Reynolds-Number Turbulence Near Walls

This paper presents a set of second-order closure models for low-Reynolds-number turbulence near the wall. Existing closure models for the Reynolds-stress equations were modified to show proper near-wall behavior. A dissipation-rate equation for the turbulent kinetic energy is also reformulated. The proposed models satisfy realizability and will not produce unphysical behavior. Fully developed channel flows are used for model testing. The equations are solved for the mean velocity, the Reynolds stresses, and the dissipation rate of the turbulent kinetic energy. The calculations are compared with both direct numerical simulations and with measurements. It is shown that the present models perform well in predicting the behavior of the turbulence near a wall. Significant improvements over previous models in predicting the components of the Reynolds stress tensor are obtained in the present models.

Shih, T.-H.↗
NASA NTRS1999

Geodynamic Effects of Ocean Tides: Progress and Problems

Satellite altimetry, particularly Topex/Poseidon, has markedly improved our knowledge of global tides, thereby allowing significant progress on some longstanding problems in geodynamics. This paper reviews some of that progress. Emphasis is given to global-scale problems, particularly those falling within the mandate of the new IERS Special Bureau for Tides: angular momentum, gravitational field, geocenter motion. For this discussion I use primarily the new ocean tide solutions GOT99.2, CSR4.0, and TPXO.4 (for which G. Egbert has computed inverse-theoretic error estimates), and I concentrate on new results in angular momentum and gravity and their solid-earth implications. One example is a new estimate of the effective tidal Q at the M_2 frequency, based on combining these ocean models with tidal estimates from satellite laser ranging. Three especially intractable problems are also addressed: (1) determining long-period tides in the Arctic [large unknown effect on the inertia tensor, particularly for Mf]; (2) determining the global psi_l tide [large unknown effect on interpretations of gravimetry for the near-diurnal free wobble]; and (3) determining radiational tides [large unknown temporal variations at important frequencies]. Problems (2) and (3) are related.

Richard, Ray↗
NASA NTRS1999

On the Alignment of Strain, Vorticity and Scalar Gradient in Turbulent, Buoyant, Nonpremixed Flames

The alignment of vorticity and scalar gradient with the eigendirections of the rate of strain tensor is investigated in turbulent buoyant nonpremixed horizontal and vertical flames. The uniqueness of a buoyant nonpremixed flame is that it contains regions with distinct alignment characteristics. The strain-enstrophy angle Psi is used to identify these regions. Examination of the vorticity field and the vorticity production in these different regions indicates that Psi and consequently the alignment properties near the flame surface identified by the mixture fraction band F approximately equals F(sub st) differ from those in the fuel region, F > F(sub st) and the oxidizer region, F < F(sub st). The F approximately equals F(sub st) band shows strain-dominance resulting in vorticity/alpha alignment while F > F(sub st) (and F < F(sub st) for the vertical flame) band(s) show(s) vorticity/beta alignment. The implication of this result is that the scalar dissipation, epsilon(sub F), attains its maximum value always near F approximately equals F(sub st). These results are also discussed within the framework of recent dynamical results [Galanti et al., Nonlinearity 10, 1675 (1997)] suggesting that the Navier-Stokes equations evolved towards an attracting solution. It is shown that the properties of such an attracting solution are also consistent with our results of buoyant turbulent nonpremixed flames.

Boratav, O. N.↗
NASA NTRS1998

Incompressible Turbulent Wing-Body Junction Flow

The overall objective of this study is to contribute to the optimized design of fan bypass systems in advanced turbofan engines. Increasing the engine bypass ratios have provided a major boost in engine performance improvement over the last fifty years. An engine with high bypass ratio (11-16:1) such as the Advanced Ducted Propulsion (ADP) is being developed and is expected to provide an additional 25% improvement in overall efficiency over the early turbofans. Such significant improvements in overall efficiency would reduce the cost per seat mile, which is a major government and Industry challenge for the 21th century. The research is part of the Advanced Subsonic Technology (AST) program that involves a NASA, U.S. Industry and FAA partnership with the goal of a safe and highly productive global air transportation system. The immediate objective of the study is to perform numerical simulation of duct-strut interactions to elucidate the loss mechanisms associated with this configuration that is typical of advanced turbofan engines such as ADP. However, at present experimental data for a duct-strut configuration are not available. Thus, as a first step a wing-body junction flow would be studied and is the specific objective of the present study. At the outset it is to be recognized that while duct-strut interaction flow is similar to that of wing-body junction flows, there are some differences owing to the presence of a wall at both ends of the strut. Likewise, some differences are due to the sheared inflow (as opposed to a uniform inflow) velocity profile. It is however expected that some features of a wing-body junction flow would persist. Next, some of the salient aspects of the complex flow near a wing-body junction, as revealed by various studies reported in the literature will be reviewed. One of the principle characteristics of the juncture flow, is the presence of the mean flow components in a plane perpendicular to the direction of the oncoming free-stream flow. The lateral curvature of the wing/strat causes the oncoming turbulent layer to skew about am axis (x-axis) parallel to the plane (xz-plane) of the mean shear. This is the principle mechanism for the generation of secondary flow. Such skew-induced secondary flows are slow to be attenuated by Reynolds stresses. Additional contribution to the generation of secondary flow comes from anisotropies in Reynolds stresses. Upstream of the strut, the mean-vorticity is directed span wise (along the y-direction). The presence of secondary flow in the vicinity of the strut causes the vorticity to stretch around the obstacle in a horse-shoe shape, with each leg having a vorticity of the opposite sense. The blockage effect of the strut imposes a severe adverse pressure gradient on the oncoming turbulent shear layer, causing boundary layer separation ahead of the leading edge, resulting in a vortex that rolls up and flows downstream into the juncture region. The separation vortices trailing in the wake of the wing can alter the lift or drag characteristics of the surfaces downstream of the wing-body juncture. Likewise, on submarines, the wake flow behind the appendage can degrade the performance of the propeller located downstream. The complex nature of this flow is caused by the presence of all six components of Reynolds stresses. Devenport and Simpson report that in the vicinity of the horse-shoe vortex there is intense recirculation with turbulent stresses being much larger than those normally observed in turbulent flows. These features contribute to making this flow a challenge to predict numerically. Some of the past studies provide useful insights into this flow that would guide our numerical efforts. In measurements reported by Shabaka and Bradshaw, the eddy viscosity tensor is seen to be non-isotropic and has negative components in certain regions. In an effort to evaluate the closure assumptions of various turbulence models, Devenport and Simpson used their own extensive measurements in juncture flows around the nose of a wing-body junction. Measured values of mean-velocity and/or turbulence kinetic energy was used to predict the magnitude of the shear stress vector. Algebraic stress models performed the best followed by Cebeci-Smith eddy viscosity model. The flow is reported to be dominated by a pressure field produced by the wing and the velocity field generated by the horseshoe vortex that is wrapped around the junction between the wing and wall. Kubendran et al. conclude from an experimental study that the shape of leading edge of the wing as characterized by its slenderness ratio is a major factor in determining the flow fields in the juncture region. The more thinner the leading edge of the juncture, the weaker the horseshoe vortex is. Also, with a slender leading edge, the secondary flow in the juncture would be due mainly to the cross-stream gradients of Reynolds stresses rather than due to a lateral skewing of the shear layer.

Krishnamurthy, R.↗
NASA NTRS1998

Quantity, Revisited: An Object-Oriented Reusable Class

"Quantity", a prototype implementation of an object-oriented class, was developed for two reasons: to help engineers and scientists manipulate the many types of quantities encountered during routine analysis, and to create a reusable software component to for large domain-specific applications. From being used as a stand-alone application to being incorporated into an existing computational mechanics toolkit, "Quantity" appears to be a useful and powerful object. "Quantity" has been designed to maintain the full engineering meaning of values with respect to units and coordinate systems. A value is a scalar, vector, tensor, or matrix, each of which is composed of Value Components, each of which may be an integer, floating point number, fuzzy number, etc., and its associated physical unit. Operations such as coordinate transformation and arithmetic operations are handled by member functions of "Quantity". The prototype has successfully tested such characteristics as maintaining a numeric value, an associated unit, and an annotation. In this paper we further explore the design of "Quantity", with particular attention to coordinate systems.

Funston, Monica Gayle↗
NASA NTRS1995

Determination of Eros Physical Parameters for Near Earth Asteroid Rendezvous Orbit Phase Navigation

Navigation of the orbit phase of the Near Earth steroid Rendezvous (NEAR) mission will re,quire determination of certain physical parameters describing the size, shape, gravity field, attitude and inertial properties of Eros. Prior to launch, little was known about Eros except for its orbit which could be determined with high precision from ground based telescope observations. Radar bounce and light curve data provided a rough estimate of Eros shape and a fairly good estimate of the pole, prime meridian and spin rate. However, the determination of the NEAR spacecraft orbit requires a high precision model of Eros's physical parameters and the ground based data provides only marginal a priori information. Eros is the principal source of perturbations of the spacecraft's trajectory and the principal source of data for determining the orbit. The initial orbit determination strategy is therefore concerned with developing a precise model of Eros. The original plan for Eros orbital operations was to execute a series of rendezvous burns beginning on December 20,1998 and insert into a close Eros orbit in January 1999. As a result of an unplanned termination of the rendezvous burn on December 20, 1998, the NEAR spacecraft continued on its high velocity approach trajectory and passed within 3900 km of Eros on December 23, 1998. The planned rendezvous burn was delayed until January 3, 1999 which resulted in the spacecraft being placed on a trajectory that slowly returns to Eros with a subsequent delay of close Eros orbital operations until February 2001. The flyby of Eros provided a brief glimpse and allowed for a crude estimate of the pole, prime meridian and mass of Eros. More importantly for navigation, orbit determination software was executed in the landmark tracking mode to determine the spacecraft orbit and a preliminary shape and landmark data base has been obtained. The flyby also provided an opportunity to test orbit determination operational procedures that will be used in February of 2001. The initial attitude and spin rate of Eros, as well as estimates of reference landmark locations, are obtained from images of the asteroid. These initial estimates are used as a priori values for a more precise refinement of these parameters by the orbit determination software which combines optical measurements with Doppler tracking data to obtain solutions for the required parameters. As the spacecraft is maneuvered; closer to the asteroid, estimates of spacecraft state, asteroid attitude, solar pressure, landmark locations and Eros physical parameters including mass, moments of inertia and gravity harmonics are determined with increasing precision. The determination of the elements of the inertia tensor of the asteroid is critical to spacecraft orbit determination and prediction of the asteroid attitude. The moments of inertia about the principal axes are also of scientific interest since they provide some insight into the internal mass distribution. Determination of the principal axes moments of inertia will depend on observing free precession in the asteroid's attitude dynamics. Gravity harmonics are in themselves of interest to science. When compared with the asteroid shape, some insight may be obtained into Eros' internal structure. The location of the center of mass derived from the first degree harmonic coefficients give a direct indication of overall mass distribution. The second degree harmonic coefficients relate to the radial distribution of mass. Higher degree harmonics may be compared with surface features to gain additional insight into mass distribution. In this paper, estimates of Eros physical parameters obtained from the December 23,1998 flyby will be presented. This new knowledge will be applied to simplification of Eros orbital operations in February of 2001. The resulting revision to the orbit determination strategy will also be discussed.

Miller, J. K.↗
NASA NTRS2000

Turbulence Model Predictions of Strongly Curved Flow in a U-Duct

The ability of three types of turbulence models to accurately predict the effects of curvature on the flow in a U-duct is studied. An explicit algebraic stress model performs slightly better than one- or two-equation linear eddy viscosity models, although it is necessary to fully account for the variation of the production-to-dissipation-rate ratio in the algebraic stress model formulation. In their original formulations, none of these turbulence models fully captures the suppressed turbulence near the convex wall, whereas a full Reynolds stress model does. Some of the underlying assumptions used in the development of algebraic stress models are investigated and compared with the computed flowfield from the full Reynolds stress model. Through this analysis, the assumption of Reynolds stress anisotropy equilibrium used in the algebraic stress model formulation is found to be incorrect in regions of strong curvature. By the accounting for the local variation of the principal axes of the strain rate tensor, the explicit algebraic stress model correctly predicts the suppressed turbulence in the outer part of the boundary layer near the convex wall.

Rumsey, Christopher L.↗