Search NASA⌕ Search

SEARCH · Search NASA

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 415 records · Page 23

Modeling the pressure-strain correlation of turbulence: An invariant dynamical systems approach

The modeling of the pressure-strain correlation of turbulence is examined from a basic theoretical standpoint with a view toward developing improved second-order closure models. Invariance considerations along with elementary dynamical systems theory are used in the analysis of the standard hierarchy of closure models. In these commonly used models, the pressure-strain correlation is assumed to be a linear function of the mean velocity gradients with coefficients that depend algebraically on the anisotropy tensor. It is proven that for plane homogeneous turbulent flows the equilibrium structure of this hierarchy of models is encapsulated by a relatively simple model which is only quadratically nonlinear in the anisotropy tensor. This new quadratic model - the SSG model - is shown to outperform the Launder, Reece, and Rodi model (as well as more recent models that have a considerably more complex nonlinear structure) in a variety of homogeneous turbulent flows. Some deficiencies still remain for the description of rotating turbulent shear flows that are intrinsic to this general hierarchy of models and, hence, cannot be overcome by the mere introduction of more complex nonlinearities. It is thus argued that the recent trend of adding substantially more complex nonlinear terms containing the anisotropy tensor may be of questionable value in the modeling of the pressure-strain correlation. Possible alternative approaches are discussed briefly.

Speziale, Charles G.↗

A contribution toward rational modeling of the pressure-strain-rate correlation

A novel method of obtaining an analytical expression of the 'linear part' of the pressure-strain-rate tensor in terms of the anisotropy tensor of the Reynolds stresses has been developed, where the coefficients of the seven independent tensor terms are functions of the invariants of the Reynolds-stress anisotropy. The coefficients are evaluated up to fourth order in the anisotropy of the Reynolds stresses to provide guidance for development of a turbulence model.

Lee, Moon Joo↗

Modelling the pressure-strain correlation of turbulence - An invariant dynamical systems approach

The modeling of the pressure-strain correlation of turbulence is examined from a basic theoretical standpoint with a view toward developing improved second-order closure models. Invariance considerations along with elementary dynamical systems theory are used in the analysis of the standard hierarchy of closure models. In these commonly used models, the pressure-strain correlation is assumed to be a linear function of the mean velocity gradients with coefficients that depend algebraically on the anisotropy tensor. It is proven that for plane homogeneous turbulent flows the equilibrium structure of this hierarchy of models is encapsulated by a relatively simple model which is only quadratically nonlinear in the anisotropy tensor. This new quadratic model - the SSG model - is shown to outperform the Launder, Reece, and Rodi model (as well as more recent models that have a considerably more complex nonlinear structure) in a variety of homogeneous turbulent flows. Some deficiencies still remain for the description of rotating turbulent shear flows that are intrinsic to this general hierarchy of models and, hence, cannot be overcome by the mere introduction of more complex nonlinearities. It is thus argued that the recent trend of adding substantially more complex nonlinear terms containing the anisotropy tensor may be of questionable value in the modeling of the pressure-strain correlation. Possible alternative approaches are discussed briefly.

Speziale, Charles G.↗

Topology of three-dimensional, variable density flows

This paper is concerned with the interpretation of unsteady, variable-density flow fields. The topology of the flow is determined by finding critical points and identifying the character of local solution trajectories. The time evolution of the flow is studied by following the paths of the critical points in the three-dimensional space of invariants of the local deformations tensor. The methodology can be applied to any smooth vector field and its associated gradient tensor including the vorticity and pressure gradient fields. This approach provides a framework for describing the geometry of complex flow patterns. Concisely summarizing that geometry in the space of invariants of the local gradient tensor may be a useful way of gaining insight into time-dependent processes described by large computational data bases. Applications to the descriptions of a flickering diffusion flame and a compressible wake are discussed.

Cantwell, Brian↗

A critical comparison of second order closures with direct numerical simulation of homogeneous turbulence

Recently, several second order closure models have been proposed for closing the second moment equations, in which the velocity-pressure gradient (and scalar-pressure gradient) tensor and the dissipation rate tensor are two of the most important terms. In the literature, these correlation tensors are usually decomposed into a so called rapid term and a return-to-isotropy term. Models of these terms have been used in global flow calculations together with other modeled terms. However, their individual behavior in different flows have not been fully examined because they are un-measurable in the laboratory. Recently, the development of direct numerical simulation (DNS) of turbulence has given us the opportunity to do this kind of study. With the direct numerical simulation, we may use the solution to exactly calculate the values of these correlation terms and then directly compare them with the values from their modeled formulations (models). Here, we make direct comparisons of five representative rapid models and eight return-to-isotropy models using the DNS data of forty five homogeneous flows which were done by Rogers et al. (1986) and Lee et al. (1985). The purpose of these direct comparisons is to explore the performance of these models in different flows and identify the ones which give the best performance. The modeling procedure, model constraints, and the various evaluated models are described. The detailed results of the direct comparisons are discussed, and a few concluding remarks on turbulence models are given.

Shih, Tsan-Hsing↗

Turbulence modeling

The performance of existing two-equation eddy viscosity models was examined. An effort was made to develop better models for near-wall turbulence using direct numerical simulations of plane channel and boundary layer flows. The asymptotic near-wall behavior of turbulence was used to examine the problems of current second order closure models and develop new models with the correct near-wall behavior. Rapid Distortion Theory was used to analytically study the effects of mean deformation on turbulence, obtain analytical solutions for the spectrum tensor, Reynolds stress tensor, anisotropy tensor and its invariants, which can be used in the turbulence model development. The potential of the renormalization group theory in turbulence modeling was studied, as well as compressible turbulent flows, and modeling of bypass transition.

Shih, Tsan-Hsing↗

Cosmic microwave background probes models of inflation

Inflation creates both scalar (density) and tensor (gravity wave) metric perturbations. We find that the tensor-mode contribution to the cosmic microwave background anisotropy on large-angular scales can only exceed that of the scalar mode in models where the spectrum of perturbations deviates significantly from scale invariance. If the tensor mode dominates at large-angular scales, then the value of DeltaT/T predicted on 1 deg is less than if the scalar mode dominates, and, for cold-dark-matter models, bias factors greater than 1 can be made consistent with Cosmic Background Explorer (COBE) DMR results.

Davis, Richard L.↗

Critical comparison of second-order closures with direct numerical simulations of homogeneous turbulence

Recently, several models have been proposed for closing the second-moment equations, in which the velocity-pressure gradient tensor and the dissipation rate tensor are two of the most important terms. In the literature, these correlation tensors are usually decomposed into a so-called rapid term and a return-to-isotropy term. Models of these terms have been used in global flow calculations together with other modeled terms. However, their individual behaviors in different flows have not been fully examined because they are unmeasurable in the laboratory. Recently, the development of direct numerical simulation (DNS) of turbulence has given us the possibility to do this kind of study. With direct numerical simulation, we may use the solution to calculate exactly the values of these correlation terms and then directly compare them with the values from their modeled formulations. In this paper, we make direct comparisons of five representative rapid models and eight return-to-isotropy models using the DNS data of 45 homogeneous flows, which were done by Rogers et al. (1986) and Lee and Reynolds (1985). The purpose of these direct comparisons is to explore the performance of these models in different flows and identify the ones that give the best performance. The paper also describes the modeling procedure, model constraints, and the various evaluated models. The detailed results of the direct comparisons are discussed, and a few concluding remarks on turbulence models are given.

Shih, Tsan-Hsing↗

Search for subgrid scale parameterization by projection pursuit regression

The dependence of subgrid-scale stresses on variables of the resolved field is studied using direct numerical simulations of isotropic turbulence, homogeneous shear flow, and channel flow. The projection pursuit algorithm, a promising new regression tool for high-dimensional data, is used to systematically search through a large collection of resolved variables, such as components of the strain rate, vorticity, velocity gradients at neighboring grid points, etc. For the case of isotropic turbulence, the search algorithm recovers the linear dependence on the rate of strain (which is necessary to transfer energy to subgrid scales) but is unable to determine any other more complex relationship. For shear flows, however, new systematic relations beyond eddy viscosity are found. For the homogeneous shear flow, the results suggest that products of the mean rotation rate tensor with both the fluctuating strain rate and fluctuating rotation rate tensors are important quantities in parameterizing the subgrid-scale stresses. A model incorporating these terms is proposed. When evaluated with direct numerical simulation data, this model significantly increases the correlation between the modeled and exact stresses, as compared with the Smagorinsky model. In the case of channel flow, the stresses are found to correlate with products of the fluctuating strain and rotation rate tensors. The mean rates of rotation or strain do not appear to be important in this case, and the model determined for homogeneous shear flow does not perform well when tested with channel flow data. Many questions remain about the physical mechanisms underlying these findings, about possible Reynolds number dependence, and, given the low level of correlations, about their impact on modeling. Nevertheless, demonstration of the existence of causal relations between sgs stresses and large-scale characteristics of turbulent shear flows, in addition to those necessary for energy transfer, provides important insight into the relation between scales in turbulent flows.

Meneveau, C.↗

Anisotropic eddy viscosity models

A general discussion on the structure of the eddy viscosity tensor in anisotropic flows is presented. The systematic use of tensor symmetries and flow symmetries is shown to reduce drastically the number of independent parameters needed to describe the rank 4 eddy viscosity tensor. The possibility of using Onsager symmetries for simplifying further the eddy viscosity is discussed explicitly for the axisymmetric geometry.

Carati, D.↗

Nonguiding Center Motion and Substorm Effects in the Magnetotail

Thick and thin models of the middle magnetotail were developed using a consistent orbit tracing technique. It was found that currents carried near the equator by groups of ions with anisotropic distribution functions are not well approximated by the guiding center expressions. The guiding center equations fail primarily because the calculated pressure tensor is not magnetic field aligned. The pressure tensor becomes field aligned as one moves away from the equator, but here there is a small region in which the guiding center equations remain inadequate because the two perpendicular components of the pressure tensor are unequal. The significance of nonguiding center motion to substorm processes then was examined. One mechanism that may disrupt a thin cross-tail current sheet involves field changes that cause ions to begin following chaotic orbits. The lowest-altitude chaotic region, characterized by an adiabaticity parameter kappa approx. equal to 0.8, is especially important. The average cross-tail particle drift is slow, and we were unable to generate a thin current sheet using such ions. Therefore, any process that tends to create a thin current sheet in a region with kappa approaching 0.8 may cause the cross-tail current to get so low that it becomes insufficient to support the lobes. A different limit may be important in resonant orbit regions of a thin current sheet because particles reach a maximum cross-tail drift velocity. If the number of ions per unit length decreases as the tail is stretched, this part of the plasma sheet also may become unable to carry the cross-tail current needed to support the lobes. Thin sheets are needed for both resonant and chaotic orbit mechanisms because the distribution function must be highly structured. A description of current continuity is included to show how field aligned currents can evolve during the transition from a two-dimensional (2-D) to a 3-D configuration.

Kaufmann, Richard L.↗

Finite-Strain Fractional-Order Viscoelastic (FOV) Material Models and Numerical Methods for Solving Them

Fraction-order viscoelastic (FOV) material models have been proposed and studied in 1D since the 1930's, and were extended into three dimensions in the 1970's under the assumption of infinitesimal straining. It was not until 1997 that Drozdov introduced the first finite-strain FOV constitutive equations. In our presentation, we shall continue in this tradition by extending the standard, FOV, fluid and solid, material models introduced in 1971 by Caputo and Mainardi into 3D constitutive formula applicable for finite-strain analyses. To achieve this, we generalize both the convected and co-rotational derivatives of tensor fields to fractional order. This is accomplished by defining them first as body tensor fields and then mapping them into space as objective Cartesian tensor fields. Constitutive equations are constructed using both variants for fractional rate, and their responses are contrasted in simple shear. After five years of research and development, we now possess a basic suite of numerical tools necessary to study finite-strain FOV constitutive equations and their iterative refinement into a mature collection of material models. Numerical methods still need to be developed for efficiently solving fraction al-order integrals, derivatives, and differential equations in a finite element setting where such constitutive formulae would need to be solved at each Gauss point in each element of a finite model, which can number into the millions in today's analysis.

Freed, Alan D.↗

Geometrical approach to neural net control of movements and posture

In one approach to modeling brain function, sensorimotor integration is described as geometrical mapping among coordinates of non-orthogonal frames that are intrinsic to the system; in such a case sensors represent (covariant) afferents and motor effectors represent (contravariant) motor efferents. The neuronal networks that perform such a function are viewed as general tensor transformations among different expressions and metric tensors determining the geometry of neural functional spaces. Although the non-orthogonality of a coordinate system does not impose a specific geometry on the space, this "Tensor Network Theory of brain function" allows for the possibility that the geometry is non-Euclidean. It is suggested that investigation of the non-Euclidean nature of the geometry is the key to understanding brain function and to interpreting neuronal network function. This paper outlines three contemporary applications of such a theoretical modeling approach. The first is the analysis and interpretation of multi-electrode recordings. The internal geometries of neural networks controlling external behavior of the skeletomuscle system is experimentally determinable using such multi-unit recordings. The second application of this geometrical approach to brain theory is modeling the control of posture and movement. A preliminary simulation study has been conducted with the aim of understanding the control of balance in a standing human. The model appears to unify postural control strategies that have previously been considered to be independent of each other. Third, this paper emphasizes the importance of the geometrical approach for the design and fabrication of neurocomputers that could be used in functional neuromuscular stimulation (FNS) for replacing lost motor control.

Review↗

Attitude Ground System (AGS) For The Magnetospheric Multi-Scale (MMS) Mission

The Magnetospheric Multiscale (MMS) mission is a Solar-Terrestrial Probe mission consisting of four identically instrumented spin-stabilized spacecraft flying in an adjustable pyramid-like formation around the Earth. The formation of the MMS spacecraft allows for three-dimensional study of the phenomenon of magnetic reconnection, which is the primary objective of the mission. The MMS spacecraft were launched early on March 13, 2015 GMT. Due to the challenging and very constricted attitude and orbit requirements for performing the science, as well as the need to maintain the spacecraft formation, multiple ground functionalities were designed to support the mission. These functionalities were incorporated into a ground system known as the Attitude Ground System (AGS). Various AGS configurations have been used widely to support a variety of three-axis-stabilized and spin-stabilized spacecraft missions within the NASA Goddard Space Flight Center (GSFC). The original MMS operational concept required the AGS to perform highly accurate predictions of the effects of environmental disturbances on the spacecraft orientation and to plan the attitude maneuvers necessary to stay within the science attitude tolerance. The orbit adjustment requirements for formation control drove the need also to perform calibrations that have never been done before in support of NASA GSFC missions. The MMS mission required support analysts to provide fast and accurately calibrated values of the inertia tensor, center of mass, and accelerometer bias for each MMS spacecraft. During early design of the AGS functionalities, a Kalman filter for estimating the attitude, body rates, center of mass, and accelerometer bias, using only star tracker and accelerometer measurements, was heavily analyzed. A set of six distinct filters was evaluated and considered for estimating the spacecraft attitude and body rates using star tracker data only. Four of the six filters are closely related and were compared during support of the Time History of Events and Macroscale Interactions during Substorms (THEMIS) and Space Technology-5 (ST-5) missions. These analyses exposed high dependency and sensitivity on the knowledge of the spacecraft inertia tensor for both body rates and accelerometer bias estimation. The conclusion of the analysis led to the design of an inertia tensor calibration technique using only star tracker data. The second most important result of the analysis was the design of two separate Kalman filters to estimate the spacecraft attitude and body rates and the accelerometer bias instead of a single combined filter. In this paper, the calibration results of the mass properties, as well as the performance of the spacecraft attitude and body rates filters using flight data are presented and compared against the mission requirements.

Estimation↗

Multigrid Preconditioning for a Space-Time Spectral-Element Discontinuous-Galerkin Solver

In this work we examine a multigrid preconditioning approach in the context of a high- order tensor-product discontinuous-Galerkin spectral-element solver. We couple multigrid ideas together with memory lean and efficient tensor-product preconditioned matrix-free smoothers. Block ILU(0)-preconditioned GMRES smoothers are employed on the coarsest spaces. The performance is evaluated on nonlinear problems arising from unsteady scale- resolving solutions of the Navier-Stokes equations: separated low-Mach unsteady ow over an airfoil from laminar to turbulent ow. A reduction in the number of ne space iterations is observed, which proves the efficiency of the approach in terms of preconditioning the linear systems, however this gain was not reflected in the CPU time. Finally, the preconditioner is successfully applied to problems characterized by stiff source terms such as the set of RANS equations, where the simple tensor product preconditioner fails. Theoretical justification about the findings is reported and future work is outlined.

Franciolini, Matteo↗

Solar Sail Torque Model Characterization for the Near Earth Asteroid Scout Mission

Near Earth Asteroid Scout (NEA Scout) was a mission to test solar sail propulsion for orbital transfer from cislunar space to flyby and image an asteroid. Had it succeeded, one of the mission goals was to characterize the solar torque on the sail to ensure successful attitude control for the orbit transfer and imaging the asteroid. The simulation used to develop the flight attitude control software uses the generalized model for solar sails, a tensor equation of the forces and torques on sails of arbitrary shape. Rios-Reyes and Scheeres developed a general process to update the torque tensor coefficients using estimates of sail torque over a range of directions to the sun. Their process was adapted and implemented for the specific case of NEA Scout using spacecraft telemetry collected during sail characterization maneuvers in combination with simulation models and parameters. The NEA Scout maneuvers were limited to the operating range of the mission and constraints of the control hardware and allowed safe testing of each attitude before proceeding to the next. The NEA Scout reaction wheel speeds are used to measure accumulated momentum, while the Active Mass Translator (AMT) position is used to subtract out the torque from the center of mass crossed with the sail force and isolate the torque from only the sail shape. The process was tested by running attitude control simulations of the characterization maneuvers, generating simulated telemetry, estimating the solar torques, then using a least squares estimating the solar torque coefficients using least-squares and then performing a least-squares fit to the solar torque tensor coefficients. These estimated coefficients were tested by evaluating the solar torques under the same conditions as the simulated telemetry and comparing to the true simulated torques. Solar force model updates can be performed separately by observing the effect of the sail on the trajectory, and the torque model can be refined using those solar force updates. This process met the needs of the NEA Scout mission and can be adapted to characterize the solar torque for other missions with different sails.

solar sail↗

Solar Sail Torque Model Characterization for the Near Earth Asteroid Scout Mission

Near Earth Asteroid Scout (NEA Scout) was a mission to test solar sail propulsion for orbital transfer from cislunar space to flyby and image an asteroid. One of the goals of the mission was to characterize the solar torque on the sail to ensure successful attitude control for the orbit transfer and imaging the asteroid. The simulation used to develop the flight attitude control software uses the generalized model for solar sails, a tensor equation of the forces and torques on sails of arbitrary shape. Rios-Reyes and Scheeres developed a general process to update the torque tensor coefficients using estimates of sail torque over a range of directions to the sun. Their process was adapted and implemented for the specific case of NEA Scout using spacecraft telemetry collected during sail characterization maneuvers in combination with simulation models and parameters. The NEA Scout maneuvers were limited to the operating range of the mission and constraints of the control hardware and allowed safe testing of each attitude before proceeding to the next. The NEA Scout reaction wheel speeds are used to measure accumulated momentum, while the Active Mass Translator (AMT) position is used to subtract out the torque from the center of mass crossed with the sail force and isolate the torque from only the sail shape. The process was tested by running attitude control simulations of the characterization maneuvers, generating simulated telemetry, estimating the solar torques, then using a least squares estimating the solar torque coefficients using least-squares and then performing a least-squares fit to the solar torque tensor coefficients. These estimated coefficients were tested by evaluating the solar torques under the same conditions as the simulated telemetry and comparing to the true simulated torques. Solar force model updates can be performed separately by observing the effect of the sail on the trajectory, and the torque model can be refined using those solar force updates. This process met the needs of the NEA Scout mission and can be adapted to characterize the solar torque for other missions with different sails.

solar sail↗

First constraints on causal sources of primordial gravitational waves from BICEP/Keck, SPTpol, SPT-3G, Planck and WMAP $B$-mode data

Non-inflationary sources of gravitational waves in the early Universe generically predict causality-limited tensor power spectra at low frequencies. We report the first-ever constraints on such sources based on cosmic microwave background (CMB) $B$-mode polarization measurements. Using data from BICEP/Keck, SPTpol, SPT-3G, Planck, and WMAP, we constrain the amplitude of an early causal tensor (ECT) power spectrum parameterized by $r_{ect}$, the ratio of causal tensor power to total scalar power at $k~=~0.01$ Mpc$^{-1}$, and obtain a 95% CL upper limit of $r_{ect}<$ 0.0077. Since $r_{ect}$ can easily be related to the parameters of a given theory, our bound robustly constrains a broad class of well-motivated gravitational wave sources in the early universe, including first-order cosmological phase transitions, enhanced small-scale density perturbations, and various topological defects. Finally, we translate our limit into a bound on the present-day energy density in gravitational waves at ultra-low frequencies otherwise inaccessible to traditional gravitational wave detection strategies, including pulsar timing arrays, interferometers, and resonant cavities.

Zebrowski, Jessica A. [Fermilab; Chicago U., KICP]↗