Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixed finite element”

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 433 records · Page 24

Data-assisted combustion simulations with dynamic submodel assignment using random forests

This investigation outlines a data-assisted approach that employs random forest classifiers for local and dynamic submodel assignment in turbulent-combustion simulations. This method is demonstrated in simulations of a single-element GOX/GCH4 rocket combustor; a priori as well as a posteriori assessments are conducted to (i) evaluate the accuracy and adjustability of the classifier for targeting different quantities of interest (QoIs), and (ii) assess improvements, resulting from the data-assisted combustion model assignment, in predicting target QoIs during simulation runtime. Results from the a priori study show that random forests, trained with local flow properties as input variables and combustion model errors as training labels, assign three different combustion models – finite-rate chemistry (FRC), flamelet progress variable (FPV) model, and inert mixing (IM) – with reasonable classification performance even when targeting multiple QoIs. Applications in a posteriori studies demonstrate improved predictions from data-assisted simulations, in temperature and CO mass fraction, when compared with monolithic FPV calculations. An additional a posteriori data-assisted simulation of a modified configuration demonstrates that the present approach can be successfully applied to different configurations, as long as thermophysical behavior can be represented by the training data. Furthermore, these results demonstrate that this data-driven framework holds promise for dynamic combustion submodel assignments in reacting flow simulations.

42 ENGINEERING↗

Lunar and Planetary Science Conference, 17th, Houston, TX, Mar. 17-21, 1986, Proceedings. Part 2

The topics discussed in this volume include lunar endogenic rocks and processes; lunar regoliths and breccias; the terrestrial planets; shergottites; primitive materials, exposure, and atmospheres; and impacts and crater tectonics. Papers are presented on the petrology and geochemistry of alkali gabbronorites from lunar breccia 67975; the formation of Apollo 17 orange and black glass beads; mixing levels, the Apennine Front soil component, and compositional trends in the Apollo 15 soils; the meteorite component of Apollo 16 noritic impact melt breccias; and constraints on the lithospheric structure of Venus from mechanical models and tectonic surface features. Consideration is also given to a fractal interpretation of topography and geoid spectra on the earth, moon, Venus, and Mars; rare earth patterns in shergottite phosphates and residues; nuclide production by primary cosmic-ray protons; gas chromatographic instrumentation for the analysis of aerosols and gases in Titan's atmosphere; and finite-element models of non-Newtonian crater relaxation.

Ryder, Graham↗

Surface flaw reliability analysis of ceramic components with the SCARE finite element postprocessor program

The SCARE (Structural Ceramics Analysis and Reliability Evaluation) computer program on statistical fast fracture reliability analysis with quadratic elements for volume distributed imperfections is enhanced to include the use of linear finite elements and the capability of designing against concurrent surface flaw induced ceramic component failure. The SCARE code is presently coupled as a postprocessor to the MSC/NASTRAN general purpose, finite element analysis program. The improved version now includes the Weibull and Batdorf statistical failure theories for both surface and volume flaw based reliability analysis. The program uses the two-parameter Weibull fracture strength cumulative failure probability distribution model with the principle of independent action for poly-axial stress states, and Batdorf's shear-sensitive as well as shear-insensitive statistical theories. The shear-sensitive surface crack configurations include the Griffith crack and Griffith notch geometries, using the total critical coplanar strain energy release rate criterion to predict mixed-mode fracture. Weibull material parameters based on both surface and volume flaw induced fracture can also be calculated from modulus of rupture bar tests, using the least squares method with known specimen geometry and grouped fracture data. The statistical fast fracture theories for surface flaw induced failure, along with selected input and output formats and options, are summarized. An example problem to demonstrate various features of the program is included.

Gyekenyesi, John P.↗

Surface flaw reliability analysis of ceramic components with the SCARE finite element postprocessor program

The SCARE (Structural Ceramics Analysis and Reliability Evaluation) computer program on statistical fast fracture reliability analysis with quadratic elements for volume distributed imperfections is enhanced to include the use of linear finite elements and the capability of designing against concurrent surface flaw induced ceramic component failure. The SCARE code is presently coupled as a postprocessor to the MSC/NASTRAN general purpose, finite element analysis program. The improved version now includes the Weibull and Batdorf statistical failure theories for both surface and volume flaw based reliability analysis. The program uses the two-parameter Weibull fracture strength cumulative failure probability distribution model with the principle of independent action for poly-axial stress states, and Batdorf's shear-sensitive as well as shear-insensitive statistical theories. The shear-sensitive surface crack configurations include the Griffith crack and Griffith notch geometries, using the total critical coplanar strain energy release rate criterion to predict mixed-mode fracture. Weibull material parameters based on both surface and volume flaw induced fracture can also be calculated from modulus of rupture bar tests, using the least squares method with known specimen geometry and grouped fracture data. The statistical fast fracture theories for surface flaw induced failure, along with selected input and output formats and options, are summarized. An example problem to demonstrate various features of the program is included.

Gyekenyesi, J. P.↗

Probabilistic and Possibilistic Analyses of the Strength of a Bonded Joint

The effects of uncertainties on the predicted strength of a single lap shear joint are examined. Probabilistic and possibilistic methods are used to account for uncertainties. A total of ten variables are assumed to be random, with normal distributions. Both Monte Carlo Simulation and the First Order Reliability Method are used to determine the probability of failure. Triangular membership functions with upper and lower bounds located at plus or minus three standard deviations are used to model uncertainty in the possibilistic analysis. The alpha cut (or vertex) method is used to evaluate the possibility of failure. Linear and geometrically nonlinear finite element analyses are used calculate the response of the joint; fracture in the adhesive and material strength failure in the strap are used to evaluate its strength. Although probabilistic and possibilistic analyses provide significantly more information than do conventional deterministic analyses, they are computationally expensive. A novel scaling approach is developed and used to substantially reduce the computational cost of the probabilistic and possibilistic analyses. The possibilistic approach for treating uncertainties appears to be viable during the conceptual and preliminary design stages when limited data are available and high accuracies are not needed. However, this viability is mixed with several cautions that are discussed herein.

Stroud, W. Jefferson↗

A linear shock cell model for jets of arbitrary exit geometry

The shock cell structures of single supersonic non-ideally expanded jets with arbitrary exit geometry are studied. Both vortex sheets and realistic mean profiles are considered for the jet shear layer. The boundary element method is used to predict the shock spacing and screech tones in a vortex sheet model of a single jet. This formulation enables the calculations to be performed only on the vortex sheet. This permits the efficient and convenient study of complicated jet geometries. Results are given for circular, elliptic and rectangular jets and the results are compared with analysis and experiment. The agreement between the predictions and measurements is very good but depends on the assumptions made to predict the geometry of the fully expanded jet. A finite diffference technique is used to examine the effect of finite mixing layer thickness for a single jet. The finite thickness of the mixing layer is found to decrease the shock spacing by approximately 20 percent over the length of the jet potential core.

Morris, P. J.↗

TPSAS-NF1676L-13122-DND

An approach for assessing the delamination propagation and growth capabilities in commercial finite element codes was developed and demonstrated for the Virtual Crack Closure Technique (VCCT) implementations in ABAQUS?. The Double Cantilever Beam (DCB) specimen was chosen as an example. First, benchmark results to assess delamination propagation capabilities under static loading were created using models simulating specimens with different delamination lengths. Second, starting from an initially straight front, the delamination was allowed to propagate based on the algorithms implemented in the commercial finite element software. The load-displacement relationship obtained from the propagation analysis results and the benchmark results were compared. Good agreements could be achieved by selecting the appropriate input parameters, which were determined in an iterative procedure. A benchmark case to assess the fatigue growth prediction capabilities was also created. Input control parameters were varied to study the effect on the computed delamination growth. The benchmarks enabled the selection of the appropriate input parameters that yielded good agreements between the results obtained from the automated analysis and the benchmark results. Overall, the results are encouraging but further assessment for mixed-mode delamination is required.

Ronald Krueger↗

Subsurface Stress Fields In Single Crystal (Anisotropic) Contacts

Single crystal superalloy turbine blades used in high pressure turbomachinery are subject to conditions of high temperature, triaxial steady and alternating stresses, fretting stresses in the blade attachment and damper contact locations, and exposure to high-pressure hydrogen. The blades are also subjected to extreme variations in temperature during start-up and shutdown transients. The most prevalent HCF failure modes observed in these blades during operation include crystallographic crack initiation/propagation on octahedral planes, and noncrystallographic initiation with crystallographic growth. Numerous cases of crack initiation and crack propagation at the blade leading edge tip, blade attachment regions, and damper contact locations have been documented. Understanding crack initiation/propagation under mixed-mode loading conditions is critical for establishing a systematic procedure for evaluating HCF life of single crystal turbine blades. This paper presents analytical and numerical techniques for evaluating two and three dimensional subsurface stress fields in anisotropic contacts. The subsurface stress results are required for evaluating contact fatigue life at damper contacts and dovetail attachment regions in single crystal nickel-base superalloy turbine blades. An analytical procedure is , presented, for evaluating the subsurface stresses in the elastic half-space, using a complex potential method outlined by Lekhnitskii. Numerical results are presented for cylindrical and spherical anisotropic contacts, using finite element analysis. Effects of crystal orientation on stress response and fatigue life are examined.

Arakere, Nagaraj K.↗

High frequency dynamic engine simulation

A digital computer simulation of a mixed flow, twin spool turbofan engine was assembled to evaluate and improve the dynamic characteristics of the engine simulation to disturbance frequencies of at least 100 Hz. One dimensional forms of the dynamic mass, momentum and energy equations were used to model the engine. A TF30 engine was simulated so that dynamic characteristics could be evaluated against results obtained from testing of the TF30 engine at the NASA Lewis Research Center. Dynamic characteristics of the engine simulation were improved by modifying the compression system model. Modifications to the compression system model were established by investigating the influence of size and number of finite dynamic elements. Based on the results of this program, high frequency engine simulations using finite dynamic elements can be assembled so that the engine dynamic configuration is optimum with respect to dynamic characteristics and computer execution time. Resizing of the compression systems finite elements improved the dynamic characteristics of the engine simulation but showed that additional refinements are required to obtain close agreement simulation and actual engine dynamic characteristics.

Schuerman, J. A.↗

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↗

Subsurface Stress Fields in FCC Single Crystal Anisotropic Contacts

Single crystal superalloy turbine blades used in high pressure turbomachinery are subject to conditions of high temperature, triaxial steady and alternating stresses, fretting stresses in the blade attachment and damper contact locations, and exposure to high-pressure hydrogen. The blades are also subjected to extreme variations in temperature during start-up and shutdown transients. The most prevalent high cycle fatigue (HCF) failure modes observed in these blades during operation include crystallographic crack initiation/propagation on octahedral planes, and non-crystallographic initiation with crystallographic growth. Numerous cases of crack initiation and crack propagation at the blade leading edge tip, blade attachment regions, and damper contact locations have been documented. Understanding crack initiation/propagation under mixed-mode loading conditions is critical for establishing a systematic procedure for evaluating HCF life of single crystal turbine blades. This paper presents analytical and numerical techniques for evaluating two and three dimensional subsurface stress fields in anisotropic contacts. The subsurface stress results are required for evaluating contact fatigue life at damper contacts and dovetail attachment regions in single crystal nickel-base superalloy turbine blades. An analytical procedure is presented for evaluating the subsurface stresses in the elastic half-space, based on the adaptation of a stress function method outlined by Lekhnitskii. Numerical results are presented for cylindrical and spherical anisotropic contacts, using finite element analysis (FEA). Effects of crystal orientation on stress response and fatigue life are examined. Obtaining accurate subsurface stress results for anisotropic single crystal contact problems require extremely refined three-dimensional (3-D) finite element grids, especially in the edge of contact region. Obtaining resolved shear stresses (RSS) on the principal slip planes also involves considerable post-processing work. For these reasons it is very advantageous to develop analytical solution schemes for subsurface stresses, whenever possible.

Arakere, Nagaraj K.↗

Thermal and dynamical regimes of single- and two-phase magmatic flow in dikes

The coupling between thermal and dynamical regimes of single- and two-phase magmatic flow in dikes, due to temperature-dependent viscosity and dissipation, was investigated using finite element calculations of magma flow in dikelike channels with length-to-width ratios of 1000:1 or more. Solutions of the steady state equations governing magma flow are obtained for a variety of conditions ranging from idealized plane-parallel models to cases involving nonparallel geometry and two-phase flows. The implications of the numerical simulations for the dynamics of flow in a dike-reservoir system and the consequences of dike entrance conditions on magmatic storage are discussed. Consideration is also given to an unmixing/self-lubrication mechanism which may be important for the lubrication of silicic magmas rising to the earth's surface in mixed magma ascent scenarios, which naturally segregates magma mixtures of two components with differing viscosities to minimize the driving pressure gradient.

Carrigan, Charles R.↗

Venus tectonic styles and crustal differentiation

Two of the most important constraints are known from Pioneer Venus data: the lack of a system of spreading rises, indicating distributed deformation rather than plate tectonics; and the high gravity/topography ratio, indicating the absence of an asthenosphere. In addition, the high depth/diameter ratios of craters on Venus indicate that Venus probably has no more crust than Earth. The problems of the character of tectonics and crustal formation and recycling are closely coupled. Venus appears to lack a recycling mechanism as effective as subduction, but may also have a low rate of crustal differentiation because of a mantle convection pattern that is more distributed, less concentrated, than Earth's. Distributed convection, coupled with the nonlinear dependence of volcanism on heat flow, would lead to much less magmatism, despite only moderately less heat flow, compared to Earth. The plausible reason for this difference in convective style is the absence of water in the upper mantle of Venus. We have applied finite element modeling to problems of the interaction of mantle convection and crust on Venus. The main emphasis has been on the tectonic evolution of Ishtar Terra, as the consequence of convergent mantle flow. The early stage evolution is primarily mechanical, with crust being piled up on the down-stream side. Then the downflow migrates away from the center. In the later stages, after more than 100 m.y., thermal effects develop due to the insulating influence of the thickened crust. An important feature of this modeling is the entrainment of some crustal material in downflows. An important general theme in both convergent and divergent flows is that of mixing vs. stratification. Models of multicomponent solid-state flow obtain that lower-density crustal material can be entrained and recycled, provided that the ration of low-density to high-density material is small enough (as in subducted slabs on Earth). The same considerations should apply in upflows; a small percent of partial melt may be carried along with its matrix and never escape to the surface. Models that assume melt automatically rising to the crust and no entrainment or other mechanism of recycling lower-density material obtain oscillatory behavior, because it takes a long time for heat to build up enough to overcome a Mg-rich low-density residuum. However, these models develop much thicker crust than consistent with estimates from crater depth/diameter ratios.

Kaula, W. M.↗

Modeling composite and fluffy grains: The effects of porosity

Recent studies of interplanetary and interstellar dust provide evidence that cosmic dust grains are fluffy, composite objects, highlighting the need for models of the electromagnetic scattering by these grains. Effective medium theory (EMT) with Mie-type series solutions has been used to explore the effects of porosity which would be important in composite dust particles. While this indirect approach is both flexible and computationally efficient, it is not necessarily a good approximation. The need for EMT and its rather restrictive assumptions may be circumvented through a direct computation of the scattering properties via finite element methods, such as the discrete dipole approximation (DDA). Recently, the utility of the DDA method has been advanced significantly through improvements in theory, in numerical algorithms, and in computer hardware. Extensive calculations with the DDA method are used here to examine more directly the effects of porosity. A particular emphasis is placed upon developing a valid methodology. For both solid and porous targets we establish both numerical and physical convergence properties over the range of size parameter that is required for our study. DDA cross sections for grains with a range of porosity are compared to those computed by the EMT/series expansion technique to examine the applicability of several mixing rules, including two extensions of the Bruggeman rule. We show that for particles with Rayleigh vacuum inclusions, the extension proposed by Rouleau & Martin is quite successful. We also investigate the effects of larger, non-Rayleigh vacuum inclusions for various levels of porosity and find that they can be significant.

Wolff, Michael J.↗

Powdered aluminum and oxygen rocket propellants: Subscale combustion experiments

Aluminum combined with oxygen has been proposed as a potential lunar in situ propellant for ascent/descent and return missions for future lunar exploration. Engine concepts proposed to use this propellant have not previously been demonstrated, and the impact on performance from combustion and two-phase flow losses could only be estimated. Therefore, combustion tests were performed for aluminum and aluminum/magnesium alloy powders with oxygen in subscale heat-sink rocket engine hardware. The metal powder was pneumatically injected, with a small amount of nitrogen, through the center orifice of a single element O-F-O triplet injector. Gaseous oxygen impinged on the fuel stream. Hot-fire tests of aluminum/oxygen were performed over a mixture ratio range of 0.5 to 3.0, and at a chamber pressure of approximately 480 kPa (70 psia). The theoretical performance of the propellants was analyzed over a mixture ratio range of 0.5 to 5.0. In the theoretical predictions the ideal one-dimensional equilibrium rocket performance was reduced by loss mechanisms including finite rate kinetics, two-dimensional divergence losses, and boundary layer losses. Lower than predicted characteristic velocity and specific impulse performance efficiencies were achieved in the hot-fire tests, and this was attributed to poor mixing of the propellants and two-phase flow effects. Several tests with aluminum/9.8 percent magnesium alloy powder did not indicate any advantage over the pure aluminum fuel.

Meyer, Mike L.↗

Coupled Fluid-Structure Interaction Analysis of Solid Rocket Motor with Flexible Inhibitors

Flexible inhibitors are generally used in solid rocket motors (SRMs) as a means to control the burning of propellant. Vortices generated by the flow of propellant around the flexible inhibitors have been identified as a driving source of instabilities that can lead to thrust oscillations in launch vehicles. Potential coupling between the SRM thrust oscillations and structural vibration modes is an important risk factor in launch vehicle design. As a means to predict and better understand these phenomena, a multidisciplinary simulation capability that couples the NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This capability is crucial to the development of NASA's new space launch system (SLS). This paper summarizes the efforts in applying the coupled software to demonstrate and investigate fluid-structure interaction (FSI) phenomena between pressure waves and flexible inhibitors inside reusable solid rocket motors (RSRMs). The features of the fluid and structural solvers are described in detail, and the coupling methodology and interfacial continuity requirements are then presented in a general Eulerian-Lagrangian framework. The simulations presented herein utilize production level CFD with hybrid RANS/LES turbulence modeling and grid resolution in excess of 80 million cells. The fluid domain in the SRM is discretized using a general mixed polyhedral unstructured mesh, while full 3D shell elements are utilized in the structural domain for the flexible inhibitors. Verifications against analytical solutions for a structural model under a steady uniform pressure condition and under dynamic modal analysis show excellent agreement in terms of displacement distribution and eigenmode frequencies. The preliminary coupled results indicate that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor. This insight could have profound implications for SRM and flexible inhibitor designs for current and future launch vehicles including SLS.

Yang, H. Q.↗

In Search of a Time Efficient Approach to Crack and Delamination Growth Predictions in Composites

Analysis benchmarking was used to assess the accuracy and time efficiency of algorithms suitable for automated delamination growth analysis. First, the Floating Node Method (FNM) was introduced and its combination with a simple exponential growth law (Paris Law) and Virtual Crack Closure technique (VCCT) was discussed. Implementation of the method into a user element (UEL) in Abaqus/Standard(Registered TradeMark) was also presented. For the assessment of growth prediction capabilities, an existing benchmark case based on the Double Cantilever Beam (DCB) specimen was briefly summarized. Additionally, the development of new benchmark cases based on the Mixed-Mode Bending (MMB) specimen to assess the growth prediction capabilities under mixed-mode I/II conditions was discussed in detail. A comparison was presented, in which the benchmark cases were used to assess the existing low-cycle fatigue analysis tool in Abaqus/Standard(Registered TradeMark) in comparison to the FNM-VCCT fatigue growth analysis implementation. The low-cycle fatigue analysis tool in Abaqus/Standard(Registered TradeMark) was able to yield results that were in good agreement with the DCB benchmark example. Results for the MMB benchmark cases, however, only captured the trend correctly. The user element (FNM-VCCT) always yielded results that were in excellent agreement with all benchmark cases, at a fraction of the analysis time. The ability to assess the implementation of two methods in one finite element code illustrated the value of establishing benchmark solutions.

Krueger, Ronald↗

An analytical and experimental investigation of edge delamination in laminates subjected to tension, bending, and torsion

An integrated two-dimensional finite element was developed to calculate interlaminar stresses and strain energy release rates for the study of delamination in composite laminates subjected to uniaxial tension, bending, and torsion loads. Addressed are the formulation, implementation, and verification of the model. Parametric studies were conducted on the effect of Poisson's ratio mismatch between plies and the stacking sequence on interlaminar stress, and on the effect of delamination opening height and delamination length, due to bending, on strain energy release rate for various laminates. A comparison of strain energy release rates in all-graphite and graphite/glass hybrid laminates is included. The preliminary results of laminates subjected to torsion are also included. Fatigue tension tests were conducted on Mode 1 and mixed mode edge-delamination coupons to establish the relationship between fatigue load vs. onset of delamination cycle. The effect on the fatigue delamination onset of different frequencies (1 and 5 Hz) was investigated for glass, graphite,and their hybrid laminates. Although a 20 percent increase in the static onset-of-delamination strength and a 10 percent increase in ultimate strength resulted from hybridizing the all-graphite laminate with a 90 deg glass ply, the fatigue onset is lower in the hybrid laminate than in the all-graphite laminate.

Chan, Wen S.↗