Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element problems”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. The transient finite-element analysis is done using TACS, and the cell-discharge battery model uses OpenMDAO and dymos. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization↗
NASA NTRS1996

Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation

We study parallel two-level overlapping Schwarz algorithms for solving nonlinear finite element problems, in particular, for the full potential equation of aerodynamics discretized in two dimensions with bilinear elements. The overall algorithm, Newton-Krylov-Schwarz (NKS), employs an inexact finite-difference Newton method and a Krylov space iterative method, with a two-level overlapping Schwarz method as a preconditioner. We demonstrate that NKS, combined with a density upwinding continuation strategy for problems with weak shocks, is robust and, economical for this class of mixed elliptic-hyperbolic nonlinear partial differential equations, with proper specification of several parameters. We study upwinding parameters, inner convergence tolerance, coarse grid density, subdomain overlap, and the level of fill-in in the incomplete factorization, and report their effect on numerical convergence rate, overall execution time, and parallel efficiency on a distributed-memory parallel computer.

Cai, Xiao-Chuan↗
NASA NTRS1983

Finite-element methods for steady solidification problems

Four Galerkin finite-element methods are tested for solving the free-boundary problem that describes steady solidification. The formulations differ in the solution method used to account for the unknown shape of the melt/solid interface, in the interphase condition (either balance of heat flux or equilibrium of temperature) distinguished for locating the interface, and in the technique used for solving the systems of algebraic equations that result from the finite-element approximations. Methods that use the melting point isotherm to locate the melt/solid interface are found more accurate and efficient than formulations based on the interfacial energy balance. Solution by a Galerkin-Newton algorithm of the free-boundary problem transformed to a fixed domain is most efficient when the field problem in each phase is made nonlinear by including radiation from the melt and solid to the surroundings.

Ettouney, H. M.↗
NASA NTRS1993

Coupling finite and boundary element methods for 2-D elasticity problems

A finite element-boundary element (FE-BE) coupling method for two-dimensional elasticity problems is developed based on a weighted residual variational method in which a portion of the domain of interest is modeled by FEs and the remainder of the region by BEs. The performance of the FE-BE coupling method is demonstrated via applications to a simple 'patch test' problem and three-crack problems. The method passed the patch tests for various modeling configurations and yielded accurate strain energy release rates for the crack problems studied.

Krishnamurthy, T.↗
NASA NTRS1993

Solution of free-boundary problems using finite-element/Newton methods and locally refined grids - Application to analysis of solidification microstructure

A new method is presented for the solution of free-boundary problems using Lagrangian finite element approximations defined on locally refined grids. The formulation allows for direct transition from coarse to fine grids without introducing non-conforming basis functions. The calculation of elemental stiffness matrices and residual vectors are unaffected by changes in the refinement level, which are accounted for in the loading of elemental data to the global stiffness matrix and residual vector. This technique for local mesh refinement is combined with recently developed mapping methods and Newton's method to form an efficient algorithm for the solution of free-boundary problems, as demonstrated here by sample calculations of cellular interfacial microstructure during directional solidification of a binary alloy.

Tsiveriotis, K.↗
NASA NTRS1982

The SPAR thermal analyzer: Present and future

The SPAR thermal analyzer, a system of finite-element processors for performing steady-state and transient thermal analyses, is described. The processors communicate with each other through the SPAR random access data base. As each processor is executed, all pertinent source data is extracted from the data base and results are stored in the data base. Steady state temperature distributions are determined by a direct solution method for linear problems and a modified Newton-Raphson method for nonlinear problems. An explicit and several implicit methods are available for the solution of transient heat transfer problems. Finite element plotting capability is available for model checkout and verification.

Marlowe, M. B.↗
NASA NTRS1983

Potential of minicomputer/array-processor system for nonlinear finite-element analysis

The potential of using a minicomputer/array-processor system for the efficient solution of large-scale, nonlinear, finite-element problems is studied. A Prime 750 is used as the host computer, and a software simulator residing on the Prime is employed to assess the performance of the Floating Point Systems AP-120B array processor. Major hardware characteristics of the system such as virtual memory and parallel and pipeline processing are reviewed, and the interplay between various hardware components is examined. Effective use of the minicomputer/array-processor system for nonlinear analysis requires the following: (1) proper selection of the computational procedure and the capability to vectorize the numerical algorithms; (2) reduction of input-output operations; and (3) overlapping host and array-processor operations. A detailed discussion is given of techniques to accomplish each of these tasks. Two benchmark problems with 1715 and 3230 degrees of freedom, respectively, are selected to measure the anticipated gain in speed obtained by using the proposed algorithms on the array processor.

Strohkorb, G. A.↗
NASA NTRS1984

Finite element computation on nearest neighbor connected machines

Research aimed at faster, more cost effective parallel machines and algorithms for improving designer productivity with finite element computations is discussed. A set of 8 boards, containing 4 nearest neighbor connected arrays of commercially available floating point chips and substantial memory, are inserted into a commercially available machine. One-tenth Mflop (64 bit operation) processors provide an 89% efficiency when solving the equations arising in a finite element problem for a single variable regular grid of size 40 by 40 by 40. This is approximately 15 to 20 times faster than a much more expensive machine such as a VAX 11/780 used in double precision. The efficiency falls off as faster or more processors are envisaged because communication times become dominant. A novel successive overrelaxation algorithm which uses cyclic reduction in order to permit data transfer and computation to overlap in time is proposed.

Mcaulay, A. D.↗
NASA NTRS1984

Finite elements and the method of conjugate gradients on a concurrent processor

An algorithm for the iterative solution of finite element problems on a concurrent processor is presented. The method of conjugate gradients is used to solve the system of matrix equations, which is distributed among the processors of a MIMD computer according to an element-based spatial decomposition. This algorithm is implemented in a two-dimensional elastostatics program on the Caltech Hypercube concurrent processor. The results of tests on up to 32 processors show nearly linear concurrent speedup, with efficiencies over 90% for sufficiently large problems.

Lyzenga, G. A.↗
NASA NTRS1986

Error-source effects in a high-accuracy optical finite-element processor

High-accuracy optical linear algebra processors are addressed with attention to three new aspects. These include: their application to the solution of finite-element problems; the first error-source models for component errors in such processors; and the first analysis of error sources in such processors.

Taylor, B. K.↗
NASA NTRS1985

Finite elements and the method of conjugate gradients on a concurrent processor

An algorithm for the iterative solution of finite element problems on a concurrent processor is presented. The method of conjugate gradients is used to solve the system of matrix equations, which is distributed among the processors of a MIMD computer according to an element-based spatial decomposition. This algorithm is implemented in a two-dimensional elastostatics program on the Caltech Hypercube concurrent processor. The results of tests on up to 32 processors show nearly linear concurrent speedup, with efficiencies over 90 percent for sufficiently large problems.

Lyzenga, G. A.↗
NASA NTRS1988

Direct finite element solution on an optical laboratory matrix-vector processor

The first optical laboratory system results employing a direct LU decomposition solution of a system of linear algebraic equations are presented for a finite element problem solution. This also represents the first laboratory demonstration of the use of sign-magnitude negative number representation as well as new bit partitioning techniques to increase the accuracy of an optical encoded processor beyond the number of bit channels available.

Casasent, David↗
NASA NTRS1993

Transient Finite Element Computations on a Variable Transputer System

A parallel program to analyze transient finite element problems was written and implemented on a system of transputer processors. The program uses the explicit time integration algorithm which eliminates the need for equation solving, making it more suitable for parallel computations. An interprocessor communication scheme was developed for arbitrary two dimensional grid processor configurations. Several 3-D problems were analyzed on a system with a small number of processors.

Smolinski, Patrick J.↗
NASA NTRS2018

Use of Simple Continuum Solutions in Finite Element Alternating Method for Fracture Problems

The performance of the finite element alternating (FEAM) method for two-dimensional crack problems is studied with respect to a polynomial pressure distribution fitted to the crack face stresses. The FEAM alternates between the analytical solution of crack in an infinite plate subjected to arbitrary polynomial distribution and a finite element solution of an uncracked body to satisfy the required boundary conditions in the crack problem. In this paper, the FEAM is applied to embedded crack and edge crack problems. For embedded crack problems, all of the constant, linear, and quadratic ( N=0,1, or 2, respectively) pressure distributions yield very accurate results with this algorithm with 4 to 5 iterations. The edge crack problems, on the other hand, require much higher order polynomials distributions (N=5 to 6) to yield accurate solutions. For slant edge crack problems, the mode-I stress-intensity factors have better accuracy than the mode-II stress-intensity factors for the same convergence tolerance.

Krishnamurthy, T.↗
NASA NTRS1992

Coupling finite and boundary element methods for two-dimensional potential problems

A finite element-boundary element (FE-BE) coupling method based on a weighted residual variational method is presented for potential problems, governed by either the Laplace or the Poisson equations. In this method, a portion of the domain of interest is modeled by finite elements (FE) and the remainder of the region by boundary elements (BE). Because the BE fundamental solutions are valid for infinite domains, a procedure that limits the effects of the BE fundamental solution to a small region adjacent to the FE region, called the transition region (TR), is developed. This procedure involves a judicious choice of functions called the transition (T) functions that have unit values on the BE-TR interface and zero values on the FE-TR interface. The present FE-BE coupling algorithm is shown to be independent of the extent of the transition region and the choice of the transition functions. Therefore, transition regions that extend to only one layer of elements between FE and BE regions and the use of simple linear transition functions work well.

Krishnamurthy, T.↗
NASA NTRS1980

The finite element duct eigenvalue problem - An improved formulation with Hermitian elements and no-flow condensation

Hermitian elements are used in a finite element solution for the eigenvalue problem in lined ducts with flow. These elements give significantly greater accuracy for reduced dimensionality when compared with Lagrangian elements. Spurious mode generation associated with the Lagrangian formulation is eliminated. A dramatic improvement in the ratio of the number of reliable eigenvalues to the total number of computed eigenvalues is effected by the use of a condensation scheme based on the no-flow eigenvectors. Results are presented for two dimensional and axisymmetric ducts. In the axisymmetric case good resolution is obtained even for high order, high frequency modes by the use of continuously graded meshes.

Astley, R. J.↗