Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Element Model”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A predictive strain-gradient model with no undetermined constants or length scales

A general meso–scale (GM) crystal plasticity (CP) model was developed that accounts for lower-order (strain hardening) and higher-order (internal stress) effects of geometrically necessary dislocations (GNDs). It is predictive: no arbitrary parameters or length scales were invoked and no ad hoc numerical techniques were employed. It uses general stress field equations for GND content and a novel harmonization technique to enforce consistency of elastic long-range singular defect fields with applied elastic-plastic fields. The model facilitates implementation in commercial finite element programs without requiring special elements, special boundary conditions, or access to element shape functions. GM simulations confirmed, with improved accuracy, previously published predictions of the Hall-Petch effect, Bauschinger effect, and anelasticity. Previously unpredicted phenomena were also predicted: anelasticity and hysteresis for single Ta crystals and strain-hardening stagnation. The internal stresses (higher-order effect) dominate at large length scales, while at small length scales, the GND density hardening (lower-order effect) dominates. As a result, GM predicts that strain heterogeneity and consequent GND internal stresses are important factors in anelasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

BISON validation of FeCrAl cladding mechanical failure during simulated reactivity-initiated accident conditions*

Here, a reactivity-initiated accident (RIA) is a postulated design basis accident in light water reactors (LWRs) in which a rapid reactivity insertion induces a fission rate increase and a fuel pellet temperature rise. During RIA, the fuel pellet thermally expands and may cause pellet-cladding mechanical interaction (PCMI). Separate effects PCMI tests were performed on C26M FeCrAl cladding tube samples, introducing biaxial stress via a well understood modified burst test (MBT) system. A high-speed camera in the MBT system captured the projections, covering a 360° view of the cladding deformation and enabling a digital image correlation (DIC) method to quantify the surface strains with high fidelity. Representative hot zero-power RIA mechanical loading conditions were applied to the sample, and the test duration ranged from 20 to 500 ms at an average temperature of 573 K. BISON finite element–based fuel performance code modeling was performed against the high-fidelity DIC data produced from the MBTs. Validation calculations were conducted with 2D models and systematically compared with test data of the burst time, burst pressure, burst hoop strains, and hoop strain rates. Based on the behaviors from the separate effects test, BISON calculations satisfactorily predicted the cladding deformation behaviors. Sensitivity analysis was conducted to identify highly influential mechanical properties responsible for the cladding failure behavior during the MBT experiments. The results highlight the significance of the cladding's mechanical strength in governing cladding strain, followed by the significance of the mechanical interactions between the pellet and the cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Material-dependent photon ionizing radiation effects in Si and GaAs PIN diodes: A numerical investigation

We present a finite-element drift-diffusion-Poisson model in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework to compare the radiation response of silicon (Si) and gallium arsenide (GaAs) PIN diodes under high-energy photon irradiation. The model solves coupled carrier continuity and Poisson’s equations with Shockley-Read-Hall recombination, and is verified against standard analytical J-V behavior. Using a simplified 1D geometry with ideal Ohmic contacts, we quantify device response under forward and reverse bias with a 100 MeV photon flux. Under forward bias, Si exhibits markedly greater radiation sensitivity than GaAs, including larger increases in current density, stronger local field and carrier-product perturbations, and higher recombination. Under reverse bias, GaAs shows larger radiation-induced photocurrent and broader current-density peaks near junctions, indicating an advantage for photodetection. Integrated steady-state recombination is consistently higher in Si across voltages. Under periodic photon pulses, GaAs produces higher-amplitude photoresponse and settles more rapidly than Si. These results highlight material-dependent trade-offs for radiation-tolerant, high-speed optoelectronics and provide guidance for selecting PIN architectures in aerospace, nuclear, and high-energy physics environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Isomerization kinetics of azobenzene crosslinked to a liquid crystal polymer network

Azobenzene is an important chromophore having numerous applications that are a direct implication of its ability to change shape upon irradiation with a suitable wavelength. Here, in this study, photoisomerization kinetics of a diacrylate azobenzene molecule (A3MA) cross-linked to a liquid crystal polymer network (LCN) is investigated. Isomerization experiments were conducted on planar LCN thin films with varying cross-linking densities for a fixed azobenzene concentration. The experiments were conducted at approximately 26 °C, which is significantly below each film’s glass transition temperature. The forward trans-cis and reverse cis-trans isomerization slows down in crosslinked azobenzene liquid crystal polymer networks (ALCNs) as compared to a solution and an uncrosslinked polymer mixture. The observed slowdown is due to the reduction in the free volume and increased steric hindrances with an increase in crosslinking. A single-term first-order kinetics rate equation is known to describe the isomerization kinetics in solutions or uncrosslinked mixtures while the same form fails to represent the reaction kinetics in moderate to highly cross-linked ALCNs. Our results for the first time indicate that a two-term first-order reaction rate is an accurate way to describe the forward trans-cis and reverse thermal cis-trans isomerization in crosslinked azobenzene liquid crystal polymer networks. In addition, the effect of film thickness and incident intensities on the reaction kinetics is investigated. Finally, the two-term isomerization kinetics is incorporated into a finite element photo-chemo-mechanical model to reveal the consequent slowdown in the actuation response of ALCN films. The current findings will be of particular interest to theoreticians and experimentalists alike, for it will help predict and understand accurately the responses of azobenzene-based coatings, and actuators and may have potential applications in underwater soft robotics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Thermomechanical response of thickly tamped targets and diamond anvil cells under pulsed hard x-ray irradiation

In the laboratory study of extreme conditions of temperature and density, the exposure of matter to high intensity radiation sources has been of central importance. Here, we interrogate the performance of multi-layered targets in experiments involving high intensity, hard x-ray irradiation, motivated by the advent of extremely high brightness hard x-ray sources, such as free electron lasers and 4th-generation synchrotron facilities. Intense hard x-ray beams can deliver significant energy in targets having thick x-ray transparent layers (tampers) around samples of interest for the study of novel states of matter and materials’ dynamics. Heated-state lifetimes in such targets can approach the microsecond level, regardless of radiation pulse duration, enabling the exploration of conditions of local thermal and thermodynamic equilibrium at extreme temperature in solid density matter. The thermal and mechanical responses of such thick layered targets following x-ray heating, including hydrodynamic relaxation and heat flow on picosecond to millisecond timescales, are modeled using radiation hydrocode simulation, finite element analysis, and thermodynamic calculations. Assessing the potential for target survival over one or more exposures and resistance to damage arising from heating and resulting mechanical stresses, this study doubles as an investigation into the performance of diamond anvil high pressure cells under high x-ray fluences. Long used in conjunction with synchrotron x-ray radiation and high power optical lasers, the strong confinement afforded by such cells suggests novel applications at emerging high intensity x-ray facilities and new routes to studying thermodynamic equilibrium states of warm, very dense matter.

Meza-Galvez, J. (ORCID:000000032699849X)↗
DOE OSTI2024

Mechanically activated and deactivated ion transport across nanopores with heterogeneous surface charge distributions

To mimic the intricate and adaptive functionalities of biological ion channels, electrohydrodynamic ion transport has been studied extensively, albeit mostly, across uniformly charged nanochannels. Here, we analyze the ion transport under coupled electric field and pressure across heterogeneously charged nanopores with oppositely charged sections on their lateral surface. We only consider such pores with symmetric hourglass-like and cylindrical shapes to focus on the effects of the non-uniform surface charge distribution. Finite-element simulations of a continuum model demonstrate that a pressure applied in either direction of the pore-axis equally suppresses or amplifies the ionic conductance, depending on the electric field polarity, by distorting the quasi-static distribution of ions in the pore. The resulting anomalous mechanical deactivation and activation of ionic current under opposite voltage biases exhibit the functional modularity of our setup, while their intensities are highly tunable, substantially greater than those of analogous behaviors in other nanochannels, and fundamentally correlated to ionic current rectification (ICR) in our pores. A detailed study of ICR subsequently reveals its counterintuitive non-monotonous variations, in the pores, with the magnitude of applied voltage and the pore length, that can help optimize their diode-like behavior. We further illustrate that while the hourglass-shaped nanopores yield the more efficient mechanical suppressors of ion transport, their cylindrical analogs are the superior rectifiers and mechanical amplifiers of ion conduction. Therefore, this article provides a blueprint for the strategic design of nanofluidic circuits to attain a robust, modular, and tunable control of ion transport under external electrical and mechanical stimuli.

Physics↗
DOE OSTI2023

A Reconfigurable Planar Capacitor Sensor for Thru the Jacket Dielectric Constant Measurement of Cable Insulation

Planar capacitor sensors have been proposed to measure the aging-related dielectric constant variation of the jacket or insulation materials of cables. Although capacitor sensors have been used to monitor the aging of cable insulation by removing the jacket, through the jacket measurement of cable insulation dielectric constant variation is difficult because once the sensor is placed on the jacket it would either measure the aging of the jacket or that of the jacket and insulation together. A reconfigurable sensor along with innovative cable-specific models could be used to determine the jacket and insulation aging using the same sensor. For the proposed novel reconfigurable sensor, the electric field penetration depth is changed by activating and deactivating PIN diode switches. In one instance, the sensor measures the dielectric constant variation of the jacket, while in the next, it measures the dielectric constant variation of both the jacket and the insulation. Here, by leveraging previously developed capacitance estimation models from large-scale finite element simulations for a specific cable, these two sets of measurement data are used to evaluate the aging-related dielectric constant variation of the insulation. Experimental data from aged cable specimens are presented to demonstrate the operation of the sensor.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

NEML2 - THE NEW ENGINEERING MATERIAL MODEL LIBRARY, VERSION 2

SF-22-157 NEML2 is an extension of an earlier open source code developed at ANL (SF-xxxx) called NEML. These libraries provide tools for constructing and using models of engineering materials, with a focus on high temperature materials for nuclear reactor applications. The software include tools for setting up material models, calibrating models, storing models for later use in other simulations, and running simple "material point" simulations of the material behavior. The main difference between NEML and NEML2 is the new software has been complete restructured to run on GPUs and other accelerators in order to improve the performance of methods used to train constitutive models against data and to improve the performance of the final models in production environments, like finite element software. While over time we expect NEML2 to become more-or-less feature complete with NEML we are filing a new software disclosure because (1) the library has been complete restructured and (2) we want to host NEML2 in a separate repository.

MESSNER, MARK↗
DOE OSTI2021

Rapid simulations of hyperspectral near-field images of three-dimensional heterogeneous surfaces

The scattering-type scanning near-field optical microscope (s-SNOM) has emerged as a powerful tool for resolving nanoscale inhomogeneities in laterally heterogeneous samples. However, most analytical models used to predict the scattering near-field signals are assuming homogenous landscapes (bulk materials), resulting in inconsistencies when applied to samples with more complex configurations. In this work, we combine the point-dipole model (PDM) to the finite-element method (FEM) to account for the lateral and vertical heterogeneities while keeping the computation time manageable. Full images, spectra, or hyperspectral line profiles can be simulated by calculating the self-consistent dipole radiation demodulated at higher harmonics of the tip oscillation, mimicking real experimental procedures. Using this formalism, we clarify several important yet puzzling experimental observations in near-field images on samples with rich typography and complex material compositions, heterostructures of two-dimensional material flakes, and plasmonic antennas. The developed method serves as a basis for future investigations of nano-systems with nontrivial topography.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Preliminary design analysis workflow for Division 5 HHA-3200 requirements for graphite core components

This report presents a design analysis workflow for graphite core components and assemblies, based on the design rules of ASME Boiler Pressure and Vessel Code, Section III, Division 5, Article HHA-3000. The workflow contains three stages: developing the design of the graphite core component, modeling the component with the finite element software MOOSE, and assessing if the component passes/fails the criteria of the HHA-3000 design rules. Since the design rules use probabilistic metrics specifically established to evaluate brittle materials, we developed a python library that performs all the statistical calculations necessary for the evaluations of the HHA-3000 criteria.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Error Analysis on Numerical Integration Algorithms in a Hypoelasticity Framework

This report determines local truncation errors for common stress integration algorithms used in explicit finite element codes with hypoelastic material models. The hypoelastic integration algorithms in question utilize an operator splitting procedure in a rotation neutralized configuration, where the stress response is determined from de- coupling the total deformation into rotational and strain dependent components. This document analyzes the error in evolving the stress given a one-step time increment Δt and compares the errors associated with both the rotational and strain components of the operator splitting method. A slight modification to a traditional algorithm is proposed and studied, where the rate of deformation is appropriately rotated from the midstep configuration at t n+1/2 to the end step configuration at t n+1 before the constitutive evaluation. The proposed modification either completely eliminates the error associated with the rotation rate or is of the same order of magnitude as the original algorithm for the three test cases considered in this report. These cases consist of an unaxial stretch with a constant true strain rate with a rigid body rotation, an uniaxial stretch with a constant engineering strain rate with a rigid body rotation, and a simple shear deformation. All three cases are compared to a closed form solution, and in almost every test case the alternative algorithm yields the most accurate one-step local truncation error.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Specimen Size Artifacts Associated with a Glovebox Deployable Laser Flash Diffusivity Measurement System

To accelerate nuclear fuel qualification and deployment efforts, Oak Ridge National Laboratory has led the world in developing vehicles for accelerating burnup accumulation in nuclear fuel by testing miniature samples in the High Flux Isotope Reactor. However, a challenge remains to establish parity between post irradiation characterization methods for miniature fuel specimens in comparison to conventional fuel elements. This work aims to identify artifacts associated with the thermal analysis of miniature disk specimens in comparison to samples of standardized dimensions through a combination of experimental and computational tools. Using a newly procured Netzsch 427 laser flash diffusivity system, preliminary data on mild steel specimens were collected, and finite element techniques were used to model heat transfer using 1D conditions to assess issues associated with pulse width and data collection frequency effects. These preliminary investigations found that significant variability in measured thermal diffusivity exist within datasets. These variations are not necessarily attributed to limitations with respect to data collection rate; nor are they associated with limitations on pulse width. Instead, these artifacts may be a result of limitations in detector intensity. Future work aims to design a focusing lens to better measure the smaller surface area of miniature specimens analyzed using the laser flash analysis system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

OpenFUSIONToolkit

The Open FUSION Toolkit (OFT) is a suite of modeling tools, and their underlying finite element framework, for problems in plasma and fusion research and engineering in arbitrary 2D and 3D geometries.

Hansen, Chris [Columbia Univ., New York, NY (Unite↗
DOE OSTI2021

BISON As-run AGR-3/4 Irradiation Test Predictions

BISON, a nuclear fuel performance application built using the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element library, was used to model the Advanced Gas Reactor (AGR)-3/4 irradiation test using as-run physics and thermal hydraulics data. The AGR-3/4 test consists of the combined third and fourth planned irradiations of the AGR Fuel Development and Qualification Program. The AGR-3/4 test train consists of twelve separate and independently controlled and monitored capsules. Each capsule contains four compacts filled with both uranium oxycarbide (UCO) unaltered “driver” fuel particles and UCO designed-to-fail (DTF) fuel particles. The DTF fraction was specified to be 1×10-2. This report documents the calculations performed to predict the failure probability of tristructural isotropic (TRISO) coated driver fuel particles during the AGR-3/4 experiment on a single compact. This report will demonstrate the capabilities of BISON to model the complex AGR-3/4 irradiation test and identify further development needed to capture the fuel particle failure probability and source term on every compact for further comparison from post-irradiation examination (PIE) data. The calculations include the modeling of the AGR-3/4 irradiation that occurred from December 2011 to April 2014 in the Advanced Test Reactor (ATR) over a total of 10 ATR cycles including seven normal cycles, one low power cycle, one unplanned outage cycle, and one Power Axial Locator Mechanism (PALM) cycle.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Reducing warpage in a hybrid large-scale additive manufacturing and compression molding process

In recent years, a hybrid manufacturing process, developed by combining extrusion-based large-scale additive manufacturing (AM) and compression molding (CM) techniques, has shown promising outcomes for producing structurally functional parts. The process can be used with both short fiber-reinforced composites and neat polymers and hence, even multi-material parts can be manufactured easily. This process offers the advantages of structural enhancement by having a desired fiber orientation using a large-scale AM process, as well as rapid manufacturing capability using a CM process. In the large-scale AM process, the alignment of fibers in the deposition direction enables significant improvement in the mechanical properties of the manufactured parts. However, the anisotropy resulting from the directional arrangement of fibers also introduces challenges related to warpage in the produced parts. This study aims to identify the causes of warpage and propose strategies to mitigate it. The research involves the use of preforms manufactured through the large-scale AM, which are then combined with the neat resin for CM manufacturing. A finite element-based numerical simulation model is developed, employing a sequentially coupled thermomechanical approach. Through a parametric study using the simulation models, optimization of printing direction and preform geometry is performed to minimize warpage. This contributes to the advancement and wider adoption of AM/CM hybrid manufacturing to produce structurally functional parts.

Jo, Eonyeon↗
DOE OSTI2025

Ceramography and Thermal Modelling of Irradiated TRISO Particles using BISON

This study presents a comprehensive analysis of irradiated TRISO fuel particles from the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program conducted at the Advanced Test Reactor (ATR) as part of my internship. Failure Analysis using high-resolution microscopy was performed to characterize mechanisms in TRISO particle layers, including the kernel, buffer, inner pyrolytic carbon (IPyC), silicon carbide (SiC), and outer pyrolytic carbon (OPyC) coatings. Failure modes were systematically classified. Thermal conductivity properties of each TRISO layer were evaluated through experimental measurements and computational models for use in BISON finite element code. Experimental thermal conductivity values ranged from 0.5 W/m·K for the buffer layer to 168 W/m·K for SiC, while BISON simulations provided averaged values of 0.57 W/m·K (buffer), 4.0 W/m·K (PyC layers), and 15.95 W/m·K (SiC). Temperature-dependent correlations from PARFUME were implemented for kernel thermal properties.

BISON↗
DOE OSTI2020

A novel framework for elastoplastic behaviour of anisotropic solids

A general framework for developing nonlinear hyperelastic/plastic constitutive laws for anisotropic solids experiencing large strains and strain rates has been developed. The proposed framework does not rely on the “a priori” known strain energy function, but instead introduces a physical decomposition of the material element into seven physically independent stress bearing mechanisms, each of which has a constitutive law in terms of internal moments described by a scalar function of a single variable. Furthermore, the model has been encoded into a combined finite-discrete element method and tested against static geomechanical test data. The numerical validation experiments show the model can reproduce plastic anisotropic behavior in both biaxial and uniaxial loading of a geomaterial.

58 GEOSCIENCES↗