Search NASA⌕ Search

SEARCH · Search NASA

Results for “Integral Equations 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

The solution of linear systems of equations with a structural analysis code on the NAS CRAY-2

Two methods for solving linear systems of equations on the NAS Cray-2 are described. One is a direct method; the other is an iterative method. Both methods exploit the architecture of the Cray-2, particularly the vectorization, and are aimed at structural analysis applications. To demonstrate and evaluate the methods, they were installed in a finite element structural analysis code denoted the Computational Structural Mechanics (CSM) Testbed. A description of the techniques used to integrate the two solvers into the Testbed is given. Storage schemes, memory requirements, operation counts, and reformatting procedures are discussed. Finally, results from the new methods are compared with results from the initial Testbed sparse Choleski equation solver for three structural analysis problems. The new direct solvers described achieve the highest computational rates of the methods compared. The new iterative methods are not able to achieve as high computation rates as the vectorized direct solvers but are best for well conditioned problems which require fewer iterations to converge to the solution.

Poole, Eugene L.↗
NASA NTRS1989

Semidiscrete Galerkin modelling of compressible viscous flow past a circular cone at incidence

A numerical study of the laminar and compressible boundary layer, about a circular cone in a supersonic free stream, is presented. It is thought that if accurate and efficient numerical schemes can be produced to solve the boundary layer equations, they can be joined to numerical codes that solve the inviscid outer flow. The combination of these numerical codes is competitive with the accurate, but computationally expensive, Navier-Stokes schemes. The primary goal is to develop a finite element method for the calculation of 3-D compressible laminar boundary layer about a yawed cone. The proposed method can, in principle, be extended to apply to the 3-D boundary layer of pointed bodies of arbitrary cross section. The 3-D boundary layer equations governing supersonic free stream flow about a cone are examined. The 3-D partial differential equations are reduced to 2-D integral equations by applying the Howarth, Mangler, Crocco transformations, a linear relation between viscosity, and a Blasius-type of similarity variable. This is equivalent to a Dorodnitsyn-type formulation. The reduced equations are independent of density and curvature effects, and resemble the weak form of the 2-D incompressible boundary layer equations in Cartesian coordinates. In addition the coordinate normal to the wall has been stretched, which reduces the gradients across the layer and provides high resolution near the surface. Utilizing the parabolic nature of the boundary layer equations, a finite element method is applied to the Dorodnitsyn formulation. The formulation is presented in a Petrov-Galerkin finite element form and discretized across the layer using linear interpolation functions. The finite element discretization yields a system of ordinary differential equations in the circumferential direction. The circumferential derivatives are solved by an implicit and noniterative finite difference marching scheme. Solutions are presented for a 15 deg half angle cone at angles of attack of 5 and 10 deg. The numerical solutions assume a laminar boundary layer with free stream Mach number of 7. Results include circumferential distribution of skin friction and surface heat transfer, and cross flow velocity distributions across the layer.

Meade, Andrew James, Jr.↗
DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗
NASA NTRS1990

Integrated fluid-thermal-structural analysis using adaptive unstructured meshes

An integrated fluid-thermal-structural analysis approach using adaptive unstructured meshes is presented. The approach combines the finite-element method and an adaptive remeshing technique to solve the Navier-tokes equations for high-speed compressible flow, the energy equation for the thermal response of the structure, and the quasi-static equilibrium equations for the structural response. The analysis solution procedure and the adaptive unstructured remeshing technique are described. The effectiveness of the approach is evaluated with three application studies. The adaptive unstructured remeshing procedure and finite-element solution algorithms combine to yield increased accuracy and effiency over standard structured meshes.

Dechaumphai, Pramote↗
NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 1: Theoretical developments and applications

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present first part of the report, we focus on the theoretical developments, and discussions of the results of numerical-performance studies using the integration schemes for GVIPS and NAV models.

Saleeb, Atef F.↗
NASA NTRS1981

Transonic flow computations in cascades using finite element method

Analysis of transonic flow through a cascade of airfoils is investigated using the finite element method. Development of a computational grid suitable for complex flow structures and different types of boundary conditions is presented. An efficient pseudo-time integration scheme is developed for the solution of equations. Modeling of the shock and the convergence characteristics of the developed scheme are discussed. Numerical results include a 45 deg staggered cascade of NACA 0012 airfoils with inlet flow Mach number of 0.8 and angles of attack 1, 0, and 1 deg.

Akay, H. U.↗
NASA NTRS1982

Algorithmic aspects of transient heat transfer problems in structures

It is noted that the application of finite element or finite difference techniques to the solution of transient heat transfer problems in structures often results in a stiff system of ordinary differential equations. Such systems are usually handled most efficiently by implicit integration techniques which require the solution of large and sparse systems of algebraic equations. The assembly and solution of these systems using the incomplete Cholesky conjugate gradient algorithm is examined. Several examples are used to demonstrate the advantage of the algorithm over other techniques.

Haftka, R. T.↗
NASA NTRS1985

Application of closed-form solutions to a mesh point field in silicon solar cells

A computer simulation method is discussed that provides for equivalent simulation accuracy, but that exhibits significantly lower CPU running time per bias point compared to other techniques. This new method is applied to a mesh point field as is customary in numerical integration (NI) techniques. The assumption of a linear approximation for the dependent variable, which is typically used in the finite difference and finite element NI methods, is not required. Instead, the set of device transport equations is applied to, and the closed-form solutions obtained for, each mesh point. The mesh point field is generated so that the coefficients in the set of transport equations exhibit small changes between adjacent mesh points. Application of this method to high-efficiency silicon solar cells is described; and the method by which Auger recombination, ambipolar considerations, built-in and induced electric fields, bandgap narrowing, carrier confinement, and carrier diffusivities are treated. Bandgap narrowing has been investigated using Fermi-Dirac statistics, and these results show that bandgap narrowing is more pronounced and that it is temperature-dependent in contrast to the results based on Boltzmann statistics.

Lamorte, M. F.↗
NASA NTRS1989

Fluid-thermal-structural interaction of aerodynamically heated leading edges

A two-dimensional finite element approach is presented for the integrated fluid-thermal-structural analysis of aerodynamically heated leading edges. The approach is combined with an adaptive unstructured remeshing technique to solve the Navier-Stokes equations for high speed compressible flow, the energy equation for the structure thermal response, and the quasi-static equilibrium equations for the structural response. Coupling and interaction between the three disciplines are demonstrated using two applications for high speed flow over a cylinder and a simulated engine leading edge verification test.

Dechaumphai, Pramote↗
NASA NTRS1991

Calculating Transient Vibrations Of Coupled Substructures

Finite-element, numerical-integration method for estimating transient vibrational response of structure composed of coupled substructures entails less computation. Responses of substructures to external forces and to forces of interaction between substructures computed. Equations of motion of each substructure solved independently. Applicable to any number and configuration of linearly responding substructures and to both determinate and indeterminate boundary conditions at interfaces between substructures. Also applicable to changing interface boundary conditions.

Admire, J. R.↗
NASA NTRS1991

Time-accurate incompressible flow computations with quadrilateral velocity-pressure elements

Quadrilateral velocity-pressure elements with constant and linear pressure interpolations are examined in the context of time-accurate finite element computation of unsteady incompressible flows. These elements involve streamline-upwind/Petrov-Galerkin stabilization and are implemented in conjunction with the one-step and multi-step temporal integration of the Navier-Stokes equations. The two test cases chosen for the performance evaluation of the formulations are the standing vortex problem and flow past a circular cylinder at Reynolds number 100.

Tezduyar, T. E.↗
NASA NTRS1999

Passively Damped Laminated Piezoelectric Shell Structures with Integrated Electric Networks

Multi-field mechanics are presented for curvilinear piezoelectric laminates interfaced with distributed passive electric components. The equations of motion for laminated piezoelectric shell structures with embedded passive electric networks are directly formulated and solved using a finite element methodology. The modal damping and frequencies of the piezoelectric shell are calculated from the poles of the system. Experimental and numerical results are presented for the modal damping and frequency of composite beams with a resistively shunted piezoceramic patch. The modal damping and frequency of plates, cylindrical shells and cylindrical composite blades with piezoelectric-resistor layers are predicted. Both analytical and experimental studies illustrate a unique dependence of modal damping and frequencies on the shunting resistance and show the effect of structural shape and curvature on piezoelectric damping.

Saravanos, Dimitris A.↗
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↗
NASA NTRS1989

A finite element approach for large motion dynamic analysis of multibody structures in space

A three-dimensional finite element formulation for modeling the transient dynamics of constrained multibody space sructures with truss-like configurations is presented. Convected coordinate systems are used to define rigid-body motion of individual elements in the system. These systems are located at one end of each element and are oriented such that one axis passes through the other end of the element. Deformation of each element, relative to its convected coordinate system, is defined by cubic flexural shape functions as used in finite element methods of structural analysis. The formulation is oriented toward joint dominated structures and places the generalized coordinates at the joint. A transformation matrix is derived to integrate joint degree-of-freedom into the equations of motion of the element. Based on the derivation, a general-purpose code LATDYN (Large Angle Transient DYNamics) was developed. Two examples are presented to illustrate the application of the code. For the spin-up of a flexible beam, results are compared with existing solutions available in the literature. For the deployment of one bay of a deployable space truss (the Minimast), results are verified by the geometric knowledge of the system and converged solution of a successively refined model.

Chang, Che-Wei↗
NASA NTRS1989

Dynamic analysis of flexible multibody structures

A system composed of several interconnected elastic components that may experience large angular motion relative to each other during operation is referred to as a flexible multibody structure. Several formulations were proposed for the determination of the dynamic response of controlled flexible multibody structures. In general, these formulations consist of superposing elastic deformations of the component body onto the large rigid body motion of the component. It was shown that this particular methodology for combining linear structural deformations with nonlinear kinematics can lead to erroneous response predictions when either the beam member is very flexible or the rotational speed is high. In addition, previous formulations introduce constraint equations to define the interrelations among system components. This approach increases the number of equations that must be solved, and may result in contraint violation when numerical error accumulates during the integration process. In order to overcome the difficulties, a new approach was suggested. The approach is essentially a finite element formulation which takes advantage of the fact that many multibody structures are joint dominated. The Large Angle Transient Dynamic Analysis (LATDYN) program for clarity of documentation, ease of use, user friendliness, modeling generality, and accuracy of results was evaluated. This required gaining a working familiarity with the code and performing several case studies.

Hernried, Alan G.↗
NASA NTRS1990

Near-wall turbulence model and its application to fully developed turbulent channel and pipe flows

A near-wall turbulence model and its incorporation into a multiple-timescale turbulence model are presented. The near-wall turbulence model is obtained from a k-equation turbulence model and a near-wall analysis. In the method, the equations for the conservation of mass, momentum, and turbulent kinetic energy are integrated up to the wall, and the energy transfer and the dissipation rates inside the near-wall layer are obtained from algebraic equations. Fully developed turbulent channel and pipe flows are solved using a finite element method. The computational results compare favorably with experimental data. It is also shown that the turbulence model can resolve the overshoot phenomena of the turbulent kinetic energy and the dissipation rate in the region very close to the wall.

Kim, S.-W.↗
NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 2: Algorithmic developments and implementation

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Li, Wei↗
DOE OSTI2020

Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids. LDRD Final Report

This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.

97 MATHEMATICS AND COMPUTING↗