Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element methods”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A parametric finite element study for determining burst strength of thin and thick-walled pressure vessels

To accurately predict the burst strength of both thin and thick-walled pressure vessels (PVs), a parametric study of PV burst strength was performed for a wide range of vessel geometries and materials using elastic-plastic finite element analysis (FEA). A valid FEA model was established through a detailed study of 2D versus 3D FEA models, the critical stress failure criterion versus the limit load criteria, and the thick-wall effect on the FEA simulations. Here, the results show that the stresses and strains at the mean diameter, rather than outside diameter, determines a more accurate burst strength for both thin and thick-walled PVs. On this basis, a parametrized FEA script using the ABAQUS Python application programming interface (API) was used to create a large database of PV burst strengths for a variety of vessel geometries and materials, demonstrating that Python scripting is a powerful technique for performing parametric studies or generating large databases. From the FEA results, using the regression method, a new burst pressure model was developed as a function of the vessel geometry (D/t ratio) and material properties (UTS and n). As validated by a large number of full-scale burst test data, the proposed burst model can very accurately predict the burst strength for both thin and thick-walled PVs.

42 ENGINEERING↗
DOE OSTI2026

ML-Based Pebble Power Reconstruction for Pebble Bed Reactor Analysis

Pebble power reconstruction has been explored to complement the conventional homogenized modeling approach in pebble bed reactor (PBR) analysis, as detailed heterogeneous geometry calculations are computationally expensive. The random distribution of pebble fuels within the core challenges the application of conventional pin power reconstruction methods. To address this, we introduce a machine learning approach based on the transformer model, composed of encoder and decoder layers, to estimate the flux and power form functions for reconstructing individual pebble neutron fluxes and powers. The homogeneous neutron flux distribution within each spectral zone (SZ) is obtained from finite element solutions of global diffusion or transport calculations. Verification tests demonstrate that the trained transformer model accurately predicts power form functions over a range of conditions, including variations in pebble enrichment, location, type, SZ size, and burnup. In particular, verification using a three-dimensional PBR benchmark with burned pebbles shows good agreement in heterogeneous pebble power distributions between Griffin and Serpent. These results highlight the potential of applying conventional pin power reconstruction approaches to PBR cores with randomly distributed pebbles.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Experimentally Validated Computational Modeling of Creep and Creep-Cracking for Nuclear Concrete Structures

In a Nuclear Power Plant, one of the most important components is the concrete nuclear reactor cavity, which serves both a structural and protective function as the biological radiation shield. Given that creep has been identified as a major knowledge gap in the assessment of nuclear structures (NUREG/CR-7153), this work helps to further the understanding of creep behavior of massive concrete containment structures for decades to enable safe and long-term operation of these facilities. This project has developed a robust, experimentally validated model to predict creep in nuclear concrete structures for up 60 years using short-term creep data thereby enabling a longer service life of critical facilities and early detection of structural failure. The work presented in this report is a pairing of computational and experimental methods. For the first time, the time temperature superposition (TTS) principle was successfully used to generate a uniaxial creep compliance master curve to predict mortar creep response for up to 22,500 days (nearly 60 years) at a reference temperature of 20°C. These data were used as input into finite element analysis (FEA) codes that use highly realistic random, 3D concrete microstructures from reconstructed coarse limestone aggregates. Finite element analysis performed provides the ability to quickly upscale mortar viscoelastic behavior to long-term concrete creep/relaxation data. A master creep compliance curve, constructed from the TTS principle, spanning 27 years, was used to validate two and a half decades of simulated concrete creep. Concurrently, three different simulated wall specimens were designed to mimic the behavior of post-tensioned concrete nuclear containment facility vessel walls over time as a result of concrete creep. The specimens were designed with different thicknesses, transverse and longitudinal reinforcement ratios, and level of post-tensioning stress. Each specimen contained various instrumentation to measure internal concrete temperature, concrete strain, and post-tensioning strain hourly for over 3 years. The concrete creep model developed in this project, based on the FEA concrete simulations, was applied to simulate the structural-scale experiments of prestressed concrete walls conducted in this project using the Grizzly code. These models can represent the effects of reinforcing and prestressing. Although there are some discrepancies with the experimental data, the model can predict the overall trends of the creep response in these experiments. One of these experimental models was also applied to an extended time to demonstrate how the findings from this study can be used to predict the behavior of actual structures of interest that have been in service for extended periods of time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Comments on “Failure analysis of corroded hydrogen-blended natural gas pipelines based on finite element analysis and genetic algorithm-back propagation neural network” [262 (2025) 111174]

This is a brief commentary paper to highlight and discuss the determination of hydrogen concentration in pipeline steel, effect of hydrogen embrittlement (HE) on the mechanical properties of the material, burst strength of corroded pipelines using finite element analysis (FEA) simulations, and curve-fit models for assessing remaining strength of X80 corroded pipelines for transporting hydrogen blended natural gas. Recently, Xie et al. [1] proposed a methodology to quantify the impact of HE on material properties and numerically determined burst pressure of X80 corroded pipelines. However, their HE quantification overestimated the degradation of tensile strength for hydrogen blending ratios beyond the original data range, and their FEA results of burst pressure are nonconservative. This work thus recharacterized the hydrogen concentration in the steel pipeline and the effect of HE on tensile strength, and then redetermined burst pressures for a set of typical corrosion defect cases considered by Xie et al. [1] based on an experimentally validated FEA modelling method. With the new FEA results, two empirical corrosion models were proposed for X80 corroded pipelines for hydrogen service. At zero hydrogen blending ratio, the novel empirical models predict burst pressures to be consistent with the industry-accepted corrosion models. Furthermore, both the numerical simulation method and the novel corrosion models are significant contributions to the pipeline industry and the hydrogen community. Application of these results will enhance the safety, reliability, and integrity of natural gas pipelines when used to transport hydrogen.

Burst pressure prediction↗
DOE OSTI2024

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Clifford Circuit Optimization with Templates and Symbolic Pauli Gates

The Clifford group is a finite subgroup of the unitary group generated by the Hadamard, the CNOT, and the Phase gates. This group plays a prominent role in quantum error correction, randomized benchmarking protocols, and the study of entanglement. Here we consider the problem of finding a short quantum circuit implementing a given Clifford group element. Our methods aim to minimize the entangling gate count assuming all-to-all qubit connectivity. First, we consider circuit optimization based on template matching and design Clifford-specific templates that leverage the ability to factor out Pauli and SWAP gates. Second, we introduce a symbolic peephole optimization method. It works by projecting the full circuit onto a small subset of qubits and optimally recompiling the projected subcircuit via dynamic programming. CNOT gates coupling the chosen subset of qubits with the remaining qubits are expressed using symbolic Pauli gates. Software implementation of these methods finds circuits that are only 0.2% away from optimal for 6 qubits and reduces the two-qubit gate count in circuits with up to 64 qubits by 64.7% on average, compared with the Aaronson-Gottesman canonical form.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Structure of a seeded palladium nanoparticle and its dynamics during the hydride phase transformation

Palladium absorbs large volumetric quantities of hydrogen at room temperature and ambient pressure, making the palladium hydride system a promising candidate for hydrogen storage. Here, we use Bragg coherent diffraction imaging to map the strain associated with defects in three dimensions before and during the hydride phase transformation of an individual octahedral palladium nanoparticle, synthesized using a seed-mediated approach. The displacement distribution imaging unveils the location of the seed nanoparticle in the final nanocrystal. By comparing our experimental results with a finite-element model, we verify that the seed nanoparticle causes a characteristic displacement distribution of the larger nanocrystal. During the hydrogen exposure, the hydride phase is predominantly formed on one tip of the octahedra, where there is a high number of lower coordinated Pd atoms. Our experimental and theoretical results provide an unambiguous method for future structure optimization of seed-mediated nanoparticle growth and in the design of palladium-based hydrogen storage systems.

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

Resonance ultrasound prediction of residual stress within a hybrid layer for additively manufactured samples

Hybrid additive manufacturing (AM) involves secondary processes or energy sources to alter specified locations within the build volume. Each hybrid step can refine the grain size, increase dislocation density, or modify residual stresses. Typically, the changes in mechanical properties are not confined within a single layer but have a compounding effect on preceding layers. Existing methods of measuring AM residual stress are limited in terms of their sensitivity, or they are destructive measurements. We propose using resonant ultrasound spectroscopy (RUS) to measure the residual stress in hybrid-AM components noninvasively, based on changes to the resonances, compared to a stress-free component. In this paper, we use finite element models to simulate residual stress in hybrid-AM components and to examine the sensitivity of RUS measurements in terms of frequency shifts and mode shapes with respect to single hybrid layers. Then, the RUS results are used to predict stress for a layer at a known location with unknown stress. Here, the approach highlights the capabilities of RUS to address an AM characterization challenge.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Framework for Compressing Unstructured Scientific Data via Serialization

We present a general framework for compressing unstructured scientific data with known local connectivity. A common application is simulation data defined on arbitrary finite element meshes. The framework employs a greedy topology preserving reordering of original nodes which allows for seamless integration into existing data processing pipelines. This reordering process depends solely on mesh connectivity and can be performed offline for optimal efficiency. However, the algorithm’s greedy nature also supports on-the-fly implementation. The proposed method is compatible with any compression algorithm that leverages spatial correlations within the data. The effectiveness of this approach is demonstrated on a large-scale real dataset using several compression methods, including MGARD, SZ, and ZFP.

Reshniak, Viktor [ORNL] (ORCID:0000000315454462)↗
DOE OSTI2022

Characterization of anisotropic mechanical properties of polymer composites from a hybrid additive manufacturing-compression molding process using x-ray computer tomography

In this study, anisotropic stiffness tensors were reconstructed based on fiber orientation distributions obtained from X-ray computer tomography (xCT). A preform was manufactured via a big area additive manufacturing (BAAM) system with carbon fiber (CF) filled acrylonitrile butadiene styrene (ABS). The tailored preform from additive manufacturing (AM) was used in the compression molding (CM) process to produce a low-void high-performance thermoplastic composite panel. An xCT technique was employed to detect the fiber orientations in CF/ABS composites manufactured via three different methods: AM from BAAM, extrusion compression molding (ECM), and AM-CM. The anisotropic stiffness tensor was obtained from the composite panel manufactured via the three manufacturing methods (AM, ECM, and AMCM). A micromechanics theory was used to obtain the orthotropic stiffness tensors of the composite panels and compared with the experimental values. The predicted stiffness tensors of AM and AM-CM composite panels were used to study the deformation characteristics of a steering wheel during airbag deployment by performing finite element analysis (FEA). The approach developed in this study can be utilized for evaluating high-performance composites.

Pokkalla, Deepak↗
DOE OSTI2021

Spectral signatures of transverse optical modes in semiconductor nanowires

In order to build nanophotonic devices, it is important to understand and ultimately control the optical mode structure within potential components such as nanoscale waveguides. However, experimental characterization of such modes in the optical regime is difficult due to the nanoscale dimensions of such components and the perturbations that would be induced by a near-field probe. Here, we demonstrate a probe-free, far-field method to characterize the optical modes within GaN nanowires (NWs) based on a novel off-axis scanning confocal microscope system. Using this microscope, we observe spectral signatures resulting from lateral leakage of waveguide modes when they exceed their respective cutoff limits. We identify these modes within hyperspectral images using an analytical model coupled with finite element simulations. The model can also be used to predict the spectral signatures for given geometrical and dielectric parameters, which enabled us to deduce the transverse dimension of the NW from hyperspectral images with an accuracy of ∼ <#comment/> 30 nm .

Wang, Maoji↗
DOE OSTI2021

Bound Constrained Partial DifferentialEquation Inverse Problem Solution by theSemi-Smooth Newton Method

We present the mathematical derivation, software implementation details, and computational results for a semi-smooth Newton method applied to two inverse problems governed by partial differential equations with bound constraints. The two problems share mathematical structural similarities to density-based topology optimization problems. The semi-smooth Newton method provides a mesh independent solution computation for the two test problems. A key step is that the complementarity part of the necessary optimality conditions are reformulated with the use of a complementarity functionφsuch that the complementarity conditions are satisfied if and only if a zero of a nonsmooth function has been obtained. The modular finite element package MFEM is utilized for the software implementation. In addition we constructed a matrix-free Operator to enable the use of efficient Krylov subspace IterativeSolver of MFEM for the solution of our two target problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Finite Element Analysis of the TRUST Nonlinear Dynamics Testbed

This paper builds on prior work conducted within the Los Alamos National Laboratory (LANL) Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) program. Specifically, it builds on finite element (FE) modeling efforts for the TRUST program’s Nonlinear Dynamics (ND) testbed. Historically, the FE model for the ND testbed has exclusively utilized a Lanczos eigensolver that linearly extracts the system’s natural frequencies. This paper investigates the Abaqus 2024’s explicit dynamic solver, which implements a central difference explicit solver. The central difference method used in Abaqus can capture nonlinear material responses in structural dynamic simulations making it suitable for the ND FE model.

42 ENGINEERING↗
DOE OSTI2022

Geomechanical and Hydrogeological Evaluation of a Shallow Hydraulic Fracture at the Devine Fracture Pilot Site, Medina County, Texas

UT-Austin’s Devine Fracture Pilot Site, 50 miles southwest of San Antonio, Texas, has been targeted for a comprehensive, multidisciplinary development of fracture diagnostic techniques that are cross-validated by ground-truth data acquisition near a recently created, 175-ft-deep, horizontal hydraulic fracture (Ahmadian et al. 2018 Demonstration of proof of concept of electromagnetic geophysical methods for high resolution illumination of induced fracture networks. In Proceedings of the SPE Hydraulic Fracturing Technology Conference and Exhibition, The Woodlands, Texas, USA, 23–25 January 2018. SPE-189858-MS.). To evaluate the fracture diagnostic methods at this site, we conducted injection tests with a predefined volumetric flow-rate profile, resembling a diagnostic fracture injection test on September 2020. Subsequently, we developed hydrogeological and geomechanical models based on flow-rate and bottomhole-pressure measurements. History-matching efforts using a simplified layer-cake hydrogeological model resulted in the field-scale formation permeability of 9.87 × 10 –15 m 2 (10 mD) and Darcy-scale fracture permeability. The analysis of the bottomhole pressure and injection-rate history showed that (1) the newly created horizontal fracture was closed adjacent to the injection well pre-injection and (2) the initial pump-pressure increase at a nominal volumetric injection rate led to near-well fracture reopening, fluid conductivity increase, and abrupt injection-rate increase. To overcome hydrogeological-model limitations of predicting fracture reopening throughout injection, we extended the modeling to a finite-element, poroelastic analysis of horizontal-fracture growth using a cohesive-zone model. Using this fracture-reopening model, we improved the history match of the transient-pressure response during the experiment by adjusting the hydromechanical properties. Furthermore, post-injection pressure transient analyses helped reduce uncertainty in the overburden-stress gradient, and the initial hydraulic-fracturing simulation verified the plausibility of achieving the surveyed propped fracture area.

02 PETROLEUM↗
DOE OSTI2021

Determining failure properties of as-received and hydrided unirradiated Zircaloy-4 from ring compression tests

An integrated experimental and computational approach was developed to determine the plastic strain at fracture and to investigate the fracture energy of as-received and hydrided unirradiated Zircaloy-4 nuclear cladding tubes in the hoop direction at room temperature from ring compression testing (RCT). This work builds on previous methods in the literature that were developed to obtain the mechanical properties of nuclear cladding before fracture (Young’s modulus, yield stress, and strain-hardening parameters) from experimental RCT results and extends the characterization of this material class to include material failure properties. Furthermore, a ductile damage approach was developed and implemented in a validated finite element model to predict material failure, evaluate fracture energy, and quantify the plastic strain at fracture initiation. The hydrogen content of the Zircaloy-4 specimens was varied up to 610 wppm when the failure became brittle. As expected, increasing hydrogen content caused embrittlement of the Zircaloy-4 samples tested and the plastic strain at fracture and fracture energy to decrease.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Active learning with multifidelity modeling for efficient rare event simulation

Here, while multifidelity modeling provides a cost-effective way to conduct uncertainty quantification with computationally expensive models, much greater efficiency can be achieved by adaptively deciding the number of required high-fidelity (HF) simulations, depending on the type and complexity of the problem and the desired accuracy in the results. We propose a framework for active learning with multifidelity modeling emphasizing the efficient estimation of rare events. Our framework works by fusing a low-fidelity (LF) prediction with an HF-inferred correction, filtering the corrected LF prediction to decide whether to call the high-fidelity model, and for enhanced subsequent accuracy, adapting the correction for the LF prediction after every HF model call. The framework does not make any assumptions as to the LF model type or its correlations with the HF model. In addition, for improved robustness when estimating smaller failure probabilities, we propose using dynamic active learning functions that decide when to call the HF model. We demonstrate our framework using several academic case studies (including some high-dimensional problems) and two finite element model case studies: estimating Navier-Stokes velocities using the Stokes approximation and estimating stresses in a transversely isotropic model subjected to displacements via a coarsely meshed isotropic model. Across these case studies, not only did the proposed framework estimate the failure probabilities accurately, but compared with either Monte Carlo or a standard variance reduction method, it also required only a small fraction of the calls to the HF model.

42 ENGINEERING↗
DOE OSTI2022

Surface coating by mechanofusion modulates bulk charging pathways and battery performance of Ni-rich layered cathodes

Ni-rich layered oxides as high-capacity battery cathodes suffer from degradation at high voltages. Here, we utilize a dry surface modification method, mechanofusion (MF), to achieve enhanced battery stability. The simplicity, high yield, and flexibility make it cost-effective and highly attractive for processing at the industrial scale. The underlying mechanisms responsible for performance improvement are unveiled by a systematic study combining multiple probes, e.g., 3D nano-tomography, spectroscopic imaging, in situ synchrotron diffraction, and finite element analysis (FEA). MF affects the bulk crystallography by introducing partially disordered structure, microstrain, and local lattice variation. Furthermore, the crack initiation and propagation pattern during delithiation are regulated and the overall mechanical fracture is reduced after such surface coating. We validate that MF can alter the bulk charging pathways. Such a synergic effect between surface modification and bulk charge distribution is fundamentally important for designing next-generation battery cathode materials.

25 ENERGY STORAGE↗
DOE OSTI2021

Thermal fluctuations and carrier localization induced by dynamic disorder in MAPbI 3 described by first-principles based tight-binding model

Halide perovskites are strongly influenced by large amplitude anharmonic lattice fluctuations at room temperature. We develop a tight binding model for dynamically disordered MAPbI$_3$ based on density functional theory (DFT) calculations to calculate electronic structure for finite temperature crystal structures at the length scale of thermal disorder and carrier localization. The model predicts individual Hamiltonian matrix elements and band structures with high accuracy, owing to the inclusion of additional matrix elements and descriptors for non-Coulombic interactions. We apply this model to electronic structure at length and time scales inaccessible to first principles methods, finding an increase in band gap, carrier mass, and the sub-picosecond fluctuations in these quantities with increasing temperature as well as the onset of carrier localization in large supercells induced by thermal disorder at 300 K. Finally, we identify the length scale $L^*= 5$ nm as the onset of localization in the electronic structure, associated with associated with decreasing band edge fluctuations, increasing carrier mass, and Rashba splitting approaching zero.

36 MATERIALS SCIENCE↗