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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Derivation of generalized transition/boundary conditions for planar multiple-layer structures

Infinite-order generalized impedance boundary conditions (GIBCs) and generalized sheet transition conditions (GSTCs) for planar multilayer configurations are developed via the Taylor series expansion method. The conditions are derived in a matrix product form where each matrix corresponds to a specific layer. An overall composite boundary/transition condition is obtained by making finite-order approximations to the elements of each matrix for the cases of 'low-contrast' and 'high-contrast' material layers. The accuracy of the truncated boundary conditions is examined by comparing their implied reflection and transmission coefficients with the corresponding exact coefficients. Design curves are also given which relate the maximum order of the conditions required to simulate a coating or layer of specific thickness and contrast. Expressions are then derived for the reflection and transmission coefficients of the GIBC/GSTC sheets, and these are compared to exact coefficients to demonstrate the validity of the derived GIBCs/GSTCs.

Ricoy, M. A.↗
NASA NTRS1983

Progress on a generalized coordinates tensor product finite element 3DPNS algorithm for subsonic

A generalized coordinates form of the penalty finite element algorithm for the 3-dimensional parabolic Navier-Stokes equations for turbulent subsonic flows was derived. This algorithm formulation requires only three distinct hypermatrices and is applicable using any boundary fitted coordinate transformation procedure. The tensor matrix product approximation to the Jacobian of the Newton linear algebra matrix statement was also derived. Tne Newton algorithm was restructured to replace large sparse matrix solution procedures with grid sweeping using alpha-block tridiagonal matrices, where alpha equals the number of dependent variables. Numerical experiments were conducted and the resultant data gives guidance on potentially preferred tensor product constructions for the penalty finite element 3DPNS algorithm.

Baker, A. J.↗
NASA NTRS2010

Hybrid FEA/SEA Assessment for an Orthogrid Cylindrical Panel Section and Periodic Subsystem Modeling Evaluation

In the lower frequency range, where particular boundary conditions can make a significant difference to panel response characteristics Statistical Energy Analysis (SEA) has never been the analytical tool of choice. In addition to boundary condition effects, SEA is not well suited in frequency bands where no modes or less than a few modes exist. The advent of the Hybrid Module has enabled integration of Finite Element Analysis to expand and enhance the capability for response calculations within VA One into the lower frequency range. Exploration of several additional modeling approaches was completed for the cylindrical orthogrid panel test article that was examined in Reference 1. Comparison of the new analytical response predictions with the measured response data from ground test and the pure SEA results from the reference will be presented. One approach that is considered promising is the periodic subsystem capability. Initially, a detailed FEM of just one region of the test article is defined. After evaluating this small region using symmetric boundary conditions, the FEM may be expanded to determine the properties of the entire system using similar connected regions that map over the entire test article. Another approach is the direct use of a very detailed finite element model of the entire panel, explicitly modeling pocket and rib details of the structure. A third approach is to approximate localized structure geometry details with a smeared property generalization using a PCOMP (NASTRAN card used to define layered composite structures) to define skin layer and ribbed layer for the orthogrid panel. The authors expect to demonstrate that the integrated Hybrid/FEM approach increases confidence in response prediction in the lower frequency range (for example from 20-300 Hz for the test article under consideration). In addition the strength and weakness of each additional approach will be highlighted and compared to those reported with those reported in an earlier paper

Smith, Andrew M.↗
NASA NTRS1976

A study of turbulent flow between parallel plates by a statistical method

Turbulent Couette flow between parallel plates was studied from a statistical mechanics approach utilizing a model equation, similar to the Boltzmann equation of kinetic theory, which was proposed by Lundgren from the velocity distribution of fluid elements. Solutions to this equation are obtained numerically, employing the discrete ordinate method and finite differences. Two types of boundary conditions on the distribution function are considered, and the results of the calculations are compared to available experimental data. The research establishes that Lundgren's equation provides a very good description of turbulence for the flow situation considered and that it offers an analytical tool for further study of more complex turbulent flows. The present work also indicates that modelling of the boundary conditions is an area where further study is required.

Srinivasan, R.↗
NASA NTRS1981

An adaptive, multi-level method for elliptic boundary value problems

Subroutine PLTMG is a FORTRAN program for solving self-adjoint elliptic boundary value problems in general regions of R-squared. It is based on a piecewise linear triangle finite element method, an adaptive grid refinement procedure, and a multi-level iterative method to solve the resulting sets of linear equations. This paper describes the method and presents some numerical results and comparisons.

Bank, R. E.↗
NASA NTRS1982

Noise transmission into semicylindrical enclosures through discretely stiffened curved panels

An analytical study of sound transmission into semicylindrical enclosures through discretely stiffened curved elastic panels is presented. The transmitted sound is estimated by solving the acoustic wave equation for the interior acoustic field, a Galerkin-like method being used. This solution is then coupled to the vibration of the stiffened panels. The response characteristics of these panels are determined by using a modal analysis where the modes are obtained by the finite element-strip method. Numerical results include spectra of the interior sound pressure due to white noise, turbulent boundary layer and propeller noise inputs.

Chang, M. T.↗
NASA NTRS1984

On accuracy conditions for the numerical computation of waves

The Helmholtz equation (Delta + K(2)n(2))u = f with a variable index of refraction n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. Such problems can be solved numerically by first truncating the given unbounded domain and imposing a suitable outgoing radiation condition on an artificial boundary and then solving the resulting problem on the bounded domain by direct discretization (for example, using a finite element method). In practical applications, the mesh size h and the wave number K, are not independent but are constrained by the accuracy of the desired computation. It will be shown that the number of points per wavelength, measured by (Kh)(-1), is not sufficient to determine the accuracy of a given discretization. For example, the quantity K(3)h(2) is shown to determine the accuracy in the L(2) norm for a second-order discretization method applied to several propagation models.

Bayliss, A.↗
NASA NTRS1985

On accuracy conditions for the numerical computation of waves

The Helmholtz equation (Delta + K(2)n(2))u = f with a variable index of refraction n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. Such problems can be solved numerically by first truncating the given unbounded domain and imposing a suitable outgoing radiation condition on an artificial boundary and then solving the resulting problem on the bounded domain by direct discretization (for example, using a finite element method). In practical applications, the mesh size h and the wave number K, are not independent but are constrained by the accuracy of the desired computation. It will be shown that the number of points per wavelength, measured by (Kh)(-1), is not sufficient to determine the accuracy of a given discretization. For example, the quantity K(3)h(2) is shown to determine the accuracy in the L(2) norm for a second-order discretization method applied to several propagation models.

Bayliss, A.↗
NASA NTRS1987

Stress concentration around a part-through hole weakening a laminated plates

A quasi-three-dimensional-elasticity type theory, based on the assumptions of transverse inextensibility and layerwise constant shear-angle, is presented, in the context of an assumed quadratic displacement potential energy approach, for analyzing an edge-loaded laminated plate weakened by the presence of a part-through hole. Numerical results, obtained using the C-deg-type triangular finite element, indicate the existence of severe cross-sectional warping in the vicinity of the hole and plate boundaries. Furthermore, the three-dimensional nature of the stress concentration factor in the neighborhood of the hole boundary is clearly exhibited. Besides, very high stress concentration factors are found in the layer weakened by the part-through hole. The numerical results presented should serve as bench-mark solutions for future comparisons.

Chaudhuri, Reaz A.↗
NASA NTRS1989

Probalistic structural analysis methods for select space propulsion system structural components

A summary of the status of this five-year project which is now in its third year of research and development is presented. The goal of the project is the development of several methodologies for probabilistic structural modeling. Probabilistic structural modeling consists of stochastic models of material properties, part geometries, boundary conditions, as well as loading conditions. The current presentation focuses on one methodology - coupling of an advanced finite element structural analysis code with probabilistic modeling strategies. The essential algorithm developments for combining the finite element and probabilistic analysis methods are reported. The validity of the resulting probabilistic structural analysis method is confirmed through a series of test problems with exact results based on Monte Carlo simulations. Additionally, the applicability of the method to a Space Propulsion System (a turbine blade) is demonstrated for static stresses.

Cruse, T. A.↗
NASA NTRS1988

Efficient vibration mode analysis of aircraft with multiple external store configurations

A coupling method for efficient vibration mode analysis of aircraft with multiple external store configurations is presented. A set of low-frequency vibration modes, including rigid-body modes, represent the aircraft. Each external store is represented by its vibration modes with clamped boundary conditions, and by its rigid-body inertial properties. The aircraft modes are obtained from a finite-element model loaded by dummy rigid external stores with fictitious masses. The coupling procedure unloads the dummy stores and loads the actual stores instead. The analytical development is presented, the effects of the fictitious mass magnitudes are discussed, and a numerical example is given for a combat aircraft with external wing stores. Comparison with vibration modes obtained by a direct (full-size) eigensolution shows very accurate coupling results. Once the aircraft and stores data bases are constructed, the computer time for analyzing any external store configuration is two to three orders of magnitude less than that of a direct solution.

Karpel, M.↗
NASA NTRS1990

Coupled BE/FE/BE approach for scattering from fluid-filled structures

NASHUA is a coupled finite element/boundary element capability built around NASTRAN for calculating the low frequency far-field acoustic pressure field radiated or scattered by an arbitrary, submerged, three-dimensional, elastic structure subjected to either internal time-harmonic mechanical loads or external time-harmonic incident loadings. Described here are the formulation and use of NASHUA for solving such structural acoustics problems when the structure is fluid-filled. NASTRAN is used to generate the structural finite element model and to perform most of the required matrix operations. Both fluid domains are modeled using the boundary element capability in NASHUA, whose matrix formulation (and the associated NASTRAN DMAP) for evacuated structures can be used with suitable interpretation of the matrix definitions. After computing surface pressures and normal velocities, far-field pressures are evaluated using an asymptotic form of the Helmholtz exterior integral equation. The proposed numerical approach is validated by comparing the acoustic field scattered from a submerged fluid-filled spherical thin shell to that obtained with a series solution, which is also derived here.

Everstine, Gordon C.↗
NASA NTRS1991

Compressible convection in a viscous Venusian mantle

Finite element simulations of axisymmetric spherical shell compressible convection were carried out to investigate the effect of various surface boundary conditions in a Venusian mantle, employing a thermal expansivity alpha which decreased with depth, a uniform viscosity an order of magnitude greater than the earth's, and zero and chondric quantities of internal heating. As long as hot plumes from the core-mantle boundary were strong, the convection pattern was typical of that for variable alpha flow; that is, it was characterized by steady upflowing regions, unsteady collections of downflowing plumes, and large aspect ratio cells. Increases in the internal heating or the temperature T0 at the top of the convecting layer weakened the hot plumes and, therefore, decreased the width of the cells. A rigid surface increased the internal temperature and also decreased the width of convection cells. Extensive regions of subadiabaticity were found in the mantle. These results are compared with those for fully three-dimensional convection under similar conditions (Schubert et al., 1990).

Leitch, A. M.↗
NASA NTRS1993

A nonconforming multigrid method using conforming subspaces

For second-order elliptic boundary value problems, we develop a nonconforming multigrid method using the coarser-grid correction on the conforming finite element subspaces. The convergence proof with an arbitrary number of smoothing steps for nu-cycle is presented.

Lee, Chang Ock↗
NASA NTRS1994

Compositional vs. thermal buoyancy and the evolution of subducted lithosphere

We formulate 2-D Cartesian finite element models that explore the fate of compositionally defined lithosphere as it encounters a viscosity increase at the boundary between the upper and lower mantle. Subducted lithosphere is represented as a cold, stiff, layered composite of denser eclogite underlain by more buoyant harzburgite. Slabs impinging on a lower mantle 30 and 100 times more viscous than the upper mantle thicken and fold strongly as they penetrate the lower mantle. Approximately a factor of two thickening occurs via pure shear just above the discontinuity, with additional enhancement due to folding by over a factor of two. No separation of the individual slab components occurs at the discontinuity, and direct comparison with models in which compositional buoyancy is explicitly ignored indicates that slab evolution is largely controlled by the thermal buoyancy. These results are at odds with hypotheses about slab evolution in which the compositional buoyancy contributions lead to component separation and the formation of slab megaliths or a compositionally layered upper mantle.

Gaherty, James B.↗
NASA NTRS1996

Development of a Higher Order Laminate Theory for Modeling Composites with Induced Strain Actuators

A refined higher order plate theory is developed to investigate the actuation mechanism of piezoelectric materials surface bonded or embedded in composite laminates. The current analysis uses a displacement field which accurately accounts for transverse shear stresses. Some higher order terms are identified by using the conditions that shear stresses vanish at all free surfaces. Therefore, all boundary conditions for displacements and stresses are satisfied in the present theory. The analysis is implemented using the finite element method which provides a convenient means to construct a numerical solution due to the discrete nature of the actuators. The higher order theory is computationally less expensive than a full three dimensional analysis. The theory is also shown to agree well with published experimental results. Numerical examples are presented for composite plates with thicknesses ranging from thin to very thick.

Chattopadhyay, Aditi↗
NASA NTRS2001

Thermal Analysis of the MC1 Engine Turbopump

The MC1 Engine turbopump supplied the propellants to the main injector. The turbopump consisted of four parts; lox pump, interpropellant seal package (IPS), RP pump and turbine. The thermal analysis was divided into two 2D finite element models; Housing or stationary parts and rotor or rotating parts. Both models were analyzed at the same boundary conditions using SINDA. The housing model consisted of, lox pump housing, ips housing, RP housing, turbine inlet housing, turbine housing, exit guide vane, heat shield and both bearing outer races. The rotor model consisted of the lox impeller; lox end bearing and id race, RP impeller, and RP bearing and id race, shaft and turbine disk. The objectives of the analysis were to: (1) verified the original design and recommend modifications to it, (2) submitted a thermal environment to support the structural analysis, (3) support the component and engine test program. and (4) to support the X34 vehicle program.

Roman, Jose↗
NASA NTRS2002

Thermal Analysis of the MCI Engine Turbopump

The MCI Engine turbopump supplied the propellants to the main injector. The turbopump consisted of four parts; lox pump, interpropellant seal package (IPS), RP pump and turbine. The thermal analysis was divided into two 2D finite element models; Housing or stationary parts and rotor or rotating parts. Both models were analyzed at the same boundary conditions using SINDA. The housing model consisted of; lox pump housing, ips housing, RP housing, turbine inlet housing, turbine housing, exit guide vane, heat shield and both bearing outer races. The rotor model consisted of the lox impeller; lox end bearing and id race, RP impeller, and RP bearing and id race, shaft and turbine disk. The objectives of the analysis were to (1) verified the original design and recommend modifications to it, (2) submitted a thermal environment to support the structural analysis, (3) support the component and engine test program and (4) to support the X34 vehicle program.

Roman, Jose↗