SEARCH · Search NASA

Results for “Complete equation of state”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 2: Algorithmic developments and implementation

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Li, Wei↗

NASA NTRS1990

PTA en route noise measurements

A long-range advanced turboprop en route noise database was obtained with weather, tracking, and onboard measurements. In-flight noise directivity measurements were made. Data repeatability within a test day was excellent. Day-to-day variability existed and is not completely understood and therefore not predicted. Comparison of a two-dimensional ray-tracing propagation model with the ensemble average ground-measured data was good; however, as stated above, the day-to-day data variability was not completely predicted. Future research will include looking at alternative propagation models. Three-dimensional ray tracing, fast field program, and the parabolic equation are possibilities. The effect of turbulence needs to be accessed.

Willshire, William L., Jr.↗

NASA NTRS1993

Application and development of the Schwinger multichannel scattering theory and the partial differential equation theory of electron-molecule scattering

One version of the multichannel theory for electron-target scattering based on the Schwinger variational principle, the SMC method, requires the introduction of a projection parameter. The role of the projection parameter a is investigated and it is shown that the principal-value operator in the SMC equation is Hermitian regardless of the value of a as long as it is real and nonzero. In a basis that is properly orthonormalizable, the matrix representation of this operator is also Hermitian. The use of such basis is consistent with the Schwinger variational principle because the Lippmann-Schwinger equation automatically builds in the correct boundary conditions. Otherwise, an auxiliary condition needs to be introduced, and Takatsuka and McKoy's original value of a is one of the three possible ways to achieve Hermiticity. In all cases but one, a can be uncoupled from the Hermiticity condition and becomes a free parameter. An equation for a based on the variational stability of the scattering amplitude is derived; its solution has an interesting property that the scattering amplitude from a converged SMC calculation is independent of the choice of a even though the SMC operator itself is a-dependent. This property provides a sensitive test of the convergence of the calculation. For a static-exchange calculation, the convergence requirement only depends on the completeness of the one-electron basis, but for a general multichannel case, the a-invariance in the scattering amplitude requires both the one-electron basis and the N plus 1-electron basis to be complete. The role of a in the SMC equation and the convergence property are illustrated using two examples: e-CO elastic scattering in the static-exchange approximation, and a two-state treatment of the e-H2 Chi(sup 1)Sigma(sub g)(+) yields b(sup 3)Sigma(sub u)(+) excitation.

Weatherford, Charles A.↗

NASA NTRS2018

OSIRIS-REx Encounters Bennu: Initial Assessment from the Approach Phase

The OSIRIS-REx spacecraft launched on September 8, 2016, on a seven-year journey to return samples from asteroid (101955) Bennu. This presentation summarizes the scientific results from the Approach and Preliminary Survey phases. Bennu observations are set to begin on August 17, 2018,when the asteroid is bright enough for detection by the PolyCam. PolyCam and MapCam collect data to survey the asteroid environment for any hazards and characterize the asteroid point-source photometric properties. Resolved images acquired during final approach, starting in late October 2018, allow the creation of a shape model using stereophotoclinometry (SPC), needed by both the navigation team and science planners. The OVIRS and OTES spectrometers characterize the point- source spectral properties over a full rotation period, providing a first look at any features and thermophysical properties. TAGSAM is released from the launch container and deployed into the sampling configuration then returned to the stow position.Preliminary Survey follows the Approach Phase in early December 2018. This phase consists of a series of hyperbolic trajectories that cross over the North and South poles and the equator of Bennu at a close-approach distance of 7 km. Images from these Preliminary Survey passes provide data to complete the 75-cm resolution SPC global shape model and solve for the rotation state. Once the shape model is complete, the asteroid coordinate system is defined for co-registration of all data products. These higher-resolution images also constrain the photometric properties and allow for an initial assessment of the geology. In Preliminary Survey the team also obtains the first OLA data, providing a measure of the surface topography. OVIRS and OTES collect data as "ride-along" instruments, with the spacecraft pointing driven by imaging constraints. These data provide a first look at the spectral variation across the surface of Bennu. Radio science measurements, combined with altimetry and imagery, determine Bennu's mass, a prerequisite to placing the spacecraft into orbit in late December 2018. Together, data from the Approach and Preliminary Survey phases set the stage for the extensive mapping planned for 2019. These dates are the baseline plan. Any contingency or unexpected discovery may change this mission profile.

Lauretta, Dante S.↗

NASA NTRS2006

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗

NASA NTRS1972

Effect of support flexibilty and damping on the dynamic response of a single mass flexible rotor in elastic bearings

A steady state analysis of the shaft and the bearing housing motion was made by assuming synchronous precession of the system. The conditions under which the support system would act as a dynamic vibration absorber at the rotor critical speed were studied; plots of the rotor and support amplitudes, phase angles, and forces transmitted were evaluated by the computer, and the performance curves were automatically plotted by a CalComp plotter unit. Curves are presented on the optimization of the support housing characteristics to attenuate the rotor unbalance response over the entire rotor speed range. The complete transient motion including rotor unbalance was examined by integrating the equations of motion numerically using a modified fourth order Runge-Kutta procedure, and the resulting whirl orbits were plotted by the CalComp plotter unit. The results of the transient analysis are discussed with regards to the design optimization procedure derived from the steady-state analysis.

Kirk, R. G.↗

NASA NTRS1975

Numerical computation of viscous flows on the lee side of blunt shapes flying at supersonic speeds

A numerical method for solving the parabolic approximation to the steady-state compressible Navier-Stokes equations is examined. The approximation neglects only the streamwise gradients of shear stress. An implicit finite difference method is used which advances the solution downstream from an initial data surface and determines the complete viscous-inviscid flow between the body and bow shock wave. It is necessary that the inviscid portion of the flow field be supersonic. Crossflow separation is determined as part of the solution. The method is applied to a 15 deg sphere-cone at 15 deg angle of attack, and the results are compared with an inviscid method-of-characteristics calculation.

Rakich, J. V.↗

NASA NTRS2004

On the Inversion for Mass (Re)Distribution from Global (Time-Variable) Gravity Field

The well-known non-uniqueness of the gravitational inverse problem states the following: The external gravity field, even if completely and exactly known, cannot Uniquely determine the density distribution of the body that produces the gravity field. This is an intrinsic property of a field that obeys the Laplace equation, as already treated in mathematical as well as geophysical literature. In this paper we provide conceptual insight by examining the problem in terms of spherical harmonic expansion of the global gravity field. By comparing the multipoles and the moments of the density function, we show that in 3-S the degree of knowledge deficiency in trying to inversely recover the density distribution from external gravity field is (n+l)(n+2)/2 - (2n+l) = n(n-1)/2 for each harmonic degree n. On the other hand, on a 2-D spherical shell we show via a simple relationship that the inverse solution of the surface density distribution is unique. The latter applies quite readily in the inversion of time-variable gravity signals (such as those observed by the GRACE space mission) where the sources over a wide range of the scales largely come from the Earth's Surface.

Chao, Benjamin F.↗

NASA NTRS2000

Resin Film Infusion (RFI) Process Modeling for Large Transport Aircraft Wing Structures

This investigation completed the verification of a three-dimensional resin transfer molding/resin film infusion (RTM/RFI) process simulation model. The model incorporates resin flow through an anisotropic carbon fiber preform, cure kinetics of the resin, and heat transfer within the preform/tool assembly. The computer model can predict the flow front location, resin pressure distribution, and thermal profiles in the modeled part. The formulation for the flow model is given using the finite element/control volume (FE/CV) technique based on Darcy's Law of creeping flow through a porous media. The FE/CV technique is a numerically efficient method for finding the flow front location and the fluid pressure. The heat transfer model is based on the three-dimensional, transient heat conduction equation, including heat generation. Boundary conditions include specified temperature and convection. The code was designed with a modular approach so the flow and/or the thermal module may be turned on or off as desired. Both models are solved sequentially in a quasi-steady state fashion. A mesh refinement study was completed on a one-element thick model to determine the recommended size of elements that would result in a converged model for a typical RFI analysis. Guidelines are established for checking the convergence of a model, and the recommended element sizes are listed. Several experiments were conducted and computer simulations of the experiments were run to verify the simulation model. Isothermal, non-reacting flow in a T-stiffened section was simulated to verify the flow module. Predicted infiltration times were within 12-20% of measured times. The predicted pressures were approximately 50% of the measured pressures. A study was performed to attempt to explain the difference in pressures. Non-isothermal experiments with a reactive resin were modeled to verify the thermal module and the resin model. Two panels were manufactured using the RFI process. One was a stepped panel and the other was a panel with two 'T' stiffeners. The difference between the predicted infiltration times and the experimental times was 4% to 23%.

Loos, Alfred C.↗

NASA NTRS1984

The Stability of the Low Degree 5 Minute Solar Oscillations

The decay rate for many of the low degree p modes observed as 5 minute oscillations of the Sun is discussed and decay rates for modes ranging from radial (l=0) to the nonradial ones with l=5 for overtones 10 through 28. Parameters needed for he solar model are given and the hydrogen mass fraction composition structure is shown as well as the structure given by Christensen-Dalsgaard (1982) for an evolved solar model. The special equation of state and opacity table with X=0.74 for the hydrogen mass fraction in the outer 0.40 of the mass needs slightly more hydrogen in the central regions than obtained by Christensen-Dalsgaard in order to give a complete and consistent model. The difference in helium production between these two models is about 10%, meaning that the total energy radiated by the Sun during its lifetime thus far agrees satisfactorily with accurately calculated evolution sequences.

Kidman, R. B.↗

NASA NTRS1987

Finite Element Flux-Corrected Transport (FEM-FCT) for the Euler and Navier-Stokes equations

A high resolution finite element method for the solution of problems involving high speed compressible flows is presented. The method uses the concepts of flux-corrected transport and is presented in a form which is suitable for implementation on completely unstructured triangular or tetrahedral meshes. Transient and steady state examples are solved to illustrate the performance of the algorithm.

Loehner, Rainald↗

NASA NTRS1987

Finite element flux-corrected transport (FEM-FCT) for the Euler and Navier-Stokes equations

A high resolution finite element method for the solution of problems involving high speed compressible flows is presented. The method uses the concepts of flux-corrected transport and is presented in a form which is suitable for implementation on completely unstructured triangular or tetrahedral meshes. Transient and steady-state examples are solved to illustrate the performance of the algorithm.

Lohner, Rainald↗

NASA NTRS1998

Ultrafast Modulation and Switching of Quantum-Well Lasers using Terahertz Fields

Modulation and switching of semiconductor lasers are important for laser-based information technology. Typically the speed of modulation and switching is limited by interband processes such as stimulated and spontaneous recombinations which occur on a nanosecond time scale. This is why the diode laser modulation has been restricted to tens of GHz. Modulation at higher speed is highly desirable as the information technology enters into the so-called tera-era. In this paper, we study the possibility of utilizing THz-field-induced plasma heating to modulate quantum-well lasers. This is a timely study since, with the advancement of THz solid-state sources and free-electron lasers, THz physics and related technology is currently coming out of its infancy. The investigation of interplaying THz and optical fields is also of intruiging fundamental interest. First, we introduce theoretical plasma heating results for the quantum-well optical amplifier in the presense of an intense half-cycle THz pulse. The heated carrier distributions are then utilized to calculate the THz-pulse-induced change in refractive index and gain profile. Since the electron-hole-plasma is heated using intraband transitions, we circumvent the usual complications due to an overall change in density, and the nonlinear recovery is governed solely by the carrier-LO-phonon interactions, typically 5 ps for a complete recovery. This procedure implies THz and sub-THz switching and recovery rates, respectively; using either gain modulation or index modulation. Plasma heating via steady-state THz fields is also studied. Finally, numerical simulation of a coupled set of equations to investigate the THz modulation based on a simplified model for quantum-well lasers is presented. Our results show that a semiconductor laser can be modulated at up to 1 THz with little distortion with a THz field amplitude at the order of a few kV/cm. Laser responses to a change in THz frequency will be shown. Constraints, practicalities, and applications will be discussed.

Ning, Cun-Zheng↗

NASA NTRS1974

Method of Characteristics Calculations and Computer Code for Materials with Arbitrary Equations of State and Using Orthogonal Polynomial Least Square Surface Fits

A numerical scheme using the method of characteristics to calculate the flow properties and pressures behind decaying shock waves for materials under hypervelocity impact is developed. Time-consuming double interpolation subroutines are replaced by a technique based on orthogonal polynomial least square surface fits. Typical calculated results are given and compared with the double interpolation results. The complete computer program is included.

Chang, T. S.↗

NASA NTRS1995

The history of re-connection and the concept of the solar wind plasma with relatively small electrical conductivity

Petschek's 're-connection' model, aspiring to be universal, treated as a boundary problem meets unresolvable difficulties connected with impossibility to specify correctly boundary and initial conditions. This problem was incorrectly formulated. Hence, ineradicable logarithmic singularities occurred on the boundary surfaces. Attempts to eliminate them by incorporating the finite electrical conductivity are incorrect. This should lead to the change in the equation type, boundary condition type and in consequence to the change in solutions. Besides, the slow mode shocks cannot be driven by small internal source. As an alternative a new plasma concept is suggested. The state of fully ionized plasma in space depends completely on the entropy of the plasma heating source and on the process in which plasma is involved. The presumptive source of the solar wind creation - the induction electric field of the solar origin - has very low entropy. The state of plasma should be very far from the thermodynamic equilibrium. Debye's screening is not complete. The excitation of the powerful resonant self-consistent electric fields in plasma provides low electric conductivity. The MHD problems should be treated in frameworks of dissipative theories.

Chertkov, A. D.↗

NASA NTRS1988

A Mathematical Formulation of the SCOLE Control Problem. Part 2: Optimal Compensator Design

The study initiated in Part 1 of this report is concluded and optimal feedback control (compensator) design for stability augmentation is considered, following the mathematical formulation developed in Part 1. Co-located (rate) sensors and (force and moment) actuators are assumed, and allowing for both sensor and actuator noise, stabilization is formulated as a stochastic regulator problem. Specializing the general theory developed by the author, a complete, closed form solution (believed to be new with this report) is obtained, taking advantage of the fact that the inherent structural damping is light. In particular, it is possible to solve in closed form the associated infinite-dimensional steady-state Riccati equations. The SCOLE model involves associated partial differential equations in a single space variable, but the compensator design theory developed is far more general since it is given in the abstract wave equation formulation. The results thus hold for any multibody system so long as the basic model is linear.

Balakrishnan, A. V.↗

NASA NTRS1994

Crystal growth and fluid mechanics problems in directional solidification

Broadly speaking, our efforts have been concentrated in two aspects of directional solidification: (A) a more complete theoretical understanding of convection effects in a Bridgman apparatus; and (B) a clear understanding of scalings of various features of dendritic crystal growth in the sensitive limit of small capillary effects. For studies that fall within class A, the principal objectives are as follows: (A1) Derive analytical formulas for segregation, interfacial shape and fluid velocities in mathematically amenable asymptotic limits. (A2) Numerically verify and extend asymptotic results to other ranges of parameter space with a view to a broader physical understanding of the general trends. With respect to studies that fall within class B, the principal objectives include answering the following questions about dendritic crystal growth: (B1) Are there unsteady dendrite solutions in 2-D to the completely nonlinear time evolving equations in the small surface tension limit with only a locally steady tip region with well defined tip radius and velocity? Is anisotropy in surface tension necessary for the existence of such solutions as it is for a true steady state needle crystal? How does the size of such a local region depend on capillary effects, anisotropy and undercooling? (B2) How do the different control parameters affect the nonlinear amplification of tip noise and dendritic side branch coarsening?

Tanveer, Saleh↗

NASA NTRS1987

Dynamics formulations for the real-time simulation of constrained motion

The Space Shuttle program has relied heavily on simulation throughout all phases of development and operation. Real-time, man-in-the-loop simulation has served the NASA manned space flight program by providing the means to evaluate systems design and integrated systems performance in a simulated flight environment as well as provide a means to train flight crews. New challenges are presented by the development and operation of a permanently manned space station. The assembly of the space station, the transferral of payloads and the use of the space station manipulator to berth the Orbiter are operations critical to the success of the space station. All these operations are examples of constrained motion among the bodies associated with the Orbiter and space station system. The state of the art of formulating the governing dynamical equations of motion for constrained systems are described. The uses of the two basic problems in multibody dynamics are discussed. The most efficient formulations of the equations of motion are addressed from the point of view of completeness. The issues surrounding incorporating the constraints into the equation of motion are presented.

Kelly, Frederick A.↗