Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
DOE OSTI2020

Full circle mechanical dynamic characterization including experimental modal analysis and finite element analysis

During operation, it was observed that a specific mechanical system experienced undesirable vibration and it became necessary to understand and mitigate this phenomenon. This document investigates the tools, methodology, and results of the dynamic characterization of the system. The characterization makes use of the experimental modal analysis (EMA) methods of single input multiple output (SIMO) and single input single output (SISO). The validity of the theory of reciprocity is confirmed to minimize measurement error, cost, and time of repeat testing. Finite element analysis (FEA) is used in choosing transducer and modal impact locations to adequately characterize the system. Single degree of freedom (SDOF) and multiple degree of freedom (MDOF) curve fitting is used to fully characterize the system’s mode shapes and natural frequencies. The EMA characterization results are used to modify and validate the FEA model so that FEA can be used to model potential structural modifications to the system to mitigate the undesirable vibration. Structural modifications are chosen, implemented, and their effectiveness is quantified using EMA. Finally, a qualitative evaluation of the methodology of FEA validation by EMA and tuning of the model to match the experimental results is discussed.

42 ENGINEERING↗
DOE OSTI2023

Computing the Central Charge of the 3D Ising CFT Using Quantum Finite Elements

The 3D Ising conformal field theory (CFT) describes different physical systems, such as uniaxial magnets or fluids, at their critical points. In absence of an analytical solution for the 3D Ising model, the scaling dimensions and operator product expansion (OPE) coefficients characterizing this CFT must be determined numerically. The currently most-cited values for these quantities have been obtained from the conformal bootstrap, while lattice calculations have so far only produced reliable results for the scaling dimensions involved in calculating the critical exponents. Using Quantum Finite Elements to investigate critical \(\phi^4\)-theory on \(\mathbb{R}\times\mathbb{S}^2\), we show in this work that it is possible to extract scaling dimensions and OPE coefficients of the 3D Ising CFT by fitting the lattice four-point function with expectations from the operator product expansion for the radially quantized CFT and extrapolating to the continuum limit. This way, we have for the first time been able to use Monte Carlo simulations to compute the central charge of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Time reversal in a phononic crystal using finite-element modeling and an equivalent circuit model

A phononic crystal acts as a dispersive medium with a phase speed that is lower than the bulk wave speed at frequencies below the resonance of a single resonator. Time reversal is used to compensate for the phase shifts caused by individual resonators as the waves enter the medium and enable focusing of acoustic waves among the crystal. An equivalent circuit, which can predict the dispersion and attenuation of the crystal model, is shown and compared to a full-wave finite-element simulation in frequency and time. The phase shift due to a single resonator is also depicted.

Kingsley, Adam D. (ORCID:0000000226611788)↗
DOE OSTI2024

Finite Element Modeling of the Phase Change in Thermally-Grown SiO 2 in SiC Systems for Gas Turbines

The operating lifetimes of SiC-based components in combustion environments are directly linked to the adhesion of the protective environmental barrier coating (EBC) layer. One of the major known failure modes for EBCs is the formation of a thick SiO 2 thermally grown oxide (TGO), which decreases coating adhesion and encourages eventual coating spallation. The effect of the TGO thickness under Yb 2 Si 2 O 7 EBCs on silicon carbide was investigated using finite element models (FEMs) with various interfacial architectures and SiO 2 TGO thicknesses. Further, the FEMs incorporated a user-defined material to simulate the volume contraction of the TGO during the silica phase transformation from β-cristobalite to α-cristobalite upon cooling from the stress-free state at 1350°C to room temperature. Systems with and without a silicon bond coating intermediary layer were assessed. It was shown that the TGO phase transformation stress (1.6–1.7 GPa) dominated the increase in stress in the TGO and EBC layers. Furthermore, it was found that stress increase in the TGO was independent of TGO thickness and interface geometry. These results indicate that stabilization of the TGO to mitigate the phase transformation could dramatically improve the performance of SiC-base components with EBCs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Deriving Transmissibility Functions from Finite Elements for Specifications

This work explores deriving transmissibility functions for a missile from a measured location at the base of the fairing to a desired location within the payload. A pressure on the outside of the fairing and the rocket motor’s excitation creates an acceleration at a measured location and a desired location. Typically, the desired location is not measured. In fact, it is typical that the payload may change, but measured acceleration at the base of the fairing is generally similar to previous test flights. Given this knowledge, it is desired to use a finite-element model to create a transmissibility function which relates acceleration from the previous test flight’s measured location at the base of the fairing to acceleration at a location in the new payload. Four methods are explored for deriving this transmissibility, with the goal of finding an appropriate transmissibility when both the pressure and rocket motor excitation are equally present. In conclusion, these methods are assessed using transient results from a simple example problem, and it is found that one of the methods gives good agreement with the transient results for the full range of loads considered.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2023

Extend an innovative HPC-Compatible Multiple Temporal-spatial Resolution Concurrent Finite Element Modeling Approach to Guide Laser Powder Bed Fusion Additive

Laser power bed fusing (PBF) additive manufacturing is a key enabling technology to manufacture highly complex and integrated automotive structures. However, the geometric complexity of PBF-AM technique also leads to highly non-uniform heating and cooling rate in the manufactured part, which may cause flaw formation and produce excessive and nonuniform residual stresses, which increase quality uncertainties and manufacture issues, leading to increases in cost and energy consumption in the form of rejected parts. In this research project, we developed an innovative Multi-Spatial-Temporal-Resolution Finite Element (MUST-FE) method and completed the corresponding high performance computation (HPC) platform-based in-house code, which enables high accuracy prediction of temperature and residual stress fields for component-scale PBF-AM manufacture in efficient computation time. The MUST-FE model is calibrated and validated with a “2D pad” AlSi10Mg experiments by matching the melt pool shape and dimension, and with a “XY-cross” AlSi10Mg experiment by matching the thermal distortion and residual stress. The innovative multi-resolution and concurrent modeling approach adopted in this code ensures accuracy and computational efficiency, which will enable energy-efficient and high-yield, low-cost manufacturing of optimized, qualifiable automotive structures and contribute towards reaching technical targets outlined in AMO’s Program Plan to develop additive manufacturing systems that deliver consistently reliable parts with predictable properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Constraint energy minimizing generalized multiscale finite element method for multi-continuum Richards equations

In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced order methods are demanded. Combining these techniques with nonlinearity, we will consider in this paper a dual-continuum model which is generalized as a multi-continuum model for a coupled system of nonlinear Richards equations as unsaturated flows, in complex heterogeneous fractured porous media; and we will solve it by a novel multiscale approach utilizing the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). In particular, such a nonlinear system will be discretized in time and then linearized by Picard iteration (whose global convergence is proved theoretically). Subsequently, we tackle the resulting linearized equations by the CEM-GMsFEM and obtain proper offline multiscale basis functions to span the multiscale space (which contains the pressure solution). More specifically, we first introduce two new sources of samples, and the GMsFEM is used over each coarse block to build local auxiliary multiscale basis functions via solving local spectral problems, that are crucial for detecting high-contrast channels. Second, per oversampled coarse region, local multiscale basis functions are created through the CEM as constrainedly minimizing an energy functional. Various numerical tests for our approach reveal that the error converges with the coarse-grid size and that only few oversampling layers as well as basis functions are needed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Coordinate transformation and construction of finite element mesh in a diverted tokamak geometry

A coordinate transformation technique between straight magnetic field line coordinate system (Ψ, θ) and Cartesian coordinate system (R, Z) is presented employing a Solov'ev solution of the Grad-Shafranov equation. Employing the equilibrium solution, the poloidal magnetic flux Ψ(R, Z) of a diverted tokamak, magnetic field line equation is solved computationally to find curves of constant poloidal angle θ, which provides us with explicit relations R = R(Ψ, θ) and Z = Z(Ψ, θ). Correspondingly, conversion from one coordinate to the other along particle trajectories in the vicinity of separatrix is demonstrated. Based on the magnetic structure, a finite element mesh is generated in a diverted tokamak geometry to solve Poisson's equation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Prediction of Substrate Distortion in Powder DED Using Finite Element Model

Abstract Leveraging the capabilities of both additive and subtractive manufacturing to repair tooling and other high value items can result in significant cost and time saving opportunities. However, when dimensional accuracy of the repaired component must be held within a specified tolerance, distortion of the component induced during repair can lead to an out of specification final component. This work presents simulation of a novel method to mitigate part distortion using sacrificial support structures. A blown powder directed energy deposition system was utilized for deposition of sacrificial support structures and clad surface depositions. Sacrificial support structures were deposited on the underside of the substrate with varying geometries prior to clad surface depositions. A finite element model was developed and verified using experimental results. The simulation and experimental control results were found to match within 6%. It has been shown that simulation can be used to predict effects of sacrificial support structures on substrate distortion in a blown powder directed energy deposition process.

Burl, Alan G.↗
DOE OSTI2020

Understanding the deformation mechanisms in Ni-based superalloys with using crystal plasticity finite element method

Nl-based superalloy is considered as a good candidate due to its excellent resistance to elevated temperature deformation for long term period application. Understanding the deformation and failure mechanisms of Ni-Based superalloys is very helpful for providing design guidelines for processing Ni-based superalloys. Experimental characterization indicates that the deformation mechanisms of Ni based superalloy is strongly microstructure dependent. Besides, damage transform from the void nucleation to the macro cracks by voids growth leading to the failure of the Ni-based superalloys are also showing strong microstructure sensitivity. Therefore, this work focuses on the prediction and comprehension of the deformation and void growth behavior in Ni based superalloy at different working conditions via crystal plasticity finite element modeling and simulation. Physically based crystal plasticity frameworks were developed for newly Ni- based superalloy Haynes 282. It was found that dislocation shearing through the precipitates were acting as the main contributor to the strength of Haynes 282 at room temperature and 815°C. Our analysis of the creeping behavior of Haynes 282 exhibited that resistance of general climb replaced by the resistance induced by the deposited climb dislocation density. In addition, in the study of void growth behavior, our simulation results demonstrated that as the main loading axis perpendicular to the grain boundary (GB), voids grow more slowly on tilt GBs in bicrystals than those in single and bicrystals with twist GBs. And tilt GBs would promote the void grow into irregular shape.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A Novel Partitioned Approach for Reduced Order Model—Finite Element Model (ROM-FEM) and ROM-ROM Coupling

Partitioned methods allow one to build a simulation capability for coupled problems by reusing existing single-component codes. In so doing, partitioned methods can shorten code development and validation times for multiphysics and multiscale applications. In this work, we consider a scenario in which one or more of the “codes” being coupled are projection-based reduced order models (ROMs), introduced to lower the computational cost associated with a particular component. We simulate this scenario by considering a model interface problem that is discretized independently on two non-overlapping subdomains. Here we then formulate a partitioned scheme for this problem that allows the coupling between a ROM “code” for one of the subdomains with a finite element model (FEM) or ROM “code” for the other subdomain. The ROM “codes” are constructed by performing proper orthogonal decomposition (POD) on a snapshot ensemble to obtain a low-dimensional reduced order basis, followed by a Galerkin projection onto this basis. The ROM and/or FEM “codes” on each subdomain are then coupled using a Lagrange multiplier representing the interface flux. To partition the resulting monolithic problem, we first eliminate the flux through a dual Schur complement. Application of an explicit time integration scheme to the transformed monolithic problem decouples the subdomain equations, allowing their independent solution for the next time step. We show numerical results that demonstrate the proposed method’s efficacy in achieving both ROM-FEM and ROM-ROM coupling.

de Castro, Amy↗
DOE OSTI2020

Evaluation of the Performance of the SAVY Container Subjected to Various Impact Events using Finite Element Methods

The SAVY container is a Safety Class Design Feature utilized both inside engineered features and outside engineered features (e.g. gloveboxes) at TA-55 Plutonium Facility Building 4. These containers were certified as low damage ratio containers by demonstrating their ability to retain material after exposure to impacts and elevated temperatures resulting from a design basis seismic event. Original testing of the SAVY did not evaluate the effects of various facility specific insults to the SAVY while within a DOT Type A Container, caused either by design basis accident or an operational mishap. These insults include: 1. Impact of a structural I-beam (drum upright & drum on side); 2. Impact of a section of Positive Pressure Chilled Water Pipe (PPCWS); 3. Forklift fork tine impact; 4. Fuel truck impact. Abaqus Explicit Dynamics Finite Element Analysis was performed on a nested SAVY/DOT Type 7A container configuration to assess the structural integrity of the SAVY under the accident conditions listed.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2023

Determining Stress in Metallic Conducting Layers of Microelectronics Devices Using High Resolution Electron Backscatter Diffraction and Finite Element Analysis

Delayed failure due to stress voiding is a concern with some aging microelectronics, as these voids can grow large enough to cause an open circuit. Local measurements of stress in the metallic layers are crucial to understanding and predicting this failure, but such measurements are complicated by the fact that exposing the aluminum conducting lines will relieve most of their stress. In this study, we instead mechanically thin the device substrate and measure distortions on the thinned surface using high resolution electron backscatter diffraction (HREBSD). These measurements are then related to the stresses in the metallic layers through elastic simulations. This study found that in legacy components that had no obvious voids, the stresses were comparable to the theoretical stresses at the time of manufacture (≈300 MPa). Distortion fields in the substrate were also determined around known voids, which may be directly compared to stress voiding models. In conclusion, the technique presented here for stress determination, HREBSD coupled with finite element analysis to infer subsurface stresses, is a valuable tool for assessing failure in layered microelectronics devices.

HREBSD↗
DOE OSTI2023

Modeling the First Hydrogen Direct Reduction Pilot Reactor for Ironmaking in the USA Using Finite Element Analysis and Its Validation Using Pilot Plant Trial Data

Direct reduction of hematite pellets with hydrogen (H 2 ) was used to produce directly reduced iron (DRI) in a pilot scale reactor at a pellet feed rate of 21.4 kg/h. At a steady state, operational parameters of the pilot plant (gas recycling rate and inlet temperature) along with key reactor output parameters, the pellet metallization, and the internal temperature profile of the reactor were reported for two scenarios with high recycle and low recycle rate of H 2 . Scenario 1, with a high recycle rate of 400 L/min H 2 along with external heating of 870 °C, gave an average metallization of 91.8%, while Scenario 2, with low recycle rate of 100 L/min H 2 and external heating of 850 °C gave a metallization of 67.8% due to the higher moles of H 2 available for reduction and the external energy required for the endothermic reduction reaction in Scenario 1 as compared with Scenario 2. Finite element analysis was used to build a model of the shaft reactor, which was validated against the metallization and internal temperature profile data. The average metallization values predicted by the model were very close to the metallization values obtained from the pilot plant samples, with 90.9% average metallization for Scenario 1 and 65.6% average metallization for Scenario 2. The internal temperature profiles in the lower region of the reactor obtained from the model were very close to these pilot plant data, with a maximum difference of 52.7 °C and 67.6 °C for Scenarios 1 and 2, respectively. The pilot plant reactor model was used extensively in the commissioning of the pilot plant and to predict the startup outcomes for a given set of operating parameters.

08 HYDROGEN↗
DOE OSTI2023

Evaluation of Surface Crack Formation in Photovoltaic Backsheets Using Fragmentation and Finite Element Simulations: Preprint

Backsheet cracking is among the most commonly observed degradation modes of photovoltaic (PV) modules in the field. Cracks can reduce the ability of backsheets to fulfil their functions, for example, protection of the modules from the environment or electrical insulation. This work presents an evaluation of the degradation and cracking propensity of two backsheets during accelerated IEC TS 62788-7-2 A3 (International Electrotechnical Commission) aging with the fragmentation test: a co-extruded polyamide backsheet (AAA), and a laminated multilayer backsheet with a polyethylene terephthalate core and outer layer and ethylene vinyl acetate inner layer (PPE). Results show the surface embrittlement of the AAA outer layer during exposure. A longer time of exposure causes cracks to form at lower strains during stretching and creates deeper cracks. Accordingly, Young's modulus of the outer layer increases, as measured by cross-sectional nanoindentation. PPE exhibits cracking after exposure as well. While no similar increase of modulus or crack depth can be observed, the outer layer of PPE exhibits more obvious signs of erosion during exposure, including progressive morphological changes and thickness losses. A finite element model to simulate surface crack formation is devised, based on initially zero-thickness decohesion elements. Decohesion criteria define a critical stress, at which these elements grow, and cracks can begin to form. These criteria are obtained via parameter optimization by comparison between simulation and experiment. The model is used to interpret crack formation in both backsheets.

backsheet↗
DOE OSTI2023

Using the Stix finite element RF code to investigate operation optimization of the ICRF antenna on Alcator C-Mod

Abstract As the Ion Cyclotron Radio Frequency range (ICRF) heating becomes more favorable in fusion devices, the urgency of predicting and mitigating impurity generation that arises from it becomes more pressing. In the ICRF regime, rectified Radio Frequency (RF) sheaths are known to form at antenna and material edges that influence negative effects like sputtering and a decrease in heating efficiency. Methods to mitigate the formation of these RF sheaths through RF image currents cancellation have been experimentally studied. A power-phasing scan done on Alcator C-Mod in which the amount of power on the two inner straps ( P in ) versus the total 4 straps ( P tot ) was varied showed a minimization of enhanced potentials between P in / P tot ∼ 0.7–0.9 while impurities were minimized for P in / P tot ∼ 0.5–0.8. New capabilities in the realm of representing the RF sheath numerically now allow for these experiments to be simulated. Given the size of the sheath relative to the scale of the device, it can be approximated as a Boundary Condition (BC). A new parallelized cold-plasma wave equation solver called Stix implements a non-linear sheath impedance model BC formulated by Myra et al (2015 Phys. Plasmas 22 062507) through the method of finite elements using the MFEM library [ http://mfem.org ]. It is seen that Stix shows qualitative agreement with the measured C-Mod enhanced potentials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗