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,423 records · Page 79

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Finite difference time domain modeling of steady state scattering from jet engines with moving turbine blades

The approach chosen to model steady state scattering from jet engines with moving turbine blades is based upon the Finite Difference Time Domain (FDTD) method. The FDTD method is a numerical electromagnetic program based upon the direct solution in the time domain of Maxwell's time dependent curl equations throughout a volume. One of the strengths of this method is the ability to model objects with complicated shape and/or material composition. General time domain functions may be used as source excitations. For example, a plane wave excitation may be specified as a pulse containing many frequencies and at any incidence angle to the scatterer. A best fit to the scatterer is accomplished using cubical cells in the standard cartesian implementation of the FDTD method. The material composition of the scatterer is determined by specifying its electrical properties at each cell on the scatterer. Thus, the FDTD method is a suitable choice for problems with complex geometries evaluated at multiple frequencies. It is assumed that the reader is familiar with the FDTD method.

Ryan, Deirdre A.↗

NASA NTRS1992

The effects of uneven fiber spacing on thermal residual stresses in a unidirectional SCS-6/Ti-15-3 laminate

High residual stresses develop in SCS-6/Ti-15-3 composites during cooldown from the fabrication temperature; these residual stresses can effect the mechanical and physical properties of the composite. Discrete fiber-matrix finite element models were used to study the residual stresses due to the temperature change during the fabrication process, including the effects of uneven fiber spacing, the free surface, and increased fiber volume fractions. To accurately model the effects of the free surface, it is only necessary to model one fiber through the thickness. Below the first ply, the analysis predicts stress distributions that are identical to the infinite array predictions. For uneven fiber space less than 0.042 mm in an interior ply, the maximum hoop stress was predicted to occur between fibers within a ply and to increase as the fiber spacing decreased. The maximum hoop stress correlated well with the observed radial cracking between fibers. For the case of touching fibers, the analysis predicted tensile radial stresses at the fiber-matrix debonding during the fabrication cooldown. Identical trends were predicted for uneven fiber spacing in surface plies with slightly greater values of maximum stresses. The analysis predicted matrix yielding to occur upon cooldown when the edge-to-edge fiber spacing was less than or equal to 0.022 mm. The stress distributions preidcted for increasing fiber volume fractions were similar to those predicted for decreasing the fiber spacing for two adjacent fibers within a ply.

Bigelow, Catherine A.↗

NASA NTRS1999

Conformal Cryogenic Tank Trade Study for Reusable Launch Vehicles

Future reusable launch vehicles may be lifting bodies with non-circular cross section like the proposed Lockheed-Martin VentureStar(tm). Current designs for the cryogenic tanks of these vehicles are dual-lobed and quad-lobed tanks which are packaged more efficiently than circular tanks, but still have low packaging efficiencies with large gaps existing between the vehicle outer mold line and the outer surfaces of the tanks. In this study, tanks that conform to the outer mold line of a non-circular vehicle were investigated. Four structural concepts for conformal cryogenic tanks and a quad-lobed tank concept were optimized for minimum weight designs. The conformal tank concepts included a sandwich tank stiffened with axial tension webs, a sandwich tank stiffened with transverse tension webs, a sandwich tank stiffened with rings and tension ties, and a sandwich tank stiffened with orthogrid stiffeners and tension ties. For each concept, geometric parameters (such as ring frame spacing, the number and spacing of tension ties or webs, and tank corner radius) and internal pressure loads were varied and the structure was optimized using a finite-element-based optimization procedure. Theoretical volumetric weights were calculated by dividing the weight of the barrel section of the tank concept and its associated frames, webs and tension ties by the volume it circumscribes. This paper describes the four conformal tank concepts and the design assumptions utilized in their optimization. The conformal tank optimization results included theoretical weights, trends and comparisons between the concepts, are also presented, along with results from the optimization of a quad-lobed tank. Also, the effects of minimum gauge values and non-optimum weights on the weight of the optimized structure are described in this paper.

Rivers, H. Kevin↗

NASA NTRS1975

A finite element formulation for evaluation of crack blunting effects in elasto-plastic solids

The sharp crack model allows the use of continuum analysis, linear elasticity, as a basic for predicting a micromechanical process, fracture, by providing a characterization of loading conditions affecting a very small volume of material. The model is described in detail.

Osias, J. R.↗

NASA NTRS1996

Measurements and Modeling of Soot Formation and Radiation in Microgravity Jet Diffusion Flames

This is a computational and experimental study for soot formation and radiative heat transfer in jet diffusion flames under normal gravity (1-g) and microgravity (0-g) conditions. Instantaneous soot volume fraction maps are measured using a full-field imaging absorption technique developed by the authors. A compact, self-contained drop rig is used for microgravity experiments in the 2.2-second drop tower facility at NASA Lewis Research Center. On modeling, we have coupled flame structure and soot formation models with detailed radiation transfer calculations. Favre-averaged boundary layer equations with a k-e-g turbulence model are used to predict the flow field, and a conserved scalar approach with an assumed Beta-pdf are used to predict gaseous species mole fraction. Scalar transport equations are used to describe soot volume fraction and number density distributions, with formation and oxidation terms modeled by one-step rate equations and thermophoretic effects included. An energy equation is included to couple flame structure and radiation analyses through iterations, neglecting turbulence-radiation interactions. The YIX solution for a finite cylindrical enclosure is used for radiative heat transfer calculations. The spectral absorption coefficient for soot aggregates is calculated from the Rayleigh solution using complex refractive index data from a Drude- Lorentz model. The exponential-wide-band model is used to calculate the spectral absorption coefficient for H20 and C02. It is shown that when compared to results from true spectral integration, the Rosseland mean absorption coefficient can provide reasonably accurate predictions for the type of flames studied. The soot formation model proposed by Moss, Syed, and Stewart seems to produce better fits to experimental data and more physically sound than the simpler model by Khan et al. Predicted soot volume fraction and temperature results agree well with published data for a normal gravity co-flow laminar flames and turbulent jet flames. Predicted soot volume fraction results also agree with our data for 1-g and 0-g laminar jet names as well as 1-g turbulent jet flames.

Ku, Jerry C.↗

NASA NTRS1993

An efficient formulation of the coupled finite element-integral equation technique for solving large 3D scattering problems

It is often desirable to calculate the electromagnetic fields inside and about a complicated system of scattering bodies, as well as in their far-field region. The finite element method (FE) is well suited to solving the interior problem, but the domain has to be limited to a manageable size. At the truncation of the FE mesh one can either impose approximate (absorbing) boundary conditions or set up an integral equation (IE) for the fields scattered from the bodies. The latter approach is preferable since it results in higher accuracy. Hence, the two techniques can be successfully combined by introducing a surface that encloses the scatterers, applying a FE model to the inner volume and setting up an IE for the tangential fields components on the surface. Here the continuity of the tangential fields is used bo obtain a consistent solution. A few coupled FE-IE methods have recently appeared in the literature. The approach presented here has the advantage of using edge-based finite elements, a type of finite elements with degrees of freedom associated with edges of the mesh. Because of their properties, they are better suited than the conventional node based elements to represent electromagnetic fields, particularly when inhomogeneous regions are modeled, since the node based elements impose an unnatural continuity of all field components across boundaries of mesh elements. Additionally, our approach is well suited to handle large size problems and lends itself to code parallelization. We will discuss the salient features that make our approach very efficient from the standpoint of numerical computation, and the fields and RCS of a few objects are illustrated as examples.

Cwik, T.↗

NASA NTRS1995

Technical report series on global modeling and data assimilation. Volume 5: Documentation of the AIRES/GEOS dynamical core, version 2

A detailed description of the numerical formulation of Version 2 of the ARIES/GEOS 'dynamical core' is presented. This code is a nearly 'plug-compatible' dynamics for use in atmospheric general circulation models (GCMs). It is a finite difference model on a staggered latitude-longitude C-grid. It uses second-order differences for all terms except the advection of vorticity by the rotation part of the flow, which is done at fourth-order accuracy. This dynamical core is currently being used in the climate (ARIES) and data assimilation (GEOS) GCMs at Goddard.

Suarez, Max J.↗

NASA NTRS1986

A scattering model for defoliated vegetation

A scattering model for defoliated vegetation is conceived as a layer of dielectric, finite-length cylinders with specified size and orientation distributions above an irregular ground surface. The scattering phase matrix of a single cylinder is computed, then the radiative transfer technique is applied to link volume scattering from vegetation to surface scattering from the soil surface. Polarized and depolarized scattering are computed and the effects of the cylinder size and orientation distributions are illustrated. It is found that size and orientation distributions have significant effects on the backscattered signal. The model is compared with scattering from defoliated trees and agricultural crops.

Karam, M. A.↗

NASA NTRS1989

Computational simulation of damping in composite structures

A computational methodology is developed for the prediction of passive damping in composite structures. The method involves multiple levels of damping modeling by integrating micromechanics, laminate, and structural damping theories. The effects of temperature and moisture on structural damping are included. The simulation of damping in the structural level is accomplished with finite-element discretization. Applications are performed on graphite/epoxy composite beams, plates, and shells to illustrate the methodology. Additional parametric studies demonstrate the variation of structural modal damping with respect to ply angles, fiber volume ratio, and temperature.

Saravanos, D. A.↗

NASA NTRS1991

Computational simulation of damping in composite structures

Saravanos, D. A.↗

NASA NTRS1991

Design of an essentially non-oscillatory reconstruction procedure on finite-element type meshes

An essentially non-oscillatory reconstruction for functions defined on finite-element type meshes was designed. Two related problems are studied: the interpolation of possibly unsmooth multivariate functions on arbitrary meshes and the reconstruction of a function from its average in the control volumes surrounding the nodes of the mesh. Concerning the first problem, we have studied the behavior of the highest coefficients of the Lagrange interpolation function which may admit discontinuities of locally regular curves. This enables us to choose the best stencil for the interpolation. The choice of the smallest possible number of stencils is addressed. Concerning the reconstruction problem, because of the very nature of the mesh, the only method that may work is the so called reconstruction via deconvolution method. Unfortunately, it is well suited only for regular meshes as we show, but we also show how to overcome this difficulty. The global method has the expected order of accuracy but is conservative up to a high order quadrature formula only. Some numerical examples are given which demonstrate the efficiency of the method.

Abgrall, R.↗

NASA NTRS1992

Design of an essentially non-oscillatory reconstruction procedure in finite-element type meshes

An essentially non oscillatory reconstruction for functions defined on finite element type meshes is designed. Two related problems are studied: the interpolation of possibly unsmooth multivariate functions on arbitary meshes and the reconstruction of a function from its averages in the control volumes surrounding the nodes of the mesh. Concerning the first problem, the behavior of the highest coefficients of two polynomial interpolations of a function that may admit discontinuities of locally regular curves is studied: the Lagrange interpolation and an approximation such that the mean of the polynomial on any control volume is equal to that of the function to be approximated. This enables the best stencil for the approximation to be chosen. The choice of the smallest possible number of stencils is addressed. Concerning the reconstruction problem, two methods were studied: one based on an adaptation of the so called reconstruction via deconvolution method to irregular meshes and one that lies on the approximation on the mean as defined above. The first method is conservative up to a quadrature formula and the second one is exactly conservative. The two methods have the expected order of accuracy, but the second one is much less expensive than the first one. Some numerical examples are given which demonstrate the efficiency of the reconstruction.

Abgrall, Remi↗

NASA NTRS2005

Lifetime Reliability Prediction of Ceramic Structures Under Transient Thermomechanical Loads

An analytical methodology is developed to predict the probability of survival (reliability) of ceramic components subjected to harsh thermomechanical loads that can vary with time (transient reliability analysis). This capability enables more accurate prediction of ceramic component integrity against fracture in situations such as turbine startup and shutdown, operational vibrations, atmospheric reentry, or other rapid heating or cooling situations (thermal shock). The transient reliability analysis methodology developed herein incorporates the following features: fast-fracture transient analysis (reliability analysis without slow crack growth, SCG); transient analysis with SCG (reliability analysis with time-dependent damage due to SCG); a computationally efficient algorithm to compute the reliability for components subjected to repeated transient loading (block loading); cyclic fatigue modeling using a combined SCG and Walker fatigue law; proof testing for transient loads; and Weibull and fatigue parameters that are allowed to vary with temperature or time. Component-to-component variation in strength (stochastic strength response) is accounted for with the Weibull distribution, and either the principle of independent action or the Batdorf theory is used to predict the effect of multiaxial stresses on reliability. The reliability analysis can be performed either as a function of the component surface (for surface-distributed flaws) or component volume (for volume-distributed flaws). The transient reliability analysis capability has been added to the NASA CARES/ Life (Ceramic Analysis and Reliability Evaluation of Structures/Life) code. CARES/Life was also updated to interface with commercially available finite element analysis software, such as ANSYS, when used to model the effects of transient load histories. Examples are provided to demonstrate the features of the methodology as implemented in the CARES/Life program.

Nemeth, Noel N.↗

NASA NTRS2017

Non-Axisymmetric Inflatable Pressure Structure (NAIPS) Full-Scale Pressure Test

Inflatable space structures have the potential to significantly reduce the required launch volume for large pressure vessels required for exploration applications including habitats, airlocks and tankage. In addition, mass savings can be achieved via the use of high specific strength softgoods materials, and the reduced design penalty from launching the structure in a densely packaged state. Large inclusions however, such as hatches, induce a high mass penalty at the interfaces with the softgoods and in the added rigid structure while reducing the packaging efficiency. A novel, Non-Axisymmetric Inflatable Pressure Structure (NAIPS) was designed and recently tested at NASA Langley Research Center to demonstrate an elongated inflatable architecture that could provide areas of low stress along a principal axis in the surface. These low stress zones will allow the integration of a flexible linear seal that substantially reduces the added mass and volume of a heritage rigid hatch structure. This paper describes the test of the first full-scale engineering demonstration unit (EDU) of the NAIPS geometry and a comparison of the results to finite element analysis.

Jones, Thomas C.↗

NASA NTRS1979

A program to compute three-dimensional subsonic unsteady aerodynamic characteristics using the doublet lattice method, L216 (DUBFLEX). Volume 2: Supplemental system design and maintenance document

The information necessary for execution of the digital computer program L216 on the CDC 6600 is described. L216 characteristics are based on the doublet lattice method. Arbitrary aerodynamic configurations may be represented with combinations of nonplanar lifting surfaces composed of finite constant pressure panel elements, and axially summetric slender bodies composed of constant pressure line elements. Program input consists of configuration geometry, aerodynamic parameters, and modal data; output includes element geometry, pressure difference distributions, integrated aerodynamic coefficients, stability derivatives, generalized aerodynamic forces, and aerodynamic influence coefficient matrices. Optionally, modal data may be input on magnetic field (tape or disk), and certain geometric and aerodynamic output may be saved for subsequent use.

Harrison, B. A.↗

NASA NTRS1996

Neutron diffraction measurements and modeling of residual strains in metal matrix composites

Neutron diffraction measurements at room temperature are used to characterize the residual strains in tungsten fiber-reinforced copper matrix, tungsten fiber-reinforced Kanthal matrix, and diamond particulate-reinforced copper matrix composites. Results of finite element modeling are compared with the neutron diffraction data. In tungsten/Kanthal composites, the fibers are in compression, the matrix is in tension, and the thermal residual strains are a strong function of the volume fraction of fibers. In copper matrix composites, the matrix is in tension and the stresses are independent of the volume fraction of tungsten fibers or diamond particles and the assumed stress free temperature because of the low yield strength of the matrix phase.

Saigal, A.↗

NASA NTRS1975

The NATA code; theory and analysis. Volume 2: User's manual

The NATA code is a computer program for calculating quasi-one-dimensional gas flow in axisymmetric nozzles and rectangular channels, primarily to describe conditions in electric archeated wind tunnels. The program provides solutions based on frozen chemistry, chemical equilibrium, and nonequilibrium flow with finite reaction rates. The shear and heat flux on the nozzle wall are calculated and boundary layer displacement effects on the inviscid flow are taken into account. The program contains compiled-in thermochemical, chemical kinetic and transport cross section data for high-temperature air, CO2-N2-Ar mixtures, helium, and argon. It calculates stagnation conditions on axisymmetric or two-dimensional models and conditions on the flat surface of a blunt wedge. Included in the report are: definitions of the inputs and outputs; precoded data on gas models, reactions, thermodynamic and transport properties of species, and nozzle geometries; explanations of diagnostic outputs and code abort conditions; test problems; and a user's manual for an auxiliary program (NOZFIT) used to set up analytical curvefits to nozzle profiles.

Bade, W. L.↗

NASA NTRS1969

ELAS: A general-purpose computer program for the equilibrium problems of linear structures. Volume 2: Documentation of the program

A general purpose digital computer program for the in-core solution of linear equilibrium problems of structural mechanics is documented. The program requires minimum input for the description of the problem. The solution is obtained by means of the displacement method and the finite element technique. Almost any geometry and structure may be handled because of the availability of linear, triangular, quadrilateral, tetrahedral, hexahedral, conical, triangular torus, and quadrilateral torus elements. The assumption of piecewise linear deflection distribution insures monotonic convergence of the deflections from the stiffer side with decreasing mesh size. The stresses are provided by the best-fit strain tensors in the least squares at the mesh points where the deflections are given. The selection of local coordinate systems whenever necessary is automatic. The core memory is used by means of dynamic memory allocation, an optional mesh-point relabelling scheme and imposition of the boundary conditions during the assembly time.

Utku, S.↗