Results for “FINITE ELEMENT METHOD” — Page 43

DOE OSTI2021

SAM User’s Guide

The System Analysis Module (SAM) is a modern system analysis tool being developed at Argonne National Laboratory for advanced non-LWR safety analysis. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for Sodium-cooled Fast Reactors (SFR), Lead-cooled Fast Reactors (LFR), and Molten Salt Reactors (MSR) or Fluoride-cooled High-temperature Reactors (FHR). SAM takes advantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented application framework (MOOSE), and its underlying meshing and finite-element library (libMesh) and linear and non-linear solvers (PETSc), to leverage the modern advanced software environments and numerical methods. This document provides a user’s guide, which will help users understand the input description and core capabilities of the SAM code. A brief overview of the code is presented, as well as how to obtain and run it. The input syntax for various parts of the code is provided. Additionally, a number of example problems, starting with simple unit component problems to problems with increasing complexity, are provided. Because the code is still under active development, this SAM User’s Guide will evolve with periodic updates.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

SAM User's Guide

The System Analysis Module (SAM) is a modern system analysis tool being developed at Argonne National Laboratory for advanced non-LWR safety analysis. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for Sodium-cooled Fast Reactors (SFR), Lead-cooled Fast Reactors (LFR), and Molten Salt Reactors (MSR) or Fluoride-cooled High-temperature Reactors (FHR). SAM takes advantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented application framework (MOOSE), and its underlying meshing and finite-element library (libMesh) and linear and non-linear solvers (PETSc), to leverage the modern advanced software environments and numerical methods. This document provides a user’s guide, which will help users understand the input description and core capabilities of the SAM code. A brief overview of the code is presented, as well as how to obtain and run it. The input syntax for various parts of the code is provided. Additionally, a number of example problems, starting with simple unit component problems to problems with increasing complexity, are provided. Because the code is still under active development, this SAM User’s Guide will evolve with periodic updates.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Hybrid Basis and Multi-Center Grid Method for Strong-Field Processes

We present a time-dependent framework that combines a hybrid basis, consisting of Gaussian-type orbitals (GTOs) and finite-element discrete-variable representation (FEDVR) functions, with a multicenter grid to simulate strong-field and attosecond dynamics in atoms and molecules. The method incorporates the construction of the orthonormal hybrid basis, the evaluation of electronic integrals, a unitary time-propagation scheme, and the extraction of optical and photoelectron observables. Its accuracy and robustness are benchmarked on one-electron systems such as atomic hydrogen and the dihydrogen cation (H$^+_2$) through comparisons with essentially-exact reference results for bound-state energies, high-harmonic generation spectra, photoionization cross sections, and photoelectron momentum distributions. This work establishes the groundwork for its integration with quantum-chemistry methods, which is already operational but will be detailed in future work, thereby enabling ab initio simulations of correlated polyatomic systems in intense ultrafast laser fields.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Non-linear Wire Rope Isolator Model to Enhance Transportation Simulation of Fragile Equipment

Abstract Introduction Wire rope isolators are often used as means of vibration isolation for the transportation of fragile machine equipment. When deployed as suspension systems in a transportation scenario, a simple spring–damper model is often used to predict the dynamic loads that act on the structure being transported. However, this simple model has proved to be too simplistic to be used in the development of an accurate numerical model that can predict the vibration levels experienced by the inner components of fragile equipment during its transportation. This paper describes the experimental tests conducted on a wire rope isolator used for the transportation of the prototype SSR1 cryomodule. Materials and methods A hysteretic Bouc–Wen model has been used to analytically describe the force–deformation relationship of the wire rope isolators. The developed model of the isolator has been implemented in a larger model to simulate the actual transportation of the prototype SSR1 cryomodule, a section of the new PIP-II linear accelerator under construction at Fermilab. A series of multibody dynamic simulations with rigid and flexible components was used to numerically determine the acceleration of some critical components. Results and conclusions An actual experimental transportation was simulated using two numerical models: the developed Bouc–Wen model and a conventional spring–damper model. It is shown how the Bouc–Wen formulation of the isolator characteristics drastically improves the correspondence between experimental and simulated results if compared to a spring–damper model, especially in the range of 0–30 Hz which is the most critical for transportation problems.

42 ENGINEERING↗

DOE OSTI2022

Enhanced resonant ultrasound spectroscopy for measurement of the elastic properties of multi-material systems

Understanding the elastic properties of materials is critical for their safe incorporation and predictable performance. Current methods of bulk elastic characterization often have notable limitations for in situ structural applications, with usage restricted to simple geometries and material distributions. To address these existing issues, this study sought to expand the capabilities of resonant ultrasound spectroscopy (RUS), an established nondestructive evaluation method, to include the characterization of isotropic multi-material samples. In this work, finite-element-based RUS analysis consisted of numerical simulations and experimental testing of composite samples comprised of material pairs with varying elasticity and density contrasts. Utilizing genetic algorithm inversion and mode matching, our results demonstrate that elastic properties of multi-material samples can be reliably identified within several percent of known or nominal values using a minimum number of identified resonance modes, given sample mass is held consistent. The accurate recovery of material properties for composite samples of varying material similarity and geometry expands the pool of viable samples for RUS and advances the method towards in situ inspection and evaluation.

36 MATERIALS SCIENCE↗

DOE OSTI2021

An FFT-based approach for Bloch wave analysis: application to polycrystals

A method based on the Fast Fourier Transform is proposed to obtain the dispersion relation of acoustic waves in heterogeneous periodic media with arbitrary microstructures. The microstructure is explicitly considered using a voxelized Representative Volume Element (RVE). The dispersion diagram is obtained solving an eigenvalue problem for Bloch waves in Fourier space. To this aim, two linear operators representing stiffness and mass are defined through the use of differential operators in Fourier space. The smallest eigenvalues are obtained using the implicitly restarted Lanczos and the subspace iteration methods, and the required inverse of the stiffness operator is done using the conjugate gradient with a preconditioner. The method is used to study the propagation of acoustic waves in elastic polycrystals, showing the strong effect of crystal anistropy and polycrystaline texture on the propagation. It is shown that the method combines the simplicity of classical Fourier series analysis with the versatility of Finite Elements to account for complex geometries proving an efficient and general approach which allows the use of large RVEs in 3D.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Multi-element flow-driven spectral chaos (ME-FSC) method for uncertainty quantification of dynamical systems

The flow-driven spectral chaos (FSC) is a recently developed method for tracking and quantifying uncertainties in the long-time response of stochastic dynamical systems using the spectral approach. The method uses a novel concept called enriched stochastic flow maps as a means to construct an evolving finite-dimensional random function space that is both accurate and computationally efficient in time. In this paper, we present a multi-element version of the FSC method (the ME-FSC method for short) to tackle (mainly) those dynamical systems that are inherently discontinuous over the probability space. In ME-FSC, the random domain is partitioned into several elements, and then the problem is solved separately on each random element using the FSC method. Subsequently, results are aggregated to compute the probability moments of interest using the law of total probability. To demonstrate the effectiveness of the ME-FSC method in dealing with discontinuities and long-time integration of stochastic dynamical systems, four representative numerical examples are presented in this paper, including the Van-der-Pol oscillator problem and the Kraichnan-Orszag three-mode problem. Results show that the ME-FSC method is capable of solving problems that have strong nonlinear dependencies over the probability space, both reliably and at low computational cost.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Fourier Analysis Method for Wireless Power Transfer Coil Design

Wireless power transfer (WPT) systems for Electric Vehicles (EVs) are designed to meet specifications such as stray field, power transfer, efficiency, and ground clearance. Typical design approaches include iterative analysis of predetermined coil geometries to identify candidates that meet these constraints. This work instead directly generates WPT coil shapes and magnetic fields to meet specifications and constraints through the optimization of Fourier basis function coefficients. The proposed Fourier Analysis Method (FAM) applies to arbitrary planar coil geometries and does not rely on iterative finite-element analysis (FEA) simulations. This flexibility allows for rapid design evaluation across a larger range of coil geometries and design specifications. A prototype coil is built to compare FAM outputs to experimental measurements and FEA simulations. The FAM is then used to illustrate the tradeoff of coil current and stray field for a given power level showing that the method is capable of generating optimized coil shapes to meet arbitrary field constraints.

Foote, Andrew↗

DOE OSTI2021

Ball-on-ring test validation for equibiaxial flexural strength testing of engineered ceramics

The validation of a ball-on-ring, equibiaxial flexural strength method to obtain the transverse rupture strength (TRS) of right cylindrical ceramic specimens was performed in this study. Validation of the test method was achieved using commercially available engineered high purity alumina disks and finite element (FE) model analysis. The validated fixture was then used to obtain the TRS and Weibull statistical analysis of MgO-partially stabilized zirconia (MSZ) and Y 2 O 3 -partially stabilized zirconia (YSZ) ceramic disks. TRS data for alumina, MSZ, and YSZ agreed with the TRS values reported in the literature. A statistically relevant number of samples (N > 30) for each material were tested to allow for a Weibull statistical analysis. Weibull parameters for these materials were within the expected values for engineered ceramics. The characteristic strength for alumina, MSZ, and YSZ were determined to be 289, 786, and 814 MPa, respectively. The Weibull modulus was determined between 10 and 25 for each material, which is typical of engineered ceramics. In addition, FE model results were in close agreement with experimental fracture values for the three ceramic materials tested in this study.

36 MATERIALS SCIENCE↗

DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗

DOE OSTI2021

High-throughput bend-strengths of ultra-small polysilicon MEMS components

The strength distribution of polysilicon bend specimens, approximately 10 μm in size, is measured using a high-throughput microelectromechanical system fabrication and testing method. The distribution is predicted from reference tests on tensile specimens and finite element analysis of the bend specimen geometry incorporated into a stochastic extreme-value strength framework. Agreement between experiment and prediction suggests that the ultra-small specimens may be at the limit of extreme-value scaling and contain only one strength-controlling flaw/specimen.

Cook, Robert F. (ORCID:0000000304228881)↗

DOE OSTI2020

System, method, and computer program for creating geometry-compliant lattice structures

A system and method of creating a shape-conforming lattice structure for a part formed via additive manufacturing. The method includes receiving a computer model of the part and generating a finite element mesh. A lattice structure including a number of lattice cellular components may also be generated. Some of the mesh elements of the finite element mesh may be deformed so that the finite element mesh conforms to the overall shape of the part. The lattice structure may then be deformed so that the lattice structure has a cellular periodicity corresponding to the finite elements of the finite element mesh. In this way, the part retains the benefits of its overall shape and the benefits of lattice features without introducing structural weak points, directional stresses, and other structural deficiencies.

36 MATERIALS SCIENCE↗

DOE OSTI2021

System, method, and computer program for creating geometry-compliant lattice structures

A system and method of creating a shape-conforming lattice structure for a part formed via additive manufacturing. The method includes receiving a computer model of the part and generating a finite element mesh. A lattice structure including a number of lattice cellular components may also be generated. Some of the mesh elements of the finite element mesh may be deformed so that the finite element mesh conforms to the overall shape of the part. The lattice structure may then be deformed so that the lattice structure has a cellular periodicity corresponding to the finite elements of the finite element mesh. In this way, the part retains the benefits of its overall shape and the benefits of lattice features without introducing structural weak points, directional stresses, and other structural deficiencies.

Vernon, Gregory John↗

DOE OSTI2024

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗

DOE OSTI2024

System, method, and computer program for creating geometry-compliant lattice structures

A system and method of creating a shape-conforming lattice structure for a part formed via additive manufacturing. The method includes receiving a computer model of the part and generating a finite element mesh. A lattice structure including a number of lattice cellular components may also be generated. Some of the mesh elements of the finite element mesh may be deformed so that the finite element mesh conforms to the overall shape of the part. The lattice structure may then be deformed so that the lattice structure has a cellular periodicity corresponding to the finite elements of the finite element mesh. In this way, the part retains the benefits of its overall shape and the benefits of lattice features without introducing structural weak points, directional stresses, and other structural deficiencies.

Vernon, Gregory John↗

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 OSTI2022

Deformation twinning and detwinning in extruded Mg-4Al: In-situ experiment and crystal plasticity simulation

Deformation twinning and detwinning in extruded Mg-4Al were investigated using in-situ SEM-DIC experiments and crystal plasticity finite element (CPFE) simulation. In this study, the in-situ SEM-DIC method was used to provide a unique set of data including twin/detwin characteristics and twin area fraction in addition to strain maps. A statistical analysis of the activation of twin variants and twin area fraction during both twinning and detwinning was conducted. A strong correlation was found between twin growth/shrinkage and the Schmid Factor (SF) for individual twin variants Higher twin SF during loading and unloading led to higher twin growth and shrinkage, respectively. However, after the applied compressive strain was removed, the pattern of the twin area fractions of the residual twin variants versus their nominal SFs did not follow the trend observed at the maximum compressive strain. Using a systematic methodology and an advanced twin/detwin model, the PRISMS-Plasticity CPFE simulation was calibrated using experimentally determined stress versus strain and twin area fraction versus strain information. A comprehensive evaluation of the CPFE model was conducted to determine its ability to capture the statistics of twin variants activation and twin area fraction. CPFE accurately captured the statistical aspects of both twinning and detwinning. Finally, it also predicted the first dominant twin variant for 47.5% of the grains and at least one of the two dominant twin variants for 80% of the grains at maximum compressive strain.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Designing Cladded Components for High Temperature Nuclear Service: Part I—Analysis Methods

Use of corrosion-resistant cladding can greatly extend the design life of structural components in many advanced reactor systems. However, there are currently no ASME design rules for cladded components to guard against creep-fatigue failure and ratcheting strain accumulation in elevated temperature nuclear service. This paper, presented in two parts, addresses this gap by proposing a design strategy for cladded components that do not require long-term testing of clad materials. The proposed approach relies on approximate design analysis methods for two types of clad materials—a soft clad that creeps much faster than and has lower yield stress than the class A base material and a hard clad that creeps much slower than and has higher yield stress than the class A base material. Part I discusses the approximate analysis strategies for the clad materials—treat a soft clad as perfectly compliant and a hard clad as linear elastic—and Part II develops a complete set of design rules for each of the two types of cladded components. In conclusion, finite element analyses of representative high temperature reactor components show that the proposed design analysis methods can bound the design quantities in soft cladded components and approximate the design quantities in hard cladded components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗