Engineering Papers⌕ Search

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.

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

At least 1,333 records · Page 74

Analysis of Composite Skin-Stiffener Debond Specimens Using Volume Elements and a Shell/3D Modeling Technique

The debonding of a skin/stringer specimen subjected to tension was studied using three-dimensional volume element modeling and computational fracture mechanics. Mixed mode strain energy release rates were calculated from finite element results using the virtual crack closure technique. The simulations revealed an increase in total energy release rate in the immediate vicinity of the free edges of the specimen. Correlation of the computed mixed-mode strain energy release rates along the delamination front contour with a two-dimensional mixed-mode interlaminar fracture criterion suggested that in spite of peak total energy release rates at the free edge the delamination would not advance at the edges first. The qualitative prediction of the shape of the delamination front was confirmed by X-ray photographs of a specimen taken during testing. The good correlation between prediction based on analysis and experiment demonstrated the efficiency of a mixed-mode failure analysis for the investigation of skin/stiffener separation due to delamination in the adherents. The application of a shell/3D modeling technique for the simulation of skin/stringer debond in a specimen subjected to three-point bending is also demonstrated. The global structure was modeled with shell elements. A local three-dimensional model, extending to about three specimen thicknesses on either side of the delamination front was used to capture the details of the damaged section. Computed total strain energy release rates and mixed-mode ratios obtained from shell/3D simulations were in good agreement with results obtained from full solid models. The good correlations of the results demonstrated the effectiveness of the shell/3D modeling technique for the investigation of skin/stiffener separation due to delamination in the adherents.

Krueger, Ronald↗

Constitutive Modeling of Piezoelectric Polymer Composites

A new modeling approach is proposed for predicting the bulk electromechanical properties of piezoelectric composites. The proposed model offers the same level of convenience as the well-known Mori-Tanaka method. In addition, it is shown to yield predicted properties that are, in most cases, more accurate or equally as accurate as the Mori-Tanaka scheme. In particular, the proposed method is used to determine the electromechanical properties of four piezoelectric polymer composite materials as a function of inclusion volume fraction. The predicted properties are compared to those calculated using the Mori-Tanaka and finite element methods.

Odegard, Gregory M.↗

Modification of codes NUALGAM and BREMRAD, Volume 1

The NUGAM2 code predicts forward and backward angular energy differential and integrated distributions for gamma photons and fluorescent radiation emerging from finite laminar transport media. It determines buildup and albedo data for scientific research and engineering purposes; it also predicts the emission characteristics of finite radioisotope sources. The results are shown to be in very good agreement with available published data. The code predicts data for many situations in which no published data is available in the energy range up to 5 MeV. The NUGAM3 code predicts the pulse height response of inorganic (NaI and CsI) scintillation detectors to gamma photons. Because it allows the scintillator to be clad and mounted on a photomultiplier as in the experimental or industrial application, it is a more practical and thus useful code than others previously reported. Results are in excellent agreement with published Monte Carlo and experimental data in the energy range up to 4.5 MeV.

Steyn, J. J.↗

An evaluation of a coupled microstructural approach for the analysis of functionally graded composites via the finite-element method

A comparison is presented between the predictions of the finite-element analysis and a recently developed higher-order theory for functionally graded materials subjected to a thorough-thickness temperature gradient. In contrast to existing micromechanical theories that utilize classical (i.e., uncoupled) homogenization schemes to calculate micro-level and macro-level stress and displacement fields in materials with uniform or nonuniform fiber spacing (i.e., functionally graded materials), the new theory explicitly couples the microstructural details with the macrostructure of the composite. Previous thermo-elastic analysis has demonstrated that such coupling is necessary when: the temperature gradient is large with respect to the dimension of the reinforcement; the characteristic dimension of the reinforcement is large relative to the global dimensions of the composite and the number of reinforcing fibers or inclusions is small. In these circumstances, the standard micromechanical analyses based on the concept of the representative volume element used to determine average composite properties produce questionable results. The comparison between the predictions of the finite-element method and the higher-order theory presented herein establish the theory's accuracy in predicting thermal and stress fields within composites with a finite number of fibers in the thickness direction subjected to a thorough-thickness thermal gradient.

Pindera, Marek-Jerzy↗

Probabilistic Structural Analysis Methods for select space propulsion system components (PSAM). Volume 3: Literature surveys and technical reports

The technical effort and computer code developed during the first year are summarized. Several formulations for Probabilistic Finite Element Analysis (PFEA) are described with emphasis on the selected formulation. The strategies being implemented in the first-version computer code to perform linear, elastic PFEA is described. The results of a series of select Space Shuttle Main Engine (SSME) component surveys are presented. These results identify the critical components and provide the information necessary for probabilistic structural analysis.

Source record↗

The 3.3K thrust chamber life prediction

The cause of low-cycle fatigue failure of the 3.3K Thrust Chamber was investigated. This thrust chamber typifies the current trend in rocket nozzle design which calls for high performance coupled with weight and volume limitations as well as the requirement of reusability. The analysis was performed with the BOPACE finite-element computer program which provides capability to determine viscoplastic response of a structure subjected to cyclic thermal and mechanical loading. Results are presented which show the critical region for low-cycle fatigue and the history of strain within that region for each thermo-mechanical loading cycle in the 3.3K thrust chamber. The predicted behavior was used to evaluate the low-cycle fatigue life near the throat plane of the chamber. The results show that BOPACE provides an extremely accurate prediction of structural behavior; the critical region was identified and the life determined from computed strains was within 154 cycles of the observed failure at 1013 cycles.

Armstrong, W. H.↗

A stratospheric chemical instability

The equations which determine partitioning of Cl(x) in steady state have multiple (three) solutions under conditions which might arise in the high-latitude winter stratosphere. Two of these solutions are stable, one is unstable, to infinitesimal perturbations. The relative stability of solutions is examined by subjecting the system to finite perturbations. The more stable solution is found to eliminate the less stable when semi-infinite volumes of the two solutions are placed in contact. The high-ClO, low NO2 solution is more stable under most conditions. Transitions from less to more stable states are slow in winter but may occur more rapidly when the seasonal variation of insolation is taken into account.

Fox, J. L.↗

Evaluation of 3 numerical methods for propulsion integration studies on transonic transport configurations

An effort has been undertaken at the NASA Langley Research Center to assess the capabilities of available computational methods for use in propulsion integration design studies of transonic transport aircraft, particularly of pylon/nacelle combinations which exhibit essentially no interference drag. The three computer codes selected represent state-of-the-art computational methods for analyzing complex configurations at subsonic and transonic flight conditions. These are: EULER, a finitie volume solution of the Euler equation; VSAERO, a panel solution of the Laplace equation; and PPW, a finite difference solution of the small disturbance transonic equations. In general, all three codes have certain capabilities that allow them to be of some value in predicting the flows about transport configurations, but all have limitations. Until more accurate methods are available, careful application and interpretation of the results of these codes are needed.

Yaros, S. F.↗

Integrated support structure for GASCAN 2

The focus of the Worcester Polytechnic Institute (WPI) Advanced Space Design Program was the preliminary design of the Integrated Support Structure for GASCAN II, a Get Away Special canister donated by the MITRE Corporation. Two teams of three students each worked on the support structure. There was a structural design team and a thermal design team. The structure will carry three experiments also undergoing preliminary design this year, the mu-gravity Ignition Experiment, the Rotational Flow in Low Gravity Experiment, and the Ionospheric Properties and Propagation Experiment. The structural design team was responsible for the layout of the GASCAN and the preliminary design of the structure itself. They produced the physical interface specifications defining the baseline weights and volumes for the equipment and produced layout drawings of the system. The team produced static and modal finite element analysis of the structure using ANSYS. The thermal design team was responsible for the power and timing requirements of the payload and for the identification and preliminary analysis of potential thermal problems. The team produced the power, timing, and energy interface specifications and assisted in the development of the specification of the battery pack. The thermal parameters of each experiment were cataloged and the experiments were subjected to worst case heat transfer scenarios.

Source record↗

A geometric representation scheme suitable for shape optimization

A geometric representation scheme is outlined which utilizes the natural design variable concept. A base configuration with distinct topological features is created. This configuration is then deformed to define components with similar topology but different geometry. The values of the deforming loads are the geometric entities used in the shape representation. The representation can be used for all geometric design studies; it is demonstrated here for structural optimization. This technique can be used in parametric design studies, where the system response is defined as functions of geometric entities. It can also be used in shape optimization, where the geometric entities of an original design are modified to maximize performance and satisfy constraints. Two example problems are provided. A cantilever beam is elongated to meet new design specifications and then optimized to reduce volume and satisfy stress constraints. A similar optimization problem is presented for an automobile crankshaft section. The finite element method is used to perform the analyses.

Tortorelli, Daniel A.↗

Bistable Deployable Composite Booms With Parabolic Cross-Sections

This paper investigates how stable equilibrium states in the extended and coiled configurations can be predicted in thin-shelled composite booms with parabolic cross-sections. These conic shapes potentially offer greater stiffness properties when compared to circular cross-sections, which is critical for improving the load bearing performance of deployed booms. Inducing bistability through the choice of composite shell layups in parabolic booms would allow for controllable self-deployment due to a less energetic coiled state when compared to monostable booms. An inextensional analytical model is used to predict the stable coiled diameters of tape spring and Collapsible Tubular Mast (CTM) booms with parabolic cross-sections. The parabolic section is discretized into circular segments using biarc spline interpolation, which allows them to be integrated into the strain energy minimization procedure used to obtain the equilibrium states. When the parabolic booms are parametrically compared against circular booms with identical layups, flattened height, and mass, the former are found to generally have better stiffness performance while being less efficient in stowed volume as evidenced by larger coiled diameters. Analytical coiled diameters and their strain energy are verified with finite element simulations for optimal parabolic tape spring and CTM booms.

Andrew J. Lee↗

Finite Element Analysis of a Deployable Space Structure

The dynamic characteristics of the Solar Array Flight Experiment (SAFE) structure during deployment and retraction are investigated. The SAFE structure consists of a deployable mast with an attached solar blanket designed with accordion type folds to permit packaging in a small volume. The planar form of the blanket geometry during deployment is maintained by a blanket tension/guidewire system. Structurally, the mast is modeled as an Euler beam column with inplane and out of plane bending and finite torsional stiffness. For out of plane motion, the blanket is modeled as a distributed mass uniformly supported by the three guidewires. For inplane motion the blanket displacements are assumed to vary linearly from the mast base to the mast tip. The mathematical model uses a virtual work formulation, required because the axial loading on the mast is nonconservative, combined with assumed beam modes to derive the differential equations of motion. Consideration of the time dependent boundary conditions results in an infinite set of ODE with time dependent coefficients. Finally, correlation of mast tip accelerations to mast base bending moments for specified modal motions are indicated.

Weeks, G. E.↗

Coupled Structural, Thermal, Phase-change and Electromagnetic Analysis for Superconductors, Volume 2

Described are 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 electromag 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 thermel and mechanical effects, and (4) computer implementation and solution of the superconductivity transition problem. 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 theoretical development is described in two volumes. Volume 1 describes mostly formulation specific problems. Volume 2 describes generalization of those formulations.

Felippa, C. A.↗

Dynamic characteristics of a hydrostatic gas bearing driven by oscillating exhaust pressure

Vibration of a statically loaded, inherently compensated hydrostatic journal bearing due to oscillating exhaust pressure is investigated. Both angular and radial vibration modes are analyzed. The time-dependent Reynolds equation governing the pressure distribution between the oscillating journal and sleeve is solved together with the journal equation of motion to obtain the response characteristics of the bearing. The Reynolds equation and the equation of motion are simplified by applying regular perturbation theory for small displacements. The numerical solutions of the perturbation equations are obtained by discretizing the pressure field using finite-difference aproximations with a discrete, nonuniform line-source model which excludes effects due to feeding hole volume. An iterative scheme is used to simultaneously satisfy the equations of motion for the journal. The results presented include Bode plots of bearing-oscillation gain and phase for a particular bearing configuration for various combinations of parameters over a range of frequencies, including the resonant frequency.

Watkins, C. B.↗

Simulation of Two-Fluid Flows by the Least-Squares Finite Element Method Using a Continuum Surface Tension Model

In this paper a numerical procedure for simulating two-fluid flows is presented. This procedure is based on the Volume of Fluid (VOF) method proposed by Hirt and Nichols and the continuum surface force (CSF) model developed by Brackbill, et al. In the VOF method fluids of different properties are identified through the use of a continuous field variable (color function). The color function assigns a unique constant (color) to each fluid. The interfaces between different fluids are distinct due to sharp gradients of the color function. The evolution of the interfaces is captured by solving the convective equation of the color function. The CSF model is used as a means to treat surface tension effect at the interfaces. Here a modified version of the CSF model, proposed by Jacqmin, is used to calculate the tension force. In the modified version, the force term is obtained by calculating the divergence of a stress tensor defined by the gradient of the color function. In its analytical form, this stress formulation is equivalent to the original CSF model. Numerically, however, the use of the stress formulation has some advantages over the original CSF model, as it bypasses the difficulty in approximating the curvatures of the interfaces. The least-squares finite element method (LSFEM) is used to discretize the governing equation systems. The LSFEM has proven to be effective in solving incompressible Navier-Stokes equations and pure convection equations, making it an ideal candidate for the present applications. The LSFEM handles all the equations in a unified manner without any additional special treatment such as upwinding or artificial dissipation. Various bench mark tests have been carried out for both two dimensional planar and axisymmetric flows, including a dam breaking, oscillating and stationary bubbles and a conical liquid sheet in a pressure swirl atomizer.

Wu, Jie↗

Rapid Prototyping Integrated With Nondestructive Evaluation and Finite Element Analysis

Most reverse engineering approaches involve imaging or digitizing an object then creating a computerized reconstruction that can be integrated, in three dimensions, into a particular design environment. Rapid prototyping (RP) refers to the practical ability to build high-quality physical prototypes directly from computer aided design (CAD) files. Using rapid prototyping, full-scale models or patterns can be built using a variety of materials in a fraction of the time required by more traditional prototyping techniques (refs. 1 and 2). Many software packages have been developed and are being designed to tackle the reverse engineering and rapid prototyping issues just mentioned. For example, image processing and three-dimensional reconstruction visualization software such as Velocity2 (ref. 3) are being used to carry out the construction process of three-dimensional volume models and the subsequent generation of a stereolithography file that is suitable for CAD applications. Producing three-dimensional models of objects from computed tomography (CT) scans is becoming a valuable nondestructive evaluation methodology (ref. 4). Real components can be rendered and subjected to temperature and stress tests using structural engineering software codes. For this to be achieved, accurate high-resolution images have to be obtained via CT scans and then processed, converted into a traditional file format, and translated into finite element models. Prototyping a three-dimensional volume of a composite structure by reading in a series of two-dimensional images generated via CT and by using and integrating commercial software (e.g. Velocity2, MSC/PATRAN (ref. 5), and Hypermesh (ref. 6)) is being applied successfully at the NASA Glenn Research Center. The building process from structural modeling to the analysis level is outlined in reference 7. Subsequently, a stress analysis of a composite cooling panel under combined thermomechanical loading conditions was performed to validate this process.

Abdul-Aziz, Ali↗

Electromagnetic scattering from a layer of finite length, randomly oriented, dielectric, circular cylinders over a rough interface with application to vegetation

A scattering model for defoliated vegetation is developed by treating a layer of defoliated vegetation as a collection of randomly oriented dielectric cylinders of finite length over an irregular ground surface. Both polarized and depolarized backscattering are computed and their behavior versus the volume fraction, the incidence angle, the frequency, the angular distribution and the cylinder size are illustrated. It is found that both the angular distribution and the cylinder size have significant effects on the backscattered signal. The present theory is compared with measurements from defoliated vegetations.

Karam, M. A.↗

Hollow conical jet models for SS 433 - A paradigm lost?

A precessing jet such as that in SS 433 may be approximated as an axisymmetric flow, if the precession time is short by comparison to the propagation time. A series of simulations has been conducted for precessing jets using an axisymmetric finite-difference hydrodynamics code. Examinations are made of hollow cylindrical jets, which lack the complication of a growing interior volume, and conical jets, which model the behavior of a precessing jet propagating on the surface of its precession cone.

Kochanek, Christopher S.↗