Search NASA⌕ Search

SEARCH · Search NASA

Results for “extended finite element method”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

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.↗
NASA NTRS2013

Global Response of the Space Shuttle External Tank with the Presence of Intertank Stringer Cracks and Radius Blocks

After propellant was loaded into the external tank (ET), the November 5, 2010 launch of Space Shuttle mission STS-133 was scrubbed due to a gaseous hydrogen leak located in a vent line near the ground umbilical and ET connection. Subsequent visual inspections identified cracks in the sprayed-on foam insulation in the forward end of the ET intertank segment, adjacent to the liquid oxygen (LOX) tank, as shown in Figure 1. These cracks necessitated repair of the foam due to debris concerns that violated launch constraints. As part of the repair process, the affected foam was removed to reveal cracks in the underlying external hat stiffeners on the intertank, as shown in Figure 2. Ultimately, five stiffeners were discovered to be cracked adjacent to the LOX tank. As the managing center for the ET Project, NASA Marshall Space Flight Center (MSFC) coordinated failure investigation and repair activities among multiple organizations, which included the ET prime contractor (Lockheed Martin Space Systems Michoud Operations), the Space Shuttle Program Office at the NASA Johnson Space Center (JSC), the NASA Kennedy Space Center (KSC), and the NASA Engineering and Safety Center (NESC). STS-133 utilized the external tank designated as ET-137. Many aspects of the investigation have been reported previously in Refs. 1-7, which focus on the root cause of the failures, the flight readiness rationale and the local analyses of the stringer failures and repair. This paper summarizes the global analyses that were conducted on ET-137 as part of the NESC effort during the investigation, which was conducted primarily to determine if the repairs that were introduced to the stringers would alter the global response of the ET. In the process of the investigation, a new STAGS tabular input capability was developed to more easily introduce the aerodynamic pressure loads using a method that could easily be extended to incorporate finite element property data such as skin and stiffener thicknesses and beam cross-sectional properties.

Lovejoy, Andrew E.↗
NASA NTRS2003

Hierarchical Strategy for Rapid Analysis Environment

A new philosophy is developed wherein the hierarchical definition of data is made use of in creating a better environment to conduct analyses of practical problems. This system can be adapted to conduct virtually any type of analysis, since this philosophy is not bound to any specific kind of analysis. It provides a framework to manage different models and its results and more importantly, the interaction between the different models. Thus, it is ideal for many types of finite element analyses like globalAoca1 analysis and those that involve multiple scales and fields. The system developed during the course of this work is just a demonstrator of the basic concepts. A complete implementation of this strategy could potentially make a major impact on the way analyses are conducted. It could considerably reduce the time frame required to conduct the analysis of real-life problems by efficient management of the data involved and reducing the human effort involved. It also helps in better decision making because of more ways to interpret the results. The strategy has been currently implemented for structural analysis, but with more work it could be extended to other fields of science when the finite element method is used to solve the differential equations numerically. This report details the work that has been done during the course of this project and its achievements and results. The following section discusses the meaning of the word hierarchical and the different references to the term in the literature. It talks about the development of the finite element method, its different versions and how hierarchy has been used to improve the methodology. The next section describes the hierarchical philosophy in detail and explains the different concepts and terms associated with it. It goes on to describe the implementation and the features of the demonstrator. A couple of problems are analyzed using the demonstrator program to show the working of the system. The two problems considered are two dimensional plane stress analysis problems. The results are compared with those obtained using conventional analysis. The different challenges faced during the development of this system are discussed. Finally, we conclude with suggestions for future work to add more features and extend it to a wider range of problems.

Whitcomb, John↗
DOE OSTI2016

Exploration of efficient reduced-order modeling and a posteriori error estimation

Efficient algorithms are considered for the computation of a reduced-order model based on the proper orthogonal decomposition methodology for the solution of parameterized elliptic partial differential equations. The method relies on partitioning the parameter space into subdomains based on the properties of the solution space and then forming a reduced basis for each of the subdomains. This yields more efficient offline and online stages for the proper orthogonal decomposition method. We extend these ideas for inexpensive adjoint based a posteriori error estimation of both the expensive finite element method solutions and the reduced-order model solutions, for a single and multiple quantities of interest. Finally, various numerical results indicate the efficacy of the approach.

42 ENGINEERING↗
NASA NTRS1992

Coupled finite-difference/finite-element approach for wing-body aeroelasticity

Computational methods using finite-difference approaches for fluids and finite-element approaches for structures have individually advanced to solve almost full-aircraft configurations. However, coupled approaches to solve fluid/structural interaction problems are still in their early stages of development, particularly for complex geometries using complete equations such as the Euler/Navier-Stokes equations. Earlier work demonstrated the success of coupling finite-difference and finite-element methods for simple wing configurations using the Euler/Navier-Stokes equations. In this paper, the same approach is extended for general wing-body configurations. The structural properties are represented by beam-type finite elements. The flow is modeled using the Euler/Navier-Stokes equations. A general procedure to fully couple structural finite-element boundary conditions with fluid finite-difference boundary conditions is developed for wing-body configurations. Computations are made using moving grids that adapt to wing-body structural deformations. Results are illustrated for a typical wing-body configuration.

Guruswamy, Guru P.↗
DOE OSTI2021

Discrete Green’s functions and spectral graph theory for computationally efficient thermal modeling

Here, this work concerns solutions of the heat equation with the spectral graph method, for which the temperature is defined at discrete points in the domain and the spatial relationship among the points is described by a graph. The heat equation on the graph is solved using matrix techniques involving the eigenvectors and eigenvalues of the Laplacian matrix. The spectral graph approach precludes the computationally intensive meshing and numerous time-integration steps of the finite element method. In the present work, the spectral graph method is extended to include heat loss at the boundaries with a generalized boundary condition, and physics-based edge weights are introduced which simplify the calibration process. From this approach a discrete Green’s function is defined which allows for solutions under a variety of heating conditions including: space-varying initial conditions; time-and-space varying internal heating; and, time-and-space-varying heating at boundaries of type 1 (Dirichlet), type 2 (Neumann) and type 3 (Robin). Results are provided for benchmark heat transfer problems in one spatial dimension and in three spatial dimensions, and verification is provided by comparison with exact analytical solutions and finite difference solutions. The spectral graph method converges within 0.4% error of the analytical solution. The practical utility of the approach is demonstrated by thermal simulation of a multilayer additive manufacturing process. The spectral graph results are compared to experimentally-obtained temperature data for two metal parts, with error less than 5% of the experimental measurements, with computation time less than one minute on a desktop computer.

36 MATERIALS SCIENCE↗
DOE OSTI2019

High-order gradients with the shifted boundary method: An embedded enriched mixed formulation for elliptic PDEs

Here, we propose an extension of the embedded boundary method known as “shifted boundary method” to elliptic diffusion equations in mixed form (e.g., Darcy flow, heat diffusion problems with rough coefficients, etc.). Our aim is to obtain an improved formulation that, for linear finite elements, is at least second-order accurate for both flux and primary variable, when either Dirichlet or Neumann boundary conditions are applied. Following previous work of Nishikawa and Mazaheri in the context of residual distribution methods, we consider the mixed form of the diffusion equation (i.e., with Darcy-type operators), and introduce an enrichment of the primary variable. This enrichment is obtained exploiting the relation between the primary variable and the flux variable, which is explicitly available at nodes in the mixed formulation. The proposed enrichment mimics a formally quadratic pressure approximation, although only nodal unknowns are stored, similar to a linear finite element approximation. We consider both continuous and discontinuous finite element approximations and present two approaches: a non-symmetric enrichment, which, as in the original references, only improves the consistency of the overall method; and a symmetric enrichment, which enables a full error analysis in the classical finite element context. Combined with the shifted boundary method, these two approaches are extended to high-order embedded computations, and enable the approximation of both primary and flux (gradient) variables with second-order accuracy, independently on the type of boundary conditions applied. We also show that the primary variable is third-order accurate, when pure Dirichlet boundary conditions are embedded.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Mechanical and Microstructural Effects of Thermal Aging on Cast Duplex Stainless Steels by Experiment and Finite Element Method

Cast duplex stainless steel piping in light water nuclear reactors expe- rience thermal aging embrittlement during operational service. Interest in extending the operational life to 80 years requires an increased understanding of the microstructural evolution and corresponding changes in mechanical behavior. We analyze the evolution of the microstructure during thermal aging of cast CF-3 and CF-8 stainless steels using electron microscopy and atom probe tomography. The evolution of the mechanical properties is measured concurrently by mechanical methods such as tensile tests, Charpy V-notch tests, and instrumented nanoinden- tation. A microstructure-based finite element method model is developed and uti- lized in conjunction with the characterization results in order to correlate the local stress-strain effects in the microstructure with the bulk measurements. This work is supported by the DOE Nuclear Energy University Programs (NEUP), contract number DE-NE0000724.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Newmark local time stepping on high-performance computing architectures

In multi-scale complex media, finite element meshes often require areas of local refinement, creating small elements that can dramatically reduce the global time-step for wave-propagation problems due to the CFL condition. Local time stepping (LTS) algorithms allow an explicit time-stepping scheme to adapt the time-step to the element size, allowing near-optimal time-steps everywhere in the mesh. We develop an efficient multilevel LTS-Newmark scheme and implement it in a widely used continuous finite element seismic wave-propagation package. In particular, we extend the standard LTS formulation with adaptations to continuous finite element methods that can be implemented very efficiently with very strong element-size contrasts (more than 100x). Capable of running on large CPU and GPU clusters, we present both synthetic validation examples and large scale, realistic application examples to demonstrate the performance and applicability of the method and implementation on thousands of CPU cores and hundreds of GPUs.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Analytical modelling of laminated composites

This paper presents a unique approach to the analysis of thick laminated composites by presenting two simple finite element methods. The first uses the predictor corrector technique to extend the simple Mindlin-type element to achieve greater accuracy, and the second develops a new least squares element which can approximate a C(1) continuous element. The least squares element has the capability to incorporate a simplified higher-order basis into a piecewise continuous displacement field creating an accurate, yet computationally simple, element. These two methods have the potential to significantly upgrade analysis methods with little additional computational cost.

Kouri, Jeffrey V.↗
NASA NTRS1993

NASA Langley developments in response calculations needed for failure and life prediction

NASA Langley developments in response calculations needed for failure and life predictions are discussed. Topics covered include: structural failure analysis in concurrent engineering; accuracy of independent regional modeling demonstrated on classical example; functional interface method accurately joins incompatible finite element models; interface method for insertion of local detail modeling extended to curve pressurized fuselage window panel; interface concept for joining structural regions; motivation for coupled 2D-3D analysis; compression panel with discontinuous stiffener coupled 2D-3D model and axial surface strains at the middle of the hat stiffener; use of adaptive refinement with multiple methods; adaptive mesh refinement; and studies on quantity effect of bow-type initial imperfections on reliability of stiffened panels.

Housner, Jerrold M.↗
DOE OSTI2020

NRC Multiphysics Analysis Capability Deployment FY2020: Part 3

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Sup- port for Advanced Non-Light Water Reactors.” The tasks completed for this report are: Task2c: Explicit modeling of pebble transient temperature response. In this simulation, the 400 MWth Pebble-Bed Modular Reactor (PBMR) design, PBMR- 400, experiences a 20-second power ramp from 100% to 150% power. This is followed by a similar reduction in the power back to 100%. Several multiscale pebble coupling approaches are tested with one pebble per mesh element in the active core region. The results show good conservation behavior and the stability of the coupling.; Extended scope part 1: An assessment of the computational efficiency of the Discontinuous Finite Element Method (DFEM) heat transfer solver shows good scalability. The DFEM solver is a factor of 4 more expensive in solution time than the Finite Element Method (FEM) solver for heat transfer problems due to the increased number of degrees of freedom. Nonetheless, the DFEM approach provides the user with the flexibility to model gap heat transfer problems.; Extended scope part 2: The GapHeatTransferInterfaceMaterial was improved to give the user increased flexibility with the modeling of heat transfer through gaps with the DFEM solver. A number of gap parameters can now be coupled both through functions and variables.; Extended scope part 3: Demonstration of how the gap width between hexagonal fuel cells can be calculated during a heat-up transient and used in the GapHeatTransferInterface model. A full-domain DFEM model with gap expansion is coupled to a SubApp that models the thermal expansion of the base plate. The results show the expected physical behavior, although have not been fully bench-marked at this point in time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1997

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗
DOE OSTI2018

High-Order Multi-Material ALE Hydrodynamics

Here, we present a new approach for multi-material arbitrary Lagrangian--Eulerian (ALE) hydrodynamics simulations based on high-order finite elements posed on high-order curvilinear meshes. The method builds on and extends our previous work in the Lagrangian and remap phases of ALE, and depends critically on a functional perspective that enables subzonal physics and material modeling. Curvilinear mesh relaxation is based on node movement, which is determined through the solution of an elliptic equation. The remap phase is posed in terms of advecting state variables between two meshes over a fictitious time interval. The resulting advection equation is solved by a discontinuous Galerkin (DG) formulation, combined with a customized Flux Corrected Transport (FCT) type algorithm. Because conservative fields are remapped, additional synchronization steps are introduced to preserve bounds with respect to primal fields. These steps include modification of the low-order FCT solutions, definition of conservative FCT fluxes based on primal field bounds, and monotone transitions between primal and conservative fields.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Unit cell geometry of multiaxial preforms for structural composites

The objective of this study is to investigate the yarn geometry of multiaxial preforms. The importance of multiaxial preforms for structural composites is well recognized by the industry but, to exploit their full potential, engineering design rules must be established. This study is a step in that direction. In this work the preform geometry for knitted and braided preforms was studied by making a range of well designed samples and studying them by photo microscopy. The structural geometry of the preforms is related to the processing parameters. Based on solid modeling and B-spline methodology a software package is developed. This computer code enables real time structural representations of complex fiber architecture based on the rule of preform manufacturing. The code has the capability of zooming and section plotting. These capabilities provide a powerful means to study the effect of processing variables on the preform geometry. the code also can be extended to an auto mesh generator for downstream structural analysis using finite element method. This report is organized into six sections. In the first section the scope and background of this work is elaborated. In section two the unit cell geometries of braided and multi-axial warp knitted preforms is discussed. The theoretical frame work of yarn path modeling and solid modeling is presented in section three. The thin section microscopy carried out to observe the structural geometry of the preforms is the subject in section four. The structural geometry is related to the processing parameters in section five. Section six documents the implementation of the modeling techniques into the computer code MP-CAD. A user manual for the software is also presented here. The source codes and published papers are listed in the Appendices.

Ko, Frank↗
NASA NTRS2001

A 3-D Magnetic Analysis of a Stirling Convertor Linear Alternator Under Load

The NASA Glenn Research Center (GRC), the Department of Energy (DOE), and the Stirling Technology Company (STC) are developing Stirling convertors for Stirling Radioisotope Power Systems (SRPS) to provide electrical power for future NASA deep space missions. STC is developing the 55-We Technology Demonstration Convertor (TDC) under contract to DOE. Of critical importance to the successful development of the Stirling convertor for space power applications is the development of a lightweight and highly efficient linear alternator. This paper presents a 3-dimensional finite element method (FEM) approach for evaluating Stirling convertor linear alternators. The model extends a magnetostatic analysis previously reported at the 35th Intersociety Energy Conversion Engineering Conference (IECEC) to include the effects of the load current. STC's 55-We linear alternator design was selected to validate the model. Spatial plots of magnetic field strength (H) are presented in the region of the exciting permanent magnets. The margin for permanent magnet demagnetization is calculated at the expected magnet operating temperature for the near earth environment and for various average magnet temperatures. These thermal conditions were selected to represent a worst-case condition for the planned deep space missions. This paper presents plots that identify regions of high H where the potential to alter the magnetic moment of the magnets exists.

Geng, Steven M.↗
NASA NTRS1990

Application of the probabilistic approximate analysis method to a turbopump blade analysis

An eigenvalue analysis of a typical space propulsion system turbopump blade is presented using an approximate probabilistic analysis methodology. The methodology was developed originally to investigate the feasibility of computing probabilistic structural response using closed-form approximate models. This paper extends the methodology to structures for which simple closed-form solutions do not exist. The finite element method will be used for this demonstration, but the concepts apply to any numerical method. The results agree with detailed analysis results and indicate the usefulness of using a probabilistic approximate analysis in determining efficient solution strategies.

Thacker, B. H.↗