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.
A Condensed Constrained Nonconforming Mortar-Based Approach for Preconditioning Finite Element Discretization Problems
This paper presents and studies an approach for constructing auxiliary space preconditioners for finite element problems using a constrained nonconforming reformulation that is based on a proposed modified version of the mortar method. The well-known mortar finite element discretization method is modified to admit a local structure, providing an element-by-element or subdomain-by-subdomain assembly property. This is achieved via the introduction of additional trace finite element spaces and degrees of freedom (unknowns) associated with the interfaces between adjacent elements or subdomains. The resulting nonconforming formulation and a reduced-via-static-condensation Schur complement form on the interfaces are used in the construction of auxiliary space preconditioners for a given conforming finite element discretization problem. Overall, the properties of these preconditioners are studied and their performance is illustrated on model second order scalar elliptic problems utilizing high order elements.
Validation of high-fidelity ion cyclotron range of frequencies antenna coupling simulations in full 3D geometry against experiments in the ASDEX Upgrade tokamak
We address the validation of finite elements ion cyclotron range of frequencies (ICRF) antenna coupling simulation against experiments performed in the ASDEX Upgrade tokamak. Measurements of the loading resistance in ICRF-heated, magnetically-perturbed 3D plasma discharges are compared against numerical predictions of the RAPLICASOL code. To this end, the 3D induction field and the 3D density profile are modeled by concatenating the PARVMEC, BMW and EMC3-EIRENE codes. The 3D density is input to RAPLICASOL, where full-wave simulations are performed on a finite element mesh retaining full 3D geometry in the ICRF antenna model and the plasma description. The results are further compared with RAPLICASOL simulations using a 1D density profile as measured at the outboard midplane in the same experiments. We find that simulations using a 1D density profile overestimate the change in loading resistance by a factor of ~197 -248%, while simulations using the full 3D density profile are in agreement with experiments within a factor ~16%.
Radiation-induced bowing of SiC/SiC composites under neutron flux gradients—integral experimental data for model validation
Here, the radiation-induced swelling of SiC and its composites, including strong dependencies on temperature and dose, can drive significant lateral bowing in the presence of temperature and/or dose gradients. In recent years, simulations have been performed to assess the extent of bowing in SiC composite light-water reactor (LWR) fuel cladding and boiling water reactor (BWR) channel boxes. However, to date, no integral experimental data exist to validate these models. This work provides the first experimental bowing evaluation of three ∼380 mm long SiC composite specimens irradiated under varying neutron dose gradients (∼50°C–60°C, 0.03–0.06 dpa): two tubes (∼9.8 mm diameter) and a miniature BWR channel box (∼30 mm square). The measured radiation-induced length swelling (∼0.3%–0.7% linear) was consistently 10%–21% higher than values obtained from 3D finite element structural analyses with inputs from 3D radiation transport calculations. This discrepancy could be at least partially explained by differences in dose rate (∼10 -8 dpa/s) compared to the literature data (∼10-6 dpa/s) used to establish the dose-to-swelling correlations in the model. Nevertheless, the modeled bowing magnitudes (<2 mm) obtained from finite element analyses and simple analytical equations were within the bounds of the experimental measurements for all specimens. With improved confidence in the ability to predict the structural response and measure the macroscopic deformations, future experiments will target transient bowing under neutron flux gradients at representative LWR temperatures and assess whether grid spacers can mitigate the tens of millimeters of bowing that would otherwise be expected in ∼4 m long LWR components.
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.
Verification and Regression Testing of a Physically Stabilized Layered Solid Element Formulation in DYNA3D/Paradyn
Recent modeling of filament-wound composite structures drove the need for a new element type within the finite element code DYNA3D/Paradyn. This new layered solid element was implemented and designed to capture the kinematics and constitutive be havior of various lamina layers defined with arbitrary orientations and volume fractions within the element to accurately model a laminated composite material. The new el ement uses a single integration point in each of the lamina layers defined to capture the constitutive behavior, and a novel physical stabilization routine is used to prevent hourglassing while mitigating shear and volumetric locking. This report documents the verification testing and the regression tests added to the DYNA3D/Paradyn software quality assurance test suite to assess the proper implementation of the new layered solid element. There were a total of 27 tests added the DYNA3D test suite, which consist of problems using the layered solid element with a single layer or multiple layers for isotropic and orthotropic material models. These tests are a combination of simple kinematically driven patch tests, beam bending tests, plate bending tests, and more complicated problems. Ultimately, the testing done verifies that the element behaves as expected and is implemented correctly.
Crystal plasticity finite element simulation of lattice rotation and x-ray diffraction during laser shock compression of tantalum
Here we present a crystal plasticity model tailored for high-pressure, high-strain-rate conditions that uses a multiscale treatment of dislocation-based slip kinetics. We use this model to analyze the pronounced plasticity-induced lattice rotations observed in shock-compressed polycrystalline tantalum via in situ x-ray diffraction. By making direct comparisons between experimentally measured and simulated texture evolution, we can explain how the details of the underlying slip kinetics control the degree of lattice rotation that ensues. Specifically, we show that only the highly nonlinear kinetics caused by dislocation nucleation can explain the magnitude of the rotation observed under shock compression. We demonstrate a good fit between our crystal plasticity model and x-ray diffraction data and exploit the data to quantify the dislocation nucleation rates that are otherwise poorly constrained by experiment in the dynamic compression regime.
Islet: interpolation semi-Lagrangian element-based transport
Abstract. Advection of trace species, or tracers, also called tracer transport, in models of the atmosphere and other physical domains is an important and potentially computationally expensive part of a model's dynamical core. Semi-Lagrangian (SL) advection methods are efficient because they permit a time step much larger than the advective stability limit for explicit Eulerian methods without requiring the solution of a globally coupled system of equations as implicit Eulerian methods do. Thus, to reduce the computational expense of tracer transport, dynamical cores often use SL methods to advect tracers. The class of interpolation semi-Lagrangian (ISL) methods contains potentially extremely efficient SL methods. We describe a finite-element ISL transport method that we call the interpolation semi-Lagrangian element-based transport (Islet) method, such as for use with atmosphere models discretized using the spectral element method. The Islet method uses three grids that share an element grid: a dynamics grid supporting, for example, the Gauss–Legendre–Lobatto basis of degree three; a physics parameterizations grid with a configurable number of finite-volume subcells per element; and a tracer grid supporting use of Islet bases with particular basis again configurable. This method provides extremely accurate tracer transport and excellent diagnostic values in a number of verification problems.
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.
A machine-learning-aided data recovery approach for predicting multi-material thermal behaviors in advanced test reactor capsules
Instrumented experiments conducted at test reactors are essential to the deployment of new advanced reactor systems. Designing new experiments and generating data on specific reactor conditions require significant investments in terms of both time and cost. Finite element analysis software can be used to create high-fidelity models of experiment environments in order to support the actual experiments, but computation time remains a concern in terms of applying outcomes to real-time usage of data (e.g., a digital twin [DT]). Here, the present research proposes a machine-learning (ML) aided approach to making temperature and displacement predictions based on the thickness of the outer gas gap on the experimental capsule used for in-pile demonstration of a novel new thermal conductivity probe in the Advanced Test Reactor (ATR). This capsule consisted of U10Zr fuel, a rodlet, sodium, and inner and outer capsules. Gas gaps existed between the fuel and the rodlet, and between the inner and the outer capsule. The learning data pertained to an experimental capsule's radial distributions of temperature and displacement, as obtained based on Abaqus and the physical features. For the first step of ML sequence, the temperature was predicted using three positional parameters. Next, the displacement was predicted using seven additional parameters. Each physical feature was normalized in order to be both nondimensional and standardized. The temperature and displacement predictions showed good agreement with the simulation results in all cases involving interpolation and extrapolation. Furthermore, data similarity enhancement increased the similarity between the training and the target data, thereby increasing the predictive accuracy of the ML models. In certain extrapolation cases involving limited original ML model accuracy, data similarity enhancement and data recovery was able to somewhat improve this accuracy.
pyNuMAD v.0.1
SAND2024-08606O The pyNuMAD software is used for managing wind turbine blade model data. pyNuMAD specializes in defining the geometry, materials, and boundary conditions for structural analysis of wind turbine blades. This includes loading in data files and providing an interface for users to make updates to the model. The software also features meshing functionality, which takes the blade model and creates a shell or brick mesh for use in finite element analysis. A typical user workflow might be: load in blade information from a yaml file, make adjustments to the blade properties, update the blade based on the adjustments, create a mesh of the blade, export this blade to another software for structural analysis. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
Sample problems for Section III, Division 5 design by inelastic analysis of Grade 91 components
This report works through two sample design analyses of representative Grade 91 components using the ASME Boiler & Pressure Vessel Code Section III, Division 5, Appendix HBB-T rules for the strain and deformation and creep-fatigue design criteria using a design by inelastic analysis approach. The component geometries and loadings were selected to span a wide variety of potential applications. The sample components are an endplate geometry, which is fairly representative of typical pressure vessel construction, and a core block, representing a critical component in a heat pipe microreactor, which has a complicated geometry and complex thermal stress history. These sample problems serve several functions. First, they are complete worked examples of how to apply the current design by inelastic analysis rules. Second, they serve as test cases to evaluate the current rules and provide feedback on how they could be optimized to deal with complicated component geometries or to be better automated in modern finite element analysis software. Finally, the analysis here uses an inelastic constitutive model for Grade 91 developed by Argonne National Laboratory and proposed for inclusion in the ASME Boiler & Pressure Vessel Code Section III, Division 5 rules as a reference constitutive model. Using this constitutive model in the sample problem analyses tests the model in the context of the ASME Code rules. The final conclusions of this report are that the material model and the current Code rules are satisfactory, but that the Code rules could be optimized to simplify automation and reduce over conservatism. Future work should address the specific optimization topics discussed in this report.
Advancements Toward ASME Nuclear Code Case for Compact Heat Exchangers
Our research team proposes to advance the state of the ASME section III code (nuclear service) for Compact Heat Exchangers (CHX). This work will improve the technical state of CHXs and lay the foundation necessary for these heat exchangers to be certified for use in nuclear service. During the course of this work, we will advance the understanding of the performance, integrity, and lifetime of the CHXs for use in any industrial application, making their use more attractive and accessible to the industry. We will do this by developing qualification and inspection procedures that utilize Non-destructive evaluation (NDE) and advanced in-service inspection techniques, with insight from the industrial utility leader EPRI. We have enlisted colleagues at MPR Associates (MPR), an elite nuclear code consulting firm, who are experts on the ASME section III code and who, with input from members of the ASME section III committee, will direct the testing and help develop a series of documents that define the rules and regulations for use of the CHX. Colleagues at North Carolina State University (NCSU) and Oregon State University (OSU) will conduct extensive tensile, creep, and fatigue experiments on diffusion bonded samples (manufactured by US-based Vacuum Process engineering) along with modeling using the elastic perfectly plastic assumptions and comprehensive full inelastic finite element analysis (FEA). This work will allow analysis by design and confidence in the strength of different internal structures. To ensure industry acceptance and long term confidence, team members at the University of Wisconsin–Madison (UW), University of Michigan (UM), Georgia Tech (GT), and the University of Idaho (UI) will extensively test prototypic heat exchangers manufactured by US-based manufactures CompRex and Vacuum Process Engineering (VPE), a leader in the development of advanced CHX. This testing will include the use of various working fluids (salt, sodium, helium, and sCO2) to evaluate operational issues as well as structural integrity under the most severe conditions. Post-test analysis of the tested CHXs coupled with pre/in-service/post NDE (ultrasonic and radiography) led by the Electric Power Research Institute (EPRI) will be incorporated into the development of the rules and regulations for their use in nuclear service.
Percheron Power Archimedes Screw Turbine Analysis (Abstract)
The development of the Percheron Power (Percheron) composite Archimedes screw turbine (AST) is based upon employing newly developed theoretical mathematical models for optimizing the water volume carried through the AST for a given diameter, combined with minimizing efficiency losses due to bypass flow and friction. The key feasibility issues associated with producing the composite blades will be to understand the strength, fatigue, and wear performance of the AST blades and associated components, and how these factors change along the length of the turbine. To get to this understanding of material property requirements, Percheron needs free-surface computational fluid dynamics (CFD) models of the AST designs and to perform interactive modeling between the turbine blade optimization and component design based on finite element method stress analysis (FEA). The workscope and analyses requested by Percheron in the Small Business Voucher (SBV) program are well aligned with PNNL capability and recent work experience in performing CFD and FEA analysis for hydropower and marine energy turbines. The CFD analysis will provide estimates of power, efficiency, and surface forces on the blades. The forces on the blade will be transferred to the FEA analysis that will be used to evaluate the material stresses. The CFD and FEA analyses will be repeated for different AST designs provided to PNNL by Percheron. The work conducted in this SBV project will complement the existing Percheron FOA project by providing better defined parameters for design and material selection for the turbine unit prior to manufacture of initial test units.
Optimal Membrane Cascade Design for Critical Mineral Recovery through Logic-based Superstructure Optimization
In this work, we extend the superstructure model proposed by Wamble et al. (2022) that considers feed input locations, recycling strategies, split fractions, stage numbers, and membrane area. We include the total number of stages as a decision variable, which might be particularly useful when there is cost as- sociated with adding additional stages. We propose a Generalized Disjunctive Programming (GDP) superstructure model that integrates all the design variables of the system. We also investigate the scalability of the model by varying the number of stages and the number of finite elements per stage to determine the impact on recovery and solution time.
Advanced System Thermal Fluids Solver Development for SAM
This work summarizes a feasibility study on testing numerical algorithms that are suitable and efficient for advanced system analysis code development under the mutli-physics framework, MOOSE. The key is the implementation of a high-order one-dimensional staggered-grid finite volume method (SG-FVM), and its direct interaction with the linear/nonlinear solver, PETSc. Leveraging the existing capabilities of the SAM code, significant code coverages were established in the finite volume method code. This in turn allows for a suite of test problems with different problem sizes and levels of complexity to be used to quantify the performance improvement of the finite volume method code. As evidently shown in this study, the implemented SG-FVM demonstrated superior performance improvement against a direct finite element method implementation through MOOSE for the wide range of selected problems. On two computer systems, the speedup was observed to be significant, with at least one order of magnitude of solving time reduction. In addition, for a complex reactor model, transient simulation was performed using the finite volume method code, the results of which agree very well with the reference results from the finite element method code. Overall, this study demonstrates a successful feasibility study on the proposed numerical algorithms and software structure to support advanced system analysis tool development. In this work, short-term priority development and testing items were identified, and long-term code adoption and integration plans were made for the eventual deployment of the finite volume method in the SAM code.
Machine tool cross beam design, fabrication, and testing using metal big area additive manufacturing
This paper describes the application of metal Big Area Additive Manufacturing (mBAAM) to the fabrication of a machine tool cross beam. The replacement of a traditional box design weldment with a new design printed by wire arc additive manufacturing using the MedUSA system at Oak Ridge National Laboratory (ORNL) is detailed. This requires a new design strategy based on the unique mBAAM capabilities. The intent of the new design is to reduce mass, while maintaining the dynamic stiffness. To compare the two designs, the natural frequencies and mode shapes are measured using impact testing and predicted using finite element analysis. It is confirmed that the printed structure dynamics agreed with the numerical model predictions, which demonstrates that it is feasible to model a large-scale mBAAM part and understand its behavior prior to printing. Another notable outcome of this study is that the significant residual stress and distortion in the print indicate that knowledge gaps remain for widespread implementation of mBAAM.
Benchmarking Microscale Ductility Measurements (Final Report of the Project DE-NE0008799)
Conventional macroscale experimentation is generally considered to be straightforward with few limitations. Conversely, micro/nanoscale experimentation presents numerous challenges in loading device design, sample preparation and handling, as well as accurate understanding of grain size and local texture effects on recorded measurements. Despite these challenges, nanopillar compression, MEMs based micro-tension, and nanoindentation approaches have been able to provide fundamental contributions to the understanding of material behavior at small lengthscales. However, the overarching shortcoming of these micro/nanoscale experimentation approaches, is the inability to directly translate measurements evaluated at the nm and µm length scales (e.g., hardness) to macroscale tensile material behavior (i.e., elastic modulus, yield strength, and ductility). The objectives of the proposed study are, 1) to establish best practices for obtaining tensile microscale ductility measurements, and 2) to validate methodologies to for comparing microscale ductility measurements to macroscale ductility measurements. In order to achieve these objectives, a multi-lengthscale, multi-temperature testing protocol and simulation framework are executed first on copper as a model material to validate the following approach, and second on reactor grade Zircaloy-2. Experiments are conducted on specimens extracted from the same test piece to ensure nominally identical grain size and texture from specimens to specimen. Motivated by the need to isolate the contribution of size-effects on obtained mechanical property measurements, specimens are manufactured with thicknesses at the micro- (1-10 µm), meso- (10-100's µm), and macroscales (sub-sized ASTM E8). In-situ full-field deformation techniques (scanning electron microscopy (SEM) grid methods and optical DIC) are incorporated into testing at each specimen length-scale to capture plasticity localization and evolution. Experimental testing for all specimens is conducted at both room temperature and elevated temperatures to probe the role of thermal activation on plastic deformation accommodation processes. Simulation efforts focus on examining the mechanical behavior of microscale specimens using a finite element approach with explicitly resolved grain morphologies, and an embedded crystal plasticity model. The cost-efficient implementation method allows for the modeling of a statistically significant number of both real (i.e., digital twin) and generated microstructures to obtain an understanding of the interrelationships between specimen microstructure and geometric variables (grain size, texture, specimen geometry, etc.) on microscale mechanical behavior.