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,261 records · Page 70

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Nonlinear mesomechanics of composites with periodic microstructure

This work is concerned with modeling the mechanical deformation or constitutive behavior of composites comprised of a periodic microstructure under small displacement conditions at elevated temperature. A mesomechanics approach is adopted which relates the microimechanical behavior of the heterogeneous composite with its in-service macroscopic behavior. Two different methods, one based on a Fourier series approach and the other on a Green's function approach, are used in modeling the micromechanical behavior of the composite material. Although the constitutive formulations are based on a micromechanical approach, it should be stressed that the resulting equations are volume averaged to produce overall effective constitutive relations which relate the bulk, volume averaged, stress increment to the bulk, volume averaged, strain increment. As such, they are macromodels which can be used directly in nonlinear finite element programs such as MARC, ANSYS and ABAQUS or in boundary element programs such as BEST3D. In developing the volume averaged or efective macromodels from the micromechanical models, both approaches will require the evaluation of volume integrals containing the spatially varying strain distributions throughout the composite material. By assuming that the strain distributions are spatially constant within each constituent phase-or within a given subvolume within each constituent phase-of the composite material, the volume integrals can be obtained in closed form. This simplified micromodel can then be volume averaged to obtain an effective macromodel suitable for use in the MARC, ANSYS and ABAQUS nonlinear finite element programs via user constitutive subroutines such as HYPELA and CMUSER. This effective macromodel can be used in a nonlinear finite element structural analysis to obtain the strain-temperature history at those points in the structure where thermomechanical cracking and damage are expected to occur, the so called damage critical points of the structure.

Walker, Kevin P.↗

NASA NTRS1993

Himalayan sandpiles

Himalayan avalanches are examined for agreement with power-law behavior between frequency of outfall and rock volume. The data support a scaling state beyond the size ranges in which corresponding sandpile experiments have suffered from finite-size effects. Large-scale data (12 orders of magnitude in rock volume) additionally suggest a possible role for the numerically predicted 'snowball effects' of large rocks triggering more net dislodging events.

Noever, David A.↗

NASA NTRS1991

Configuration optimization of space structures

The objective is to develop a computer aid for the conceptual/initial design of aerospace structures, allowing configurations and shape to be apriori design variables. The topics are presented in viewgraph form and include the following: Kikuchi's homogenization method; a classical shape design problem; homogenization method steps; a 3D mechanical component design example; forming a homogenized finite element; a 2D optimization problem; treatment of volume inequality constraint; algorithms for the volume inequality constraint; object function derivatives--taking advantage of design locality; stiffness variations; variations of potential; and schematics of the optimization problem.

Felippa, Carlos↗

NASA NTRS1994

Computer program for calculation of complex chemical equilibrium compositions and applications. Part 1: Analysis

This report presents the latest in a number of versions of chemical equilibrium and applications programs developed at the NASA Lewis Research Center over more than 40 years. These programs have changed over the years to include additional features and improved calculation techniques and to take advantage of constantly improving computer capabilities. The minimization-of-free-energy approach to chemical equilibrium calculations has been used in all versions of the program since 1967. The two principal purposes of this report are presented in two parts. The first purpose, which is accomplished here in part 1, is to present in detail a number of topics of general interest in complex equilibrium calculations. These topics include mathematical analyses and techniques for obtaining chemical equilibrium; formulas for obtaining thermodynamic and transport mixture properties and thermodynamic derivatives; criteria for inclusion of condensed phases; calculations at a triple point; inclusion of ionized species; and various applications, such as constant-pressure or constant-volume combustion, rocket performance based on either a finite- or infinite-chamber-area model, shock wave calculations, and Chapman-Jouguet detonations. The second purpose of this report, to facilitate the use of the computer code, is accomplished in part 2, entitled 'Users Manual and Program Description'. Various aspects of the computer code are discussed, and a number of examples are given to illustrate its versatility.

Gordon, Sanford↗

NASA NTRS2012

Biomedical Simulation Models of Human Auditory Processes

Detailed acoustic engineering models that explore noise propagation mechanisms associated with noise attenuation and transmission paths created when using hearing protectors such as earplugs and headsets in high noise environments. Biomedical finite element (FE) models are developed based on volume Computed Tomography scan data which provides explicit external ear, ear canal, middle ear ossicular bones and cochlea geometry. Results from these studies have enabled a greater understanding of hearing protector to flesh dynamics as well as prioritizing noise propagation mechanisms. Prioritization of noise mechanisms can form an essential framework for exploration of new design principles and methods in both earplug and earcup applications. These models are currently being used in development of a novel hearing protection evaluation system that can provide experimentally correlated psychoacoustic noise attenuation. Moreover, these FE models can be used to simulate the effects of blast related impulse noise on human auditory mechanisms and brain tissue.

Bicak, Mehmet M. A.↗

NASA NTRS2014

Three-Dimensional Imaging and Numerical Reconstruction of Graphite/Epoxy Composite Microstructure Based on Ultra-High Resolution X-Ray Computed Tomography

A combined experimental and computational study aimed at high-resolution 3D imaging, visualization, and numerical reconstruction of fiber-reinforced polymer microstructures at the fiber length scale is presented. To this end, a sample of graphite/epoxy composite was imaged at sub-micron resolution using a 3D X-ray computed tomography microscope. Next, a novel segmentation algorithm was developed, based on concepts adopted from computer vision and multi-target tracking, to detect and estimate, with high accuracy, the position of individual fibers in a volume of the imaged composite. In the current implementation, the segmentation algorithm was based on Global Nearest Neighbor data-association architecture, a Kalman filter estimator, and several novel algorithms for virtualfiber stitching, smoothing, and overlap removal. The segmentation algorithm was used on a sub-volume of the imaged composite, detecting 508 individual fibers. The segmentation data were qualitatively compared to the tomographic data, demonstrating high accuracy of the numerical reconstruction. Moreover, the data were used to quantify a) the relative distribution of individual-fiber cross sections within the imaged sub-volume, and b) the local fiber misorientation relative to the global fiber axis. Finally, the segmentation data were converted using commercially available finite element (FE) software to generate a detailed FE mesh of the composite volume. The methodology described herein demonstrates the feasibility of realizing an FE-based, virtual-testing framework for graphite/fiber composites at the constituent level.

Czabaj, M. W.↗

NASA NTRS1974

Cyclic fatigue analysis of rocket thrust chambers. Volume 1: OFHC copper chamber low cycle fatigue

A three-dimensional finite element elasto-plastic strain analysis was performed for the throat section of a regeneratively cooled rocket combustion chamber. The analysis employed the RETSCP finite element computer program. The analysis included thermal and pressure loads, and the effects of temperature dependent material properties, to determine the strain range corresponding to the chamber operating cycle. The analysis was performed for chamber configuration and operating conditions corresponding to a hydrogen-oxygen combustion chamber which was fatigue tested to failure. The computed strain range at typical chamber operating conditions was used in conjunction with oxygen-free, high-conductivity (OHFC) copper isothermal fatigue test data to predict chamber low-cycle fatigue life.

Miller, R. W.↗

NASA NTRS1974

Analytical and experimental investigation of a 1/8-scale dynamic model of the shuttle orbiter. Volume 3B: Supporting data

The NASA Structural Analysis System (NASTRAN) Model 1 finite element idealization, input data, and detailed analytical results are presented. The data presented include: substructuring analysis for normal modes, plots of member data, plots of symmetric free-free modes, plots of antisymmetric free-free modes, analysis of the wing, analysis of the cargo doors, analysis of the payload, and analysis of the orbiter.

Mason, P. W.↗

NASA NTRS1993

A study of fiber volume fraction effects in notched unidirectional SCS-6/Ti-15V-3Cr-3Al-3Sn composite

Notched unidirectional SCS-6/Ti-15-3 composite of three different fiber volume fractions (vf = 0.15, 0.37, and 0.41) was investigated for various room temperature microstructural and material properties including: fatigue crack initiation, fatigue crack growth, and fracture toughness. While the matrix hardness is similar for all fiber volume fractions, the fiber/matrix interfacial shear strength and matrix residual stress increases with fiber volume fraction. The composite fatigue crack initiation stress is shown to be matrix controlled and occurs when the net maximum matrix stress approaches the endurance limit stress of the matrix. A model is presented which includes residual stresses and presents the composite initiation stress as a function of fiber volume fraction. This model predicts a maximum composite initiation stress at vf approximately 0.15 which agrees with the experimental data. The applied composite stress levels were increased as necessary for continued crack growth. The applied Delta(K) values at crack arrest increase with fiber volume fraction by an amount better approximated using an energy based formulation rather than when scaled linear with modulus. After crack arrest, the crack growth rate exponents for vf37 and vf41 were much lower and toughness much higher, when compared to the unreinforced matrix, because of the bridged region which parades with the propagating fatigue crack. However, the vf15 material exhibited a higher crack growth rate exponent and lower toughness than the unreinforced matrix because once the bridged fibers nearest the crack mouth broke, the stress redistribution broke all bridged fibers, leaving an unbridged crack. Degraded, unbridged behavior is modeled using the residual stress state in the matrix ahead of the crack tip. Plastic zone sizes were directly measured using a metallographic technique and allow prediction of an effective matrix stress intensity which agrees with the fiber pressure model if residual stresses are considered. The sophisticated macro/micro finite element models of the 0.15 and 0.37 fiber volume fractions presented show good agreement with experimental data and the fiber pressure model when an estimated effective fiber/matrix debond length is used.

Covey, Steven J.↗

NASA NTRS1983

Volume integrals associated with the inhomegeneous Helmholtz equation. Part 2: Cylindrical region; rectangular region

Results are presented for volume integrals associated with the Helmholtz operator, nabla(2) + alpha(2), for the cases of a finite cylindrical region and a region of rectangular parallelepiped. By using appropriate Taylor series expansions and multinomial theorem, these volume integrals are obtained in series form for regions r r' and r 4', where r and r' are distances from the origin to the point of observation and source, respectively. When the wave number approaches zero, the results reduce directly to the potentials of variable densities.

Zhong, W. F.↗

NASA NTRS2002

The Speed of Axial Propagation of a Cylindrical Bubble Through a Cylindrical Vortex

Inspired by the rapid elongation of air columns injected into vortices by dolphins, we present an exact inviscid solution for the axial speed (assumed steady) of propagation of the tip of a semi-infinite cylindrical bubble along the axis of a cylindrical vortex. The bubble is assumed to be held at constant pressure by being connected to a reservoir, the lungs of the dolphin, say. For a given bubble pressure, there is a modest critical rotation rate above which steadily propagating bubbles exist. For a bubble at ambient pressure, the propagation speed of the bubble (relative to axial velocity within the vortex) varies between 0.5 and 0.6 of the maximum rotational speed of the vortex. Surprisingly, the bubble tip can propagate (almost as rapidly) even when the pressure minimum in the vortex core is greater than the bubble pressure; in this case, solutions exhibit a dimple on the nose of the bubble. A situation important for incipient vortex cavitation, and one which dolphins also demonstrate, is elongation of a free bubble, i.e., one whose internal pressure may vary. Under the assumption that the acceleration term is small (checked a posteriori), the steady solution is applied at each instant during the elongation. Three types of behavior are then possible depending on physical parameters and initial conditions: (A) Unabated elongation with slowly increasing bubble pressure, and nearly constant volume. Volume begins to decrease in the late stages. (B1) Elongation with decreasing bubble pressure. A limit point of the steady solution is encountered at a finite bubble length. (B2) Unabated elongation with decreasing bubble pressure and indefinite creation of volume. This is made possible by the existence of propagating solutions at bubble pressures below the minimum vortex pressure. As the bubble stretches, its radius initially decreases but then becomes constant; this is also observed in experiments on incipient vortex cavitation.

Shariff, Karim↗

NASA NTRS1975

Subgrid or Reynolds stress-modeling for three-dimensional turbulence computations

A review is given of recent advances in two distinct computational methods for evaluating turbulence fields, namely, statistical Reynolds stress modeling and turbulence simulation, where large eddies are followed in time. It is shown that evaluation of the mean Reynolds stresses, rather than use of a scalar eddy viscosity, permits an explanation of streamline curvature effects found in several experiments. Turbulence simulation, with a new volume averaging technique and third-order accurate finite-difference computing is shown to predict the decay of isotropic turbulence in incompressible flow with rather modest computer storage requirements, even at Reynolds numbers of aerodynamic interest.

Rubesin, M. W.↗

NASA NTRS2005

Prediction of Mechanical Properties of Polymers With Various Force Fields

The effect of force field type on the predicted elastic properties of a polyimide is examined using a multiscale modeling technique. Molecular Dynamics simulations are used to predict the atomic structure and elastic properties of the polymer by subjecting a representative volume element of the material to bulk and shear finite deformations. The elastic properties of the polyimide are determined using three force fields: AMBER, OPLS-AA, and MM3. The predicted values of Young s modulus and shear modulus of the polyimide are compared with experimental values. The results indicate that the mechanical properties of the polyimide predicted with the OPLS-AA force field most closely matched those from experiment. The results also indicate that while the complexity of the force field does not have a significant effect on the accuracy of predicted properties, small differences in the force constants and the functional form of individual terms in the force fields determine the accuracy of the force field in predicting the elastic properties of the polyimide.

Odegard, Gregory M.↗

NASA NTRS1977

DESAP 2: A structural design program with stress and buckling constraints. Volume 1: Theoretical and user's manual

DESAP 2 is described as a finite element program for computer-automated, minimum weight design of elastic structures with constraints on stresses (including local instability criteria) and buckling loads. No limits are placed on the number of load conditions for stress-constrained design, but only one of these load conditions can be chosen as the potential buckling load. A substantial portion of DESAP 2, particularly the analysis of the prebuckling state, is derived from the SOLID SAP finite element program. The stress-constrained design is based on the classical stress ratio method, which drives the design towards a fully stressed state. The constraints on the buckling load are handled by solving the appropriate optimality criterion by successive iterations. During each iteration, the element sizes determined by the stress ratio method are used as the minimum size constraints. The element subroutines are organized in a manner that permits the user to make additions and changes with a minimal programming effort. Consequently, DESAP 2 can readily be changed into a special-purpose program to handle the user's specific design requirements and failure criteria.

Kiusalaas, J.↗

NASA NTRS1982

Turbofan forced mixer-nozzle internal flowfield. Volume 3: A computer code for 3-D mixing in axisymmetric nozzles

A finite difference method is developed for making detailed predictions of three dimensional subsonic turbulent flow in turbofan lobe mixers. The governing equations are solved by a forward-marching solution procedure which corrects an inviscid potential flow solution for viscous and thermal effects, secondary flows, total pressure distortion and losses, internal flow blockage and pressure drop. Test calculations for a turbulent coaxial jet flow verify that the turbulence model performs satisfactorily for this relatively simple flow. Lobe mixer flows are presented for two geometries typical of current mixer design. These calculations included both hot and cold flow conditions, and both matched and mismatched Mach number and total pressure in the fan and turbine streams.

Kreskovsky, J. P.↗

NASA NTRS1993

Effective thermoelastic and thermal properties of unidirectional fiber-reinforced composites and their sensitivity coefficients

Three-dimensional finite element models are used to assess the accuracy of the thermoelastic and thermal properties of unidirectional fiber-reinforced composites predicted by six different micromechanical models. The six models are: simple mechanics of materials type equations, fiber substructuring model, vanishing fiber diameter model, self-consistent model, Mori-Tanaka model, and method of cells. In addition, the finite element models are used to assess the accuracy of derivatives of the effective properties, with respect to each of the constituent material properties and fiber-volume ratio, computed using the six micromechanical models. The predictions of the micromechanical and finite element models for four advanced composite material systems are compared with experimental data. The results obtained in the present study show that the predictions of the Mori-Tanaka model and the method of cells are closer to those of the finite element models than those of all the other micromechanical models.

Noor, Ahmed K.↗

NASA NTRS1992

Three-dimensional finite element simulation of intermingled-fiber hybrid composite behavior

Three-dimensional finite element methods and the intraply hybrid micromechanics equations are used to predict composite properties for a unidirectional graphite-epoxy primary composite with S-glass fibers used as hybridizing fibers. The micromechanics equations are embedded in a computer code ICAN (Integrated Composites Analyzer). The three-dimensional finite element model consists of three-by-three unit cell array, with a total fiber volume ratio of 0.54. There is a good agreement between the composite properties and microstresses obtained from both methods. The results indicate that the finite element methods and micromechanics equations can be used to obtain the properties of intermingled hybrid composites needed for analysis/design of hybrid composite structures.

Mital, Subodh K.↗

NASA NTRS1971

Transpiration and film cooling boundary layer computer program. Volume 1: Numerical solutions of the turbulent boundary layer equations with equilibrium chemistry

A finite difference turbulent boundary layer computer program has been developed. The program is primarily oriented towards the calculation of boundary layer performance losses in rocket engines; however, the solution is general, and has much broader applicability. The effects of transpiration and film cooling as well as the effect of equilibrium chemical reactions (currently restricted to the H2-O2 system) can be calculated. The turbulent transport terms are evaluated using the phenomenological mixing length - eddy viscosity concept. The equations of motion are solved using the Crank-Nicolson implicit finite difference technique. The analysis and computer program have been checked out by solving a series of both laminar and turbulent test cases and comparing the results to data or other solutions. These comparisons have shown that the program is capable of producing very satisfactory results for a wide range of flows. Further refinements to the analysis and program, especially as applied to film cooling solutions, would be aided by the acquisition of a firm data base.

Levine, J. N.↗