Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boundary 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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

On the Effects of Surface Roughness on Boundary Layer Transition

Surface roughness can influence laminar-turbulent transition in many different ways. This paper outlines selected analyses performed at the NASA Langley Research Center, ranging in speed from subsonic to hypersonic Mach numbers and highlighting the beneficial as well as adverse roles of the surface roughness in technological applications. The first theme pertains to boundary-layer tripping on the forebody of a hypersonic airbreathing configuration via a spanwise periodic array of trip elements, with the goal of understanding the physical mechanisms underlying roughness-induced transition in a high-speed boundary layer. The effect of an isolated, finite amplitude roughness element on a supersonic boundary layer is considered next. The other set of flow configurations examined herein corresponds to roughness based laminar flow control in subsonic and supersonic swept wing boundary layers. A common theme to all of the above configurations is the need to apply higher fidelity, physics based techniques to develop reliable predictions of roughness effects on laminar-turbulent transition.

Choudhari, Meelan M.↗
NASA NTRS2001

SGBEM-FEM Alternating Method for Analyzing 3D Non-Planar Cracks and Their Growth in Structural Components

The highly accurate and efficient Symmetric Galerkin Boundary Element Method (SGBEM), a Finite Element Method (FEM)-based alternating method, is proposed for analyzing three-dimensional non-planar cracks and their growth. The cracks are modeled using the symmetric Galerkin boundary element method as a distribution of displacement discontinuities, simulating an infinite medium. The finite element method only analyzes the stress for the uncracked body. The solution for the cracked structural component is determined by an iteration procedure. This process alternates between an FEM solution for the uncracked body and the SGBEM solution for a crack in an infinite body. Numerical analysis, and the Java code used, evaluate stress intensity factors and model fatigue crack growth. Examples of non-planar cracks in infinite media and planar cracks in finite bodies, as well as growth under fatigue, show the accuracy of the method.

Nikishkov, G. P.↗
NASA NTRS1984

ESCHER: An interactive mesh-generating editor for preparing finite-element input

ESCHER is an interactive mesh generation and editing program designed to help the user create a finite-element mesh, create additional input for finite-element analysis, including initial conditions, boundary conditions, and slidelines, and generate a NEUTRAL FILE that can be postprocessed for input into several finite-element codes, including ADINA, ADINAT, DYNA, NIKE, TSAAS, and ABUQUS. Two important ESCHER capabilities, interactive geometry creation and mesh archival storge are described in detail. Also described is the interactive command language and the use of interactive graphics. The archival storage and restart file is a modular, entity-based mesh data file. Modules of this file correspond to separate editing modes in the mesh editor, with data definition syntax preserved between the interactive commands and the archival storage file. Because ESCHER was expected to be highly interactive, extensive user documentation was provided in the form of an interactive HELP package.

Oakes, W. R., Jr.↗
DOE OSTI2020

The second-generation Shifted Boundary Method and its numerical analysis

Recently, the Shifted Boundary Method (SBM) was proposed within the class of unfitted (or immersed, or embedded) finite element methods. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we present enhanced variational SBM formulations for the Poisson and Stokes problems with improved flexibility and robustness. These simplified variational forms allow to relax some of the assumptions required by the mathematical proofs of stability and convergence of earlier implementations. First, we show that these new SBM implementations can be proved asymptotically stable and convergent even without the rather restrictive assumption that the inner product between the normals to the true and surrogate boundaries is positive. Second, we show that it is not necessary to introduce a stabilization term involving the tangential derivatives of the solution at Dirichlet boundaries, therefore avoiding the calibration of an additional stabilization parameter. Finally, we prove enhanced L 2 -estimates without the cumbersome assumption – of earlier proofs – that the surrogate domain is convex. Instead we rely on a conventional assumption that the boundary of the true domain is smooth, which can also be replaced by requiring convexity of the true domain. The aforementioned improvements open the way to a more general and efficient implementation of the Shifted Boundary Method, particularly in complex three-dimensional geometries. We complement these theoretical developments with numerical experiments in two and three dimensions.

42 ENGINEERING↗
NASA NTRS1985

Three-dimensional elastic-plastic finite-element analysis of fatigue crack propagation

Fatigue cracks are a major problem in designing structures subjected to cyclic loading. Cracks frequently occur in structures such as aircraft and spacecraft. The inspection intervals of many aircraft structures are based on crack-propagation lives. Therefore, improved prediction of propagation lives under flight-load conditions (variable-amplitude loading) are needed to provide more realistic design criteria for these structures. The main thrust was to develop a three-dimensional, nonlinear, elastic-plastic, finite element program capable of extending a crack and changing boundary conditions for the model under consideration. The finite-element model is composed of 8-noded (linear-strain) isoparametric elements. In the analysis, the material is assumed to be elastic-perfectly plastic. The cycle stress-strain curve for the material is shown Zienkiewicz's initial-stress method, von Mises's yield criterion, and Drucker's normality condition under small-strain assumptions are used to account for plasticity. The three-dimensional analysis is capable of extending the crack and changing boundary conditions under cyclic loading.

Goglia, G. L.↗
NASA NTRS1992

A boundary element alternating method for two-dimensional mixed-mode fracture problems

A boundary element alternating method, denoted herein as BEAM, is presented for two dimensional fracture problems. This is an iterative method which alternates between two solutions. An analytical solution for arbitrary polynomial normal and tangential pressure distributions applied to the crack faces of an embedded crack in an infinite plate is used as the fundamental solution in the alternating method. A boundary element method for an uncracked finite plate is the second solution. For problems of edge cracks a technique of utilizing finite elements with BEAM is presented to overcome the inherent singularity in boundary element stress calculation near the boundaries. Several computational aspects that make the algorithm efficient are presented. Finally, the BEAM is applied to a variety of two dimensional crack problems with different configurations and loadings to assess the validity of the method. The method gives accurate stress intensity factors with minimal computing effort.

Raju, I. S.↗
NASA NTRS2012

An Integrated Biomechanical Model for Microgravity-Induced Visual Impairment

When gravitational unloading occurs upon entry to space, astronauts experience a major shift in the distribution of their bodily fluids, with a net headward movement. Measurements have shown that intraocular pressure spikes, and there is a strong suspicion that intracranial pressure also rises. Some astronauts in both short- and long-duration spaceflight develop visual acuity changes, which may or may not reverse upon return to earth gravity. To date, of the 36 U.S. astronauts who have participated in long-duration space missions on the International Space Station, 15 crew members have developed minor to severe visual decrements and anatomical changes. These ophthalmic changes include hyperopic shift, optic nerve distension, optic disc edema, globe flattening, choroidal folds, and elevated cerebrospinal fluid pressure. In order to understand the physical mechanisms behind these phenomena, NASA is developing an integrated model that appropriately captures whole-body fluids transport through lumped-parameter models for the cerebrospinal and cardiovascular systems. This data feeds into a finite element model for the ocular globe and retrobulbar subarachnoid space through time-dependent boundary conditions. Although tissue models and finite element representations of the corneo-scleral shell, retina, choroid and optic nerve head have been integrated to study pathological conditions such as glaucoma, the retrobulbar subarachnoid space behind the eye has received much less attention. This presentation will describe the development and scientific foundation of our holistic model.

Nelson, Emily S.↗
NASA NTRS1993

A comparative study of an ABC and an artificial absorber for truncating finite element meshes

The type of mesh termination used in the context of finite element formulations plays a major role on the efficiency and accuracy of the field solution. The performance of an absorbing boundary condition (ABC) and an artificial absorber (a new concept) for terminating the finite element mesh was evaluated. This analysis is done in connection with the problem of scattering by a finite slot array in a thick ground plane. The two approximate mesh truncation schemes are compared with the exact finite element-boundary integral (FEM-BI) method in terms of accuracy and efficiency. It is demonstrated that both approximate truncation schemes yield reasonably accurate results even when the mesh is extended only 0.3 wavelengths away from the array aperture. However, the artificial absorber termination method leads to a substantially more efficient solution. Moreover, it is shown that the FEM-BI method remains quite competitive with the FEM-artificial absorber method when the FFT is used for computing the matrix-vector products in the iterative solution algorithm. These conclusions are indeed surprising and of major importance in electromagnetic simulations based on the finite element method.

Oezdemir, T.↗
NASA NTRS1991

A locally refined rectangular grid finite element method - Application to computational fluid dynamics and computational physics

The present FEM technique addresses both linear and nonlinear boundary value problems encountered in computational physics by handling general three-dimensional regions, boundary conditions, and material properties. The box finite elements used are defined by a Cartesian grid independent of the boundary definition, and local refinements proceed by dividing a given box element into eight subelements. Discretization employs trilinear approximations on the box elements; special element stiffness matrices are included for boxes cut by any boundary surface. Illustrative results are presented for representative aerodynamics problems involving up to 400,000 elements.

Young, David P.↗
NASA NTRS2018

Cyclic Symmetry Finite Element Forced Response Analysis of a Distortion-Tolerant Fan with Boundary Layer Ingestion

Accurate prediction of the blade vibration stress is required to determine overall durability of fan blade design under Boundary Layer Ingestion (BLI) distorted flow environments. Traditional single blade modeling technique is incapable of representing accurate modeling for the entire rotor blade system subject to complex dynamic loading behaviors and vibrations in distorted flow conditions. A particular objective of our work was to develop a high-fidelity full-rotor aeromechanics analysis capability for a system subjected to a distorted inlet flow by applying cyclic symmetry finite element modeling methodology. This reduction modeling method allows computationally very efficient analysis using a small periodic section of the full rotor blade system. Experimental testing by the use of the 8-foot by 6-foot Supersonic Wind Tunnel Test facility at NASA Glenn Research Center was also carried out for the system designated as the Boundary Layer Ingesting Inlet/Distortion-Tolerant Fan (BLI2DTF) technology development. The results obtained from the present numerical modeling technique were evaluated with those of the wind tunnel experimental test, toward establishing a computationally efficient aeromechanics analysis modeling tool facilitating for analyses of the full rotor blade systems subjected to a distorted inlet flow conditions. Fairly good correlations were achieved hence our computational modeling techniques were fully demonstrated. The analysis result showed that the safety margin requirement set in the BLI2DTF fan blade design provided a sufficient margin with respect to the operating speed range.

forced response analysis↗
NASA NTRS1976

The Relation of Finite Element and Finite Difference Methods

Finite element and finite difference methods are examined in order to bring out their relationship. It is shown that both methods use two types of discrete representations of continuous functions. They differ in that finite difference methods emphasize the discretization of independent variable, while finite element methods emphasize the discretization of dependent variable (referred to as functional approximations). An important point is that finite element methods use global piecewise functional approximations, while finite difference methods normally use local functional approximations. A general conclusion is that finite element methods are best designed to handle complex boundaries, while finite difference methods are superior for complex equations. It is also shown that finite volume difference methods possess many of the advantages attributed to finite element methods.

Vinokur, M.↗
NASA NTRS1985

Analysis of a finite composite plate with smooth rigid pin

An elastic continuum method of analysis developed for a smooth rigid pin in a finite composite plate by using complex potentials that satisfy the governing orthotropic differential equation is presented. The plate is subjected to uniaxial loading and the nonlinear problems of partial contact in the pin-plate are studied. Three types of pin fits are examined: interference, push, and clearance. The problems of stresses due to pure interference, load contact behavior, and stresses variations are addressed. The effects of the rectangular geometry and the orthotropic property of the plate are demonstrated through numerical results. It is found that the use of large hole diameter/plate width ratios have deleterious effects and the use of interference has beneficial effects. The solutions obtained will help in the development of special finite elements for fastener joints. Finally, the boundary conditions on the finite outer boundaries are found to be satisfied by a successive integration technique.

Naidu, A. C. B.↗
NASA NTRS1995

Thermal behavior spiral bevel gears

An experimental and analytical study of the thermal behavior of spiral bevel gears is presented. Experimental data were taken using thermocoupled test hardware and an infrared microscope. Many operational parameters were varied to investigate their effects on the thermal behavior. The data taken were also used to validate the boundary conditions applied to the analytical model. A finite element-based solution sequence was developed. The three-dimensional model was developed based on the manufacturing process for these gears. Contact between the meshing gears was found using tooth contact analysis to describe the location, curvatures, orientations, and surface velocities. This information was then used in a three-dimensional Hertzian contact analysis to predict contact ellipse size and maximum pressure. From these results, an estimate of the heat flux magnitude and the location on the finite element model was made. The finite element model used time-averaged boundary conditions to permit the solution to attain steady state in a computationally efficient manner.Then time- and position-varying boundary conditions were applied to the model to analyze the cyclic heating and cooling due to the gears meshing and transferring heat to the surroundings, respectively. The model was run in this mode until the temperature behavior stabilized. The transient flash temperature on the surface was therefore described. The analysis can be used to predict the overall expected thermal behavior of spiral bevel gears. The experimental and analytical results were compared for this study and also with a limited number of other studies. The experimental and analytical results attained in the current study were basically within 10% of each other for the cases compared. The experimental comparison was for bulk thermocouple locations and data taken with an infrared microscope. The results of a limited number of other studies were compared with those obtained herein and predicted the same basic behavior.

Handschuh, Robert F.↗
NASA NTRS1992

Development of 3D electromagnetic modeling tools for airborne vehicles

The main goal of this project is to develop methodologies for scattering by airborne composite vehicles. Although our primary focus continues to be the development of a general purpose code for analyzing the entire structure as a single unit, a number of other tasks are also pursued in parallel with this effort. These tasks are important in testing the overall approach and in developing suitable models for materials coatings, junctions and, more generally, in assessing the effectiveness of the various parts comprising the final code. Here, we briefly discuss our progress on the five different tasks which were pursued during this period. Our progress on each of these tasks is described in the detailed reports (listed at the end of this report) and the memoranda included. The first task described below is, of course, the core of this project and deals with the development of the overall code. Undoubtedly, it is the outcome of the research which was funded by NASA-Ames and the Navy over the past three years. During this year we developed the first finite element code for scattering by structures of arbitrary shape and composition. The code employs a new absorbing boundary condition which allows termination of the finite element mesh only 0.3 lambda from the outer surface of the target. This leads to a remarkable reduction of the mesh size and is a unique feature of the code. Other unique features of this code include capabilities to model resistive sheets, impedance sheets and anisotropic materials. This last capability is the latest feature of the code and is still under development. The code has been extensively validated for a number of composite geometries and some examples are given. The validation of the code is still in progress for anisotropic and larger non-metallic geometries and cavities. The developed finite element code is based on a Galerkin's formulation and employs edge-based tetrahedral elements for discretizing the dielectric sections and the region between the target and the outer mesh termination boundary (ATB). This boundary is placed in conformity with the target's outer surface, thus resulting in additional reduction of the unknown count.

Volakis, John L.↗
NASA NTRS2012

Vibration Propagation of Gear Dynamics in a Gear-Bearing-Housing System Using Mathematical Modeling and Finite Element Analysis

Vibration and noise caused by gear dynamics at the meshing teeth propagate through power transmission components to the surrounding environment. This study is devoted to developing computational tools to investigate the vibro-acoustic propagation of gear dynamics through a gearbox using different bearings. Detailed finite element/contact mechanics and boundary element models of the gear/bearing/housing system are established to compute the system vibration and noise propagation. Both vibration and acoustic models are validated by experiments including the vibration modal testing and sound field measurements. The effectiveness of each bearing type to disrupt vibration propagation is speed-dependent. Housing plays an important role in noise radiation .It, however, has limited effects on gear dynamics. Bearings are critical components in drivetrains. Accurate modeling of rolling element bearings is essential to assess vibration and noise of drivetrain systems. This study also seeks to fully describe the vibro-acoustic propagation of gear dynamics through a power-transmission system using rolling element and fluid film wave bearings. Fluid film wave bearings, which have higher damping than rolling element bearings, could offer an energy dissipation mechanism that reduces the gearbox noise. The effectiveness of each bearing type to disrupt vibration propagation in explored using multi-body computational models. These models include gears, shafts, rolling element and fluid film wave bearings, and the housing. Radiated noise is mapped from the gearbox surface to surrounding environment. The effectiveness of rolling element and fluid film wave bearings in breaking the vibro-acoustic propagation path from the gear to the housing is investigated.

Parker, Robert G.↗
NASA NTRS1990

Task reports on developing techniques for scattering by 3D composite structures and to generate new solutions in diffraction theory using higher order boundary conditions

There are two tasks described in this report. First, an extension of a two dimensional formulation is presented for a three dimensional body of revolution. With the introduction of a Fourier expansion of the vector electric and magnetic fields, a coupled two dimensional system is generated and solved via the finite element method. An exact boundary condition is employed to terminate the mesh and the fast fourier transformation is used to evaluate the boundary integrals for low O(n) memory demand when an iterative solution algorithm is used. Second, the diffraction by a material discontinuity in a thick dielectric/ferrite layer is considered by modeling the layer as a distributed current sheet obeying generalized sheet transition conditions (GSTC's).

Volakis, John L.↗
NASA NTRS1974

COMOC: Three dimensional boundary region variant, programmer's manual

The three-dimensional boundary region variant of the COMOC computer program system solves the partial differential equation system governing certain three-dimensional flows of a viscous, heat conducting, multiple-species, compressible fluid including combustion. The solution is established in physical variables, using a finite element algorithm for the boundary value portion of the problem description in combination with an explicit marching technique for the initial value character. The computational lattice may be arbitrarily nonregular, and boundary condition constraints are readily applied. The theoretical foundation of the algorithm, a detailed description on the construction and operation of the program, and instructions on utilization of the many features of the code are presented.

Orzechowski, J. A.↗
NASA NTRS1989

Probabilistic structural analysis methods for select space propulsion system components

The Probabilistic Structural Analysis Methods (PSAM) project developed at the Southwest Research Institute integrates state-of-the-art structural analysis techniques with probability theory for the design and analysis of complex large-scale engineering structures. An advanced efficient software system (NESSUS) capable of performing complex probabilistic analysis has been developed. NESSUS contains a number of software components to perform probabilistic analysis of structures. These components include: an expert system, a probabilistic finite element code, a probabilistic boundary element code and a fast probability integrator. The NESSUS software system is shown. An expert system is included to capture and utilize PSAM knowledge and experience. NESSUS/EXPERT is an interactive menu-driven expert system that provides information to assist in the use of the probabilistic finite element code NESSUS/FEM and the fast probability integrator (FPI). The expert system menu structure is summarized. The NESSUS system contains a state-of-the-art nonlinear probabilistic finite element code, NESSUS/FEM, to determine the structural response and sensitivities. A broad range of analysis capabilities and an extensive element library is present.

Millwater, H. R.↗