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.

SourceResearch connectionFrom yearThrough year

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

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Small rocket flowfield diagnostic chambers

Instrumented and optically-accessible rocket chambers are being developed to be used for diagnostics of small rocket (less than 440 N thrust level) flowfields. These chambers are being tested to gather local fluid dynamic and thermodynamic flowfield data over a range of test conditions. This flowfield database is being used to better understand mixing and heat transfer phenomena in small rockets, influence the numerical modeling of small rocket flowfields, and characterize small rocket components. The diagnostic chamber designs include: a chamber design for gathering wall temperature profiles to be used as boundary conditions in a finite element heat flux model; a chamber design for gathering inner wall temperature and static pressure profiles; and optically-accessible chamber designs, to be used with a suite of laser-based diagnostics for gathering local species concentration, temperature, density, and velocity profiles. These chambers were run with gaseous hydrogen/gaseous oxygen (GH2/GO2) propellants, while subsequent versions will be run on liquid oxygen/hydrocarbon (LOX/HC) propellants. The purpose, design, and initial test results of these small rocket flowfield diagnostic chambers are summarized.

Morren, Sybil↗

NASA NTRS1979

Application of the method of integral relations to unsteady fluid flow problems with shocks

A mixed method using both the method of integral relations and the finite difference technique is developed for the solution of unsteady flow problems. The integral relations method is based upon a chosen interpolating function dependent only on the time domain. The resulting local semidiscrete finite element equations obtained are assembled into a global form. The spatial derivatives at the nodes are replaced by finite difference operators and the discretized nonlinear algebraic system is solved by an iterative scheme. Solutions are obtained for the one-dimensional gasdynamics equation, the one-dimensional wave equation and Burger's model of turbulence. Agreement with other numerical and analytical solutions is excellent for the gasdynamics problem and satisfactory for the Burger equation in cases of small viscous effects.

Wadia, A. R.↗

NASA NTRS1977

Finite element analysis of aeroacoustic jet-flap flows

A computational analysis was performed on the steady, turbulent aerodynamic flowfields associated with a jet-blown flap. For regions devoid of flow separation, a parabolic approximation to the governing time-averaged Navier-Stokes equations was applied. Numerical results are presented for the symmetry plane flow of a slot-nozzle planar jet flap geometry, including prediction of flowfield evolution within the secondary mixing region immediately downstream of the trailing edge. Using a two equation turbulence kinetic energy closure model, rapid generation and decay of large spatial gradients in mean and correlated fluctuating velocity components within the immediate wake region were predicted. Modifications to the turbulent flow structure, as induced by porous surface treatment of the flap, were evaluated. The recirculating flow within a representative discrete slot in the surface was evaluated, using the two dimensional, time-averaged Navier-Stokes equations.

Baker, A. J.↗

NASA NTRS2004

Penalty-Based Interface Technology for Prediction of Delamination Growth in Laminated Structures

An effective interface element technology has been developed for connecting and simulating crack growth between independently modeled finite element subdomains (e.g., composite plies). This method has been developed using penalty constraints and allows coupling of finite element models whose nodes do not necessarily coincide along their common interface. Additionally, the present formulation leads to a computational approach that is very efficient and completely compatible with existing commercial software. The present interface element has been implemented in the commercial finite element code ABAQUS as a User Element Subroutine (UEL), making it easy to test the approach for a wide range of problems. The interface element technology has been formulated to simulate delamination growth in composite laminates. Thanks to its special features, the interface element approach makes it possible to release portions of the interface surface whose length is smaller than that of the finite elements. In addition, the penalty parameter can vary within the interface element, allowing the damage model to be applied to a desired fraction of the interface between the two meshes. Results for double cantilever beam DCB, end-loaded split (ELS) and fixed-ratio mixed mode (FRMM) specimens are presented. These results are compared to measured data to assess the ability of the present damage model to simulate crack growth.

Averill, Ronald C.↗

NASA NTRS1993

Solution of mixed convection heat transfer from isothermal in-line fins

Transient and steady state combined natural and forced convective flows over two in-line finite thickness fins (louvers) in a vertical channel are numerically solved using two methods. The first method of solution is based on the 'Simple Arbitrary Lagrangian Eulerian' (SALE) technique which incorporates mainly two computational phases: (1) a Lagrangian phase in which the velocity field is updated by the effects of all forces, and (2) an Eulerian phase that executes all advective fluxes of mass, momentum and energy. The second method of solution uses the finite element code entitled FIDAP. In the first part, comparison of the results by FIDAP, SALE, and available experimental work were done and discussed for steady state forced convection over louvered fins. Good agreements were deduced between the three sets of results especially for the flow over a single fin. In the second part and in the absence of experimental literature, the numerical predictions were extended to the transient transports and to the opposing flow where pressure drop is reversed. Results are presented and discussed for heat transfer and pressure drop in assisting and opposing mixed convection flows.

Khalilollahi, Amir↗

NASA NTRS1993

Small hydrogen/oxygen rocket flowfield behavior from heat flux measurements

The mixing and heat transfer phenomena in small rocket flow fields with fuel film cooling is not well understood. An instrumented, water-cooled chamber with a gaseous hydrogen/gaseous oxygen injector was used to gather steady-state inner and outer wall temperature profiles. The chamber was tested at 414 kPa (60 psia) chamber pressure, from mixture ratios of 3.41 to 8.36. Sixty percent of the fuel was used for film cooling. These temperature profiles were used as boundary conditions in a finite element analysis program, MSC/NASTRAN, to calculate the local radial and axial heat fluxes in the chamber wall. The normal heat fluxes were then calculated and used as a diagnostic of the rocket's flow field behavior. The normal heat fluxes determined were on the order of 1.0 to 3.0 MW/meters squared (0.6 to 1.8 Btu/sec-inches squared). In the cases where mixture ratio was 5 or above, there was a sharp local heat flux maximum in the barrel section of the chamber. This local maximum seems to indicate a reduction or breakdown of the fuel film cooling layer, possibly due to increased mixing in the shear layer between the film and core flows. However, the flow was thought to be completely laminar, as the throat Reynolds numbers were below 50,000 for all the cases. The increased mixing in the shear layer in the higher mixture ratio cases appeared not to be due to the transition of the flow from laminar to turbulent, but rather due to increased reactions between the hydrogen film and oxidizer-rich core flows.

Reed, Brian D.↗

NASA NTRS1999

Experimentally Determined Crack Location and Mixed Mode Stress Intensity Factors

This presentation first examines five strain field equations: the two modes of the Westergaard equations, and a mixed mode; the generalized coordinate system, and the significance testing of coefficients. It also discusses the weighted multiple linear least squares method, which shows why it is difficult to find the crack tip. The presentation then discusses the experiments to examine stress intensity. Lastly, comparison of finite element method results and experimental results is reviewed.

Zachary, Loren W.↗

NASA NTRS2004

Quantifying Square Membrane Wrinkle Behavior Using MITC Shell Elements

For future membrane based structures, quantified predictions of membrane wrinkling behavior in terms of amplitude, angle and wavelength are needed to optimize the efficiency and integrity of such structures, as well as their associated control systems. For numerical analyses performed in the past, limitations on the accuracy of membrane distortion simulations have often been related to the assumptions made while using finite elements. Specifically, this work demonstrates that critical assumptions include: effects of gravity. supposed initial or boundary conditions, and the type of element used to model the membrane. In this work, a 0.2 square meter membrane is treated as a structural material with non-negligible bending stiffness. Mixed Interpolation of Tensorial Components (MTTC) shell elements are used to simulate wrinkling behavior due to a constant applied in-plane shear load. Membrane thickness, gravity effects, and initial imperfections with respect to flatness were varied in numerous nonlinear analysis cases. Significant findings include notable variations in wrinkle modes for thickness in the range of 50 microns to 1000 microns, which also depend on the presence of an applied gravity field. However, it is revealed that relationships between overall strain energy density for cases with differing initial conditions are independent of assumed initial con&tions. In addition, analysis results indicate that the relationship between amplitude scale (W/t) and structural scale (L/t) is linear in the presence of a gravity field.

Jacobson, Mindy B.↗

NASA NTRS2005

Skin-Stiffener Debond Prediction Based on Computational Fracture Analysis

Interlaminar fracture mechanics has proven useful for characterizing the onset of delaminations in composites and has been used with limited success primarily to investigate onset in fracture toughness specimens and laboratory size coupon type specimens. Future acceptance of the methodology by industry and certification authorities however, requires the successful demonstration of the methodology on structural level. For this purpose a panel was selected that is reinforced with stringers. Shear loading causes the panel to buckle and the resulting out-of-plane deformations initiate skin/stringer separation at the location of an embedded defect. For finite element analysis, the panel and surrounding load fixture were modeled with shell elements. A small section of the stringer foot and the panel in the vicinity of the embedded defect were modeled with a local 3D solid model. Across the width of the stringer foot the mixed-mode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with the mixed-mode failure criterion of the graphite/epoxy material. For small applied loads the failure index is well below one across the entire width. With increasing load the failure index approaches one first near the edge of the stringer foot from which delamination is expected to grow. With increasing delamination lengths the buckling pattern of the panel changes and the failure index increases which suggests that rapid delamination growth from the initial defect is to be expected.

Krueger, Ronald↗

NASA NTRS2006

Analysis of Composite Panel-Stiffener Debonding Using a Shell/3D Modeling Technique

Krueger, Ronald↗

NASA NTRS1979

Interlaminar crack growth in fiber reinforced composites during fatigue

This paper presents an investigation of interlaminar crack growth behavior in fiber-reinforced composites subjected to fatigue loading. In the experimental phase of the study, interlaminar crack propagation rates and mechanisms were determined for the cases of various geometries, laminate parameters and cyclic stress levels. An advanced singular hybrid-stress finite element method was used in conjunction with the experimental results to examine the local crack-tip behavior and to characterize the crack propagation during fatigue. Results elucidate the basic nature of the cyclic delamination damage and relate the interlaminar crack growth rate to the range of mixed-mode crack-tip stress intensity factors. The study provides fundamental insight into the problem, reveals several important features of the interlaminar fatigue failure, and should be of practical importance in selection, testing and design of composite materials.

Wang, S. S.↗

NASA NTRS1981

Interlaminar crack growth in fiber reinforced composites during fatigue, part 3

Interlaminar crack growth behavior in fiber-reinforced composites subjected to fatigue loading was investigated experimentally and theoretically. In the experimental phase, inter-laminar crack propagation rates and mechanisms were determined for the cases of various geometries, laminate parameters and cyclic stress levels. A singular hybrid-stress finite element method was used in conjuction with the experimental results to examine the local crack-tip behavior and to characterize the crack propagation during fatigue. Results elucidate the basic nature of the cyclic delamination damage, and relate the interlaminar crack growth rate to the range of mixed-mode crack-tip stress intensity factors. The results show that crack growth rates are directly related to the range of the mixed-mode cyclic stress intensity factors by a power law relationship.

Wang, S. S.↗

NASA NTRS1987

Shock structure in jets of arbitrary exit geometry

The shock structure and hydrodynamic stability of single and twin 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. Results are given for circular, elliptic and rectangular jets and the results compared with analysis and experiment. A finite difference technique is used to examine the effect of finite mixing layer thickness for a single jet. The hydrodynamic stability of twin circular jets is also considered. The variation of the stability characteristics are presented for the vortex sheet case for several separation distances. The formulation of the problem to account for finite thickness effects is also given.

Morris, P. J.↗

NASA NTRS1999

A Mixed-Mode I/II Fracture Criterion and Its Application in Crack Growth Predictions

A crack tip opening displacement (CTOD)-based, mixed mode fracture criterion is developed for predicting the onset and direction of crack growth. The criterion postulates that crack growth occurs in either the Mode I or Mode II direction, depending on whether the maximum in either the opening or the shear component of CTOD, measured at a specified distance behind the crack tip, attains a critical value. For crack growth direction prediction, the proposed CTOD criterion is shown to be equivalent to seven commonly used crack growth criteria under linearly elastic and asymptotic conditions. Under elastic-plastic conditions the CTOD criterion's prediction of the dependence of the crack growth direction on the crack-up mode mixity is in excellent agreement with the Arcan test results. Furthermore, the CTOD criterion correctly predicts the existence of a crack growth transition from mode I to mode II as the mode mixity approaches the mode II loading condition. The proposed CTOD criterion has been implemented in finite element crack growth simulation codes Z1P2DL and FRANC2DL to predict the crack growth paths in (a) a modified Arcan test specimen and fixture made of AL 2024-T34 and (b) a double cantilever beam (DCB) specimen made of AL 7050. A series of crack growth simulations have been carried out for the crack growth tests in the Arcan and DCB specimens and the results further demonstrate the applicability of the mixed mode CTOD fracture criterion crack growth predictions and residual strength analyses for airframe materials.

Sutton, Michael A.↗

NASA NTRS2016

An Enriched Shell Finite Element for Progressive Damage Simulation in Composite Laminates

A formulation is presented for an enriched shell nite element capable of progressive damage simulation in composite laminates. The element uses a discrete adaptive splitting approach for damage representation that allows for a straightforward model creation procedure based on an initially low delity mesh. The enriched element is veri ed for Mode I, Mode II, and mixed Mode I/II delamination simulation using numerical benchmark data. Experimental validation is performed using test data from a delamination-migration experiment. Good correlation was found between the enriched shell element model results and the numerical and experimental data sets. The work presented in this paper is meant to serve as a rst milestone in the enriched element's development with an ultimate goal of simulating three-dimensional progressive damage processes in multidirectional laminates.

McElroy, Mark W.↗

NASA NTRS1993

A new approach to mixed H2/H infinity controller synthesis using gradient-based parameter optimization methods

In the past few years, the mixed H(sub 2)/H-infinity control problem has been the object of much research interest since it allows the incorporation of robust stability into the LQG framework. The general mixed H(sub 2)/H-infinity design problem has yet to be solved analytically. Numerous schemes have considered upper bounds for the H(sub 2)-performance criterion and/or imposed restrictive constraints on the class of systems under investigation. Furthermore, many modern control applications rely on dynamic models obtained from finite-element analysis and thus involve high-order plant models. Hence the capability to design low-order (fixed-order) controllers is of great importance. In this research a new design method was developed that optimizes the exact H(sub 2)-norm of a certain subsystem subject to robust stability in terms of H-infinity constraints and a minimal number of system assumptions. The derived algorithm is based on a differentiable scalar time-domain penalty function to represent the H-infinity constraints in the overall optimization. The scheme is capable of handling multiple plant conditions and hence multiple performance criteria and H-infinity constraints and incorporates additional constraints such as fixed-order and/or fixed structure controllers. The defined penalty function is applicable to any constraint that is expressible in form of a real symmetric matrix-inequity.

Ly, Uy-Loi↗

NASA NTRSDate not supplied

Application of a Detached Eddy Simulation Approach with Finite-Rate Chemistry to Mars-Relevant Retropropulsion Operating Environments

Human-scale Mars vehicles will require retropropulsion for descent and landing, replacing heritage supersonic parachute systems with an extended phase of powered flight. Due to the limitations of terrestrial testing in Mars-relevant conditions, design and analysis will increasingly rely on computational modeling and simulation. This paper provides an overview of a computational campaign investigating the aerodynamics of a Mars lander concept along various points on a powered descent trajectory including supersonic, transonic, and subsonic conditions using finite-rate chemistry. Simulations using unstructured grids containing billions of elements are performed at scale using thousands of Graphics Processing Units, enabling run-times of a few days for each simulation presented. At each freestream condition, significant minor species concentrations are observed external to the nozzles in the large mixing region upstream of the vehicle. While the flowfields are highly non-stationary in all cases, the mean integrated forces and moments on the vehicle remain small in comparison to the deceleration provided through retropropulsion.

Ashley Korzun↗

NASA NTRS1994

Technology Benefit Estimator (T/BEST): User's Manual

The Technology Benefit Estimator (T/BEST) system is a formal method to assess advanced technologies and quantify the benefit contributions for prioritization. T/BEST may be used to provide guidelines to identify and prioritize high payoff research areas, help manage research and limited resources, show the link between advanced concepts and the bottom line, i.e., accrued benefit and value, and to communicate credibly the benefits of research. The T/BEST software computer program is specifically designed to estimating benefits, and benefit sensitivities, of introducing new technologies into existing propulsion systems. Key engine cycle, structural, fluid, mission and cost analysis modules are used to provide a framework for interfacing with advanced technologies. An open-ended, modular approach is used to allow for modification and addition of both key and advanced technology modules. T/BEST has a hierarchical framework that yields varying levels of benefit estimation accuracy that are dependent on the degree of input detail available. This hierarchical feature permits rapid estimation of technology benefits even when the technology is at the conceptual stage. As knowledge of the technology details increases the accuracy of the benefit analysis increases. Included in T/BEST's framework are correlations developed from a statistical data base that is relied upon if there is insufficient information given in a particular area, e.g., fuel capacity or aircraft landing weight. Statistical predictions are not required if these data are specified in the mission requirements. The engine cycle, structural fluid, cost, noise, and emissions analyses interact with the default or user material and component libraries to yield estimates of specific global benefits: range, speed, thrust, capacity, component life, noise, emissions, specific fuel consumption, component and engine weights, pre-certification test, mission performance engine cost, direct operating cost, life cycle cost, manufacturing cost, development cost, risk, and development time. Currently, T/BEST operates on stand-alone or networked workstations, and uses a UNIX shell or script to control the operation of interfaced FORTRAN based analyses. T/BEST's interface structure works equally well with non-FORTRAN or mixed software analysis. This interface structure is designed to maintain the integrity of the expert's analyses by interfacing with expert's existing input and output files. Parameter input and output data (e.g., number of blades, hub diameters, etc.) are passed via T/BEST's neutral file, while copious data (e.g., finite element models, profiles, etc.) are passed via file pointers that point to the expert's analyses output files. In order to make the communications between the T/BEST's neutral file and attached analyses codes simple, only two software commands, PUT and GET, are required. This simplicity permits easy access to all input and output variables contained within the neutral file. Both public domain and proprietary analyses codes may be attached with a minimal amount of effort, while maintaining full data and analysis integrity, and security. T/BESt's sotware framework, status, beginner-to-expert operation, interface architecture, analysis module addition, and key analysis modules are discussed. Representative examples of T/BEST benefit analyses are shown.

Generazio, Edward R.↗