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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2015

Analytical Modeling of a Novel Transverse Flux Machine for Direct Drive Wind Turbine Applications

This paper presents a nonlinear analytical model of a novel double sided flux concentrating Transverse Flux Machine (TFM) based on the Magnetic Equivalent Circuit (MEC) model. The analytical model uses a series-parallel combination of flux tubes to predict the flux paths through different parts of the machine including air gaps, permanent magnets (PM), stator, and rotor. The two-dimensional MEC model approximates the complex three-dimensional flux paths of the TFM and includes the effects of magnetic saturation. The model is capable of adapting to any geometry which makes it a good alternative for evaluating prospective designs of TFM as compared to finite element solvers which are numerically intensive and require more computation time. A single phase, 1 kW, 400 rpm machine is analytically modeled and its resulting flux distribution, no-load EMF and torque, verified with Finite Element Analysis (FEA). The results are found to be in agreement with less than 5% error, while reducing the computation time by 25 times.

17 WIND ENERGY↗

NASA NTRS1986

Finite element solver for 3-D compressible viscous flows

The space shuttle main engine (SSME) has extremely complex internal flow structure. The geometry of the flow domain is three-dimensional with complicated topology. The flow is compressible, viscous, and turbulent with large gradients in flow quantities and regions of recirculations. The analysis of the flow field in SSME involves several tedious steps. One is the geometrical modeling of the particular zone of the SSME being studied. Accessing the geometry definition, digitalizing it, and developing surface interpolations suitable for an interior grid generator require considerable amount of manual labor. There are several types of grid generators available with some general-purpose finite element programs. An efficient and robust computational scheme for solving 3D Navier-Stokes equations has to be implemented. Post processing software has to be adapted to visualize and analyze the computed 3D flow field. The progress made in a project to develop software for the analysis of the flow is discussed. The technical approach to the development of the finite element scheme and the relaxation procedure are discussed. The three dimensional finite element code for the compressible Navier-Stokes equations is listed.

Reddy, K. C.↗

NASA NTRS1994

An assessment of the adaptive unstructured tetrahedral grid, Euler Flow Solver Code FELISA

A three-dimensional solution-adaptive Euler flow solver for unstructured tetrahedral meshes is assessed, and the accuracy and efficiency of the method for predicting sonic boom pressure signatures about simple generic models are demonstrated. Comparison of computational and wind tunnel data and enhancement of numerical solutions by means of grid adaptivity are discussed. The mesh generation is based on the advancing front technique. The FELISA code consists of two solvers, the Taylor-Galerkin and the Runge-Kutta-Galerkin schemes, both of which are spacially discretized by the usual Galerkin weighted residual finite-element methods but with different explicit time-marching schemes to steady state. The solution-adaptive grid procedure is based on either remeshing or mesh refinement techniques. An alternative geometry adaptive procedure is also incorporated.

Djomehri, M. Jahed↗

NASA NTRS1996

Finite Element Analysis of Geodesically Stiffened Cylindrical Composite Shells Using a Layerwise Theory

Layerwise finite element analyses of geodesically stiffened cylindrical shells are presented. The layerwise laminate theory of Reddy (LWTR) is developed and adapted to circular cylindrical shells. The Ritz variational method is used to develop an analytical approach for studying the buckling of simply supported geodesically stiffened shells with discrete stiffeners. This method utilizes a Lagrange multiplier technique to attach the stiffeners to the shell. The development of the layerwise shells couples a one-dimensional finite element through the thickness with a Navier solution that satisfies the boundary conditions. The buckling results from the Ritz discrete analytical method are compared with smeared buckling results and with NASA Testbed finite element results. The development of layerwise shell and beam finite elements is presented and these elements are used to perform the displacement field, stress, and first-ply failure analyses. The layerwise shell elements are used to model the shell skin and the layerwise beam elements are used to model the stiffeners. This arrangement allows the beam stiffeners to be assembled directly into the global stiffness matrix. A series of analytical studies are made to compare the response of geodesically stiffened shells as a function of loading, shell geometry, shell radii, shell laminate thickness, stiffener height, and geometric nonlinearity. Comparisons of the structural response of geodesically stiffened shells, axial and ring stiffened shells, and unstiffened shells are provided. In addition, interlaminar stress results near the stiffener intersection are presented. First-ply failure analyses for geodesically stiffened shells utilizing the Tsai-Wu failure criterion are presented for a few selected cases.

Gerhard, Craig Steven↗

DOE OSTI2014

The use of discrete harmonics in direct multi-scale embedding of polycrystal plasticity

In this work, we describe an approach for directly embedding polycrystal plasticity models in component scale calculations, with an emphasis on computational tractability. Previously, we have employed adaptive sampling to mitigate the computational cost of direct embedding, achieving two or more orders of magnitude in wall-clock speedup compared to more traditional approaches. However, in our previous work the crystal orientation distribution function (crystallographic texture) was not allowed to evolve significantly. Here we discuss an approach that allows for evolving texture by employing discrete harmonics, effectively decoupling considerations related to accuracy of integrals in the homogenization from those related to adequate representation of the evolving texture. We discuss the basic behaviors and convergence of the new polycrystal plasticity framework. Specific applications focus on the deformation of titanium, including the effects of twinning. Overall, the discrete harmonic based framework offers an attractive path forward for computationally efficient multi-scale embedding of polycrystal plasticity.

36 MATERIALS SCIENCE↗

DOE OSTI2014

Percept

The Percept software package is a collection of libraries and executables that provide tools for verifying computer simulations of engineering components and systems. Percept is useful for simulations using the finite element or finite volume methods on unstructured meshes. Percept includes API's for adaptive mesh refinement, geometry representation, the method of manufactured solutions, analysis of convergence including the convergence of vibrational eigenmodes, and metrics for analyzing the difference between fields represented on two different overlapping unstructured grids.

Carnes, Brian↗

NASA NTRS1985

Probabilistic Finite Element: Variational Theory

The goal of this research is to provide techniques which are cost-effective and enable the engineer to evaluate the effect of uncertainties in complex finite element models. Embedding the probabilistic aspects in a variational formulation is a natural approach. In addition, a variational approach to probabilistic finite elements enables it to be incorporated within standard finite element methodologies. Therefore, once the procedures are developed, they can easily be adapted to existing general purpose programs. Furthermore, the variational basis for these methods enables them to be adapted to a wide variety of structural elements and to provide a consistent basis for incorporating probabilistic features in many aspects of the structural problem. Tasks concluded include the theoretical development of probabilistic variational equations for structural dynamics, the development of efficient numerical algorithms for probabilistic sensitivity displacement and stress analysis, and integration of methodologies into a pilot computer code.

Belytschko, T.↗

NASA NTRS1988

Finite element (MARC) solution technologies for viscoplastic analyses

A need for development of realistic constitutive models for structural components operating at high temperatures, accompanied by appropriate solution technologies for stress/life analyses of these components is studied. Viscoplastic models provide a better description of inelastic behavior of materials, but their mathematical structure is very complex. The highly nonlinear and stiff nature of the constitutive equations makes analytical solutions difficult. Therefore, suitable solution, finite element or other numerical, technologies must be developed to make these models adaptable for better and rational designs of components. NASA-Lewis has developed several solution technologies and successfully applied them to the solution of a number of uniaxial and multiaxial problems. Some of these solution technologies are described along with the models and representative results. The solution technologies developed and presented encompass a wide range of models, such as, isotropic, anisotropic, metal matrix composites, and single crystal models.

Arya, V. K.↗

NASA NTRS1987

Isoparametric finite element methods to estimate discontinuous coefficients in two-dimensional elliptic equations

A summary of findings is presented on the problem of estimating discontinuous coefficients, including the locations of interfaces at which discontinuities occur, for elliptic equations on a two-dimensional domain. Finite-element approximations are constructed that are based on an underlying grid that adaptively moves to follow locations of interfaces as they are iteratively updated in the numerical estimation process.

Lamm, Patricia K.↗

NASA NTRS1990

Effects of turbulence mixing, variable properties, and vaporization on spray droplet combustion

Combustion of liquid fuels in the form of spray droplets is simulated numerically. Various vaporization models are examined as to their performance in finite element calculations involving a turbulent flow field. The Eulerian coordinate for the gas and Lagrangian coordinate for the liquid spray droplets are coupled through source terms being updated in the equations of continuity, momentum, and energy. The k-epsilon and modified eddy breakup models are used for simulating turbulent spray combustion flow field. Numerical results for the droplet trajectories, droplet heating, recirculation characteristics, and effects of evaporation models are evaluated. It is also shown that the finite element method is advantageous in dealing with complex geometries, complex boundary conditions, adaptive unstructured grids.

Kim, Y. M.↗

NASA NTRS2006

Strip Yield Model Numerical Application to Different Geometries and Loading Conditions

A new numerical method based on the strip-yield analysis approach was developed for calculating the Crack Tip Opening Displacement (CTOD). This approach can be applied for different crack configurations having infinite and finite geometries, and arbitrary applied loading conditions. The new technique adapts the boundary element / dislocation density method to obtain crack-face opening displacements at any point on a crack, and succeeds by obtaining requisite values as a series of definite integrals, the functional parts of each being evaluated exactly in a closed form.

Hatamleh, Omar↗

DOE OSTI2014

Percept

Source record↗

DOE OSTI2022

The Schwarz Alternating Method for the Seamless Coupling of Nonlinear Reduced Order Models and Full Order Models

Projection-based model order reduction allows for the parsimonious representation of full order models (FOMs), typically obtained through the discretization of a set of partial differential equations (PDEs) using conventional techniques (e.g., finite element, finite volume, finite difference methods) where the discretization may contain a very large number of degrees of freedom. As a result of this more compact representation, the resulting projection-based reduced order models (ROMs) can achieve considerable computational speedups, which are especially useful in real-time or multi-query analyses. One known deficiency of projection-based ROMs is that they can suffer from a lack of robustness, stability and accuracy, especially in the predictive regime, which ultimately limits their useful application. Another research gap that has prevented the widespread adoption of ROMs within the modeling and simulation community is the lack of theoretical and algorithmic foundations necessary for the “plug-and-play” integration of these models into existing multi-scale and multi-physics frameworks. This paper describes a new methodology that has the potential to address both of the aforementioned deficiencies by coupling projection-based ROMs with each other as well as with conventional FOMs by means of the Schwarz alternating method [41]. Leveraging recent work that adapted the Schwarz alternating method to enable consistent and concurrent multiscale coupling of finite element FOMs in solid mechanics [35, 36], we present a new extension of the Schwarz framework that enables FOM-ROM and ROM-ROM coupling, following a domain decomposition of the physical geometry on which a PDE is posed. In order to maintain efficiency and achieve computation speed-ups, we employ hyper-reduction via the Energy-Conserving Sampling and Weighting (ECSW) approach [13]. We evaluate the proposed coupling approach in the reproductive as well as in the predictive regime on a canonical test case that involves the dynamic propagation of a traveling wave in a nonlinear hyper-elastic material.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1985

Finite element methods for high speed flows

An explicit finite element based solution procedure for solving the equations of compressible viscous high speed flow is presented. The method uses domain splitting to advance the solution with different timesteps on different portions of the mesh. For steady inviscid flows, adaptive mesh refinement procedures are successfully employed to enhance the definition of discontinuities. Preliminary ideas on the application of adaptive mesh refinement to the solution of problems involving steady viscous flow are presented. Sample timings are given for the performance of the finite element code on modern supercomputers.

Loehner, R.↗

NASA NTRS1997

COMET-AR User's Manual: COmputational MEchanics Testbed with Adaptive Refinement

The COMET-AR User's Manual provides a reference manual for the Computational Structural Mechanics Testbed with Adaptive Refinement (COMET-AR), a software system developed jointly by Lockheed Palo Alto Research Laboratory and NASA Langley Research Center under contract NAS1-18444. The COMET-AR system is an extended version of an earlier finite element based structural analysis system called COMET, also developed by Lockheed and NASA. The primary extensions are the adaptive mesh refinement capabilities and a new "object-like" database interface that makes COMET-AR easier to extend further. This User's Manual provides a detailed description of the user interface to COMET-AR from the viewpoint of a structural analyst.

Moas, E.↗

NASA NTRS1986

Cyclic creep analysis from elastic finite-element solutions

A uniaxial approach was developed for calculating cyclic creep and stress relaxation at the critical location of a structure subjected to cyclic thermomechanical loading. This approach was incorporated into a simplified analytical procedure for predicting the stress-strain history at a crack initiation site for life prediction purposes. An elastic finite-element solution for the problem was used as input for the simplified procedure. The creep analysis includes a self-adaptive time incrementing scheme. Cumulative creep is the sum of the initial creep, the recovery from the stress relaxation and the incremental creep. The simplified analysis was exercised for four cases involving a benchmark notched plate problem. Comparisons were made with elastic-plastic-creep solutions for these cases using the MARC nonlinear finite-element computer code.

Kaufman, A.↗

DOE OSTI2019

A dual mesh method with adaptivity for stress-constrained topology optimization

This paper is concerned with the topological optimization of elastic structures, with the goal of minimizing the compliance and/or mass of the structure, subject to a stress constraint. It is well known that depending upon the geometry and the loading conditions, the stress field can exhibit singularities, if these singularities are not adequately resolved, the topological optimization process will be ineffective. For computational efficiency, adaptive mesh refinement is required to adequately resolve the stress field. This poses a challenge for the traditional Solid Isotropic Material with Penalization (SIMP) method that employs a one-to-one correspondence between the finite element mesh and the optimization design variables because as the mesh is refined, the optimization process must somehow be re-started with a new set of design variables. The proposed solution here is a dual mesh approach, one finite element mesh defines the material distribution, a second finite element mesh is used for the computation of the displacement and stress fields. This allows for stress-based adaptive mesh refinement without modifying the definition of the optimization design variables. A second benefit of this dual mesh approach is that there is no need to apply a filter to the design variables to enforce a length scale, instead the length scale is determined by the design mesh. This reduces the number of design variables and allows the designer to apply spatially varying length scale if desired. The efficacy of this dual mesh approach is established via several stress-constrained topology optimization problems.

36 MATERIALS SCIENCE↗

NASA NTRS2013

ISSM: Ice Sheet System Model

In order to have the capability to use satellite data from its own missions to inform future sea-level rise projections, JPL needed a full-fledged ice-sheet/iceshelf flow model, capable of modeling the mass balance of Antarctica and Greenland into the near future. ISSM was developed with such a goal in mind, as a massively parallelized, multi-purpose finite-element framework dedicated to ice-sheet modeling. ISSM features unstructured meshes (Tria in 2D, and Penta in 3D) along with corresponding finite elements for both types of meshes. Each finite element can carry out diagnostic, prognostic, transient, thermal 3D, surface, and bed slope simulations. Anisotropic meshing enables adaptation of meshes to a certain metric, and the 2D Shelfy-Stream, 3D Blatter/Pattyn, and 3D Full-Stokes formulations capture the bulk of the ice-flow physics. These elements can be coupled together, based on the Arlequin method, so that on a large scale model such as Antarctica, each type of finite element is used in the most efficient manner. For each finite element referenced above, ISSM implements an adjoint. This adjoint can be used to carry out model inversions of unknown model parameters, typically ice rheology and basal drag at the ice/bedrock interface, using a metric such as the observed InSAR surface velocity. This data assimilation capability is crucial to allow spinning up of ice flow models using available satellite data. ISSM relies on the PETSc library for its vectors, matrices, and solvers. This allows ISSM to run efficiently on any parallel platform, whether shared or distrib- ISSM: Ice Sheet System Model NASA's Jet Propulsion Laboratory, Pasadena, California uted. It can run on the largest clusters, and is fully scalable. This allows ISSM to tackle models the size of continents. ISSM is embedded into MATLAB and Python, both open scientific platforms. This improves its outreach within the science community. It is entirely written in C/C++, which gives it flexibility in its design, and the power/speed that C/C++ allows. ISSM is svn (subversion) hosted, on a JPL repository, to facilitate its development and maintenance. ISSM can also model propagation of rifts using contact mechanics and mesh splitting, and can interface to the Dakota software. To carry out sensitivity analysis, mesh partitioning algorithms are available, based on the Scotch, Chaco, and Metis partitioners that ensure equal area mesh partitions can be done, which are then usable for sampling and local reliability methods.

Larour, Eric↗