Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element methods”

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 505 records · Page 28

The stability of pseudospectral-Chebyshev methods

The pseudospectral-Chebyshev methods are shown to be convergent in variable coefficient problems and, in some cases, hyperbolic problems. The analysis demonstrates that the rate of convergence is greater for finite difference methods or the finite element method. For a single first-order hyperbolic equation, the method is seen as remaining stable even when the coefficient changes sign, although in this case it is specified that care must be taken to have adequate spatial resolution. It is noted that this fact, combined with the fact that collocation methods are easy to apply in the nonlinear case, shows that the pseudospectral method is in general preferable to the Galerkin or Tau methods.

Gottlieb, D.↗

Megamechanics research Consortium

The mechanical behavior of large space lattice systems is investigated. Continuum modeling methods such as the Load Correction Method (LCM) and finite element operations aid in solving structural analysis problems. Static analysis of frames and trusses, dynamic analysis of trusses, analysis of thermal and mechanical loads on frames and trusses, and load corrections for three dimensional complex geometries is possible by using the LCM. The finite element method and the LCM are used to extend Saint-Venant's principle via the Boussinesq problem into the regime of large repetitive lattice systems.

Broome, T. H., Jr.↗

An advanced probabilistic structural analysis method for implicit performance functions

In probabilistic structural analysis, the performance or response functions usually are implicitly defined and must be solved by numerical analysis methods such as finite element methods. In such cases, the most commonly used probabilistic analysis tool is the mean-based, second-moment method which provides only the first two statistical moments. This paper presents a generalized advanced mean value (AMV) method which is capable of establishing the distributions to provide additional information for reliability design. The method requires slightly more computations than the second-moment method but is highly efficient relative to the other alternative methods. In particular, the examples show that the AMV method can be used to solve problems involving non-monotonic functions that result in truncated distributions.

Wu, Y.-T.↗

A new procedure for calculating contact stresses in gear teeth

A numerical procedure for evaluating and monitoring contact stresses in meshing gear teeth is discussed. The procedure is intended to extend the range of applicability and to improve the accuracy of gear contact stress analysis. The procedure is based upon fundamental solution from the theory of elasticity. It is an iterative numerical procedure. The method is believed to have distinct advantages over the classical Hertz method, the finite-element method, and over existing approaches with the boundary element method. Unlike many classical contact stress analyses, friction effects and sliding are included. Slipping and sticking in the contact region are studied. Several examples are discussed. The results are in agreement with classical results. Applications are presented for spur gears.

Somprakit, Paisan↗

Modeling of crack bridging in a unidirectional metal matrix composite

The effective fatigue crack driving force and crack opening profiles were determined analytically for fatigue tested unidirectional composite specimens exhibiting fiber bridging. The crack closure pressure due to bridging was modeled using two approaches; the fiber pressure model and the shear lag model. For both closure models, the Bueckner weight function method and the finite element method were used to calculate crack opening displacements and the crack driving force. The predicted near crack tip opening profile agreed well with the experimentally measured profiles for single edge notch SCS-6/Ti-15-3 metal matrix composite specimens. The numerically determined effective crack driving force, Delta K(sup eff), was calculated using both models to correlate the measure crack growth rate in the composite. The calculated Delta K(sup eff) from both models accounted for the crack bridging by showing a good agreement between the measured fatigue crack growth rates of the bridged composite and that of unreinforced, unbridged titanium matrix alloy specimens.

Ghosn, Louis J.↗

Monitoring of contact stresses in advanced propulsion systems

This paper discusses a new and innovative numerical procedure for evaluating and monitoring contact stresses in the mating surfaces of propulsion systems. The procedure is intended to extend the range of applicability and the accuracy of contact stress analyses. Originally proposed by Burton Paul, et. al., (1981) the procedure is based upon fundamental solutions from the theory of elasticity. It is an iterative numerical procedure. It is believed that the method has distinct advantages over the classical Hertz method, the finite-element method and over existing approaches with the boundary element method. Several examples are discussed. Favorable comparisons are obtained with classical results. Applications with mating element of propulsion systems are discussed.

Somprakit, Paisan↗

Modeling of crack bridging in a unidirectional metal matrix composite

The effective fatigue crack driving force and crack opening profiles were determined analytically for fatigue tested unidirectional composite specimens exhibiting fiber bridging. The crack closure pressure due to bridging was modeled using two approaches: the fiber pressure model and the shear lag model. For both closure models, the Bueckner weight function method and the finite element method were used to calculate crack opening displacements and the crack driving force. The predicted near crack tip opening profile agreed well with the experimentally measured profiles for single edge notch SCS-6/Ti-15-3 metal matrix composite specimens. The numerically determined effective crack driving force, Delta K(eff), was calculated using both models to correlate the measure crack growth rate in the composite. The calculated Delta K(eff) from both models accounted for the crack bridging by showing a good agreement between the measured fatigue crack growth rates of the bridged composite and that of unreinforced, unbridged titanium matrix alloy specimens.

Ghosn, Louis J.↗

Contact stress analysis of spiral bevel gears using nonlinear finite element static analysis

A procedure is presented for performing three-dimensional stress analysis of spiral bevel gears in mesh using the finite element method. The procedure involves generating a finite element model by solving equations that identify tooth surface coordinates. Coordinate transformations are used to orientate the gear and pinion for gear meshing. Contact boundary conditions are simulated with gap elements. A solution technique for correct orientation of the gap elements is given. Example models and results are presented.

Bibel, G. D.↗

Theoretical and numerical difficulties in 3-D vector potential methods in finite element magnetostatic computations

This paper describes the results of application of three well known 3D magnetic vector potential (MVP) based finite element formulations for computation of magnetostatic fields in electrical devices. The three methods were identically applied to three practical examples, the first of which contains only one medium (free space), while the second and third examples contained a mix of free space and iron. The first of these methods is based on the unconstrained curl-curl of the MVP, while the second and third methods are predicated upon constraining the divergence of the MVP 10 zero (Coulomb's Gauge). It was found that the two latter methods cease to give useful and meaningful results when the global solution region contains a mix of media of high and low permeabilities. Furthermore, it was found that their results do not achieve the intended zero constraint on the divergence of the MVP.

Demerdash, N. A.↗

Contact stress analysis of spiral bevel gears using nonlinear finite element static analysis

A procedure is presented for performing three-dimensional stress analysis of spiral bevel gears in mesh using the finite element method. The procedure involves generating a finite element model by solving equations that identify tooth surface coordinates. Coordinate transformations are used to orientate the gear and pinion for gear meshing. Contact boundary conditions are simulated with gap elements. A solution technique for correct orientation of the gap elements is given. Example models and results are presented.

Bibel, G. D.↗

Hybrid equations of motion for flexible multibody systems using quasi-coordinates

A variety of engineering systems, such as automobiles, aircraft, rotorcraft, robots, spacecraft, etc., can be modeled as flexible multibody systems. The individual flexible bodies are in general characterized by distributed parameters. In most earlier investigations they were approximated by some spatial discretization procedure, such as the classical Rayleigh-Ritz method or the finite element method. This paper presents a mathematical formulation for distributed-parameter multibody systems consisting of a set of hybrid (ordinary and partial) differential equations of motion in terms of quasi-coordinates. Moreover, the equations for the elastic motions include rotatory inertia and shear deformation effects. The hybrid set is cast in state form, thus making it suitable for control design.

Meirovitch, L.↗

Theoretical and Experimental Study of Microstrip-to-Slot Line Uniplanar Transition

Recent advances in MMCI technology make it possible to construct transitions from CPW-to-microstrip with via hole, microstrip-to-slot line and microshield line-to-CPW all of which have potential applications in the feed network of antennas. In this study we investigate the characteristics of the microstrip-to-slot line uniplanar transition using the finite element methods (FEM) and finite difference time domain (FDTD) techniques, and compared the theoretical results with the measurements. In both cases, the results agree with the measurements within a few percent.

Yook, Jong-Gwan↗

Unified CFD Methods Via Flowfield-Dependent Variation Theory

This paper addresses the flowfield-dependent variation (FDV) methods in which complex physical phenomena are taken into account in the final form of partial differential equations to be solved so that finite difference methods (FDM) or finite element methods (FEM) themselves will not dictate the physics, but rather are no more than simply the options how to discretize between adjacent nodal points or within an element. The variation parameters introduced in the formulation are calculated from the current flowfield based on changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkohler numbers between adjacent nodal points, which play many significant roles such as adjusting the governing equations (hyperbolic, parabolic, and/or e!liptic), resolving various physical phenomena, and controlling the accuracy and stability of the numerical solution. The theory is verified by a number of example problems addressing the physical implications of the variation parameters which resemble the flowfield itself, shock capturing mechanism, transitions and interactions between inviscid/viscous, compressibility/incompressibility, and laminar/turbulent flows.

Chung, T. J.↗

Contact Stress Analysis of Spiral Bevel Gears Using Finite Element Analysis

A procedure is presented for performing three-dimensional stress analysis of spiral bevel gears in mesh using the finite element method. The procedure involves generating a finite element model by solving equations that identify tooth surface coordinates. Coordinate transformations are used to orientate the gear and pinion for gear meshing. Contact boundary conditions are simulated with gap elements. A solution technique for correct orientation of the gap elements is given. Example models and results are presented.

Bibel, G. D.↗

How Certain Physical Considerations Impact Aerostructural Wing Optimization

Wing design optimization has been studied extensively and is of continued interest as optimization tools are developed and become more accessible. In each of these studies, certain assumptions and simplifications are made to make the design problem tractable. However, it is difficult to find systematic studies in which several considerations are added or removed one at a time to study how much impact they have. In this work, we examine how certain physical considerations (viscous drag, wave drag, thrust loads, and inertial relief from structural, fuel, and engine masses), impact the aerostructural optimization results for three distinct aircraft wings. The goal is to help develop a rough idea of how important these physical considerations are. We do this using gradient-based optimization and a multidisciplinary design optimization framework, OpenMDAO. We use the open-source tool OpenAeroStruct that couples a vortex lattice method to a finite element method. We establish a baseline aerostructural design optimization problem then perform a series of optimizations, each with one physical consideration removed from the baseline case. We find that depending on the size of the aircraft and flight conditions, the importance of some of these physical considerations varies considerably whereas the importance of others do not. Specifically, the optimal designs change radically without proper viscous and wave drag considerations and smaller aircraft with more distributed propulsion are more affected by the inclusion of engine loads.

Aeropropulsion↗

Pronghorn Porous Media Model Validation with Pressure Drop Measurements

The verification and validation (V&V) of Pronghorn is imperative to assert its accuracy when predicting the fluid velocity, temperature, and pressure in high temperature gas-cooled reactors. Pronghorn is a coarse-mesh, intermediate-fidelity, and multidimensional thermal-hydraulics (TH) code developed by the Idaho National Laboratory (INL). New pebble bed experiments are used to observe the details of the fluid motion and pressure drop in the porous bed under the reactor normal operation. This paper focuses on the validation of the Pronghorn compressible and incompressible Navier-Stokes equations using the pressure drop measurements performed at the engineering-scale pebble bed facility at the Texas A&M university (TAMU). Various pressure drop correlations and porosity functions are implemented in both Pronghorn and STAR-CCM+ to compare the pressure drop due to the combined viscous and inertial resistances in the porous bed. The correlations accounting for the near-wall effect are also utilized to observe if the pressure drop estimates can be improved. Pronghorn porous media models predict the pressure drop well relative to the STAR-CCM+ simulation results and 1D correlations, and both the finite element method (FEM) and finite volume method (FVM) perform accurately. Pronghorn models are also validated with the experimental measurements given the different Reynolds number ranges and specific aspect ratios. The likelihood of the statistical significance between the pressure drop measurements and specific correlations or simulations is low provided that the overlap of their confidence intervals is more than the half of a single arm. Several validation metrics are reasonable in regard to the similar studies from other literature. The precise average pebble bed porosity estimation has much impact on the pressure drop, and the Foumeny and Montillet (dense packing) models carry out the accurate pressure drop prediction by considering the near-wall effect.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗

Practical model reduction for flexible structures

The most accepted methods for developing models of flexible structures for the analysis, design and simulation of control systems are finite element methods and extraction from experimental data. Using either of these methods, models are often produced with modes that do not significantly add to the fidelity of the model. This paper presents techniques for eliminating these modes. Particular attention is paid to multiple-input, multiple-output systems. First, control system models developed for flexible structures, using finite element methods and experimental data, are briefly discussed. the shortcomings of using models with unnecessary and/or residual modes are delineated. Then, two techniques for reducing the order of models are presented; the first is applicable to single-input, single-output systems and the second for multiple-input, multiple-output systems. Finally, both are illustrated using model data from the proposed NASA Shuttle-C.

Mitchell, Jerrel R.↗