Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

A finite element study of incompressible flows past oscillating cylinders and aerofoils

We present our numerical results for certain unsteady flows past oscillating cylinders and aerofoils. The computations are based on the stabilized space-time finite element formulation. The implicit equation systems resulting from the space-time finite element discretizations are solved using iterative solution techniques. One of the problems studied is flow past a cylinder which is forced to oscillate in the horizontal direction. In this case we observe a change from an unsymmetric mode of vortex shedding to a symmetric one. An extensive study was carried out for the case in which a cylinder is mounted on lightly damped springs and allowed to oscillate in the vertical direction. In this case the motion of the cylinder needs to be determined as part of the solution, and under certain conditions this motion changes the vortex-shedding pattern of the flow field significantly. This nonlinear fluid-structure interaction exhibits certain interesting behavior such as 'lock-in' and 'hysteresis', which are in good agreement with the laboratory experiments carried out by other researchers in the past. Preliminary results for flow past a pitching aerofoil are also presented.

Mittal, S.↗
NASA NTRS1990

Large scale nonlinear numerical optimal control for finite element models of flexible structures

This paper discusses the development of large scale numerical optimal control algorithms for nonlinear systems and their application to finite element models of structures. This work is based on our expansion of the optimal control algorithm (DDP) in the following steps: improvement of convergence for initial policies in non-convex regions, development of a numerically accurate penalty function method approach for constrained DDP problems, and parallel processing on supercomputers. The expanded constrained DDP algorithm was applied to the control of a four-bay, two dimensional truss with 12 soft members, which generates geometric nonlinearities. Using an explicit finite element model to describe the structural system requires 32 state variables and 10,000 time steps. Our numerical results indicate that for constrained or unconstrained structural problems with nonlinear dynamics, the results obtained by our expanded constrained DDP are significantly better than those obtained using linear-quadratic feedback control.

Shoemaker, Christine A.↗
NASA NTRS2005

Finite Element Analysis of Adaptive-Stiffening and Shape-Control SMA Hybrid Composites

Shape memory alloy hybrid composites with adaptive-stiffening or morphing functions are simulated using finite element analysis. The composite structure is a laminated fiber-polymer composite beam with embedded SMA ribbons at various positions with respect to the neutral axis of the beam. Adaptive stiffening or morphing is activated via selective resistance heating of the SMA ribbons or uniform thermal loads on the beam. The thermomechanical behavior of these composites was simulated in ABAQUS using user-defined SMA elements. The examples demonstrate the usefulness of the methods for the design and simulation of SMA hybrid composites. Keywords: shape memory alloys, Nitinol, ABAQUS, finite element analysis, post-buckling control, shape control, deflection control, adaptive stiffening, morphing, constitutive modeling, user element

Gao, Xiujie↗
DOE OSTI2018

In-situ neutron diffraction and crystal plasticity finite element modeling to study the kinematic stability of retained austenite in bearing steels

This work integrates in-situ neutron diffraction and crystal plasticity finite element modeling to study the kinematic stability of retained austenite in high carbon bearing steels. The presence of a kinematically metastable retained austenite in bearing steels can significantly affect the macro-mechanical and micro-mechanical material response. Mechanical characterization of metastable austenite is a critical component in accurately capturing the micro-mechanical behavior under typical application loads. Traditional mechanical characterization techniques are unable to discretely quantify the micro-mechanical response of the austenite, and as a result, the computational predictions rely heavily on trial and error or qualitative descriptions of the austenite phase. In order to overcome this, in the present work, we use in-situ neutron diffraction of a uniaxial tension test of an A485 Grade 1 bearing steel specimen. The mechanical response determined from the neutron diffraction analysis was incorporated into a hybrid crystal plasticity finite element model that accounts for the martensite's crystal plasticity and the stress-assisted transformation from austenite to martensite in bearing steels. Here, the modeling response was used to estimate the single crystal elastic constants of the austenite and martensite phases. Finally, the results show that using in-situ neutron diffraction, coupled with a crystal plasticity model, can successfully predict both the micro-mechanical and macro-mechanical responses of bearing steels while accounting for the martensitic transformation of the retained austenite.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2010

Finite element modeling of the deformation of magnetoelastic film

Recently a new class of biocompatible elastic polymers loaded with small ferrous particles, a magnetoelastic polymer, has been developed. This engineered material is formed into a thin film using spin casting. An applied magnetic field will deform the film. The magnetic deformation of this film has many possible applications, particularly in microfluidic pumps and pressure regulators. In this paper a finite element method suitable for the transient simulation of arbitrarily shaped three-dimensional magnetoelastic polymers subjected to time-varying magnetic fields is developed. The approach is similar to that employed in finite elment magnetohydrodynamic simulations, the key difference is a more complex hyperelastic material model. In order to confirm the validity of the approach, finite element solutions for an axially symmetric thin film are compared to an analytical solution based on the membrane (infinitely thin) approximation. For this particular problem the two approaches give qualitatively similar results and converge as the film thickness approaches zero.

36 MATERIALS SCIENCE↗
NASA NTRS1977

PLANS; a finite element program for nonlinear analysis of structures. Volume 2: User's manual

The PLANS system, rather than being one comprehensive computer program, is a collection of finite element programs used for the nonlinear analysis of structures. This collection of programs evolved and is based on the organizational philosophy in which classes of analyses are treated individually based on the physical problem class to be analyzed. Each of the independent finite element computer programs of PLANS, with an associated element library, can be individually loaded and used to solve the problem class of interest. A number of programs have been developed for material nonlinear behavior alone and for combined geometric and material nonlinear behavior. The usage, capabilities, and element libraries of the current programs include: (1) plastic analysis of built-up structures where bending and membrane effects are significant, (2) three dimensional elastic-plastic analysis, (3) plastic analysis of bodies of revolution, and (4) material and geometric nonlinear analysis of built-up structures.

Pifko, A.↗
NASA NTRS1993

Elastic-plastic mixed-iterative finite element analysis: Implementation and performance assessment

An elastic-plastic algorithm based on Von Mises and associative flow criteria is implemented in MHOST-a mixed iterative finite element analysis computer program developed by NASA Lewis Research Center. The performance of the resulting elastic-plastic mixed-iterative analysis is examined through a set of convergence studies. Membrane and bending behaviors of 4-node quadrilateral shell finite elements are tested for elastic-plastic performance. Generally, the membrane results are excellent, indicating the implementation of elastic-plastic mixed-iterative analysis is appropriate.

Sutjahjo, Edhi↗
DOE OSTI2024

A Low-Rank QTT-based Finite Element Method for Elasticity Problems

We present an efficient and robust numerical algorithm for solving the linear elasticity problem that combines the Quantized Tensor Train format and a domain partitioning strategy. This approach makes it possible to solve the linear elasticity problem on a computational domain that is more general than a square. By integrating Z-ordering and subdomain concatenation, our method substantially decreases memory usage and achieves a notable reduction in rank compared to established Finite Element implementations like the FEniCS platform. This efficiency is maintained while still guaranteeing exponential convergence with respect to the number of degrees of freedom. This performance gain, however, requires a fundamental rethinking of how core finite element operations are implemented. This includes changes to mesh discretization, node and degree of freedom ordering, stiffness matrix and internal nodal force assembly, and the execution of algebraic matrix-vector operations. In this work, we discuss all these aspects in detail and assess the method’s performance in the numerical approximation of three representative test cases.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Generating Finite-Element Models Of Composite Materials

Program starts at micromechanical level, from simple inputs supplied by user. COMGEN, COmposite Model GENerator, is interactive FORTRAN program used to create wide variety of finite-element models of continuous-fiber composite materials at micromechanical level. Quickly generates batch or "session files" to be submitted to finite-element preprocessor and postprocessor program, PATRAN. COMGEN requires PATRAN to complete model.

Melis, M. E.↗
NASA NTRSDate not supplied

Finite Element Modeling for Compression Strength Prediction of Honeycomb Cores with Geometric Imperfections Measured using X-ray CT Imaging

Aluminum honeycomb cores have been used extensively in composite sandwich panels due to having high bending rigidity while maintaining low density. During manufacturing of the aluminum honeycomb, the thin metal walls are susceptible to imperfections that deviate from an ideal honeycomb shape. Computational tools to quantify imperfections and investigate their effects on quasi-static compression and impact response can inform design criteria to construct safer and lighter launch vehicle structures. This paper presents the results of a study on finite element modeling of metallic honeycomb cores (HCC) with geometric imperfections measured using X-ray Computer Tomography (CT), to predict the compression response. Finite element (FE) models are constructed with measured imperfections with appropriate boundary conditions. It is shown that the average behavior of larger 10x11 cell models can be predicted by sampling single cell models from the larger domain. Amplification of imperfections for cells at the center of sampled specimens are mode switched from the larger out-of-plane imperfection to the mode shape of the adjacent cells resulting in strengthening of the center specimens.

aluminum honeycomb↗
DOE OSTI2013

Three dimensional finite element methods: Their role in the design of DC accelerator systems

High Voltage Engineering has designed, built and tested a 2 MV dual irradiation system that will be applied for radiation damage studies and ion beam material modification. The system consists of two independent accelerators which support simultaneous proton and electron irradiation (energy range 100 keV - 2 MeV) of target sizes of up to 300 Multiplication-Sign 300 mm{sup 2}. Three dimensional finite element methods were used in the design of various parts of the system. The electrostatic solver was used to quantify essential parameters of the solid-state power supply generating the DC high voltage. The magnetostatic solver and ray tracing were used to optimize the electron/ion beam transport. Close agreement between design and measurements of the accelerator characteristics as well as beam performance indicate the usefulness of three dimensional finite element methods during accelerator system design.

36 MATERIALS SCIENCE↗
DOE OSTI2018

The stabilized mixed finite element scheme of elasticity problem

In this paper, we consider a new mixed finite element scheme of the elasticity problem in two and three dimensions, and the new scheme can impose Neumann boundary condition directly. A new stabilized method is proposed for this scheme, in which the equal order linear element pair is employed to approximate the stress and displacement, and an abstract operator is used to characterize the lack of the inf-sup condition. The new stabilized method is locking free. Numerical results show the excellent stability and accuracy of the new method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A discrete-continuum coupled finite element modelling approach for fibre reinforced concrete

Highlights: • An easy-to-implement but effective meso-scale modelling approach is developed for FRC. • It is capable of simulating and interpreting all main failure mechanisms in FRC. • It is critically validated by experiments with single or multiple random fibres under various loading conditions. • It can be used to optimise key design factors of FRC at meso/material scale by parametric studies. - Abstract: Fibre reinforced concrete (FRC) exhibits complicated failure modes such as fibre breakage, mortar cracking, crushing and spalling, fibre-mortar interfacial debonding, depending on many material properties, geometric dimensions, boundary and loading conditions. Most existing numerical models are unable to reproduce these failure modes that may occur simultaneously or sequentially in a specimen, mainly due to difficulties in generating finite element meshes with a large number of randomly-oriented fibres. Herein we develop a discrete-continuum coupled finite element modelling approach for FRC materials capable of effectively simulating all the major failure modes. The continuum damaged plasticity model is used to simulate damage and fracture behaviour of the mortar, while debonding of fibre-mortar interfaces is modelled by nonlinear cohesive interfacial elements. Unique techniques are devised to generate conforming meshes between fibres and the surrounding mortar so that the randomly-oriented fibres are easily modelled. The modelling approach is validated by simulating single fibre pullout tests with different inclination angles, notched and non-notched direct tensile tests and three-point bending beam tests with randomly-distributed multiple fibres.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Galerkin finite element method for time-fractional stochastic diffusion equations

In this paper, Galerkin finite element method for solving the time-fractional stochastic diffusion equations with multiplicative noise is proposed and investigated. The pathwise regularity properties of solutions to the semidiscrete Galerkin approximations are demonstrated and the convergence of optimal rates are derived. And also we construct the fully discrete scheme which is based on the approximations of the Mittag–Leffler function and analyze the error estimates of convergence in L{sub 2}-norm space. Finally, numerical results are conducted to confirm our theoretical findings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Hybrid Finite Element--Spectral Method for the Fractional Laplacian: Approximation Theory and Efficient Solver

Here, a numerical scheme is presented for approximating fractional order Poisson problems in two and three dimensions. The scheme is based on reformulating the original problem posed over $\Omega$ on the extruded domain $\mathcal{C}=\Omega\times[0,\infty)$ following. The resulting degenerate elliptic integer order PDE is then approximated using a hybrid FEM-spectral scheme. Finite elements are used in the direction parallel to the problem domain $\Omega$, and an appropriate spectral method is used in the extruded direction. The spectral part of the scheme requires that we approximate the true eigenvalues of the integer order Laplacian over $\Omega$. We derive an a priori error estimate which takes account of the error arising from using an approximation in place of the true eigenvalues. We further present a strategy for choosing approximations of the eigenvalues based on Weyl's law and finite element discretizations of the eigenvalue problem. The system of linear algebraic equations arising from the hybrid FEM-spectral scheme is decomposed into blocks which can be solved effectively using standard iterative solvers such as multigrid and conjugate gradient. Numerical examples in two and three dimensions suggest that the approach is quasi-optimal in terms of complexity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A single crystal plasticity finite element formulation with embedded deformation twins

Deformation twinning is an important plastic deformation mechanism in some polycrystalline metals such as titanium and magnesium. In this paper, we present a novel crystal plasticity finite element framework that accounts for deformation twinning explicitly, in addition to crystallographic slip. Within this computational framework, deformation twins are treated as weak discontinuities embedded within individual finite elements, such that a jump in the velocity gradient field is introduced (via the discretized gradient operator) between the twinned and untwinned crystalline regions, taking into account compatibility and traction continuity conditions at the interface between these two regions. The deformation gradient is multiplicatively split into elastic and plastic parts in the untwinned region, as is customary in finite-deformation crystal plasticity formulations. A different multiplicative decomposition of the deformation gradient into elastic, plastic (slip), and twinning parts is adopted in the twinned region, allowing deformation twinning to be accounted for as an additional mode of plastic deformation. Here, a stochastic model is used to predict twin nucleation at grain boundaries, and the evolution of the length and thickness of the twinned region under subsequent deformation is taken into account.

36 MATERIALS SCIENCE↗