Search NASA⌕ Search

SEARCH · Search NASA

Results for “Adaptive finite elements”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Adaptive quadrilateral and triangular finite-element scheme for compressible flows

The development of an adaptive mesh refinement procedure for analyzing high-speed compressible flows using the finite-element method is described. This new adaptation procedure, which uses both quadrilateral and triangular elements, was implemented with two explicit finite-element algorithms - the two-step Taylor-Galerkin and the multistep Galerkin-Runge-Kutta schemes. A von Neumann stability analysis and a rotating 'cosine hill'problem demonstrate the instability of the Taylor-Galerkin scheme when coupled with the adaptation procedure. For the same adaptive refinement scheme, the Galerkin-Runge-Kutta procedure yields stable solutions within its explicit stability limit. The utility of this new adaptation procedure for the prediction of compressible flow features is illustrated for inviscid problems involving strong shock interactions at hypersonic speeds.

Ramakrishnan, R.↗
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↗
NASA NTRS2019

Adjoint-Based Anisotropic Mesh Adaptation for a Stabilized Finite-Element Flow Solver

An adjoint solver is implemented in the FUN3D stabilized finite-element flow solver. The adjoint solution is used to generate anisotropic, adapted meshes to control error in scalar out- put functionals, such as lift or drag coefficient. The inviscid and turbulent adjoints are verified with a finite-difference derivative approximation and can be used in design optimization in addition to mesh adaptation. The adjoint capability represents an extension of previous research using the FUN3D stabilized finite-element capability with metric-based mesh adaptation and interpolation-based error estimates to generate highly anisotropic adapted meshes for turbulent flows. In the present work, a metric-based approach is again utilized, where the adjoint and the primal solutions both contribute to the generation of a metric tensor field that is subsequently used to produce the required anisotropic mesh for each adaptation cycle. Adaptive results are then shown for an inviscid supersonic flow over a diamond airfoil, inviscid transonic flow over the ONERA M6 wing, and viscous laminar flow over the NACA 0012 airfoil, all using drag force as the output functional. Adjoint-based adaptation is compared with a multiscale solution-based approach that controls the L norm of Mach number interpolation error to demonstrate the effectiveness and effciency of the adjoint-based adaptive mesh technology.

Balan, Aravind↗
NASA NTRS1989

Adaptive numerical solutions of the Euler equations in 3D using finite elements

The development of an adaptive mesh solution for a flow involving shock interaction on a swept cylinder and an initial solution for a flow past a complex fighter configuration is reported. The finite element solution algorithm, the mesh generation, and the adaptivity of the solution are described. Sample results for the flow past an F-18 configuration at Mach 0.9 and alpha of 3 deg and for shock interaction on a swept cylinder at Mach 8.04 are summarized.

Peraire, J.↗
NASA NTRS2005

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

The usage of shape memory materials has extended rapidly to many fields, including medical devices, actuators, composites, structures and MEMS devices. For these various applications, shape memory alloys (SMAs) are available in various forms: bulk, wire, ribbon, thin film, and porous. In this work, the focus is on SMA hybrid composites with adaptive-stiffening or morphing functions. These composites are created by using SMA ribbons or wires embedded in a polymeric based composite panel/beam. Adaptive stiffening or morphing is activated via selective resistance heating or uniform thermal loads. To simulate the thermomechanical behavior of these composites, a SMA model was implemented using ABAQUS user element interface and finite element simulations of the systems were studied. Several examples are presented which show that the implemented model can be a very useful design and simulation tool for SMA hybrid composites.

Gao, Xiu-Jie↗
NASA NTRS1985

Advances and trends in structures and dynamics; Proceedings of the Symposium, Washington, DC, October 22-25, 1984

Among the topics discussed are developments in structural engineering hardware and software, computation for fracture mechanics, trends in numerical analysis and parallel algorithms, mechanics of materials, advances in finite element methods, composite materials and structures, determinations of random motion and dynamic response, optimization theory, automotive tire modeling methods and contact problems, the damping and control of aircraft structures, and advanced structural applications. Specific topics covered include structural design expert systems, the evaluation of finite element system architectures, systolic arrays for finite element analyses, nonlinear finite element computations, hierarchical boundary elements, adaptive substructuring techniques in elastoplastic finite element analyses, automatic tracking of crack propagation, a theory of rate-dependent plasticity, the torsional stability of nonlinear eccentric structures, a computation method for fluid-structure interaction, the seismic analysis of three-dimensional soil-structure interaction, a stress analysis for a composite sandwich panel, toughness criterion identification for unidirectional composite laminates, the modeling of submerged cable dynamics, and damping synthesis for flexible spacecraft structures.

Noor, A. K.↗
NASA NTRS1986

An adaptive characteristic Petrov-Galerkin finite element method for convection-dominated linear and nonlinear parabolic problems in one space variable

The present adaptive FEM technique for convection-dominated problems is based on a Petrov-Galerkin scheme for spatial approximation, whose typical time-step employs test functions chosen to yield an approximate solution coinciding with the exact solutions at the finite element grid nodes. The derivation of truly local a posteriori error estimates is made possible by this procedure, which is also shown to be a very effective solver by the numerical examples presented.

Demkowicz, L.↗
NASA NTRS1988

An adaptive remeshing method for finite element thermal analysis

A finite element remeshing approach that makes use of quadrilateral and triangular elements is described. The approach uses the solution on a previous mesh to create a new mesh. Meshes are completely unstructured with highly refined elements in regions of steep gradients and larger elements where gradients are smaller. Studies of convergence rates for heat conduction problems with exact solutions show that for problems with highly localized solution variations, the remeshing approach gives smaller solution errors with fewer unknowns than refinement of uniform, structured meshes.

Thornton, Earl A.↗
NASA NTRS1990

Unstructured finite element mesh generation and adaptive procedures for CFD

A method is described for generating unstructured meshes of triangles or tetrahedra for computational domains of complex geometrical shape. To illustrate the power of the approach, it is applied to the solution of flows past several complete aircraft configurations. The advocated approach allows for the natural incorporation of mesh adaptivity and this is demonstrated for both inviscid and viscous computations in two and three dimensions.

Peraire, J.↗
NASA NTRS2004

Stress Distribution During Deformation of Polycrystalline Aluminum by Molecular-Dynamics and Finite-Element Modeling

In this paper, a multiscale modelling strategy is used to study the effect of grain-boundary sliding on stress localization in a polycrystalline microstructure with an uneven distribution of grain size. The development of the molecular dynamics (MD) analysis used to interrogate idealized grain microstructures with various types of grain boundaries and the multiscale modelling strategies for modelling large systems of grains is discussed. Both molecular-dynamics and finite-element (FE) simulations for idealized polycrystalline models of identical geometry are presented with the purpose of demonstrating the effectiveness of the adapted finite-element method using cohesive zone models to reproduce grain-boundary sliding and its effect on the stress distribution in a polycrystalline metal. The yield properties of the grain-boundary interface, used in the FE simulations, are extracted from a MD simulation on a bicrystal. The models allow for the study of the load transfer between adjacent grains of very different size through grain-boundary sliding during deformation. A large-scale FE simulation of 100 grains of a typical microstructure is then presented to reveal that the stress distribution due to grain-boundary sliding during uniform tensile strain can lead to stress localization of two to three times the background stress, thus suggesting a significant effect on the failure properties of the metal.

Yamakov, V.↗
NASA NTRS1999

Crack Turning and Arrest Mechanisms for Integral Structure

In the course of several years of research efforts to predict crack turning and flapping in aircraft fuselage structures and other problems related to crack turning, the 2nd order maximum tangential stress theory has been identified as the theory most capable of predicting the observed test results. This theory requires knowledge of a material specific characteristic length, and also a computation of the stress intensity factors and the T-stress, or second order term in the asymptotic stress field in the vicinity of the crack tip. A characteristic length, r(sub c), is proposed for ductile materials pertaining to the onset of plastic instability, as opposed to the void spacing theories espoused by previous investigators. For the plane stress case, an approximate estimate of r(sub c), is obtained from the asymptotic field for strain hardening materials given by Hutchinson, Rice and Rosengren (HRR). A previous study using of high order finite element methods to calculate T-stresses by contour integrals resulted in extremely high accuracy values obtained for selected test specimen geometries, and a theoretical error estimation parameter was defined. In the present study, it is shown that a large portion of the error in finite element computations of both K and T are systematic, and can be corrected after the initial solution if the finite element implementation utilizes a similar crack tip discretization scheme for all problems. This scheme is applied for two-dimensional problems to a both a p-version finite element code, showing that sufficiently accurate values of both K(sub I) and T can be obtained with fairly low order elements if correction is used. T-stress correction coefficients are also developed for the singular crack tip rosette utilized in the adaptive mesh finite element code FRANC2D, and shown to reduce the error in the computed T-stress significantly. Stress intensity factor correction was not attempted for FRANC2D because it employs a highly accurate quarter-point scheme to obtain stress intensity factors.

Pettit, Richard↗
DOE OSTI2016

Calculation of reinforced-concrete frame strength under a simultaneous static cross section load and a column lateral impact

When designing buildings with reinforced concrete that are planned to resist dynamic loads it is necessary to calculate this structural behavior under operational static and emergency impact and blast loads. Calculations of the structures under shock-wave loads can be performed by solving dynamic equations that do not consider static loads. Due to this fact the calculation of reinforced concrete frame under a simultaneous static and dynamic load in full 3d settings becomes a very non trivial and resource consuming problem. This problem can be split into two tasks. The first one is a shock-wave problem that can be solved using software package RANET-3, which allows solving the problem using finite elements method adapted for dynamic task. This method calculates strain-stress state of the material and its dynamic destruction, which is considered as growth and consolidation of micro defects under loading. On the second step the results of the first step are taken as input parameters for quasi static calculation of simultaneous static and dynamic load using finite elements method in AMP Civil Engineering-11.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Adaptive-Grid Methods for Phase Field Models of Microstructure Development

In this work the authors show how the phase field model can be solved in a computationally efficient manner that opens a new large-scale simulational window on solidification physics. Our method uses a finite element, adaptive-grid formulation, and exploits the fact that the phase and temperature fields vary significantly only near the interface. We illustrate how our method allows efficient simulation of phase-field models in very large systems, and verify the predictions of solvability theory at intermediate undercooling. We then present new results at low undercoolings that suggest that solvability theory may not give the correct tip speed in that regime. We model solidification using the phase-field model used by Karma and Rappel.

Provatas, Nikolas↗
DOE OSTI2013

Minimising the error in eigenvalue calculations involving the Boltzmann transport equation using goal-based adaptivity on unstructured meshes

This article presents a method for goal-based anisotropic adaptive methods for the finite element method applied to the Boltzmann transport equation. The neutron multiplication factor, k{sub eff}, is used as the goal of the adaptive procedure. The anisotropic adaptive algorithm requires error measures for k{sub eff} with directional dependence. General error estimators are derived for any given functional of the flux and applied to k{sub eff} to acquire the driving force for the adaptive procedure. The error estimators require the solution of an appropriately formed dual equation. Forward and dual error indicators are calculated by weighting the Hessian of each solution with the dual and forward residual respectively. The Hessian is used as an approximation of the interpolation error in the solution which gives rise to the directional dependence. The two indicators are combined to form a single error metric that is used to adapt the finite element mesh. The residual is approximated using a novel technique arising from the sub-grid scale finite element discretisation. Two adaptive routes are demonstrated: (i) a single mesh is used to solve all energy groups, and (ii) a different mesh is used to solve each energy group. The second method aims to capture the benefit from representing the flux from each energy group on a specifically optimised mesh. The k{sub eff} goal-based adaptive method was applied to three examples which illustrate the superior accuracy in criticality problems that can be obtained.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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↗
DOE OSTI2020

A matrix-free approach for finite-strain hyperelastic problems using geometric multigrid

This paper investigates matrix-free algorithms for problems in quasi-static finite-strain hyperelasticity. Iterative solvers with matrix-free operator evaluation have emerged as an attractive alternative to sparse matrices in the fluid dynamics and wave propagation communities because they significantly reduce the memory traffic, the limiting factor in classical finite element solvers. Specifically, we study different matrix-free realizations of the finite element tangent operator and determine whether generalized methods of incorporating complex constitutive behavior might be feasible. Furthermore, in order to improve the convergence behavior of iterative solvers, we also propose a method by which to construct level tangent operators and employ them to define a geometric multigrid preconditioner. Additionally, the performance of the matrix-free operator and the geometric multigrid preconditioner is compared to the matrix-based implementation with an algebraic multigrid (AMG) preconditioner on a single node for a representative numerical example of a heterogeneous hyperelastic material in two and three dimensions. We find that matrix-free methods for finite-strain solid mechanics are very promising, outperforming linear matrix-based schemes by two to five times, and that it is possible to develop numerically efficient implementations that are independent of the hyperelastic constitutive law.

42 ENGINEERING↗