Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element analysis”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Finite-element methodology for thermal analysis of convectively cooled structures

A finite-element method for steady-state thermal analysis of convectively cooled structures is presented. The method is based on representing the coolant passages by finite elements with fluid bulk temperature nodes and fluid/structure interface modes. Four finite elements are described: two basic elements (mass transport and surface convection) and two integrated elements for applications to discrete tube and plate-fin convectively cooled structural configurations. Comparative finite-element and lumped-parameter thermal analyses of several convectively cooled structures demonstrate the practicality of utilizing finite-element methodology for thermal analysis of realistic convectively cooled structures.

Source record↗
NASA NTRS1990

Development of an integrated BEM approach for hot fluid structure interaction

A comprehensive boundary element method is presented for transient thermoelastic analysis of hot section Earth-to-Orbit engine components. This time-domain formulation requires discretization of only the surface of the component, and thus provides an attractive alternative to finite element analysis for this class of problems. In addition, steep thermal gradients, which often occur near the surface, can be captured more readily since with a boundary element approach there are no shape functions to constrain the solution in the direction normal to the surface. For example, the circular disc analysis indicates the high level of accuracy that can be obtained. In fact, on the basis of reduced modeling effort and improved accuracy, it appears that the present boundary element method should be the preferred approach for general problems of transient thermoelasticity.

Dargush, G. F.↗
DOE OSTI2023

Fluid structural interaction simulation of the MP-1 plate performance irradiated in the advanced test reactor

Due to the high neutron fluxes they generate and increased heat transfer performance, plate type fuels are used in the U.S. high-performance research reactors. During irradiation, a significant amount of fission energy (i.e., ~200 MeV per atom) is released by the U-235 chain reaction that is carried away by the coolant. Through thermal hydraulic analysis, the coolant’s heat transfer capability is investigated to ensure that the desired fuel temperature can be maintained. At high temperatures, the plate undergoes elastic/plastic deformation, creep, and swelling as a result of both the temperature gradients and the fission gas production within the fuel. These effects are studied via fuel performance analysis. When the plate deformation is small enough that the coolant flow remains relatively unchanged, conducting these two types of analyses independently will suffice. But at high fuel burnups, the swollen plates may encroach into the coolant channels that separate the fuel plates from each other and cause these channels to narrow. Large reductions in channel gap size imperil cooling performance, causing fuel temperatures to rise. Further, if the plate deformation is asymmetric, the fuel centerline will shift toward one side of the channel, causing an uneven reduction in coolability. In addition, a boehmite (oxide) layer will, over time, grow on the plate surface, further obstructing the heat transfer process. To precisely predict fuel temperatures/deformation, a complete coupled analysis that considers coolant flow, heat transfer, oxide growth, elastic/plastic deformation, creep, and swelling is needed; however, this type of analysis method is not available in the literature. To fill the gap, this research developed a fluid structure interaction (FSI) approach to the fuel plate analysis, then successfully applied it to the Mini-Plate (MP)-1 experiment, which was irradiated in the Advanced Test Reactor (ATR) for both one and two cycles. The complete analysis coupled STAR-CCM+, a computational fluid dynamic (CFD) software for calculating flow, with Abaqus, a finite element analysis code for calculating plate deformation. Improvements in the results were found when comparing the fully coupled analysis to the independently conducted analyses but they were not significant due to the miniature size of the plates and the relatively short irradiation time. In the future, the fully coupled approach presented herein will be applied to full-size fuel plates with longer irradiation cycles once additional experiments become available.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Modeling the hydrodynamic impact on the tool influence function during hemispherical subaperture optical polishing

To fabricate high-precision and accurate optics relative to the optical design surface, a high level of deterministic control of material removal (i.e., the tool influence function, TIF) during subaperture tool polishing is required. In this study, a detailed analysis of the pressure distribution, which is a key component of the TIF, has been performed using finite element analysis to couple together solid mechanics and fluid dynamics. Modeling experimental parameters of recently published work reveals that, when considering tool deformation, which in turn influences the fluid film thickness between the tool and workpiece, the effective pressure profile has a flat-top distribution. This flat-top pressure profile differs from the parabolic pressure distributions predicted by Hertzian mechanics. Furthermore, the shear contribution is shown here to be a key contributor to material removal, inducing the removal at the periphery of the contact edge and even outside the generally accepted contact area. Finally, the simulated fluid velocities provide evidence of mixed-mode contact polishing, supporting recent experimental findings that also suggest that onset of hydroplaning contributions lead to material removal drop-off.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Technology Demonstration of a High-Pressure Swirl Oxy-Coal Combustor

This technical report presents the exploration of the design and prototyping of a High-Pressure Swirl Oxy-coal Combustor. Pressurized oxy-coal combustion systems have the potential to improve efficiency along with an increased carbon capture rate. Reduction of flue gas at higher pressure, smaller system size, and capital cost reductions render high-pressure oxy-coal systems particularly attractive as next-generation energy-producing systems. High-pressure oxy-coal combustion systems are a recent concept, and thus operability issues of combustor designs for such systems are not fully understood. Significant challenges exist to maintain oxy-coal combustion stability at elevated pressure and a high CO 2 diluent environment. Although a body of knowledge exists for high-pressure oxygen combustion in rocket engines (or similar applications), it is yet to be strategized how these fundamental concepts can be translated to low-temperature CO 2 diluent combustion regimes. The realization of the pressurized oxy-coal based systems requires combustor components to be designed and demonstrated for an operating pressure over 10 bar. However, pressurized oxy-coal combustor design information at this pressure range and scale relevant to validate those proposed systems is currently limited. Experimental data from MWth scale oxy-coal combustors are needed to identify the optimal trade-off between net efficiency and systems size. The proposed effort is aimed at demonstrating a 1 MWth down-fired swirl Oxy-Coal combustor and investigate the interrelation between combustor operating conditions (pressure; flame stability; flue gas recirculation ratio) and conversion efficiencies to minimize oxygen requirements. One of the key challenges is to configure burner design (i.e., swirl number and injector) and operating conditions for high-pressure oxy-coal combustion systems. These experiments differ from current systems partly due to the high theoretical flame temperature and related burner operability issues associated with oxy-combustion. An ASPEN PLUS® model study for 550 MWe TIPS and ENEL pressurized oxy-coal systems with CO 2 recirculation was performed to evaluate system design, subsystems sizing, and operating condition determination. The system analysis effort included TRL and technology gap determination of subsystems and critical components. This information was scaled to develop design requirements (design pressure and flue gas recirculation: RR Flue Gas = $\frac{m_{flue}}{m_{total}}$) for the 1 MWth combustor. The effects of a wide range of carbon dioxide recirculation ratios on the thermal efficiency of ENEL and TIPS cycles are studied. The pressure of 10 bar and 80 bar are used for ENEL and TIPS cycles, respectively. The thermal efficiency of ENEL is significantly higher than the efficiency of TIPS at a pressure of less than 10 bar. The insights from system analysis were then used to design a 1 MWth swirl oxy-coal combustor. Flame temperature analysis and material strength analysis was performed to determine the combustor thickness. The structural integrity of the combustor was validated by finite element analysis using Abacus® and Hypermesh®. Feasibility of igniters and secondary burners are investigated in successful high-pressure oxy-methane combustion. The secondary burners are designed in such a way that it can operate between 100 to 500 kW firing input. Three generations of the pintle injector were designed based on swirl numbers (S=0, 0.9, and 1.2). Key pintle injector parameters such as pintle size, pintle orifice size, spray pattern were investigated by cold flow tests. Information from these tests was used to modify injector design for smooth and successful operation. A 5 mm pintle orifice size was decided upon as the optimum size for oxy-coal operation for the combustor. Shadow sizing experiments were performed to identify the atomization rate of each injector. Different coal water slurry mixtures (30 – 50% coal by wt% in the mixture) at various total momentum ratios (TMR) were investigated for this purpose. These experiments provided decisive information to choose the best design of the injector. The injector with 1.2 swirl provided higher atomization in all cases than other designs. The mean equivalent droplet size of the jet was similar at different TMR and mixture ratios using this injector, thus making it suitable for use in most cases. Therefore, the 1.2 swirl-pintle injector was chosen for the shakedown test. The combustor and other sub-systems, including feed systems and control and data acquisition, have been manufactured, assembled, and integrated. The total system integration and installation began on July 1, 2020. The shake-down tests and initial operational capability demonstration are expected to be completed by September 30, 2020.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS1974

Specialists Meeting on Rotorcraft Dynamics, Moffett Field, Calif., February 13-15, 1974, Proceedings

Analysis of specific problems in rotorcraft dynamics. Topics include hingeless rotor theory, dynamic stall modelling, periodic systems identification, analysis of complex systems with phasing matrices, flapping stability, flap-lag dynamics at high advance ratios, finite element analysis and fuselage free-vibration characteristics, coupled rotor/frame vibration methods, gust response characteristics with unsteady stall effects, antiresonance theory, cyclic feathering motions and dynamic loads, control load envelope shaping, rotor aeroelasticity, use of Floquet theory, theory of proprotors and tilt-rotors, two-bladed teetering rotors, stability of air and ground resonance, vertical-plane pendulum absorbers, multicyclic jet-flap control, engine/frame interface dynamics, and others. The minutes of the question and answer periods following the presentations are presented in the supplement. Individual items are announced in this issue.

Source record↗
NASA NTRS1991

Inelastic Strain Analysis of Solder Joint in NASA Fatigue Specimen

The solder fatigue specimen designed by NASA-GSFC/UNISYS is analyzed in order to obtain the inelastic strain history during two different representative temperature cycles specified by UNISYS. In previous reports (dated July 25, 1990, and November 15, 1990), results were presented of the elastic-plastic and creep analysis for delta T = 31 C cycle, respectively. Subsequent results obtained during the current phase, from viscoplastic finite element analysis of the solder fatigue specimen for delta T = 113 C cycle are summarized. Some common information is repeated for self-completeness. Large-deformation continuum formulations in conjunction with a standard linear solid model is utilized for modeling the solder constitutive creep-plasticity behavior. Relevant material properties are obtained from the literature. Strain amplitudes, mean strains, and residual strains (as well as stresses) accumulated due to a representative complete temperature cycle are obtained as a result of this analysis. The partitioning between elastic strains, time-independent inelastic (plastic) strains, and time-dependent inelastic (creep) strains is also explicitly obtained for two representative cycles. Detailed plots are presented for two representative temperature cycles. This information forms an important input for fatigue damage models, when predicting the fatigue life of solder joints under thermal cycling

Dasgupta, Abhijit↗
NASA NTRS1998

Aircraft Structural Mass Property Prediction Using Conceptual-Level Structural Analysis

This paper describes a methodology that extends the use of the Equivalent LAminated Plate Solution (ELAPS) structural analysis code from conceptual-level aircraft structural analysis to conceptual-level aircraft mass property analysis. Mass property analysis in aircraft structures has historically depended upon parametric weight equations at the conceptual design level and Finite Element Analysis (FEA) at the detailed design level ELAPS allows for the modeling of detailed geometry, metallic and composite materials, and non-structural mass coupled with analytical structural sizing to produce high-fidelity mass property analyses representing fully configured vehicles early in the design process. This capability is especially valuable for unusual configuration and advanced concept development where existing parametric weight equations are inapplicable and FEA is too time consuming for conceptual design. This paper contrasts the use of ELAPS relative to empirical weight equations and FEA. ELAPS modeling techniques are described and the ELAPS-based mass property analysis process is detailed Examples of mass property stochastic calculations produced during a recent systems study are provided This study involved the analysis of three remotely piloted aircraft required to carry scientific payloads to very high altitudes at subsonic speeds. Due to the extreme nature of this high-altitude flight regime,few existing vehicle designs are available for use in performance and weight prediction. ELAPS was employed within a concurrent engineering analysis process that simultaneously produces aerodynamic, structural, and static aeroelastic results for input to aircraft performance analyses. The ELAPS models produced for each concept were also used to provide stochastic analyses of wing structural mass properties. The results of this effort indicate that ELAPS is an efficient means to conduct multidisciplinary trade studies at the conceptual design level.

Sexstone, Matthew G.↗
NASA NTRS2001

Notch Sensitivity of Woven Ceramic Matrix Composites Under Tensile Loading: An Experimental, Analytical, and Finite Element Study

The stress concentrations associated with circular notches and subjected to uniform tensile loading in woven ceramic matrix composites (CMCs) have been investigated for high-efficient turbine engine applications. The CMC's were composed of Nicalon silicon carbide woven fabric in SiNC matrix manufactured through polymer impregnation process (PIP). Several combinations of hole diameter/plate width ratios and ply orientations were considered in this study. In the first part, the stress concentrations were calculated measuring strain distributions surrounding the hole using strain gages at different locations of the specimens during the initial portion of the stress-strain curve before any microdamage developed. The stress concentration was also calculated analytically using Lekhnitskii's solution for orthotropic plates. A finite-width correction factor for anisotropic and orthotropic composite plate was considered. The stress distributions surrounding the circular hole of a CMC's plate were further studied using finite element analysis. Both solid and shell elements were considered. The experimental results were compared with both the analytical and finite element solutions. Extensive optical and scanning electron microscopic examinations were carried out for identifying the fracture behavior and failure mechanisms of both the notched and notched specimens. The stress concentration factors (SCF) determined by analytical method overpredicted the experimental results. But the numerical solution underpredicted the experimental SCF. Stress concentration factors are shown to increase with enlarged hole size and the effects of ply orientations on stress concentration factors are observed to be negligible. In all the cases, the crack initiated at the notch edge and propagated along the width towards the edge of the specimens.

Haque, A.↗
NASA NTRS2002

Application of Interface Technology in Progressive Failure Analysis of Composite Panels

A progressive failure analysis capability using interface technology is presented. The capability has been implemented in the COMET-AR finite element analysis code developed at the NASA Langley Research Center and is demonstrated on composite panels. The composite panels are analyzed for damage initiation and propagation from initial loading to final failure using a progressive failure analysis capability that includes both geometric and material nonlinearities. Progressive failure analyses are performed on conventional models and interface technology models of the composite panels. Analytical results and the computational effort of the analyses are compared for the conventional models and interface technology models. The analytical results predicted with the interface technology models are in good correlation with the analytical results using the conventional models, while significantly reducing the computational effort.

Sleight, D. W.↗
NASA NTRS2003

CARES/LIFE Ceramics Analysis and Reliability Evaluation of Structures Life Prediction Program

This manual describes the Ceramics Analysis and Reliability Evaluation of Structures Life Prediction (CARES/LIFE) computer program. The program calculates the time-dependent reliability of monolithic ceramic components subjected to thermomechanical and/or proof test loading. CARES/LIFE is an extension of the CARES (Ceramic Analysis and Reliability Evaluation of Structures) computer program. The program uses results from MSC/NASTRAN, ABAQUS, and ANSYS finite element analysis programs to evaluate component reliability due to inherent surface and/or volume type flaws. CARES/LIFE accounts for the phenomenon of subcritical crack growth (SCG) by utilizing the power law, Paris law, or Walker law. The two-parameter Weibull cumulative distribution function is used to characterize the variation in component strength. The effects of multiaxial stresses are modeled by using either the principle of independent action (PIA), the Weibull normal stress averaging method (NSA), or the Batdorf theory. Inert strength and fatigue parameters are estimated from rupture strength data of naturally flawed specimens loaded in static, dynamic, or cyclic fatigue. The probabilistic time-dependent theories used in CARES/LIFE, along with the input and output for CARES/LIFE, are described. Example problems to demonstrate various features of the program are also included.

Nemeth, Noel N.↗
NASA NTRS1998

Aircraft Structural Mass Property Prediction Using Conceptual-Level Structural Analysis

This paper describes a methodology that extends the use of the Equivalent LAminated Plate Solution (ELAPS) structural analysis code from conceptual-level aircraft structural analysis to conceptual-level aircraft mass property analysis. Mass property analysis in aircraft structures has historically depended upon parametric weight equations at the conceptual design level and Finite Element Analysis (FEA) at the detailed design level. ELAPS allows for the modeling of detailed geometry, metallic and composite materials, and non-structural mass coupled with analytical structural sizing to produce high-fidelity mass property analyses representing fully configured vehicles early in the design process. This capability is especially valuable for unusual configuration and advanced concept development where existing parametric weight equations are inapplicable and FEA is too time consuming for conceptual design. This paper contrasts the use of ELAPS relative to empirical weight equations and FEA. ELAPS modeling techniques are described and the ELAPS-based mass property analysis process is detailed. Examples of mass property stochastic calculations produced during a recent systems study are provided. This study involved the analysis of three remotely piloted aircraft required to carry scientific payloads to very high altitudes at subsonic speeds. Due to the extreme nature of this high-altitude flight regime, few existing vehicle designs are available for use in performance and weight prediction. ELAPS was employed within a concurrent engineering analysis process that simultaneously produces aerodynamic, structural, and static aeroelastic results for input to aircraft performance analyses. The ELAPS models produced for each concept were also used to provide stochastic analyses of wing structural mass properties. The results of this effort indicate that ELAPS is an efficient means to conduct multidisciplinary trade studies at the conceptual design level.

Sexstone, Matthew G.↗
NASA NTRS2006

Three-Dimensional Field Solutions for Multi-Pole Cylindrical Halbach Arrays in an Axial Orientation

This article presents three-dimensional B field solutions for the cylindrical Halbach array in an axial orientation. This arrangement has applications in the design of axial motors and passive axial magnetic bearings and couplers. The analytical model described here assumes ideal magnets with fixed and uniform magnetization. The field component functions are expressed as sums of 2-D definite integrals that are easily computed by a number of mathematical analysis software packages. The analysis is verified with sample calculations and the results are compared to equivalent results from traditional finite-element analysis (FEA). The field solutions are then approximated for use in flux linkage and induced EMF calculations in nearby stator windings by expressing the field variance with angular displacement as pure sinusoidal function whose amplitude depends on radial and axial position. The primary advantage of numerical implementation of the analytical approach presented in the article is that it lends itself more readily to parametric analysis and design tradeoffs than traditional FEA models.

Thompson, William K.↗
NASA NTRS2009

Failure Analysis of Sapphire Refractive Secondary Concentrators

Failure analysis was performed on two sapphire, refractive secondary concentrators (RSC) that failed during elevated temperature testing. Both concentrators failed from machining/handling damage on the lens face. The first concentrator, which failed during testing to 1300 C, exhibited a large r-plane twin extending from the lens through much of the cone. The second concentrator, which was an attempt to reduce temperature gradients and failed during testing to 649 C, exhibited a few small twins on the lens face. The twins were not located at the origin, but represent another mode of failure that needs to be considered in the design of sapphire components. In order to estimate the fracture stress from fractographic evidence, branching constants were measured on sapphire strength specimens. The fractographic analysis indicated radial tensile stresses of 44 to 65 MPa on the lens faces near the origins. Finite element analysis indicated similar stresses for the first RSC, but lower stresses for the second RSC. Better machining and handling might have prevented the fractures, however, temperature gradients and resultant thermal stresses need to be reduced to prevent twinning.

Salem, Jonathan A.↗
NASA NTRS2014

The Vibration Ring

The vibration ring was conceived as a driveline damping device to prevent structure-borne noise in machines. It has the appearance of a metal ring, and can be installed between any two driveline components like an ordinary mechanical spacer. Damping is achieved using a ring-shaped piezoelectric stack that is poled in the axial direction and connected to an electrical shunt circuit. Surrounding the stack is a metal structure, called the compression cage, which squeezes the stack along its poled axis when excited by radial driveline forces. The stack in turn generates electrical energy, which is either dissipated or harvested using the shunt circuit. Removing energy from the system creates a net damping effect. The vibration ring is much stiffer than traditional damping devices, which allows it to be used in a driveline without disrupting normal operation. In phase 1 of this NASA Seedling Fund project, a combination of design and analysis was used to examine the feasibility of this concept. Several designs were evaluated using solid modeling, finite element analysis, and by creating prototype hardware. Then an analytical model representing the coupled electromechanical response was formulated in closed form. The model was exercised parametrically to examine the stiffness and loss factor spectra of the vibration ring, as well as simulate its damping effect in the context of a simplified driveline model. The results of this work showed that this is a viable mechanism for driveline damping, and provided several lessons for continued development.

Asnani, Vivake M.↗
NASA NTRS1995

Effect of various approximations on predicted progressive failure in plain weave composites

Three-dimensional finite element analysis was used to simulate progressive failure of a plain weave composite subjected to in-plane extension. The loading was parallel to one of the tow directions. The effects of various characteristics of the finite element model on predicted behavior were examined. The predicted behavior was found to be sensitive to quadrature order, mesh refinement, and the material degradation model. Also the sensitivity of the predictions to the tow waviness was studied. The predicted strength decreased considerably with increased waviness. More numerical studies and comparisons with experimental data are needed to establish reliable guidelines for accurate progressive failure prediction.

Whitcomb, John↗
NASA NTRS2001

Effect of Speed (Centrifugal Load) on Gear Crack Propagation Direction

The effect of rotational speed (centrifugal force) on gear crack propagation direction was explored. Gears were analyzed using finite element analysis and linear elastic fracture mechanics. The analysis was validated with crack propagation experiments performed in a spur gear fatigue rig. The effects of speed, rim thickness, and initial crack location on gear crack propagation direction were investigated. Crack paths from the finite element method correlated well with those deduced from gear experiments. For the test gear with a backup ratio (rim thickness divided by tooth height) of nib = 0.5, cracks initiating in the tooth fillet propagated to rim fractures when run at a speed of 10,000 rpm and became tooth fractures for speeds slower than 10,000 rpm for both the experiments and anal~ sis. From additional analysis, speed had little effect on crack propagation direction except when initial crack locations were near the tooth/rim fracture transition point for a given backup ratio. When at that point, higher speeds tended to promote rim fracture while lower speeds (or neglecting centrifugal force) produced tooth fractures.

Lewicki, David G.↗