Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Elemental analysis”

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 325 records · Page 18

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

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

97 MATHEMATICS AND COMPUTING↗

Finite element analysis of unsteady incompressible flow around an oscillating obstacle of arbitrary shape.

An analytical procedure based on Navier-Stokes equations was developed for representing unsteady flow patterns around oscillating obstacles. A variational formulation of the Helmholtz vorticity equation was discretized in finite element form and integrated numerically. At each step of the numerical integration the velocity field around the obstacle was determined from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints at each time step of the numerical integration. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.

Bratanow, T.↗

Pit Stability Predictions of Additively Manufactured SS316 Surfaces Using Finite Element Analysis

Stainless steels are susceptible to localized forms of corrosion attack, such as pitting. The size and lifetime of a nucleated pit can vary, depending on a critical potential or current density criterion, which determines if the pit repassivates or continues growing. This work uses finite element method (FEM) modeling to compare the critical pit radii predicted by thermodynamic and kinetic repassivation criteria. Experimental electrochemical boundary conditions are used to capture the active pit kinetics. Geometric and environmental parameters, such as the pit shape and size (analogous to additively manufactured lack-of-fusion pores), solution concentration, and water layer thickness were considered to assess their impact on the pit repassivation criterion. The critical pit radius (the transition point from stable growth to repassivation) predicted for a hemispherical pit was larger when using the repassivation potential (E rp ) criteria, as opposed to the current density criteria (pit stability product). Including both the pit stability product and E rp into its calculations, the analytical maximum pit model predicted a critical radius two times more conservative than the FEA approach, under the conditions studied herein. The complex pits representing lack-of-fusion pores were shown to have minimal impact on the critical radius in atmospheric conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Symbolic derivation of material property matrices in finite element analysis

The principles and operation of MMAX, a symbolic-computation program which automates the process of generating property matrices for structural materials, are briefly described and illustrated with sample analyses of a rubberlike material and an elastoplastic material. MMAX is written in LISP under the symbolic finite-element generator FINGER and the general symbolic manipulator MACSYMA; it first derives the formulas required by mathematical manipulation, and then translates the formulas into FORTRAN code, adapted to the particular type of machine to be used for the numerical calculations. This approach is shown to combine efficiently the advantages of symbolic and numerical computation for engineering applications.

Tan, H. Q.↗

A REDUCED ORDER MODELING APPROACH TO PROBABILISTIC CREEP-DAMAGE PREDICTIONS IN FINITE ELEMENT ANALYSIS

This paper introduces a computationally efficient Reduced Order Modeling (ROM) approach for the probabilistic prediction of creep-damage failure. Component-level probabilistic simulations are needed to assess the reliability and safety of high-temperature components. Full-scale probabilistic creep-damage modeling in finite element (FE) approach is computationally expensive requiring many hundreds of simulations to replicate the uncertainty of component failure. To that end, ROM is proposed to minimize the elevated computational cost while controlling the loss of accuracy. It is proposed that full-scale probabilistic simulations can be completed in 1D at a reduced cost, the extremum conditions extracted, and those conditions applied for lower-cost 2D/3D probabilistic simulations of components that capture the mean and uncertainty of failure. The probabilistic Sine-hyperbolic (Sinh) model is selected which in previous work was calibrated to alloy 304 stainless steel. The Sinh model includes probability density functions (pdfs) for test condition (stress and temperature), initial damage (i.e. microstructure), and material properties uncertainty. The Sinh model is programmed into ANSYS finite element software using the USERCREEP.F material subroutine. First, the Sinh model and FE code are subject to verification and validation to ensure the accuracy of the simulations. Numerous Monte Carlo simulations are executed in a 1D model to generate probabilistic creep deformation, damage, and rupture data. This data is analyzed and the probabilistic parameters corresponding to extreme creep response are extracted. The ROM concept is applied where only the extreme conditions are applied in the 2D probabilistic prediction of a component. The probabilistic predictions between the 1D and 2D geometry is compared to assess ROM for creep. The accuracy of the probabilistic prediction employing the ROM approach will potentially reduce the time and cost of simulating complex engineering systems. Future studies will introduce multi-stage Sinh, stochasticity, and spatial uncertainty for improved prediction.

36 MATERIALS SCIENCE↗

Vibration Propagation of Gear Dynamics in a Gear-Bearing-Housing System Using Mathematical Modeling and Finite Element Analysis

Vibration and noise caused by gear dynamics at the meshing teeth propagate through power transmission components to the surrounding environment. This study is devoted to developing computational tools to investigate the vibro-acoustic propagation of gear dynamics through a gearbox using different bearings. Detailed finite element/contact mechanics and boundary element models of the gear/bearing/housing system are established to compute the system vibration and noise propagation. Both vibration and acoustic models are validated by experiments including the vibration modal testing and sound field measurements. The effectiveness of each bearing type to disrupt vibration propagation is speed-dependent. Housing plays an important role in noise radiation .It, however, has limited effects on gear dynamics. Bearings are critical components in drivetrains. Accurate modeling of rolling element bearings is essential to assess vibration and noise of drivetrain systems. This study also seeks to fully describe the vibro-acoustic propagation of gear dynamics through a power-transmission system using rolling element and fluid film wave bearings. Fluid film wave bearings, which have higher damping than rolling element bearings, could offer an energy dissipation mechanism that reduces the gearbox noise. The effectiveness of each bearing type to disrupt vibration propagation in explored using multi-body computational models. These models include gears, shafts, rolling element and fluid film wave bearings, and the housing. Radiated noise is mapped from the gearbox surface to surrounding environment. The effectiveness of rolling element and fluid film wave bearings in breaking the vibro-acoustic propagation path from the gear to the housing is investigated.

Parker, Robert G.↗

Real-time elemental analysis of liquids for process monitoring using laser-induced breakdown spectroscopy with a liquid wheel sampling approach

This article presents an engineered sampling system that used a rotating wheel to form a thin liquid layer, permitting the use of laser-induced breakdown spectroscopy (LIBS) for in situ, real-time elemental impurity quantification during liquid processing. The sampling approach was demonstrated on eight elements from across the periodic table (Na, Al, K, Ca, Ti, Sr, Mo, and Yb). Univariate and multivariate calibrations were presented for each element. The average value for percent root mean square errors of cross-validation for the multivariate models was 3.64%, highlighting the method's strong prediction accuracy. Additionally, the limits of detection for each analyte were estimated from their univariate models: Na = 0.0532, Al = 18.5, K = 0.105, Ca = 0.273, Ti = 67.7, Sr = 0.640, Mo = 22.4, and Yb = 22.9 μg mL –1 . Finally, a test in which multivariate models were used to monitor a liquid system for 80 min was performed to investigate the real-time monitoring capabilities of this liquid LIBS sampling approach. Rigorous measurements were performed to effectively predict the absence and concentrations of multiple analytes as they were spiked and diluted. This demonstration showed the feasibility of using LIBS for real-time liquid quantification models with estimated precision ≤ 8.1%. Finally, the limitations of this approach and potential future improvements are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interplanetary dust - Trace element analysis of individual particles by neutron activation

Although micrometeorites of cometary origin are thought to be the dominant component of interplanetary dust, it has never been possible to positively identify such micrometer-sized particles. Two such particles have been identified as definitely micrometeorites since their abundances of volatile and nonvolatile trace elements closely match those of primitive solar system material.

Ganapathy, R.↗

Nonlinear finite-element analysis of laminated composite shells

Two aspects of the geometrically nonlinear analysis of laminated composite shells are considered in this paper. The first pertains to the simplifications and computational advantages resulting from the use of group-theoretic methods in conjunction with computerized symbolic (algebraic) manipulation in the element development. The second aspect pertains to the combined effects of shear deformation and nonlinearity on the accuracy and convergence of two-dimensional, shear-flexible stiffness finite element models.

Noor, A. K.↗

Recent developments in finite element analysis for transonic airfoils

The prediction of aerodynamic forces in the transonic regime generally requires a flow field calculation to solve the governing non-linear mixed elliptic-hyperbolic partial differential equations. Finite difference techniques were developed to the point that design and analysis application are routine, and continual improvements are being made by various research groups. The principal limitation in extending finite difference methods to complex three-dimensional geometries is the construction of a suitable mesh system. Finite element techniques are attractive since their application to other problems have permitted irregular mesh elements to be employed. The purpose of this paper is to review the recent developments in the application of finite element methods to transonic flow problems and to report some recent results.

Hafez, M. M.↗

Evaluation of the use of a singularity element in finite element analysis of center-cracked plates

Two different methods are applied to the analyses of finite width linear elastic plates with central cracks. Both methods give displacements as a primary part of the solution. One method makes use of Fourier transforms. The second method employs a coarse mesh of triangular second-order finite elements in conjunction with a single singularity element subjected to appropriate additional constraints. The displacements obtained by these two methods are in very good agreement. The results suggest considerable potential for the use of a cracked element for related crack problems, particularly in connection with the extension to nonlinear material behavior.

Mendelson, A.↗

Finite element analysis of the effect of a non-planar solid-liquid interface on the lateral solute segregation during unidirectional solidification

The effect of solid-liquid interface shape on lateral solute segregation during steady-state unidirectional solidification of a binary mixture is calculated under the assumption of no convection in the liquid. A finite element technique is employed to compute the concentration field in the liquid and the lateral segregation in the solid with a curved boundary between the liquid and solid phases. The computational model is constructed assuming knowledge of the solid-liquid interface shape; no attempt is made to relate this shape to the thermal field. The influence of interface curvature on the lateral compositional variation is investigated over a range of system parameters including diffusivity, growth speed, distribution coefficient, and geometric factors of the system. In the limiting case of a slightly nonplanar interface, numerical results from the finite element technique are in good agreement with the analytical solutions of Coriell and Sekerka obtained by using linear theory. For the general case of highly non-planar interface shapes, the linear theory fails and the concentration field in the liquid as well as the lateral solute segregation in the solid can be calculated by using the finite element method.

Carlson, F. M.↗

Three-Dimensional Effects in the Plate Element Analysis of Stitched Textile Composites

Three-dimensional effects related to the analysis of stitched textile composites are discussed. The method of calculation is based on the virtual crack closure technique (VCCT), and models that model the upper and lower surface of the delamination or debond with two-dimensional (2D) plate elements rather than three-dimensional (3D) solid elements. The major advantages of the plate element modeling technique are a smaller model size and simpler geometric modeling. Details of the modeling of the laminated plate and the stitching are discussed.

Glaessgen, E. H.↗

Multi-Fidelity Scheme for Accelerating 4D Finite Element Analysis for Mircoreactors

Next-generation microreactors are currently being designed to be operated terrestrial and extraterrestrial for remote surface power production. These systems will provide an alternative source of carbon-free energy that is versatile and can be utilized for various applications. These applications of microreactors have prompted the usage of high-fidelity unstructured finite element (FE) based approaches to provide time-dependent solutions for multiple physics fields. Computing these high-fidelity (4-D) solutions for multiple design iterations, physics fields, and transient events requires an immense computational cost. These high-fidelity solutions are computational expensive and the cost can be reduced through the implementation of less accurate low-fidelity solutions. Unlike the current fleet of commercial nuclear reactors, these next-gen systems present challenges due to the material and physical limitations required. Due to these constraints, the brute force technique of parameterizing important system characteristics determines whether a design meets the project objective. A system such as a nuclear reactor could have thousands of design parameters that affect system performance. In order to analyze the entire parameter space of such a complex system would require millions of CPU hours and countless design iterations. This costly approach is not practical due to regulatory and budget limitations. In this proposal, an approach that utilizes hybrid high-fidelity and low-fidelity FE models to reduce the computational cost of evaluating these applications in 4-D will be presented. The accuracy of the high-fidelity model and the computational efficiency of the low-fidelity model are taken advantage of to produce a solution that closely resembles the full order high-fidelity solution. By utilizing an unconverged coarse FE mesh, operation limits such as temperature, structural loading, etc., can be evaluated in an accelerated fashion and then can be spatially interpolated onto a finer mesh. The resulting error arising from the coarse mesh can be mitigated with a discrepancy function that actively quantifies and corrects the error in the coarse solution. This discrepancy function can be periodically updated with high-fidelity calculations across the temporal domain thus requiring less iterations on the fine mesh. As a result, the computational cost can be reduced for the evaluation and design iteration of microreactors. The proposal for this research contains three sections: Section 2 provides a literature review, Section 3 outlines the methodology for the proposed multi-fidelity scheme, and Section 4 displays preliminary results of the proposed research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Probabilistic Fatigue Damage Prognosis Using a Surrogate Model Trained Via 3D Finite Element Analysis

Utilizing inverse uncertainty quantification techniques, structural health monitoring can be integrated with damage progression models to form probabilistic predictions of a structure's remaining useful life. However, damage evolution in realistic structures is physically complex. Accurately representing this behavior requires high-fidelity models which are typically computationally prohibitive. In the present work, a high-fidelity finite element model is represented by a surrogate model, reducing computation times. The new approach is used with damage diagnosis data to form a probabilistic prediction of remaining useful life for a test specimen under mixed-mode conditions.

Leser, Patrick E.↗

Finite Element Analysis of the Random Response Suppression of Composite Panels at Elevated Temperatures using Shape Memory Alloy Fibers

A feasibility study on the use of shape memory alloys (SMA) for suppression of the random response of composite panels due to acoustic loads at elevated temperatures is presented. The constitutive relations for a composite lamina with embedded SMA fibers are developed. The finite element governing equations and the solution procedures for a composite plate subjected to combined acoustic and thermal loads are presented. Solutions include: 1) Critical buckling temperature; 2) Flat panel random response; 3) Thermal postbuckling deflection; 4) Random response of a thermally buckled panel. The preliminary results demonstrate that the SMA fibers can completely eliminate the thermal postbuckling deflection and significantly reduce the random response at elevated temperatures.

Turner, Travis L.↗

Hybrid-finite-element analysis of some nonlinear and 3-dimensional problems of engineering fracture mechanics

In this paper, efficient numerical methods for the analysis of crack-closure effects on fatigue-crack-growth-rates, in plane stress situations, and for the solution of stress-intensity factors for arbitrary shaped surface flaws in pressure vessels, are presented. For the former problem, an elastic-plastic finite element procedure valid for the case of finite deformation gradients is developed and crack growth is simulated by the translation of near-crack-tip elements with embedded plastic singularities. For the latter problem, an embedded-elastic-singularity hybrid finite element method, which leads to a direct evaluation of K-factors, is employed.

Atluri, S. N.↗