SEARCH · Engineering Papers

Results for “Finite Volume”

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 1,495 records · Page 83

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Interplanetary Program to Optimize Simulated Trajectories (IPOST). Volume 2: Analytic manual

The Interplanetary Program to Optimize Space Trajectories (IPOST) is intended to support many analysis phases, from early interplanetary feasibility studies through spacecraft development and operations. The IPOST output provides information for sizing and understanding mission impacts related to propulsion, guidance, communications, sensor/actuators, payload, and other dynamic and geometric environments. IPOST models three degree of freedom trajectory events, such as launch/ascent, orbital coast, propulsive maneuvering (impulsive and finite burn), gravity assist, and atmospheric entry. Trajectory propagation is performed using a choice of Cowell, Encke, Multiconic, Onestep, or Conic methods. The user identifies a desired sequence of trajectory events, and selects which parameters are independent (controls) and dependent (targets), as well as other constraints and the cost function. Targeting and optimization is performed using the Stanford NPSOL algorithm. IPOST structure allows subproblems within a master optimization problem to aid in the general constrained parameter optimization solution. An alternate optimization method uses implicit simulation and collocation techniques.

Hong, P. E.↗

NASA NTRS2002

Efficient Reformulation of the Thermoelastic Higher-order Theory for Fgms

Functionally graded materials (FGMs) are characterized by spatially variable microstructures which are introduced to satisfy given performance requirements. The microstructural gradation gives rise to continuously or discretely changing material properties which complicate FGM analysis. Various techniques have been developed during the past several decades for analyzing traditional composites and many of these have been adapted for the analysis of FGMs. Most of the available techniques use the so-called uncoupled approach in order to analyze graded structures. These techniques ignore the effect of microstructural gradation by employing specific spatial material property variations that are either assumed or obtained by local homogenization. The higher-order theory for functionally graded materials (HOTFGM) is a coupled approach developed by Aboudi et al. (1999) which takes the effect of microstructural gradation into consideration and does not ignore the local-global interaction of the spatially variable inclusion phase(s). Despite its demonstrated utility, however, the original formulation of the higher-order theory is computationally intensive. Herein, an efficient reformulation of the original higher-order theory for two-dimensional elastic problems is developed and validated. The use of the local-global conductivity and local-global stiffness matrix approach is made in order to reduce the number of equations involved. In this approach, surface-averaged quantities are the primary variables which replace volume-averaged quantities employed in the original formulation. The reformulation decreases the size of the global conductivity and stiffness matrices by approximately sixty percent. Various thermal, mechanical, and combined thermomechanical problems are analyzed in order to validate the accuracy of the reformulated theory through comparison with analytical and finite-element solutions. The presented results illustrate the efficiency of the reformulation and its advantages in analyzing functionally graded materials.

Bansal, Yogesh↗

NASA NTRS2007

Lunar Lander Structural Design Studies at NASA Langley

The National Aeronautics and Space Administration is currently developing mission architectures, vehicle concepts and flight hardware to support the planned human return to the Moon. During Phase II of the 2006 Lunar Lander Preparatory Study, a team from the Langley Research Center was tasked with developing and refining two proposed Lander concepts. The Descent-Assisted, Split Habitat Lander concept uses a disposable braking stage to perform the lunar orbit insertion maneuver and most of the descent from lunar orbit to the surface. The second concept, the Cargo Star Horizontal Lander, carries ascent loads along its longitudinal axis, and is then rotated in flight so that its main engines (mounted perpendicular to the vehicle longitudinal axis) are correctly oriented for lunar orbit insertion and a horizontal landing. Both Landers have separate crew transport volumes and habitats for surface operations, and allow placement of large cargo elements very close to the lunar surface. As part of this study, lightweight, efficient structural configurations for these spacecraft were proposed and evaluated. Vehicle structural configurations were first developed, and preliminary structural sizing was then performed using finite element-based methods. Results of selected structural design and trade studies performed during this activity are presented and discussed.

Wu, K. Chauncey↗

NASA NTRS1992

Phase noise in pulsed Doppler lidar and limitations on achievable single-shot velocity accuracy

The smaller sampling volumes afforded by Doppler lidars compared to radars allows for spatial resolutions at and below some sheer and turbulence wind structure scale sizes. This has brought new emphasis on achieving the optimum product of wind velocity and range resolutions. Several recent studies have considered the effects of amplitude noise, reduction algorithms, and possible hardware related signal artifacts on obtainable velocity accuracy. We discuss here the limitation on this accuracy resulting from the incoherent nature and finite temporal extent of backscatter from aerosols. For a lidar return from a hard (or slab) target, the phase of the intermediate frequency (IF) signal is random and the total return energy fluctuates from shot to shot due to speckle; however, the offset from the transmitted frequency is determinable with an accuracy subject only to instrumental effects and the signal to noise ratio (SNR), the noise being determined by the LO power in the shot noise limited regime. This is not the case for a return from a media extending over a range on the order of or greater than the spatial extent of the transmitted pulse, such as from atmospheric aerosols. In this case, the phase of the IF signal will exhibit a temporal random walk like behavior. It will be uncorrelated over times greater than the pulse duration as the transmitted pulse samples non-overlapping volumes of scattering centers. Frequency analysis of the IF signal in a window similar to the transmitted pulse envelope will therefore show shot-to-shot frequency deviations on the order of the inverse pulse duration reflecting the random phase rate variations. Like speckle, these deviations arise from the incoherent nature of the scattering process and diminish if the IF signal is averaged over times greater than a single range resolution cell (here the pulse duration). Apart from limiting the high SNR performance of a Doppler lidar, this shot-to-shot variance in velocity estimates has a practical impact on lidar design parameters. In high SNR operation, for example, a lidar's efficiency in obtaining mean wind measurements is determined by its repetition rate and not pulse energy or average power. In addition, this variance puts a practical limit on the shot-to-shot hard target performance required of a lidar.

Mcnicholl, P.↗

NASA NTRS2001

Visualization of CFD Results in Immersive Virtual Environments

An object-oriented event-driven immersive virtual environment (VE) is described for the visualization of computational fluid dynamics (CFD) results. The VE incorporates the following types of primitive software objects: interface objects, support objects, geometric entities, and finite elements. The fluid domain is discretized using either a multi-block structured grid or an unstructured finite element mesh. The VE allows natural 'fly-through' visualization of the model, the CFD grid, and the model's surroundings. In order to help visualize the flow and its effects on the model, the VE incorporates the following objects: stream objects (lines, surface-restricted lines. ribbons. and volumes); colored surfaces; elevation surfaces; surface arrows; global and local iso-surfaces; vortex cores; and separation/attachment surfaces and lines. Most of these objects can be used for dynamically probing the flow. Particles and arrow animations can be displayed on top of stream objects. Primitive response quantities as well as derived quantities can be used. A recursive tree search algorithm is used for real-time point and value search in the CFD grid.

Wasfy, Tamer M.↗

NASA NTRS2004

Further Improvement in 3DGRAPE

3DGRAPE/AL:V2 denotes version 2 of the Three-Dimensional Grids About Anything by Poisson's Equation with Upgrades from Ames and Langley computer program. The preceding version, 3DGRAPE/AL, was described in Improved 3DGRAPE (ARC-14069) NASA Tech Briefs, Vol. 21, No. 5 (May 1997), page 66. These programs are so named because they generate volume grids by iteratively solving Poisson's Equation in three dimensions. The grids generated by the various versions of 3DGRAPE have been used in computational fluid dynamics (CFD). The main novel feature of 3DGRAPE/AL:V2 is the incorporation of an optional scheme in which anisotropic Lagrange-based trans-finite interpolation (ALBTFI) is coupled with exponential decay functions to compute and blend interior source terms. In the input to 3DGRAPE/AL:V2 the user can specify whether or not to invoke ALBTFI in combination with exponential-decay controls, angles, and cell size for controlling the character of grid lines. Of the known programs that solve elliptic partial differential equations for generating grids, 3DGRAPE/AL:V2 is the only code that offers a combination of speed and versatility with most options for controlling the densities and other characteristics of grids for CFD.

Alter, Stephen↗

NASA NTRS1982

Numerical boundary condition procedures and multigrid methods; Proceedings of the Symposium, NASA Ames Research Center, Moffett Field, CA, October 19-22, 1981

Papers presented in this volume provide an overview of recent work on numerical boundary condition procedures and multigrid methods. The topics discussed include implicit boundary conditions for the solution of the parabolized Navier-Stokes equations for supersonic flows; far field boundary conditions for compressible flows; and influence of boundary approximations and conditions on finite-difference solutions. Papers are also presented on fully implicit shock tracking and on the stability of two-dimensional hyperbolic initial boundary value problems for explicit and implicit schemes.

Source record↗

NASA NTRS1988

The critical complex in nucleation is not the nucleus

A theory of nucleation is presented which combines features of classical nucleation. In this theory, the nucleus is a cluster of finite size with a sharp interface and spinodal decomposition, in which a fluctuation in density or composition generates the critical complex and nucleus. Thus, the nucleus is a product of nucleation. Schematic models for the classical and the proposed routes to nucleation are compared, focusing on the relation between pressure and volume in the models. In the proposed model, the critical complex, or critical state, is between the parent state and the nucleus.

Rasmussen, Don H.↗

NASA NTRS1993

The effect of nonsymmetric pressure stiffness on the dynamic characteristics of Solid Rocket Motors

This paper discusses the effect of pressure on the dynamics of pre-stiffened structures such as the Advanced Solid Rocket Motor (ASRM). Previous work in which the stiffness terms resulting from constant pressure were derived has been extended to enable modeling of nonconstant pressure applied over nonenclosed volumes. These conditions will result in nonsymmetric terms in the global stiffness matrix which will not cancel out. Three new pressure stiffness elements incorporating these nonsymmetric terms have been implemented as dummy elements in COSMIC NASTRAN and have been tested on various simple examples as well as an existing ASRM NASTRAN finite element model. The results indicate that for all load cases of practical interest to the ASRM program, the nonsymmetric terms have very little effect on the dynamic characteristics. In addition, the pressure stiffness elements developed in the previous work which assumed constant pressure gave virtually the same results as the new elements even for problems in which the pressures are not constant. The original elements appear to work well as long as the pressure gradient across any individual element is no larger than about 0.75 psi/inch. The new elements are therefore most useful for determining the conditions under which the original pressure stiffness elements can be used.

Christensen, Eric R.↗

NASA NTRS1993

Buckling analysis of laminated thin shells in a hot environment

Results are presented of parametric studies to assess the effects of various parameters on the buckling behavior of angle-ply, laminated thin shells in a hot environment. These results were obtained by using a three-dimensional finite element analysis. An angle-ply, laminated thin shell with fiber orientation of (theta/-theta)(sub 2) was subjected to compressive mechanical loads. The laminated thin shell had a cylindrical geometry. The laminate contained T300 graphite fibers embedded in an intermediate-modulus, high-strength (IMHS) matrix. The fiber volume fraction was 55 percent and the moisture content was 2 percent. The residual stresses induced into the laminate structure during the curing were taken into account. Parametric studies were performed to examine the effect on the critical buckling load of the following parameters: cylinder length and thickness, internal hydrostatic pressure, different ply thicknesses, different temperature profiles through the thickness of the structure, and different lay up configurations and fiber volume fractions. In conjunction with these parameters the ply orientation was varied from 0 deg to 90 deg. Seven ply angles were examined: 0 deg, 15 deg, 30 deg, 45 deg, 60 deg, 75 deg, and 90 deg. The results show that the ply angle theta and the laminate thickness had significant effects on the critical buckling load. The fiber volume fraction, the fiber orientations, and the internal hydrostatic pressure had important effects on the critical buckling load. The cylinder length had a moderate influence on the buckling load. The thin shell with (theta/-theta)(sub 2) or (theta/-theta)(sub s) angle-ply laminate had better buckling-load performance than the thin shell with (theta)(sub 4) angle-ply laminate. The temperature profiles through the laminate thickness and various laminates with the different ply thicknesses has insignificant effects on the buckling behavior of the thin shells.

Gotsis, Pascal K.↗

NASA NTRS1999

Maximum Energies of Shock-Accelerated Electrons in Young Shell Supernova Remnants

Young supernova remnants (SNRs) are often assumed to be the source of cosmic rays up to energies approaching the slight steepening in the cosmic ray spectrum at around 1000 TeV, known as the "knee." We show that the observed X-ray emission of 14 radio-bright shell remnants, including all five historical shells, can be used to put limits on E(sub max), the energy at which the electron energy distribution must steepen from its slope at radio-emitting energies. Most of the remnants show thermal spectra, so any synchrotron component must fall below the observed X-ray fluxes. We obtain upper limits on E(sub max) by considering the most rapid physically plausible cutoff in the relativistic electron distribution, an exponential, which is as sharp or sharper than found in any more elaborate models. This maximally curved model then gives us the highest possible E(sub max) consistent with not exceeding observed X-rays. Our results are thus independent of particular models for the electron spectrum in SNRs. Assuming homogeneous emitting volumes with a constant magnetic field strength of 10 uG, no object could reach 1000 TeV, and only one, Kes 73, has an upper limit on E(sub max), above 100 TeV. All the other remnants have limits at or below 80 TeV. E(sub max) is probably set by the finite remnant lifetime rather than by synchrotron losses for remnants younger than a few thousand years, so that an observed electron steepening should be accompanied by steepening at the same energy for protons. More complicated, inhomogeneous models could allow higher values of E(sub max) in parts of the remnant, but the emission-weighted average value, that characteristic of typical electrons, should obey these limits. The young remnants are not expected to improve much over their remaining lives at producing the highest energy Galactic cosmic rays; if they cannot, this picture of cosmic-ray origin may need major alteration.

Reynolds, Stephen P.↗

NASA NTRS1972

Thermal elastoplastic structural analysis of non-metallic thermal protection systems

An incremental theory and numerical procedure to analyze a three-dimensional thermoelastoplastic structure subjected to high temperature, surface heat flux, and volume heat supply as well as mechanical loadings are presented. Heat conduction equations and equilibrium equations are derived by assuming a specific form of incremental free energy, entropy, stresses and heat flux together with the first and second laws of thermodynamics, von Mises yield criteria and Prandtl-Reuss flow rule. The finite element discretization using the linear isotropic three-dimensional element for the space domain and a difference operator corresponding to a linear variation of temperature within a small time increment for the time domain lead to systematic solutions of temperature distribution and displacement and stress fields. Various boundary conditions such as insulated surfaces and convection through uninsulated surface can be easily treated. To demonstrate effectiveness of the present formulation a number of example problems are presented.

Chung, T. J.↗

NASA NTRS1988

Structural analyses of engine wall cooling concepts and materials

The severe thermal environments under which hypersonic aircraft such as the National Aerospace Plane (NASP) will operate require cooling of the engine walls, especially in the combustor. A preliminary assessment is made of some candidate materials based on structural analyses for a number of convective cooling configurations. Three materials were studied: graphite/copper and tungsten/copper composite alloys with 50 percent fiber volume fractions and a wrought cobalt-base superalloy, Haynes 188. Anisotropic mechanical and thermal properties for the composites were obtained from a computer code, ICAN, which determines the composite material properties from the individual properties of the fiber and matrix materials. The structural analyses were performed by using the MARC nonlinear finite element code. Heat transfer analyses were conducted to calculate the metal temperature distributions.

Kaufman, Albert↗

NASA NTRS1997

Complex Volume Grid Generation Through the Use of Grid Reusability

This paper presents a set of surface and volume grid generation techniques which reuse existing surface and volume grids. These methods use combinations of data manipulations to reduce grid generation time, improve grid characteristics, and increase the capabilities of existing domain discretization software. The manipulation techniques utilize physical and computational domains to produce basis function on which to operate and modify grid character and smooth grids using Trans-Finite Interpolation, a vector interpolation method and parametric re-mapping technique. With these new techniques, inviscid grids can be converted to viscous grids, multiple zone grid adaption can be performed to improve CFD solver efficiency, and topological changes to improve modeling of flow fields can be done simply and quickly. Examples of these capabilities are illustrated as applied to various configurations.

Alter, Stephen J.↗

NASA NTRS1993

Coupled structural, thermal, phase-change and electromagnetic analysis for superconductors, volume 1

This research program has dealt with the theoretical development and computer implementation of reliable and efficient methods for the analysis of coupled mechanical problems that involve the interaction of mechanical, thermal, phase-change and electromagnetic subproblems. The focus application has been the modeling of superconductivity and associated quantum-state phase-change phenomena. In support of this objective the work has addressed the following issues: (1) development of variational principles for finite elements; (2) finite element modeling of the electromagnetic problem; (3) coupling of thermal and mechanical effects; and (4) computer implementation and solution of the superconductivity transition problem. The research was carried out over the period September 1988 through March 1993. The main accomplishments have been: (1) the development of the theory of parametrized and gauged variational principles; (2) the application of those principled to the construction of electromagnetic, thermal and mechanical finite elements; and (3) the coupling of electromagnetic finite elements with thermal and superconducting effects; and (4) the first detailed finite element simulations of bulk superconductors, in particular the Meissner effect and the nature of the normal conducting boundary layer. The grant has fully supported the thesis work of one doctoral student (James Schuler, who started on January 1989 and completed on January 1993), and partly supported another thesis (Carmelo Militello, who started graduate work on January 1988 completing on August 1991). Twenty-three publications have acknowledged full or part support from this grant, with 16 having appeared in archival journals and 3 in edited books or proceedings.

Felippa, Carlos A.↗

NASA NTRS1989

Comet ejection and dynamics of nonspherical dust particles and meteoroids

This paper generalizes the formalism for calculating the ejection velocity of meteoroids and dust from comets and the forces to which such objects are subject in interplanetary space, including the dust tail of comets. It is found that spheres have the smallest cross section of any geometrical figures of the same valume averaged over random orientations, so for a fixed volume and mass, both the ejection velocity and beta reaches a minimum for bodies of spherical shapes. Flakes in random orientation are ejected near 70 percent of the highest ejection velocity for any orientation. Needles in random orientation escape a comet at nearly 90 percent of their maximum velocities. Randomly oriented cylinders of finite thickness escape at lower velocities that are slightly closer to their maximum velocities. The average beta acting on spin-aligned, perfectly absorbing needles is more than half that acting on a sphere of the same material and radius.

Gustafson, Bo A. S.↗

NASA NTRS1992

International Congress on Recent Developments in Air- and Structure-Borne Sound and Vibration, 2nd, Auburn Univ., AL, Mar. 4-6, 1992, Proceedings. Vols. 1-3

The papers presented in this volume cover a variety of topics, including sound intensity, structural intensity, modal analysis and synthesis, statistical energy analysis and energy methods, passive and active damping, and boundary element methods. Attention is also given to diagnostics and condition monitoring, material characterization and nondestructive evaluation, active noise and vibration control, sound radiation and scattering, and finite element analysis.

Crocker, Malcolm J.↗

NASA NTRS1994

Technical report series on global modeling and data assimilation. Volume 2: Direct solution of the implicit formulation of fourth order horizontal diffusion for gridpoint models on the sphere

High order horizontal diffusion of the form K Delta(exp 2m) is widely used in spectral models as a means of preventing energy accumulation at the shortest resolved scales. In the spectral context, an implicit formation of such diffusion is trivial to implement. The present note describes an efficient method of implementing implicit high order diffusion in global finite difference models. The method expresses the high order diffusion equation as a sequence of equations involving Delta(exp 2). The solution is obtained by combining fast Fourier transforms in longitude with a finite difference solver for the second order ordinary differential equation in latitude. The implicit diffusion routine is suitable for use in any finite difference global model that uses a regular latitude/longitude grid. The absence of a restriction on the timestep makes it particularly suitable for use in semi-Lagrangian models. The scale selectivity of the high order diffusion gives it an advantage over the uncentering method that has been used to control computational noise in two-time-level semi-Lagrangian models.

Max J. Suarez↗