Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element 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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Uncertainty Quantification Framework for Predicting Material Response with Large Number of Parameters: Application to Creep Prediction in Ferritic-Martensitic Steels Using Combined Crystal Plasticity and Grain Boundary Models

This paper presents an uncertainty quantification (UQ) framework for the physics-based model prediction of material response with a large number of parameters. The application problem presented in this work is that of predicting creep in Grade 91 steel at 600°C. The material response is defined with a physically based microstructural model with constitutive equations emulating several observed phenomena in Grade 91 and embodied into an explicit geometry mesoscale finite element model for prior austenite grains and grain boundaries. Creep within the grains and in grain boundaries are represented by crystal plasticity for dislocation motion and a physics-based model for cavity growth and nucleation, respectively. The creep behavior of this material is influenced by several parameters, some of which have a wide range of variation based on experimental data. UQ combined with microstructural modeling can discover the core microstructural causes of experimental variability, leading to improved materials with lower variability in critical long-term material properties. In this study, we investigate the model's uncertainty to identify material properties that may be modified during production to increase creep life and analyze different components of the crystal plasticity model for improvements. For this purpose, a quantity of interest is defined as time to minimum creep rate, which correlates well to the creep failure of the material. A deep neural network model was trained and validated to be used as a surrogate for the finite element model. Then, a variance-based sensitivity analysis is performed on the surrogate model to find the Sobol indices of the input parameters in respect to the output quantity of interest. The Sobol indices are used to reduce the dimensionality of the model. Generalized polynomial chaos expansion is used on the reduced basis models to propagate the uncertainty from the input parameters to the quantity of interest using the deep neural network surrogate model. These results are benchmarked against uncertainty propagation using Monte Carlo simulations. In conclusion, the UQ performed through the reduced basis model captures almost all the uncertainty in the model with significantly fewer simulations, making it possible to perform the UQ directly via simulations with the finite element model rather than surrogate machine-learned models.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Few-nucleon matrix elements in pionless effective field theory in a finite volume

Here, pionless effective field theory in a finite volume (FVEFT $\notπ$ ) is investigated as a framework for the analysis of multinucleon spectra and matrix elements calculated in lattice QCD (LQCD). By combining FVEFT π with the stochastic variational method, the spectra of nuclei with atomic number A ∈ {2, 3} are matched to existing finite-volume LQCD calculations at heavier-than-physical quark masses corresponding to a pion mass m $\notπ$ = 806 MeV , thereby enabling infinite-volume binding energies to be determined using infinite-volume variational calculations. Based on the variational wave functions that are constructed in this approach, the finite-volume matrix elements of various local operators are computed in FVEFT $\notπ$ and matched to LQCD calculations of the corresponding QCD operators in the same volume, thereby determining the relevant one- and two-body effective field theory counterterms and enabling an extrapolation of the LQCD matrix elements to infinite volume. As examples, the scalar, tensor, and axial matrix elements are considered, in addition to the magnetic moments and the isovector longitudinal momentum fraction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Topology-driven compressive behavior of Inconel 718 lattice structures with Z-strut reinforcement fabricated by laser powder bed fusion

This study investigates the compressive deformation behavior and mechanical performance of Inconel 718 lattice structures fabricated by laser powder bed fusion (LPBF). Four unit-cell topologies—BCC, FCC, BCCZ, and FCCZ—were designed with a unit-cell size of 3 mm and fabricated under identical process conditions to isolate the effect of topology. Measured relative densities ranged from 14.65% to 17.72%. Compressive testing showed that Z-strut-reinforced topologies (BCCZ: 54.6 MPa; FCCZ: 80.2 MPa) exhibited higher strength than their unreinforced counterparts, which may be associated with mixed-mode deformation behavior enabled by the vertically aligned Z-struts. Finite element simulations and Digital Image Correlation (DIC) analysis support the observation of a transition from node-dominated deformation in BCC/FCC to mixed-mode deformation in BCCZ/FCCZ. These findings suggest that unit-cell topology is a key design variable for tailoring deformation mechanisms in LPBF lattice structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Improvements to the Modeling of the TREAT Reactor and Experiments

This paper summarizes the latest improvements and lessons learned from the modeling and simulation of the transient test reactor at Idaho National Laboratory using the MAM- MOTH reactor physics application. MAMMOTH is a MOOSE-based, Finite Element Method application that specializes in the analysis of the spatial dynamics behavior of nuclear reactors. Since early 2018 several transient tests have been conducted at TREAT, thus providing the opportunity to apply and benchmark modern modeling and simulation tools. MAMMOTH was used to provide predictions of the power coupling factor between the core and the experiment for various experiments. Even though the power coupling factor predictions agree very well with the experimental data, within the bounds of the experimental uncertainty, one shortcoming was the underprediction of the total energy deposited in the core and in the experiment. Determination of the sources for this discrepancy are ongoing, but several key problems have been identified and resolved, thus providing valuable insights for future research. This paper discusses several of these lessons learned. First, the heat capacity data for the TREAT fuel has some significant problems due to limitations of the measurement techniques used circa 1960s. The sensitivity of the peak power and the total energy deposition to various representations of the heat capacity is approximately 5%. Second, the effects of the biological shield and thermal column on the modeling of the core are non-negligible, since they affect the mean generation time and the effective reflection of neutrons back into the core, which is suspected to be important during the core heat up. Matching the reactor period resolves the fact that the reduced spatial domain used in the MAMMOTH model underpredicts the mean generation time. The neutron reflection from these regions is marginally improved with the use of an albedo boundary condition. Third, modeling of the control rod movement with a multi-scheme method is introduced and its current limitations are exposed. Fourth, we explore the effects of using a homogenized model with Superhomogenization equivalence and how that differs from fully heterogeneous simulations. Finally, the energy condensation effects for this graphite core are significant. Solutions with 10 and 26 energy groups show the benefits of using a finer coarse group structure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

MOOSE Framework Meshing Enhancements to Support Reactor Analysis

MOOSE-based physics codes require an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The use of difficult-touse, external licensed software is often required to generate high quality meshes for reactor geometries. High-fidelity geometry modeling also requires elaborate tracking of groups of elements for material property assignment and output reporting which can be considerably complex for the user to identify and maintain. Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, several meshingrelated enhancements have been developed for the MOOSE framework to address user challenges in creating finite element meshes for advanced reactor geometries. MOOSE mesh generators have been developed to mesh hexagonal geometries (pins, ducted assemblies, and cores) commonly found in liquid-metal cooled fast reactor concepts. The mesh generator used for hexagonal pin cells is generic for regular polygons and therefore may also be used for Cartesian pin cells. Hexagonal pin cells can be stitched into ducted assemblies, and assemblies can be stitched together into a core. The user may specify region ids, region names, and other preferences on the mesh. This control is useful for later material mapping in the MOOSE-based physics codes input. A capability was also developed for meshing rotating control drums including determination of material volume fractions in each mesh element as a function of time. Control drum meshes may be stitched to other hexagonal assemblies to create a core configuration. Additional mesh generators were developed that wrap around the hexagonal meshing capabilities and utilize “extra element integer” ID values on each element. In regular Cartesian or hexagonal assemblies or cores, the bookkeeping of element groups for both material assignment and output reporting can now be automated through assignment of pin, assembly, core, axial and depletion id values stored as extra element integers. The extra element tags on the mesh greatly speed the reactor analyst’s efforts to map materials to meshes, track depletion zones, and parse output such as axial pin power distributions. At the highest level, pin, assembly, and core mesh generators (with this reactor terminology) have also been developed to easily generate regular Cartesian and hexagonal cores, including axial extrusion. These reactor geometry builders call upon the previously mentioned capabilities to produce analysis-ready 3D meshes including material assignments. Open source mesh triangulation capabilities were also investigated for integration into the MOOSE framework to address the need for meshing the core periphery region which extends from the irregular outer assembly border to a cylindrical boundary. Options are limited due to licensing constraints, and the recommendation is pursue building a native MOOSE Delaunay triangulator routine with full functionality. Finally, a series of verification problems were performed with NEAMS physics tools. All developed capabilities will be available in the new open-source “Reactor” module of the MOOSE framework, which is accessible to any MOOSE-based NEAMS physics tool.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Radioactive Source Localization via Bayesian Particle Filter

In the event of a misplaced radioactive source or other emergency situation, measuring a radiation field, mapping its distribution, and determining a source location are essential tasks to ameliorating the situation. However, radiation fields may be extremely hazardous to human surveyors and minimizing received radiation doses is just as essential. Robots appear to be a potential solution to these problems. Beyond simply measuring radiation, the robot's computer processing capabilities offer a way to apply complex data analysis methods to radiation measurements in real-time. Methods which predict likely source locations can then feed this information into other processes, potentially improving path planning and enabling more efficient measurements. Given a robot mounted with a gamma-ray detector, can we: develop a methodology to account for detector performance across a wide range of source angles, distances, and photon energies? operate an autonomously navigating robot to effectively survey and characterize an area of interest? implement a data analysis method, conventionally used in measurements of motion, for source localization purposes? An open-source TurtleBot 3 robot, running Robot Operating System (ROS) on Ubuntu 16.04 LTS, was fitted with a Kromek GR1{sup R} Cadmium Zinc Telluride (CZT) solid-state gamma-ray detector. As a part of ROS, the packages OpenSlam, gmapping, and amcl were used to perform Simultaneous Localization and Mapping (SLAM), determining the robot's position and mapping the surrounding area. Data was acquired via Lidar mounted on top the TurtleBot 3. Detector Calibration Fit: The equation was fit to 365 counts of various energies, distances, and angles. A MATLAB{sup R} program was written to simulate measurements taken a robot on a random walk, with count data and positions discretized into finite element pixels. Using this program, a sample of 100 runs was performed on a map with a simulated source at the center, with a total of 200 of 2 pixels each. Similarly, multiple runs of the filter were performed on recorded robot measurement data. In both simulation and real tests, when corrected for errors (particles placed outside of bounds or on the robot, and simulation-specific errors), corresponding t-tests of predicted x and y-coordinates were within a 95% confidence interval of the actual position. For the real trial, these positions are slightly skewed right in the x-axis as the robot remained mainly to the left side of the source within the sample area. These simulations demonstrate potential validity for the usage of a particle filter as method of radioactive source localization. In the future, true real-time implementation and data fusion may further augment the performance of the robot to localize lost sources. Additionally, identification of multiple sources, determination of source types, and usage of a collimator are areas to potentially be explored.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

A Combined Computational and Mathematical Analysis of Interconnect Fatigue Potential in Photovoltaic Modules

A finite element model of a 60-cell monocrystalline silicon glass-polymer photovoltaic module was simulated with ±1.0 kPa and ±2.4 kPa loads applied to the glass to calculate the deformation under load. Cell-to-cell displacements were used to approximate interconnect strain and stress. A mathematical fatigue cycle life relation was fitted to data for the interconnect material (copper), to generate a life prediction at each interconnect location based on the local stress means, reversal extents, and amplitudes. Interconnect stress was found to be significantly asymmetric about zero despite symmetric positive and negative module loads due to laminate thickness offsets about the neutral plane and the effects of module framing. Cycle life results indicated that interconnect fatigue failure was unlikely to occur over a 30-year lifetime of conservative wind and snow load cycles since the typical cell design feature of leaving some unconstrained length between the cell edge and first solder pad increases the effective gauge length and decreases the stress levels below the material endurance limit. Follow-up analyses found that 3.6 mm and 6.4 mm were the minimum unconstrained lengths required to survive the assumed lifetime of wind and snow cycles, respectively, confirming that typical industrial module constructions with 8–15 mm unconstrained lengths should survive conservatively. Notably, large magnitude, low-cycle snow loading was consistently the limiting factor requiring a longer unconstrained interconnect length. Finally, insights and workflows from this study inform module interconnection design limits for survival against mechanical fatigue in deployment environments.

14 SOLAR ENERGY↗
DOE OSTI2025

Direct sensitivity analysis on the parameterization of crystal plasticity models

Various methods for calibrating crystal plasticity finite element (CPFE) models lead to non-unique input parameter values, which subsequently introduce uncertainty in the predicted mechanical response. Sensitivity analysis (SA) conducted on crystal plasticity models is used to identify how variability in these parameters contribute to output uncertainty. Traditional SA on CPFE parameters uses simplified surrogate models to save computational time. However, the accuracy of the surrogate models depends on the quantity of training data used, and any modeling error can propagate into the SA results, potentially affecting their reliability. In this work, the elementary effects test (EET) method, a global SA technique using direct CPFE simulations was employed, and the results obtained were compared with the First Order Second Moment (FOSM) method. ExaConstit, an open-source GPU-enabled CPFE code, was used to perform the simulations and direct SA. The EET method was accurately able to capture the non-linear effects of all the input parameters on the output and is a valuable approach for reliably attributing parameter sensitivities in CPFE models. Based on the results, efficient strategies to perform future parameter calibration and SA are discussed. Additionally, the SA trends observed in different single crystal orientations closely mirrored the activity of the slip systems.

Elementary Effects Test↗
DOE OSTI2022

Verification of MOOSE/Bison's Heat Conduction Solver Using Combined Spatiotemporal Convergence Analysis

Bison is a computational physics code that uses the finite element method to model the thermo-mechanical response of nuclear fuel. Since Bison is used to inform high-consequence decisions, it is important that its computational results are reliable and predictive. One important step in assessing the reliability and predictive capabilities of a simulation tool is the verification process, which quantifies numerical errors in a discrete solution relative to the exact solution of the mathematical model. One step in the verification process—called code verification—ensures that the implemented numerical algorithm is a faithful representation of the underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. In this paper, the code verification process is applied to spatiotemporal heat conduction problems in Bison. Simultaneous refinement of the discretization in space and time is employed to reveal any potential mistakes in the numerical algorithms for the interactions between the spatial and temporal components of the solution. For each verification problem, the correct spatial and temporal order of accuracy is demonstrated for both first- and second-order accurate finite elements and a variety of time-integration schemes. Furthermore, these results provide strong evidence that the Bison numerical algorithm for solving spatiotemporal problems reliably represents the underlying mathematical model in MOOSE. The selected test problems can also be used in other simulation tools that numerically solve for conduction or diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Stochastic nonlinear analysis of unidirectional fiber composites using image-based microstructural uncertainty quantification

Here, we present a data-driven nonlinear uncertainty quantification and propagation framework to study the microstructure-induced stochastic performance of unidirectional (UD) carbon fiber reinforced polymer (CFRP) composites. The proposed approach integrates (1) microscopic image characterization, (2) stochastic microstructure reconstruction, and (3) efficient multiscale finite element simulations enabled by self-consistent clustering (SCA) analysis. To model the complex microstructural variability, the proposed UQ methods take the non-Gaussian uncertainty sources into account through a distribution-free sampling approach leveraging nonparametric and asymptotic statistical tools. A hierarchical conditional sampling strategy enables the simultaneous sampling of multiple sources of uncertainties. Our approach provides insights into the impact of microstructural variabilities, which are shown to have an increasing impact on the nonlinear responses of UD CFRP parts under progressive compression loading and ultimately on the failure rate over time. We discover that before CFRP parts start to fail, a characteristic time period emerges with distinctive uncertainty distributions specific to the microstructure variability and the probability of failure. Identifying the failure time period is crucial to the reliability prediction, which is an essential component of CFRP design.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sockeye: A One-Dimensional, Two-Phase, Compressible Flow Heat Pipe Application

Sockeye is a heat pipe analysis application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework. The primary purpose of Sockeye is to provide a transient heat pipe simulation tool to be used in the analysis of nuclear microreactor designs. Sockeye provides the capability to perform one-dimensional, two-phase, compressible flow simulation of a heat pipe working fluid and two-dimensional, axisymmetric heat conduction for the heat pipe cladding and its surroundings. Sockeye is demonstrated against analytical solutions and experimental data from the SAFE-30 heat pipe module test.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Nanoscopic Imaging of Self-Propelled Ultrasmall Catalytic Nanomotors

Ultrasmall nanomotors (<100 nm) are highly desirable nanomachines for their size-specific advantages over their larger counterparts in applications spanning nanomedicine, directed assembly, active sensing, and environmental remediation. While there are extensive studies on motors larger than 100 nm, the design and understanding of ultrasmall nanomotors have been scant due to the lack of high-resolution imaging of their propelled motions with orientation and shape details resolved. Here, we report the imaging of the propelled motions of catalytically powered ultrasmall nanomotors─hundreds of them─at the nanometer resolution using liquid-phase transmission electron microscopy. These nanomotors are Pt nanoparticles of asymmetric shapes (“tadpoles” and “boomerangs”), which are colloidally synthesized and observed to be fueled by the catalyzed decomposition of NaBH4 in solution. Statistical analysis of the orientation and position trajectories of fueled and unfueled motors, coupled with finite element simulation, reveals that the shape asymmetry alone is sufficient to induce local chemical concentration gradient and self-diffusiophoresis to act against random Brownian motion. Our work elucidates the colloidal design and fundamental forces involved in the motions of ultrasmall nanomotors, which hold promise as active nanomachines to perform tasks in confined environments such as drug delivery and chemical sensing.

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

Toward large-scale simulation of residual stress and distortion in wire and arc additive manufacturing

In this study, we aim to advance the structural analysis of wire and arc additive manufacturing (WAAM) by considering the thermomechanical features inherent in direct energy deposition. Simulation approaches including the iterative substructure method (ISM), dynamic mesh refining method (DMRM), and graphics processing unit (GPU) based explicit finite element method (FEM) were developed for accelerating additive manufacturing stress analysis that is very time consuming by conventional numerical methods. ISM and DMRM take advantage of the strong nonlinearity phenomenon near the moving heat source by reducing the global iterations and refining the local mesh, respectively. In addition to the use of GPUs, the explicit FEM is accelerated by a time scaling technique based on the inherent strain concept. The residual stress and distortion of two large builds were analyzed, showing very consistent numerical results and good agreement with experiments. Compared with the commercial software Abaqus, the novel approaches reduced the computational cost substantially without compromising accuracy. The GPU code showed the highest computational efficiency (30~70×), while DMRM and ISM had acceleration factors of 9× and 3×, respectively. Finally, such high-fidelity modeling approaches will be very useful for building up a digital twin of WAAM to reduce development time and cost.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Measurements of the thermal conductivity of reference liquids using a modified transient hot-wire needle probe

Many power generation processes feature some type of working fluid that plays a critical role in the transfer of heat throughout the system. To optimize thermal efficiencies, an accurate characterization of the temperature-dependent thermal conductivity of the fluid is needed. Herein, this study develops a time-domain Modified Transient Hot-wire Needle Probe to overcome phenomena that have impaired traditional measurement techniques when the fluid is both corrosive and at an elevated temperature. First, the physical design of the instrument is presented to demonstrate features that limit interference from corrosion, electrical conduction, and convective and radiative heat transfer. Next, a governing equation for 1D (radial) heat conduction is fully developed using the thermal quadrupoles technique. A finite-element study is performed for validation. A sensitivity analysis is done to assess the time-dependent influence of critical parameters in the governing equation for several samples. Finally, thermal conductivity measurements are presented for seven different liquid substances by fitting the thermal model to experimental data across a range of temperatures. Reference correlations are included for comparison. The sample thermal conductivity and ambient temperature of the experimental conditions range from 0.125–16.5$\frac{W}{m·K}$ and 10–400°C, respectively. The total uncertainty, comprised of both the experimental design and fitting process, of all measurements is quantified.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

CU-BENs: A structural modeling finite element library

The present work discusses capabilities within the finite element library CU-BENs. CU-BENs focuses on applying the finite element method to structural mechanics problems encountered within the context of inverse problems and partitioned fluid–structure interaction; thus element formulations are primarily of a structural type — truss, frame, and triangular discrete Kirchhoff theory shells. CU-BENs defaults to the skyline sparse storage scheme for the system matrix, but also supports other storage schemes when using external libraries such as LAPACK and UMFPACK. CU-BENs includes built-in nonlinear solution strategies, such as the Newton Raphson method and modified spherical arc length method, that are available within static analyses as well within the context of transient dynamic analyses involving a generalized-α implementation of the Newmark implicit time integration scheme.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Condenser Enclosure Vibration Isolation Mount Selection

A mini split cooling system will be used to maintain temperature requirements in a mobile secure transport system. While the split cooling system was designed to be used in static residential or commercial applications, it was selected for this transportation application due to a unique set of security requirements. However, the system’s ability maintain reliability and survive prolonged long-term shock and vibration is a significant concern. The mitigation strategy is to select vibration isolation mounts and perform lifetime shock and vibration testing to demonstrate survivability. The goal of this study is to generate a finite element model of the system and perform modal analysis to inform selection of vibration mounts to minimize the amount of vibrational energy transferred to the split cooling system. The scope of this report is limited to study of the condensing unit only, and geometric variation of the assembly will not be allowed.

42 ENGINEERING↗
DOE OSTI2020

Discrete spherical harmonic functions for texture representation and analysis

A basis of discrete harmonic functions for efficient representation and analysis of crystallographic texture is presented. Discrete harmonics are a numerical representation of the harmonics on the sphere. A finite element formulation is utilized to calculate these orthonormal basis functions, which provides several advantageous features for quantitative texture analysis. These include high-precision numerical integration, a simple implementation of the non-negativity constraint and computational efficiency. Simple examples of pole figure and texture interpolation and of Fourier filtering using these basis sets are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2022

CAMERA: A method for cost-aware, adaptive, multifidelity, efficient reliability analysis

Estimating probability of failure in aerospace systems is a critical requirement for flight certification and qualification. Failure probability estimation involves resolving tails of probability distributions, and Monte Carlo sampling methods are intractable when expensive high-fidelity simulations have to be queried. Here, we propose a method to use models of multiple fidelities that trade accuracy for computational efficiency. Specifically, we propose the use of multifidelity Gaussian process models to efficiently fuse models at multiple fidelity, thereby offering a cheap surrogate model that emulates the original model at all fidelities. Furthermore, we propose a novel sequential acquisition function based experiment design framework that can automatically select samples from appropriate fidelity models to make predictions about quantities of interest at the highest fidelity. We use our proposed approach in an importance sampling setting and demonstrate our method on the failure level set and probability estimation on synthetic test functions and two real-world applications, namely, the reliability analysis of a gas turbine engine blade using a finite element method and a transonic aerodynamic wing test case using Reynolds-averaged Navier-Stokes equations. We show that our method predicts the failure boundary and probability more accurately and at a fraction of the computational cost compared with using just a single expensive high-fidelity model. Finally, we show that our sequential approach is guaranteed to asymptotically converge to the true failure boundary with high probability.

97 MATHEMATICS AND COMPUTING↗