Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous 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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Elastic wave suppression through additively manufactured petal lattice metamaterials

Lattice-core sandwich structure metamaterials are lightweight alternatives to monolithic materials that can present better mechanical, thermal, and energy dampening performance. Manufacturing lattice metamaterials to follow curved surfaces can pose a challenge, as the lattices rely on their geometric orientation to the substrate for their mechanical properties. This work rationally designed a lattice structure where the surface is broken up into “petals” connected to the underlying lattice, which localizes the petals’ impact response. This design opens a pathway for implementation of lattice-core sandwich structures onto complex surface geometries. These petal structures were evaluated for their energy absorption efficiency experimentally by utilizing pressure waves generated with nanosecond lasers and computationally via finite element modeling. The lattice structures exhibited a two-orders-of-magnitude decrease in transmitted pressure compared to their constituent steel at equivalent mass. Furthermore, localizing energy absorption into petal structures provided a 44% reduction in peak load compared to a continuous “single-petal” design.

Sarracino, Alex↗
NASA NTRS1990

Micromechanics analysis of space simulated thermal deformations and stresses in continuous fiber reinforced composites

Space simulated thermally induced deformations and stresses in continuous fiber reinforced composites were investigated with a micromechanics analysis. The investigation focused on two primary areas. First, available explicit expressions for predicting the effective coefficients of thermal expansion (CTEs) for a composite were compared with each other, and with a finite element (FE) analysis, developed specifically for this study. Analytical comparisons were made for a wide range of fiber/matrix systems, and predicted values were compared with experimental data. The second area of investigation focused on the determination of thermally induced stress fields in the individual constituents. Stresses predicted from the FE analysis were compared to those predicted from a closed-form solution to the composite cylinder (CC) model, for two carbon fiber/epoxy composites. A global-local formulation, combining laminated plate theory and FE analysis, was used to determine the stresses in multidirectional laminates. Thermally induced damage initiation predictions were also made.

Bowles, David E.↗
NASA NTRS1996

First-Order System Least Squares for the Stokes Equations, with Application to Linear Elasticity

Following our earlier work on general second-order scalar equations, here we develop a least-squares functional for the two- and three-dimensional Stokes equations, generalized slightly by allowing a pressure term in the continuity equation. By introducing a velocity flux variable and associated curl and trace equations, we are able to establish ellipticity in an H(exp 1) product norm appropriately weighted by the Reynolds number. This immediately yields optimal discretization error estimates for finite element spaces in this norm and optimal algebraic convergence estimates for multiplicative and additive multigrid methods applied to the resulting discrete systems. Both estimates are uniform in the Reynolds number. Moreover, our pressure-perturbed form of the generalized Stokes equations allows us to develop an analogous result for the Dirichlet problem for linear elasticity with estimates that are uniform in the Lame constants.

Cai, Z.↗
NASA NTRS2003

Simulating Fatigue Crack Growth in Spiral Bevel Pinion

This project investigates computational modeling of fatigue crack growth in spiral bevel gears. Current work is a continuation of the previous efforts made to use the Boundary Element Method (BEM) to simulate tooth-bending fatigue failure in spiral bevel gears. This report summarizes new results predicting crack trajectory and fatigue life for a spiral bevel pinion using the Finite Element Method (FEM). Predicting crack trajectories is important in determining the failure mode of a gear. Cracks propagating through the rim may result in catastrophic failure, whereas the gear may remain intact if one tooth fails and this may allow for early detection of failure. Being able to predict crack trajectories is insightful for the designer. However, predicting growth of three-dimensional arbitrary cracks is complicated due to the difficulty of creating three-dimensional models, the computing power required, and absence of closed- form solutions of the problem. Another focus of this project was performing three-dimensional contact analysis of a spiral bevel gear set incorporating cracks. These analyses were significant in determining the influence of change of tooth flexibility due to crack growth on the magnitude and location of contact loads. This is an important concern since change in contact loads might lead to differences in SIFs and therefore result in alteration of the crack trajectory. Contact analyses performed in this report showed the expected trend of decreasing tooth loads carried by the cracked tooth with increasing crack length. Decrease in tooth loads lead to differences between SIFs extracted from finite element contact analysis and finite element analysis with Hertz contact loads. This effect became more pronounced as the crack grew.

Ural, Ani↗
DOE OSTI2021

An efficient reconstruction algorithm for diffusion on triangular grids using the nodal discontinuous Galerkin method

High-energy-density (HED) hydrodynamics studies such as those relevant to inertial confinement fusion and astrophysics require highly disparate densities, temperatures, viscosities, and other diffusion parameters over relatively short spatial scales. This presents a challenge for high-order accurate methods to effectively resolve the hydrodynamics at these scales, particularly in the presence of highly disparate diffusion. A significant volume of engineering and physics applications use an unstructured discontinuous Galerkin (DG) method developed based on the finite element mesh generation and algorithmic framework. This work discusses the application of an affine reconstructed nodal DG method for unstructured grids of triangles. Solving the diffusion terms in the DG method is non-trivial due to the solution representations being piecewise continuous. Hence, the diffusive flux is not defined on the interface of elements. The proposed numerical approach reconstructs a smooth solution in a parallelogram that is enclosed by the quadrilateral formed by two adjacent triangle elements. The interface between these two triangles is the diagonal of the enclosed parallelogram. Similar to triangles, the mapping of parallelograms from a physical domain to a reference domain is an affine mapping, which is necessary for an accurate and efficient implementation of the numerical algorithm. Thus, all computations can still be performed on the reference domain, which promotes efficiency in computation and storage. This reconstruction does not make assumptions on choice of polynomial basis. Reconstructed DG algorithms have previously been developed for modal implementations of the convection–diffusion equations. However, to the best of the authors’ knowledge, this is the first practical guideline that has been proposed for applying the reconstructed algorithm on a nodal discontinuous Galerkin method with a focus on accuracy and efficiency. As a result, the algorithm is demonstrated on a number of benchmark cases as well as a challenging substantive problem in HED hydrodynamics with highly disparate diffusion parameters.

Computational efficiency↗
DOE OSTI2025

Multi-phase Wireless Power Transfer with High Power Density Inductive Coils for Electric Drone Charging

Wireless charging of unmanned ground vehi-cles and aircraft has been proposed to increase charging reliability and security, allow for autonomous functionality, and either reduce battery size or increase continuous flight time. This paper proposes a three-phase Litz wire primary and a two-phase PCB secondary for high secondary-side power density considering misalignment tolerances, surface and volumetric power density, and coil sizing. Elec-tromagnetic 3D finite element analysis (FEA) simulations are conducted to study variation in mutual inductance and coupling coefficient with different secondary coil sizes and number of turns, horizontal and vertical misalignment between the primary and secondary, and a combination of both for a set of designs with the same rated power. Results indicate that receiver coils of reduced size relative to the primary can deliver rated power for increased power density and higher tolerance to misalignment.

Gastineau, Lucas [University of Kentucky]↗
NASA NTRS1996

Coupled mixed-field laminate theory and finite element for smart piezoelectric composite shell structures

Mechanics for the analysis of laminated composite shells with piezoelectric actuators and sensors are presented. A new mixed-field laminate theory for piezoelectric shells is formulated in curvilinear coordinates which combines single-layer assumptions for the displacements and a layerwise representation for the electric potential. The resultant coupled governing equations for curvilinear piezoelectric laminates are described. Structural mechanics are subsequently developed and an 8-node finite-element is formulated for the static and dynamic analysis of adaptive composite structures of general laminations containing piezoelectric layers. Evaluations of the method and comparisons with reported results are presented for laminated piezoelectric-composite plates, a closed cylindrical shell with a continuous piezoceramic layer and a laminated composite semi-circular cantilever shell with discrete cylindrical piezoelectric actuators and/or sensors.

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

Stirling Technology Development at NASA GRC

The Department of Energy, Stirling Technology Company (STC), and NASA Glenn Research Center (NASA Glenn) are developing a free-piston Stirling convertor for a high efficiency Stirling Radioisotope Generator (SRG) for NASA Space Science missions. The SRG is being developed for multimission use, including providing electric power for unmanned Mars rovers and deep space missions. NASA Glenn is conducting an in-house technology project to assist in developing the convertor for space qualification and mission implementation. Recent testing of 55-We Technology Demonstration Convertors (TDCs) built by STC includes mapping of a second pair of TDCs, single TDC testing, and TDC electromagnetic interference and electromagnetic compatibility characterization on a nonmagnetic test stand. Launch environment tests of a single TDC without its pressure vessel to better understand the convertor internal structural dynamics and of dual-opposed TDCs with several engineering mounting structures with different natural frequencies have recently been completed. A preliminary life assessment has been completed for the TDC heater head, and creep testing of the IN718 material to be used for the flight convertors is underway. Long-term magnet aging tests are continuing to characterize any potential aging in the strength or demagnetization resistance of the magnets used in the linear alternator (LA). Evaluations are now beginning on key organic materials used in the LA and piston/rod surface coatings. NASA Glenn is also conducting finite element analyses for the LA, in part to look at the demagnetization margin on the permanent magnets. The world's first known integrated test of a dynamic power system with electric propulsion was achieved at NASA Glenn when a Hall-effect thruster was successfully operated with a free-piston Stirling power source. Cleveland State University is developing a multidimensional Stirling computational fluid dynamics code to significantly improve Stirling loss predictions and assist in identifying convertor areas for further improvements. This paper will update the status and results for these efforts.

Thieme, Lanny G.↗
NASA NTRS2002

Stirling Technology Development at NASA GRC

The Department of Energy, Stirling Technology Company (STC), and NASA Glenn Research Center (NASA Glenn) are developing a free-piston Stirling convertor for a high-efficiency Stirling Radioisotope Generator (SRG) for NASA Space Science missions. The SRG is being developed for multimission use, including providing electric power for unmanned Mars rovers and deep space missions. NASA Glenn is conducting an in-house technology project to assist in developing the convertor for space qualification and mission implementation. Recent testing, of 55-We Technology Demonstration Convertors (TDC's) built by STC includes mapping, of a second pair of TDC's, single TDC testing, and TDC electromagnetic interference and electromagnetic compatibility characterization on a nonmagnetic test stand. Launch environment tests of a single TDC without its pressure vessel to better understand the convertor internal structural dynamics and of dual-opposed TDC's with several engineering mounting structures with different natural frequencies have recently been completed. A preliminary life assessment has been completed for the TDC heater head, and creep testing of the IN718 material to be used for the flight convertors is underway. Long-term magnet aging tests are continuing to characterize any potential aging in the strength or demagnetization resistance of the magnets used in the linear alternator (LA). Evaluations are now beginning on key organic materials used in the LA and piston/rod surface coatings. NASA Glenn is also conducting finite element analyses for the LA, in part to look at the demagnetization margin on the permanent magnets. The world's first known integrated test of a dynamic power system with electric propulsion was achieved at NASA Glenn when a Hall-effect thruster was successfully operated with a free-piston Stirling power source. Cleveland State University is developing a multidimensional Stirling computational fluid dynamics code to significantly improve Stirling loss predictions and assist in identifying convertor areas for further improvements. This paper will update the status and results for these efforts.

Thieme, Lanny G.↗
NASA NTRS1996

Two-Dimensional Dynamic Simulation of a Continuous Foil Bearing

In this paper, the two dimensional(radial and circumferential) transient Navier-Stokes equations are used to solve the hydrodynamic problem in conjunction with the time dependent motion of the journal, and the deformable, spring supported foil. The elastic deformation of the foil and its supports are simulated by a finite element model. The time-dependent Navier-Stokes formulation is used to solve for the interaction between the fluid lubricant, the motion of the journal and the deformable foil boundary. The steady state, the quasi-transient and the full transient dynamic simulation of the foil-fluid journal interaction are examined on a comparative basis. For the steady state simulation, the fluid lubricant pressures are evaluated for a particular journal position, by means of an iterative scheme until convergence is achieved in both the fluid pressures and the corresponding foil deformation. For the quasi-transient case, the transient motion of the journal is calculated using a numerical integration scheme for the velocity and displacement of the journal. The deformation of the foil is evaluated through numerical iteration in feedback mode with the fluid film pressure generated by the journal motion until convergence at every time step is achieved. For the full transient simulation, a parallel real-time integration scheme is used to evaluate simultaneously the new journal position and the new deformed shape of the foil at each time step. The pressure of the fluid lubricant is iterated jointly with the corresponding journal position and the deformed foil geometry until convergence is achieved. A variable time-stepping Newmark-Beta integration procedure is used to evaluate the transient dynamics at each time step of the bearing.

Braun, M. Jack↗
DOE OSTI2021

Vibration-Based Sensor Design: A Grey-Box Approach

Knowledge of the internal structure of an object or device under investigation proceeds from the basic idea of constructing its dynamic behavioral relations governed by a set of differential/algebraic equations that characterize its response. These equations can be partial differential equations leading to finite element or finite difference relations requiring a complex numerical solution on a super computer or ordinary differential equations requiring sophisticated numerical integration techniques to obtain the desired solution. Discrete dynamic systems evolving from digitized data acquisition are typically captured by sampled-data (continuous-to-discrete) representations characterized by a set of difference equations specifying the underlying system dynamics. In any case, with a mathematical description in hand, Grey-Box modeling techniques have evolved, concerned with the estimation of model parameters embedded in a prescribed set of equations (the system) governing its behavior, while capturing the underlying physical phenomenology of the problem at hand.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Development of 3D electromagnetic modeling tools for airborne vehicles

The main goal of this report is to advance the development of methodologies for scattering by airborne composite vehicles. Although the primary focus continues to be the development of a general purpose computer code for analyzing the entire structure as a single unit, a number of other tasks are also being pursued in parallel with this effort. One of these tasks discussed within is on new finite element formulations and mesh termination schemes. The goal here is to decrease computation time while retaining accuracy and geometric adaptability.The second task focuses on the application of wavelets to electromagnetics. Wavelet transformations are shown to be able to reduce a full matrix to a band matrix, thereby reducing the solutions memory requirements. Included within this document are two separate papers on finite element formulations and wavelets.

Volakis, John L.↗
NASA NTRS2013

Hypersonic Inflatable Aerodynamic Decelerator (HIAD) Technology Development Overview

The successful flight of the Inflatable Reentry Vehicle Experiment (IRVE)-3 has further demonstrated the potential value of Hypersonic Inflatable Aerodynamic Decelerator (HIAD) technology. This technology development effort is funded by NASA's Space Technology Mission Directorate (STMD) Game Changing Development Program (GCDP). This paper provides an overview of a multi-year HIAD technology development effort, detailing the projects completed to date and the additional testing planned for the future. The effort was divided into three areas: Flexible Systems Development (FSD), Mission Advanced Entry Concepts (AEC), and Flight Validation. FSD consists of a Flexible Thermal Protection Systems (FTPS) element, which is investigating high temperature materials, coatings, and additives for use in the bladder, insulator, and heat shield layers; and an Inflatable Structures (IS) element which includes manufacture and testing (laboratory and wind tunnel) of inflatable structures and their associated structural elements. AEC consists of the Mission Applications element developing concepts (including payload interfaces) for missions at multiple destinations for the purpose of demonstrating the benefits and need for the HIAD technology as well as the Next Generation Subsystems element. Ground test development has been pursued in parallel with the Flight Validation IRVE-3 flight test. A larger scale (6m diameter) HIAD inflatable structure was constructed and aerodynamically tested in the National Full-scale Aerodynamics Complex (NFAC) 40ft by 80ft test section along with a duplicate of the IRVE-3 3m article. Both the 6m and 3m articles were tested with instrumented aerodynamic covers which incorporated an array of pressure taps to capture surface pressure distribution to validate Computational Fluid Dynamics (CFD) model predictions of surface pressure distribution. The 3m article also had a duplicate IRVE-3 Thermal Protection System (TPS) to test in addition to testing with the Aerocover configuration. Both the Aerocovers and the TPS were populated with high contrast targets so that photogrammetric solutions of the loaded surface could be created. These solutions both refined the aerodynamic shape for CFD modeling and provided a deformed shape to validate structural Finite Element Analysis (FEA) models. Extensive aerothermal testing has been performed on the TPS candidates. This testing has been conducted in several facilities across the country. The majority of the testing has been conducted in the Boeing Large Core Arc Tunnel (LCAT). HIAD is continuing to mature testing methodology in this facility and is developing new test sample fixtures and control methodologies to improve understanding and quality of the environments to which the samples are subjected. Additional testing has been and continues to be performed in the NASA LaRC 8ft High Temperature Tunnel, where samples up to 2ft by 2ft are being tested over representative underlying structures incorporating construction features such as sewn seams and through-thickness quilting. With the successful completion to the IRVE-3 flight demonstration, mission planning efforts are ramping up on the development of the HIAD Earth Atmospheric Reenty Test (HEART) which will demonstrate a relevant scale vehicle in relevant environments via a large-scale aeroshell (approximately 8.5m) entering at orbital velocity (approximately 7km/sec) with an entry mass on the order of 4MT. Also, the Build to Print (BTP) hardware built as a risk mitigation for the IRVE-3 project to have a "spare" ready to go in the event of a launch vehicle delivery failure is now available for an additional sub-orbital flight experiment. Mission planning is underway to define a mission that can utilize this existing hardware and help the HIAD project further mature this technology.

Hughes, Stephen J.↗
DOE OSTI2020

Hydro-mechanical coupled modeling of hydraulic fracturing using the mesh fragmentation technique

The inclusion of special interface elements between standard finite elements is becoming a common technique to model discontinuities in porous media using the finite element method (FEM). Recently, it was proved that high aspect ratio (HAR) elements can reproduce very satisfactorily the main features of behavior observed in a single hydraulically-induced fracture. In this paper, the mesh fragmentation technique (MFT) is proposed to extend the application of HAR elements to simulate more general problems involving the formation of hydraulic fractures in rocks by introducing this type of element in-between the standard FE of a typical mesh. The proposed framework makes use of standard FE techniques. Therefore, it is relatively easy to upgrade an existing FE program for continuous porous media to deal with evolving discontinuities. The proposed approach is verified and validated against available analytical and numerical solutions. The crack patterns predicted by the MFT are also compared with the ones observed in the laboratory. The formation of multiple fractures is also analyzed. In all cases, the proposed technique was able to properly capture the main features of rock behavior subjected to hydraulic fracturing.

58 GEOSCIENCES↗
NASA NTRS2018

An Efficient Inverse Finite Element Method for Shape and Stress Sensing of Laminated Composite and Sandwich Plates and Shells

The inverse Finite Element Method (iFEM) is a revolutionary methodology for real-time reconstruction of full-field structural displacements and stresses in structures that are instrumented with strain sensors. This inverse problem is commonly referred to as shape and stress sensing, which is well-recognized as an enabling technology for structural health monitoring systems. In this study, an improved iFEM formulation is proposed for shape and stress sensing of laminated composite and sandwich plates and shells. The formulation includes the kinematics of a shear deformation plate theory known as Refined Zigzag Theory (RZT) as its baseline. The present iFEM formulation is based upon the minimization of a weighted-least-squares functional that uses the complete set of section strains of RZT. The improved iFEM methodology is applicable for shape and stress sensing of thin and moderately thick plate and shell structures involving a relatively small number of strain gauges. The main advantage of the current formulation is that highly accurate through-the-thickness distributions of displacements, strains, and stresses are attainable using an element based on simple C0-continuous displacement interpolation functions. A three-node inverseshell element, named i3-RZT, is developed. Two example problems are examined in detail: (1) a simply supported rectangular laminated composite plate and (2) a wedge structure with a hole near one of the clamped ends. For both problems, the experimental strain data are generated numerically by the direct finite element analysis using high-fidelity discretizations. These strains are then regarded as the experimental strains obtained from surface mounted strain gauges or embedded fiber Bragg grating (FBG) sensors. The numerical results demonstrate the superior capability and potential applicability of the i3- RZT/iFEM methodology for performing accurate shape and stress sensing of complex composite structures.

Kefal, Adnan↗
DOE OSTI2023

FINITE ELEMENT MODEL MESH REFINEMENT EFFECTS ON QUALIFICATION OF NUCLEAR GRADE GRAPHITE COMPONENTS

The American Society of Mechanical Engineers (ASME) provides the full and simplified design-by-analysis probabilistic assessments for determining acceptance of nuclear grade graphite core components. The assessments can be characterized by three parts: (1) a component stress distribution, often determined by a finite element (FE) model; (2) a Weibull probability density function (pdf) that characterizes the experimental tensile strength distribution; and (3) the post-processor, which combines the FE model and the Weibull strength distribution in accordance with the full and simplified assessments to determine component acceptance. It is known that the level of mesh refinement in FE models can affect the modeled component’s calculated stress distribution. Depending on the component geometry, the stress distribution may converge with sufficient refinement. It was previously unknown whether the acceptance decision resulting from the full and simplified assessments might change even with sufficient mesh refinement. This study explores that question using experimental strength results for a dog-bone geometry for two graphite grades, IG-110 and PCEA. The simplified assessment has two criteria that must be met, the first limits the combined membrane stress by the allowable stress and the second limits the peak equivalent stress by the allowable stress scaled by the ratio of flexural to tensile strength. In the application of the simplified assessment, convergence of the peak equivalent stress required extreme mesh refinement, however, the acceptance decision was not affected. It is hypothesized that more complex geometries with stress concentrations may present mesh refinement effects on the simplified assessment acceptance decision. Mesh refinement did affect the acceptance decision in the full assessment for the applied pressure loadings in this study. This work suggests component stress distribution convergence is not a sufficient criteria for POF convergence in the full assessment and that mesh refinement should continue until the POF has converged, especially where the resulting POF is bordering the SRC acceptable POF limit.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2006

A Computational and Experimental Investigation of Shear Coaxial Jet Atomization

The instability and subsequent atomization of a viscous liquid jet emanated into a high-pressure gaseous surrounding is studied both computationally and experimentally. Liquid water issued into nitrogen gas at elevated pressures is used to simulate the flow conditions in a coaxial shear injector element relevant to liquid propellant rocket engines. The theoretical analysis is based on a simplified mathematical formulation of the continuity and momentum equations in their conservative form. Numerical solutions of the governing equations subject to appropriate initial and boundary conditions are obtained via a robust finite difference scheme. The computations yield real-time evolution and subsequent breakup characteristics of the liquid jet. The experimental investigation utilizes a digital imaging technique to measure resultant drop sizes. Data were collected for liquid Reynolds number between 2,500 and 25,000, aerodynamic Weber number range of 50-500 and ambient gas pressures from 150 to 1200 psia. Comparison of the model predictions and experimental data for drop sizes at gas pressures of 150 and 300 psia reveal satisfactory agreement particularly for lower values of investigated Weber number. The present model is intended as a component of a practical tool to facilitate design and optimization of coaxial shear atomizers.

Ibrahim, Essam A.↗