Using the conformal decomposition finite element method to model the Direct Ink Write process for sinusoidal extrusion patterns
Explore the source record for details and available documents.
Kucala, Alec [Sandia National Laboratories (SNL-NM↗
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.
Multifidelity Finite Element and Surrogate Modeling of Post-Shock Miniature Springs for Optimal Experimental Design
Explore the source record for details and available documents.
Shimanek, John David [Sandia National Laboratories↗
Micro-cantilever beam experiments and modeling in porous polycrystalline UO 2
Understanding the impact of microstructure on the thermo-mechanical behavior of oxide nuclear fuels is vital to predicting their performance through multiscale models. Evaluating the mechanical properties at the sub-grain length scale is key to developing these multiscale models. In this work, 3D finite element (FE) models were constructed to simulate the micrometer-scale bending of micro-cantilever beams fabricated using porous polycrystalline uranium dioxide (UO 2 ) and tested at room temperature. Here, the results showed that the porosity and elastic anisotropy of individual grains can play a significant role in determining the effective mechanical properties of the material deduced from the tests. Specifically, the porosity had a non-negligible effect, given that the pore size was of the same order of magnitude as the dimensions of the micro-beams. Correlations between load-deflection data, pore location, and elastic properties (effective Young's modulus) were investigated using UO 2 micro-beam FE models, where pore clusters were included and placed at different locations along the length of the beam. Results indicated that the presence of pore clusters near the substrate, i.e., the clamp of the micro-cantilever beam, has the strongest effect on the load-deflection behavior, with the porosity leading to a reduction of stiffness that is the largest for any location of the pore clusters. Furthermore, it was also found that pore clusters located towards the middle of the span and close to the end of the beam have a comparatively small effect on the load-deflection behavior. Therefore, it is concluded that accurate estimates of Young's modulus can be obtained from micro-cantilever experiments after accounting for porosity on the one third of the beam length close to the clamp. This, in turn, provides an avenue to improve microscale experiments and their analysis in porous, anisotropic elastic materials.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Minimizing thickness variation in monolithic U-10Mo fuel foil and Zr interlayer during hot rolling: A microstructure-based finite element method analysis
Low-enriched uranium alloyed with 10 wt. % molybdenum (U-10Mo) has been identified as a promising alternative to highly enriched uranium fuel for the United States’ high performance research reactors. The monolithic U-10Mo fuel plate consists of a metallic U-10Mo fuel foil with a 25 µm Zr interlayer and a relatively thick cladding of aluminum alloy 6061. The Zr interlayer is typically applied during the hot co-rolling process, and this process dictates the uniformity of the Zr interlayer. Thickness variation observed in the U-10Mo and Zr interlayer has been attributed to several sources: the initial grain size of the U-10Mo castings, can materials, rolling temperature, inhomogeneous molybdenum content, and porosity in the cast U-10Mo. This thickness variation limits the ability to meet the dimensional specification; thus, a better understanding of the factors causing the nonuniform thickness is needed. In this work, we used a novel, microstructure-based finite element method to model the hot rolling process to address these concerns. Grain microstructures in U-10Mo were tessellated and explicitly considered in the finite element model. Each grain was assigned a random material property to mimic the grain strength variations induced by different grain orientations. Simulations were performed using six steel can thicknesses, four grain sizes, and with or without a Zr interlayer to investigate the influences of those variables on the thickness nonuniformity. The simulation results showed that a thinner steel can and finer U-10Mo grain size reduce thickness variations in both the U-10Mo fuel foil and Zr interlayer. The direct findings from the simulations and analysis can be used to optimize the hot rolling schedule, reduce fabrication defects, and meet the dimensional specifications. The proposed microstructure-based finite element model can be also coupled with experimental microstructure characterization data, images, and models to simulate multi-pass hot rolling.
36 MATERIALS SCIENCE↗
Advanced Permanent Magnet Generator Topologies Using Multimaterial Shape Optimization and 3D Printing: Preprint
A vast majority of utility-scale wind turbine generators in the United States are dependent on foreign- sourced rare-earth permanent magnets that are vulnerable to supply chain uncertainties. Many small wind original equipment manufacturers are motivated to pursue continuous improvements to the generator design to lower the material and production costs and improve performance by lowering cogging torque and increasing the efficiency. Traditional design and manufacturing offer limited opportunities. In this work, we demonstrate advanced design approaches for a 15-kW baseline wind turbine generator by making use of recent progress in three-dimensional (3D) printing of polymer- bonded magnets and, electrical and structural steel. We explore three methods of magnet parametrization using Bezier curves resulting in symmetric, asymmetric and multimaterial magnet designs. We employ a multiphysics approach combining parametric computer-aided design modeling, finite- element analysis and targeted sampling to identify novel designs with more opportunities for reducing rare-earth material, improving efficiency and minimizing cogging torque. The results show that asymmetric-pole design and multimaterial-pole designs offer a greater opportunity to minimize rare-earth magnet materials by up to 35% with similar performance as the baseline generator, suggesting newer opportunities with design freedom beyond traditional limits of symmetry and as allowed by 3D printing.
Bezier curves↗
Biaxial steel plated concrete constitutive models for composite structures: Implementation and validation
The Steel-plated Concrete (SC) technique is an alternative construction technique with faster onsite construction speed, reduced construction time, and increased structural performance. Aiming to predict the force transfer mechanism of SC elements, an innovative constitutive model package is proposed by implementing experimental-based biaxial steel plate concrete models into the nonlinear finite element (FE) model “Membrane Model of SC (MM-CS).” First, the formulation and implementation of the MM-SC is illustrated in detail, including the equilibrium and compatibility equations and the implementation of constitutive models. Next, various experimental data of SC members are selected and simulated using the proposed MM-SC model. In this research, two types of structures, SC panels and framed SC shear walls, and three types of loading conditions, including uniaxial compression, pure shear, and combined axial-flexural-shear tests, are analyzed respectively. Good agreements were obtained between the reported results and the FE simulation results in terms of yield capacity and ultimate capacity, proving the reliability of the implemented analytical constitutive material models and the biaxial membrane model formulations.
42 ENGINEERING↗
A layered solid finite element formulation with interlaminar enhanced displacements for the modeling of laminated composite structures
Accurate modeling of layered composite structures often requires the use of detailed finite element models which can sufficiently resolve the kinematics and material behavior within each layer of the composite. However, individually discretizing each material layer into finite elements presents a prohibitive computational expensive given the large number of thin layers comprising some laminated composites. To address these challenges, an 8-node layered solid hexahedral finite element is formulated with the aim of striking an appropriate balance between efficiency and fidelity. The element is discretized into an arbitrary number of distinct material layers, and employs reduced in-plane integration within each layer. The chosen reduced integration scheme is supplemented by a novel physical stabilization approach which includes layerwise enhancements to mitigate various forms of locking phenomena. The proposed framework additionally supports the inclusion of interlaminar enhanced displacements to better represent the kinematics of general layered composite materials. Finally, the described element formulation has been implemented in the ParaDyn finite element code, and its efficacy for modeling laminated composite structures is demonstrated on a variety of verification problems.
42 ENGINEERING↗
Simulation of a TRISO MiniFuel irradiation experiment with data-informed uncertainty quantification
An irradiation experiment using tristructural isotropic (TRISO) fuel particles and the miniature fuel (MiniFuel) irradiation vehicle was performed in Oak Ridge National Laboratory’s High Flux Isotope Reactor (HFIR) to support development of the Kairos Power fluoride salt–cooled, high-temperature reactor (KP-FHR). Here, this paper describes modeling predictions of temperatures and fuel burnup for the as-built experiment. An uncertainty quantification (UQ) analysis was performed to determine the effect of TRISO particle volume and position on the temperature predictions at various fuel heat generation rates (HGRs). This UQ study utilized fuel kernel position and volume measurements previously collected using X-ray computed tomography (XCT) techniques and Monte Carlo sampling methods to generate fuel compact cases that were then analyzed using a finite element thermal model. The UQ analysis indicated that uncertainty in calculated temperatures caused by varying TRISO particle arrangement is relatively small, even at high fuel HGR. Final predictions of particle temperatures throughout the irradiation are shown to be relevant to KP-FHR normal and off-normal operating conditions and to previous TRISO irradiation experiments. The combination of XCT with UQ analyses will inform post-irradiation examination (PIE) of the irradiated fuel compacts, and these analyses can be used to develop fuel performance models for coated particle fuel forms. Both PIE of separate-effects irradiation data and enhanced fuel performance modeling support accelerated qualification of TRISO fuels for a broad range of advanced reactor applications. The novel approach demonstrated here of measuring TRISO particle configurations with XCT methods and generating representative fuel compacts for finite element modeling and UQ analysis could be leveraged by the broader particle fuel community in the development of other TRISO fuel experiments in which these variables may have a significant impact on key outcomes.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Thermal modeling of directed energy deposition additive manufacturing using graph theory
Purpose: The purpose of this paper is to develop, apply and validate a mesh-free graph theory–based approach for rapid thermal modeling of the directed energy deposition (DED) additive manufacturing (AM) process. Design/methodology/approach: Here, the authors develop a novel mesh-free graph theory–based approach to predict the thermal history of the DED process. Subsequently, the authors validated the graph theory predicted temperature trends using experimental temperature data for DED of titanium alloy parts (Ti-6Al-4V). Temperature trends were tracked by embedding thermocouples in the substrate. The DED process was simulated using the graph theory approach, and the thermal history predictions were validated based on the data from the thermocouples. Findings: The temperature trends predicted by the graph theory approach have mean absolute percentage error of approximately 11% and root mean square error of 23°C when compared to the experimental data. Moreover, the graph theory simulation was obtained within 4 min using desktop computing resources, which is less than the build time of 25 min. By comparison, a finite element–based model required 136 min to converge to similar level of error. Research limitations/implications: This study uses data from fixed thermocouples when printing thin-wall DED parts. In the future, the authors will incorporate infrared thermal camera data from large parts. Practical implications: The DED process is particularly valuable for near-net shape manufacturing, repair and remanufacturing applications. However, DED parts are often afflicted with flaws, such as cracking and distortion. In DED, flaw formation is largely governed by the intensity and spatial distribution of heat in the part during the process, often referred to as the thermal history. Accordingly, fast and accurate thermal models to predict the thermal history are necessary to understand and preclude flaw formation. Originality/value: This paper presents a new mesh-free computational thermal modeling approach based on graph theory (network science) and applies it to DED. The approach eschews the tedious and computationally demanding meshing aspect of finite element modeling and allows rapid simulation of the thermal history in additive manufacturing. Although the graph theory has been applied to thermal modeling of laser powder bed fusion (LPBF), there are distinct phenomenological differences between DED and LPBF that necessitate substantial modifications to the graph theory approach.
42 ENGINEERING↗
Vibration of a cryomodule’s coupler during road transportation: Multibody analysis and experimental validation
The road transportation of critical components for cryogenic cryomodules is a challenging task. Indeed, the shipping and the handling of fundamental sub-assemblies can expose these components to harmful dynamic loads that can jeopardize their structural integrity. This article has the purpose of developing a finite element multibody model capable of assessing the acceleration on a test coupler for a cavity of the Single Spoke Resonator 1 subjected to a typical road transportation trip. A multibody model consisting of rigid bodies and finite element–imported flexible bodies has been created allowing a comparison of the simulated and experimental acceleration. Afterwards, an experimental on-road test, whose set-up is described in this article, has been carried out at Fermilab with instrumented components. Being able to effectively simulate and validate the dynamic effects on such sub-assemblies makes it possible, for the future, to simulate more complex structures subject to dynamic loadings, such as the entire cryomodules during shipping. Furthermore, the results from the simulations can be used to guide the design of new suspension systems suitable for the reduction mitigation of structural vibration during the road transportation of cryomodules.
42 ENGINEERING↗
RF sheath induced sputtering on Proto-MPEX. I. Sheath equivalent dielectric layer for modeling the RF sheath
The pulsed linear plasma device Prototype Material Plasma Exposure eXperiment (Proto-MPEX) uses a radio frequency (RF) helicon antenna with an aluminum nitride ceramic window for plasma production. The RF sheath created under the helicon antenna is sufficient to cause ion impact energies to be greater than the sputtering threshold of the AlN helicon window material and for impurities to be created. In this study, we investigate the RF sheath on the inner diameter of the helicon window and its impact on the impurity production rates in Proto-MPEX. Three models—a 3D COMSOL finite element RF model of the Proto-MPEX helicon region, a rectified DC sheath potential model, and the 3D Global Impurity TRansport code—are coupled together to study impurity production and transportation. This novel method of impurity generation and transport modeling spans length scales ranging from the sheath (millimeters or less) up to the full device (meters) and can be applied to other radio frequency sources and antennas in a wide range of plasma physics studies, including basic plasmas, low-temperature processing plasmas, plasma thrusters, and fusion plasmas.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Creep and Creep Fracture Modeling with Surrogate Creep Models and the Extended Finite Element Method
Alloy components in advanced nuclear reactors will be subjected to environmental conditions that could include high temperatures, irradiation, and exposure to corrosive salts. These conditions could lead to the formation of crack-like defects, which could grow over time in a mechanism known as creep crack growth (CCG). Predicting growth rates of these defects is important for assessing the safe operating life of advanced reactors. This project documents progress toward developing and testing next-generation data-driven constitutive models for deformation creep. It also documents the application of the extended finite element method in conjunction with surrogate creep models to predict CCG parameters under a variety of conditions. These important incremental developments contribute to the longer-term objective of developing microstructure-aware constitutive models that can be used for predicting creep deformation and CCG at the component scale with improved accuracy.
36 MATERIALS SCIENCE↗
A Multi-Scale Computational Platform for Predictive Modeling of Corrosion in Al-Steel Joints (Final Report)
The research team proposed to develop innovative multi-scale models to predict corrosion and the resulting mechanical performances in aluminum-steel joints. The methods of joining considered are resistance spot welding, self-piercing riveting, and rivet-welding, all suitable for mass production applications. The multi-scale models integrate high throughput first-principle calculations based on density functional theory (DFT), high throughput calculation of phase diagrams (CALPHAD) modeling, and finite element method (FEM) simulations. These models are to be validated through laboratory experiments. Furthermore, the models are available as open source so as to enable scientists and engineers in the community to adapt and contribute to the development and application. The approaches rely on the research team’s extensive experience on the prediction of properties of individual phases at finite temperatures and variable compositions through DFT calculations, and our broad expertise on dissimilar material joining and their corrosion. The proposed computational framework enables high throughput computations for improved predictions of corrosion and the associated mechanical performance in dissimilar material joints, resulting in significant reduction in computational time needed by the current state-of-the-art methods. With the participation of researchers from three universities, an auto manufacturer, two manufacturing technology/equipment suppliers, and a software developer/vendor, the interdisciplinary research team applies the technical development on both phase-based modeling and laboratory experiments into the automobile body joining processes for validation and technology demonstration. The global cost of corrosion was estimated at about 3.4% of the global GDP in 2013. By using available corrosion control practices, it is estimated a saving between 15-35% of the cost of corrosion. In the U.S., more than $276 billion is spent repairing corrosion damage. Prediction of the corrosion and its impact on performance of the dissimilar material joints is critical for reducing the massive number of the current corrosion-based recalls for automobiles. Thus, the project goal is to develop models to enable predictive maintenance and end-of-life planning of multi-metal joints with risk of corrosion under different conditions such as exposure to high temperatures in summer and salt solutions in winter, quantified through its pH. An academia-industry consortium led by the University of Michigan and including Pennsylvania State University, University of Illinois Urbana-Champaign, University of Georgia, General Motors Company, Livermore Software Technology Corporation, and Optimal Process Technologies, LLC. created multi-scale models for prediction of corrosion in aluminum-steel joint structures such of them used in vehicle subassemblies – chassis and transmission systems. Starting from the first principle calculations, the team developed mathematical and data-driven models to predict the metallic components, which are formed during joining of two metals, for example aluminum and steel - a lightweight multilateral system which is currently used in more than 60% car bodies. These models were used for simulating chemical reactions that are happening when the joining metallic components are exposed to high temperatures and different pH values. The team was able to predict how the corrosion installs on the metallic components and how they lead to a sudden failure of components in cars. Newly developed machine learning algorithms combining Science, Technology, Engineering and Math disciplines, advanced finite element simulation and experimental validations have been integrated in a platform for prediction of the corrosion evolution and prediction the failure of joints under mechanical loadings and fatigue. Moreover, based on machine learning and inverse analysis, the team proposed solutions for designing new metallic alloys less susceptible to corrosion when joining multi-material assembles. An average of 4% error compared with experiments was achieved for the most common joints that are used in vehicle subassemblies.
36 MATERIALS SCIENCE↗
Pressure Transient Analyses and Poroelastic Modeling of Hydraulic Fracture Dilation for Multiple Injections at the Devine Fracture Pilot Site
Our team has conducted electromagnetic (EM) surveys for the past six years to monitor hydraulic-fracture behavior at the Devine Fracture Pilot Site (DFPS). The sub-horizontal orientation of a shallow hydraulic fracture at the DFPS provides uniform access to the fracture area for interrogation and data collection. Ahmadian et al. (2023) suggested a possible correlation between spatiotemporal changes in the flow rate, bottomhole pressure (BHP), and the observed surface recorded electric field at the DFPS. In this paper, we present the development of poroelastic forward models and pressure transient analyses (PTAs) to support the development of a multiphysics inverse model for these EM surveys. First, we conducted PTAs of the shut-in periods after six injections out of 10 to determine the fracture closure pressure (FCP) or the overburden pressure used in a poroelastic fracture reopening model. Second, we developed a finite-element poroelastic model throughout five injection cycles to include the effect of the cumulative injected volumes due to the previous injections on current fracture dilation in the presence of highly permeable unpropped and propped zones adjacent to the cohesive layer that models fracture reopening. Fracture reopening in this poroelastic model is based on a calibrated traction-separation response using the bottomhole pressure collected in two injection campaigns in 2020 and 2022. We used the outcomes of a previous simulation study of the primary hydraulic-fracturing stimulation to define the dimension of an unpropped fracture zone ahead of the propped fracture area. The PTAs led to FCPs consistent with those obtained using the injection data collected at the DFPS in 2020. Further, these analyses showed that at later injections, the fracture closure occurred at a later time with respect to the shut-in time, inferring the effect of cumulative injected volumes in previous injections. The simulation results show that considering the propped and unpropped fracture zones improves our poroelastic model in predicting the injection-well BHP. The numerical simulation results demonstrate a significant excess pore pressure near the fracture because of the preceding formation loadings by the previous injections. The obtained fracture dilation area and fluid pressure distribution provide a basis to improve the development of a multiphysics inverse model. Furthermore, in an iteratively coupled scheme, this pressure distribution can be introduced into EM models to render a holistic view of the causative mechanisms for the surface signal anomalies.
02 PETROLEUM↗
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↗
Multiscale Modeling of the Mechanical Response of Silicon Carbide Composite Within the Accelerated Fuel Qualification Framework
The accelerated fuel qualification (AFQ) framework has been used for the initial development of multiscale modeling of silicon carbide (SiC) fiber reinforced composite (SiC-SiC). The AFQ framework provides a methodology to leverage physics-informed multiscale modeling along with a reduced set of empirical test data to reduce the time and cost of licensing and qualification of new nuclear fuel systems while maintaining the overall nuclear power plant safety case. SiC-SiC is being proposed for in-core applications, most notably fuel cladding, for current and next-generation nuclear reactors because of its high temperature stability, irradiation tolerance, and ability to withstand many accident conditions. As these composites exhibit multiscale architectures and complex microstructure-based fracture mechanics, it is an appealing use case for the AFQ methodology. While the end goal of this work is a single multiscale model that can be used for predictive in-core performance, current focus is on the individual various length scale models. Four individual models have been initially developed from microscale to engineering system level to capture key physics-based effects across different length scales. These models include a microscale homogenized tow model, a mesoscale fast Fourier transform–based weave model that integrates the homogenized tow model, a mesoscale finite element–based weave model, and a system-level BISON fuel performance model. Results of these models have undergone an initial comparison with separate-effects test data showing a good match to experimental results. By using the AFQ framework during model development, several near-term benefits have been secured including a reduction in development time for the SiC-SiC cladding, more targeted irradiation testing, and a better understanding of uncertainty.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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↗
PRISMS-PF: A general framework for phase-field modeling with a matrix-free finite element method
Abstract A new phase-field modeling framework with an emphasis on performance, flexibility, and ease of use is presented. Foremost among the strategies employed to fulfill these objectives are the use of a matrix-free finite element method and a modular, application-centric code structure. This approach is implemented in the new open-source PRISMS-PF framework. Its performance is enabled by the combination of a matrix-free variant of the finite element method with adaptive mesh refinement, explicit time integration, and multilevel parallelism. Benchmark testing with a particle growth problem shows PRISMS-PF with adaptive mesh refinement and higher-order elements to be up to 12 times faster than a finite difference code employing a second-order-accurate spatial discretization and first-order-accurate explicit time integration. Furthermore, for a two-dimensional solidification benchmark problem, the performance of PRISMS-PF meets or exceeds that of phase-field frameworks that focus on implicit/semi-implicit time stepping, even though the benchmark problem’s small computational size reduces the scalability advantage of explicit time-integration schemes. PRISMS-PF supports an arbitrary number of coupled governing equations. The code structure simplifies the modification of these governing equations by separating their definition from the implementation of the numerical methods used to solve them. As part of its modular design, the framework includes functionality for nucleation and polycrystalline systems available in any application to further broaden the phenomena that can be used to study. The versatility of this approach is demonstrated with examples from several common types of phase-field simulations, including coarsening subsequent to spinodal decomposition, solidification, precipitation, grain growth, and corrosion.
36 MATERIALS SCIENCE↗