Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tensors”

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 1,117 records · Page 62

[Two-phase Interfaces in Weak External Fields]

Our aim has been that of understanding from first principles the behavior of two-phase interfaces in the absence of gravitational constraints. This is fundamental to our ability to deal with the fluid structures that abound in the real biological, chemical, and physical world. A substantial effort was mounted to determine how familiar hydrodynamic concepts have to be modified and interpreted to make them appropriate to the multi-level structure alluded to above. This was primarily in the context of the microscopic symmetric pressure tensor, which was, for the first time, expressed in the invaluable density functional format, and the used to follow the predictions of popular microscopic models of the energetics of interfacial systems. In the course of these investigations, the previous murky relation between pressure tensor and thermodynamics was completely clarified. The process of extending thermodynamic information to interfacial dynamics was initiated along two paths. One was from the viewpoint of an inertialess lattice gas, resulting in the surprising conclusion that at this level, all transport is governed by precisely the thermodynamic free energy, albeit with a non-trivial effective particle mobility. The other aimed at understanding the fashion in which slow macroscopic motions, accounted for by a time-varying microscopic energy, generate effective hydrodynamic parameters. By examining a solvable model system, it was found that all current procedures for doing so are deficient, and suitable alleviation suggested. The major effect of this project was to set the stage for the analysis of the substantial dynamical regimes in which extensive equilibrium information provides the dominant background. This produces a smooth junction to the models of Araki and Munakata, Giacomin and Lebowitz, and Oxtoby. It is also crucial to our understanding of the complex interfacial equilibrium configurations required for intermediate stages of two-phase separation, for which the "phase-field" techniques we have been developing are uniquely effective. And it puts us within striking range of the physical experiments that can provide incisive tests of the theoretical framework.

Percus, J. K.↗

Integrated Computational System for Aerodynamic Steering and Visualization

In February of 1994, an effort from the Fluid Dynamics and Information Sciences Divisions at NASA Ames Research Center with McDonnel Douglas Aerospace Company and Stanford University was initiated to develop, demonstrate, validate and disseminate automated software for numerical aerodynamic simulation. The goal of the initiative was to develop a tri-discipline approach encompassing CFD, Intelligent Systems, and Automated Flow Feature Recognition to improve the utility of CFD in the design cycle. This approach would then be represented through an intelligent computational system which could accept an engineer's definition of a problem and construct an optimal and reliable CFD solution. Stanford University's role focused on developing technologies that advance visualization capabilities for analysis of CFD data, extract specific flow features useful for the design process, and compare CFD data with experimental data. During the years 1995-1997, Stanford University focused on developing techniques in the area of tensor visualization and flow feature extraction. Software libraries were created enabling feature extraction and exploration of tensor fields. As a proof of concept, a prototype system called the Integrated Computational System (ICS) was developed to demonstrate CFD design cycle. The current research effort focuses on finding a quantitative comparison of general vector fields based on topological features. Since the method relies on topological information, grid matching and vector alignment is not needed in the comparison. This is often a problem with many data comparison techniques. In addition, since only topology based information is stored and compared for each field, there is a significant compression of information that enables large databases to be quickly searched. This report will (1) briefly review the technologies developed during 1995-1997 (2) describe current technologies in the area of comparison techniques, (4) describe the theory of our new method researched during the grant year (5) summarize a few of the results and finally (6) discuss work within the last 6 months that are direct extensions from the grant.

Hesselink, Lambertus↗

Co-Seismic Mass Dislocation and Its Effect on Earth's Rotation and Gravity

Mantle processes often involve large-scale mass transport, ranging from mantle convection, tectonic motions, glacial isostatic adjustment, to tides, atmospheric and oceanic loadings, volcanism and seismicity. On very short time scale of less than an hour, co-seismic event, apart from the "shaking" that is the earthquake, leaves behind permanent (step-function-like) dislocations in the crust and mantle. This redistribution of mass changes the Earth's inertia tensor (and hence Earth's rotation in both length-of-day and polar motion), and the gravity field (in terms of spherical harmonic Stokes coefficients). The question is whether these effects are large enough to be of any significance. In this paper we report updated calculation results. The calculation uses the normal mode summation scheme, applied to 15,814 major earthquakes that occurred during 1976-1998, according to source mechanism solutions given by the Harvard Central Moment Tensor catalog. Compared to the truly large ones earlier in the century, the earthquakes we study are individually all too small to have left any discernible signature in geodetic records of Earth rotation or global gravity field. However, their collective effects continue to exhibit an extremely strong statistical tendencies. For example, earthquakes conspire to decrease J(sub 2) and J(sub 22) while shortening LOD, resulting in a rounder and more compact Earth. Strong tendency is also seen in the earthquakes trying to "nudge" the Earth rotation pole towards about 140 degree E, roughly opposite to the observed polar drift direction. The geophysical significance and implications will be further studied.

Chao, Benjamin F.↗

Coseismic Excitation of the Earth's Polar Motion

Apart from the "shaking" near the epicenter that is the earthquake, a seismic event creates a permanent field of dislocation in the entire Earth. This redistribution of mass changes (slightly) the Earth's inertia tensor; and the Earth's rotation will change in accordance with the conservation of angular momentum. Similar to this seismic excitation of Earth rotation variations, the same mass redistribution causes (slight) changes in the Earth's gravitational field expressible in terms of changes in the Stokes coefficients of its harmonic expansion. In this paper, we give a historical background of the subject and discuss the related physics; we then compute the geodynamic effects caused by earthquakes based on normal-mode summation scheme. The effects are computed using the centroid moment tensor (CMT) solutions for 15,814 major earthquakes from Jan., 1977, through Feb., 1999, as provided in the Harvard CMT catalog. The computational results further strengthens these findings and conclusions: (i) the strong tendency for earthquakes to make the Earth rounder and more compact (however slightly) continues; (ii) so does the trend in the seismic "nudging" of the rotation pole toward the general direction of approx. 140 E, roughly opposite to that of the observed polar drift, but two orders of magnitude smaller in drift speed.

Chao, B. F.↗

Scale-Similar Models for Large-Eddy Simulations

Scale-similar models employ multiple filtering operations to identify the smallest resolved scales, which have been shown to be the most active in the interaction with the unresolved subgrid scales. They do not assume that the principal axes of the strain-rate tensor are aligned with those of the subgrid-scale stress (SGS) tensor, and allow the explicit calculation of the SGS energy. They can provide backscatter in a numerically stable and physically realistic manner, and predict SGS stresses in regions that are well correlated with the locations where large Reynolds stress occurs. In this paper, eddy viscosity and mixed models, which include an eddy-viscosity part as well as a scale-similar contribution, are applied to the simulation of two flows, a high Reynolds number plane channel flow, and a three-dimensional, nonequilibrium flow. The results show that simulations without models or with the Smagorinsky model are unable to predict nonequilibrium effects. Dynamic models provide an improvement of the results: the adjustment of the coefficient results in more accurate prediction of the perturbation from equilibrium. The Lagrangian-ensemble approach [Meneveau et al., J. Fluid Mech. 319, 353 (1996)] is found to be very beneficial. Models that included a scale-similar term and a dissipative one, as well as the Lagrangian ensemble averaging, gave results in the best agreement with the direct simulation and experimental data.

Sarghini, F.↗

Co-Seismic Mass Dislocation and its Effect on Earth's Rotation and Gravity

Mantle processes often involve large-scale mass transport, ranging from mantle convection, tectonic motions, glacial isostatic adjustment, to tides, atmospheric and oceanic loadings, volcanism and seismicity. On very short time scale of less than an hour, co-seismic event, apart from the shaking that is the earthquake, leaves behind permanent (step-function-like) dislocations in the crust and mantle. This redistribution of mass changes the Earth's inertia tensor (and hence Earth's rotation in both length-of-day and polar motion), and the gravity field (in terms of spherical harmonic Stokes coefficients). The question is whether these effects are large enough to be of any significance. In this paper we report updated calculation results based on Chao & Gross (1987). The calculation uses the normal mode summation scheme, applied to nearly twenty thousand major earthquakes that occurred during 1976-2002, according to source mechanism solutions given by the Harvard Central Moment Tensor catalog. Compared to the truly large ones earlier in the century, the earthquakes we study are individually all too small to have left any discernible signature in geodetic records of Earth rotation or global gravity field. However, their collective effects continue to exhibit an extremely strong statistical tendencies. For example, earthquakes conspire to decrease J2 and J22 while shortening LOD, resulting in a rounder and more compact Earth. Strong tendency is also seen in the earthquakes trying to nudge the Earth rotation pole towards approximately 140 degrees E, roughly opposite to the observed polar drift direction. The geophysical significance and implications will be further studied.

Chao, B. F.↗

Macroscopic Models of Radiative Transfer as Applied to Computation of the Radiation Field in the Solar Atmosphere

In this paper a new three dimensional half-moment model for radiative transfer is presented for a gray medium. It describes the evolution of the zeroth and first directional half moments of the radiative intensity. The closure is provided, similarly to Dubroca and Klar, by the maximum entropy concept. This work generalizes that model to three dimensions. The model presented here (the derivation being done in Ripoll and Wray, called the M(sup 1/2)(sub 1) model, is a hyperbolic system consisting of a total of eight equations in three dimensions, four equations for each direction. Each half model has the classical form of a macroscopic moment model in which the pressure tensor is constructed from the well-known Eddington tensor with a particular Eddington factor. Moreover, different source and border terms occur. The latter introduce couplings between the macroscopic and microscopic quantities and between the + and - streams, through the intensity in the plane perpendicular to the flux. The main theoretical application of the half moment model, treated in this paper, is its reduction to a full moment model, called M(sup +)(sub 1), for the particular but important case of a hot, opaque source radiating in a cold transparent (or semi-transparent) medium for very specific applications, such as stellar interiors or atmospheres, or combustion problems. The structure of the paper is as follows. In section 2, the model M(sup 1/2)(sub 1) is presented. In section 3, for the particular case of a hot, opaque source radiating into a cold medium, the half moment model is reduced to the M(sup +)(sub 1) model. In section 4, we first solve a simple and academic problem to validate the models, followed by a simplified solar atmosphere.

Ripoll, J.-F.↗

Parameter Estimation in Atmospheric Data Sets

In this study the structure tensor technique is used to estimate dynamical parameters in atmospheric data sets. The structure tensor is a common tool for estimating motion in image sequences. This technique can be extended to estimate other dynamical parameters such as diffusion constants or exponential decay rates. A general mathematical framework was developed for the direct estimation of the physical parameters that govern the underlying processes from image sequences. This estimation technique can be adapted to the specific physical problem under investigation, so it can be used in a variety of applications in trace gas, aerosol, and cloud remote sensing. As a test scenario this technique will be applied to modeled dust data. In this case vertically integrated dust concentrations were used to derive wind information. Those results can be compared to the wind vector fields which served as input to the model. Based on this analysis, a method to compute atmospheric data parameter fields will be presented. .

Wenig, Mark↗

A new continuum model for suspensions of gyrotactic micro-organisms

A new continuum model is formulated for dilute suspensions of swimming micro-organisms with asymmetric mass distributions. Account is taken of randomness in a cell's swimming direction, p, by postulating that the probability density function for p satisfies a Fokker-Planck equation analogous to that obtained for colloid suspensions in the presence of rotational Brownian motion. The deterministic torques on a cell, viscous and gravitational, are balanced by diffusion, represented by an isotropic rotary diffusivity Dr, which is unknown a priori, but presumably reflects stochastic influences on the cell's internal workings. When the Fokker-Planck equation is solved, macroscopic quantities such as the average cell velocity Vc, the particle diffusivity tensor D and the effective stress tensor sigma can be computed; Vc and D are required in the cell conservation equation, and sigma in the momentum equation. The Fokker-Planck equation contains two dimensionless parameters, lambda and epsilon; lambda is the ratio of the rotary diffusion time Dr-1 to the torque relaxation time B (balancing gravitational and viscous torques), while epsilon is a scale for the local vorticity or strain rate made dimensionless with B. In this paper we solve the Fokker-Planck equation exactly for epsilon = 0 (lambda arbitrary) and also obtain the first-order solution for small epsilon. Using experimental data on Vc and D obtained with the swimming alga, Chlamydomonas nivalis, in the absence of bulk flow, the epsilon = 0 results can be used to estimate the value of lambda for that species (lambda approximately 2.2; Dr approximately 0.13 s-1). The continuum model for small epsilon is then used to reanalyse the instability of a uniform suspension, previously investigated by Pedley, Hill & Kessler (1988). The only qualitatively different result is that there no longer seem to be circumstances in which disturbances with a non-zero vertical wavenumber are more unstable than purely horizontal disturbances. On the way, it is demonstrated that the only significant contribution to sigma, other than the basic Newtonian stress, is that derived from the stresslets associated with the cells' intrinsic swimming motions.

NASA Discipline Number 28-20↗

Co-Seismic Mass Displacement and its Effect on Earth's Rotation and Gravity

Mantle processes often involve large-scale mass transport, ranging from mantle convection, tectonic motions, glacial isostatic adjustment, to tides, atmospheric and oceanic loadings, volcanism and seismicity. On very short time scale of less than an hour, co-seismic event, apart from the "shaking" that is the earthquake, leaves behind permanent (step-function-like) displacements in the crust and mantle. This redistribution of mass changes the Earth's inertia tensor (and hence Earth's rotation in both length-of-day and polar motion), and the gravity field. The question is whether these effects are large enough to be of any significance. In this paper we report updated calculation results based on Chao & Gross. The calculation uses the normal mode summation scheme, applied to over twenty thousand major earthquakes that occurred during 1976-2002, according to source mechanism solutions given by the Harvard Centroid Moment Tensor catalog. Compared to the truly large ones earlier in the century, the earthquakes we study are individually all too small to have left any discernible signature in geodetic records of Earth rotation or global gravity field. However, their collective effects continue to exhibit an extremely strong statistical tendencies, conspiring to decrease J2 and J22 while shortening LOD, resulting in a rounder and more compact Earth. Strong tendency is also seen in the earthquakes trying to "nudge" the Earth rotation pole towards approx. 140 deg.E, roughly opposite to the observed polar drift direction. Currently, the Gravity Recovery And Climate Experiment (GRACE) is measuring the time-variable gravity to high degree and order with unprecedented accuracy. Our results show that great earthquakes such as the 1960 Chilean or 1964 Alaskan events cause gravitational field changes that are large enough to be detected by GRACE.

Chao, B. F.↗

Taking Measure of the Milky Way

We intend to use SIM to make definitive measurements of fundamental structural and dynamical parameters of the Milky Way. The important niche in dynamical parameter space afforded by SIM can be exploited to resolve, with unprecedented precision, a number of classical problems of Galactic astronomy. In addition, we have developed new tests of the Galactic mass distribution specifically designed for data with the special properties of SIM products. Our proposed suite of experiments will utilize the SIM Astrometric Grid as well as complementary observations of star clusters and other strategically-selected, distant "test particles" for a definitive characterization of the major components (bulge, disk, halo, satellite system) of the Milky Way. Specifically, our goals will be: 1) The determination of two fundamental parameters that play a central role in virtually every problem in Galactic astronomy, namely (a) the solar distance to the center of the Milky Way, R(sub 0); (b) the solar angular velocity around the Galactic: center, omega(sub 0). 2) The measurement of fundamental dynamical properties of the Milky Way, among them (a) the pattern speed of the central bar (b) the rotation field and velocity-dispersion tensor in the disk (c) the kinematics (mean rotational velocity and velocity dispersion tensor) of the halo as a function of position 3. The definition of the mass distribution of the Galaxy, which is dominated by the presence of dark matter. We intend to measure (a) the relative contribution of the disk and halo to the gravitational potential (b) the local volume and surface mass density of the disk (c) the shape, mass and extent of the dark halo of the Milky Way out to 250 kpc.

Majewski, Steven R.↗

A Thermodynamically Consistent Damage Model for Advanced Composites

A continuum damage model for the prediction of damage onset and structural collapse of structures manufactured in fiber-reinforced plastic laminates is proposed. The principal damage mechanisms occurring in the longitudinal and transverse directions of a ply are represented by a damage tensor that is fixed in space. Crack closure under load reversal effects are taken into account using damage variables established as a function of the sign of the components of the stress tensor. Damage activation functions based on the LaRC04 failure criteria are used to predict the different damage mechanisms occurring at the ply level. The constitutive damage model is implemented in a finite element code. The objectivity of the numerical model is assured by regularizing the dissipated energy at a material point using Bazant's Crack Band Model. To verify the accuracy of the approach, analyses of coupon specimens were performed, and the numerical predictions were compared with experimental data.

Maimi, Pere↗

Particle Demagnetization in Collisionless Magnetic Reconnection

The dissipation mechanism of magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. In this presentation, we present analytical theory results, as well as 2.5 and three-dimensional PIC simulations of guide field magnetic reconnection. We will show that diffusion region scale sizes in moderate and large guide field cases are determined by electron Larmor radii, and that analytical estimates of diffusion region dimensions need to include description of the heat flux tensor. The dominant electron dissipation process appears to be based on thermal electron inertia, expressed through nongyrotropic electron pressure tensors. We will argue that this process remains viable in three dimensions by means of a detailed comparison of high resolution particle-in-cell simulations.

Hesse, Michael↗

How Does Collisionless Magnetic Reconnection Work in the Presence of a Guide Magnetic Field?

The dissipation mechanism of guide field magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. On the other hand, a second set of studies emphasizes the role of wave-particle interactions in providing anomalous resistivity in the diffusion region. In this presentation, we present analytical theory results, as well as 2.5 and three-dimensional PIC simulations of guide-field magnetic reconnection. We will show that diffusion region scale sizes in moderate and large guide field cases are determined by electron Larmor radii, and that analytical estimates of diffusion region dimensions need to include description of the heat flux tensor. The dominant electron dissipation process appears to be based on thermal electron inertia, expressed through nongyrotropic electron pressure tensors. We will argue that this process remains viable in three dimensions by means of a detailed comparison of high resolution particle-in-cell simulations.

Hesse, Michael↗

MASCOT - MATLAB Stability and Control Toolbox

MASCOT software was created to provide the conceptual aircraft designer accurate predictions of air vehicle stability and control characteristics. The code takes as input mass property data in the form of an inertia tensor, aerodynamic loading data, and propulsion (i.e. thrust) loading data. Using fundamental non-linear equations of motion, MASCOT then calculates vehicle trim and static stability data for any desired flight condition. Common predefined flight conditions are included. The predefined flight conditions include six horizontal and six landing rotation conditions with varying options for engine out, crosswind and sideslip, plus three takeoff rotation conditions. Results are displayed through a unique graphical interface developed to provide stability and control information to the conceptual design engineers using a qualitative scale indicating whether the vehicle has acceptable, marginal, or unacceptable static stability characteristics. This software allows the user to prescribe the vehicle s CG location, mass, and inertia tensor so that any loading configuration between empty weight and maximum take-off weight can be analyzed. The required geometric and aerodynamic data as well as mass and inertia properties may be entered directly, passed through data files, or come from external programs such as Vehicle Sketch Pad (VSP). The current version of MASCOT has been tested with VSP used to compute the required data, which is then passed directly into the program. In VSP, the vehicle geometry is created and manipulated. The aerodynamic coefficients, stability and control derivatives, are calculated using VorLax, which is now available directly within VSP. MASCOT has been written exclusively using the technical computing language MATLAB . This innovation is able to bridge the gap between low-fidelity conceptual design and higher-fidelity stability and control analysis. This new tool enables the conceptual design engineer to include detailed static stability and trim constraints in the conceptual design loop. The unique graphical interface developed for this tool presents the stability data in a format that is understandable by the conceptual designer, yet also provides the detailed quantitative results if desired.

Kenny, Sean↗

Magnetospheric Multiscale (MMS) Mission Attitude Ground System Design

This paper presents an overview of the attitude ground system (AGS) currently under development for the Magnetospheric Multiscale (MMS) mission. The primary responsibilities for the MMS AGS are definitive attitude determination, validation of the onboard attitude filter, and computation of certain parameters needed to improve maneuver performance. For these purposes, the ground support utilities include attitude and rate estimation for validation of the onboard estimates, sensor calibration, inertia tensor calibration, accelerometer bias estimation, center of mass estimation, and production of a definitive attitude history for use by the science teams. Much of the AGS functionality already exists in utilities used at NASA's Goddard Space Flight Center with support heritage from many other missions, but new utilities are being created specifically for the MMS mission, such as for the inertia tensor, accelerometer bias, and center of mass estimation. Algorithms and test results for all the major AGS subsystems are presented here.

Sedlak, Joseph E.↗

Scale Dependence of Magnetic Helicity in the Solar Wind

We determine the magnetic helicity, along with the magnetic energy, at high latitudes using data from the Ulysses mission. The data set spans the time period from 1993 to 1996. The basic assumption of the analysis is that the solar wind is homogeneous. Because the solar wind speed is high, we follow the approach first pioneered by Matthaeus et al. by which, under the assumption of spatial homogeneity, one can use Fourier transforms of the magnetic field time series to construct one-dimensional spectra of the magnetic energy and magnetic helicity under the assumption that the Taylor frozen-in-flow hypothesis is valid. That is a well-satisfied assumption for the data used in this study. The magnetic helicity derives from the skew-symmetric terms of the three-dimensional magnetic correlation tensor, while the symmetric terms of the tensor are used to determine the magnetic energy spectrum. Our results show a sign change of magnetic helicity at wavenumber k approximately equal to 2AU(sup -1) (or frequency nu approximately equal to 2 microHz) at distances below 2.8AU and at k approximately equal to 30AU(sup -1) (or nu approximately equal to 25 microHz) at larger distances. At small scales the magnetic helicity is positive at northern heliographic latitudes and negative at southern latitudes. The positive magnetic helicity at small scales is argued to be the result of turbulent diffusion reversing the sign relative to what is seen at small scales at the solar surface. Furthermore, the magnetic helicity declines toward solar minimum in 1996. The magnetic helicity flux integrated separately over one hemisphere amounts to about 10(sup 45) Mx(sup 2) cycle(sup -1) at large scales and to a three times lower value at smaller scales.

Brandenburg, Axel↗

Simulations of Magnetic Reconnection - Kinetic Mechanisms Underlying the Fluid Description of Ions

Because of its ability to transfer the energy stored in magnetic field together with the breaking of the flux freezing constraint, magnetic reconnection is considered as one of the most important phenomena in plasma physics. When it happens in a collision less environment such as the terrestrial magnetosphere, it should a priori be modelled with in the framework of kinetic physics. The evidence of kinetic features has incidentally for a long time, been shown by researchers with the help of both numerical simulations and satellite observations. However, most of our understanding of the process comes from the more intuitive fluid interpretation with simple closure hypothesis which do not include kinetic effects. To what extent are these two separate descriptions of the same phenomenon related? What is the role of kinetic effects in the averaged/fluid dynamics of reconnection? This thesis addresses these questions for the proton population in the particular case of anti parallel merging with the help of 2D Hybrid simulations. We show that one can not assume, as is usually done, that the acceleration of the proton flow is only due to the Laplace force. Our results show, for symmetric and asymmetric connection, the importance of the pressure force, opposed to the electric one on the separatrices, in the decoupling region. In the symmetric case, we emphasize the kinetic origin of this force by analyzing the proton distribution functions and explain their structure by studying the underlying particle dynamics. Protons, as individual particles, are shown to bounce in the electric potential well created by the Hall effect. The spatial divergence of this well results in a mixing in phase space responsible for the observed structure of the pressure tensor. A detailed energy budget analysis confirms the role of the pressure force for the acceleration; but, contrary to what is sometimes assumed, it also reveals that the major part of the incoming Poynting flux is transferred to the thermal energy flux rather than to the convective kinetic energy flux, although the latter is generally supposed dominant. In the symmetric case, we propose the pressure tensor to be an additional proxy of the ion decoupling region in satellite data and verify this suggestion by studying a reconnection event encountered by the Cluster spacecrafts. Finally, the last part of this thesis is devoted to the study of the kinetic structure of asymmetric tangential current sheets where connection can develop. This theoretical part consists in finding a steady state solution to the Vlasov-Maxwell system for the protons in such a configuration. We present the theory and its first confrontation to numerical tests.

Aunai, icolas↗