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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Thermostructural design of a carbon-carbon heatshield for a Jovian entry

The thermostructural response of three candidate carbon-carbon composites for the Jovian entry probe heatshield was investigated. The analysis for the three materials, Sandia Felt, Carbitex 700, and SAI 4-D weave carbon-carbon was conducted using a dual finite element approach which involved heat conduction as well as the structural response. A receding boundary due to ablation and inertial loads encountered by the probe were included. Severe cracking, circumferential and radial, and interlaminar shear failure was observed during the radiative heating pulse for the Sandia Felt and Carbitex 700 materials, respectively. The 4-D weave material showed no failures over the entire entry.

White, M. J.↗
NASA NTRS1984

The numerical solution of the Helmholtz equation for wave propagation problems in underwater acoustics

The Helmholtz Equation (-delta-K(2)n(2))u=0 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. A numerical algorithm was developed and a computer code implemented that can effectively solve this equation in the intermediate frequency range. The equation is discretized using the finite element method, thus allowing for the modeling of complicated geometrices (including interfaces) and complicated boundary conditions. A global radiation boundary condition is imposed at the far field boundary that is exact for an arbitrary number of propagating modes. The resulting large, non-selfadjoint system of linear equations with indefinite symmetric part is solved using the preconditioned conjugate gradient method applied to the normal equations. A new preconditioner is developed based on the multigrid method. This preconditioner is vectorizable and is extremely effective over a wide range of frequencies provided the number of grid levels is reduced for large frequencies. A heuristic argument is given that indicates the superior convergence properties of this preconditioner.

Bayliss, A.↗
NASA NTRS1985

The dynamics and control of large flexible space structures, 8

A development of the in plane open loop rotational equations of motion for the proposed Spacecraft Control Laboratory Experiment (SCOLE) in orbit configuration is presented based on an Eulerian formulation. The mast is considered to be a flexible beam connected to the (rigid) shuttle and the reflector. Frequencies and mode shapes are obtained for the mast vibrational appendage modes (assumed to be decoupled) for different boundary conditions based on continuum approaches and also preliminary results are obtained using a finite element representation of the mast reflector system. The linearized rotational in plane equation is characterized by periodic coefficients and open loop system stability can be examined with an application of the Floquet theorem. Numerical results are presented to illustrate the potential instability associated with actuator time delays even for delays which represent only a small fraction of the natural period of oscillation of the modes contained in the open loop model of the system. When plant and measurement noise effects are added to the previously designed deterministic model of the hoop column system, it is seen that both the system transient and steady state performance are degraded. Mission requirements can be satisfied by appropriate assignment of cost function weighting elements and changes in the ratio of plant noise to measurement noise.

Bainum, P. M.↗
NASA NTRS1986

Survey of the status of finite element methods for partial differential equations

The finite element methods (FEM) have proved to be a powerful technique for the solution of boundary value problems associated with partial differential equations of either elliptic, parabolic, or hyperbolic type. They also have a good potential for utilization on parallel computers particularly in relation to the concept of domain decomposition. This report is intended as an introduction to the FEM for the nonspecialist. It contains a survey which is totally nonexhaustive, and it also contains as an illustration, a report on some new results concerning two specific applications, namely a free boundary fluid-structure interaction problem and the Euler equations for inviscid flows.

Temam, Roger↗
NASA NTRS1994

Software For Three-Dimensional Stress And Thermal Analyses

BEST3D is advanced engineering software system for three-dimensional thermal and stress analyses, particularly of components of hot sections of gas-turbine engines. Utilizes boundary element method, offering, in many situations, more accuracy, efficiency, and ease of use than finite element method. Performs engineering analyses of following types: elastic, heat transfer, plastic, forced vibration, free vibration, and transient elastodynamic. Written in FORTRAN 77.

Banerjee, P. K.↗
NASA NTRS1996

A Thermostructural Analysis of a Diboride Composite Leading Edge

In an effort to support the design of zirconium diboride composite leading edges for hypersonic vehicles, a finite element model (FEM) of a prototype leading edge was created and finite element analysis (FEA) was employed to assess its thermal and structural response to aerothermal boundary conditions. Unidirectional material properties for the structural components of the leading edge, a continuous fiber reinforced diboride composite, were computed with COSTAR. These properties agree well with those experimentally measured. To verify the analytical approach taken with COSMOS/M, an independent FEA of one of the leading edge assembly components was also done with COSTAR. Good agreement was obtained between the two codes. Both showed that a unidirectional lay-up had the best margin of safety for a simple loading case. Both located the maximum stress in the same region and ply. The magnitudes agreed within 4 percent. Trajectory based aerothermal heating was then applied to the leading edge assembly FEM created with COSMOS/M to determine steady state temperature response, displacement, stresses, and contact forces due to thermal expansion and thermal strains. Results show that the leading edge stagnation line temperature reached 4700 F. The maximum computed failure index for the laminated composite components peaks at 4.2, and is located at the bolt flange in layer 2 of the side bracket. The temperature gradient in the tip causes a compressive stress of 279 ksi along its width and substantial tensile stresses within its depth.

Kowalski, Tom↗
NASA NTRS2012

A New Approach in Force-Limited Vibration Testing of Flight Hardware

The force-limited vibration test approaches discussed in NASA-7004C were developed to reduce overtesting associated with base shake vibration tests of aerospace hardware where the interface responses are excited coherently. This handbook outlines several different methods of specifying the force limits. The rationale for force limiting is based on the disparity between the impedances of typical aerospace mounting structures and the large impedances of vibration test shakers when the interfaces in general are coherently excited. Among these approaches, the semi-empirical method is presently the most widely used method to derive the force limits. The inclusion of the incoherent excitation of the aerospace structures at mounting interfaces has not been accounted for in the past and provides the basis for more realistic force limits for qualifying the hardware using shaker testing. In this paper current methods for defining the force limiting specifications discussed in the NASA handbook are reviewed using data from a series of acoustic and vibration tests. A new approach based on considering the incoherent excitation of the structural mounting interfaces using acoustic test data is also discussed. It is believed that the new approach provides much more realistic force limits that may further remove conservatism inherent in shaker vibration testing not accounted for by methods discussed in the NASA handbook. A discussion on using FEM/BEM analysis to obtain realistic force limits for flight hardware is provided.

force limiting responses↗
NASA NTRS2023

Acoustics of a Supersonic Mach 1.4 Axisymmetric Spike Inlet

An inlet for a supersonic commercial aircraft is analyzed for noise radiation in low-speed configuration, with auxiliary inlet open. Duct mode amplitudes generated by the fan are determined by post-processing a URANS simulation and become a boundary condition to the inlet. Far-field sound is modeled using commercial acoustics finite element software. The effect of inlet geometry and acoustic liners is simulated and discussed. Absolute levels for inlet tone noise are predicted.

Supersonics↗
DOE OSTI2025

Integrating Crack Detection and Pipe Shape Optimization for Enhanced Sewage System Durability

Crack detection in underground reinforced concrete pipes has been essential in determining the state of stormwater infrastructure. Detection models have been implemented for detecting cracks and other defects in pipes using CCTV footage for stormwater drainage systems. In addition, Finite element models have been used to determine optimum shapes and pipe thickness for different boundary conditions such as header pipes in power plants. The concept of shape optimization emerges as a crucial factor in power plant design and operation, with the potential to maximize performance while minimizing the use of materials. Shape optimization not only enhances efficiency but also contributes to reducing the environmental footprint. This paper discusses the integration of both topics by using the cracks detected in underground pipes as boundary conditions for shape optimization of the pipes. A machine learning model has been developed which uses limited data for training and outlines the location of detected cracks. A shape optimization methodology is proposed in which ANSYS modules are used to analyze fluid flow and then optimize the shape of the pipe. The crack detection model developed has been applied to a crack detected in lab setting and machine learning model used has an accuracy of 98% using a random forest algorithm.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

The Influence of Carbon Fiber Composite Specimen Design Parameters on Artificial Lightning Strike Current Dissipation and Material Thermal Damage

Previous artificial lightning strike direct effect research has examined a broad range of specimen design parameters. No works have studied how such specimen design parameters and electrical boundary conditions impact the dissipation of electric current flow through individual plies. This article assesses the influence of carbon fiber composite specimen design parameters (design parameters = specimen size, shape, and stacking sequence) and electrical boundary conditions on the dissipation of current and the spread of damage resulting from Joule heating. Thermal-electric finite element (FE) modelling is used and laboratory scale (<1 m long) and aircraft scale (>1 m long) models are generated in which laminated ply current dissipation is predicted, considering a fixed artificial lightning current waveform. The simulation results establish a positive correlation between the current exiting the specimen from a given ply and the amount of thermal damage in that ply. The results also establish that the distance to ground, from the strike location to the zero potential boundary conditions (ground), is the controlling factor which dictates the electric current dissipation in each ply. Significantly, this distance to ground is dependent on each of the specimen shape, dimensions, stacking sequence, and location of ground boundary conditions. Therefore, it is not possible to decouple current dissipation and damage from specimen design and boundary condition setup. However, it is possible to define a specimen size for a given specimen shape, stacking sequence, and waveform which limit the influence of specimen dimensions on the resulting current distribution and damage. For a rectangular specimen design which appears in literature multiple times, as 100 × 150 mm and with a stacking sequence of [45/0/-45/90] 4s , a specimen design of greater than 300 × 200 mm is required to limit the influence of specimen dimensions on current distribution and damage.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Boundary identification for 2-D parabolic problems arising in thermal testing of materials

Problems on the identification of two-dimensional spatial domains arising in the detection and characterization of structural flaws in materials are considered. For a thermal diffusion system with external boundary input, observations of the temperature on the surface are used in an output least square approach. Parameter estimation techniques based on the method of mappings are discussed, and approximation schemes are developed based on a finite-element Galerkin approach. Theoretical convergence results for computational techniques are given, and the results are applied to the identification of two kinds of boundary shapes.

Banks, H. T.↗
NASA NTRS1984

Interactive modeling using the PATRAN-G program

PATRAN-G is an interactive finite element generating and graphics display program. The program is currently executable on the VAX or PRIME computers and is compatible with 19 commercially available graphics terminals. The PATRAN program can be used in two modes: generation of a finite element model, or postprocessing of output data from a finite element analysis program. The finite element generating mode is divided into two phases. Phase I is the geometry construction. During this phase, the geometric boundaries of the model are defined and the material properties are input. Mirror imaging, translation and scaling of lines, surfaces, and volumes, and rotation of regions permit the geometry to be described efficiently. Points may be defined in alternate rectangular, cylindrical, or spherical coordinate systems. Also, PATRAN-G can calculate the intersection of lines or surfaces. This capability allows a user to easily model, for example, the intersection of a cylinder which cuts through a sphere at an angle. Once the geometry of the model has been described in terms of points, lines, surfaces, and volumes. Phase II is used to generate the computer finite element model.

Stalnaker, W. A.↗
NASA NTRS1991

An h-p Taylor-Galerkin finite element method for compressible Euler equations

An extension of the familiar Taylor-Galerkin method to arbitrary h-p spatial approximations is proposed. Boundary conditions are analyzed, and a linear stability result for arbitrary meshes is given, showing the unconditional stability for the parameter of implicitness alpha not less than 0.5. The wedge and blunt body problems are solved with both linear, quadratic, and cubic elements and h-adaptivity, showing the feasibility of higher orders of approximation for problems with shocks.

Demkowicz, L.↗
NASA NTRS1996

First-Order System Least-Squares for Second-Order Elliptic Problems with Discontinuous Coefficients

The first-order system least-squares methodology represents an alternative to standard mixed finite element methods. Among its advantages is the fact that the finite element spaces approximating the pressure and flux variables are not restricted by the inf-sup condition and that the least-squares functional itself serves as an appropriate error measure. This paper studies the first-order system least-squares approach for scalar second-order elliptic boundary value problems with discontinuous coefficients. Ellipticity of an appropriately scaled least-squares bilinear form of the size of the jumps in the coefficients leading to adequate finite element approximation results. The occurrence of singularities at interface corners and cross-points is discussed. and a weighted least-squares functional is introduced to handle such cases. Numerical experiments are presented for two test problems to illustrate the performance of this approach.

Manteuffel, Thomas A.↗
NASA NTRS1982

Finite element analysis of the effect of a non-planar solid-liquid interface on the lateral solute segregation during unidirectional solidification

The effect of solid-liquid interface shape on lateral solute segregation during steady-state unidirectional solidification of a binary mixture is calculated under the assumption of no convection in the liquid. A finite element technique is employed to compute the concentration field in the liquid and the lateral segregation in the solid with a curved boundary between the liquid and solid phases. The computational model is constructed assuming knowledge of the solid-liquid interface shape; no attempt is made to relate this shape to the thermal field. The influence of interface curvature on the lateral compositional variation is investigated over a range of system parameters including diffusivity, growth speed, distribution coefficient, and geometric factors of the system. In the limiting case of a slightly nonplanar interface, numerical results from the finite element technique are in good agreement with the analytical solutions of Coriell and Sekerka obtained by using linear theory. For the general case of highly non-planar interface shapes, the linear theory fails and the concentration field in the liquid as well as the lateral solute segregation in the solid can be calculated by using the finite element method.

Carlson, F. M.↗
NASA NTRS1993

An efficient formulation of the coupled finite element-integral equation technique for solving large 3D scattering problems

It is often desirable to calculate the electromagnetic fields inside and about a complicated system of scattering bodies, as well as in their far-field region. The finite element method (FE) is well suited to solving the interior problem, but the domain has to be limited to a manageable size. At the truncation of the FE mesh one can either impose approximate (absorbing) boundary conditions or set up an integral equation (IE) for the fields scattered from the bodies. The latter approach is preferable since it results in higher accuracy. Hence, the two techniques can be successfully combined by introducing a surface that encloses the scatterers, applying a FE model to the inner volume and setting up an IE for the tangential fields components on the surface. Here the continuity of the tangential fields is used bo obtain a consistent solution. A few coupled FE-IE methods have recently appeared in the literature. The approach presented here has the advantage of using edge-based finite elements, a type of finite elements with degrees of freedom associated with edges of the mesh. Because of their properties, they are better suited than the conventional node based elements to represent electromagnetic fields, particularly when inhomogeneous regions are modeled, since the node based elements impose an unnatural continuity of all field components across boundaries of mesh elements. Additionally, our approach is well suited to handle large size problems and lends itself to code parallelization. We will discuss the salient features that make our approach very efficient from the standpoint of numerical computation, and the fields and RCS of a few objects are illustrated as examples.

Cwik, T.↗
DOE OSTI2022

Crystal plasticity modeling of strain-induced martensitic transformations to predict strain rate and temperature sensitive behavior of 304 L steels: Applications to tension, compression, torsion, and impact

This paper advances crystallographically-based Olson-Cohen (direct γ → α’) and deformation mechanism (indirect γ→ε→α’) phase transformation models for predicting strain-induced austenite to martensite transformation. Here, the advanced transformation models enable predictions of not only strain-path sensitive, but also of strain-rate and temperature sensitive deformation of polycrystalline stainless steels (SSs). The deformation of constituent grains in SSs is modeled as a combination of anisotropic elasticity, crystallographic slip, and phase transformation, while the hardening is based on the evolution of dislocation density and explicit shifts in phase fractions. Such grain-scale deformation is implemented within the meso-scale elasto-plastic self-consistent (EPSC) homogenization model, which is coupled with the implicit finite element (FE) method to provide a constitutive response at each FE integration point for solving boundary value problems at the macro-scale. Parameters pertaining to the hardening and transformation models within FEEPSC are calibrated and validated on a suite of data including flow curves and phase fractions for monotonic compression, tension, and torsion as a function of strain-rate and temperature for wrought and additively manufactured (AM) SS304L. To illustrate the potential and accuracy of the integrated multi-level FE-EPSC simulation framework, geometry, mechanical response, phase fractions, and texture evolution are simulated during gas-gun impact deformation of a cylinder and quasi-static tension of a notched specimen made of AM SS304L. Details of the simulation framework, comparison between experimental and simulation results, and insights from the results are presented and discussed.

304L steels↗
DOE OSTI2024

Alloying effects on deformation induced microstructure evolution in copper

In this work, we investigated the effects of alloying elements on plastic deformation and microstructure evolution in polycrystalline copper (Cu) and Cu alloyed with 1 wt. % lead (Cu-1 %Pb). These materials were selected due to the size mismatch between Cu and Pb, with the latter forming precipitates at grain boundaries. Multi-modal characterization techniques, including neutron diffraction, electron backscatter diffraction (EBSD), and transmission electron microscopy (TEM), along with finite element simulations were employed to study the deformation behavior across multiple length scales. While both Cu and Cu-1%Pb exhibited similar macroscale response and final deformation textures, both dislocation line profile analysis and TEM revealed increased dislocation density in deformed Cu-1%Pb specimens. The presence of lead precipitates also significantly affected local plastic deformation during compression, with their influence diminishing with increasing strain. These results demonstrate the complex relationships between alloying elements, plastic deformation, microstructural evolution, and material behavior under load. The insights gained from this multi-scale and multi-technique approach contribute to the fundamental understanding of microstructural evolution in immiscible alloys and are valuable for tailoring the properties of structural materials for specific engineering applications.

36 MATERIALS SCIENCE↗