Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Triangle elements”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Triangular spectral elements for incompressible fluid flow

We discuss the use of triangular elements in the spectral element method for direct simulation of incompressible flow. Triangles provide much greater geometric flexibility and are better conditioned and more accurate than quadrilateral elements when small angles arise. We employ a family of tensor product algorithms for triangles, allowing triangular elements to be handled with comparable arithmetic complexity to quadrilateral elements. The triangular discretizations are applied to the Poisson equation and are validated. The triangular discretizations are then applied to the incompressible Navier-Stokes equations, and a laminar channel flow solution is given. The new triangular spectral elements can be combined with standard quadrilateral elements, yielding a general and flexible high order method for complex geometries in two dimensions.

Mavriplis, Catherine↗
NASA NTRS2015

Computations of Disturbance Amplification Behind Isolated Roughness Elements and Comparison with Measurements

Computations are performed to study laminar-turbulent transition due to isolated roughness elements in boundary layers at Mach 3.5 and 5.95, with an emphasis on flow configurations for which experimental measurements from low disturbance wind tunnels are available. The Mach 3.5 case corresponds to a roughness element with right-triangle planform with hypotenuse that is inclined at 45 degrees with respect to the oncoming stream, presenting an obstacle with spanwise asymmetry. The Mach 5.95 case corresponds to a circular roughness element along the nozzle wall of the Purdue BAMQT wind tunnel facility. In both cases, the mean flow distortion due to the roughness element is characterized by long-lived streamwise streaks in the roughness wake, which can support instability modes that did not exist in the absence of the roughness element. The linear amplification characteristics of the wake flow are examined towards the eventual goal of developing linear growth correlations for the onset of transition.

Choudhari, Meelan↗
NASA NTRS1978

Finite element analysis of periodic transonic flow problems

Flow about an oscillating thin airfoil in a transonic stream was considered. It was assumed that the flow field can be decomposed into a mean flow plus a periodic perturbation. On the surface of the airfoil the usual Neumman conditions are imposed. Two computer programs were written, both using linear basis functions over triangles for the finite element space. The first program uses a banded Gaussian elimination solver to solve the matrix problem, while the second uses an iterative technique, namely SOR. The only results obtained are for an oscillating flat plate.

Fix, G. J.↗
NASA NTRS2005

Generating Solid Models from Topographical Data

A method of generating solid models of terrain involves the conversion of topographical data into a form useable by a rapid-prototyping (RP) machine. The method was developed to enable the use of the RP machine to make solid models of Martian terrain from Mars Orbiter laser-altimeter topographical data. The method is equally applicable to the generation of models of the terrains of other astronomical bodies, including other planets, asteroids, and Earth. Topographical data describe a terrain in terms of a set of three-dimensional coordinates [e.g., Cartesian (x,y,z) or polar (latitude, longitude, radius) coordinates] of points or nodes on the terrain surface. The input data for the RP machines are required to provide a three-dimensional description, not of a single surface, but of a volume in this case, a ground volume that underlies the terrain surface. The description is required to be in the form of triangular elements that connect the nodes of all the surfaces and that completely bound the volume, with no open areas, no overlap of triangles, and no extraneous geometric elements. The software used in the present model-generation method was written in IDL - an advanced programming language that affords a number of tools, including subroutines that triangularize surfaces. The software creates a volume from the topographical surface data by adding sides to the edges of the terrain surface and joining the sides with a bottom surface. Each of the sides is triangularized by use of IDL subroutines, and then the software searches for extraneous elements and removes them. Topographical data are usually presented in a grid corresponding to polar coordinates, so that a model generated from such data is equivalent to a topographical map in Mercator projection. However an RP machine is fully capable of including the curvature of a planetary body in a model that it makes. Therefore, the software also offers a capability to transform the topographical data to a projection onto a surface having a curvature corresponding to that of the surface of the modeled planet.

Keller, John W.↗
NASA NTRS1973

New plate and shell elements for NASTRAN

A new higher order triangular plate-bending finite element is presented which possesses high accuracy for practical mesh subdivisions and which uses only translations and rotations as grid point degrees of freedom. The element has 18 degrees of freedom, the transverse displacement and two rotations at the vertices and mid-side grid points of the triangle. The transverse displacement within the element is approximated by a quintic polynomial; the bending strains thus vary cubically within the element. Transverse shear flexibility is taken into account in the stiffness formulation. Two examples of static and dynamic analysis are included to show the behavior of the element.

Narayanaswami, R.↗
NASA NTRS1988

A mixed formulaton of C(0)-linear triangular plate/shell element - The role of edge shear constraints

A simple and effective linear C(0)-tiangular element for plates and shells is developed on the basis of the Hellinger-Reissner mixed variational principle with independently assumed stress and displacement fields. Two main features are emphasized in this development concerning the assumed transverse shear stress field. First, this assumption results in edge-type penalty constraints in the thin plate/shell regimes. Second, the element is used in the form of a four-triangle, cross-diagonal macroelement mesh. As a result, the element is shown to be free from shear locking in thin plate/shell applications, and it satisfies all the appropriate patch tests required for Kirchhoff plate models. In addition, the element exhibits good overall convergence properties in a variety of test problems for plates and shells. Finally, from the computational standpoint, the element is very efficient since all stiffness component matrices are derived in explicit forms.

Saleeb, A. F.↗
NASA NTRS1971

Gravity flow of powder in a lunar environment. Part 2: Analysis of flow initiation

A small displacement-small strain finite element technique utilizing the constant strain triangle and incremental constitutive equations for elasticplastic (media nonhardening and obeying a Coulomb yield condition) was applied to the analysis of gravity flow initiation. This was done in a V-shaped hopper containing a powder under lunar environmental conditions. Three methods of loading were examined. Of the three, the method of computing the initial state of stress in a filled hopper prior to drawdown, by adding material to the hopper layer by layer, was the best. Results of the analysis of a typical hopper problem show that the initial state of stress, the elastic moduli, and the strength parameters have an important influence on material response subsequent to the opening of the hopper outlet.

Pariseau, W. G.↗
NASA NTRS2008

Sub-Surface Windscreen for the Measurement of Outdoor Infrasound

A windscreen has been developed that features two advantages favorable for the measurement of outdoor infrasound. First, the sub-surface location, with the top of the windscreen flush with the ground surface, minimizes the mean velocity of the impinging wind. Secondly, the windscreen material (closed cell polyurethane foam) has a sufficiently low acoustic impedance (222 times that of air) and wall thickness (0.0127 m) to provide a transmission coefficient of nearly unity over the infrasonic frequency range (0-20 Hz). The windscreen, a tightly-sealed box having internal dimensions of 0.3048 x 0.3048 x 0.3556 m, contains a microphone, preamplifier, and a cable feed thru to an external power supply. Provisions are made for rain drainage and seismic isolation. A three-element array, configured as an equilateral triangle with 30.48 m spacing and operating continuously in the field, periodically receives highly coherent signals attributed to emissions from atmospheric turbulence. The time delays between infrasonic signals received at the microphones permit determination of the bearing and elevation of the sources, which correlate well with locations of pilot reports (PIREPS) within a 320 km radius about the array. The test results are interpreted to yield spectral information on infrasonic emissions from clear air turbulence.

Shams, Qamar A.↗
NASA NTRS1990

Shock interference studies on a circular cylinder at Mach 16

This paper presents the results of a computational study on the shock interference problems on a cylindrical body typical of an engine inlet cowl leading edge at a nominal Mach number of 16. The two-dimensional Navier-Stokes equations are solved assuming the flow to be in chemical and thermal equilibrium and using a finite element method. The algorithm employs a cell-centered fully implicit upwind scheme. Adaptively generated unstructured meshes of triangles and quadrilaterals are employed. Under certain conditions the finite element code resulted in oscillatory solutions for shock interference at Mach 16. Some of the causes of the unsteady behavior are identified, and to the extent possible, such situations are avoided in the present application. Two shock interference conditions involving a Type IV (supersonic jet) and a Type III (attaching shear layer) are solved. The results are compared with available experimental data and reasonable agreement is observed. A semi-empirical method is also used to estimate the maximum surface heat flux and static pressure.

Prabhu, Ramadas K.↗
NASA NTRS1977

NASTRAN finite element idealization study

The investigation of the effects of variations of mesh refinement and mesh pattern were conducted using a basic rectangular mesh pattern. When employing the constant strain TRMEM element, the basic rectangular pattern was subdivided into triangles. This subdivision employs two different triangular patterns and allows results to be obtained which demonstrate the effect of modelling bias. Errors in tip deflection, direct stress, and shearing stress as a function of mesh size and element aspect ratio were obtained as well as mid-span stress distributions. All problems were solved on an IBM 360/95 computer using MacNeal-Schwendler Version MSC-38 Rigid Format-1 of the NASTRAN computer program. While NASTRAN uses double precision arithmetic for the solution of the global equations for displacements, subsequent computations to obtain element stresses are carried out in single precision. This suggests that some improvement in stress recovery might be expected when using the higher precision DCD machines.

Case, W. R.↗
NASA NTRS1973

Finite Element Method for Thermal Analysis

A two- and three-dimensional, finite-element thermal-analysis program which handles conduction with internal heat generation, convection, radiation, specified flux, and specified temperature boundary conditions is presented. Elements used in the program are the triangle and tetrahedron for two- and three-dimensional analysis, respectively. The theory used in the program is developed, and several sample problems demonstrating the capability and reliability of the program are presented. A guide to using the program, description of the input cards, and program listing are included.

Heuser, J.↗
NASA NTRS1971

Finite element solution of lubrication problems

A variational formulation of the transient lubrication problem is presented and the corresponding finite element equations derived for three and six point triangles, and, four and eight point quadrilaterals. Test solutions for a one dimensional slider bearing used in validating the computer program are given. Utility of the method is demonstrated by a solution of the shrouded step bearing.

Reddi, M. M.↗
NASA NTRS1973

The finite element method in shell stability analysis.

A development of the finite element method for thin shell instability analysis is presented, covering three principal aspects: (1) representation of shell geometry, (2) representation of element behavior, and (3) algorithmic tools for solution of the large-order systems of nonlinear algebraic equations which characterize various phases of shell instability. Two shell elements are described, an arbitrary quadrilateral and a triangle, and numerical results are presented for two widely-employed comparison problems for linear (stable) analysis. Two shell problems which include instability effects are also solved.-

Gallagher, R. H.↗
NASA NTRS1992

Membrane triangles with corner drilling freedoms. III - Implementation and performance evaluation

This paper completes a three-part series on the formulation of 3-node, 9-dof membrane triangles with corner drilling freedoms based on parametrized variational principles. The first four sections cover element implementation details including determination of optimal parameters and treatment of distributed loads. Then three elements of this type, labeled ALL, FF and EFF-ANDES, are tested on standard plane stress problems. ALL represents numerically integrated versions of Allman's 1988 triangle; FF is based on the free formulation triangle presented by Bergan and Felippa in 1985; and EFF-ANDES represent two different formulations of the optimal triangle derived in Parts I and II. The numerical studies indicate that the ALL, FF and EFF-ANDES elements are comparable in accuracy for elements of unitary aspect ratios. The ALL elements are found to stiffen rapidly in inplane bending for high aspect ratios, whereas the FF and EFF elements maintain accuracy. The EFF and ANDES implementations have a moderate edge in formation speed over the FF.

Felippa, Carlos A.↗
NASA NTRS1980

A two dimensional mesh verification algorithm

A two dimensional mesh verification algorithm is discussed which assumes that a finite element mesh is represented in a program by lists of data, i.e., vertex coordinates, element incidences, boundary data, which describe a collection of triangles. A set of conditions are placed on the triangles (i.e., on the list data) which ensure that the triangles tile a region and which also can be verified by the algorithm.

Simpson, R. B.↗
NASA NTRS1988

Three-dimensional axisymmetric ball model for the PC-TAP program

Spherical-geometry solids' thermal modeling is presently undertaken with a three-dimensional FEM method that can be used in conjunction with the IBM-PC Thermal Analyzer Program (TAP). The elements of the model's nodal network are spherical triangles having approximately equivalent volumes for transient cases, thereby decreasing computer run times. Attention is given to results for the three test cases to which the model has been applied, which are compared with results from analytical closed-form solutions. Absolute temperature errors not greater than 1 percent are indicated by the comparisons.

Pidcoke, L. H.↗
NASA NTRS1992

Spatial and temporal adaptive procedures for the unsteady aerodynamic analysis of airfoils using unstructured meshes

An algorithm which combines spatial and temporal adaption for the time integration of the two dimensional Euler equations on unstructured meshes of triangles is presented. Spatial adaption involves mesh enrichment to add elements in high gradient regions of the flow and mesh coarsening to remove elements where they are no longer needed. Temporal adaption is a time accurate, local time stepping procedure which integrates the flow equations in each cell according to the local numerical stability constraint. The flow solver utilizes a four stage Runge-Kutta time integration scheme with an upwind flux-split spatial discretization. Results obtained using spatial and temporal adaption indicate that highly accurate solutions can be obtained with a significant savings of computing time over global time stepping.

Hooker, John R.↗
NASA NTRS1992

Spatial and temporal adaptive procedures for the unsteady aerodynamic analysis of airfoils using unstructured meshes

An algorithm which combines spatial and temporal adaption for the time integration of the two-dimensional Euler equations on unstructured meshes of triangles is presented. Spatial adaption involves mesh enrichment to add elements in high gradient regions of the flow and mesh coarsening to remove elements where they are no longer needed. Temporal adaption is a time accurate, local time stepping procedure which integrates the flow equations in each cell according to the local numerical stability constraint. The flow solver utilizes a four-stage Runge-Kutta time integration scheme with an upwind flux-split spatial discretization. Results obtained using spatial and temporal adaption indicate that highly accurate solutions can be obtained with a significant savings of computing time over global time stepping.

Hooker, John R.↗